基于Jurkevich法Cyg X-1光变周期特性研究

针对天文观测数据不等间隔的特点,利用Jurkevich方法对黑洞候选体X射线双星Cyg X-1近16年的Rossi X-ray Timing Explorer(RXTE)All-Sky Monitor(ASM)数据进行了频谱分析,并结合Kidger周期存在可信度的判据,发现Cyg X-1存在363±11.5 d的光变周期,但没有发现420 d的光变周期,而且150 d的光变周期可能仅仅是51±1.4 d光变周期的观测表象。     

黑洞候选体X射线双星CygX-1长期谱变化的研究

利用Rossi X-ray Timing Explorer(RXTE)卫星从1996年1月5日到2005年5月1日的all-sky monitor(ASM,1.5-12KeV)资料,详细分析了黑洞候选体X射线双星Cyg X-1在A-band(1.5-3KeV),B-band(3-5KeV)和C-band(5-12KeV)三个能段上光子计数率与谱线硬度比HR2(5-12KeV/3-5KeV)的相关性.利用1天时标的ASM资料研究发现: (1)Cyg X-1处于软态时,在A-band上光子计数率与硬度比HR2呈现反相关性,而在B-band和C-band上呈现正相关性.当该天体处于硬态时,在A、B、C三个能段上光子计数率与硬度比HR2都是反相关的.(2)Cyg X-1处于硬态和软态时,硬度比HR2与HR1都是正相关的.此外,还分析了ASM上"Dwell by Dwell"资料,得到了十分有趣的结果,即: (1)MJD=52600-52760期间(此时Cyg X-1处于硬态),A-band和B-band上的光子计数率与HR2是反相关性的,而在C-band上却呈现出较强的正相关性. (2)硬态时,硬度比HR2与HR1存在明显的反相关性.     

黑洞候选体X射线双星Cyg X-1光变曲线的研究

利用Rossi X-ray Timing Explorer(RXTE)卫星从1996年1月到2005年5月的公共资料,研究了黑洞候选体X射线双星Cyg X-1光变曲线的周期性,并且发现了一些有趣的周期特性(T=1．0±0．2天和T=18．0±3．0天)．在硬态时,T=1．0±0．2天和T=18．0±3．0天的周期特性同时出现在它的光变曲线之中,但在软态时仅有T=1．0±0．2天的周期出现．采用一种新的方法(即微分法),同时为了检验微分法的正确性,还利用传统的方法(即快速傅立叶变换和叠加法)分析了相同的资料,并得到了与微分法相同的结果．另外也对Cyg X-1的轨道周期进行了研究,结果发现T=5．6天的轨道周期不仅出现在硬态,同时也出现在软态,只是在硬态时比在软态时明显．     

粘滞对黑洞吸积盘的截断的影响

采用吸积盘-冕模型研究了粘滞对黑洞X射线双星和低光度活动星系核(LLAGN)中吸积盘的截断以及黑洞X射线双星中高低态转变的影响.以前的分析表明,冕的结构对粘滞的大小非常敏感.因此详细计算了一系列粘滞系数情况下冕的结构.为了便于与观测比较,将数值计算结果进行解析拟合得到最大蒸发率和粘滞系数a的关系,M/MEDD≈1.08a3.35;截断半径和粘滞系数的关系,R/Rs≈36.11a-1.94.这些结果可以用来解释光谱态的高低态的转换和截断半径的变化.并将这些结果应用到几个黑洞X射线双星XTE J118 480,GX 339-4,以及活动星系核NGC 4636中.     

天鹅座X-1硬X射线的能量谱

1985年9月22日,我们用空间硬X射线望远镜HAPI-2,在我国的高空气球上,对天鹅座X-1进行了观测,望远镜几何面积为140cm~2,观测能区20—200keV,空间定向精度±0.2°。获得的硬X射线能谱陡,总辐射能流LX～1.0×1~(37)erg/sec,低于常规态。本文将介绍望远镜的性能和观测过程,能谱分析的结果,并对当时天鹅座X-1可能的“态”进行讨论。     

Circinus X-1 在转换期的能谱和功率谱

Circinus X-1是小质量X射线双星(LMXB),它的能谱和时变十分复杂而独特。在不同的观测时期,它的硬度强度图(HID)和双色图(CCD)显示了不同的形状。本文选择了Circinus X-1正处于X射线流量从高向低转换期间的观测数据,通过对这一时期HID和相应能谱和时变性质的讨论,并将所得结果与1997年观测到的高亮度期间的能谱和时变性质(Shirey 1999)进行对比研究,进一步找到源的强度变化对其X射线辐射性质的影响。     

XTE J1550-564 2001年“迷你爆发”的X射线能谱研究

黑洞X射线暂现源的迷你爆发是一类峰值光度较低、持续时间较短的爆发.由于观测数据较少,其物理机制仍不清楚.利用RXTE (Rossi X-ray Timing Explorer)卫星从2001年1月28日到3月14日的数据,研究了黑洞X射线暂现源XTE J1550–564 2001年迷你爆发的X射线能谱特性.发现在本次迷你爆发中, XTE J1550–564的X射线能谱可以用幂律谱很好地拟合.整个爆发的硬度强度图不是标准的q型,而是一直保持在最右侧.此外,还分析了此次爆发谱指数Γ与未吸收的2–10 keV能段的X射线流量F_(2–10 keV)的相关性,发现Γ-F_(2–10 keV)呈反相关关系,且谱指数Γ∈[1.35, 1.72].上述结果表明2001年这次爆发一直处于低/硬态,它的X射线辐射主要来自于辐射低效的吸积模式,如ADAF(Advection-Dominated Accretion Flow).     

Cyg X-3X光子场对γ射线的吸收效应

本文详细讨论了X射线双星Cyg X-3 的X光子场对γ射线的吸收效应.计算结果表明:当Cyg x-3 的X光子发射场处于高态时,从致密中子星附近发射的能量为10~2—10~4MeV的γ射线约有60％被其吸收,吸收过程产生的正负电子的逆Compton 散射会使出射能量10—70MeV的γ射线强度平均增加约50％.计算的γ射线能谱与实验观测结果基本相符.     

黑洞天体活动的基本面研究

<正>活动星系核(中心黑洞质量M_(BH)～10~6-10~(10)M_⊙)和黑洞X射线双星(M_(BH)～10M_⊙)普遍被认为具有相似的中心引擎:黑洞、吸积盘和喷流.类似的中心引擎、质量相却差如此之大(6–9个数量级)的两类黑洞系统是否具有相似的物理仍不清楚.本文围绕不同尺度黑洞天体的物理性质和观测特性展开,主要研究了不同尺度黑洞天体活动的基本面关系以及黑洞X射线双星的能谱演化.     

银河系甚弱X射线暂现源的研究进展

银河系内的X射线双星暂现源大多是吸积的黑洞或中子星系统。最近十几年以来,还探测到几十颗具有极低爆发光度(2～10 keV的光度约10~(27)～10~(29) J·s~(-1))的X射线暂现源,称为甚弱X射线暂现源(very faint X-ray transient,简称VFXT)。根据观测和理论分析,VFXT的长期平均物质吸积率约小于10~(-13)～10~(-12) M_⊙·a~(-1),这样低的吸积率无法用传统的X射线双星演化理论解释。首先总结了VFXT的观测特征,指出其族群多样性的特征。评述了现有可能的机制,并指出还需要更多的观测和理论研究来揭示这类奇特的暂现源的本质。     

On the origin of the solar system and the exceptional position of the Sun in the Galaxy

The solar system's position in the Galaxy is an exclusive one, since the Sun is close to the corotation circle, which is the place where the angular velocity of the galactic differential rotation is equal to that of density waves displaying as spiral arms. Each galaxy contains only one corotation circle; therefore, it is an exceptional place. In the Galaxy, the deviation of the Sun from the corotation is very small — it is equal to ΔR/R _⊙≈0.03, where ΔR=R _c ?R _⊙,R _c is the corotation distance from the galactic center andR _⊙ is the Sun's distance from the galactic center. The special conditions of the Sun's position in the Galaxy explain the origin of the fundamental cosmogony timescalesT ₁≈4.6×10⁹ yr,T ₂?10⁸ yr,T ₃?10⁶ yr detected by the radioactive decay of various nuclides. The timescaleT ₁ (the solar system's ‘lifetime’) is the protosolar cloud lifetime in a space between the galactic spiral arms. The timescaleT ₂ is the presolar cloud lifetime in a spiral arm.T ₃ is a timescale of hydrodynamical processes of a cloud-wave interaction. The possibility of the natural explanation of the cosmogony timescales by the unified process (on condition that the Sun is near the state of corotation) can become an argument in favour of the fact that the nearness to the corotation is necessary for the formation of systems similar to the Solar system. If the special position of the Sun is not incidental, then the corotation circles of our Galaxy, as well as those of other galaxies, are just regions where situations similar to ours are likely to be found.     

Perturbations by the oblateness of the Earth and by the planets in the motion of the Moon

Perturbations in the motion of the Moon are computed for the effect by the oblateness of the Earth and for the indirect effect of planets. Based on Delaunay's analytical solution of the main problem, the computations are performed by a method of Fourier series operation. The effect of the oblateness of the Earth is obtained to the second order, partly adopting an analytical evaluation. Both in longitude and latitude are found a few terms whose coefficient differs from the current lunar ephemeris based on Brown's theory by about 0.01. While, concerning the indirect effect of planets, several periodic terms in the current ephemeris seem to have errors reaching 0.05.As for the secular variations of and due to the figure of the Earth and the indirect effect of planets, the newly-computed values agree within 1/cy with Brown's results reduced to the same values of the parameters. Further, the accelerations in the mean longitude, and caused by the secular changes in the eccentricity of the Earth's orbite and in the obliquity of the ecliptic are obtained. The comparison with Brown shows an agreement within 0.3/cy² for the former cause and 0.02/cy² for the latter. An error is found in the argument of the principal term for the perturbations due to the ecliptic motion in the current ephemeris.Proceedings of the Conference on Analytical Methods and Ephemerides: Theory and Observations of the Moon and Planets. Facultés universitaires Notre Dame de la Paix, Namur, Belgium, 28–31 July, 1980.     

On the model of the accumulation of the moon compatible with the data on the composition and the age of lunar rocks

It is suggested that the overall early melting of the lunar surface is not necessary for the explanation of facts and that the structure of highlands is more complicated than a solidified anorthositic ‘plot’. The early heating of the interior of the Moon up to 1000K is really needed for the subsequent thermal history with the maximum melting 3.5 × 10⁹ yr ago, to give the observed ages for mare basalts. This may be considered as an indication that the Moon during the accumulation retained a portion of its gravitational energy converted into heat, which may occur only at rapid processes. A rapid (t < 10³ yr) accretion of the Moon from the circumterrestrial swarm of small particles would give necessary temperature, but it is not compatible with the characteristic time 10⁸ yr of the replenishment of this swarm which is the same as the time-scale of the accumulation of the Earth. It is shown that there were conditions in the circumterrestial swarm for the formation at a first stage of a few large protomoons. Their number and position is evaluated from the simple formal laws of the growth of satellites in the vicinity of a planet. Such ‘systems’ of protomoons are compared with the observed multiple systems, and the conclusion is reached that there could have been not more than 2–3 large protomoons with the Earth. The tidal evolution of protomoon orbits was short not only for the present value of the tidal phase-lag but also for a considerably smaller value. The coalescence of protomoons into a single Moon had to occur before the formation of the observed relief on the Moon. If we accept the age 3.9 × 10⁹ yr for the excavation of the Imbrium basin and ascribe the latter to the impact of an Earth satellite, this collision had to be roughly at 30R, whereR is the radius of the Earth, because the Moon at that time had to be somewhere at this distance. Therefore, the protomoons had to be orbiting inside 20–25R, and their coalescence had to occur more than 4.0x10⁹ yr ago. The energy release at coalescence is equivalent to several hundred degrees and even 1000 K. The process is very rapid (of the order of one hour). Therefore, the model is valid for the initial conditions of the Moon.     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

Accounting for the viscosity of the fluid core in the problem of the physical libration of the moon

A two-component theoretical model of the physical libration of the Moon in longitude is constructed with account taken of the viscosity of the core. In the new version, a hydrodynamic problem of motion of a fluid filling a solid rotating shell is solved. It is found that surfaces of equal angular velocity are spherical, and a velocity field of the fluid core of the Moon is described by elementary functions. A distribution of the internal pressure in the core is found. An angular momentum exchange between the fluid core and solid mantle is described by a third-order differential equation with a right-hand side. The roots of a characteristic equation are studied and the stability of rotation is proved. A libration angle as a function of time is found using the derived solution of the differential equation. Limiting cases of infinitely large and infinitely small viscosity are considered and an effect of lag of a libration phase from a phase of action of an external moment of forces is ascertained. This makes it possible to estimate the viscosity and sizes of the lunar fluid core from data of observations.     

On some problems concerning the calculation of the form of the magnetopause and the electric field on the magnetopause in the framework of the continuum medium model

In order to understand the reason of the existence of the electric field in the magnetosphere, and for the theoretical evaluation of its value, it is necessary to find the solution of the problem of determination of the magnetosphere boundary form in the frameworks of the continuum medium model which takes into account part of the magnetospheric plasma movement in supporting the magnetospheric boundary equilibrium. A number of problems for finding the distribution of the pressure, the density, the magnetic field and the electric field on the particular tangential discontinuity is considered in the case when the form of discontinuity is set (the direct problem) and a number of problems for finding the form of the discontinuity and the distribution of the above-mentioned physical quantities on the discontinuity is considered when the law of the change of the external pressure along the boundary is set (for example, with the help of the approximate Newton equation). The problem which is considered here, which deals with the calculation of the boundary form and with the calculation of the distribution of the corresponding physical quantities on the discontinuity of the 1st kind for the compressible fluid with the magnetic field with field lines which are perpendicular to the plane of the flow in question, concerns the last sort of problems. The comparison of the results of the calculation with the data in the equatorial cross-section of the magnetosphere demonstrates that the calculated form of the boundary, the value of the velocity of the return flow and the value of the electric field on the magnetopause, agree satisfactorily with the observational data.