大质量黑洞吸积盘的径向自引力

本文研究了大质量黑洞吸积盘的自引力,用薄盘位形上积分的方法计算了吸积盘自引力的径向与垂向分量,着重讨论了径向自引力。主要结果为:对于大质量黑洞(M～10~8—10~(10)M_⊙)吸积盘,在(R/R_g)～10~5—10~4的距离上,径向自引力会超过中心天体引力。在这个距离上,吸积盘的动力学结构完全不同于开普勒盘。提出了径向自引力不稳定扰动作为一种能源机制。本文还得到吸积盘自引力与中心天体引力量级比较的两个判据,并由此得到大质量黑洞吸积盘外半径的近似解析估计。本文结果可用于类星体和星系核吸积盘。     

吸积盘随机振荡光度功率谱密度中低频准周期振荡现象的研究

采用含有频率涨落噪声和指数形式关联随机力作用的广义朗之万(Langevin)方程模型描述黑洞吸积盘的垂向振荡,推导出吸积盘随机振荡光度功率谱密度的解析表达式,并讨论了系统参数对功率谱密度中低频准周期振荡(Low Frequency Quasi-Periodic Oscillations,LFQPOs)现象的影响。研究结果发现选取合适的系统参数时,功率谱密度曲线上出现了一个基频和一个二次谐频的共振双峰低频准周期振荡,基频峰对应的中心频率为吸积盘振荡的特征频率;随机力关联时间决定了基频峰的高度和宽度,频率噪声强度和粘滞阻尼只对二次谐频峰产生影响。结果说明吸积盘的随机振荡模型可以作为低频准周期振荡起源的一种解释。     

几何薄光学厚吸积盘谱线轮廓及成像研究

假设位于黑洞赤道面上做圆形轨道运动的吸积盘是几何薄、光学厚的.利用光子追踪法计算在Kerr度规下的光子运动轨迹,通过数值计算研究薄吸积盘的相对论谱线轮廓及成像.在大角度观测时,吸积盘下表面的光子对谱线轮廓及成像的影响是显著的.     

射电类星体B3 1637+436A双峰轮廓宽发射线的光变研究

通过对比间隔17yr (静止系12.5yr)的Sloan Digital Sky Survey (SDSS)、Keck LRIS (Low-Resolution Imaging Spectrograph)和Baryon Oscillation Spectroscopic Survey (BOSS)的光学光谱, 对宽线双峰射电类星体B3 1637+436A的吸积盘性质做了初步研究.该类星体的SDSS (2001年5月21日)和Keck LRIS (2003年6月28\lk 日)的连续谱没有明显差异, 而BOSS (2018月6月17日)连续谱在V波段比SDSS暗了约1.2星等.这3条光谱的Hα宽发射线均表现出明显的双峰轮廓. 和连续谱一样, 宽线双峰的发射线轮廓也被认为是来自于吸积盘.通过上述3条光谱的Hα宽线双峰的盘模型拟合发现: BOSS光谱的Hα宽线双峰可以用单个盘模型很好地拟合,谱线发射区距离中心黑洞约900--3000引力半径. 而SDSS和Keck LRIS的Hα双峰轮廓需要用两个盘模型拟合, 对应的吸积盘发射区有两个区域,其中一个发射区(``外盘'')与BOSS类似, 另一个发射区(``内盘'')位于约400--900引力半径范围内, 到中心黑洞的距离显著小于外盘. 结合连续谱的光变特征,推测内盘消失是造成SDSS/Keck LRIS光谱和BOSS光谱在连续谱和发射线轮廓上同时发生光变的主要原因.     

射电宁静类星体紫外/光学光变的混合模型研究Ⅱ:来自SDSS DR7的的数据支持

到目前为止射电宁静类星体的光变起源还不清楚,其中一个可能的机制是吸积盘-再辐射模型。通过对斯隆数字巡天第7次释放数据(Sloan Digital Sky Survey data release 7, SDSS DR7)中789个经过多历元观测并且表现出变亮变蓝现象的射电宁静类星体进行分析,发现黑洞质量较小类星体的综合相对光变谱(残余谱/综合谱)的变亮变蓝幅度更大;其次,对于比较长的光变时标,尤其是在紫外波段,类星体变亮变蓝的幅度也更大。吸积盘-再辐射模型对以上两个主要结果中相对光变谱都拟合得很好,为该模型提供了进一步的支持。     

吸积盘内区向外区的辐射谱

本文计算了吸积盘光学厚的内区向光学薄的外区发射的辐射谱，在中介频谱区辐射谱不同于谱指数为1／3的幂律谱，所得结果对更精细的研究吸积盘的连续谱、发射线以及自辐照模型可能有应用的前景。     

径移主导吸积盘的含时数值模拟

黑洞吸积被认为是活动星系核的能量来源,而研究黑洞吸积最主要的难题是要解释吸积气体的角动量是如何转移出去的.黑洞吸积盘理论提出通过粘滞力矩转移角动量.径移主导吸积盘模型(advection-dominated accretion flow,简称ADAF)是几何厚、光学薄的低吸积率吸积盘模型,它在解释低光度活动星系核的一些观测现象时获得了很大的成功.近年来,随着计算机技术的迅猛发展,数值模拟在天体物理研究中获得了广泛的应用,也取得了突破性的进展.主要工作是用数值模拟方法研究ADAF整体解作为初始条件模拟黑洞的吸积过程.     

黑洞吸积盘理论进展

在已知的四种黑洞吸积盘模型中，Shapiro-Lightman-Eardley(SLE)盘是不稳定的，细（Slim)盘研究得还不够，较为成功的是最早建立的Shakura-Sunyaev盘（SSD）和近年成为热点的径移主导吸积流（ADAF）。SSD和ADAF看来分别适用于吸积流中离黑洞较远和较近的区域，故二者的结合即ADAF＋SSD模型有望对黑洞吸积流作出较为完整的描述，但这个结合模型也还有不少未解决的问题。     

含电子-正电子对的双温吸积盘

本文研究包含正负电子对的双温吸积盘,正负电子对是通过非热的级联过程产生的.结果表明,对于轫致辐射盘和轫致辐射康普顿化盘,都存在一个临界吸积率.当吸积率小于临界吸积率时,盘结构有唯一的稳态解,大于临界吸积率则没有稳态解.我们求得了轫致辐射康普顿化盘的临界吸积率与α和ε_e的关系曲线.对于轫致辐射盘,电子对的影响可以忽略;而对于轫致辐射康普顿化盘,电子对的密度可达10~(-2)～10~3,电子对对盘的结构和光谱都有很大的影响.     

AGN吸积盘辐射区半径的多方法研究

通过收集数据,用标准α盘模型法、短时标光变法、连续波混响滞后法、微引力透镜法分析了活动星系核吸积盘辐射区半径,并对4种方法进行分析讨论。结果表明:(1)将耀变体和类星体参数比较,发现两个子类的中心黑洞质量、光度和吸积盘辐射区半径都没有明显差别。选择高光度源进行研究,出现了两个子类以上参数无明显区别的情况;(2)短时标光变法得到的吸积盘辐射区半径比和黑洞质量没有明显的相关性;(3)连续波混响滞后法得到的吸积盘辐射区半径比对黑洞质量的依赖性高于标准α盘模型法,对于无法确定中心黑洞质量的源,可以利用连续波混响滞后法得到吸积盘辐射区半径;(4)微引力透镜法得到的吸积盘辐射区半径比和黑洞质量没有明显的相关性,且该方法只适用于存在微引力透镜效应的类星体;(5)这些方法从观测上证明了吸积盘辐射区半径与黑洞质量存在相关性,验证了标准模型成立。为进一步在观测上寻找这4种方法所需要的源提供了理论指导,对活动星系核的吸积盘辐射区半径研究有重要意义。     

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.     

Calculation of the global system of the field-aligned currents on the dayside of the magnetosphere

Using the well-known equation for the normal component of the current which exist near the tangential discontinuity in the plasma in the case of the frozen-in magnetic field, and supposing that the current closes in the ionosphere in the auroral oval in the region 1, one calculates and compares with the data of observations the dependence of the density of the field-aligned current at the level of the ionosphere on the local time.     

On the systematic errors of the determinations of the absolute orientation of the meridian marks

The classical method of determination of the absolute azimuth (or Bessel's parameter n) can secure sufficiently precision for RA from observations of stars at high geographical latitudes during polar night only.     

Rotation of the elastic Earth: the role of the angular-velocity-dependence of the elasticity-caused perturbation

We calculate the so-called convective term, which shows up in the expression for the angular velocity of the elastic Earth, within the Andoyer formalism. The term emerges due to the fact that the elasticity-caused perturbation depends not only on the instantaneous orientation of the Earth but also on its instantaneous angular velocity. We demonstrate that this term makes a considerable contribution into the overall angular velocity. At the same time the convective term turns out to be automatically included into the correction to the nutation series due to the elasticity, if the series is defined by the perturbation of the figure axis (and not of the rotational axis) in accordance with the current IAU resolution. Hence it is not necessary to take the effect of the convective term into consideration in the perturbation of the elastic Earth as far as the nutation is related to the motion of the figure axis.     

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the role of the motion of the photospheric footpoints of the magnetic field lines

By means of a simple relation between the velocity v of the fluid particle and the velocity v_f of the photospheric footpoint of the magnetic field line v_z and B_z being respectively the components of v and the magnetic field B normal to the photospheric surface, it is shown formally that through the phtospheric surface the transport of all the quantities attributed to the magnetic field, such as the magnetic flux, the magnetic energy and the helicity, is independent of v_z, and v_f is the only kinematical quantity on which the transport depends. In addition, in the neighborhood of the neutral line the velocity v_l of the moving curve of constant B_z is found to be equal approximately to the component of v or v_f in the direction of v_l. Since v_l can be measured or extimated, so can the components of v and v_f near the neutral line.     

Determination of the temperature of the solar corona from the spectrum of the electron-scattering continuum

When the K-corona is formed by the scattering of photospheric radiation from free electrons, the Fraunhofer lines are greatly broadened by the thermal motions of the hot electrons. This paper discusses the possibility of measuring the coronal electron temperature from the residual depressions in the K-coronal spectrum. If the ratio of the intensities at 4100 Å and 3900 Å can be measured to an accuracy of ±1%, the coronal temperature can be inferred to an accuracy of ±0.2 MK. The temperature of a coronal inhomogeneity may also be measured by this method, provided the position angle is known.Now at Fraunhofer Institute, Freiburg, Germany.