by Dauser, T. and Wilms, J. and Reynolds, C. S. and Brenneman, L. W.
7 pages, 6 figures; accepted by MNRAS for Publication

We present an extended scheme for the calculation of the profiles of emission lines from accretion discs around rotating black holes. The scheme includes discs with angular momenta which are parallel and antiparallel with respect to the black hole’s angular momentum, as both configurations are assumed to be stable (King et al., 2005). We discuss line shapes for such discs and present a code for modelling observational data with this scheme in X-ray data analysis programs. Based on a Green’s function approach, an arbitrary radius dependence of the disc emissivity and arbitrary limb darkening laws can be easily taken into account, while the amount of precomputed data is significantly reduced with respect to other available models.

by Wang, J. -M. and Yan, C. -S. and Gao, H. -Q. and Hu, C. and Li, Y. -R. and Zhang, S.
emulateapj.sty, 5 page, 4 figures (in press)

Self-gravitating accretion disks collapse to star-forming(SF) regions extending to the inner edge of the dusty torus in active galactic nuclei (AGNs). A full set of equations including feedback of star formation is given to describe the dynamics of the regions. We explore the role of supernovae explosion (SNexp), acting to excite turbulent viscosity, in the transportation of angular momentum in the regions within 1pc scale. We find that accretion disks with typical rates in AGNs can be driven by SNexp in the regions and metals are produced spontaneously. The present model predicts a metallicity–luminosity relationship consistent with that observed in AGNs. As relics of SF regions, a ring (or belt) consisting of old stars remains for every episode of supermassive black hole activity. We suggest that multiple stellar rings with random directions interact and form a nuclear star cluster after episodes driven by star formation.

by Fu, Wen and Lai, Dong
15 pages, 10 figures, submitted to MNRAS

Hot accretion tori around a compact object are known to be susceptible to a global hydrodynamical instability, the so-called Papaloizou-Pringle (PP) instability, arising from the interaction of non-axisymmetric waves across the corotation radius, where the wave pattern speed matches the fluid rotation rate. However, accretion tori produced in various astrophysical situations (e.g., collapsars and neutron star binary mergers) are likely to be highly magnetized. We study the effect of magnetic fields on the PP instability in incompressible tori with various magnetic strengths and structures. In general, toroidal magnetic fields have significant effects on the PP instability: For thin tori (with the fractional width relative to the outer torus radius much less than unity), the instability is suppressed at large field strengths with the corresponding toroidal Alfven speed $latex v_{A\phi}\go 0.2r\Omega$ (where $latex \Omega$ is the flow rotation rate). For thicker tori (with the fractional width of order 0.4 or larger), which are hydrodynamically stable, the instability sets in for sufficiently strong magnetic fields (with $latex v_{A\phi}\go 0.2 r\Omega$). Our results suggest that highly magnetized accretion tori may be subjected to global instability even when it is stable against the usual magneto-rotational instability.

by Schnittman, Jeremy D.
10 pages, 4 figures, submitted to ApJ; comments welcome

We calculate the location and stability of the L_4 and L_5 Lagrange equilibrium points in the circular restricted three-body problem as the binary system evolves via gravitational radiation losses. Relative to the purely Newtonian case, we find that the L_4 equilibrium point moves towards the secondary mass and becomes slightly less stable, while the L_5 point moves away from the secondary and gains in stability. We discuss a number of astrophysical applications of these results, in particular as a mechanism for producing electromagnetic counterparts to gravitational-wave signals.

by Tanaka, Takamitsu and Haiman, Zoltan and Menou, Kristen
27 pages, 5 figures, submitted to ApJ

The coalescence of a supermassive black hole binary (SMBHB) is thought to be accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of the surrounding circumbinary gas disk after the merger. It has been proposed that once the merger has been detected in gravitational waves (GWs) by LISA, follow-up EM searches for this afterglow can help identify the EM counterpart of the LISA source. Here we study whether the afterglows may be sufficiently bright and numerous to be detectable in EM surveys alone. The viscous afterglow, which lasts for years to decades for SMBHBs in LISA’s sensitivity window, is characterized by rapid increases in both the bolometric luminosity and in the spectral hardness of the source. If quasar activity is triggered by the same major galaxy mergers that produce SMBHBs, then the afterglow could be interpreted as a signature of the birth of a quasar. Using an idealized model for the post-merger viscous spreading of the circumbinary disk and the resulting light curve, and using the observed luminosity function of quasars as a proxy for the SMBHB merger rate, we delineate the survey requirements for identifying such birthing quasars. If circumbinary disks have a high disk surface density and viscosity, an all-sky soft X-ray survey with a sensitivity of ~10%/yr. If >1% of the X-ray emission is reprocessed into optical frequencies, birthing quasars could also be identified in optical transient surveys such as the LSST. Distinguishing a birthing quasar from other variable sources may be facilitated by the monotonic hardening of its spectrum, but will likely remain challenging. This reinforces the notion that joint EM-plus-GW observations offer the best prospects for identifying the EM signatures of SMBHB mergers.

by Shang, Cien and Haiman, Zoltan
16 pages, 9 tables, 12 figures, submitted to MNRAS

The gravitational waves (GWs) emitted by inspiraling binary black holes, expected to be detected by the Laser Interferometer Space Antenna (LISA), could be used to determine the luminosity distance to these sources with the unprecedented precision of <~ 1%. We study cosmological parameter constraints from such standard sirens, in the presence of gravitational lensing by large-scale structure. Lensing introduces magnification with a probability distribution function (PDF) whose shape is highly skewed and depends on cosmological parameters. We use Monte-Carlo simulations to generate mock samples of standard sirens, including a small intrinsic scatter, as well as the additional, larger scatter from lensing, in their inferred distances. We derive constraints on cosmological parameters, by simultaneously fitting the mean and the distribution of the residuals on the distance vs redshift (d_L – z) Hubble diagram. We find that for standard sirens at redshift z ~ 1, the sensitivity to a single cosmological parameter, such as the matter density Omega_m, or the dark energy equation of state w, is ~ 50%-80% tighter when the skewed lensing PDF is used, compared to the sensitivity derived from a Gaussian PDF with the same variance. When these two parameters are constrained simultaneously, the skewness yields a further enhanced improvement (by ~ 120%), owing to the correlation between the parameters. The sensitivity to the amplitude of the matter power spectrum, sigma_8 from the cosmological dependence of the PDF alone, however, is ~ 20% worse than that from the Gaussian PDF. At higher redshifts, the PDF resembles a Gaussian more closely, and the effects of the skewness become less prominent. These results highlight the importance of obtaining an accurate and reliable PDF of the lensing convergence, in order to realize the full potential of standard sirens as cosmological probes.

by Lesur, G. and Ogilvie, G. I.
6 pages, 5 figures, accepted in MNRAS

The mechanism of angular momentum transport in accretion discs has long been debated. Although the magnetorotational instability appears to be a promising process, poorly ionized regions of accretion discs may not undergo this instability. In this letter, we revisit the possibility of transporting angular momentum by turbulent thermal convection. Using high-resolution spectral methods, we show that strongly turbulent convection can drive outward angular momentum transport at a rate that is, under certain conditions, compatible with observations of discs. We find however that the angular momentum transport is always much weaker than the vertical heat transport. These results indicate that convection might be another way to explain global disc evolution, provided that a sufficiently unstable vertical temperature profile can be maintained.

by Zanotti, Olindo and Rezzolla, Luciano and Del Zanna, Luca and Palenzuela, Carlos
17 pages, 11 figures, submitted to MNRAS, movies available at http://numrel.aei.mpg.de/Visualisations/Archive/BinaryBlackHoles/EMCounterparts/EMCounterparts.html

We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. More specifically, we perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a “shock detector”, thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. In contrast with what done in similar works, we here question the estimates of the bremsstrahlung luminosity when computed without properly taking into account the radiation transfer, thus yielding cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic analogue of the isothermal evolution of \citet{Corrales2009}, that the luminosity produced can reach a peak value above $latex L \simeq 10^{43} {\rm erg/s} $ at about $latex \sim 20 {\rm d}$ after the merger of a binary with total mass $latex M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system.

by Sorathia, Kareem A. and Reynolds, Christopher S. and Armitage, Philip J.
8 Pages, 7 Figures ApJ, In Press

We analyze a suite of global magnetohydrodynamic (MHD) accretion disk simulations in order to determine whether scaling laws for turbulence driven by the magnetorotational instability, discovered via local shearing box studies, are globally robust. The simulations model geometrically-thin disks with zero net magnetic flux and no explicit resistivity or viscosity. We show that the local Maxwell stress is correlated with the self-generated local vertical magnetic field in a manner that is similar to that found in local simulations. Moreover, local patches of vertical field are strong enough to stimulate and control the strength of angular momentum transport across much of the disk. We demonstrate the importance of magnetic linkages (through the low-density corona) between different regions of the disk in determining the local field, and suggest a new convergence requirement for global simulations — the vertical extent of the corona must be fully captured and resolved. Finally, we examine the temporal convergence of the average stress, and show that an initial long-term secular drift in the local flux-stress relation dies away on a time scale that is consistent with turbulent mixing of the initial magnetic field.

by Curir, Anna and de Romeri, Valentina and Murante, Giuseppe
10 pages, 8 figures, accepted for pubblication in “Astrophysics and Space Science”

The bar formation is still an open problem in modern astrophysics. In this paper we present numerical simulation performed with the aim of analyzing the growth of the bar instability inside stellar-gaseous disks, where the star formation is triggered, and a central black hole is present. The aim of this paper is to point out the impact of such a central massive black hole on the growth of the bar. We use N-body-SPH simulations of the same isolated disk-to-halo mass systems harboring black holes with different initial masses and different energy feedback on the surrounding gas. We compare the results of these simulations with the one of the same disk without black hole in its center. We make the same comparison (disk with and without black hole) for a stellar disk in a fully cosmological scenario. A stellar bar, lasting 10 Gyrs, is present in all our simulations. The central black hole mass has in general a mild effect on the ellipticity of the bar but it is never able to destroy it. The black holes grow in different way according their initial mass and their feedback efficiency, the final values of the velocity dispersions and of the black hole masses are near to the phenomenological constraints.

by King, Andrew
invited review, IAU Symposium 267, Co-Evolution of Central Black Holes and Galaxies, B.M. Peterson, R.S. Somerville, and T. Storchi-Bergmann, eds typos in eq (2.2) corrected

I review accretion and outflow in active galactic nuclei. Accretion appears to occur in a series of very small–scale, chaotic events, whose gas flows have no correlation with the large–scale structure of the galaxy or with each other. The accreting gas has extremely low specific angular momentum and probably represents only a small fraction of the gas involved in a galaxy merger, which may be the underlying driver.

Eddington accretion episodes in AGN must be common in order for the supermassive black holes to grow. I show that they produce winds with velocities $latex v \sim 0.1c$ and ionization parameters implying the presence of resonance lines of helium– and hydrogenlike iron. The wind creates a strong cooling shock as it interacts with the interstellar medium of the host galaxy, and this cooling region may be observable in an inverse Compton continuum and lower–excitation emission lines associated with lower velocities. The shell of matter swept up by the shocked wind stalls unless the black hole mass has reached the value $latex M_{\sigma}$ implied by the $latex M – \sigma$ relation. Once this mass is reached, further black hole growth is prevented. If the shocked gas did not cool as asserted above, the resulting (`energy-driven’) outflow would imply a far smaller SMBH mass than actually observed. Minor accretion events with small gas fractions can produce galaxy-wide outflows, including fossil outflows in galaxies where there is little current AGN activity.

by Dutan, Ioana
12 pages, 10 figures, submitted to MNRAS

We present an alternative model for launching relativistic jets in active galactic nuclei (AGN) from an accreting Kerr black hole (BH) by converting the accretion disc energy into jet energy, when the rotational energy of the BH is transferred to the inner disc by closed magnetic field lines which connects the BH to the disc (BH-disc magnetic connection). In this way, the available disc energy is increased by the BH rotational energy. We assume that the BH may undergo recurring episodes of its activity with: (i) a first phase when accretion power dominates, and (ii) a second phase when BH spin-down power dominates. In both cases the jet is driven by a low-luminosity, (geometrically) thin accretion disc, as the disc energy is used to launch the jet. We use the general relativistic conservation laws to calculate the mass flow rate into the jets, the launching power of the jets, and the angular momentum transported by the jets. We consider BHs with a spin parameter $latex a_* \geqslant 0.95$, so that the jets are launched from the region inside of the BH ergosphere. The angular momentum removed from the accretion disc is carried away by the disc particles that ultimately form the jets. As far as the BH is concerned, it can (i) spin up by accreting matter and (ii) spin down due to the magnetic counter-acting torque on the BH. We found that a stationary state of the BH ($latex a_* = $ const) can be reached if the mass accretion rate is larger than $latex \dot{m} \sim 0.001$. The maximum value of the BH spin parameter depends on $latex \dot{m}$ being less but close to 0.9982 (Thorne’s model). In addition, the maximum AGN lifetime can be much longer than $latex \sim 10^{7}$ yr when using the BH spin-down power. This result is consistent with the estimation of the maximum AGN lifetime when the AGN output power is provided by the Blandford–Znajek mechanism.

by Noble, Scott C. and Krolik, Julian H. and Hawley, John F.
Accepted for publication in ApJ, 52 pages, 38 figures, AASTEX. High-resolution versions can be found at the following links: http://ccrg.rit.edu/~scn/papers/schwarzstress.ps, http://ccrg.rit.edu/~scn/papers/schwarzstress.pdf

We explore the parameter dependence of inner disk stress in black hole accretion by contrasting the results of a number of simulations, all employing 3-d general relativistic MHD in a Schwarzschild spacetime. Five of these simulations were performed with the intrinsically conservative code HARM3D, which allows careful regulation of the disk aspect ratio, H/R; our simulations span a range in H/R from 0.06 to 0.17. We contrast these simulations with two previously reported simulations in a Schwarzschild spacetime in order to investigate possible dependence of the inner disk stress on magnetic topology. In all cases, much care was devoted to technical issues: ensuring adequate resolution and azimuthal extent, and averaging only over those time-periods when the accretion flow is in approximate inflow equilibrium. We find that the time-averaged radial-dependence of fluid-frame electromagnetic stress is almost completely independent of both disk thickness and poloidal magnetic topology. It rises smoothly inward at all radii (exhibiting no feature associated with the ISCO) until just outside the event horizon, where the stress plummets to zero. Reynolds stress can also be significant near the ISCO and in the plunging region; the magnitude of this stress, however, depends on both disk thickness and magnetic topology. The two stresses combine to make the net angular momentum accreted per unit rest-mass 7-15% less than the angular momentum of the ISCO.

by Hobbs, Alexander and Nayakshin, Sergei and Power, Chris and King, Andrew
19 pages, 17 figures. Submitted to MNRAS

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. However, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, inside bulges the gas velocity distribution can be randomised by the velocity kicks due to feedback from star formation. Here we perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of “ballistic” accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. This prevents the efficient gas mixing that was found in the simulations without turbulence, and allows a fraction of gas to impact the innermost boundary of the simulations directly. Using the ballistic approximation, we derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this “ballistic” mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars. We therefore argue that star formation in bulges can be a strong catalyst for SMBH accretion.

by Yuan, Feng and Bu, Defu
9 pages, 9 figures; submitted to MNRAS

The most important finding of two-dimensional hydrodynamical simulations of hot accretion flows is that the flow is convectively unstable, because of its inward increase of entropy. As a result, the profile of the mass accretion rate is a function of radius, i.e., only a small fraction of accretion gas available at the outer boundary can finally fall onto the black hole, while the rest is lost in the convective outflows. Radiation is usually neglected in these simulations. When the radiative cooling becomes more and more important, the entropy will increase slower inward. The entropy can even decrease when the radiation becomes stronger than the viscous heating, i.e, the flow enters into the luminous hot accretion flow regime. In the present paper, we investigate the convective instability and correspondingly the profile of accretion rate in the presence of strong radiative cooling by performing two-dimensional hydrodynamical numerical simulation. This problem is important because the profile of the mass accretion rate determines the observational appearance of accretion flows, the growth of black hole, and the evolution of black hole spin. We find that the flow is still strongly convectively unstable, and the radial profile of accretion rate changes little compared to the case of non-radiative flow. This is because the gradient of entropy in the gravitational direction still increases inward although the gradient of entropy decreases.

by Lattanzi, Massimiliano and Montani, Giovanni
6 pages, 4 figures. To appear in Europhysics Letters

We provide a new analysis of the system of partial differential equations describing the radial and vertical equilibria of the plasma in accretion disks. In particular, we show that the partial differential system can be separated once a definite, oscillatory (or hyperbolic) form for the radial dependence of the relevant physical quantities is assumed. The system is thus reduced to an ordinary differential system in the vertical dimensionless coordinate. The resulting equations can be integrated analytically in the limit of small magnetic pressure. We complete our analysis with a direct numerical integration of the more general case. The main result is that a ring-like density profile (i.e., radial oscillations in the mass density) can appear even in the limit of small magnetic pressure.

by Mayer, Lucio and Kazantzidis, Stelios and Escala, Andres and Callegari, Simone
26 pages, 4 Figures, submitted to Nature (includes Supplementary Information)

Observations of distant bright quasars suggest that billion solar mass supermassive black holes (SMBHs) were already in place less than a billion years after the Big Bang. Models in which light black hole seeds form by the collapse of primordial metal-free stars cannot explain their rapid appearance due to inefficient gas accretion. Alternatively, these black holes may form by direct collapse of gas at the center of protogalaxies. However, this requires metal-free gas that does not cool efficiently and thus is not turned into stars, in contrast with the rapid metal enrichment of protogalaxies. Here we use a numerical simulation to show that mergers between massive protogalaxies naturally produce the required central gas accumulation with no need to suppress star formation. Merger-driven gas inflows produce an unstable, massive nuclear gas disk. Within the disk a second gas inflow accumulates more than 100 million solar masses of gas in a sub-parsec scale cloud in one hundred thousand years. The cloud undergoes gravitational collapse, which eventually leads to the formation of a massive black hole. The black hole can grow to a billion solar masses in less than a billion years by accreting gas from the surrounding disk.

by Hopkins, Philip F. and Quataert, Eliot
15 figures, 30 pages. Submitted to MNRAS. Movies of the simulations described here can be found at http://www.cfa.harvard.edu/~phopkins/Site/Movies_zoom.html

We use multi-scale SPH simulations to follow the inflow of gas from galactic scales to <0.1pc, where the gas begins to resemble a traditional Keplerian accretion disk. The key ingredients are gas, stars, black holes (BHs), self-gravity, star formation, and stellar feedback. We use ~100 simulations to survey a large parameter space of galaxy properties and subgrid models for the ISM physics. We generate initial conditions for our simulations of galactic nuclei (<~300pc) using galaxy scale simulations, including both major mergers and isolated bar-(un)stable disk galaxies. For sufficiently gas-rich, disk-dominated systems, a series of gravitational instabilities generates large accretion rates of up to 1-10 M_sun/yr onto the BH (at <~10pc, our simulations resemble the ‘bars within bars’ model, but the gas exhibits diverse morphologies, including spirals, rings, clumps, and bars; their duty cycle is modest, complicating attempts to correlate BH accretion with nuclear morphology. At ~1-10pc, the gravitational potential becomes dominated by the BH and bar-like modes are no longer present. However, the gas becomes unstable to a standing, eccentric disk or a single-armed spiral mode (m=1), driving the gas to sub-pc scales. Proper treatment of this mode requires including star formation and the self-gravity of both the stars and gas. We predict correlations between BHAR and SFR at different galactic nuclei: nuclear SF is more tightly coupled to AGN activity, but correlations exist at all scales.

by Shcherbakov, Roman V. and Baganoff, Frederick K.
3 pages, 2 figures; to be published in “The Galactic Center: a Window to the Nuclear Activity of Disk Galaxies”, ed. Mark Morris et al., Astron. Soc. Pacific, 2010

We propose a 2-temperature radial dynamical model of plasma flow near Sgr A* and fit the bremsstrahlung emission to extensive quiescent X-Ray Chandra data. The model extends from several arcseconds to black hole (BH) gravitational radius, describing the outer accretion flow together with the infalling region. The model incorporates electron heat conduction, relativistic heat capacity of particles and feeding by stellar winds. Stellar winds from each star are considered separately as sources of mass, momentum and energy. Self-consistent search for the stagnation and sonic points is performed. Most of gas is found to outflow from the region. The accretion rate is limited to below 1% of Bondi rate due to the effect of thermal conduction enhanced by entropy production in a turbulent flow. The X-Ray brightness profile proves too steep near the BH, thus a synchrotron self-Compton point source is inferred with luminosity L=3×10^32erg/s. We fit the sub-mm emission from the inner flow, thus aiming at a single model of Sgr A* accretion suitable at any radius.

by Palenzuela, Carlos and Lehner, Luis and Yoshida, Shin
12 pages

In addition to producing loud gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary’s dynamics induce a variability in possible electromagnetically induced emissions as well as an enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

by Zhang, Dong and Dai, Z. G.
61 pages, 14 figures, 4 tables, submitted to ApJ

(Abridged) The hyperaccreting neutron star or magnetar disks cooled via neutrino emission can be a candidate of gamma-ray burst (GRB) central engines. The strong field $latex \geq10^{15}-10^{16}$ G of the magnetar can play a significant role in affecting the disk properties and even lead to the funnel accretion process. We investigate the effects of strong fields on the disks around magnetars, and discuss implications of such accreting magnetar systems for GRB and GRB-like events. We discuss quantum effects of the strong fields on the disk, and use the MHD conservation equations to describe the behavior of the disk flow coupled with a large scale field, which is generated by the star-disk interaction. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity, and change the electron fraction and degeneracy state significantly. A magnetized disk is always viscously stable outside the Alfv\’{e}n radius, but will be thermally unstable near the Alfv\’{e}n radius where the magnetic field plays a more important role in transferring the angular momentum and heating the disk than the viscous stress. The funnel accretion process will be only important for an extremely strong field, which creates a magnetosphere inside the Alfv\’{e}n radius and truncates the plane disk. Because of higher temperature and more concentrated neutrino emission of the magnetar surface ring-like belt region covered by funnel accretion, the neutrino annihilation rate from the accreting magnetars can be much higher than that from accreting neutron stars without fields. Furthermore, the neutrino annihilation mechanism and the magnetically-driven pulsar wind from the magnetar surface can work together to generate and feed an ultra-relativistic jet along the stellar magnetic poles.

by Krolik, Julian H.
14 pages. Ap J, in press

Although recent work in numerical relativity has made tremendous strides in quantifying the gravitational wave luminosity of black hole mergers, very little is known about the electromagnetic luminosity that might occur in immediate conjunction with these events. We show that whenever the heat deposited in the gas near a pair of merging black holes is proportional to its total mass, and the surface density of the gas in the immediate vicinity is greater than the (quite small) amount necessary to make it optically thick, the characteristic scale of the luminosity emitted in direct association with the merger is the Eddington luminosity independent of the gas mass. The duration of the photon signal is proportional to the gas mass, and is generally rather longer than the merger event. At somewhat larger distances, dissipation associated with realigning the gas orbits to the new spin orientation of the black hole can supplement dissipation of the energy gained from orbital adjustment to the mass lost in gravitational radiation; these two heat sources can combine to augment the electromagnetic radiation over longer timescales.

by Bode, Tanja and Haas, Roland and Bogdanovic, Tamara and Laguna, Pablo and Shoemaker, Deirdre
19 pages, 14 figures

Coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures from coalescence events of supermassive black holes are the next observational grand challenge. Such detections will provide the means to study cosmological evolution and accretion processes associated with these gargantuan compact objects. More generally, the observations will enable testing general relativity in the strong, nonlinear regime and will provide independent cosmological measurements to high precision. Understanding the conditions under which coincidences of EM and GW signatures arise during supermassive black hole mergers is therefore of paramount importance. As an essential step towards this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning supermassive black hole binaries in a gas cloud. We find that variable EM signatures correlated with GWs can arise in merging systems as a consequence of shocks and accretion combined with the effect of relativistic beaming. The most striking EM variability is observed for systems where spins are aligned with the orbital axis and where orbiting black holes form a stable set of density wakes, but all systems exhibit some characteristic signatures that can be utilized in searches for EM counterparts. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, calculated luminosities imply that they may be identified by EM searches to z = 1, while lower mass systems and binaries immersed in low density ambient gas can only be detected in the local universe.

by Shen, Yue and Loeb, Abraham
11 pages. Submitted to MNRAS

Double-peaked broad lines in Active Galactic Nuclei (AGNs) may indicate the existence of a supermassive black hole (SMBH) binary whose two broad line regions (BLRs) contribute together to the line profile. An alternative interpretation of the double-peaked broad line feature is a disk origin for the line emission. We calculate the expected broad line profiles for a SMBH binary with various separations, using simple BLR models. Under reasonable assumptions that both BLRs are illuminated by the two black holes (BHs) and that the ionizing flux at the BLR position is roughly constant, we confirm the emergence of double-peaked features and radial velocity drifts of the two peaks due to the binary orbital motion when the two BHs are close enough such that the light-of-sight orbital velocity difference is larger than the FWHM of individual broad components. However, when the two BHs are even closer such that the two BLRs are no longer distinct, the line profile becomes more complex and there are no coherent radial velocity drifts in the two peaks with time. We discuss the temporal variations of the broad line profile for binary SMBHs and highlight the different behaviors of reverberation mapping in the binary and disk emitter scenarios.

by McWilliams, Sean T. and Thorpe, James Ira and Baker, John G. and Kelly, Bernard J.
16 pages, 9 figures, submitted to Phys. Rev. D

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin).

by Giri, Kinsuk and Chakrabarti, Sandip K. and Samanta, Madan M. and Ryu, Dongsu
19 pages, 13 figures, 1 Table MNRAS (In press)

We study the accretion processes on a black hole by numerical simulation. We use a grid based finite difference code for this purpose. We scan the parameter space spanned by the specific energy and the angular momentum and compare the time-dependent solutions with those obtained from theoretical considerations. We found several important results (a) The time dependent flow behaves close to a constant height model flow in the pre-shock region and a flow with vertical equilibrium in the post-shock region. (c) The infall time scale in the post-shock region is several times higher than the free-fall time scale. (b) There are two discontinuities in the flow, one being just outside of the inner sonic point. Turbulence plays a major role in determining the locations of these discontinuities. (d) The two discontinuities oscillate with two different frequencies and behave as a coupled harmonic oscillator. A Fourier analysis of the variation of the outer shock location indicates higher power at the lower frequency and lower power at the higher frequency. The opposite is true when the analysis of the inner shock is made. These behaviours will have implications in the spectral and timing properties of black hole candidates.

by Tanaka, Takamitsu and Menou, Kristen
46 pages, 8 figures. Submitted to ApJ

The Laser Interferometer Space Antenna (LISA) will detect gravitational wave signals from coalescing pairs of massive black holes in the total mass range (10^5 – 10^7)/Msol out to cosmological distances. Identifying and monitoring the electromagnetic counterparts of these events would enable cosmological studies and offer new probes of gas physics around well-characterized massive black holes. Milosavljevic & Phinney (2005) proposed that a circumbinary disk around a binary of mass ~10^6 Msol will emit an accretion-powered X-ray afterglow approximately one decade after the gravitational wave event. We revisit this scenario by using Green’s function solutions to calculate the temporal viscous evolution and the corresponding electromagnetic signature of the circumbinary disk. Our calculations suggest that an electromagnetic counterpart may become observable as a rapidly brightening source soon after the merger, i.e. several years earlier than previously thought. The afterglow can reach super-Eddington luminosities without violating the local Eddington flux limit. It is emitted in the soft X-ray by the innermost circumbinary disk, but it may be partially reprocessed at optical and infrared frequencies. We also find that the spreading disk becomes increasingly geometrically thick close to the central object as it evolves, indicating that the innermost flow could become advective and radiatively inefficient, and generate a powerful outflow. We conclude that the mergers of massive black holes detected by LISA offer unique opportunities for monitoring on humanly tractable timescales the viscous evolution of accretion flows and the emergence of outflows around massive black holes with precisely known masses, spins and orientations.

by Farris, Brian D. and Liu, Yuk Tung and Shapiro, Stuart L.
33 pages, 24 Figures, two tables. Submitted to PRD

Merging supermassive black hole-black hole (BHBH) binaries produced in galaxy mergers are promising sources of detectable gravitational waves. If such a merger takes place in a gaseous environment, there is a possibility of a simultaneous detection of electromagnetic and gravitational radiation, as the stirring, shock heating and accretion of the gas may produce variability and enhancements in the electromagnetic flux. Such a simultaneous detection can provide a wealth of opportunities to study gravitational physics, accretion physics, and cosmology. We investigate this scenario by performing fully general relativistic, hydrodynamic simulations of merging, equal-mass, nonspinning BHBH binaries embedded in gas clouds. We evolve the metric using the BSSN formulation with standard moving puncture gauge conditions and handle the hydrodynamics via a high-resolution shock-capturing (HRSC) scheme. We consider both “binary Bondi accretion” in which the binary is at rest relative to the ambient gas cloud, as well as “binary Bondi-Hoyle-Lyttleton accretion” in which the binary moves relative to the gas cloud. The gas cloud is assumed to be homogeneous far from the binary and governed by a \Gamma-law equation of state. We vary \Gamma between 4/3 and 5/3. For each simulation, we compute the gas flow and accretion rate and estimate the electromagnetic luminosity due to bremsstrahlung and synchrotron emission. We find evidence for significant enhancements in both the accretion rate and luminosity over values for a single black hole of the same mass as the binary. We estimate that this luminosity enhancement should be detectable by LSST for a 10^6 M_sun binary in a hot gas cloud of density n~10/cm^3 and temperature T~10^6 K at z=1, reaching a maximum of L~3×10^43 erg/s, with the emission peaking in the visible band.

by Shapiro, Stuart L.
9 pages, 5 figures; submitted to PRD

Tidal torques from a binary black hole (BHBH) empty out the central regions in any circumbinary gaseous accretion disk. The balance between tidal torques and viscosity maintain the inner edge of the disk at a radius r ~ 1.5a — 2a, where a is the binary semimajor axis. Eventually, the inspiraling binary decouples from disk and merges, leaving behind a central hollow (”donut hole”) in the disk orbiting the remnant black hole. We present a simple, time-dependent, Newtonian calculation that follows the secular (viscous) evolution of the disk as it fills up the hollow down to the black hole innermost stable circular orbit and then relaxes to stationary equilibrium. We use our model to calculate the electromagnetic radiation (”afterglow”) spectrum emitted during this transient accretion epoch. Observing the temporal increase in the total electromagnetic flux and the hardening of the spectrum as the donut hole fills may help confirm a BHBH merger detected by a gravitational wave interferometer. We show how the very existence of the initial hollow can lead to super-Eddington accretion during this secular phase if the rate is not very far below Eddington prior to decoupling. Our model, though highly idealized, may be useful in establishing some of the key parameters, thermal emission features and scalings that characterize this transient. It can serve as a guide in the design and calibration of future radiation-magnetohydrodynamic simulations in general relativity.

by Bhattacharya, Debbijoy and Ghosh, Shubhrangshu and Mukhopadhyay, Banibrata
26 pages including 7 figures; submitted in ApJ

The mechanism by which outflows and plausible jets are driven from black hole systems, still remains observationally elusive. Notwithstanding, several observational evidences and deeper theoretical insights reveal that accretion and outflow/jet are strongly correlated. Here, we model an advective disk-outflow coupled dynamics, incorporating explicitly the vertical flux. Inter-connecting dynamics of outflow and accretion essentially upholds the conservation laws. We investigate the properties of the disk-outflow surface and its strong dependence on the rotation parameter of the black hole. The energetics of disk-outflow strongly depend on mass, accretion rate and spin of the black holes. The model clearly shows that the outflow power extracted from the disk increases strongly with the spin of the black hole, inferring that the power of the observed astrophysical jets has a proportional correspondence with the spin of the central object. In case of blazars (BL Lacs and Flat Spectrum Radio Quasars), most of their emission are believed to be originated from their jets. It is observed that BL Lacs are relatively low luminous than Flat Spectrum Radio Quasars (FSRQs). The luminosity might be linked to the power of the jet, which in turn reflects that the nuclear regions of the BL Lac objects have a relatively low spinning black hole compared to that in the case of FSRQ. If the extreme gravity is the source to power strong outflows and jets, then spin of the black hole, perhaps, might be the fundamental parameter to account for the observed astrophysical processes in an accretion powered system.

by Dotti, M. and Volonteri, M. and Perego, A. and Colpi, M. and Ruszkowski, M. and Haardt, F.
11 pages, 3 figures. Accepted for publication in MNRAS

Using high resolution hydrodynamical simulations, we explore the spin evolution of massive dual black holes orbiting inside a circumnuclear disc, relic of a gas-rich galaxy merger. The black holes spiral inwards from initially eccentric co or counter-rotating coplanar orbits relative to the disc’s rotation, and accrete gas that is carrying a net angular momentum. As the black hole mass grows, its spin changes in strength and direction due to its gravito-magnetic coupling with the small-scale accretion disc. We find that the black hole spins loose memory of their initial orientation, as accretion torques suffice to align the spins with the angular momentum of their orbit on a short timescale (<1-2 Myr). A residual off-set in the spin direction relative to the orbital angular momentum remains, at the level of <10 degrees for the case of a cold disc, and <30 degrees for a warmer disc. Alignment in a cooler disc is more effective due to the higher coherence of the accretion flow near each black hole that reflects the large-scale coherence of the disc’s rotation. If the massive black holes coalesce preserving the spin directions set after formation of a Keplerian binary, the relic black hole resulting from their coalescence receives a relatively small gravitational recoil. The distribution of recoil velocities inferred from a simulated sample of massive black hole binaries has median <70 km/s much smaller than the median resulting from an isotropic distribution of spins.

by Fragile, P. Chris
5 pages, 4 figures, accepted to ApJ Letters

One of the primary means of determining the spin of an astrophysical black hole is by actually measuring the inner radius of a surrounding accretion disk and using that to infer the spin. By comparing a number of different estimates of the inner radius from simulations of tilted accretion disks with differing black-hole spins, we show that such a procedure can give quite wrong answers. Over the range 0 <= a/M <= 0.9, we find that, for moderately thick disks (H/r ~ 0.2) with modest tilt (15 degrees), the inner radius is nearly independent of spin. This result is likely dependent on tilt, such that for larger tilts, it may even be that the inner radius would increase with increasing spin. In the opposite limit, we confirm through numerical simulations of untilted disks that, in the limit of zero tilt, the inner radius recovers approximately the expected dependence on spin.

by Anderson, Matthew and Lehner, Luis and Megevand, Miguel and Neilsen, David
9 pages, 6 figures

We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. To this end, we study radiation transfer from the system incorporating bremsstrahlung emission from a Maxwellian plasma and absorption given by Kramer’s opacity law modified to incorporate blackbody effects. We employ this model in the radiation transfer integration to compute the luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations from a transparent limit (in which the emissivity is integrated over the computational domain and over all frequencies) and with a simple thermal emission model. We find close agreement between the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for electromagnetic signals above $latex 10^{14}$ Hz. For lower frequencies, we find a self-eclipsing behavior in the disk, resulting in a strong intensity variability connected to the orbital period of the disk.

by Rajesh, S. R. and Mukhopadhyay, Banibrata
25 pages including 22 figures; to appear in MNRAS

We investigate the viscous two temperature accretion discs around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes selfconsistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms, for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter 0.998. It is found that the flow, during its infall from the Keplerian to sub-Keplerian transition region to the black hole event horizon, passes through various phases of advection — general advective paradigm to radiatively inefficient phase and vice versa. Hence the flow governs much lower electron temperature ~10^8-10^{9.5} K, in the range of accretion rate in Eddington units 0.01 <~ \mdot <~ 100, compared to the hot protons of temperature ~ 10^{10.2} – 10^{11.8}K. Therefore, the solution may potentially explain the hard X-rays and \gamma-rays emitted from AGNs and X-ray binaries. We then show that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the under-fed AGNs and quasars (e.g. Sgr A^*) to ultra-luminous X-ray sources (e.g. SS433), at different combinations of input parameters such as mass accretion rate, ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67).

by Li, Guang-Xing and Yuan, Ye-Fei and Cao, Xinwu
15 pages, 13 figures

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are investigated in this work. The resulting emergent spectra are calculated using relativistic ray-tracing method. The angular dependence of the disk luminosity, the effects of the heat advection and the effect of disk thickness on the estimation of the black hole spin and accretion rate are discussed. The improvements compared to previous works are that we use self-consistent disk equations and we consider the disk self-shadowing effect. We find that at moderate accretion rate, with inclusion of the heat advection effect, radiation trapped in the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to the emergent spectra. At high accretion rate, large inclination and large black hole spin, both the disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these effects will influence the measurement of the black hole spin based on the spectral fitting and influence the angular dependence of luminosity.

by Corrales, Lia R. and Haiman, Zoltán and MacFadyen, Andrew
16 pages with 14 figures, submitted to MNRAS

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern.

by Rossi, Elena M. and Lodato, G. and Armitage, P. J. and Pringle, J. E. and King, A. R.
16 pages, accepted by MNRAS. Animations of the simulations are available at http://jilawww.colorado.edu/~pja/recoil.html

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.

by Bogdanovic, Tamara and Eracleous, Michael and Sigurdsson, Steinn
8 pages, 8 figures, submitted to the New Astronomy Reviews as a part of the SCSLSA-7 proceedings

Detection of electromagnetic (EM) counterparts of pre-coalescence binaries has very important implications for our understanding of the evolution of these systems as well as the associated accretion physics. In addition, a combination of EM and gravitational wave signatures observed from coalescing supermassive black hole binaries (SBHBs) would provide independent measurements of redshift and luminosity distance, thus allowing for high precision cosmological measurements. However, a statistically significant sample of these objects is yet to be attained and finding them observationally has proven to be a difficult task. Here we discuss existing observational evidence and how further advancements in the theoretical understanding of observational signatures of SBHBs before and after the coalescence can help in future searches.

by Loeb, Abraham
3 pages

Mergers of gas-rich galaxies lead to black hole binaries that coalesce as a result of dynamical friction on the ambient gas. Once the binary tightens to <10^3 Schwarzschild radii, its merger is driven by the emission of gravitational waves (GWs). We show that this transition occurs generically at orbital periods of ~1-10 years and an orbital velocity V of a few thousand km/s, with a very weak dependence on the supply rate of gas (V proportional to Mdot^{1/8}). Therefore, as binaries enter their GW-dominated inspiral, they inevitably induce large periodic shifts in the broad emission lines of any associated quasar(s). The probability of finding a binary in tighter configurations scales as V^{-8} owing to their much shorter lifetimes. Systematic monitoring of the broad emission lines of quasars on timescales of months to decades can set a lower limit on the expected rate of GW sources for LISA.

by Alig, Christian and Burkert, Andreas and Johansson, Peter H. and Schartmann, Marc
10 pages, 6 figures, submitted to MNRAS

We simulate clouds in the Galactic Centre (GC) crossing over the black hole in parts and present this as a possible formation mechanism for the observed stellar disks in the GC through the redistribution of angular momentum by colliding material with opposite angular momentum. A parameter study using six high resolution simulations of an isothermal cloud of constant density falling onto the black hole and crossing over it in parts demonstrates that this mechanism is able to reproduce the observed disk properties in the GC. The evolution of the ensuing accretion disks is highly non-linear with the redistribution of the angular momentum through dissipative processes being a dominant effect. We analyse the resulting Toomre unstable, eccentric gaseous disk and show that this already yields a good comparison with the observed stellar disk size and eccentricity in the GC. The best simulation results in an outer radius of 1 pc, a mass of 10$latex ^4$ M$latex _{\sun}$ and an eccentricity of 0.24 for the Toomre unstable disk, which compares well with the observations.

by van Meter, James R. and Wise, John H. and Miller, M. Coleman and Reynolds, Christopher S. and Centrella, Joan M. and Baker, John G. and Boggs, William D. and Kelly, Bernard J. and McWilliams, Sean T.
5 pages, 4 figures

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

by Shapiro, Charles and Bacon, David and Hendry, Martin and Hoyle, Ben
9 pages, 4 figures, submitted to MNRAS

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.

by Garofalo, David
Phys.RevD accepted

We study the spin dependence of accretion onto rotating Kerr black holes using analytic techniques. In its linear regime, angular momentum transport in MHD turbulent accretion flow involves the generation of radial magnetic field connecting plasma in a differentially rotating flow. We take a first principles approach, highlighting the constraint that limits the generation and amplification of radial magnetic fields, stemming from the transfer of energy from mechanical to magnetic form. Because the energy transferred in magnetic form is ultimately constrained by gravitational potential energy or Killing energy, the spin-dependence of the latter allows us to derive spin-dependent constraints on the success of the accreting plasma to expel its angular momentum and accrete. We find an inverse relationship between this ability and black hole spin. If this radial magnetic field generation forms the basis for angular momentum transfer in accretion flows, accretion rates involving Kerr black holes are expected to be lower as the black hole spin increases in the prograde sense.

by Dai, Lixin and Fuerst, Steven V. and Blandford, Roger
16 pages, 11 figures, submitted to MNRAS; corrected typos

We present simulated results of quasi-periodic flares generated by the inelastic collisions of a star bound to a super-massive black hole (SMBH) and its attendant accretion disc. We show that the behavior of the quasi-periodicity is affected by the mass and spin of the black hole and the orbital elements of the stellar orbit. We also evaluate the possibility of extracting useful information on these parameters and verifying the character of the Kerr metric from such quasi-periodic signals. Comparisons are made with the observed optical outbursts of OJ287, infrared flares from the Galactic center and X-ray variability in RE J1034+396.

by Chang, P. and Strubbe, L. E. and Menou, K. and Quataert, E.
11 pages, 6 figures, submitted to MNRAS

Using a one-dimensional height integrated model, we calculate the evolution of an unequal mass binary black hole with a coplanar gas disk that contains a gap due to the presence of the secondary black hole. Viscous evolution of the outer circumbinary disk initially hardens the binary, while the inner disk drains onto the primary (central) black hole. As long as the inner disk remains cool and thin at low $latex \dot{M}_{\rm ext}$ (rather than becoming hot and geometrically thick), the mass of the inner disk reaches an asymptotic mass typically $latex \sim 10^{-3} – 10^{-4} M_{\odot}$. Once the semimajor axis shrinks below a critical value, angular momentum losses from gravitational waves dominate over viscous transport in hardening the binary. The inner disk then no longer responds viscously to the inspiraling black holes. Instead, tidal interactions with the secondary rapidly drive the inner disk into the primary. Tidal and viscous dissipation in the inner disk lead to a late time brightening in luminosity $latex L\propto t_{\rm minus}^{5/4}$, where $latex t_{\rm minus}$ is the time prior to the final merger. This late time brightening peaks $latex \sim 1$ day prior to the final merger at $latex \sim 0.1 L_{\rm Edd}$. This behavior is relatively robust because of self regulation in the coupled viscous-gravitational evolution of such binary systems. It constitutes a unique electromagnetic signature of a binary supermassive black hole merger and may allow the host galaxy to be identified if used in conjunction with the Laser Interferometric Space Antenna (LISA) localization.

by van Putten, Maurice H. P. M.
Invited review, in Orleans School On The Sciences Of The Universe: Mass Motion (Springer-Verlag), to appear

Ultra-high energy cosmic rays (UHECRs) and gamma-ray bursts (GRBs) are the most exceptional nonthermal transient events, that appear to be associated with black holes. Here, we describe radiation mechanisms induced by turbulent flows around rapidly rotating black holes: high-energy emissions from a relativistic capillary effect along the black hole spin-axis and low-energy emissions by catalytic conversion of spin-energy. High-energy emissions arise, concurrently, in photons and, upstream of an outgoing Alfv\’en front, in ionic contaminants by linear acceleration. The latter develop into ultra-high energy cosmic rays (UHECRs) about the Greisen-Zatsepin-Kuzmin (GZK) threshold in low-luminosity, intermittent active galactic nuclei. These may include Seyfert galaxies and Cen A suggested by detections of UHECRs by the Pierre Auger Observatory and, for the latter, also of Very High Energy (VHE) gamma-rays by the High Energy Stereoscopic System (HESS). Nearly complete spin-down of stellar mass black holes is common to collapsars and mergers of neutron stars with another neutron star or companion black hole. Thus, long GRBs from rotating black holes explain events with and without supernovae and a diversity in their X-ray afterglows. Their intrinsic exponential decay is remarkably consistent with the average of 600 light curves of long GRBs, whose total output agrees with observed peak and true energies in gamma-rays. We conclude that long GRBs are spin-powered. Gravitational radiation from turbulent flows in SgrA* might be of interest to the planned Laser Interferometric Space Antenna (LISA) and, for stellar mass black holes in GRBs, should be detectable by LIGO-Virgo. Long GRBs from naked inner engines produced in mergers produce long-duration radio-burst that may be seen in all-sky surveys by the Low Frequency Array (LOFAR).

by Palenzuela, Carlos and Anderson, Matthew and Lehner, Luis and Liebling, Steven L. and Neilsen, David
4 pages, 4 figures

In addition to producing gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary’s dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

by Prince, T. A. and Team, for the LISA International Science
Science White Paper submitted to the Astro2010 Decadal Survey

This Astro2010 science white paper provides an overview of the opportunities in low-frequency gravitational-wave astronomy, a new field that is poised to make significant advances. While discussing the broad context of gravitational-wave astronomy, this paper concentrates on the low-frequency region (10^(-5) to 1 Hz), a frequency range abundantly populated in strong sources of gravitational waves including massive black hole mergers, ultra-compact stellar-mass galactic binaries, and capture of compact objects by massive black holes in the nuclei of galaxies.

by Phinney, E. S.
Science White Paper submitted to the Astro2010 Decadal Survey

The principal goal of this whitepaper is not so much to demonstrate that gravitational wave detectors like LIGO and LISA will help answer many central questions in astronomy and astrophysics, but to make the case that they can help answer a far greater range of questions if we prepare to make the (sometimes substantial) effort to identify electromagnetic counterparts to the gravitational wave sources.

by Bloom, Joshua S. and Holz, Daniel E. and Hughes, Scott A. and Menou, Kristen and Adams, Allan and Anderson, Scott F. and Becker, Andy and Bower, Geoffrey C. and Brandt, Niel and Cobb, Bethany and Cook, Kem and Corsi, Alessandra and Covino, Stefano and Fox, Derek and Fruchter, Andrew and Fryer, Chris and Grindlay, Jonathan and Hartmann, Dieter and Haiman, Zoltan and Kocsis, Bence and Jones, Lynne and Loeb, Abraham and Marka, Szabolcs and Metzger, Brian and Nakar, Ehud and Nissanke, Samaya and Perley, Daniel A. and Piran, Tsvi and Poznanski, Dovi and Prince, Tom and Schnittman, Jeremy and Soderberg, Alicia and Strauss, Michael and Shawhan, Peter S. and Shoemaker, David H. and Sievers, Jonathan and Stubbs, Christopher and Tagliaferri, Gianpiero and Ubertini, Pietro and Wozniak, Przemyslaw
7 pages (plus one coverpage), submitted to the US Astro2010 Decadal Survey. This is a living document, with updates expected to be posted to this archive. Those interested in contributing should contact J. S. Bloom

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?