HI 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 .Comment: 20 pages, 7 figures. Accepted on the Astrophysical Journa

Accretion vs colliding wind models for the gamma-ray binary LS I +61 303: an assessment

LS I +61 303 is a puzzling Be/X-ray binary with variable gamma-ray emission at up TeV energies. The nature of the compact object and the origin of the high-energy emission are unclear. One family of models invokes particle acceleration in shocks from the collision between the B-star wind and a relativistic pulsar wind, while another centers on a relativistic jet powered by accretion. Recent high-resolution radio observations showing a putative "cometary tail" pointing away from the Be star near periastron have been cited as support for the pulsar-wind model. We wish here to carry out a quantitative assessment of these competing models for this extraordinary source. We apply a 3D SPH code for dynamical simulations of both the pulsar-wind-interaction and accretion-jet models. The former yields a description of the shape of the wind-wind interaction surface. The latter provides an estimation of the accretion rate. The results allow critical evaluation of how the two distinct models confront the data in various wavebands under a range of conditions. When one accounts for the 3D dynamical wind interaction under realistic constraints for the relative strength of the B-star and pulsar winds, the resulting form of the interaction front does not match the putative "cometary tail" claimed from radio observations. On the other hand, dynamical simulations of the accretion-jet model indicate that the orbital phase variation of accretion power includes a secondary broad peak well away from periastron, thus providing a plausible way to explain the observed TeV gamma ray emission toward apastron. We conclude that the colliding-wind model is not clearly established for LS I +61 303, while the accretion-jet model can reproduce many key characteristics of the observed TeV gamma-ray emission.Comment: Accepted for publication in A&A. The resolution of the figures is lower than in the journal paper to minimize file sizes. Seven pages, 5 figure

Gas driven massive black hole binaries: signatures in the nHz gravitational wave background

Pulsar timing arrays (PTAs) measure nHz frequency gravitational waves (GWs) generated by orbiting massive black hole binaries (MBHBs) with periods between 0.1-10 yr. Previous studies on the nHz GW background assumed that the inspiral is purely driven by GWs. However, torques generated by a gaseous disk can shrink the binary much more efficiently than GW emission, reducing the number of binaries at these separations. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of MBHBs through physically distinct regions of the disk, until GWs take over their evolution. We extract MBHB cosmological merger rates from the Millennium simulation, generate Monte Carlo realizations of a population of gas driven binaries, and calculate the corresponding GW amplitudes of the most luminous individual binaries and the stochastic GW background. For steady state alpha-disks with alpha>0.1 we find that the nHz GW background can be significantly modified. The number of resolvable binaries is however not changed by the presence of gas; we predict 1-10 individually resolvable sources to stand above the noise for a 1-50 ns timing precision. Gas driven migration reduces predominantly the number of small total mass or unequal mass ratio binaries, which leads to the attenuation of the mean stochastic GW--background, but increases the detection significance of individually resolvable binaries. The results are sensitive to the model of binary--disk interaction. The GW background is not attenuated significantly for time-dependent models of Ivanov, Papaloizou, & Polnarev (1999).Comment: Accepted by MNRAS, 15 pages, 8 figure

Properties of Accretion Flows Around Coalescing Supermassive Black Holes

What are the properties of accretion flows in the vicinity of coalescing supermassive black holes (SBHs)? The answer to this question has direct implications for the feasibility of coincident detections of electromagnetic (EM) and gravitational wave (GW) signals from coalescences. Such detections are considered to be the next observational grand challenge that will enable testing general relativity in the strong, nonlinear regime and improve our understanding of evolution and growth of these massive compact objects. In this paper we review the properties of the environment of coalescing binaries in the context of the circumbinary disk and hot, radiatively inefficient accretion flow models and use them to mark the extent of the parameter space spanned by this problem. We report the results from an ongoing, general relativistic, hydrodynamical study of the inspiral and merger of black holes, motivated by the latter scenario. We find that correlated EM+GW oscillations can arise during the inspiral phase followed by the gradual rise and subsequent drop-off in the light curve at the time of coalescence. While there are indications that the latter EM signature is a more robust one, a detection of either signal coincidentally with GWs would be a convincing evidence for an impending SBH binary coalescence. The observability of an EM counterpart in the hot accretion flow scenario depends on the details of a model. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, upper limits on luminosity imply that they may be identified by EM searches out to z~0.1-1. However, given the radiatively inefficient nature of the gas flow, we speculate that a majority of massive binaries may appear as low luminosity AGN in the local universe.Comment: Revised version accepted to Class. Quantum Grav. for proceedings of 8th LISA Symposium. 15 pages, 3 figures, includes changes suggested in referee report

Massive binary black holes in galactic nuclei and their path to coalescence

Massive binary black holes form at the centre of galaxies that experience a merger episode. They are expected to coalesce into a larger black hole, following the emission of gravitational waves. Coalescing massive binary black holes are among the loudest sources of gravitational waves in the Universe, and the detection of these events is at the frontier of contemporary astrophysics. Understanding the black hole binary formation path and dynamics in galaxy mergers is therefore mandatory. A key question poses: during a merger, will the black holes descend over time on closer orbits, form a Keplerian binary and coalesce shortly after? Here we review progress on the fate of black holes in both major and minor mergers of galaxies, either gas-free or gas-rich, in smooth and clumpy circum-nuclear discs after a galactic merger, and in circum-binary discs present on the smallest scales inside the relic nucleus.Comment: Accepted for publication in Space Science Reviews. To appear in hard cover in the Space Sciences Series of ISSI "The Physics of Accretion onto Black Holes" (Springer Publisher

Modeling high-energy light curves of the PSR B1259-63/LS 2883 binary based on 3-D SPH simulations

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.Comment: 32 pages, 6 figues. Accepted for publication in Ap

The pre-outburst flare of the A 0535+26 August/September 2005 outburst

We study the spectral and temporal behavior of the High Mass X-ray Binary A 0535+26 during a `pre-outburst flare' which took place ~5 d before the peak of a normal (type I) outburst in August/September 2005. We compare the studied behavior with that observed during the outburst. We analyse RXTE observations that monitored A 0535+26 during the outburst. We complete spectral and timing analyses of the data. We study the evolution of the pulse period, present energy-dependent pulse profiles both at the initial pre-outburst flare and close to outburst maximum, and measure how the cyclotron resonance-scattering feature (hereafter CRSF) evolves. We present three main results: a constant period P=103.3960(5)s is measured until periastron passage, followed by a spin-up with a decreasing period derivative of Pdot=(-1.69+/-0.04)x10^(-8)s/s at MJD 53618, and P remains constant again at the end of the main outburst. The spin-up provides evidence for the existence of an accretion disk during the normal outburst. We measure a CRSF energy of Ecyc~50kev during the pre-outburst flare, and Ecyc~46kev during the main outburst. The pulse shape, which varies significantly during both pre-outburst flare and main outburst, evolves strongly with photon energy.Comment: 4 pages, 4 figures, accepted for publication in A&A Letters. To be published in parallel to Postnov et al. 200

Multimessenger astronomy with pulsar timing and X-ray observations of massive black hole binaries

We demonstrate that very massive (>10^8\msun), cosmologically nearby (z<1) black hole binaries (MBHBs), which are primary targets for ongoing and upcoming pulsar timing arrays (PTAs), are particularly appealing multimessenger carriers. According to current models for massive black hole formation and evolution, the planned Square Kilometer Array (SKA) will collect gravitational wave signals from thousands of such massive systems, being able to individually resolve and locate in the sky several of them (maybe up to a hundred). By employing a standard model for the evolution of MBHBs in circumbinary discs, with the aid of dedicated numerical simulations, we characterize the gas-binary interplay, identifying possible electromagnetic signatures of the PTA sources. We concentrate our investigation on two particularly promising scenarios in the high energy domain, namely, the detection of X-ray periodic variability and of double broad K\alpha iron lines. Up to several hundreds of periodic X-ray sources with a flux >10^-13 erg s^-1 cm^-2 will be in the reach of upcoming X-ray observatories. Double relativistic K\alpha lines may be observable in a handful of low redshift (z<0.3) sources by proposed deep X-ray probes, such as Athena. (Abridged)Comment: 19 pages, 11 figures, submitted to MNRAS, minor revision of the reference lis

Viscous propagation of mass flow variability in accretion discs

We study mass flow rate through a disc resulting from a varying mass supply rate. Variable mass supply rate occurs, e.g., during disc state transitions, and in interacting eccentric binaries. It is, however, damped by the viscosity of the disc. Here, we calculate this damping in detail. We derive an analytical description of the propagation of the flow rate using the solution of Lynden-Bell & Pringle, in which the disc is assumed to extend to infinity. In particular, we derive the accretion-rate Green's function, and its Fourier transform, which gives the fractional damping at a given variability frequency. We then compare this model to that of a finite disc with the mass supply at its outer edge. We find significant differences with respect to the infinite disc solution, which we find to overestimate the viscous damping. In particular, the asymptotic form of the Green's function is power-law for the infinite disc and exponential for the finite one. We then find a simple fitting form for the latter, and also calculate its Fourier transform. In general, the damping becomes very strong when the viscous time at the outer edge of the disc becomes longer than the modulation time scale. We apply our results to a number of astrophysical systems. We find the effect is much stronger in low-mass X-ray binaries, where the disc size is comparable to that of the Roche lobe, than in high-mass binaries, where the wind-fed disc can have a much smaller size.Comment: MNRAS, in pres

Quasi-periodic flares in EXO 2030+375 observed with INTEGRAL

Context: Episodic flaring activity is a common feature of X-ray pulsars in HMXBs. In some Be/X-ray binaries flares were observed in quiescence or prior to outbursts. EXO 2030+375 is a Be/X-ray binary showing "normal" outbursts almost every ~46 days, near periastron passage of the orbital revolution. Some of these outbursts were occasionally monitored with the INTEGRAL observatory. Aims: The INTEGRAL data revealed strong quasi-periodic flaring activity during the rising part of one of the system's outburst. Such activity has previously been observed in EXO 2030+375 only once, in 1985 with EXOSAT. (Some indications of single flares have also been observed with other satellites.) Methods: We present the analysis of the flaring behavior of the source based on INTEGRAL data and compare it with the flares observed in EXO 2030+375 in 1985. Results: Based on the observational properties of the flares, we argue that the instability at the inner edge of the accretion disk is the most probable cause of the flaring activity.Comment: Accepted for publication in A&A Lette