Quasar Clustering in Cosmological Hydrodynamic Simulations: Evidence for mergers

We examine the clustering properties of a population of quasars drawn from fully hydrodynamic cosmological simulations that directly follow black hole growth. We find that the black hole correlation function is best described by two distinct components: contributions from BH pairs occupying the same dark matter halo ('1-halo term') which dominate at scales below 300 kpc/h, and contributions from BHs occupying separate halos ('2-halo term') which dominate at larger scales. From the 2-halo BH term we find a typical host halo mass for faint-end quasars (those probed in our simulation volumes) ranging from 10^11 to a few 10^12 solar masses from z=5 to z=1 respectively (consistent with the mean halo host mass). The BH correlation function shows a luminosity dependence as a function of redshift, though weak enough to be consistent with observational constraints. At small scales, the high resolution of our simulations allows us to probe the 1-halo clustering in detail, finding that the 1-halo term follows an approximate power law, lacking the characteristic decrease in slope at small scales found in 1-halo terms for galaxies and dark matter. We show that this difference is a direct result of a boost in the small-scale quasar bias caused by galaxies hosting multiple quasars (1-subhalo term) following a merger event, typically between a large central subgroup and a smaller, satellite subgroup hosting a relatively small black hole. We show that our predicted small-scale excess caused by such mergers is in good agreement with both the slope and amplitude indicated by recent small-scale measurements. Finally, we note the excess to be a strong function of halo mass, such that the observed excess is well matched by the multiple black holes of intermediate mass (10^7-10^8 solar masses) found in hosts of 4-8*10^11 solar masses, a range well probed by our simulations.Comment: 12 pages, 10 figures. Submitted to MNRA

Giant Clumps in Simulated High-z Galaxies: Properties, Evolution and Dependence on Feedback

We study the evolution of giant clumps in high-z disc galaxies using AMR cosmological simulations at redshifts z=6-1. Our sample consists of 34 galaxies, of halo masses 10^{11}-10^{12}M_s at z=2, run with and without radiation pressure (RP) feedback from young stars. While RP has little effect on the sizes and global stability of discs, it reduces the amount of star-forming gas by a factor of ~2, leading to a decrease in stellar mass by a similar factor by z~2. Both samples undergo violent disc instability (VDI) and form giant clumps of masses 10^7-10^9M_s at a similar rate, though RP significantly reduces the number of long-lived clumps. When RP is (not) included, clumps with circular velocity <40(20)km/s, baryonic surface density <200(100)M_s/pc^2 and baryonic mass <10^{8.2}(10^{7.3})M_s are short-lived, disrupted in a few free-fall times. The more massive and dense clumps survive and migrate toward the disc centre over a few disc orbital times. In the RP simulations, the distribution of clump masses and star-formation rates (SFRs) normalized to their host disc is very similar at all redshifts. They exhibit a truncated power-law with a slope slightly shallower than -2. Short-lived clumps preferentially have young stellar ages, low masses, high gas fractions and specific SFRs (sSFR), and they tend to populate the outer disc. The sSFR of massive, long-lived clumps declines with age as they migrate towards the disc centre, producing gradients in mass, stellar age, gas fraction, sSFR and metallicity that distinguish them from short-lived clumps. Ex situ mergers make up ~37% of the mass in clumps and ~29% of the SFR. They are more massive and with older stellar ages than the in situ clumps, especially near the disc edge. Roughly half the galaxies at redshifts z=4-1 are clumpy over a wide range of stellar mass, with clumps accounting for ~3-30% of the SFR but ~0.1-3% of the stellar mass.Comment: Accepted to MNRAS. Updated measurements of outflows from clumps and new study of clump virial parameters. Other small change

The Halo Occupation Distribution of Active Galactic Nuclei

Using a fully cosmological hydrodynamic simulation that self-consistently incorporates the growth and feedback of supermassive black holes and the physics of galaxy formation, we examine the effects of environmental factors (e.g., local gas density, black hole feedback) on the halo occupation distribution of low luminosity active galactic nuclei (AGN). We decompose the mean occupation function into central and satellite contribution and compute the conditional luminosity functions (CLF). The CLF of the central AGN follows a log-normal distribution with the mean increasing and scatter decreasing with increasing redshifts. We analyze the light curves of individual AGN and show that the peak luminosity of the AGN has a tighter correlation with halo mass compared to instantaneous luminosity. We also compute the CLF of satellite AGN at a given central AGN luminosity. We do not see any significant correlation between the number of satellites with the luminosity of the central AGN at a fixed halo mass. We also show that for a sample of AGN with luminosity above 10^42 ergs/s the mean occupation function can be modeled as a softened step function for central AGN and a power law for the satellite population. The radial distribution of AGN inside halos follows a power law at all redshifts with a mean index of -2.33 +/- 0.08. Incorporating the environmental dependence of supermassive black hole accretion and feedback, our formalism provides a theoretical tool for interpreting current and future measurements of AGN clustering.Comment: 14 pages, 11 figures, 2 Tables (Matches the MNRAS accepted version

The z=5 Quasar Luminosity Function from SDSS Stripe 82

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.Comment: 29 pages, 22 figures, ApJ accepted (updated to published version

Faint-end Quasar Luminosity Functions from Cosmological Hydrodynamic Simulations

We investigate the predictions for the faint-end quasar luminosity function (QLF) and its evolution using fully cosmological hydrodynamic simulations which self-consistently follow star formation, black hole growth and associated feedback processes. We find remarkably good agreement between predicted and observed faint end of the optical and X-ray QLFs (the bright end is not accessible in our simulated volumes) at z < 2. At higher redshifts our simulations tend to overestimate the QLF at the faintest luminosities. We show that although the low (high) luminosity ranges of the faint-end QLF are dominated by low (high) mass black holes, a wide range of black hole masses still contributes to any given luminosity range. This is consistent with the complex lightcurves of black holes resulting from the detailed hydrodynamics followed in the simulations. Consistent with the results on the QLFs, we find good agreement for the evolution of the comoving number density (in optical, soft and hard X-ray bands) of AGN for luminosities above 10^43 erg/s. However, the luminosity density evolution from the simulation appears to imply a peak at higher redshift than constrained from hard X-ray data (but not in optical). Our predicted excess at the faintest fluxes at z >= 2 does not lead to an overestimate to the total X-ray background and its contribution is at most a factor of two larger than the unresolved fraction of the 2-8 keV background. Even though this could be explained by some yet undetected, perhaps heavily obscured faint quasar population, we show that our predictions for the faint sources at high redshifts (which are dominated by the low mass black holes) in the simulations are likely affected by resolution effects.Comment: 12 pages, 9 figures; submitted and reviewed by MNRA

Morphological evolution of supermassive black hole merger hosts and multimessenger signatures

With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overall rate of mergers, and the correlation between merging black holes and their host galaxies. We find that these mergers occur in typical galaxies along the $M_{\rm{BH}}-M_*$ relation, and that between LISA and PTAs we expect to probe the full range of galaxy masses. As galaxy mergers can trigger increased star formation, we find that galaxies hosting low-mass black hole mergers tend to show a slight increase in star formation rates compared to a mass-matched sample. However, high-mass merger hosts have typical star formation rates, due to a combination of low gas fractions and powerful AGN feedback. Although minor black hole mergers do not correlate with disturbed morphologies, major mergers (especially at high-masses) tend to show morphological evidence of recent galaxy mergers which survives for ~500 Myr. This is on the same scale as the infall/hardening time of the merging black holes, suggesting that electromagnetic followups to gravitational wave signals may not be able to observe this correlation. We further find that incorporating a realistic timescale delay for the black hole mergers could shift the distribution of merger masses toward higher-masses, decreasing the rate of LISA detections while increasing the rate of PTA detections.Comment: 14 pages, 17 figures. Published in MNRA

The MassiveBlack-II simulation: the evolution of haloes and galaxies to z 0

(Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts.Comment: 26 pages, 25 figures. Submitted to MNRAS. High-resolution version and MBII galaxy catalogs can be found at http://mbii.phys.cmu.edu/data