OJ 287的引力波辐射研究

类星体OJ 287的研究已经有100多年的历史,它呈现周期性双峰爆发现象,爆发周期为12 yr,比较好地解释观测的模型是双黑洞模型,即次黑洞绕主黑洞运动,并撞击主黑洞吸积盘从而引起爆发.此模型合理地解释了OJ 287的光变曲线并正确预言了未来爆发时间,这间接证明了广义相对论进动效应以及引力波的存在.星系中心大质量黑洞是一类重要的引力波源,而精确确定了其内部组分运动学方程的星系非常少,由于双黑洞模型提供了精确的黑洞运动轨道,故可以在这个运动轨道基础上研究其引力波辐射.在已有工作基础上,采用后牛顿近似方法首次得到了引力波辐射功率和波形随时间的演化关系,根据目前引力波探测设备IPTA(International Pulsar Timing Array)的进展情况,未来十几年内对OJ 287的引力波直接探测将成为可能.     

估算有中子星作伴星的黑洞质量

在本文中，我们讨论了由中子星和黑洞构成的密近双星系统的引力效应。首先，我们清楚地解释了引力辐射存在引起的轨道圆化，并估算了圆化的时间尺度。求解了围绕史瓦西黑洞作圆轨道运动的矢量进动方程，得到了一个精确解。由进动周期和中心黑洞的质量的关系，我们得到了估算中子星伴随的黑洞质量的方法。随后我们讨论了以中子星脉冲信号随自旋进动变化来寻找黑洞的方法。     

估算有中子星作伴星的黑洞质量

在本文中，我们讨论了由中子星和黑洞构成的密近双星系统的引力效应，首先，我们清楚地解释了引力辐射存在引起的轨道圆化，并估算了圆化的时间尺度，求解了围绕史瓦西黑洞作圆轨道运动的矢量进动方程，得到了一个精确解，由进动周期和中心黑洞的质量的关系，我们得到了估算中子星伴随的黑洞质量的方法，随后我们讨论了以中子星脉冲信号随自旋进动变化来寻找黑洞的方法。     

双星系统中大质量比致密天体质量四极矩对轨道频率的影响分析

极端质量比旋进系统是空间引力波探测器最重要的波源之一。对引力波的探测需要高精度波形模版。当前主流的极端质量比旋进系统引力波计算模型中,人们一般将小质量天体当作试验粒子进行计算,而忽略了其结构及自身引力对背景引力场的影响。利用Mathisson-Papapetrou-Dixon方程研究延展体在弯曲时空中的运动,以及小天体自旋和质量多极矩对引力波信号识别产生的影响。结果表明,质量比在10?6-10?4范围的旋进系统,其自旋达到很大时,自旋对延展体的轨道运动有不可忽略的影响;在质量比10?4-10?2区间内,需要考虑中心黑洞潮汐作用导致的白矮星形变;在质量比大于10?4,且白矮星自旋很大时,其自旋产生的形变会对小天体轨道运动产生不可忽略的影响。大质量黑洞潮汐作用导致的恒星级黑洞或中子星产生的形变可以忽略,中子星和黑洞的自旋会对轨道运动产生不可忽略的影响,而自旋产生的四极矩对轨道运动不产生影响。     

OJ 287,0716+714与BL Lacertae的光变分析

OJ 287是存在着剧烈活动的低峰频BL Lac天体,其低频段的能谱与另两个TeV BL Lac天体(0716+714和BL Lacertae)在低频段的能谱很相似,但是切仑科夫望远镜却没能探测到它的TeV射线.利用这3个天体的观测数据,比较它们在22 GHz、37 GHz和B波段的最小光变周期及延迟的异同,进一步寻找没有观测到OJ 287的TeV伽马射线的可能原因.分析结果显示:(1)最小光变周期方面,OJ 287在37 GHz和B波段的周期偏小,在22 GHz,OJ 287与0716+714的结果相当,但与BL Lacertae相比要小很多,OJ 287的周期更短表明其活动性更强,却没有探测到来自OJ 287的TeV伽马射线,这表明OJ 287在TeV波段的辐射与这3个低能波段最小光变周期之间可能没有联系;(2)延迟方面,OJ 287在B波段相对于37 GHz的延迟要长于0716+714,短于BL Lacertae;在37 GHz相对于22 GHz的延迟要短于0716+714,而BL Lacertae在37 GHz相对于22 GHz的时延为负值,表明22 GHz要超前于37 GHz.通过对延迟的比较分析,并没有发现OJ 287与0716+714和BL Lacertae之间存在明显的差异;从能谱来看,很可能是由于OJ 287在TeV波段的能谱较陡造成切仑科夫望远镜没有探测到来自OJ 287的伽马辐射,但TeV能段较陡的能谱对低能段光变的影响目前还不是很清楚.     

自旋致密双星后牛顿哈密顿系统轨道动力学数值研究

由中子星或黑洞构成的旋转致密双星后牛顿哈密顿系统属于相对论二体问题,该系统不但含有丰富的共振和混沌等动力学现象,而且成为探测引力波的理想天然波源.引力体的轨道动力学性质会在引力波中得到反映.因此,实际天体的混沌性既可能是对引力波探测的挑战,又可望是获得观测效应的机遇.本学位论文正是在这样的国际学术氛围下数值研究旋转致密双星后牛顿保守哈密顿动力学问题.基于最小二乘法原理我们构造了单和双标度因子等几种流形改正方法,分别对旋转致密双星后牛     

R-N场中雷达回波延迟实验的后后牛顿修正

从球对称引力场中光子在赤道面上的运动微分方程出发,采用后牛顿近似方法,讨论了Reissner-Nordstr(o|¨)m度规场中雷达回波延迟实验的后后牛顿修正.通过计算给出了后牛顿以及后后牛顿的修正结果,当荷电量Q为零时,该结果适用于电中性天体,其后牛顿部分的修正与经典文献一致.文中采用的近似展开方法同样适用于讨论更高阶的修正.雷达回波的高阶修正对目前展开的高精度空间引力试验有着重要的参考价值.     

短周期彗星表面水冰升华分布研究

彗星是太阳系遗留的原始星子,研究彗星彗核的演化对理解太阳系其他天体的形成和演化历史具有重要意义.在太阳的辐射作用下,彗星携带的挥发性成分会发生升华,并带动尘埃运动,造成彗核物质的损失.因此,彗核的升华活动对其表面形貌甚至整体形状演化都会产生影响.从IAU (International Astronomical Union) MPC (Minor Planet Center)获取轨道数据,并考虑了彗核的自转以及进动,利用MONET (Mass lossdriven shape evolution model)形状演化模型对短周期彗星做数值模拟,计算得到了短周期彗星1P/Halley、9P/Tempel 1、 19P/Borrelly、 67P/C-G (Churyumov-Gerasimenko)、 81P/Wild 2和103P/Hartley 2在一个轨道周期内的太阳辐射能量以及表面侵蚀深度的分布,结合其动力学参数讨论了自转、进动和公转等特性对其表面水冰升华分布的影响以及造成南北侵蚀差异的可能性.     

有吸积盘时径向运动云的谱线轮廓

在活动星系核的黑洞模型中,黑洞附近运动的电离云所发出的光线将要受到多普勒频移、引力红移和光线偏折等效应的影响,如果黑洞周围还有吸积盘存在,电离云的辐射还可能被吸积盘遮挡。本文全面地考虑了这些效应,利用光子输运方程方法,给出了在Schwarzchild度规中径向运动电离云发出的谱线轮廓的精确解,并发现在某些情况下将出现不对称的双峰结构。     

粒子在标量张量理论的二阶后闵可夫斯基近似弱引力场中的引力偏转

标量张量理论是目前为止最成功的相对论性引力理论之一,经受住了已有实验和观测的检验.随着实验和观测精度的提高,理论将得到进一步的检验.为此将该理论的可观测效应计算到二阶后牛顿近似非常必要.标量张量理论的二阶后牛顿近似结果已经由Xie等人给出.为得到有质量粒子在弱引力场中的偏转情况,利用二阶后闵可夫斯基近似下的度规解得到了粒子的测地线方程,并通过迭代的方法求得在非束缚条件下,粒子在弱引力场中的轨迹的解析解.进一步利用此结果计算了粒子在弱引力场中的偏转角(与光线的偏转角不同)并与前人的结果进行了比较.     

The Orbit Computation of Binary Black Hole OJ287—A 3.5-th Order Post-Newtonian Approximation

The massive binary black hole OJ287 is a source of intense gravita- tional radiation. To detect the signal of its gravitational waves, a knowledge of the signal waveform will be of great help, and this is mainly determined by the orbital motion of the binary. For this, we carry out a detailed calculation on the orbital motion of OJ287, using the post-Newtonian (PN) approximation up to the 3.5th order within the framework of general relativity. Our result is one order higher than the previous work made by others. As in the process of radiation, there is a time delay from the instance when the secondary black hole impacts on the accretion disk of the primary to the moment of the optical outburst. This time delay has to be taken into consideration when we try to fit the calculating orbit with the observed times of outbursts. Adopting a linear relation between the time delay and the impact distance as an empirical model, we fit the cal- culating orbit with the recent 7 outbursts of OJ287, and obtain the solution of its orbital motion, as well as its averaged orbital parameters. By analyzing the result of 3.5 PN order calculation of the binary system, we find some interesting features. In the early period, the rate of precession of the secondary black hole increases, while in the late period approximate to merging, the rate of precession attains its maximum. Afterwards it diminishes, and finally becomes negative. At present we cannot determine whether this behavior is due to the insuffcient accuracy of the 3.5-th order approximation. For the term of dissipative radiation in the equation of motion, the coeffcients of the 2.5 and 3.5 PN orders possess opposite signs. This implies that the 3.5-th order term represents the absorption of energy from outside. However, the sum of the 2.5-th order and 3.5-th order terms still behaves as radiating gravitational waves outward, the rate of energy variation of the system is negative. The calculated result of this work may be useful for more accurate calculations of the gravitational radiation of OJ287.     

A Research on the Gravitational Wave Radiation of OJ 287

The research on quasar OJ 287 has lasted over 100 years. OJ 287 exhibits the phenomenon of periodic two-peak outbursts with the eruptive period of 12 years. Observations are rather well interpreted with the black hole binary model. In this model, the secondary black hole moves around the primary black hole and crashes against the accretion disk of the primary black hole, causing outbursts. This model reasonably explains the light curves of OJ 287 and correctly predicts the time of future outbursts. These indirectly justify the precessional effect of general relativity and the existence of gravitational waves. The massive black hole in the center of galaxy is an important kind of gravitational wave source. However, the number of the galaxies with precisely determined kinematical equations of inner components is quite small. The precise kinematic orbits of black holes are provided by the black hole binary model, so the radiation of gravitational waves can be studied on the basis of these kinematic orbits. Based on the existing work, the evolutionary relations of the radiation power and waveform of gravitational waves with time are first derived by using the post-Newtonian approximation method. According to the current progress of the detection equipment of gravitational waves, i.e., IPTA (International Pulsar Timing Array), the direct detection of gravitational waves from OJ 287 may be possible within the future decade.     

Measuring black hole spin in OJ287

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of ${1.84\cdot 10^{10}M_\odot}$ , the secondary mass ${1.40\cdot 10^{8} M_\odot}$ , major axis precession rate 39°.1 per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is 0.28 ± 0.01 (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the post-Newtonian equations of motion does not exactly follow Einstein’s theory: a parameter q is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein’s value q = 1. This supports the no-hair theorem of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein’s value with the accuracy of 2% (1 sigma). There is a possibility of improving the accuracy of both quantities using the exact timing of the outburst on 21 July 2019. Because of closeness of OJ287 to the Sun (8–12°), the observations would be best carried out by a telescope in space.     

A new model for the periodic outbursts of the BL Lac object OJ287

Recent observations provide strong evidence for the BL Lac object OJ287 exhibiting a 11.6±0.5 yr periodicity with a double-peaked maxima in its optical flux variations. Several models have been proposed for the optical behavior. The 2005 November outburst in OJ287 gives us a surprising result since calculation based on the periodicity was predicting such an outburst in late 2006. Here we suggest a new model, it can not only explain the optical quasi-periodic behavior, but also the radio flares behavior which is and simultaneous with the optical flares. We propose that OJ287 is a binary pair of super-massive black holes, both of them creating a jet. The quasi-periodic double peaks would be due to the relativistic beaming effect on the emission coming from the double helix jets. We used “core flares” to explain the large lags between light curves at different frequencies, and the assumption of two jets appear to be merged with each other partly in the radio frequency emitting regions provides a viable interpretation that we can see only a broad maximum which contains two radio flares that we cannot distinguish.     

Laser Interferometer Space Antenna double black holes: dynamics in gaseous nuclear discs

We study the inspiral of double black holes, with masses in the Laser Interferometer Space Antenna ( LISA ) window of detectability, orbiting inside a massive circumnuclear, rotationally supported gaseous disc. Using high-resolution smoothed particle hydrodynamics simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they corotate with the gaseous disc. As a consequence, the massive black holes bind forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits, the initial eccentricity (if present) does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed (≃1 pc). In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits; eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, the active galactic nucleus activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance scales of tens of parsecs.     

Gravitational signals due to tidal interactions between white dwarfs and black holes

In this paper, we compute the gravitational signal emitted when a white dwarf moves around a black hole on a closed or open orbit using the affine-model approach. We compare the orbital and the tidal contributions to the signal, assuming that the star moves in a safe region where, although very close to the black hole, the strength of the tidal interaction is insufficient to provoke the stellar disruption. We show that for all considered orbits the tidal signal presents sharp peaks corresponding to the excitation of the non-radial oscillation modes of the star, the amplitude of which depends on how deep the star penetrates the black hole tidal radius and on the type of orbit. Further structure is added to the emitted signal by the coupling between the orbital and the tidal motions.     

The problem of few black holes

The problem of few black holes becomes important in multiple mergers of galaxies. If supermassive black holes in centres of galaxies are common, then interaction of three or four supermassive black holes should also be common. The merger of two galaxies with one black hole each produces a semi-stable black hole binary system. Subsequent mergers of galaxies with their own central black holes produces dynamical few-body evolution in which mergers of black holes occur. According to our numerical simulations this evolution typically ends when only one or two black holes remain and, in the latter case, they are ejected in opposite directions from the center of the galaxy. Even when we pick the initial black hole masses at random from a wide distribution, the two black hole ejections happen rather symmetrically. Sometimes the final masses differ considerably in which case only the lighter black hole is ejected. This is caused by the potential barrier of the galaxy itself which prevents the heavy slowly moving black hole flying out of the galaxy. We discuss OJ287 as a possible example of a multiple black hole system.     

Supermassive black-hole binary systems in active galactic nuclei

An analysis of periodic components of flux variability was carried out based on the long-term monitoring of the nuclei of active galaxies 3C454.3, 1633+382, and 3C120, performed in the Crimean Astrophysical Observatory from 1985 to 2008 at 22.2 and 36.8 GHz. Long-period components of the variability (12–14 yrs) were detected and interpreted in terms of the precessional motion of the central body in binary systems. Short-period components (1.5–3 yrs) related to the orbital periods for the motions of the central supermassive black holes were also detected. We concluded that the brightest active galaxies observed as nonstationary sources in a wide range of wavelengths are binary systems of supermassive black holes at the stage of evolution close to coalescence. For the supposed binary black-hole systems, the masses of the central objects and their companions, the orbital radii of the companions, and the coalescence times were determined. The ratios of the masses in the binary systems in all cases proved to be less than ten, pointing to a strong gravitational effect of the companion on the central black hole. The velocities of the central body motion proved to be high, approximately 1000 km/s. This fact should be accounted for in the calculations of the rate of accretion onto the central body. The orbital radii of the companions fall into a narrow range between 4 × 10¹⁶ cm and 6 × 10¹⁶ cm, demonstrating a strong dependence of the masses of the binary systems on the orbital sizes and the energy loss for the gravitational radiation. Within the orbit of the companion during its motion through the accretion disk, a high temperature of surrounding gas is achieved. The high density of the medium, 10⁹–10¹⁰ cm^?3, combined with the magnetic field and shock waves propagating in the accretion disk, develop the conditions for powerful energy release in the directed jets.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.