Search

We present an analysis searching for dual AGN among 62 high-redshift ($2.5 < z < 3.5$) X-ray sources selected from publicly available deep Chandra fields. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis tool that calculates the Bayes factor for whether a given archival Chandra AGN is more likely a single or dual point source. We find no strong evidence for dual AGN in any individual source in our sample. We then increase our sensitivity to search for dual AGN across the sample by comparing our measured distribution of Bayes factors to that expected from a sample composed entirely of single point sources, and again find no evidence for dual AGN in the observed sample distribution. Although our analysis utilizes one of the largest Chandra catalogs of high-$z$ X-ray point sources available to study, the findings remain limited by the modest number of sources observed at the highest spatial resolution with Chandra and the typical count rates of the detected sources. Our non-detection allows us to place an upper-limit on the X-ray dual AGN fraction between $2.5

We investigate how the properties of massive black hole binaries influence the observed properties of core galaxies. We compare the observed trend in stellar mass deficit as a function of total stellar mass in the core galaxy with predicted trends in IllustrisTNG. We calculate mass deficits in simulated galaxies by applying sub-grid, post-processing physics based on the results of literature N-body experiments. We find the median value of the posterior distribution for the minimum binary mass ratio capable of creating a core is 0.7. For the gas mass fraction above which a core is erased we find a median value of 0.6. Thus low mass ratio binaries do not contribute to core formation and gas-rich mergers must lead to star formation within the nucleus, ultimately erasing the core. Such constraints have important implications for the overall massive black hole binary population, black hole-galaxy co-evolution, and the origin of the gravitational wave background., Comment: 10 pages, 12 figures, Accepted by MNRAS

We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3\sigma$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4\sigma$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal., Comment: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA., Comment: 93 pages with a lovely cover page thanks to Bernard Kelly and Elizabeth Ferrara

We present results from a multi-wavelength analysis searching for multiple AGN systems in nearby (z<0.077) triple galaxy mergers. Combining archival Chandra, SDSS, WISE, and VLA observations, we quantify the rate of nearby triple AGN, as well as investigate possible connections between SMBH accretion and merger environments. Analyzing the multi-wavelength observations of 7 triple galaxy mergers, we find that 1 triple merger has a single AGN (NGC 3341); we discover, for the first time, 4 likely dual AGN (SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356-1016); we confirm one triple AGN system, SDSS J0849+1114; and 1 triple merger in our sample remains ambiguous (SDSS J0858+1822). Analyzing the WISE data, we find a trend of increasing N_H (associated with the primary AGN) as a function of increasing W1-W2 color, reflecting that the motions of gas and dust are coupled in merging environments, where large amount of both can be funneled into the active central region during mergers. Additionally, we find that the one triple AGN system in our sample has the highest levels of N_H and W1-W2 color, while the dual AGN candidates all have lower levels; these results are consistent with theoretical merger simulations that suggest higher levels of nuclear gas are more likely to activate AGN in mergers., Comment: 19 pages, 3 figures, accepted to ApJ

We present results from our X-ray analysis of the first systematic search for triple AGN in nearby (z<0.077) triple galaxy mergers. We analyze archival Chandra observations of 7 triple galaxy mergers with BAYMAX (Bayesian Analysis of Multiple AGN in X-rays), fitting each observation with single, dual, and triple X-ray point source models. In doing so, we conclude that 1 triple merger has one X-ray point source (SDSS J0858+1822, although it's unlikely to be an AGN); 5 triple mergers are likely composed of two X-ray point sources (NGC 3341, SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356$-$1016); and one system is composed of three X-ray point sources (SDSS J0849+1114). By fitting the individual X-ray spectra of each point source, we analyze the 2-7 keV luminosities as well as the levels of obscuration associated with each potential AGN. We find that 4/5 dual X-ray point source systems have primary and secondary point sources with bright X-ray luminosities (L_2-7 kev >10^40 erg s^-1), possibly associated with 4 new undetected dual AGN. The dual and triple point source systems are found to have physical separations between 3-9 kpc and flux ratios between 2x10^-3 - 0.84. A multi-wavelength analysis to determine the origin of the X-ray point sources discovered in this work is presented in our companion paper (Foord et al. 2020c)., Comment: 27 pages, 8 figures, accepted to ApJ

We use Chandra X-ray observations to look for evidence of a recoiling black hole from the brightest cluster galaxy in Abell 2261 (A2261-BCG). A2261-BCG is a strong candidate for a recoiling black hole because of its large, flat stellar core, revealed by Hubble Space Telescope imaging observations. We took 100-ksec observations with Chandra and combined it with 35 ksec of archival observations to look for low-level accretion onto a black hole of expected mass $M\sim10^{10}\ M_{\scriptscriptstyle \odot}$ that could possibly be located in one of four off-center stellar knots near the galaxy's center or else in the optical center of the galaxy or in the location of radio emission. We found no X-ray emission arising from a point source in excess of the cluster gas and can place limits on the accretion of any black hole in the central region to a 2-7 keV flux below $4.3 \times 10^{-16}\ \mathrm{erg\ s^{-1}\ cm^{-2}}$, corresponding to a bolometric Eddington fraction of about $10^{-6}$. Thus there is either no $10^{10}\ M_{\scriptscriptstyle \odot}$ black hole in the core of A2261-BCG, or it is accreting at a low level. We also discuss the morphology of the X-ray emitting gas in the cluster and how its asymmetry is consistent with a large dynamic event., Comment: Accepted to AAS Journals. 12 pages of main text and figures, plus references and appendix tables and figures

We present the first results from BAYMAX (Bayesian AnalYsis of Multiple AGN in X-rays), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual AGNs is to use X-ray observations, only sources that are widely separated relative to the instrument PSF are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandra's angular resolution (<1''). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ~0.3'' from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find the probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0.3'' to be very low. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, and studying the dual AGN population where past spatial resolution limits have prevented systematic analyses., Comment: 16 pages, 5 figures, accepted to ApJ

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist., Comment: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA)., Comment: Submission to the 2020-2030 Astronomy and Astrophysics Decadal Survey (Astro2020), 10 pages, 3 figures

We have examined a sample of 13 sub-Eddington supermassive black holes hosted by galaxies spanning a variety of morphological classifications to further understand the empirical fundamental plane of black hole activity. This plane describes black holes from stellar-mass to supermassive and relates the mass of an accreting black hole and its radio and X-ray luminosities. A key factor in studying the fundamental plane is the turnover frequency, the frequency at which the radio continuum emission becomes optically thin. We measured this turnover frequency using new VLA observations combined, when necessary, with archival Chandra observations. Radio observations are in the range of 5--40 GHz across four frequency bands in B-configuration, giving high spatial resolution to focus on the core emission. We use Markov Chain Monte Carlo methods to fit the continuum emission in order to find the turnover frequency. After testing for correlations, the turnover frequency does not display a significant dependence on either mass or mass accretion rate, indicating that more complicated physics than simple scaling and optical depth effects are at play, as has been suggested by recent theoretical work., Comment: 10 pages, 4 figures, 2 tables, accepted for publication in ApJ

We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations., Comment: 9 pages, 2 figures, accepted to ApJ

We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done., Comment: 12 pages, 9 figures, submitted to ApJ

Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$\sigma$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs., Comment: ApJ accepted, 12 pages, 8 Figures, 4 Tables

Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations., Comment: Published in the Oct 22 2015 issue of Nature

Accreting black holes are observed to launch relativistic, collimated jets of matter and radiation. In some sources, discrete ejections have been detected with highly relativistic velocities. These particular sources typically have very high mass accretion rates, while sources lower knot velocities are predominantly associated with black holes with relatively low mass accretion rates. We quantify this behavior by examining knot velocity with respect to X-ray luminosity, a proxy for mass accretion rate onto the black hole. We find a positive correlation between the mass-scaled X-ray luminosity and jet knot velocity. In addition, we find evidence that the jet velocity is also a function of polar angle, supporting the "spine-sheath" model of jet production. Our results reveal a fundamental aspect of how accretion shapes mechanical feedback from black holes into their host environments., Comment: Accepted for publication in ApJ letters, 5 pages, 2 figures

We put active galactic nuclei (AGNs) with low-mass black holes on the fundamental plane of black hole accretion---the plane that relates X-ray emission, radio emission, and mass of an accreting black hole---to test whether or not the relation is universal for both stellar-mass and supermassive black holes. We use new Chandra X-ray and Very Large Array radio observations of a sample of black holes with masses less than $10^{6.3}\ M_{\scriptscriptstyle \odot}$, which have the best leverage for determining whether supermassive black holes and stellar-mass black holes belong on the same plane. Our results suggest that the two different classes of black holes both belong on the same relation. These results allow us to conclude that the fundamental plane is suitable for use in estimating supermassive black hole masses smaller than $\sim 10^7\ M_{\scriptscriptstyle \odot}$, in testing for intermediate-mass black holes, and in estimating masses at high accretion rates., Comment: ApJL accepted. 5 pages

Simultaneous observations of X-rays and radio luminosities have been well studied in accreting stellar-mass black holes. These observations are performed in order to understand how mass accretion rates and jetted outflows are linked in these individual systems. Such contemporaneous studies in supermassive black holes (SMBH) are harder to perform, as viscous times scale linearly with mass. However, as NGC 4395 is the lowest known mass Seyfert galaxy, we have used it to examine the simultaneous X-ray (Swift) and radio (VLA) correlation in a SMBH in a reasonably timed observing campaign. We find that the intrinsic X-ray variability is stronger than the radio variability, and that the fluxes are only weakly or tentatively coupled, similar to prior results obtained in NGC 4051. If the corona and the base of the jet are one and the same, this may suggest that the corona in radio-quiet AGN filters disk variations, only transferring the strongest and/or most sustained variations into the jet. Further, when both NGC 4395 and NGC 4051 are placed on the stellar-mass L_X-L_R plane, they appear to reside on the steeper L_X-L_R track. This suggests that SMBHs also follow two distinct tracks just as stellar-mass black holes do, and supports the idea that the same physical disk-jet mechanisms are at play across the mass scale., Comment: 9 pages, 10 figures, Accepted to ApJ Letters

We examine the role of spin in launching jets from compact objects across the mass scale. Our work includes a total of 37 Seyferts, 11 stellar-mass black holes, and 13 neutron stars. We find that when the Seyfert reflection lines are modeled with Gaussian line features (a crude proxy for inner disk radius and therefore spin), only a slight inverse correlation is found between the Doppler-corrected radio luminosity at 5 GHz (a proxy for jet power) and line width. When the Seyfert reflection features are fit with relativistically-blurred disk reflection models that measure spin, there is a tentative positive correlation between the Doppler-corrected radio luminosity and the spin measurement. Further, when we include stellar-mass black holes in the sample, to examine the effects across the mass scale, we find a slightly stronger correlation with radio luminosity per unit mass and spin, at a marginal significance (2.3 sigma confidence level). Finally, when we include neutron stars, in order to probe lower spin values, we find a positive correlation (3.3 sigma confidence level) between radio luminosity per unit mass and spin. Although tentative, these results suggest that spin may have a role in determining the jet luminosity. In addition, we find a slightly more significant correlation (4.4 sigma confidence level) between radio luminosity per Bolometric luminosity and spin, using our entire sample of black holes and neutrons stars. Again, although tentative, these relations point to the possibility that the mass accretion rate, i.e. Bolometric luminosity, is also important in determining the jet luminosity, in addition to spin. Our analysis suggests that mass accretion rate and disk or coronal magnetic field strength may be the "throttle" in these compact systems, to which the Eddington limit and spin may set the maximum jet luminosity that can be achieved., Comment: 14 pages, 13 Figures, ApJ Accepted

We present X-ray and radio observations of the new Galactic supernova remnant (SNR) G306.3-0.9, recently discovered by Swift. Chandra imaging reveals a complex morphology, dominated by a bright shock. The X-ray spectrum is broadly consistent with a young SNR in the Sedov phase, implying an age of 2500 yr for a distance of 8 kpc, plausibly identifying this as one of the 20 youngest Galactic SNRs. Australia Telescope Compact Array (ATCA) imaging reveals a prominent ridge of radio emission that correlates with the X-ray emission. We find a flux density of ~ 160 mJy at 1 GHz, which is the lowest radio flux recorded for a Galactic SNR to date. The remnant is also detected at 24microns, indicating the presence of irradiated warm dust. The data reveal no compelling evidence for the presence of a compact stellar remnant., Comment: 6 pages, 4 figures, 1 table, formatted using emulateapj.cls. Accepted for publication in ApJ

We calculate the observable signature of a black hole accretion disk with a gap or hole created by a secondary black hole embedded in the disk. We find that for an interesting range of parameters of black hole masses (~10^6 to 10^9 solar masses), orbital separation (~1 AU to ~0.1 pc), and gap width (10 to 190 disk scale heights), the missing thermal emission from a gap manifests itself in an observable decrement in the spectral energy distribution. We present observational diagnostics in terms of power-law forms that can be fit to line-free regions in AGN spectra or in fluxes from sequences of broad filters. Most interestingly, the change in slope in the broken power-law is almost entirely dependent on the width of gap in the accretion disk, which in turn is uniquely determined by mass ratio of the black holes, such that it scales roughly as q^(5/12). Thus one can use spectral observations of the continuum of bright active galactic nuclei to infer not only the presence of a closely separated black hole binary but also the mass ratio. When the black hole merger opens an entire hole (or cavity) in the inner disk, the broad band SED of the AGN or quasar may serve as a diagnostic. Such sources should be especially luminous in optical bands but intrinsically faint in X-rays (i.e., not merely obscured). We briefly note that viable candidates may have already been identified, though extant detailed modeling of those with high quality data have not yet revealed an inner cavity., Comment: Accepted by ApJ

We present a study of the mechanical power generated by both winds and jets across the black hole mass scale. We begin with the study of ionized X-ray winds and present a uniform analysis using Chandra grating spectra. The high quality grating spectra facilitate the characterization of the outflow velocity, ionization and column density of the absorbing gas. We find that the kinetic power of the winds scales with increasing bolometric luminosity as log(L_wind) \propto (1.58 \pm 0.07) log(L_Bol). This means that SMBH may be more efficient than stellar-mass black holes in launching winds. In addition, the simplicity of the scaling may suggest common driving mechanisms across the mass scale. For comparison, we next examine jet production, estimating jet power based on the energy required to inflate local bubbles. The jet relation is log(L_Jet)\propto (1.18\pm0.24) log(L_Bol). The energetics of the bubble associated with Cygnus X-1 are particularly difficult to determine, and the bubble could be a background SNR. If we exclude Cygnus X-1, then the jets follow a consistent relation to the winds within errors but with a higher normalization, log(L_Jet) \propto (1.34 \pm 0.50) log(L_Bol). The formal consistency in the wind and jet scaling relations suggests that a common launching mechanism may drive both flows; magnetic processes are viable possibilities. We also examine winds with especially high velocities, v > 0.01c. These ultra-fast outflows tend to resemble the jets more than the winds, indicating we may be observing a regime in which winds become jets. This study allows for the total power from black hole accretion, both mechanical and radiative, to be characterized in a simple manner and suggests a possible connection between winds and jets. Finally, we find at low Eddington fractions, the jet power is dominant, and at high Eddington fractions the wind power is dominant., Comment: 24 pages, 10 figures, Accepted to ApJ on 16 Nov 2012

We present Chandra observations of 12 galaxies that contain supermassive black holes with dynamical mass measurements. Each galaxy was observed for 30 ksec and resulted in a total of 68 point source detections in the target galaxies including supermassive black hole sources, ultraluminous X-ray sources, and extragalactic X-ray binaries. Based on our fits of the X-ray spectra, we report fluxes, luminosities, Eddington ratios, and slope of the power-law spectrum. Normalized to the Eddington luminosity, the 2--10 keV band X-ray luminosities of the SMBH sources range from $10^{-8}$ to $10^{-6}$, and the power-law slopes are centered at $\sim2$ with a slight trend towards steeper (softer) slopes at smaller Eddington fractions, implying a change in the physical processes responsible for their emission at low accretion rates. We find 20 ULX candidates, of which six are likely ($>90%$ chance) to be true ULXs. The most promising ULX candidate has an isotropic luminosity in the 0.3--10 keV band of $1.0_{-0.3}^{+0.6} \times 10^{40}$ erg/s., Comment: Accepted by ApJ. 16 pages, 8 figures, 5 tables

We present HST STIS observations of the galaxy NGC 4382 (M85) and axisymmetric models of the galaxy to determine mass-to-light ration (M/L, V-band) and central black hole mass (M_BH). We find M/L = 3.74 +/- 0.1 (solar units) and M_BH = 1.3 (+5.2, -1.2) \times 10^7 M_sun at an assumed distance of 17.9 Mpc, consistent with no black hole. The upper limit, M_BH < 9.6 \times 10^7 M_sun (2{\sigma}) or M_BH < 1.4 \times 10^8 M_sun (3{\sigma}) is consistent with the current M-{\sigma} relation, which predicts M_BH = 8.8 \times 10^7 M_sun at {\sigma}_e = 182 km/s, but low for the current M-L relation, which predicts M_BH = 7.8 \times 10^8 M_sun at L_V = 8.9 \times 10^10 L_sun,V. HST images show the nucleus to be double, suggesting the presence of a nuclear eccentric stellar disk, in analogy to the Tremaine disk in M31. This conclusion is supported by the HST velocity dispersion profile. Despite the presence of this non-axisymmetric feature and evidence of a recent merger, we conclude that the reliability of our black hole mass determination is not hindered. The inferred low black hole mass may explain the lack of nuclear activity., Comment: 10 pages, 8 figures, accepted by ApJ

We examine the possibility that the observed relation between black-hole mass and host-galaxy stellar velocity dispersion (the M-sigma relation) is biased by an observational selection effect, the difficulty of detecting a black hole whose sphere of influence is smaller than the telescope resolution. In particular, we critically investigate recent claims that the M-sigma relation only represents the upper limit to a broad distribution of black-hole masses in galaxies of a given velocity dispersion. We find that this hypothesis can be rejected at a high confidence level, at least for the early-type galaxies with relatively high velocity dispersions (median 268 km/s) that comprise most of our sample. We also describe a general procedure for incorporating observational selection effects in estimates of the properties of the M-sigma relation. Applying this procedure we find results that are consistent with earlier estimates that did not account for selection effects, although with larger error bars. In particular, (i) the width of the M-sigma relation is not significantly increased; (ii) the slope and normalization of the M-sigma relation are not significantly changed; (iii) most or all luminous early-type galaxies contain central black holes at zero redshift. Our results may not apply to late-type or small galaxies, which are not well-represented in our sample., Comment: 8 pages, ApJ accepted

In this paper, we examine if the properties of central black holes in galactic nuclei correlate with their host dark matter halo. We analyze the entire sample of galaxies where black hole mass, velocity dispersion, sigma, and asymptotic circular velocity, Vc, have all been measured. We fit M-sigma and M-Vc to a power law, and find that in both relationships the scatter and slope are similar. This model-independent analysis suggests that although the black hole masses are not uniquely determined by dark matter halo mass, when considered for the current sample as a whole, the M-Vc correlation may be as strong (or as weak) as M-sigma. Although the data are sparse, there appears to be more scatter in the correlation for both sigma and Vc at the low--mass end. This is not unexpected given our current understanding of galaxy and black hole assembly. In fact, there are several compelling reasons that account for this: (i) SMBH formation is likely less efficient in low-mass galaxies with large angular momentum content; (ii) SMBH growth is less efficient in low-mass disk galaxies that have not experienced major mergers; (iii) dynamical effects, such as gravitational recoil, increase scatter preferentially at the low-mass end. Therefore, the recent observational claim of the absence of central SMBHs in bulgeless, low mass galaxies, or deviations from the correlations defined by high-mass black holes in large galaxies today is, in fact, predicated by current models of black hole growth. We show how this arises as a direct consequence of the coupling between dark matter halos and central black holes at the earliest epochs., Comment: ApJ, in press

We present the stellar kinematics in the central 2" of the luminous elliptical galaxy M87 (NGC 4486), using laser adaptive optics to feed the Gemini telescope integral-field spectrograph, NIFS. The velocity dispersion rises to 480 km/s at 0.2". We combine these data with extensive stellar kinematics out to large radii to derive a black-hole mass equal to (6.6+-0.4)x10^9 Msun, using orbit-based axisymmetric models and including only the NIFS data in the central region. Including previously-reported ground-based data in the central region drops the uncertainty to 0.25x10^9 Msun with no change in the best-fit mass; however, we rely on the values derived from the NIFS-only data in the central region in order to limit systematic differences. The best-fit model shows a significant increase in the tangential velocity anisotropy of stars orbiting in the central region with decreasing radius; similar to that seen in the centers of other core galaxies. The black-hole mass is insensitive to the inclusion of a dark halo in the models --- the high angular-resolution provided by the adaptive optics breaks the degeneracy between black-hole mass and stellar mass-to-light ratio. The present black-hole mass is in excellent agreement with the Gebhardt & Thomas value, implying that the dark halo must be included when the kinematic influence of the black hole is poorly resolved. This degeneracy implies that the black-hole masses of luminous core galaxies, where this effect is important, may need to be re-evaluated. The present value exceeds the prediction of the black hole-dispersion and black hole-luminosity relations, both of which predict about 1x10^9 Msun for M87, by close to twice the intrinsic scatter in the relations. The high-end of the black hole correlations may be poorly determined at present., Comment: 12 pages, accepted for publication in the Astrophysical Journal

We report on the results of a simultaneous monitoring campaign employing eight Chandra X-ray (0.5-10 keV) and six VLA/EVLA (8.4 GHz) radio observations of NGC 4051 over seven months. Evidence for compact jets is observed in the 8.4 GHz radio band; This builds on mounting evidence that jet production may be prevalent even in radio-quiet Seyferts. Assuming comparatively negligible local diffuse emission in the nucleus, the results also demonstrate an inverse correlation of L_radio proportional to L_X-ray ^(-0.72+/-0.04) . Current research linking the mass of supermassive black holes and stellar-mass black holes in the "low/hard" state to X-ray luminosities and radio luminosities suggest a "fundamental plane of accretion onto black holes" that has a positive correlation of L_radio proportional to L_X-ray^(0.67+/-0.12) . Our simultaneous results differ from this relation by more than 11 sigma, indicating that a separate mode of accretion and ejection may operate in this system. A review of the literature shows that the inverse correlation seen in NGC 4051 is seen in three other black hole systems, all of which accrete at near 10% of their Eddington luminosity, perhaps suggesting a distinct mode of disk-jet coupling at high Eddington fractions. We discuss our results in the context of disk and jets in black holes and accretion across the black hole mass scale., Comment: 12 pages, 9 figures

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)., Comment: Accepted for publication in MNRAS

We report on recently derived improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on ~50 M_BH measurements and ~20 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. The scatter was found to be significantly larger than estimated in most previous studies. The large scatter requires revision of the local black hole mass function, and it implies that there may be substantial selection bias in studies of the evolution of the M-sigma and M-L relations. When only considering ellipticals, the scatter appears to decrease. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We also report on the effects on the fits of culling the sample according to the resolution of the black hole's sphere of influence., Comment: 6 pages. To appear in the proceedings of the IAU Symposium 267, Co-evolution of Central Black Holes and Galaxies

Black hole accretion and jet production are areas of intensive study in astrophysics. Recent work has found a relation between radio luminosity, X-ray luminosity, and black hole mass. With the assumption that radio and X-ray luminosity are suitable proxies for jet power and accretion power, respectively, a broad fundamental connection between accretion and jet production is implied. In an effort to refine these links and enhance their power, we have explored the above relations exclusively among black holes with direct, dynamical mass-measurements. This approach not only eliminates systematic errors incurred through the use of secondary mass measurements, but also effectively restricts the range of distances considered to a volume-limited sample. Further, we have exclusively used archival data from the Chandra X-ray Observatory to best isolate nuclear sources. We find log(L_R) = (4.80 +/- 0.24) + (0.78 +/- 0.27) log(M_BH) + (0.67 +/- 0.12) log(L_X), in broad agreement with prior efforts. Owing to the nature of our sample, the plane can be turned into an effective mass predictor. When the full sample is considered, masses are predicted less accurately than with the well-known M-sigma relation. If obscured AGN are excluded, the plane is potentially a better predictor than other scaling measures., Comment: Accepted by ApJ. 13 pages, 6 figures, updated to include changes to match accepted version

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals., Comment: 27 pages, 18 figures, 7 tables, ApJ accepted

We report five new measurements of central black hole masses based on STIS and WFPC2 observations with the Hubble Space Telescope and on axisymmetric, three-integral, Schwarzschild orbit-library kinematic models. We selected a sample of galaxies within a narrow range in velocity dispersion that cover a range of galaxy parameters (including Hubble type and core/power-law surface density profile) where we expected to be able to resolve the galaxy's sphere of influence based on the predicted value of the black hole mass from the M-sigma relation. We find masses in units of 10^8 solar masses for the following galaxies: NGC 3585, M_BH = 3.4 (+1.5, -0.6); NGC 3607, M_BH = 1.2 (+0.4, -0.4); NGC 4026, M_BH = 2.1 (+0.7, -0.4); and NGC 5576, M_BH = 1.8 (+0.3, -0.4), all significantly excluding M_BH = 0. For NGC 3945, M_BH = 0.09 (+0.17, -0.21), which is significantly below predictions from M-sigma and M-L relations and consistent with M_BH = 0, though the presence of a double bar in this galaxy may present problems for our axisymmetric code., Comment: 15 pages, 21 figures, 13 tables. Accepted by ApJ

We investigate the distribution of massive black holes (MBHs) in the Virgo cluster. Observations suggest that AGN activity is widespread in massive galaxies (M>1e10 solar masses), while at lower galaxy masses star clusters are more abundant, which might imply a limited presence of central black holes in these galaxy-mass regimes. We explore if this possible threshold in MBH hosting, is linked to nature, nurture, or a mixture of both. The nature scenario arises naturally in hierarchical cosmologies, as MBH formation mechanisms typically are efficient in biased systems, which would later evolve into massive galaxies. Nurture, in the guise of MBH ejections following MBH mergers, provides an additional mechanism that is more effective for low mass, satellite galaxies. The combination of inefficient formation, and lower retention of MBHs, leads to the natural explanation of the distribution of compact massive ob jects in Virgo galaxies. If MBHs arrive to the correlation with the host mass and velocity dispersion during merger-triggered accretion episodes, sustained tidal stripping of the host galaxies creates a population of MBHs which lie above the expected scaling between the holes and their host mass, suggesting a possible environmental dependence., Comment: MNRAS letters

During the inspiral and merger of a binary black hole, gravitational radiation is emitted anisotropically due to asymmetries in the merger configuration. This anisotropic radiation leads to a gravitational wave kick, or recoil velocity, as large as ~ 4000 km/sec. We investigate the effect gravitational recoil has on the retention of intermediate mass black holes (IMBH) within Galactic globular clusters. Assuming that our current understanding of IMBH-formation is correct and yields an IMBH-seed in every globular cluster, we find a significant problem retaining low mass IMBHs (1000 $\Msun$) in the typical merger-rich globular cluster environment. Given a uniform black hole spin distribution and orientation and a Kroupa IMF, we find that at most 3% of the globular clusters can retain an IMBH larger than 1000 $\Msun$ today. For a population of black holes that better approximates mass loss from winds and supernovae, we find that 16% of globulars can retain an IMBH larger than 1000 $\Msun$. Our calculations show that if there are black holes of mass $M > 60 \Msun$ in a cluster, repeated IMBH-BH encounters will eventually eject a 1000 $\Msun$ IMBH with greater than 30% probability. As a consequence, a large population of rogue black holes may exist in our Milky Way halo. We discuss the dynamical implications of this subpopulation, and its possible connection to ultraluminous X-ray sources (ULXs)., Comment: 8 pages, 6 figures. Pick up high resolution version of figure 1 at http://gravity.psu.edu/~kellyhb/kickfigures. Submitted to Astrophysical Journal

We present numerical three-body experiments that include the effects of gravitational radiation reaction by using equations of motion that include the 2.5-order post-Newtonian force terms, which are the leading order terms of energy loss from gravitational waves. We simulate binary-single interactions and show that close approach cross sections for three 1 solar mass objects are unchanged from the purely Newtonian dynamics except for close approaches smaller than 1.0e-5 times the initial semimajor axis of the binary. We also present cross sections for mergers resulting from gravitational radiation during three-body encounters for a range of binary semimajor axes and mass ratios including those of interest for intermediate-mass black holes (IMBHs). Building on previous work, we simulate sequences of high-mass-ratio three-body encounters that include the effects of gravitational radiation. The simulations show that the binaries merge with extremely high eccentricity such that when the gravitational waves are detectable by LISA, most of the binaries will have eccentricities e > 0.9 though all will have circularized by the time they are detectable by LIGO. We also investigate the implications for the formation and growth of IMBHs and find that the inclusion of gravitational waves during the encounter results in roughly half as many black holes ejected from the host cluster for each black hole accreted onto the growing IMBH., Comment: 34 pages, 14 figures, minor corrections to match version accepted by ApJ

We present results of numerical simulations of sequences of binary-single scattering events of black holes in dense stellar environments. The simulations cover a wide range of mass ratios from equal mass objects to 1000:10:10 solar masses and compare purely Newtonian simulations to simulations in which Newtonian encounters are interspersed with gravitational wave emission from the binary. In both cases, the sequence is terminated when the binary's merger time due to gravitational radiation is less than the arrival time of the next interloper. We find that black hole binaries typically merge with a very high eccentricity (0.93 < e < 0.95 pure Newtonian; 0.85 < e < 0.90 with gravitational wave emission) and that adding gravitational wave emission decreases the time to harden a binary until merger by ~ 30% to 40%. We discuss the implications of this work for the formation of intermediate-mass black holes and gravitational wave detection., Comment: 28 pages including 9 figures, submitted to ApJ

Evidence has been mounting for the existence of black holes with masses from 10^2 to 10^4 M_Solar associated with stellar clusters. Such intermediate-mass black holes (IMBHs) will encounter other black holes in the dense cores of these clusters. The binaries produced in these interactions will be perturbed by other objects as well thus changing the orbital characteristics of the binaries. These binaries and their subsequent mergers due to gravitational radiation are important sources of gravitational waves. We present the results of numerical simulations of high mass ratio encounters, which help clarify the interactions of intermediate-mass black holes in globular clusters and help determine what types of detectable gravitational wave signatures are likely., Comment: 4 pages, 3 figures to appear in the proceedings of The Astrophysics of Gravitational Wave Sources, College Park, MD, 24-26 April 2003