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 Haiman, Zoltán
10 pages, 5 figures. Review contribution to the Proceedings of “The First Stars and Galaxies: Challenges for the Next Decade”, Austin, TX, March 8-11, 2010

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may arise by the combination of Eddington-limited growth and mergers of stellar-mass seed BHs left behind by the first generation of metal-free stars, or by the rapid direct collapse of gas in rare special environments where the gas can avoid fragmenting into stars. In this contribution, I review these two competing scenarios. I also briefly mention some more exotic ideas and how the different models may be distinguished in the future by LISA and other instruments.

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 Graham, Alister W. and Onken, Christopher A. and Athanassoula, E. and Combes, Francoise
17 pages. Submitted to MNRAS on July 19

[Abridged] We present an updated and improved M_bh-sigma diagram containing 64 galaxies for which M_bh measurements (not just upper limits) are available. Due to new and increased black hole masses at the high-mass end, and a better representation of barred galaxies at the low-mass end, the “classical” (all morphological type) M_bh-sigma relation for predicting black hole masses is log(M_bh/M_Sun) = (5.13+/-0.34)log[sigma/200] + (8.13+/-0.05), with an rms scatter of 0.43 dex. Modifying the regression analysis to correct for a hitherto over-looked sample bias in which black holes with masses <10^6 M_Sun are not (yet) detectable, the relation steepens further to give log(M_bh/M_Sun) = (5.95+/-0.44)log[sigma/200] + (8.15+/-0.06). We have also updated the “barless'' and “elliptical-only'' M_bh-sigma relations introduced by Graham and Hu in 2008 due to the offset nature of barred/disc galaxies. These relations have a total scatter as low as 0.34 dex and currently define the upper envelope of points in the M_bh-sigma diagram. These relations also have a slope consistent with the value 5, in agreement with the prediction by Silk & Rees based on feedback from massive black holes in bulges built by monolithic-collapse.

Using updated virial products and velocity dispersions from 28 active galactic nuclei, we determine that the optimal scaling factor f – which brings their virial products in line with the 64 directly measured black hole masses – is 2.8^{+0.7}_{-0.5}. This is roughly half the value reported by Onken et al. and Woo et al., and consequently halves the mass estimates of most high-redshift quasars. We have explored the results after separating the samples into barred and non-barred galaxies, and we have also developed a preliminary corrective term to the velocity dispersion based on bar dynamics.

by Decarli, Roberto and Dotti, Massimo and Montuori, Carmen and Liimets, Tiina and Ederoclite, Alessandro
2 figures, 1 table. Accepted for publication in ApJ Letters

We report the discovery of peculiar features in the optical spectrum of 4C+22.25, a flat spectrum radio quasar at z=0.4183 observed in the SDSS and in a dedicated spectroscopic follow-up from the Nordic Optical Telescope. The Hbeta and Halpha lines show broad profiles (FWHM~12,000 km/s), faint fluxes and extreme offsets (Delta v=8,700+/-1,300 km/s) with respect to the narrow emission lines. These features show no significant variation in a time lag of ~3.1 yr (rest frame). We rule out possible interpretations based on the superposition of two sources or on recoiling black holes, and we discuss the virtues and limitations of a massive black hole binary scenario.

by Greene, J. E. and Peng, C. Y. and Kim, M. and Kuo, C. Y. and Braatz, J. A. and Impellizzeri, C. M. V. and Condon, J. J. and Lo, K. Y. and Henkel, C. and Reid, M. J.
21 pages, 14 figures, accepted for publication in the Astrophysical Journal

The black hole (BH)-bulge correlations have greatly influenced the last decade of effort to understand galaxy evolution. Current knowledge of these correlations is limited predominantly to high BH masses (M_BH> 10^8 M_sun) that can be measured using direct stellar, gas, and maser kinematics. These objects, however, do not represent the demographics of more typical L< L* galaxies. This study transcends prior limitations to probe BHs that are an order of magnitude lower in mass, using BH mass measurements derived from the dynamics of H_2O megamasers in circumnuclear disks. The masers trace the Keplerian rotation of circumnuclear molecular disks starting at radii of a few tenths of a pc from the central BH. Modeling of the rotation curves, presented by Kuo et al. (2010), yields BH masses with exquisite precision. We present stellar velocity dispersion measurements for a sample of nine megamaser disk galaxies based on long-slit observations using the B&C spectrograph on the Dupont telescope and the DIS spectrograph on the 3.5m telescope at Apache Point. We also perform bulge-to-disk decomposition of a subset of five of these galaxies with SDSS imaging. The maser galaxies as a group fall below the M_BH-sigma* relation defined by elliptical galaxies. We show, now with very precise BH mass measurements, that the low-scatter power-law relation between M_BH and sigma* seen in elliptical galaxies is not universal. The elliptical galaxy M_BH-sigma* relation cannot be used to derive the BH mass function at low mass or the zeropoint for active BH masses. The processes (perhaps BH self-regulation or minor merging) that operate at higher mass have not effectively established an M_BH-sigma* relation in this low-mass regime.

by Garfinkle, David and Pretorius, Frans and Yunes, Nicolas
4 pages, no figures, submitted to Rapid Communications

We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice-versa.

by Key, Joey Shapiro and Cornish, Neil J.
11 pages, 19 figures

The Laser Interferometer Space Antenna (LISA) is designed to detect gravitational wave signals from astrophysical sources, including those from coalescing binary systems of compact objects such as black holes. Colliding galaxies have central black holes that sink to the center of the merged galaxy and begin to orbit one another and emit gravitational waves. Some galaxy evolution models predict that the binary black hole system will enter the LISA band with significant orbital eccentricity, while other models suggest that the orbits will already have circularized. Using a full seventeen parameter waveform model that includes the effects of orbital eccentricity, spin precession and higher harmonics, we investigate how well the source parameters can be inferred from simulated LISA data. Defining the reference eccentricity as the value one year before merger, we find that for typical LISA sources, it will be possible to measure the eccentricity to an accuracy of parts in a thousand. The accuracy with which the eccentricity can be measured depends only very weakly on the eccentricity, making it possible to distinguish circular orbits from those with very small eccentricities. LISA measurements of the orbital eccentricity can provide strong constraints on theories of galaxy mergers in the early universe.

by Kelly, Brandon C. and Vestergaard, Marianne and Fan, Xiaohui and Hopkins, Philip and Hernquist, Lars and Siemiginowska, Aneta
Accepted by ApJ, 25 pages (emulateapj), 11 figures

We present an estimate of the black hole mass function (BHMF) of broad line quasars (BLQSOs) that self-consistently corrects for incompleteness and the statistical uncertainty in the mass estimates, based on a sample of 9886 quasars at 1 < z 1 it is highly incomplete at M_BH < 10^9 M_Sun and L / L_Edd 1, where the BLQSO phase occurs at the end of a fueling event when black hole feedback unbinds the accreting gas, halting the accretion flow.

by Park, KwangHo and Ricotti, Massimo
14 pages, 12 figures, submitted to ApJ

We study the effect of radiative feedback on accretion onto intermediate mass black holes (IMBHs) using the hydrodynamical code ZEUS-MP with a radiative transfer algorithm. In this paper, the first of a series, we assume accretion from a uniformly dense gas with zero angular momentum. Our 1D and 2D simulations explore how X-ray and UV radiation emitted near the black hole regulates the gas supply from large scales. Both 1D and 2D simulations show similar accretion rate and period between peaks in accretion, meaning that the hydro-instabilities that develop in 2D simulations do not affect the mean flow properties. We present a suite of simulations exploring accretion across a large parameter space, including different radiative efficiencies and radiation spectra, black hole masses, density and temperature, $latex T_\infty$, of the neighboring gas. In agreement with previous studies we find regular oscillatory behavior of the accretion rate, with duty cycle $latex \sim 7%$, mean accretion rate 3-6% $latex (T_{\infty}/10^4 {\rm K})^{2.5}$ of the Bondi rate and peak accretion $latex \sim 10$ times the mean. We derive parametric formulas for the period between bursts, the mean accretion rate and the peak luminosity of the bursts and thus provide a formulation of how feedback regulated accretion operates. The temperature profile of the hot ionized gas is crucial in determining the accretion rate, while the period of the bursts is proportional to the mean size of the Str\”{o}mgren sphere. We also find that softer spectrum of radiation produces higher accretion rate. This study is a first step to model the growth of seed black holes in the early universe and to make a prediction of the number and the luminosity of ultra-luminous X-ray sources in galaxies produced by IMBHs accreting from the interstellar medium.

by Schneider, Justus and Amaro-Seoane, Pau and Spurzem, Rainer
Submitted to MNRAS, comments welcome

Dense stellar systems such as globular clusters, galactic nuclei and nuclear star clusters are ideal loci to study stellar dynamics due to the very high densities reached, usually a million times higher than in the solar neighborhood; they are unique laboratories to study processes related to relaxation. There are a number of different techniques to model the global evolution of such a system. In statistical models we assume that relaxation is the result of a large number of two-body gravitational encounters with a net local effect. We present two moment models that are based on the collisional Boltzmann equation. By taking moments of the Boltzmann equation one obtains an infinite set of differential moment equations where the equation for the moment of order $latex n$ contains moments of order $latex n+1$. In our models we assume spherical symmetry but we do not require dynamical equilibrium. We truncate the infinite set of moment equations at order $latex n=4$ for the first model and at order $latex n=5$ for the second model. The collisional terms on the right-hand side of the moment equations account for two-body relaxation and are computed by means of the Rosenbluth potentials. We complete the set of moment equations with closure relations which constrain the degree of anisotropy of our model by expressing moments of order $latex n+1$ by moments of order $latex n$. The accuracy of this approach relies on the number of moments included from the infinite series. Since both models include fourth order moments we can study mechanisms in more detail that increase or decrease the number of high velocity stars. The resulting model allows us to derive a velocity distribution function, with unprecedented accuracy, compared to previous moment models.

by Robinson, Andrew and Young, Stuart and Axon, David J. and Kharb, Preeti and Smith, James E.
18 pages, 4 figures, accepted for publication in the Astrophysical Journal as a Letter

We report spectropolarimetric observations of the quasar E1821+643 (z=0.297), which suggest that it may be an example of gravitational recoil due to anisotropic emission of gravitational waves following the merger of a supermassive black hole (SMBH) binary. In total flux, the broad Balmer lines are redshifted by ~1000 km/s relative to the narrow lines and have highly red asymmetric profiles, whereas in polarized flux the broad H_alpha line exhibits a blueshift of similar magnitude and a strong blue asymmetry. We show that these observations are consistent with a scattering model in which the broad-line region has two components, moving with different bulk velocities away from the observer and towards a scattering region at rest in the host galaxy. If the high velocity system is identified as gas bound to the SMBH, this implies that the SMBH is itself moving with a velocity ~2100 km/s relative to the host galaxy. We discuss some implications of the recoil hypothesis and also briefly consider whether our observations can be explained in terms of scattering of broad-line emission originating from the active component of an SMBH binary, or from an outflowing wind.

by Sesana, A.
15 pages, 9 figures, accepted for publication in the Astrophysical Journal

We construct evolutionary tracks for massive black hole binaries (MBHBs) embedded in a surrounding distribution of stars. The dynamics of the binary is evolved by taking into account the erosion of the central stellar cusp bound to the massive black holes, the scattering of unbound stars feeding the binary loss cone, and the emission of gravitational waves (GWs). Stellar dynamics is treated in a hybrid fashion by coupling the results of numerical 3-body scattering experiments of bound and unbound stars to an analytical framework for the evolution of the stellar density distribution and for the efficiency of the binary loss cone refilling. Our main focus is on the behaviour of the binary eccentricity, in the attempt of addressing its importance in the merger process and its possible impact for GW detection with the planned Laser Interferometer Space Antenna ({\it LISA}), and ongoing and forthcoming pulsar timing array (PTA) campaigns. We produce a family of evolutionary tracks extensively sampling the relevant parameters of the system which are the binary mass, mass ratio and initial eccentricity, the slope of the stellar density distribution, its normalization and the efficiency of loss cone refilling. We find that, in general, stellar dynamics causes a dramatic increase of the MBHB eccentricity, especially for initially already mildly eccentric and/or unequal mass binaries. When applied to standard MBHB population models, our results predict eccentricities in the ranges $latex 10^{-3}-0.2$ and $latex 0.03-0.3$ for sources detectable by {\it LISA} and PTA respectively. Such figures may have a significant impact on the signal modelling, on source detection, and on the development of parameter estimation algorithms.

by Porto, Rafael A.
25 pages, 4 figures, revtex4

Using effective field theory (EFT) techniques we calculate the next-to-leading order (NLO) spin-orbit contributions to the gravitational potential of inspiralling compact binaries. We use the covariant spin supplementarity condition (SSC), and explicitly prove the equivalence with previous results by Faye et al. in arXiv:gr-qc/0605139. We also show that the direct application of the Newton-Wigner SSC at the level of the action leads to the correct dynamics using a canonical (Dirac) algebra. This paper then completes the calculation of the necessary spin dynamics within the EFT formalism that will be used in a separate paper to compute the spin contributions to the energy flux and phase evolution to NLO.

by Schneider, Raffaella and Marassi, Stefania and Ferrari, Valeria
10 pages, 9 figures, proceedings of the GWDAW 10 Conference, submitted to Class. & Quantum Grav

Astrophysical sources emit gravitational waves in a large variety of processes occurred since the beginning of star and galaxy formation. These waves permeate our high redshift Universe, and form a background which is the result of the superposition of different components, each associated to a specific astrophysical process. Each component has different spectral properties and features that it is important to investigate in view of a possible, future detection. In this contribution, we will review recent theoretical predictions for backgrounds produced by extragalactic sources and discuss their detectability with current and future gravitational wave observatories.

by Sethi, Shiv K. and Haiman, Zoltán and Pandey, Kanhaiya
submitted to ApJ, 5 emulateapj pages and 5 figures

It has been proposed that primordial gas in early dark matter halos, with virial temperatures above 10^4 K, can avoid fragmentation and undergo rapid collapse, possibly resulting in a supermassive black hole (SMBH). This requires the gas to avoid cooling and to remain at temperatures near T=10^4 K. We show that this condition can be satisfied in the presence of a sufficiently strong primordial magnetic field, which heats the collapsing gas via ambipolar diffusion. If the field has a strength above B = 3.6 (comoving) nG, the collapsing gas is kept warm (T=10^4K) until it reaches the critical density n_crit=10^3 cm^{-3} at which the roto-vibrational states of H_2 approach local thermodynamic equilibrium. H_2-cooling then remains inefficient, and the gas temperature stays near 10^4K, even as it continues to collapse to higher densities. The critical magnetic field strength required to permanently suppress H_2-cooling is somewhat higher than upper limit of approx. 2 nG from the cosmic microwave background (CMB). However, it can be realized in the rare (2-3)-sigma regions of the spatially fluctuating B-field; these regions contain a sufficient number of halos to account for the z=6 quasar BHs.

by Gergely, László Á. and Biermann, Peter L. and Caramete, Laurenţiu I.
11 pages, 2 figures

During post-Newtonian evolution of a compact binary, a mass ratio different from one provides a second small parameter, which can lead to unexpected results. We present a statistics of supermassive black hole candidates, which enables us first to derive their mass distribution, then to establish a logarithmically even probability of the mass ratios at their encounter. In the mass ratio range (1/30,1/3) of supermassive black hole mergers representing 40% of all possible cases, the combined effect of spin-orbit precession and gravitational radiation leads to a spin-flip of the dominant spin during the inspiral phase of the merger. This provides a mechanism for explaining a large set of observations on X-shaped radio galaxies. In another 40%, with mass ratios (1/30,1/1000) a spin-flip never happens, while in the remaining 20% of mergers with mass ratios (1/3,1) it may occur during the plunge. We analyze the magnitude of the spin-flip angle occurring during the inspiral as function of the mass ratio and original relative orientation of the spin and orbital angular momentum. We also derive a formula for the final spin at the end of the inspiral in this mass ratio range.

by Batcheldor, D. and Robinson, A. and Axon, D. J. and Perlman, E. S. and Merritt, D.
ApJ Letters accepted

Isophotal analysis of M87, using data from the Advanced Camera for Surveys, reveals a projected displacement of 6.8 +/- 0.8 pc (~ 0.1 arcsec) between the nuclear point source (presumed to be the location of the supermassive black hole, SMBH) and the photo-center of the galaxy. The displacement is along a position angle of 307 +/- 17 degrees and is consistent with the jet axis. This suggests the active SMBH in M87 does not currently reside at the galaxy center of mass, but is displaced in the counter-jet direction. Possible explanations for the displacement include orbital motion of an SMBH binary, gravitational perturbations due to massive objects (e.g., globular clusters), acceleration by an asymmetric or intrinsically one-sided jet, and gravitational recoil resulting from the coalescence of an SMBH binary. The displacement direction favors the latter two mechanisms. However, jet asymmetry is only viable, at the observed accretion rate, for a jet age of >0.1 Gyr and if the galaxy restoring force is negligible. This could be the case in the low density core of M87. A moderate recoil ~1 Myr ago might explain the disturbed nature of the nuclear gas disk, could be aligned with the jet axis, and can produce the observed offset. Alternatively, the displacement could be due to residual oscillations resulting from a large recoil that occurred in the aftermath of a major merger any time in the last 10 Gyr.

by Klein, Antoine and Jetzer, Philippe
7 pages, 1 figure; Accepted for publication in Phys. Rev. D

We compute here the spin-orbit and spin-spin couplings needed for an accurate computation of the phasing of gravitational waves emitted by comparable-mass binaries on eccentric orbits at the second post-Newtonian (PN) order. We use a quasi-Keplerian parametrization of the orbit free of divergencies in the zero eccentricity limit. We find that spin-spin couplings induce a residual eccentricity for coalescing binaries at 2PN, of the order of $latex 10^{-4}$-$latex 10^{-3}$ for supermassive black hole binaries in the LISA band. Spin-orbit precession also induces a non-trivial pattern in the evolution of the eccentricity, which could help to reduce the errors on the determination of the eccentricity and spins in a gravitational wave measurement.

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 Hirata, Christopher M. and Holz, Daniel E. and Cutler, Curt
11 pages, 6 figures, to be submitted to PRD

Gravitational wave sources are a promising cosmological standard candle because their intrinsic luminosities are determined by fundamental physics (and are insensitive to dust extinction). They are, however, affected by weak lensing magnification due to the gravitational lensing from structures along the line of sight. This lensing is a source of uncertainty in the distance determination, even in the limit of perfect standard candle measurements. It is commonly believed that the uncertainty in the distance to an ensemble of gravitational wave sources is limited by the standard deviation of the lensing magnification distribution divided by the square root of the number of sources. Here we show that by exploiting the non-Gaussian nature of the lensing magnification distribution, we can improve this distance determination, typically by a factor of 2–3; we provide a fitting formula for the effective distance accuracy as a function of redshift for sources where the lensing noise dominates.

by Bekki, Kenji and Graham, Alister W.
15 page, 5 figures, accepted in ApJL

Giant elliptical galaxies, believed to be built from the merger of lesser galaxies, are known to house a massive black hole at their center rather than a compact star cluster. If low- and intermediate-mass galaxies do indeed partake in the hierarchical merger scenario, then one needs to explain why their dense nuclear star clusters are not preserved in merger events. A valuable clue may the recent revelation that nuclear star clusters and massive black holes frequently co-exist in intermediate mass bulges and elliptical galaxies. In an effort to understand the physical mechanism responsible for the disappearance of nuclear star clusters, we have numerically investigated the evolution of merging star clusters with seed black holes. Using black holes that are 1-5% of their host nuclear cluster mass, we reveal how their binary coalescence during a merger dynamically heats the newly wed star cluster, expanding it, significantly lowering its central stellar density, and thus making it susceptible to tidal destruction during galaxy merging. Moreover, this mechanism provides a pathway to explain the observed reduction in the nucleus-to-galaxy stellar mass ratio as one proceeds from dwarf to giant elliptical galaxies.

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 Silk, Joe and Nusser, Adi

Force balance considerations put a limit on the rate of AGN radiation momentum output, $latex L/c$, capable of driving galactic superwinds. We show that this condition is insufficient: black holes obeying the observed $latex \mbh -\sigma $ relation cannot supply enough energy in radiation which can drive the gas out by pressure alone. The shortfall is by up to an order of magnitude in most, but not all, cases. We propose that outflow-triggering of star formation by enhancing the intercloud medium turbulent pressure and squeezing clouds can supply the necessary boost, and suggest possible tests of this hypothesis. We further point out that the time-scales for Bondi accretion and for orbital decay of merging clumps by dynamical friction in the nuclear disk around a central black hole both follow a similar scaling with mass, favoring the most massive black holes, but the latter process is up to two orders of magnitude more rapid at $latex z\gtsim 10.$ The combination of accretion and coalescence results in earlier formation of more massive black holes, and, in particular, can account for the masses of the black holes inferred to power AGN at $latex z\sim 6.$

by Burkert, Andreas and Tremaine, Scott
13 pages, 3 figures, submitted to ApJ

Elliptical, lenticular, and early-type spiral galaxies show a remarkably tight power-law correlation between the mass M_BH of their central supermassive black hole (SMBH) and the number N_GC of globular clusters: M_BH=m*N_GC^(1.11+/-0.04) with m=1.3*10^5 solar masses. Thus, to a good approximation the SMBH mass is the same as the total mass of the globular clusters. Based on a limited sample of 13 galaxies, this relation appears to be a better predictor of SMBH mass (rms scatter 0.2 dex) than the M_BH-sigma relation between SMBH mass and velocity dispersion sigma. The small scatter reflects the fact that galaxies with high globular cluster specific frequency S_N tend to harbor SMBHs that are more massive than expected from the M_BH-sigma relation. A possible explanation is that both large black-hole masses and large globular cluster populations are associated with recent major mergers.

by Kesden, Michael and Sperhake, Ulrich and Berti, Emanuele
7 pages, 4 figures, submitted to ApJL

Numerical-relativity simulations indicate that the black hole produced in a binary merger can recoil with a velocity up to v_max ~ 4,000 km/s with respect to the center of mass of the initial binary. This challenges the paradigm that most galaxies form through hierarchical mergers, yet retain supermassive black holes at their centers despite having escape velocities much less than v_max. Interaction with a circumbinary disk can align the binary black hole spins with their orbital angular momentum, reducing the recoil velocity of the final black hole produced in the subsequent merger. However, the effectiveness of this alignment depends on highly uncertain accretion flows near the binary black holes. In this Letter, we show that if the spin S_1 of the more massive binary black hole is even partially aligned with the orbital angular momentum L, relativistic spin precession on sub-parsec scales can align the binary black hole spins with each other. This alignment significantly reduces the recoil velocity even in the absence of gas. For example, if the angle between S_1 and L at large separations is 10 degrees while the second spin S_2 is isotropically distributed, the spin alignment discussed in this paper reduces the median recoil from 864 km/s to 273 km/s for maximally spinning black holes with a mass ratio of 9/11. This reduction will greatly increase the fraction of galaxies retaining their supermassive black holes.

by Krumholz, Mark R. and Burkert, Andreas
14 pages, 3 figures, emulateapj format, submitted to ApJ

We present a first-principles derivation of the evolution equations describing a thin axisymmetric disk of gas and stars with an arbitrary rotation curve that is kept in a state of marginal gravitational instability and energy equilibrium due to the balance between energy released by accretion and energy lost due to decay of turbulence. Unlike previous analyses of this problem, our results do not depend on an assumed model for the rate of mass and angular momentum transport due to gravitational instability, or on an order-of-magnitude energy equilibrium argument. Instead, we self-consistently determine the position- and time-dependent transport rates from the fluid dynamical equations. We show that there is a steady-state configuration for disks dominated by gravitational instability, and for disks in this state we analytically determine the velocity dispersion, surface density, and rates of mass and angular momentum transport as a function of the gas mass fraction, the rotation curve, and the rate of external accretion onto the disk edge. We show that disks that are initially out of steady state will evolve into it on timescales comparable to the orbital period if the accretion rate is high. Finally, we discuss the implications of these results for the structure of disks in a broad range of environments, including high redshift galaxies, the outer gaseous disks of local galaxies, and accretion disks around protostars.

by Volonteri, Marta
To appear in The Astronomy and Astrophysics Review. The final publication is available at http://www.springerlink.com

Evidence shows that massive black holes reside in most local galaxies. Studies have also established a number of relations between the MBH mass and properties of the host galaxy such as bulge mass and velocity dispersion. These results suggest that central MBHs, while much less massive than the host (~ 0.1%), are linked to the evolution of galactic structure. In hierarchical cosmologies, a single big galaxy today can be traced back to the stage when it was split up in hundreds of smaller components. Did MBH seeds form with the same efficiency in small proto-galaxies, or did their formation had to await the buildup of substantial galaxies with deeper potential wells? I briefly review here some of the physical processes that are conducive to the evolution of the massive black hole population. I will discuss black hole formation processes for `seed’ black holes that are likely to place at early cosmic epochs, and possible observational tests of these scenarios.

by Liu, Xin and Greene, Jenny E. and Shen, Yue and Strauss, Michael A.
Submitted to ApJL, Mar 10, 2010

We report the discovery of four kpc-scale binary AGNs. These objects were originally selected from the Sloan Digital Sky Survey based on double-peaked [O III] 4959,5007 emission lines in their fiber spectra. The double peaks could result from pairing active supermassive black holes (SMBHs) in a galaxy merger, or could be due to bulk motions of narrow-line region gas around a single SMBH. Deep near-infrared (NIR) images and optical slit spectra obtained from the Magellan 6.5 m and the APO 3.5 m telescopes strongly support the binary SMBH scenario for the four objects. In each system, the NIR images reveal tidal features and double stellar bulges with a projected separation of several kpc, while optical slit spectra show two Seyfert 2 nuclei spatially coincident with the stellar bulges, with line-of-sight velocity offsets of a few hundred km/s. These objects were drawn from a sample of only 43 objects, demonstrating the efficiency of this technique to find kpc-scale binary AGNs.

by Seth, Anil C. and Cappellari, Michele and Neumayer, Nadine and Caldwell, Nelson and Bastian, Nate and Olsen, Knut and Blum, Robert D. and Debattista, Victor P. and McDermid, Richard and Puzia, Thomas and Stephens, Andrew
To appear in ApJ. 21 pages, 19 figures, higher resolution PDF available at: http://www.cfa.harvard.edu/~aseth/ngc404paper.pdf

We examine the nuclear morphology, kinematics, and stellar populations in nearby S0 galaxy NGC 404 using a combination of adaptive optics assisted near-IR integral-field spectroscopy, optical spectroscopy, and HST imaging. These observations enable study of the NGC 404 nucleus at a level of detail possible only in the nearest galaxies. The surface brightness profile suggests the presence of three components, a bulge, a nuclear star cluster, and a central light excess within the cluster at radii <3 pc. These components have distinct kinematics with modest rotation seen in the nuclear star cluster and counter-rotation seen in the central excess. Molecular hydrogen emission traces a disk with rotation nearly orthogonal to that of the stars. The stellar populations of the three components are also distinct, with half of the mass of the nuclear star cluster having ages of ~1 Gyr (perhaps resulting from a galaxy merger), while the bulge is dominated by much older stars. Dynamical modeling of the stellar kinematics gives a total nuclear star cluster mass of 1.1×10^7 Msol. Dynamical detection of a possible intermediate mass black hole is hindered by uncertainties in the central stellar mass profile. Assuming a constant mass-to-light ratio, the stellar dynamical modeling suggests a black hole mass of <1×10^5 Msol, while the molecular hydrogen gas kinematics are best fit by a black hole with mass of 4.5×10^5 Msol. Unresolved and possibly variable dust emission in the near-infrared and AGN-like molecular hydrogen emission line ratios do suggest the presence of an accreting black hole in this nearby LINER galaxy.

by Volonteri, Marta
10 pages. To appear in the proceedings of the conference “Accretion and ejection in AGN: a global view” (Como, 22-26 June 2009)

Black hole spins affect the efficiency of the “classical” accretion processes, hence the radiative output from quasars. Spins also determine how much energy is extractable from the hole itself. Recently it became clear that massive black hole spins also affect the retention of black holes in galaxies, be cause of the impulsive “gravitational recoil”, up to thousands km/s, due to anisotropic emission of gravitational waves at merger. I discuss here the evolution of massive black hole spins along the cosmic history, due to the combination of mergers and accretion events. I describe recent simulations of accreting black holes in merger remnants, and discuss the implication for the spins of black holes in quasars.

by Farrell, S. A. and Servillat, M. and Oates, S. R. and Heywood, I. and Godet, O. and Webb, N. A. and Barret, D.
4 pages, 2 figures. Accepted 11th of Feb 2010. Contributed talk to appear in Proceedings of “X-ray Astronomy 2009: Present Status, Multi-Wavelength Approach and Future Perspectives”, Bologna, Italy, September 7-11, 2009, AIP, eds. A. Comastri, M. Cappi, and L. Angelini

The brightest Ultra-Luminous X-ray source HLX-1 in the galaxy ESO 243-49 currently provides strong evidence for the existence of intermediate mass black holes. Here we present the latest multi-wavelength results on this intriguing source in X-ray, UV and radio bands. We have refined the X-ray position to sub-arcsecond accuracy. We also report the detection of UV emission that could indicate ongoing star formation in the region around HLX-1. The lack of detectable radio emission at the X-ray position strengthens the argument against a background AGN.

by Holley-Bockelmann, Kelly and Micic, Miroslav and Sigurdsson, Steinn and Rubbo, Louis
18 pages, 10 figures, accepted in the Astrophysical Journal

We calculate the gravitational wave signal from the growth of 10 million solar mass supermassive black holes (SMBH) from the remnants of Population III stars. The assembly of these lower mass black holes is particularly important because observing SMBHs in this mass range is one of the primary science goals for the Laser Interferometer Space Antenna (LISA), a planned NASA/ESA mission to detect gravitational waves. We use high resolution cosmological N-body simulations to track the merger history of the host dark matter halos, and model the growth of the SMBHs with a semi-analytic approach that combines dynamical friction, gas accretion, and feedback. We find that the most common source in the LISA band from our volume consists of mergers between intermediate mass black holes and SMBHs at redshifts less than 2.

This type of high mass ratio merger has not been widely considered in the gravitational wave community; detection and characterization of this signal will likely require a different technique than is used for SMBH mergers or extreme mass ratio inspirals. We find that the event rate of this new LISA source depends on prescriptions for gas accretion onto the black hole as well as an accurate model of the dynamics on a galaxy scale; our best estimate yields about 40 sources with a signal-to-noise ratio greater than 30 occur within a volume like the Local Group during SMBH assembly — extrapolated over the volume of the universe yields roughly 500 observed events over 10 years, although the accuracy of this rate is affected by cosmic variance.

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 Callegari, S. and Kazantzidis, S. and Mayer, L. and Colpi, M. and Bellovary, J. M. and Quinn, T. and Wadsley, J.
5 pages, 4 figures, submitted to ApJ

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the evolution of massive black hole (MBH) pairs during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events occurring around a predicted peak of MBH pair formation at a redshift of $latex z \sim 3$, in the sensitivity window of the Laser Interferometer Space Antenna. Owing to strong tidal torques acting on its host and orbital circularization inside the disk of the primary galaxy, the companion MBH undergoes distinct episodes of enhanced accretion which cause an increase of the initial 1:10 mass ratio of the MBHs. We also find that the efficiency of MBH pair formation in the nuclei of the remnants correlates with the final mass ratio of the pair itself, so that MBH pairs with larger mass ratios are produced more effectively and promptly. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the final mass ratios of the resulting MBH pairs can be as large as 1:2, suggesting that minor galaxy mergers can give rise to MBH pairs with major mass ratios. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their host merging galaxies, but are a consequence of the complex interplay between accretion and merger dynamics.

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 Batcheldor, D.
ApJ Letters, accepted

The observed relation between supermassive black hole (SMBH) mass (M) and bulge stellar velocity dispersion (Sigma) is described by log(M) = alpha + beta*log(Sigma/200 km/s). As this relation has important implications for models of galaxy and SMBH formation and evolution, there continues to be great interest in adding to the M catalog. The “sphere of influence” (r) argument uses spatial resolution to exclude some M estimates and pre-select additional galaxies for further SMBH studies. This Letter quantifies the effects of applying the r argument to a population of galaxies and SMBHs that do not follow the M-Sigma relation. All galaxies with known values of Sigma, closer than 100 Mpc, are given a random M and selected when r is spatially resolved. These random SMBHs produce an M-Sigma relation of alpha=8.3, beta=4.0, consistent with observed values. Consequently, future proposed M estimates should not be justified solely on the basis of resolving r. This Letter shows the observed M-Sigma relation may simply be a result of available spatial resolution. However, it also implies the observed M-Sigma relation defines an upper limit. This potentially provides valuable new insight into the processes of galaxy and SMBH formation and evolution.

by Feoli, A. and Mancini, L. and Marulli, F. and Bergh, S. van den
16 pages, 5 figures, to be published in a special issue of General Relativity and Gravitation

The relation between the mass of supermassive black holes located in the center of the host galaxies and the kinetic energy of random motions of the corresponding bulges can be reinterpreted as an age-temperature diagram for galaxies. This relation fits the experimental data better than the M_bh-M_G, M_bh-L_G, and M_bh-sigma laws. The validity of this statement has been confirmed by using three samples extracted from different catalogues of galaxies. In the framework of the LambdaCDM cosmology our relation has been compared with the predictions of two galaxy formation models based on the Millennium Simulation.

by Neumayer, Nadine
8 pages, accepted for publication in PASA, contribution for “The Many Faces of Centaurus A” conference in Sydney, 2009

At less than 4 Mpc distance the radio galaxy NGC 5128 (Centaurus A) is the prime example to study the supermassive black hole and its influence on the environment in great detail. To model and understand the feeding and feedback mechanisms one needs an accurate determination of the mass of the supermassive black hole. The aim of this review is to give an overview of the recent studies that have been dedicated to measure the black hole mass in Centaurus A from both gas and stellar kinematics. It shows how the advancement in observing techniques and instrumentation drive the field of black hole mass measurements and concludes that adaptive optics assisted integral field spectroscopy is the key to identify the effects of the AGN on the surrounding ionised gas. Using data from SINFONI at the ESO Very Large Telescope, the best-fit black hole mass is M_BH=4.5 +1.7/-1.0 x 10^7 Msolar (from H_2 kinematics) and M_BH= (5.5 +/- 3.0) x 10^7 Msolar (from stellar kinematics; both with 3 sigma errors). This is one of the cleanest gas vs star comparison of a M_BH determination, and brings Centaurus A into agreement with the M_BH-sigma relation.

by Van Wassenhove, S. and Volonteri, M. and Walker, M. G. and Gair, J. R.
Submitted to MNRAS on November 30, 2009

As massive black holes (MBHs) grow from lower-mass seeds, it is natural to expect that a leftover population of progenitor MBHs should also exist in the present universe. Dwarf galaxies undergo a quiet merger history, and as a result, we expect that dwarfs observed in the local Universe retain some `memory’ of the original seed mass distribution. Consequently, the properties of MBHs in nearby dwarf galaxies may provide clean indicators of the efficiency of MBH formation. In order to examine the properties of MBHs in dwarf galaxies, we evolve different MBH populations within a Milky Way halo from high-redshift to today. We consider two plausible MBH formation mechanisms: `massive seeds’ formed via gas-dynamical instabilities and a Population III remnant seed model. `Massive seeds’ have larger masses than PopIII remnants, but form in rarer hosts. We dynamically evolve all halos merging with the central system, taking into consideration how the interaction modifies the satellites, stripping their outer mass layers. We compute different properties of the MBH population hosted in these satellites. We find that for the most part MBHs retain the original mass, thus providing a clear indication of what the properties of the seeds were. We derive the black hole occupation fraction (BHOF) of the satellite population at z=0. MBHs generated as `massive seeds’ have large masses that would favour their identification, but their typical BHOF is always below 40 per cent and decreases to less than per cent for observed dwarf galaxy sizes. In contrast, Population III remnants have a higher BHOF, but their masses have not grown much since formation, inhibiting their detection.

by Merritt, David
To appear in “The Galactic Center: A Window on the Nuclear Environment of Disk Galaxies”, ed. Mark Morris, Daniel Q. Wang and Feng Yuan

The distribution of late-type (old) stars in the inner parsec of the Milky Way is very different than expected for a relaxed population around a supermassive black hole. Instead of a density cusp, there is a 0.5 pc core. This article discusses what sorts of dynamical models might explain this “conundrum of old age.” A straightforward interpretation is that the nucleus is unrelaxed, and that the distribution of the old giants reflects the distribution of fainter stars and stellar remnants generally in the core. On the other hand, a density cusp could be present in the unobserved populations, and the deficit of bright giants could be a result of interactions with these objects. At the present time, no model is clearly preferred.

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 Devecchi, B. and Volonteri, M. and Colpi, M. and Haardt, F.
11 pages, 4 figures, submitted to MNRAS

Galactic nuclei host central massive objects either in the form of supermassive black holes or nuclear stellar clusters. Recent investigations have shown that both components co-exist in at least a few galaxies. In this paper we explore the possibility of a connection between nuclear star clusters and black holes that establishes at the moment of their formation. We here model the evolution of high redshift discs, hosted in dark matter halos with virial temperatures 10^4 K, whose gas has been polluted with metals just above the critical metallicity for fragmentation. A nuclear cluster forms as a result of a central starburst from gas inflowing from the unstable disc. The nuclear stellar cluster provides a suitable environment for the formation of a black hole seed, ensuing from runaway collisions among the most massive stars. Typical masses for the nuclear stellar clusters at the time of black hole formation (z~10) are inthe range 10^4-10^6 solar masses and have half mass radii < 0.5 pc. The black holes forming in these dense, high redshift clusters can have masses in the range ~300-2000 solar masses.

by Peterson, Bradley M.
10 pages, 2 figures. To be published in the Proceedings of IAU Symposium 267 “Co-Evolution of Central Black Holes and Galaxies”

We review briefly direct and indirect methods of measuring the masses of black holes in galactic nuclei, and then focus attention on supermassive black holes in active nuclei, with special attention to results from reverberation mapping and their limitations. We find that the intrinsic scatter in the relationship between the AGN luminosity and the broad-line region size is very small, ~0.11 dex, comparable to the uncertainties in the better reverberation measurements. We also find that the relationship between reverberation-based black hole masses and host-galaxy bulge luminosities also seems to have surprisingly little intrinsic scatter, ~0.17 dex. We note, however, that there are still potential systematics that could affect the overall mass calibration at the level of a factor of a few.

by Hayasaki, Kimitake and Ueda, Yoshihiro and Isobe, Naoki
9 pages, 11 figures, submitted to Publications of the Astronomical Society of Japan

If the activity of active galactic nuclei (AGNs) is predominantly induced by major galaxy mergers, then a significant fraction of AGNs should harbor binary massive black holes in their centers. We study the mass function of binary massive black holes in nearby AGNs based on the theory of evolution of binary massive black holes interacting with ambient gaseous disks proposed by Hayasaki (2009). The timescale of orbital decay is estimated as the order of $latex 10^8 yr$, being independent of the black hole mass but only dependent on the mass ratio and Eddington ratio. This makes it possible for any binary massive black holes to merge within a Hubble time. We find that $latex 1.3%–1.7%$ of the total number of nearby AGNs have close, binary massive black-holes with orbital period less than ten-years, detectable with on-going highly sensitive X-ray monitors such as Monitor of All-sky X-ray Image and/or Swift/Burst Alert Telescope. Close binaries with total black-hole masses of $latex 10^{6.5-7}M_sun$ are the most frequent in massive binary black-hole populations of nearby AGNs.

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 Volonteri, Marta and Gultekin, Kayhan and Dotti, Massimo
Submitted to MNRAS

We assess the influence of massive black hole (MBH) ejections from galaxy centres, due to the gravitational radiation recoil, along the cosmic merger history of the MBH population. We discuss the ‘danger’ of the recoil for MBHs as a function of different MBH spin/orbit configurations and of the host halo cosmic bias, and on how that reflects on the ‘occupation fraction’ of MBHs. We assess ejection probabilities for mergers occurring in a gas-poor environment, where the MBH binary coalescence is driven by stellar dynamical processes, and the spin/orbit configuration is expected to be isotropically distributed. We contrast this case with the ‘aligned’ case. The latter is the most realistic situation for ‘wet’, gas-rich mergers, which are the expectation for high-redshift galaxies. We find that if all halos at z>5-7 host a MBH, the probability of the Milky Way (or similar size galaxy) to host a MBH today is less than 50%, unless MBHs form continuously in galaxies. The ‘occupation fraction’ of MBHs, intimately related to halo bias and MBH formation efficiency, plays a crucial role in increasing the retention fraction. Small halos, with shallow potential wells and low escape velocities, have a high ejection probability, but the MBH merger rate is very low along their galaxy formation merger hierarchy: MBH formation processes are likely inefficient in such shallow potential wells. Recoils can decrease the overall frequency of MBHs in small galaxies to ~60%, while they have little effect on the frequency of MBHs in large galaxies (at most a 20% effect).

by Trenkel, Christian and Kemble, Steve and Bevis, Neil and Magueijo, Joao
Twin paper to arXiv:0912.0710

We suggest that LISA Pathfinder could be used to subject TEVES, and in particular the non-relativistic MOND phenomenology it incorporates, to a direct, controlled experimental test, in just a few years’ time. The basic concept is to fly LISA Pathfinder through the region around the Sun-Earth saddle point, following its nominal mission, in order to look for anomalous gravity gradients. We examine various strategies to reach the saddle point, and conclude that the preferred strategy, resulting in relatively short transfer times of order one year, probably involves a lunar fly-by. We present robust estimates of the MOND gravity gradients that LISA Pathfinder should be exposed to, and conclude that if the gradiometer on-board the spacecraft achieves its nominal performance, these gradients will not just be detected, but measured and characterised in some detail, should they exist. Conversely, given the large predicted signal based on standard assumptions, a null result would most likely spell the end of TEVES/MOND.

by Valtonen, M. and Mikkola, S. and Lehto, H. J. and Hyvönen, T. and Nilsson, K. and Merritt, D. and Gopakumar, A. and Rampadarath, H. and Hudec, R. and Basta, M. and Saunders, R.
12 pages, 4 figures, IAU261

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of $latex 1.84\cdot 10^{10} M_\odot$, the secondary mass $latex 1.46\cdot 10^{8} M_\odot$, major axis precession rate $latex 39^\circ.1$ per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is $latex 0.28\:\pm\:0.01$ (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the Post Newtonian equations of motion does not exactly follow Einstein’s theory: a parameter $latex q$ is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein’s value $latex q = 1$. This supports the $latex no-hair theorem$ of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein’s value with the accuracy of 2% (1 sigma).