by Hayasaki, Kimitake and Yagi, Kent and Tanaka, Takahiro and Mineshige, Shin
10 pages, 2 figures, submitted to ApJ

When binary black holes are embedded in a gaseous environment, a rotating disk surrounding them, the so-called circumbinary disk, will be formed. The binary exerts a gravitational torque on the circumbinary disk and thereby the orbital angular momentum is transferred to it, while the angular momentum of the circumbinary disk is transferred to the binary through the mass accretion. The binary undergoes an orbital decay due to both the gravitational wave emission and the binary-disk interaction. This causes the phase evolution of the gravitational wave signal. The precise measurement of the gravitational wave phase thus may provide information regarding the circumbinary disk. In this paper, we assess the detectability of the signature of the binary-disk interaction using the future space-borne gravitational wave detectors such as DECIGO and BBO by the standard matched filtering analysis. We find that the effect of the circumbinary disk around binary black holes in the mass range $latex 6M_sun\le{M}\lesssim3\times10^3M_sun$ is detectable at a statistically significant level in five year observation, provided that gas accretes onto the binary at a rate greater than $latex \dot{M}\sim1.4\times10^{17} [gs^{-1}] j^{-1}(M/10M_sun)^{33/23}$ with 10% mass-to-energy conversion efficiency, where j represents the efficiency of the angular momentum transfer from the binary to the circumbinary disk. We show that $latex O(0.1)$ coalescence events are expected to occur in sufficiently dense molecular clouds in five year observation. We also point out that the circumbinary disk is detectable, even if its mass at around the inner edge is by over 10 orders of magnitude less than the binary mass.

by Blecha, Laura and Loeb, Abraham and Narayan, Ramesh
20 pages, 9 figures. Submitted to MNRAS

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN are most likely to be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)

by O’Shaughnessy, R. and Healy, J. and London, L. and Meeks, Z. and Shoemaker, D.
Submitted to PRD

The gravitational wave signature emitted from a merging binary depends on the orientation of an observer relative to the binary. Previous studies suggest that emission along the total initial or total final angular momenta leads to both the strongest and simplest signal from a precessing compact binary. In this paper we describe a concrete counterexample: a binary with $latex m_1/m_2=4$, $latex a_1=0.6 \hat{x} = -a_2$, placed in orbit in the x,y plane. We extract the gravitational wave emission along several proposed emission directions, including the initial (Newtonian) orbital angular momentum; the final (~ initial) total angular momentum; and the dominant principal axis of $latex _M$. Using several diagnostics, we show that the suggested preferred directions are not representative. For example, only for a handful of other directions (0.95). We conclude that the information available in just one direction (or mode) does not adequately encode the complexity of orientation-dependent emission for even short signals from merging black hole binaries. Future investigations of precessing, unequal-mass binaries should carefully explore and model their orientation-dependent emission.

by Agudo, Ivan and Marscher, Alan P. and Jorstad, Svetlana G. and Gomez, Jose L. and Perucho, Manel and Piner, B. Glenn and Rioja, Maria and Dodson, Richard
Accepted for Publication in The Astrophysical Journal. 11 pages, 6 figures, 4 tables. High resolution images on figure 1 and complete tables 1 and 2 may be provided on request

We present the results from an ultra-high-resolution 7mm Very Long Baseline Array (VLBA) study of the relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. Analysis of the image sequence reveals a sharp jet-position-angle swing by >100 deg. during [2004,2006], as viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our images also show for the first time a prominent erratic wobbling behavior of the innermost ~0.4mas of the jet with fluctuations in position angle of up to ~40 deg. over time scales ~2yr. This is accompanied by highly superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as possible origins of the phenomenon on the scales probed by our observations, although such processes may cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet.

by Alsing, Justin and Berti, Emanuele and Will, Clifford and Zaglauer, Helmut
19 pages, 2 figures, 2 tables

We derive the equations of motion, the periastron shift, and the gravitational radiation damping for quasicircular compact binaries in a massive variant of the Brans-Dicke theory of gravity. We also study the Shapiro time delay and the Nordtvedt effect in this theory. By comparing with recent observational data, we put bounds on the two parameters of the theory: the Brans-Dicke coupling parameter \omega_{BD} and the scalar mass m_s. We find that the most stringent bounds come from Cassini measurements of the Shapiro time delay in the Solar System, that yield a lower bound \omega_{BD}>40000 for scalar masses m_s1000 for m_s1250 for m_s<10^{-20} eV. A first estimate suggests that bounds comparable to the Shapiro time delay may come from observations of radiation damping in the eccentric white dwarf-neutron star binary PSR J1141-6545, but a quantitative prediction requires the extension of our work to eccentric orbits.

by Liu, F. K. and Wang, Dong and Chen, Xian
47 pages, 9 figures, 1 table; accepted for publication in ApJ

Numerical relativity simulations predict that coalescence of supermassive black hole (SMBH) binaries not only leads to a spin flip but also to a recoiling of the merger remnant SMBHs. In the literature, X-shaped radio sources are popularly suggested to be candidates for SMBH mergers with spin flip of jet-ejecting SMBHs. Here we investigate the spectral and spatial observational signatures of the recoiling SMBHs in radio sources undergoing black hole spin flip. Our results show that SMBHs in most spin-flip radio sources have mass ratio $latex q\ga 0.3$ with a minimum possible value $latex q_{\rm min} \simeq 0.05$. For major mergers, the remnant SMBHs can get a kick velocity as high as $latex 2100 km s^{-1}$ in the direction within an angle $latex \la 40^\circ$ relative to the spin axes of remnant SMBHs, implying that recoiling quasars are biased to be with high Doppler-shifted broad emission lines while recoiling radio galaxies are biased to large apparent spatial off-center displacements. We also calculate the distribution functions of line-of-sight velocity and apparent spatial off-center for spin-flip radio sources with different apparent jet reorientation angles. Our results show that the larger the apparent jet reorientation angle is, the larger the Doppler-shifting recoiling velocity and apparent spatial off-center displacement will be. We investigate the effects of recoiling velocity on the dust torus in spin-flip radio sources and suggest that recoiling of SMBHs would lead to “dust poor” AGNs. Finally, we collect a sample of 19 X-shaped radio objects and for each object give the probability of detecting the predicted signatures of recoiling SMBH.

by Lahav, Carmit Gordon and Meiron, Yohai and Soker, Noam
10 pages, 3 figures, 2 tables

We examine the scatter in the correlation between black hole masses and their host bulge masses, and find that the increase of the scatter with mass suggests that mergers alone cannot produce the observed correlation. A simple merger scenario of small galaxies leads to a proportionality relation between the late-time black hole and bulge masses, with intrinsic scatter increasing along the ridge line of the relation as the square root of the mass. By examining a sample of 87 galaxies, we find that the scatter increases with mass more rapidly than expected from the merger-only explanation. We do not dispute that mergers play a role, but our results favour a universal feedback mechanism that works on all scales. We discuss the possibility that the feedback mechanism that operated during galaxy formation involved the presence of a cooling flow.

by Narayan, Ramesh and McClintock, Jeffrey E.
Published online in MNRAS, November 21, 2011

We show that the 5-GHz radio flux of transient ballistic jets in black hole binaries correlates with the dimensionless black hole spin parameter a* estimated via the continuum-fitting method. The data suggest that jet power scales either as the square of a* or the square of the angular velocity of the horizon. This is the first direct evidence that jets may be powered by black hole spin energy. The observed correlation validates the continuum-fitting method of measuring spin. In addition, for those black holes that have well-sampled radio observations of ballistic jets, the correlation may be used to obtain rough estimates of their spins.

by Marx, Jay and Danzmann, Karsten and Hough, James and Kuroda, Kazuaki and McClelland, David and Mours, Benoit and Phinney, Sterl and Rowan, Sheila and Sathyaprakash, B. and Vetrano, Flavio and Vitale, Stefano and Whitcomb, Stan and Will, Clifford
116 pages. Original document in higher resolution can be found at https://gwic.ligo.org/roadmap/

Gravitational wave science is on the verge of direct observation of the waves predicted by Einstein’s General Theory of Relativity and opening the exciting new field of gravitational wave astronomy. In the coming decades, ultra-sensitive arrays of ground-based instruments and complementary spaced-based instruments will observe the gravitational wave sky, inevitably discovering entirely unexpected phenomena while providing new insight into many of the most profound astrophysical phenomena known. in July 2007 the Gravitational Wave International Committee (GWIC) initiated the development of a strategic roadmap for the field of gravitational wave science with a 30-year horizon. The goal of this roadmap is to serve the international gravitational wave community and its stakeholders as a tool for the development of capabilities and facilities needed to address the exciting scientific opportunities on the intermediate and long-term horizons.

by Barausse, Enrico and Buonanno, Alessandra and Tiec, Alexandre Le
11 pages, 2 figures

Using the main result of a companion paper, in which the binding energy of a circular-orbit non-spinning compact binary system is computed at leading-order beyond the test-particle approximation, the exact expression of the effective-one-body (EOB) metric component $latex g^\text{eff}_{tt}$ is obtained through first order in the mass ratio. Combining these results with the recent gravitational self-force calculation of the periastron advance for circular orbits in the Schwarzschild geometry, the EOB metric component $latex g^\text{eff}_{rr}$ is also determined at linear order in the mass ratio. These results assume that the mapping between the real and effective Hamiltonians at the second and third post-Newtonian (PN) orders holds at all PN orders. Our findings also confirm the advantage of resumming the PN dynamics around the test-particle limit if the goal is to obtain a flexible model that can smoothly connect the test-mass and equal-mass limits.

by Tiec, Alexandre Le and Barausse, Enrico and Buonanno, Alessandra
5 pages, 1 figure

Using the first law of binary black-hole mechanics, we compute the binding energy E and total angular momentum J of two non-spinning compact objects moving on circular orbits with frequency Omega, at leading order beyond the test-particle approximation. By minimizing E(Omega) we recover the exact frequency shift of the Schwarzschild innermost stable circular orbit induced by the conservative piece of the gravitational self-force. Comparing our results for the coordinate invariant relation E(J) to those recently obtained from numerical simulations of comparable-mass non-spinning black-hole binaries, we find a remarkably good agreement, even in the strong-field regime. Our findings confirm that the domain of validity of perturbative calculations may extend well beyond the extreme mass-ratio limit.

by Tiec, Alexandre Le and Blanchet, Luc and Whiting, Bernard F.
43 pages, 3 figures

First laws of black hole mechanics, or thermodynamics, come in a variety of different forms. In this paper, from a purely post-Newtonian (PN) analysis, we obtain a first law for binary systems of point masses moving along an exactly circular orbit. Our calculation is valid through 3PN order and includes, in addition, the contributions of logarithmic terms at 4PN and 5PN orders. This first law of binary point-particle mechanics is then derived from first principles in general relativity, and analogies are drawn with the single and binary black hole cases. Some consequences of the first law are explored for PN spacetimes. As one such consequence, a simple relation between the PN binding energy of the binary system and Detweiler’s redshift observable is established. Through it, we are able to determine with high precision the numerical values of some previously unknown high order PN coefficients in the circular-orbit binding energy. Finally, we propose new gauge invariant notions for the energy and angular momentum of a particle in a binary system.

by Comerford, Julia M. and Gerke, Brian F. and Stern, Daniel and Cooper, Michael C. and Weiner, Benjamin J. and Newman, Jeffrey A. and Harrison, Fiona and Madsen, Kristin and Barrows, R. Scott
11 pages, 8 figures, submitted to ApJ

Merger-remnant galaxies with kpc-scale separation dual active galactic nuclei (AGNs) should be widespread as a consequence of galaxy mergers and triggered gas accretion onto supermassive black holes, yet very few dual AGNs have been observed. Galaxies with double-peaked narrow AGN emission lines in the Sloan Digital Sky Survey are plausible dual AGN candidates, but their double-peaked profiles could also be the result of gas kinematics or AGN-driven outflows and jets on small or large scales. To help distinguish between these scenarios, we have obtained spatial profiles of the AGN emission via follow-up longslit spectroscopy of 81 double-peaked narrow-line AGNs in SDSS at 0.03 < z < 0.36 using Lick, Palomar, and MMT Observatories. We find that all 81 systems exhibit double AGN emission components with ~kpc projected spatial separations on the sky, which suggests that they are produced by kpc-scale dual AGNs or kpc-scale outflows, jets, or rotating gaseous disks. In addition, we find that the subsample (58%) of the objects with spatially-compact emission components may be preferentially produced by dual AGNs, while the subsample (42%) with spatially-extended emission components may be preferentially produced by AGN outflows. We also find that for 32% of the sample the two AGN emission components are preferentially aligned with the host galaxy major axis, as expected for dual AGNs orbiting in the host galaxy potential. Our results both narrow the list of possible physical mechanisms producing the double AGN components, and suggest several observational criteria for selecting the most promising dual AGN candidates from the full sample of double-peaked narrow-line AGNs. Using these criteria, we determine the 17 most compelling dual AGN candidates in our sample.

by Barausse, Enrico and Buonanno, Alessandra and Hughes, Scott A. and Khanna, Gaurav and O’Sullivan, Stephen and Pan, Yi
19 pages, 14 figures, 6 tables

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries. Although in this analysis we only consider equatorial inspirals (orbital angular momentum parallel to spin), there is no issue of principle preventing us from considering inclined binaries. We will extend this analysis to examine misaligned spin-orbit configurations in future work.

by Aranha, R. F. and Soares, I. Damião and Tonini, E. V.
18 pages, 12 Figures

We examine numerically the post-merger regime of two Schwarzschild black holes in non head-on collision. Our treatment is made in the realm of non-axisymmetric Robinson-Trautman spacetimes which are appropriate for the description of the system. Characteristic initial data for the system are constructed and the Robinson-Trautman equation is integrated using a numerical code based on the Galerkin spectral method. The collision is planar, restricted to the plane determined by the directions of the two initial colliding black holes, with the net momentum fluxes of gravitational waves confined to this plane. We evaluate the efficiency of mass-energy extraction, the total energy and momentum carried out by gravitational waves and the momentum distribution of the remnant black hole. Our analysis is based on the Bondi-Sachs four momentum conservation laws. Head-on collisions and orthogonal collisions constitute, respectively, upper and lower bounds to the power emission and to the efficiency of mass-energy extraction by gravitational waves. The momentum extraction and the pattern of the momentum fluxes, as a function of the incidence angle, are examined. The momentum extraction characterizes a regime of strong deceleration of the system. The angular pattern of gravitational wave signals is also examined. They are typically bremsstrahlung for early times emission. Gravitational waves are also emitted outside the plane of collision but this component has a zero net momentum flux. The relation between the incidence angle of collision and the exit angle of the remnant closely approximates a relation for inelastic collisions of classical particles in Newtonian dynamics.

by Dotti, M. and Sesana, A. and Decarli, R.
4 Figures. Accepted for publication in Advances in Astronomy

The study of the dynamical evolution of massive black hole pairs in mergers is crucial in the context of a hierarchical galaxy formation scenario. The timescales for the formation and the coalescence of black hole binaries are still poorly constrained, resulting in large uncertainties in the expected rate of massive black hole binaries detectable in the electromagnetic and gravitational wave spectra. Here we review the current theoretical understanding of the black hole pairing in galaxy mergers, with a particular attention to recent developments and open issues. We conclude with a review of the expected observational signatures of massive binaries, and of the candidates discussed in literature to date.

by Sesana, A.
21 pages, 4 figures, accepted for publication in Advances in Astronomy

I review our current understanding of massive black hole (MBH) formation and evolution along the cosmic history. After a brief introductory overview of the relevance of MBHs in the hierarchical structure formation paradigm, I discuss the main viable channels for seed BH formation at high redshift and for their subsequent mass growth and spin evolution. The emerging hierarchical picture, where MBHs grow through merger triggered accretion episodes, acquiring their mass while shining as quasars, is overall robust, but too simplistic to explain the diversity observed in MBH phenomenology. I briefly discuss which future observations will help to shed light on the MBH cosmic history in the near future, paying particular attention to the upcoming gravitational wave window.

by Yagi, Kent and Stein, Leo C. and Yunes, Nicolas and Tanaka, Takahiro
26 pages, 3 figures, 2 tables; submitted to PRD

We consider a general class of quantum gravity-inspired, modified gravity theories, where the Einstein-Hilbert action is extended through the addition of all terms quadratic in the curvature tensor coupled to scalar fields with standard kinetic energy. This class of theories includes Einstein-Dilaton-Gauss-Bonnet and Chern-Simons modified gravity as special cases. We analytically derive and solve the coupled field equations in the post-Newtonian approximation, assuming a comparable-mass, spinning black hole binary source in a quasi-circular, weak-field/slow-motion orbit. We find that a naive subtraction of divergent piece associated with the point-particle approximation is ill-suited to represent compact objects in these theories. Instead, we model them by appropriate effective sources built so that known strong-field solutions are reproduced in the far-field limit. In doing so, we prove that black holes in Einstein-Dilaton-Gauss-Bonnet and Chern-Simons theory can have hair, while neutron stars have no scalar monopole charge, in diametrical opposition to results in scalar-tensor theories. We then employ techniques similar to the direct integration of the relaxed Einstein equations to obtain analytic expressions for the scalar field, metric perturbation, and the associated gravitational wave luminosity measured at infinity. We find that scalar field emission mainly dominates the energy flux budget, sourcing electric-type (even-parity) dipole scalar radiation and magnetic-type (odd-parity) quadrupole scalar radiation, correcting the General Relativistic prediction at relative -1PN and 2PN orders. Such modifications lead to corrections in the emitted gravitational waves that can be mapped to the parameterized post-Einsteinian framework. Such modifications could be strongly constrained with gravitational wave observations.

by Valtonen, Mauri J. and Villforth, Carolin and Wiik, Kaj
to appear in Monthly Notices of Royal Astronomical Society

We study the variation of the optical polarization angle in the blazar OJ287 and compare it with the precessing binary black hole model with a ‘live’ accretion disk. First, a model of the variation of the jet direction is calculated, and the main parameters of the model are fixed by the long term optical brightness evolution. Then this model is compared with the variation of the parsec scale radio jet position angle in the sky. Finally, the variation of the polarization angle is calculated using the same model, but using a magnetic field configuration which is at a constant angle relative to the optical jet. It is found that the model fits the data reasonably well if the field is almost parallel to the jet axis. This may imply a steady magnetic field geometry, such as a large-scale helical field.

by Damour, Thibault and Nagar, Alessandro and Pollney, Denis and Reisswig, Christian
4 pages, 2 figures

Using accurate numerical relativity simulations of (nonspinning) black-hole binaries with mass ratios 1:1, 2:1 and 3:1 we compute the gauge invariant relation between the (reduced) binding energy $latex E$ and the (reduced) angular momentum $latex j$ of the system. We show that the relation $latex E(j)$ is an accurate diagnostic of the dynamics of a black-hole binary in a highly relativistic regime. By comparing the numerical-relativity $latex E^{\rm NR} (j)$ curve with the predictions of several analytic approximation schemes, we find that, while the usual, non-resummed post-Newtonian-expanded $latex E^{\rm PN} (j)$ relation exhibits large and growing deviations from $latex E^{\rm NR} (j)$, the prediction of the effective one-body formalism, based purely on known analytical results (without any calibration to numerical relativity), agrees strikingly well with the numerical-relativity results.

by Hergt, Steven and Steinhoff, Jan and Schaefer, Gerhard
37 pages, submitted

The present paper addresses open questions regarding the handling of the spin supplementary condition within the effective field theory approach to the post-Newtonian approximation. In particular it is shown how the spin supplementary condition can be eliminated at the level of the potential and how the dynamics can be cast into a fully reduced Hamiltonian form. Two different methods are used and compared, one based on the well-known Dirac bracket and the other based on an action principle. It is discussed how the latter approach can be used to improve the Feynman rules by formulating them in terms of reduced canonical spin variables.

by Lovelace, Geoffrey and Boyle, Michael and Scheel, Mark A. and Szilagyi, Bela
17 pages, 7 figures, submitted to Classical and Quantum Gravity

Motivated by the possibility of observing gravitational waves from merging black holes whose spins are nearly extremal (i.e., 1 in dimensionless units), we present numerical waveforms from simulations of merging black holes with the highest spins simulated to date: (1) a 25.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.97 aligned with the orbital angular momentum; and (2) a previously reported 12.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.95 anti-aligned with the orbital angular momentum. First, we consider the horizon mass and spin evolution of the new aligned-spin simulation. During the inspiral, the horizon area and spin evolve in remarkably close agreement with Alvi’s analytic predictions, and the remnant hole’s final spin agrees reasonably well with several analytic predictions. We also find that the total energy emitted by a real astrophysical system with these parameters—almost all of which is radiated during the time included in this simulation—would be 10.952% of the initial mass at infinite separation. Second, we consider the gravitational waveforms for both simulations. After estimating their uncertainties, we compare the waveforms to several post-Newtonian approximants, finding significant disagreement well before merger, although the phase of the TaylorT4 approximant happens to agree remarkably well with the numerical prediction in the aligned-spin case. We find that the post-Newtonian waveforms have sufficient uncertainty that hybridized waveforms will require far longer numerical simulations (in the absence of improved post-Newtonian waveforms) for accurate parameter estimation of low-mass binary systems.

by Li, Nan and Mao, Shude and Gao, Liang and Loeb, Abraham and Di Stefano, R.
9 pages, 5 figures

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $latex 10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $latex 10^{-3}$ is $latex \sim 0.2 \fBH$ ($latex \sim 0.05 \fBH$) in a multiply-imaged lens system, where $latex \fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

by Gualandris, Alessia and Dotti, Massimo and Sesana, Alberto
5 pages, 5 figures, submitted to MNRAS

We study the evolution of the orientation of the orbital plane of massive black hole binaries (BHBs) in rotating stellar systems in which the total angular momentum of the stellar cusp is misaligned with respect to that of the binary. We compare results from direct summation N-body simulations with predictions from a simple theoretical model. We find that the same encounters between cusp stars and the BHB that are responsible for the hardening and eccentricity evolution of the binary, lead to a reorientation of the binary orbital plane. In particular, binaries whose angular momentum is initially misaligned with respect to that of the stellar cusp tend to realign their orbital planes with the angular momentum of the cusp on a timescale of a few hardening times. This is due to angular momentum exchange between stars and the BHB during close encounters, and may have important implications for the relative orientation of host galaxies and radio jets.

by Barrows, R. Scott and Stern, Daniel and Madsen, Kristin and Harrison, Fiona and Assef, Roberto J. and Comerford, Julia M. and Cushing, Michael C. and Fassnacht, Christopher D. and Gonzalez, Anthony and Griffith, Roger and Hickox, Ryan and Kirkpatrick, J. Davy and Lagattuta, David J.
18 pages and 9 figures. Accepted for publication in The Astrophysical Journal

The X-ray source CXOXBJ142607.6+353351 (CXOJ1426+35), which was identified in a 172 ks Chandra image in the Bootes field, shows double-peaked rest-frame optical/UV emission lines, separated by 0.69″ (5.5 kpc) in the spatial dimension and by 690 km s^-1 in the velocity dimension. The high excitation lines and emission line ratios indicate both systems are ionized by an AGN continuum, and the double-peaked profile resembles that of candidate dual AGN. At a redshift of z=1.175, this source is the highest redshift candidate dual AGN yet identified. However, many sources have similar emission line profiles for which other interpretations are favored. We have analyzed the substantial archival data available in this field, as well as acquired near-infrared (NIR) adaptive optics (AO) imaging and NIR slit spectroscopy. The X-ray spectrum is hard, implying a column density of several 10^23 cm^-2. Though heavily obscured, the source is also one of the brightest in the field, with an absorption-corrected 2-10 keV luminosity of ~10^45 erg s^-1. Outflows driven by an accretion disk may produce the double-peaked lines if the central engine accretes near the Eddington limit. However, we may be seeing the narrow line regions of two AGN following a galactic merger. While the AO image reveals only a single source, a second AGN would easily be obscured by the significant extinction inferred from the X-ray data. Understanding the physical processes producing the complex emission line profiles seen in CXOJ1426+35 and related sources is important for interpreting the growing population of dual AGN candidates.

by Moesta, Philipp and Alic, Daniela and Rezzolla, Luciano and Zanotti, Olindo and Palenzuela, Carlos
4 pages, 3 figures

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.

by Shields, G. A. and Rosario, D. J. and Junkkarinen, V. and Chapman, S. C. and Bonning, E. W. and Chiba, T.
14 pages, 5 figures. Submitted to The Astrophysical Journal

We present HST and UKIRT spectra and images of the 2 kpc binary quasar LBQS 0103-2753 (z=0.858). The HST images (V- and I-band) show tidal features demonstrating that this system is a major galaxy merger in progress. A two-color composite image brings out knots of star formation along the tidal arc and elsewhere. The infrared spectrum shows that both objects are at the same redshift, and that the discrepant redshift of C IV in component A is a consequence of the BAL absorption in the spectrum of this component. LBQS 0103-2753 is one of the most closely spaced binary QSOs known, and is one of relatively few dual AGN showing confirmed broad emission lines from both components. While statistical studies of binary QSOs suggest that simultaneous fueling of both black holes during a merger may be relatively rare, LBQS 0103-2753 demonstrates that such fueling can occur at high luminosity at a late stage in the merger at nuclear spacing of only a few kpc, without severe obscuration of the nuclei.

by O’Shaughnessy, R. and Kaplan, D. L. and Sesana, A. and Kamble, A.
Submitted to ApJ. v1 original submission; v2 minor changes in response to referee

In the last few years before merger, supermassive black hole binaries will rapidly inspiral and precess in a magnetic field imposed by a surrounding circumbinary disk. Multiple simulations suggest this relative motion will convert some of the local energy to a Poynting-dominated outflow, with a luminosity 10^{43} erg/s * (B/10^4 G)^2(M/10^8 Msun)^2 (v/0.4 c)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. On top of a secular increase in power on the gravitational wave inspiral timescale, orbital motion will produce significant, detectable modulations, both on orbital periods and (if black hole spins are not aligned with the binary’s total angular momenta) spin-orbit precession timescales. Because the gravitational wave merger time increases rapidly with separation, we find vast numbers of these transients are ubiquitously predicted, unless explicitly ruled out (by low efficiency $latex \epsilon$) or obscured (by accretion geometry f_{geo}). If the fraction of Poynting flux converted to radio emission times the fraction of lines of sight accessible $latex f_{geo}$ is sufficiently large (f_{geo} \epsilon > 2\times 10^{-4} for a 1 year orbital period), at least one event is accessible to future blind surveys at a nominal 10^4 {deg}^2 with 0.5 mJy sensitivity. Our procedure generalizes to other flux-limited surveys designed to investigate EM signatures associated with many modulations produced by merging SMBH binaries.

by Popovic, Luka C.
The work was presented as an invited talk at special workshop “Spectral lines and super-massive black holes” held on June 10, 2011 as a part of activity within the frame of COST action 0905 “Black holes in a violent universe” and as a part of the 8th Serbian Conference on Spectral Line Shapes in Astrophysics.Sent to New Astronomy Review as a review paper

The broad emission spectral lines emitted from AGNs are our main probe of the geometry and physics of the broad line region (BLR) close to the SMBH. There is a group of AGNs that emits very broad and complex line profiles, showing two displaced peaks, one blueshifted and one redshifted from the systemic velocity defined by the narrow lines, or a single such peak. It has been proposed that such line shapes could indicate a supermassive binary black hole (SMB) system. We discuss here how the presence of an SMB will affect the BLRs of AGNs and what the observational consequences might be.

We review previous claims of SMBs based on broad line profiles and find that they may have non-SMB explanations as a consequence of a complex BLR structure. Because of these effects it is very hard to put limits on the number of SMBs from broad line profiles. It is still possible, however, that unusual broad line profiles in combination with other observational effects (line ratios, quasi-periodical oscillations, spectropolarimetry, etc.) could be used for SMBs detection.

Some narrow lines (e.g., [O\,III]) in some AGNs show a double-peaked profile. Such profiles can be caused by streams in the Narrow Line Region (NLR), but may also indicate the presence of a kilo-parsec scale mergers. A few objects indicated as double-peaked narrow line emitters are confirmed as kpc-scale margers, but double-peaked narrow line profiles are mostly caused by the complex NLR geometry.

We briefly discuss the expected line profile of broad Fe K$latex \alpha$ that probably originated in the accretion disk(s) around SMBs.

Finally we consider rare configurations where a SMB system might be gravitationally lensed by a foreground galaxy, and discuss the expected line profiles in these systems.

by Fabbiano, G. and Junfeng, Wang and Elvis, M. and Risaliti, G.
Preprint (not final) version of a paper to appear in Nature

The current picture of galaxy evolution advocates co-evolution of galaxies and their nuclear massive black holes (MBHs), through accretion and merging. Quasar pairs (6,000-300,000 light-years separation) exemplify the first stages of this gravitational interaction. The final stages, through binary MBHs and final collapse with gravitational wave emission, are consistent with the sub-light-year separation MBHs inferred from optical spectra and light-variability of two quasars. The double active nuclei of few nearby galaxies with disrupted morphology and intense star formation (e.g., NGC 6240 and Mkn 463; ~2,400 and ~12,000 light-years separation respectively) demonstrate the importance of major mergers of equal mass spirals in this evolution, leading to an elliptical galaxy, as in the case of the double radio nucleus (~15 light-years separation) elliptical 0402+379. Minor mergers of galaxies with a smaller companion should be a more common occurrence, evolving into spiral galaxies with active MBH pairs, but have hitherto not been seen. Here we report the presence of two active MBHs, separated by ~430 light-years, in the Seyfert galaxy NGC 3393. The regular spiral morphology and predominantly old circum-nuclear stellar population of this galaxy, and the closeness of the MBHs embedded in the bulge, suggest the result of minor merger evolution.

by Valtonen, M. J. and Mikkola, S. and Lehto, H. J. and Gopakumar, A. and Hudec, R. and Polednikova, J.
14 pages, 14 figures

We examine the ability to test the black hole no-hair theorem at the 10% level in this decade using the binary black hole in OJ287. In the test we constrain the value of the dimensionless parameter q that relates the scaled quadrupole moment and spin of the primary black hole: q2 = -q 2 . At the present we can say that q = 1 \pm 0.3 (one), in agreement with General Relativity and the no-hair theorems. We demonstrate that this result can be improved if more observational data is found in historical plate archives for the 1959 and 1971 outbursts. We also show that the predicted 2015 and 2019 outbursts will be crucial in improving the accuracy of the test. Space-based photometry is required in 2019 July due the proximity of OJ287 to the Sun at the time of the outburst. The best situation would be to carry out the photometry far from the Earth, from quite a different vantage point, in order to avoid the influence of the nearby Sun. We have considered in particular the STEREO space mission which would be ideal if it has a continuation in 2019 or LORRI on board the New Horizons mission to Pluto.

by Galaviz, Pablo
15 pages, 15 figures and 4 tables

We study the stability and chaos of three compact objects using post-Newtonian (PN) equations of motion derived from the Arnowitt-Deser-Misner-Hamiltonian formulation, where we include terms up to 2.5 PN order in the orbital part and the leading order in spin corrections. We perform numerical simulations of a hierarchical configuration of three compact bodies in which a binary system is perturbed by a third, lighter body initially far away from the binary. The relative importance of the different PN orders is examined. The basin boundary method and the computation of Lyapunov exponent are employed to analyze the stability and chaotic properties of the system. The 1 PN terms produce a small but noticeable change in the stability regions of the parameters considered. On the other hand, the inclusion of spin or gravitational radiation does not produce a significant change with respect to the inclusion of the 1 PN terms.

by Huwyler, Cédric and Klein, Antoine and Jetzer, Philippe
26 pages, 28 figures

In this paper, we compute the accuracy at which the planned space-based gravitational wave detector LISA will be able to observe deviations from General Relativity. To do so, we introduce six correction parameters that account for modified gravity in the second post-Newtonian gravitational wave phase for inspiralling supermassive black hole binaries with spin precession on quasi-circular orbits. The precession of the spins and the angular momentum modulate the gravitational waveform, resulting in additional structure which could reduce correlations in the parameter space and increase the detection accuracy of the alternative theory parameters. Also, the use of higher harmonics could create further structure and increase the time during which the signal lasts in the frequency window of LISA. In order to find error distributions for the alternative theory parameters, we use the Fisher information formalism and carry out Monte Carlo simulations for 17 different binary black hole mass configurations in the range 10^5 Msun < M < 10^8 Msun with 10^3 randomly distributed points in the parameter space each, using the full (FWF) and restricted (RWF) version of the gravitational waveform. We find that the binaries can roughly be separated into two groups: one with low (\precsim 10^7 Msun) and one with high total masses (\succsim 10^7 Msun). The RWF errors on the alternative theory parameters are two orders of magnitude higher than the FWF errors for high-mass binaries while almost comparable for low-mass binaries. Due to dilution of the available information, the accuracy of the binary parameters is reduced by factors of a few, except for the luminosity distance which is affected more seriously in the high-mass regime. As an application, we compute an optimal lower bound on the graviton mass which is increased by a factor of ~1.5 when using the FWF.

by Zanotti, Olindo
6 pages, 4 figures. Accepted by New Astronomy

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, for a binary with a total mass M~10^6 M_\odot. 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.

by Ingram, Adam and Done, Chris
11 pages, 6 figures, 2 tables

Black hole accretion flows show rapid X-ray variability. The Power Spectral Density (PSD) of this is typically fit by a phenomenological model of multiple Lorentzians for both the broad band noise and Quasi-Periodic Oscillations (QPOs). Our previous paper (Ingram & Done 2011) developed the first physical model for the PSD and fit this to observational data. This was based on the same truncated disc/hot inner flow geometry which can explain the correlated properties of the energy spectra. This assumes that the broad band noise is from propagating fluctuations in mass accretion rate within the hot flow, while the QPO is produced by global Lense-Thirring precession of the same hot flow.

Here we develop this model, making some significant improvements. Firstly we specify that the viscous frequency (equivalently, surface density) in the hot flow has the same form as that measured from numerical simulations of precessing, tilted accretion flows. Secondly, we refine the statistical techniques which we use to fit the model to the data. We re-analyse the PSD from the 1998 rise to outburst of XTE J1550-564 with our new model in order to assess the impact of these changes. We find that the derived outer radii of the hot flow (set by the inner radius of the truncated disc) are rather similar, changing from ~68-13 Rg throughout the outburst rise. However, the more physical assumptions of our new model also allow us to constrain the scale height of the flow. This decreases as the outer radius of the flow decreases, as expected from the spectral evolution. The spectrum steepens in response to the increased cooling as the as the truncation radius sweeps in, so gas pressure support for the flow decreases.

The new model, propfluc, is publically available within the xspec spectral fitting package.

by Lousto, Carlos O. and Zlochower, Yosef
4 pages, 3 figures, revtex 4

We revisit the scenario of the gravitational radiation recoil acquired by the final remnant of a black-hole-binary merger by studying a set of configurations that have components of the spin both aligned with the orbital angular momentum and in the orbital plane. We perform a series of 24 new full numerical simulations for equal-mass and equal-spin-magnitude binaries, but with different spin orientations.

We extend previous recoil fitting formulas to include nonlinear terms in the spins and successfully include both the new and known results. For this new formula the predicted maximum velocity approaches 5000km/s. More importantly, from the astrophysical point of view, it reaches this maximum for spins partially aligned with the orbital angular momentum. The optimal configuration is near an equipartition of the hangup and superkick contributions. This newly discovered contribution to the recoil leads to an important increase of the probabilities of large recoils in generic astrophysical mergers. We measure these probabilities for the case of accretion-aligned spins and find non-negligible probabilities for supermassive black hole encounters leading to recoil velocities of several thousand km/s.

by Nieuwenhuizen, Theo M.
5 pages latex

It is assumed that a galaxy starts as a dark halo of a few million Jeans clusters (JCs), each of which consists of nearly a trillion micro brown dwarfs, MACHOs of earth mass. JCs in the galaxy center heat up their MACHOs by tidal forces, which makes them expand, so that coagulation and star formation occurs. Being continuously fed by matter from bypassing JCs, the star(s) may transform into a super massive black hole. It has a fast $latex t^3$ growth during the first mega years, and a slow $latex t^{1/3}$ growth at giga years. JCs disrupted by a close encounter can provide matter for the bulge. Those that survive can be so agitated that they form stars and become globular star clusters. Thus black holes mostly arise together with galactic bulges in their own environment and are about as old as the oldest globular clusters. The age 13.2 Gyr of the star HE 1523-0901 (Frebel et al. 2007) puts forward that the Galactic halo was fully assembled at that moment. In case of merging super massive black holes the JCs passing near the galactic center provide ideal assistance to overcome the last parsec.

by Favata, Marc
25 pages, 8 figures

The nonlinear gravitational-wave memory causes a time-varying but nonoscillatory correction to the gravitational-wave polarizations. It arises from gravitational-waves that are sourced by gravitational-waves. Previous considerations of the nonlinear memory effect have focused on quasicircular binaries. Here I consider the nonlinear memory from Newtonian orbits with arbitrary eccentricity. Expressions for the waveform polarizations and spin-weighted spherical-harmonic modes are derived for elliptic, hyperbolic, parabolic, and radial orbits. In the hyperbolic, parabolic, and radial cases the nonlinear memory provides a 2.5 post-Newtonian (PN) correction to the leading-order waveforms. This is in contrast to the elliptical and quasicircular cases, where the nonlinear memory corrects the waveform at leading (0PN) order. This difference in PN order arises from the fact that the memory builds up over a short “scattering” timescale in the hyperbolic case, as opposed to a much longer radiation-reaction timescale in the elliptical case. The nonlinear memory corrections presented here complete our knowledge of the leading-order (Peters-Mathews) waveforms for elliptical orbits. These calculations are also relevant for binaries with quasicircular orbits in the present epoch which had, in the past, large eccentricities. Because the nonlinear memory depends sensitively on the past evolution of a binary, I discuss the effect of this early-time eccentricity on the value of the late-time memory in nearly-circularized binaries. I also discuss the observability of large “memory jumps” in a binary’s past that could arise from its formation in a capture process. Lastly, I provide estimates of the signal-to-noise ratio of the linear and nonlinear memories from hyperbolic and parabolic binaries.

by Filloux, Ch. and Pacheco, J. A. de Freitas and Durier, F. and de Araujo, J. C. N.
18 pages, 7 figures, to appear in the IJMPD

The coalescence history of massive black holes has been derived from cosmological simulations, in which the evolution of those objects and that of the host galaxies are followed in a consistent way. The present study indicates that supermassive black holes having masses greater than $latex \sim 10^{9} M_{\odot}$ underwent up to 500 merger events along their history. The derived coalescence rate per comoving volume and per mass interval permitted to obtain an estimate of the expected detection rate distribution of gravitational wave signals (”ring-down”) along frequencies accessible by the planned interferometers either in space (LISA) or in the ground (Einstein). For LISA, in its original configuration, a total detection rate of about $latex 15 yr^{-1}$ is predicted for events having a signal-to-noise ratio equal to 10, expected to occur mainly in the frequency range $latex 4-9 mHz$. For the Einstein gravitational wave telescope, one event each 14 months down to one event each 4 years is expected with a signal-to-noise ratio of 5, occurring mainly in the frequency interval $latex 10-20 Hz$. The detection of these gravitational signals and their distribution in frequency would be in the future an important tool able to discriminate among different scenarios explaining the origin of supermassive black holes.

by Levi, Michele
24 pages, revtex4-1, 5 figures

We calculate the next-to-next-to-leading order (NNLO) spin1-spin2 dynamics of a compact binary evaluated at fourth post-Newtonian (PN) order. We use an effective field theory (EFT) approach, and first demonstrate here the ability of the EFT approach to go at NNLO in the PN corrections of spin effects. The NNLO spin1-spin2 interaction sector includes contributions from diagrams, which are not pure spin1-spin2 diagrams, as they arise from other sectors. These diagrams contribute through the leading order spin accelerations and precessions, that should be first taken into account here. The EFT calculation is carried out in terms of the nonrelativistic gravitational (NRG) fields. The fact that the spin is derivative-coupled adds significantly to the complexity of computations. In particular, for the irreducible two-loop diagrams, which are the most complicated in this sector, irreducible two-loop tensor integrals up to order 4 are encountered. Moreover, not all of the benefits of the NRG fields apply to spin interactions, as all possible diagram topologies are realized at each order of G included. Still, the NRG fields remain advantageous, and thus there was no use of automated computations in this work. Our final result can be reduced, and a NNLO spin1-spin2 Hamiltonian can be derived from it.

by Hartung, Johannes and Steinhoff, Jan
7 pages, submitted to AdP

We present the next-to-next-to-leading order post-Newtonian (PN) spin(1)-spin(2) Hamiltonian for two self-gravitating spinning compact objects. If both objects are rapidly rotating, then the corresponding interaction is comparable in strength to a 4PN effect. The Hamiltonian is checked via the global Poincare algebra with the center-of-mass vector uniquely determined by an ansatz.

by Gualandris, Alessia and Merritt, David
22 pages, 23 figures, submitted to ApJ

We simulate mergers between galaxies containing collisionally-relaxed nuclei around massive black holes (BHs). Our galaxies contain four mass groups, representative of old stellar populations; a primary goal is to understand the distribution of stellar-mass BHs after the merger. Mergers are followed using direct-summation N-body simulations, assuming a mass ratio of 1:3 and two different orbits. Evolution of the massive BH binary is followed until its separation has shrunk by a factor of 20 below the hard-binary separation. During the galaxy merger, large cores are carved out in the stellar distribution, with radii several times the influence radius of the massive BH. Much of the pre-existing mass segregation is erased during this phase. We follow the evolution of the merged galaxies for approximately three, central relaxation times after coalescence of the massive binary; both standard, and top-heavy, mass functions are considered. The cores that were formed in the stellar distribution persist, and the distribution of the stellar-mass black holes evolves against this essentially fixed background. Even after three central relaxation times, these models look very different from the relaxed, multi-mass models that are often assumed to describe the distribution of stars and stellar remnants near a massive BH; in particular, the density of stellar BHs is much smaller than in those models. We discuss the implications of our results for the EMRI problem and for the existence of Bahcall-Wolf cusps.

by Lang, Meagan and Holley-Bockelmann, Kelly and Bogdanovic, Tamara and Amaro-Seoane, Pau and Sesana, Alberto
9 pages, 1 figure, submitted to ApJ. Comments are welcome and may be incorporated into the paper with attribution

Observations of the Galactic Center (GC) have accumulated a multitude of “forensic” evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity — this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

by Volonteri, Marta and Stark, Daniel P.
MNRAS in press

Motivated by recent observational results that focus on high redshift black holes, we explore the effect of scatter and observational biases on the ability to recover the intrinsic properties of the black hole population at high redshift. We find that scatter and selection biases can hide the intrinsic correlations between black holes and their hosts, with ‘observable’ subsamples of the whole population suggesting, on average, positive evolution even when the underlying population is characterized by no- or negative evolution. We create theoretical mass functions of black holes convolving the mass function of dark matter halos with standard relationships linking black holes with their hosts. Under these assumptions, we find that the local MBH – sigma correlation is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), overestimating the number density of all but the most massive black holes. Positive evolution or including scatter in the MBH – sigma correlation makes the discrepancy worse, as it further increases the number density of observable black holes. We notice that if the MBH – sigma correlation at z = 6 is steeper than today, then the mass function becomes shallower. This helps reproducing the mass function of z = 6 black holes proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of order 1/1000) host an active massive black hole at z = 6, or that most AGN are obscured, hindering their detection in optical surveys. Current measurements of the high redshift black hole mass function might be underestimating the density of low mass black holes if the active fraction or luminosity are a function of host or black hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper MBH – sigma relation at high redshift.

by Kelly, Bernard J. and Baker, John G. and Boggs, William D. and McWilliams, Sean T. and Centrella, Joan
19 pages, 17 figures

We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [Baker et al., Phys. Rev. D 78:044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission — an implicit rotating source (IRS). We further apply the late-time merger-ringdown model for the rotational frequency introduced in Baker et al. (2008), along with an improved amplitude model appropriate for the dominant (2,+/- 2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above ~ 150 MSun. This model may be directly applicable to gravitational-wave detection of intermediate-mass black hole mergers.

by Ohme, Frank and Hannam, Mark and Husa, Sascha
15 pages, 9 figures, PDFLaTeX

With recent advances in post-Newtonian (PN) theory and numerical relativity (NR) it has become possible to construct inspiral-merger-ringdown gravitational waveforms from coalescing compact binaries by combining both descriptions into one complete hybrid signal. It is important to estimate the error of such waveforms. Previous studies have identified the PN contribution as the dominant source of error, which can be reduced by incorporating longer NR simulations. There are two outstanding issues that make it difficult to determine the minimum simulation length necessary to produce suitably accurate hybrids: (1) the relevant criteria for a signal search is the mismatch between the true waveform and a set of model waveforms, optimized over all waveforms in the model. For discrete hybrids this optimization is not possible. (2) these calculations require that NR waveforms already exist, while ideally we would like to know the necessary length before performing the simulation. Here we overcome these difficulties by developing a general procedure that allows us to estimate hybrid mismatch errors without numerical data, and to optimize them over all physical parameters. Using this procedure we find that, contrary to some earlier studies, ~10 NR orbits before merger allow for the construction of waveform families that are accurate enough for detection in a broad range of parameters, only excluding highly spinning, unequal-mass systems. Nonspinning binaries, even with high mass-ratio (>20) are well modeled for astrophysically reasonable component masses. In addition, the parameter bias is only of the order of 1% for total mass and symmetric mass-ratio and less than 0.1 for the dimensionless spin magnitude. We take the view that similar NR waveform lengths will remain the state of the art in the Advanced detector era, and begin to assess the limits of the science that can be done with them.

by Kulkarni, Girish and Loeb, Abraham
18 pages, 14 figures, 2 tables; submitted to MNRAS; comments welcome

We examine the formation of groups of multiple supermassive black holes (SMBHs) in gas-poor galactic nuclei due to the high merger rate of galaxies at high redshifts. We calculate the relative likelihood of binary, triple, and quadruple SMBH systems, by considering the timescales for relevant processes and combining merger trees with N-body simulations for the dynamics of stars and SMBHs in galactic nuclei. Typical haloes today with mass $latex M_0\approx 10^{14}$ M$latex _\odot$ have an average mass $latex M_{z=6}=5\times 10^{11}$ M$latex _\odot$ at $latex z\sim 6$, while rare haloes with current mass $latex M_0\gtrsim 10^{15}$ M$latex _\odot$ have an average mass $latex M_{z=6}=5\times 10^{12}$ M$latex _\odot$ at that redshift. These cluster-size haloes are expected to host single galaxies at $latex z\sim 6$. We expect about 30% galaxies within haloes with present-day mass $latex M_0\approx 10^{14}$ M$latex _\odot$ to contain more than two SMBHs at redshifts $latex 2\lesssim z\lesssim 6$. For larger present-day haloes, with $latex M_0\gtrsim 10^{15}$ M$latex _\odot$, this fraction is almost 60%. The existence of multiple SMBHs at high redshifts can potentially explain the mass deficiencies observed in the cores of massive elliptical galaxies, which are up to 5 times the mass of their central BHs. Multiple SMBHs would also lead to an enhanced rate of tidal disruption of stars, modified gravitational wave signals compared to isolated BH binaries, and slingshot ejection of SMBHs from galaxies at high speeds in excess of 2000 km s$latex ^{-1}$.

by Babiuc, M. C. and Winicour, J. and Zlochower, Y.
11 pages, 7 figures

In this work, we present a work in progress towards an efficient and economical computational module which interfaces between Cauchy and characteristic evolution codes. Our goal is to provide a standardized waveform extraction tool for the numerical relativity community which will allow CCE to be readily applied to a generic Cauchy code. The tool provides a means of unambiguous comparison between the waveforms generated by evolution codes based upon different formulations of the Einstein equations and different numerical approximation.

by Nagar, Alessandro
11 pages, no figures. Submitted to Phys. Rev. D

Building on the recently computed next-to-next-to-leading order (NNLO) post-Newtonian (PN) spin-orbit Hamiltonian for spinning binaries \cite{Hartung:2011te} we extend the effective-one-body (EOB) description of the dynamics of two spinning black-holes to NNLO in the spin-orbit interaction. The calculation that is presented extends to NNLO the next-to-leading order (NLO) spin-orbit Hamiltonian computed in Ref. \cite{Damour:2008qf}. The present EOB Hamiltonian reproduces the spin-orbit coupling through NNLO in the test-particle limit case. In addition, in the case of spins parallel or antiparallel to the orbital angular momentum, when circular orbits exist, we find that the inclusion of NNLO spin-orbit terms moderates the effect of the NLO spin-orbit coupling.

by Tiec, Alexandre Le and Mroué, Abdul H. and Barack, Leor and Buonanno, Alessandra and Pfeiffer, Harald P. and Sago, Norichika and Taracchini, Andrea
5 pages, 3 figures

The general relativistic periastron advance of non-spinning black hole binaries on quasi-circular orbits has been computed using black hole perturbation theory, post-Newtonian expansions, and the effective-one-body formalism. We compare these approximations with accurate numerical relativity simulations of mass ratios 1/8 < m1/m2 m1m2/(m1+m2)^2. The effective-one-body prediction also agrees very well over the entire mass-ratio range considered.