OJ287双黑洞轨道计算:3.5阶后牛顿近似

大质量双黑洞OJ287是一个强引力辐射源.为了探测其引力波信号,需要知道波形,而这主要是由轨道运动所决定.为此,从广义相对论3.5阶后牛顿近似的运动方程出发对OJ287的轨道进行仔细研究,取大黑洞位置固定作为近似,给出了后牛顿近似下3.5阶的次黑洞轨道解,比他人2.5阶的工作高了一阶.次黑洞撞击吸积盘面到光学爆发存在时间延迟,这对于确定轨道参数有很大影响.利用径向距离与爆发时间关系的线性模型,对最近7次爆发时刻的观测值拟合,给出了更精确的OJ287双黑洞的轨道参数及其运动轨道.分析了计算结果,研究运动特征,并且发现了两个新性质:次黑洞进动在初期增加,在晚期接近并和时,进动达到最大值,然后减小并越过0而趋于负值.尚不能确定晚期的这个行为是否由3.5阶近似不够准确所造成.运动方程中耗散性的辐射项,后牛顿2.5阶和3.5阶的系数具有相反的符号.这意味着3.5阶项反而是从外界吸收能量.但2.5阶与3.5阶之和仍然是向外辐射引力波的,体系能量变化率为负.这个工作的计算结果可以用来更精确地计算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.     

Analyses of the Light Variations of OJ 287, 0716+714 and BL Lacertae

OJ 287 is a BL Lac object which exhibits intense activities of low peak-frequencies. Its energy spectrum in low frequency band is quite similar with those of two other TeV BL Lac objects (i.e., 0716+714 and BL Lacertae). However, the Cerenkov telescope did not detect its TeV rays. By using the observational data of these three heavenly bodies and comparing the discrepan- cies of their minimal periods of light variations and delays at 22 GHz, 37 GHz and B-waveband, we have further investigated the possible reason why the TeV gamma-rays of OJ 287 have not been observed. The results of analyses are as fol- lows. (1) For the minimal periods of light variations, the periods of OJ 287 at 37 GHZ and B-waveband are short. At 22 GHz the results of OJ 287 and 0716+714 are comparable, but the period of OJ 287 is much shorter in comparison with that of BL Lacertae, and this shows that its activity is more intense. However, the TeV gamma-rays of OJ 287 have not been detected, which implies that the radiation of OJ 287 in TeV waveband may have no connection with the minimal periods of light variations in these three low-energy wavebands. (2) In respect of delays, the delay of OJ 287 in the B waveband with respect to 37 GHz is longer than that of 0716+714, but shorter than that of BL Lacertae. Its delay at 37 GHz in respect to 22 GHz is shorter than that of 0716+714. Meanwhile, the delay of BL Lacertae at 37 GHz in respect to 22 GHz is negative, which implies that 22 GHz precedes 37 GHz. Via the comparison and analysis of delays, no obvious differences between OJ 287 and 0716+714 as well as BL Lacertae have been found. On the side of energy spectra, it is quite possible that due to the steep energy spectrum of OJ 287 in TeV waveband, the Cerenkov telescope has not detected the gamma radiation of OJ 287. However, nowadays it is still not clear whether the steep energy spectrum in TeV energy range has some influence on the light variations in low energy realm.     

Alignment time-scale of the microquasar GRO J1655−40

The microquasar GRO J1655−40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady-state disc warped by Lense–Thirring precession and find the time-scale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically 50 times higher in the latter case but we find that, in both the cases, the lifetime of the mass-transfer state is at most a few times the alignment time-scale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter aligned after its supernova kick but that it is most likely to be close to alignment rather than counter alignment now.     

Some secondary indications of gravitational collapse

A stellar core becomes somewhat less massive due to neutrinos radiated away during its collapse in a neutron star or a black hole. The paper deals with the hydrodynamic motion of stellar envelope induced by such a mass loss. Depending on the structure of the outer stellar layers, the motion results either in ejection of an envelope with mass and energy proper for Nova outbursts; or nearly instantaneous excitation of strong pulsations of the star; or lastly in a slow slipping away of the whole stellar envelope. These phenomena are of importance when more powerful events, like supernova outbursts presumably associated with gravitational collapse, are absent. Such secondary indications of gravitational collapse are of special interest, since they may be a single observable manifestation (besides neutrinos and gravitational waves) of massive black hole formation.     

Gravitational loss-cone instability in stellar systems with retrograde orbit precession

We study spherical and disc clusters in a near-Keplerian potential of galactic centres or massive black holes. In such a potential orbit precession is commonly retrograde, that is, the direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly-radial orbits. We show that if there is a loss-cone at low angular momentum (e.g. due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics   l ≥ 3  . If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independent of azimuthal number m . The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum.
The analysis is based on simple characteristic equations describing small perturbations in a disc or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analysing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose–Nyquist criterion for disc and spherical gravitational systems.     

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.     

Precession of the Orbit of a Planet around Stars Revolving along: Low-Eccentricity Orbits in a Binary System: Analytical Solution

In our previous paper (hereafter, paper I) we presented analytical results on the non-planar motion of a planet around a binary star for the cases of the circular orbits of the components of the binary. We found that the orbital plane of the planet (the plane containing the “unperturbed” elliptical orbit of the planet), in addition to precessing about the angular momentum of the binary, undergoes simultaneously the precession within the orbital plane. We demonstrated that the analytically calculated frequency of this additional precession is not the same as the frequency of the precession of the orbital plane about the angular momentum of the binary, though the frequencies of both precessions are of the same order of magnitude. In the present paper we extend the analytical results from paper I by relaxing the assumption that the binary is circular – by allowing for a relatively small eccentricity ε of the stars orbits in the binary. We obtain an additional, ε-dependent term in the effective potential for the motion of the planet. By analytical calculations we demonstrate that in the particular case of the planar geometry (where the planetary orbit is in the plane of the stars orbits), it leads to an additional contribution to the frequency of the precession of the planetary orbit. We show that this additional, ε-dependent contribution to the precession frequency of the planetary orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Besides, we also obtain analytical results for another type of the non-planar configuration corresponding to the linear oscillatory motion of the planet along the axis of the symmetry of the circular orbits of the stars. We show that as the absolute value of the energy increases, the period of the oscillations decreases.     

Gravitational and electromagnetic emission by magnetized coalescing binary systems

We discuss the possibility to obtain an electromagnetic emission accompanying the gravitational waves emitted in the coalescence of a compact binary system. Motivated by the existence of black hole configurations with open magnetic field lines along the rotation axis, we consider a magnetic dipole in the system, the evolution of which leads to (i) electromagnetic radiation, and (ii) a contribution to the gravitational radiation, the luminosity of both being evaluated. Starting from the observations on magnetars, we impose upper limits for both the electromagnetic emission and the contribution of the magnetic dipole to the gravitational wave emission. Adopting this model for the evolution of neutron star binaries leading to short gamma ray bursts, we compare the correction originated by the electromagnetic field to the gravitational waves emission, finding that they are comparable for particular values of the magnetic field and of the orbital radius of the binary system. Finally we calculate the electromagnetic and gravitational wave energy outputs which result comparable for some values of magnetic field and radius.     

Orbital evolution measurement of the accreting millisecond X-ray pulsar SAX J1808.4-3658

We present results from a pulse timing analysis of the accretion-powered millisecond X-ray pulsar SAX J1808.4-3658 using X-ray data obtained during four outbursts of this source. Extensive observations were made with the proportional counter array of the Rossi X-ray Timing Explorer (RXTE) during the four outbursts that occurred in 1998, 2000, 2002 and 2005. Instead of measuring the arrival times of individual pulses or the pulse arrival time delay measurement that is commonly used to determine the orbital parameters of binary pulsars, we have determined the orbital ephemeris during each observation by optimizing the pulse detection against a range of trial ephemeris values. The source exhibits a significant pulse shape variability during the outbursts. The technique used by us does not depend on the pulse profile evolution, and is therefore, different from the standard pulse timing analysis. Using 27 measurements of orbital ephemerides during the four outbursts spread over more than 7 years and more than 31,000 binary orbits, we have derived an accurate value of the orbital period of 7249.156862(5) s (MJD = 50915) and detected an orbital period derivative of (3.14 ± 0.21) × 10⁻¹² s s⁻¹. We have included a table of the 27 mid-eclipse time measurements of this source that will be valuable for further studies of the orbital evolution of the source, especially with ASTROSAT. We point out that the measured rate of orbital period evolution is considerably faster than the most commonly discussed mechanisms of orbital period evolution like mass transfer, mass loss from the companion star and gravitational wave radiation. The present time scale of orbital period change, 73 Myr is therefore likely to be a transient high value of period evolution and similar measurements during subsequent outbursts of SAX J1808.4-3658 will help us to resolve this.     

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.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

Supermassive binary black holes among cosmic gamma-ray sources

Supermassive binary black holes (SBBHs) are a natural outcome of galaxy mergers. Here we show that low-frequency (f≤10⁻⁶ Hz) quasi-periodic variability observed from cosmic blazar sources can provide substantial inductive support for the presence of close (d≲0.1 pc) SBBHs at their centers. It is argued on physical grounds that such close binary systems are likely to give rise to different (although not independent) periodicities in the radio, optical and X-ray/TeV regime, and, hence that detection of appropriate period ratios significantly corroborates the SBBH interpretation. This is illustrated for a binary model where optical longterm periodicity is related to accretion disk interactions, radio periodicity to Newtonian jet precession, and periodicities in the high energy bands to the orbital motion of the jet. We use the observed periodicities to constrain the properties for a sample of SBBH candidates including OJ 287 and AO 0235+16, and discuss the results within the context of jet activity and binary evolution.     

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能段较陡的能谱对低能段光变的影响目前还不是很清楚.     

Analytical solution for the three-dimensional motion of a circumbinary planet around a binary star

By using the method of separating rapid and slow subsystem, we obtain an analytical solution for a stable three-dimensional motion of a circumbinary planet around a binary star. We show that the motion of the planet is more complicated than it was obtained for this situation analytically by Farago and Laskar (2010). Namely, in addition to the precession of the orbital plane of the planet around the angular momentum of the binary (found by Farago and Laskar (2010)), there is simultaneously the precession of the orbital plane of the planet within the orbital plane. We show that the frequency of this additional precession is different from the frequency of the precession of the orbital plane around the angular momentum of the binary. We demonstrate that this problem is mathematically equivalent both to the problem of the motion of a satellite around an oblate planet and to the problem of a hydrogen Rydberg atom in the field of a high-frequency linearly-polarized laser radiation, thus discovering yet another connection between astrophysics and atomic physics. We point out that all of the above physical systems have a higher than geometrical symmetry, which is a counterintuitive result. In particular, it is manifested by the fact that, while the elliptical orbit of the circumbinary planet (around a binary star) or of the satellite (around an oblate planet) or of the Rydberg electron (in the laser field) undergoes simultaneously two types of the precession, the shape of the orbit does not change. The fact that a system, consisting of a circumbinary planet around a binary star, possesses the hidden symmetry should be of a general physical interest. Our analytical results could be used for benchmarking future simulations.     

Perihelion precession in binary systems: higher order corrections

Higher order corrections (up to n-th order) are obtained for the perihelion precession in binary systems like OJ287 using the Schwarzschild metric and complex integration. The corrections are performed considering the third root of the motion equation and developing the expansion in terms of \(r_{s}/ (a(1-e^{2}) )\).The results are compared with other expansions that appear in the literature giving corrections to second and third order. Finally, we simulate the shape of relativistic orbits for binary systems with different masses.