Supermassive Binary Black Hole System Research Articles

Parameter estimation for inspiraling eccentric compact binaries including pericenter precession

Inspiraling supermassive black hole binary systems with high orbital eccentricity are important sources for space-based gravitational wave observatories like the Laser Interferometer Space Antenna. Eccentricity adds orbital harmonics to the Fourier transform of the gravitational wave signal, and relativistic pericenter precession leads to a three-way splitting of each harmonic peak. We study the parameter estimation accuracy for such waveforms with different initial eccentricity, using the Fisher matrix method and a Monte Carlo sampling of the initial binary orientation. The eccentricity improves the parameter estimation by breaking degeneracies between different parameters. In particular, we find that the source localization precision improves significantly for higher-mass binaries due to eccentricity. The typical sky position errors are $\ensuremath{\sim}1\text{ }\text{ }\mathrm{deg}$ for a nonspinning, ${10}^{7}{M}_{\ensuremath{\bigodot}}$, equal-mass binary at redshift $z=1$, if the initial eccentricity 1 yr before merger is ${e}_{0}\ensuremath{\sim}0.6$. Pericenter precession does not affect the source localization accuracy significantly, but it does further improve the mass and eccentricity estimation accuracy systematically by a factor of 3--10 for masses between ${10}^{6}{M}_{\ensuremath{\bigodot}}$ and ${10}^{7}{M}_{\ensuremath{\bigodot}}$ for ${e}_{0}\ensuremath{\sim}0.3$.

THE FIRST SPECTROSCOPICALLY RESOLVED SUB-PARSEC ORBIT OF A SUPERMASSIVE BINARY BLACK HOLE

One of the most intriguing scenarios proposed to explain how active galactic nuclei are triggered involves the existence of a supermassive binary black hole system in their cores. Here we present an observational evidence for the first spectroscopically resolved sub-parsec orbit of a such system in the core of Seyfert galaxy NGC 4151. Using a method similar to those typically applied for spectroscopic binary stars we obtained radial velocity curves of the supermassive binary system, from which we calculated orbital elements and made estimates about the masses of components. Our analysis shows that periodic variations in the light and radial velocity curves can be accounted for an eccentric, sub-parsec Keplerian orbit of a 15.9-year period. The flux maximum in the lightcurve correspond to the approaching phase of a secondary component towards the observer. According to the obtained results we speculate that the periodic variations in the observed H{\alpha} line shape and flux are due to shock waves generated by the supersonic motion of the components through the surrounding medium. Given the large observational effort needed to reveal this spectroscopically resolved binary orbital motion we suggest that many such systems may exist in similar objects even if they are hard to find. Detecting more of them will provide us with insight into black hole mass growth process.

Parsec-scale jet precession in BL Lacertae (2200+420)

BL Lacertae is the prototype of the BL Lac class of active galactic nuclei, exhibiting intensive activity on parsec (pc) scales, such as intense core variability and multiple ejections of jet components. In particular, in previous works the existence of precession motions in the pc-scale jet of BL Lacertae has been suggested. In this work we revisit this issue, investigating temporal changes of the observed right ascension and declination offsets of the jet knots in terms of our relativistic jet-precession model. The seven free parameters of our precession model were optimized via a heuristic cross-entropy method, comparing the projected precession helix with the positions of the jet components on the plane of the sky and imposing constraints on their maximum and minimum superluminal velocities. Our optimized best model is compatible with a jet having a bulk velocity of 0.9824c, which is precessing with a period of about 12.1 yr in the observer's reference frame and changing its orientation in relation to the line of sight between 4 and 5 degrees, approximately. Assuming that the jet precession has its origin in a supermassive binary black hole system, we show that the 2.3-yr periodic variation in the structural position angle of the very-long-baseline interferometry (VLBI) core of BL Lacertae reported by Stirling et al. is compatible with a nutation phenomenon if the secondary black hole has a mass higher than about six times that of the primary black hole.

Observing the Dynamics of Supermassive Black Hole Binaries with Pulsar Timing Arrays

Pulsar timing arrays are a prime tool to study unexplored astrophysical regimes with gravitational waves. Here, we show that the detection of gravitational radiation from individually resolvable supermassive black hole binary systems can yield direct information about the masses and spins of the black holes, provided that the gravitational-wave-induced timing fluctuations both at the pulsar and at Earth are detected. This in turn provides a map of the nonlinear dynamics of the gravitational field and a new avenue to tackle open problems in astrophysics connected to the formation and evolution of supermassive black holes. We discuss the potential, the challenges, and the limitations of these observations.

Photocentric variability of quasars caused by variations in their inner structure: consequences forGaiameasurements

We study the photocenter position variability due to variations in the quasar inner structure. We consider variability in the accretion disk emissivity and torus structure variability due to different illumination by the central source. We discuss possible detection of these effects by Gaia. Observations of the photocenter variability in two AGNs, SDSS J121855+020002 and SDSS J162011+1724327 have been reported and discussed. With investigation of the variations in the quasar inner structure we explore how much this effect can affect the position determination and whether it can be (or not) detected with Gaia mission. We used (a) a model of a relativistic disk, including the perturbation that can increase brightness of a part of the disk, and consequently offset the photocenter position, and (b) a model of a dusty torus which absorbs and re-emits the incoming radiation from accretion disk. We estimated the value of the photocenter offset due to these two effects. We found that perturbations in the inner structure can significantly offset the photocenter. It depends on the characteristics of perturbation and accretion disk and structure of the torus. In the case of two considered QSOs the observed photocenter offsets cannot be explained by variations in the accretion disk and other effects should be considered. We discussed possibility of exploding stars very close to the AGN source, and also possibility that there are two variable sources in the center of these two AGNs that may indicate a binary super-massive black hole system on a kpc (pc) scale. The Gaia mission seems to be very perspective, not only for astrometry, but also for exploring the inner structure of AGNs. We conclude that variations in the quasar inner structure can affect the observed photocenter (up to several mas). There is a chance to observe such effect in the case of bright and low-redshifted QSOs.

Investigation of some galactic and extragalactic gravitational phenomena

Here we present a short overview of the most important results of our investigations of the following galactic and extragalactic gravitational phenomena: supermassive black holes in centers of galaxies and quasars, supermassive black hole binaries, gravitational lenses and dark matter. For the purpose of these investigations, we developed a model of a relativistic accretion disk around a supermassive black hole, based on the ray-tracing method in the Kerr metric, a model of a bright spot in an accretion disk and three different models of gravitational microlenses. All these models enabled us to study physics, spacetime geometry and effects of strong gravity in the vicinity of supermassive black holes, variability of some active galaxies and quasars, different effects in the lensed quasars with multiple images, as well as the dark matter fraction in the Universe. We also found an observational evidence for the first spectroscopically resolved sub-parsec orbit of a supermassive black hole binary system in the core of active galaxy NGC 4151. Besides, we studied applications of one potential alternative to dark matter in the form of a modified theory of gravity on Galactic scales, to explain the recently observed orbital precession of some S-stars, which are orbiting around a massive black hole at the Galactic center.

Pulsars

AbstractPulsars can be considered as the ultimate time-variable source. They show variations on time-scales ranging from nanoseconds to as long as years, and they emit over almost the entire electromagnetic spectrum. The dominant modulation is associated with the rotation period, which can vary from slighty more than a millisecond to upwards of ten seconds (if we include the magnetars). Variations on time-scales shorter than the pulse period are mostly associated with emission processes and are manifested as giant pulses, microstructure and sub-pulses (to name a few). On time-scales of a rotation to a few hundred rotations are other phenomena also associated with the emission, such as nulling, moding, drifting and intermittency.By probing these and slightly longer time-scales we find that pulsars exhibit “glitches”, which are rapid variations in spin rates. They are believed to be related to the interaction between the superfluid interior of the neutron star and the outer crust. Detailed studies of glitches can reveal much about the properties of the constituents of neutron stars—the only way to probe the physics of material at such extreme densities. Time-scales of about an hour or longer reveal that some pulsars are in binary systems, in particular the most rapidly rotating systems. Discovering and studying those binary systems provides vital clues to the evolution of massive stars, while some of the systems are also the best probes of strong-field gravity theories; the elusive pulsar-black hole binary would be the ultimate system.Pulsars are tools that allow us to probe a range of phenomena and time-scales. It is possible to measure the time of arrival of pulses from some pulsars to better than a few tens of nanoseconds over years, making them some of the most accurate clocks known. Concerning their rotation, deviations from sphericity may cause pulsars to emit gravitational waves which might then be detected by next-generation gravitational-wave detectors. Pulsars themselves can be used as the arms of a Galactic-scale gravitational-wave detector. Measuring correlated deviations in the arrival times of pulses from a number of pulsars distributed throughout the Galaxy could give rise to a direct detection of the stochastic gravitational-wave background, which is associated with the astrophysics of the early Universe—most likely from supermassive black-hole binary systems, but potentially also from cosmic strings. While they are famed for their clock-like rotational stability, some pulsars—in particular the more youthful ones—exhibit modulation in pulse arrival times, often called timing noise. It was recently demonstrated that in at least some cases this variability is deterministic and is associated with modulations in the pulsar emission properties and the spin-down rate. This breakthrough may lead to further improvements in the precision which can be achieved with pulsar timing, and enhance still further the ability to test theories of gravity directly and to make a direct detection of gravitational waves.I presented some of the history of what is known about the variations in pulsars on all these time-scales and reviewed some of the recent achievements in our understanding of the phenomena. I also highlighted how new transients associated with radio-emitting neutron stars are being discovered, and how other transient sources are being identified by the same techniques. These continued improvements have come about without new telescopes, but the next generation of very sensitive wide-field instruments will permit observational cadences which will reveal many new manifestations and will further revolutionise our understanding of this class of objects which have such high astrophysical potential.

Electromagnetic counterparts from counter-rotating relativistic kicked discs

We show the results of two dimensional general relativistic inviscid and isothermal hydrodynamical simulations comparing the behavior of co-rotating (with respect to the black hole rotation) and counter-rotating circumbinary quasi-Keplerian discs in the post merger phase of a supermassive binary black hole system. While confirming the spiral shock generation within the disc due to the combined effects of mass loss and recoil velocity of the black hole, we find that the maximum luminosity of counter-rotating discs is a factor ∼(2–12) higher than in the co-rotating case, depending on the spin of the black hole. On the other hand, the luminosity peak happens ∼10 days later with respect to the co-rotating case. Although the global dynamics of counter-rotating discs in the post merger phase of a merging event is very similar to that for co-rotating discs, an important difference has been found. In fact, increasing the spin of the central black hole produces more luminous co-rotating discs while less luminous counter-rotating ones.

MULTIPLE TIDAL DISRUPTIONS AS AN INDICATOR OF BINARY SUPERMASSIVE BLACK HOLE SYSTEMS

We find that the majority of systems hosting multiple tidal disruptions are likely to contain hard binary SMBH systems, and also show that the rates of these repeated events are high enough to be detected by LSST over its lifetime. Therefore, these multiple tidal disruption events provide a novel method to identify super-massive black hole (SMBH) binary systems with parsec to sub-parsec separations. The rates of tidal disruptions are investigated using simulations of non-interacting stars initially orbiting a primary SMBH and the potential of the model stellar cusp. The stars are then evolved forward in time and perturbed by a secondary SMBH inspiraling from the edge of the cusp to its stalling radius. We find with conservative magnitude estimates that the next generation transient survey LSST should detect multiple tidal disruptions in approximately 3 galaxies over 5 years of observation, though less conservative estimates could increase this rate by an order of magnitude.

THE FIRST KINEMATIC DETERMINATION OF MILLION-YEAR PRECESSION PERIOD OF ACTIVE GALACTIC NUCLEI

Short precession periods like 164d of SS433 can be well determined by observations of time scales longer or much longer than the precession period. However, it doesn't work for sources with precession periods of millions of years. This paper utilizes the particular morphologies of X-shaped sources, so that the 3 dimension kinematics of lobes can be obtained. Thus, for the first time, the million-year precession period of X-shaped sources by observer on the Earth can be determined elegantly: $(6.1\pm 1.5)$Myr, $(1.8\pm 0.5)$Myr, and $(3.2\pm 1.2)$Myr for 3C52, 3C223.1 and 4C12.03 respectively. The result naturally explains the asymmetry displayed in the morphology of these sources, and the effect of propagation time on the diversity of morphologies is well demonstrated. The precession period may originate from long-term effects of a binary super-massive black hole system, which is a potential source of gravitational wave radiation.

Retrograde accretion and merging supermassive black holes

We investigate whether a circumbinary gas disc can coalesce a supermassive black hole binary system in the centre of a galaxy. This is known to be problematic for a prograde disc. We show that in contrast, interaction with a retrograde circumbinary disc is considerably more effective in shrinking the binary because there are no orbital resonances. The binary directly absorbs negative angular momentum from the circumbinary disc by capturing gas into a disc around the secondary black hole, or discs around both holes if the binary mass ratio is close to unity. In many cases the binary orbit becomes eccentric, shortening the pericentre distance as the eccentricity grows. In all cases the binary coalesces once it has absorbed the angular momentum of a gas mass comparable to that of the secondary black hole. Importantly, this conclusion is unaffected even if the gas inflow rate through the disc is formally super--Eddington for either hole. The coalescence timescale is therefore always $\sim M_2/\dot M$, where $M_2$ is the secondary black hole mass and $\dot M$ the inflow rate through the circumbinary disc.

Long-term variation time scales in OJ 287

The light curve data from 1894 to 2008 are compiled for the BL Lacertae object OJ 287 from the available literature. Periodicity analysis methods (the Discrete Correlation Function-DCF, the Jurkevich method, the power spectral (Fourier) analysis, and the CLEANest method) are performed to search for possible periodicites in the light curve of OJ 287. Significance levels are given for the possible periods. The analysis results confirm the existence of the 12.2±0.6 yr time scale and show a hint of a ∼53 yr time scale. The 12.2±0.6 yr period is used as the orbital period to investigate the supermassive binary black hole system parameters.

Short-term VHE variability in blazars: PKS 2155-304

Context: The $\gamma$-ray blazar PKS 2155-304 has attracted considerable attention because of its extreme TeV variability characteristics during an exceptional flaring period in 2006. Among the observed key findings are (i) a minimum variability timescale as short as $\sim 200$ sec and (ii) highly variable TeV emission, which in the frequency interval [$10^{-4}$ Hz, $10^{-2}$ Hz] can be described by a log-normal distribution and suggests an underlying multiplicative (and not additive) process. Aims: Simultaneously accounting for these findings appears difficult within conventional approaches. Following earlier suggestions for the TeV blazar Mkn 501, we explore a possible scenario where PKS 2155-304 is supposed to harbor a supermassive binary black hole system and where the observed TeV variability is dominated by emission from the less massive black hole. Methods: We analyze the constraints on the very high energy (VHE) source imposed by the observed variability characteristics and the integrated VHE luminosity output, and discuss its implications for a binary black hole system. Results: We show that for a secondary mass of $m_{\rm BH} \sim 10^7 M_{\odot}$, fluctuations in the disk accretion rate that feed the jet could account for the observed red-noise type variability process down to frequencies of $\sim 10^{-2}$ Hz. Jet curvature induced by orbital motion, on the other hand, could further relax constraints on the intrinsic jet speeds. Conclusions: Because a binary system can lead to different (yet not independent) periodicities in different energy bands, a longterm (quasi-) periodicity analysis could offer important insights into the real nature of the central engine of PKS~2155-304.

The sensitivity of the Parkes Pulsar Timing Array to individual sources of gravitational waves

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

Delensing gravitational wave standard sirens with shear and flexion maps

Supermassive black hole binary systems (SMBHB) are standard sirens -- the gravitational wave analogue of standard candles -- and if discovered by gravitational wave detectors, they could be used as precise distance indicators. Unfortunately, gravitational lensing will randomly magnify SMBHB signals, seriously degrading any distance measurements. Using a weak lensing map of the SMBHB line of sight, we can estimate its magnification and thereby remove some uncertainty in its distance, a procedure we call "delensing." We find that delensing is significantly improved when galaxy shears are combined with flexion measurements, which reduce small-scale noise in reconstructed magnification maps. Under a Gaussian approximation, we estimate that delensing with a 2D mosaic image from an Extremely Large Telescope (ELT) could reduce distance errors by about 30-40% for a SMBHB at z=2. Including an additional wide shear map from a space survey telescope could reduce distance errors by 50%. Such improvement would make SMBHBs considerably more valuable as cosmological distance probes or as a fully independent check on existing probes.

Elliptic waveforms for inspiralling compact binaries

The inspiral of supermassive black hole binary systems with high orbital eccentricity are the most promising sources for the gravitational wave observatories. The importance of elliptic gravitational waveforms in various physical scenarios has been emphasized by several authors (Wahlquist 1987, Moreno-Garrido, Buitrago and Mediavilla 1994, Martel and Poisson 1999). Taking into account the eccentricity of the orbit in the total waveform improves the parameter estimation for these sources, as it is shown by the construction and analyzation of the Fisher information matrix. In our work we use the Fourier-Bessel analysis of the Kepler motion and the stationary phase approximation of time-depend waveforms.

THE QUASAR SDSS J1536+0441: AN UNUSUAL DOUBLE-PEAKED EMITTER

The quasar SDSS J153636.22+044127.0, exhibiting peculiar broad emission-line profiles with multiple components, was proposed as a candidate sub-parsec binary supermassive black hole system. More recently, imaging revealed two spatially distinct sources, leading some to suggest the system to be a quasar pair separated by ~5 kpc. We present Palomar and Keck optical spectra of this system from which we identify a third velocity component to the emission lines. We argue that the system is more likely an unusual member of the class of active galactic nuclei (AGNs) known as "double-peaked emitters" than a sub-parsec black hole binary or quasar pair. We find no significant velocity evolution of the two main peaks over the course of 0.95 yr, with a 3-sigma upper limit on any secular change of 70 km/s/yr. We also find that the three velocity components of the emission lines are spatially coincident to within 0.015" along the slit, apparently ruling out the double-quasar hypothesis.

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 × 1016 cm and 6 × 1016 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, 109–1010 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.

H I OBSERVATIONS OF THE SUPERMASSIVE BINARY BLACK HOLE SYSTEM IN 0402+379

We have recently discovered a supermassive binary black hole system with a projected separation between the two black holes of 7.3 parsecs in the radio galaxy 0402+379. This is the most compact supermassive binary black hole pair yet imaged by more than two orders of magnitude. We present Global VLBI observations at 1.3464 GHz of this radio galaxy, taken to improve the quality of the HI data. Two absorption lines are found toward the southern jet of the source, one redshifted by 370 +/- 10 km/s and the other blueshifted by 700 +/- 10 km/s with respect to the systemic velocity of the source, which, along with the results obtained for the opacity distribution over the source, suggests the presence of two mass clumps rotating around the central region of the source. We propose a model consisting of a geometrically thick disk, of which we only see a couple of clumps, that reproduces the velocities measured from the HI absorption profiles. These clumps rotate in circular Keplerian orbits around an axis that crosses one of the supermassive black holes of the binary system in 0402+379. We find an upper limit for the inclination angle of the twin jets of the source to the line of sight of 66 degrees, which, according to the proposed model, implies a lower limit on the central mass of ~7 x 10^8 Msun and a lower limit for the scale height of the thick disk of ~12 pc .

A candidate sub-parsec supermassive binary black hole system

The role of mergers in producing galaxies, together with the finding that most large galaxies harbour black holes in their nuclei, implies that binary supermassive black hole systems should be common. Here we report that the quasar SDSS J153636.22+044127.0 is a plausible example of such a system. This quasar shows two broad-line emission systems, separated in velocity by 3,500 km s(-1). A third system of unresolved absorption lines has an intermediate velocity. These characteristics are unique among known quasars. We interpret this object as a binary system of two black holes, having masses of 10(7.3) and 10(8.9) solar masses separated by approximately 0.1 parsec with an orbital period of approximately 100 years.