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.     

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的引力辐射.     

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.     

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 waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays

Massive black holes are key components of the assembly and evolution of cosmic structures, and a number of surveys are currently on going or planned to probe the demographics of these objects and to gain insight into the relevant physical processes. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses  ≳10⁷ M_⊙  . The whole cosmic population produces a stochastic background that could be detectable with upcoming PTAs. Sources sufficiently close and/or massive generate gravitational radiation that significantly exceeds the level of the background and could be individually resolved. We consider a wide range of massive black hole binary assembly scenarios, investigate the distribution of the main physical parameters of the sources, such as masses and redshift, and explore the consequences for PTAs observations. Depending on the specific massive black hole population model, we estimate that on average at least one resolvable source produces timing residuals in the range  ∼5–50 ns  . PTAs, and in particular the future Square Kilometre Array, can plausibly detect these unique systems, although the events are likely to be rare. These observations would naturally complement on the high-mass end of the massive black hole distribution function future surveys carried out by the Laser Interferometer Space Antenna .     

tempo2: a new pulsar timing package – III. Gravitational wave simulation

Analysis of pulsar timing data sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.     

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.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

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.     

Dark matter accretion wakes of high-redshift black holes

Anisotropic emission of gravitational waves during the merger of black holes induces a recoil velocity on the centre of mass of the binary and the final merger product can then be ejected from its host galaxy. We consider ejected black holes which stay on bound orbits around their host haloes. A recoiled black hole which moves on an almost radial orbit outside the virial radius of its central galaxy, in the cold dark matter background, reaches its apapsis in a finite time. Due to small dark matter velocity dispersion at high redshifts and also the small black hole velocity near the apapsis passage a high-density wake forms around these black hole. Gamma-ray emission can result from the enhancement of dark matter annihilation in these wakes. The diffuse high-energy gamma-ray background from the ensemble of such black holes in the Hubble volume is also evaluated.     

An axion-like scalar field environment effect on binary black hole merger

The environment, such as an accretion disk, could modify the signal of the gravitational wave from astrophysical black hole binaries. In this article, we model the matter field around intermediatemass binary black holes by means of an axion-like scalar field and investigate their joint evolution. In detail, we consider equal mass binary black holes surrounded by a shell of axion-like scalar field both in spherically symmetric and non-spherically symmetric cases, and with different strengths of the scalar field. Our result shows that the environmental scalar field could essentially modify the dynamics. Firstly,in the spherically symmetric case, with increase of the scalar field strength, the number of circular orbits for the binary black hole is reduced. This means that the scalar field could significantly accelerate the merger process. Secondly, once the scalar field strength exceeds a certain critical value, the scalar field could collapse into a third black hole with its mass being larger than that of the binary. Consequently,the new black hole that collapses from the environmental scalar field could accrete the binary promptly and the binary collides head-on with each other. In this process, there is almost no quadrupole signal produced, and, consequently, the gravitational wave is greatly suppressed. Thirdly, when the scalar field strength is relatively smaller than the critical value, the black hole orbit could develop eccentricity through accretion of the scalar field. Fourthly, during the initial stage of the inspiral, the gravitational attractive force from the axion-like scalar field could induce a sudden turn in the binary orbits, hence resulting in a transient wiggle in the gravitational waveform. Finally, in the non-spherical case, the scalar field could gravitationally attract the binary moving toward the center of mass for the scalar field and slow down the merger process.     

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.     

Variability of Blazars

Variability is one of the characteristics of blazars. The rapid variability is superposed on the long term variation. In this work, the variability on different time scales, such as intra-day (IDV), short-term (STV) and long-term (LTV) variations are presented for some sources. We also presented our own observations of some selected objects, for which the historical data were compiled for periodicity analysis using several methods. The parameters of the binary black hole system OJ 287 are determined.     

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.     

A hybrid N-body code incorporating algorithmic regularization and post-Newtonian forces

We describe a novel N -body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite- N gravitational fluctuations on the orbits of the S-stars, and inspiral of an intermediate-mass black hole into the Galactic Centre.     

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.     

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.     

The first gravitational-wave burst GW150914, as predicted by the scenario machine

The Advanced LIGO observatory recently reported (Abbott et al., 2016a) the first direct detection of gravitational waves predicted by Einstein (1916). The detection of this event was predicted in 1997 on the basis of the Scenario Machine population synthesis calculations (Lipunov et al., 1997b) Now we discuss the parameters of binary black holes and event rates predicted by different scenarios of binary evolution. We give a simple explanation of the big difference between detected black hole masses and the mean black hole masses observed in of X-ray Nova systems. The proximity of the masses of the components of GW150914 is in good agreement with the observed initial mass ratio distribution in massive binary systems, as is used in Scenario Machine calculations for massive binaries.     

Gravitational waves: The future of black hole physics

The new millennium will witness the operation of several long-baseline ground-based interferometric detectors, possibly a space-based detector too, which will make it possible to directly observe black holes by catching gravitational waves emitted by them during their formation or when they are perturbed or when a binary consisting of black holes in-spirals due to radiation reaction. Such observations will help us not only to test some of the fundamental predictions of Einstein's general relativity but will also give us the unique opportunity to map black hole spacetimes, to measure the masses and spins of black holes and their population, etc.