-
High-redshift supermassive black hole mergers in simulations with dynamical friction modelling
Authors:
Colin DeGraf,
Nianyi Chen,
Yueying Ni,
Tiziana Di Matteo,
Simeon Bird,
Michael Tremmel,
Rupert Croft
Abstract:
In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological si…
▽ More
In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass BHs, and directly incorporates dynamical friction so as to accurately model the dynamical processes which bring black holes to the galaxy center where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to 10$^{4.5}$ M$_\odot$, and remain broadly consistent with the TNG simulations at scales > 10$^6$ M$_\odot$ (the seed mass used in TNG). By resolving lower-mass black holes, the overall merger rate is ~5x higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentering scheme, reducing the high-z merger rate mass-matched mergers by a factor of ~2x. We also calculate the expected LISA Signal-to-Noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA's peak sensitivity (10$^6$ M$_\odot$) to resolve the majority of LISA's GW detections.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
Properties and Evolution of Dual and Offset AGN in the ASTRID Simulation at $z \sim 2$
Authors:
Nianyi Chen,
Tiziana Di Matteo,
Yueying Ni,
Michael Tremmel,
Colin DeGraf,
Yue Shen,
A. Miguel Holgado,
Simeon Bird,
Rupert Croft,
Yu Feng
Abstract:
We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Com…
▽ More
We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Compared with the full AGN population, duals are characterized by a low $M_\text{BH}/M_*$ ratio, a high specific star-formation rate (sSFR) of $\sim 1\,\text{Gyr}^{-1}$, and a high Eddington ratio ($\sim 0.05$, double that of single AGN). The dual AGN are formed in major galaxy mergers (typically involving $M_\text{halo}<10^{13}\,M_\odot$), with BHs that have similar masses. At small separations (when their host galaxies are in the late phase of the merger) duals become $2\sim8$ times brighter (albeit more obscured) than at larger separations. $80\%$ of these bright, close duals merge in the simulation within $\sim500\,\text{Myrs}$. Notably, the initially less-massive BH in duals frequently becomes the brighter AGN during the galaxy merger. In offset AGN, the active BH is typically $\gtrsim 10$ times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the center of massive halos ($ M_\text{ halo}\sim 10^{13-14}\,M_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $Δr\sim5\,\text{kpc}$ for a few hundred Myrs.
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
An estimate of the stochastic gravitational wave background from the MassiveBlackII simulation
Authors:
Bailey Sykes,
Hannah Middleton,
Andrew Melatos,
Tiziana Di Matteo,
Colin DeGraf,
Aklant Bhowmick
Abstract:
A population of super-massive black hole binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of $10^{-9}$--$10^{-7}$ Hz. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work we use super-massive black hole binary mergers from the Ma…
▽ More
A population of super-massive black hole binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of $10^{-9}$--$10^{-7}$ Hz. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work we use super-massive black hole binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disk. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1 yr^{-1}$ is found to be $A_{yr^{-1}} = 6.9 \times 10^{-16}$ and $A_{yr^{-1}} = 6.4 \times 10^{-16}$ in the gravitational-wave driven and stellar scattering/gas disk cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely gravitational wave driven model as frequency decreases inside the PTA frequency band.
△ Less
Submitted 10 February, 2022;
originally announced February 2022.
-
Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation
Authors:
Nianyi Chen,
Yueying Ni,
A. Miguel Holgado,
Tiziana Di Matteo,
Michael Tremmel,
Colin DeGraf,
Simeon Bird,
Rupert Croft,
Yu Feng
Abstract:
We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to…
▽ More
We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs.
△ Less
Submitted 15 December, 2021;
originally announced December 2021.
-
The ASTRID simulation: the evolution of Supermassive Black Holes
Authors:
Yueying Ni,
Tiziana Di Matteo,
Simeon Bird,
Rupert Croft,
Yu Feng,
Nianyi Chen,
Michael Tremmel,
Colin DeGraf,
Yin Li
Abstract:
We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses…
▽ More
We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers.
△ Less
Submitted 20 March, 2022; v1 submitted 27 October, 2021;
originally announced October 2021.
-
On the detectability of massive black hole merger eventsby LISA
Authors:
Samuel Banks,
Katharine Lee,
Nazanin Azimi,
Kendall Scarborough,
Nikolai Stefanov,
Indra Periwal,
Colin DeGraf,
Tiziana Di Matteo
Abstract:
The launch of space based gravitational wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays (PTAs) will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for th…
▽ More
The launch of space based gravitational wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays (PTAs) will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for the upcoming low frequency ($10^{-4}-10^{-1}$ Hz) window probed by LISA. It is important to investigate the expected MBH merger rates and associated signals, to determine how potential LISA events are affected by physics included in current models. To study this, we post-process the large population of MBHBs in the Illustris simulation to account for dynamical friction time delays associated with BH infall/inspiral. We show that merger delays associated with binary evolution have the potential to decrease the expected merger rates, with $M_{\rm{BH}} > 10^6 M_\odot$ MBHBs (the lowest mass in Illustris) decreasing from $\sim 3$ yr$^{-1}$ to $\sim 0.1 $yr$^{-1}$, and shifting the merger peak from z $\sim 2$ to $\sim 1.25$. During this time, we estimate that accretion grows the total merging mass by as much as 7x from the original mass. Importantly, however, dynamical friction associated delays (which shift the mergers toward lower-redshift and higher-masses) lead to a stronger signal/strain for the emitted GWs in the LISA band, increasing mean frequency from $10^{-3.1}$ to $10^{-3.4}-10^{-4.0}$ Hz, and mean strain from $10^{-17.2}$ to $10^{-16.3}-10^{-15.3}$. Finally, we show that after including a merger delay and associated $M_{\rm{BH}}$ growth, mergers still tend to lie on the typical $M_{\rm{BH}}-M_*$ relation, but with an increased likelihood of an undermassive black hole.
△ Less
Submitted 19 July, 2021;
originally announced July 2021.
-
Dynamical Friction Modeling of Massive Black Holes in Cosmological Simulations and Effects on Merger Rate Predictions
Authors:
Nianyi Chen,
Yueying Ni,
Michael Tremmel,
Tiziana Di Matteo,
Simeon Bird,
Colin DeGraf,
Yu Feng
Abstract:
In this work we establish and test methods for implementing dynamical friction for massive black hole pairs that form in large volume cosmological hydrodynamical simulations which include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of dynamical friction…
▽ More
In this work we establish and test methods for implementing dynamical friction for massive black hole pairs that form in large volume cosmological hydrodynamical simulations which include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of dynamical friction (DF) from collisionless particles, black holes can effectively sink close to the galaxy center, provided that the black hole's dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes' orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole's evolution. DF from gas becomes less than $1\%$ of DF from collisionless particles for BH masses $> 10^{7}$ M$_{\odot}$. Using our best DF model, we calculate the merger rate down to $z=1.1$ using an $L_{\rm box}=35$ Mpc$/h$ simulation box. We predict $\sim 2$ mergers per year for $z>1.1$ peaking at $z\sim 2$. These merger rates are within the range obtained in previous work using similar-resolution hydro-dynamical simulations. We show that the rate is enhanced by factor of $\sim 2$ when DF is taken into account in the simulations compared to the no-DF run. This is due to $>40\%$ more black holes reaching the center of their host halo when DF is added.
△ Less
Submitted 2 April, 2021; v1 submitted 31 March, 2021;
originally announced April 2021.
-
Resolving Disks & Mergers in $z\sim2$ Heavily Reddened Quasars and their Companion Galaxies with ALMA
Authors:
Manda Banerji,
Gareth C. Jones,
Stefano Carniani,
Colin DeGraf,
Jeff Wagg
Abstract:
We present sub-arcsecond resolution ALMA imaging of the CO(3-2) emission in two $z\sim2.5$ heavily reddened quasars (HRQs) - ULASJ1234+0907 and ULASJ2315+0143 - and their companion galaxies. Dynamical modeling of the resolved velocity fields enables us to constrain the molecular gas morphologies and host galaxy masses. Combining the new data with extensive multi-wavelength observations, we are abl…
▽ More
We present sub-arcsecond resolution ALMA imaging of the CO(3-2) emission in two $z\sim2.5$ heavily reddened quasars (HRQs) - ULASJ1234+0907 and ULASJ2315+0143 - and their companion galaxies. Dynamical modeling of the resolved velocity fields enables us to constrain the molecular gas morphologies and host galaxy masses. Combining the new data with extensive multi-wavelength observations, we are able to study the relative kinematics of different molecular emission lines, the molecular gas fractions and the locations of the quasars on the M$_{\rm{BH}}$-M$_{\rm{gal}}$ relation. Despite having similar black-hole properties, the two HRQs display markedly different host galaxy properties and local environments. J1234 has a very massive host, M$_{\rm{dyn}} \sim 5 \times 10^{11}$M$_\odot$ and two companion galaxies that are similarly massive located within 200 kpc of the quasar. The molecular gas fraction is low ($\sim$6%). The significant ongoing star formation in the host galaxy is entirely obscured at rest-frame UV and optical wavelengths. J2315 is resolved into a close-separation major-merger ($Δ$r=15 kpc; $Δ$v=170 km/s) with a $\sim$1:2 mass ratio. The total dynamical mass is estimated to be $\lesssim$10$^{11}$M$_\odot$ and the molecular gas fraction is high ($>$45%). A new HSC image of the galaxy shows unobscured UV-luminous star-forming regions co-incident with the extended reservoir of cold molecular gas in the merger. We use the outputs from the Illustris simulations to track the growth of such massive black holes from $z\sim6$ to the present day. While J1234 is consistent with the simulated $z\sim2$ relation, J2315 has a black hole that is over-massive relative to its host galaxy.
△ Less
Submitted 18 March, 2021;
originally announced March 2021.
-
Morphological evolution of supermassive black hole merger hosts and multimessenger signatures
Authors:
Colin DeGraf,
Debora Sijacki,
Tiziana Di Matteo,
Kelly Holley-Bockelmann,
Greg Snyder,
Volker Springel
Abstract:
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 overal…
▽ More
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.
△ Less
Submitted 28 July, 2021; v1 submitted 1 December, 2020;
originally announced December 2020.
-
Cosmological simulations of massive black hole seeds: predictions for next generation electromagnetic and gravitational wave observations
Authors:
Colin DeGraf,
Debora Sijacki
Abstract:
We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated mo…
▽ More
We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated model where black hole seeds only form in progenitor galaxies that conform to the conditions for direct collapse black hole formation. While all seed models considered here are in a broad agreement with present observational constraints on black hole populations from optical, UV and X-ray studies, we find that they lead to widely different black hole number densities and halo occupation fractions which are currently observationally unconstrained. In terms of future electromagnetic spectrum observations, the faint-end quasar luminosity function and the low-mass-end black hole-host galaxy scaling relations are very sensitive to the specific massive seed prescription. Specifically, the direct collapse model exhibits a seeding efficiency which decreases rapidly with cosmic time and produces much fewer black holes in low mass galaxies, in contrast to the original Illustris simulation. We further find that the total black hole merger rate varies by more than one order of magnitude for different seed models, with the redshift evolution of the chirp mass changing as well. Supermassive black hole merger detections with LISA and International Pulsar Timing Array may hence provide the most direct means of constraining massive black hole seed formation in the early Universe.
△ Less
Submitted 29 July, 2021; v1 submitted 26 June, 2019;
originally announced June 2019.
-
The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
Authors:
Ji-hoon Kim,
Oscar Agertz,
Romain Teyssier,
Michael J. Butler,
Daniel Ceverino,
Jun-Hwan Choi,
Robert Feldmann,
Ben W. Keller,
Alessandro Lupi,
Thomas Quinn,
Yves Revaz,
Spencer Wallace,
Nickolay Y. Gnedin,
Samuel N. Leitner,
Sijing Shen,
Britton D. Smith,
Robert Thompson,
Matthew J. Turk,
Tom Abel,
Kenza S. Arraki,
Samantha M. Benincasa,
Sukanya Chakrabarti,
Colin DeGraf,
Avishai Dekel,
Nathan J. Goldbaum
, et al. (18 additional authors not shown)
Abstract:
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st…
▽ More
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
△ Less
Submitted 8 March, 2018; v1 submitted 10 October, 2016;
originally announced October 2016.
-
Black hole clustering and duty cycles in the Illustris simulation
Authors:
Colin DeGraf,
Debora Sijacki
Abstract:
We use the high-resolution cosmological simulation Illustris to investigate the clustering of supermassive black holes across cosmic time, the link between black hole clustering and host halo masses, and the implications for black hole duty cycles. Our predicted black hole correlation length and bias match the observational data very well across the full redshift range probed. Black hole clusterin…
▽ More
We use the high-resolution cosmological simulation Illustris to investigate the clustering of supermassive black holes across cosmic time, the link between black hole clustering and host halo masses, and the implications for black hole duty cycles. Our predicted black hole correlation length and bias match the observational data very well across the full redshift range probed. Black hole clustering is strongly luminosity-dependent on small, 1-halo scales, with some moderate dependence on larger scales of a few Mpc at intermediate redshifts. We find black hole clustering to evolve only weakly with redshift, initially following the behaviour of their hosts. However below z ~ 2 black hole clustering increases faster than that of their hosts, which leads to a significant overestimate of the clustering-predicted host halo mass. The full distribution of host halo masses is very wide, including a low-mass tail extending up to an order of magnitude below the naive prediction for minimum host mass. Our black hole duty cycles follow a power-law dependence on black hole mass and decrease with redshift, and we provide accurate analytic fits to these. The increase in clustering amplitude at late times, however, means that duty cycle estimates based on black hole clustering can overestimate duty cycles substantially, by more than two orders of magnitude. We find the best agreement when the minimum host mass is assumed to be $10^{11.2} M_\odot$ , which provides an accurate measure across all redshifts and luminosity ranges probed by our simulation.
△ Less
Submitted 20 January, 2017; v1 submitted 21 September, 2016;
originally announced September 2016.
-
Giant Clumps in Simulated High-z Galaxies: Properties, Evolution and Dependence on Feedback
Authors:
Nir Mandelker,
Avishai Dekel,
Daniel Ceverino,
Colin DeGraf,
Yicheng Guo,
Joel Primack
Abstract:
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,…
▽ More
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.
△ Less
Submitted 15 September, 2016; v1 submitted 29 December, 2015;
originally announced December 2015.
-
The Confinement of Star-Forming Galaxies into a Main Sequence through Episodes of Gas Compaction, Depletion, and Replenishment
Authors:
Sandro Tacchella,
Avishai Dekel,
C. Marcella Carollo,
Daniel Ceverino,
Colin DeGraf,
Sharon Lapiner,
Nir Mandelker,
Joel R. Primack
Abstract:
Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the…
▽ More
Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.
△ Less
Submitted 8 February, 2016; v1 submitted 8 September, 2015;
originally announced September 2015.
-
Evolution of Density Profiles in High-z Galaxies: Compaction and Quenching Inside-Out
Authors:
Sandro Tacchella,
Avishai Dekel,
C. Marcella Carollo,
Daniel Ceverino,
Colin DeGraf,
Sharon Lapiner,
Nir Mandelker,
Joel R. Primack
Abstract:
Using cosmological simulations, we address the interplay between structure and star formation in high-redshift galaxies via the evolution of surface density profiles. Our sample consists of 26 galaxies evolving in the redshift range $z=7-1$, spanning the stellar mass range $(0.2-6.4)\times 10^{10}M_\odot$ at $z=2$. We recover the main trends by stacking the profiles in accordance to their evolutio…
▽ More
Using cosmological simulations, we address the interplay between structure and star formation in high-redshift galaxies via the evolution of surface density profiles. Our sample consists of 26 galaxies evolving in the redshift range $z=7-1$, spanning the stellar mass range $(0.2-6.4)\times 10^{10}M_\odot$ at $z=2$. We recover the main trends by stacking the profiles in accordance to their evolution phases. Following a wet compaction event that typically occurs when the stellar mass is $\sim10^{9.5}~M_{\odot}$ at $z\sim2-4$, the gas develops a cusp inside the effective radius, associated with a peak in star-formation rate (SFR). The SFR peak and the associated feedback, in the absence of further gas inflow to the centre, marks the onset of gas depletion from the central 1 kpc, leading to quenching of the central SFR. An extended, star-forming ring that forms by fresh gas during the central quenching process shows as a rising specific SFR (sSFR) profile, which is interpreted as inside-out quenching. Before quenching, the stellar density profile grows self-similarly, maintaining its log-log shape because the sSFR is similar at all radii. During the quenching process, the stellar density saturates to a constant value, especially in the inner 1 kpc. The stellar mass and SFR profiles deduced from observations show very similar shapes, consistent with the scenario of wet compaction leading to inside-out quenching and the subsequent saturation of a dense stellar core. We predict a cuspy gas profile during the blue nugget phase, and a gas-depleted core, sometimes surrounded by a ring, in the post-blue nugget phase.
△ Less
Submitted 21 February, 2016; v1 submitted 31 August, 2015;
originally announced September 2015.
-
Compaction and Quenching of High-z Galaxies in Cosmological Simulations: Blue and Red Nuggets
Authors:
Adi Zolotov,
Avishai Dekel,
Nir Mandelker,
Dylan Tweed,
Shigeki Inoue,
Colin DeGraf,
Daniel Ceverino,
Joel Primack,
Guillermo Barro,
Sandra M. Faber
Abstract:
We use cosmological simulations to study a characteristic evolution pattern of high redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (blue nuggets) at z~4-2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z~2. The…
▽ More
We use cosmological simulations to study a characteristic evolution pattern of high redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (blue nuggets) at z~4-2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z~2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star-formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc ($Σ_1$). Massive galaxies quench earlier, faster, and at a higher $Σ_1$ than lower-mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-radius-mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum $Σ_1$ anti-correlated with the initial spin parameter, as predicted by Dekel & Burkert (2014). The central quenching is triggered by the high SFR and stellar/supernova feedback (possibly also AGN feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing.
△ Less
Submitted 6 March, 2015; v1 submitted 15 December, 2014;
originally announced December 2014.
-
Scaling relations between black holes and their host galaxies: comparing theoretical and observational measurements, and the impact of selection effects
Authors:
Colin DeGraf,
Tiziana Di Matteo,
Tommaso Treu,
Yu Feng,
Jong-Hak Woo,
Daeseong Park
Abstract:
We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental…
▽ More
We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental of the three relations. We find Gaussian scatter about all three relations, except among the highest mass galaxies, which host more massive black holes. Below z~2 the slopes for the full population remain roughly z-independent, and only steepen by 50% by z~4. The normalization of the sigma, LV relations evolve by 0.3, 0.43 dex, while the MBH correlation does not evolve to at least z~2. Testing for selection biases, we find samples selected by MBH or M* have steeper slopes than randomly selected samples. If unaccounted for, such a selection function would find faster evolution than inferred from a randomly selected sample, as objects at the highend of the relation tend to evolve more rapidly. We find a potential bias among high-LBH subsamples (tending to reside in higher mass galaxies), but these bright-AGN exhibit no intrinsic bias relative to fainter ones in equivalent-mass hosts, nor is there a significant difference between active- and inactive-samples. Finally we characterize the evolution of individual black holes along the scaling planes. Below the local relation, black holes grow faster than their host (72% of black holes 0.3 dex below the mean relation have a MBH-M* trajectory steeper than the local relation), while those above have shallower trajectories (only 14% are steeper than local). Black holes tend to grow faster than their hosts until surpassing the local relation, at which point their growth is suppressed while their hosts continue to grow, returning them to the mean relation.
△ Less
Submitted 19 December, 2014; v1 submitted 12 December, 2014;
originally announced December 2014.
-
Black hole growth and AGN feedback under clumpy accretion
Authors:
Colin DeGraf,
Avishai Dekel,
Jared Gabor,
Frederic Bournaud
Abstract:
High-resolution simulations of supermassive black holes in isolated galaxies have suggested the importance of short (~10 Myr) episodes of rapid accretion caused by interactions between the black hole and massive dense clouds within the host. Accretion of such clouds could potentially provide the dominant source for black hole growth in high-z galaxies, but it remains unresolved in cosmological sim…
▽ More
High-resolution simulations of supermassive black holes in isolated galaxies have suggested the importance of short (~10 Myr) episodes of rapid accretion caused by interactions between the black hole and massive dense clouds within the host. Accretion of such clouds could potentially provide the dominant source for black hole growth in high-z galaxies, but it remains unresolved in cosmological simulations. Using a stochastic subgrid model calibrated by high-resolution isolated galaxy simulations, we investigate the impact that variability in black hole accretion rates has on black hole growth and the evolution of the host galaxy. We find this clumpy accretion to more efficiently fuel high-redshift black hole growth. This increased mass allows for more rapid accretion even in the absence of high-density clumps, compounding the effect and resulting in substantially faster overall black hole growth. This increased growth allows the black hole to efficiently evacuate gas from the central region of the galaxy, driving strong winds up to ~2500 km/s, producing outflows ~10x stronger than the smooth accretion case, suppressing the inflow of gas onto the host galaxy, and suppressing the star formation within the galaxy by as much as a factor of two. This suggests that the proper incorporation of variability is a key factor in the co-evolution between black holes and their hosts.
△ Less
Submitted 26 April, 2017; v1 submitted 11 December, 2014;
originally announced December 2014.
-
The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0
Authors:
Nishikanta Khandai,
Tiziana Di Matteo,
Rupert Croft,
Stephen M. Wilkins,
Yu Feng,
Evan Tucker,
Colin DeGraf,
Mao-Sheng Liu
Abstract:
(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 pro…
▽ More
(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.
△ Less
Submitted 4 February, 2014;
originally announced February 2014.
-
The z=5 Quasar Luminosity Function from SDSS Stripe 82
Authors:
Ian D. McGreer,
Linhua Jiang,
Xiaohui Fan,
Gordon T. Richards,
Michael A. Strauss,
Nicholas P. Ross,
Martin White,
Yue Shen,
Donald P. Schneider,
Adam D. Myers,
W. Niel Brandt,
Colin DeGraf,
Eilat Glikman,
Jian Ge,
Alina Streblyanska
Abstract:
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…
▽ More
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.
△ Less
Submitted 22 April, 2013; v1 submitted 18 December, 2012;
originally announced December 2012.
-
The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine
Authors:
Nicholas P. Ross,
Ian D. McGreer,
Martin White,
Gordon T. Richards,
Adam D. Myers,
Nathalie Palanque-Delabrouille,
Michael A. Strauss,
Scott F. Anderson,
Yue Shen,
W. N. Brandt,
Christophe Yeche,
Molly E. C. Swanson,
Eric Aubourg,
Stephen Bailey,
Dmitry Bizyaev,
Jo Bovy,
Howard Brewington,
J. Brinkmann,
Colin DeGraf,
Tiziana Di Matteo,
Garrett Ebelke,
Xiaohui Fan,
Jian Ge,
Elena Malanushenko,
Viktor Malanushenko
, et al. (17 additional authors not shown)
Abstract:
We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the lu…
▽ More
We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.
△ Less
Submitted 23 October, 2012;
originally announced October 2012.
-
Growth of Early Supermassive Black Holes and the High-Redshift Eddington Ratio Distribution
Authors:
Colin DeGraf,
Tiziana Di Matteo,
Nishikanta Khandai,
Rupert Croft
Abstract:
Using a new large-scale (~ 0.75 Gpc)^3 hydrodynamic cosmological simulation we investigate the growth rate of supermassive black holes in the early universe (z > 4.75). Remarkably, we find a clear peak in the typical Eddington ratio at black hole masses of 4-8 * 10^7 solar masses (typically found in halos of ~7 * 10^11 to 10^12 solar masses), independent of redshift and indicative that most of BH…
▽ More
Using a new large-scale (~ 0.75 Gpc)^3 hydrodynamic cosmological simulation we investigate the growth rate of supermassive black holes in the early universe (z > 4.75). Remarkably, we find a clear peak in the typical Eddington ratio at black hole masses of 4-8 * 10^7 solar masses (typically found in halos of ~7 * 10^11 to 10^12 solar masses), independent of redshift and indicative that most of BH growth occurs in the cold-flow dominated regime. Black hole growth is by and large regulated by the evolution of gas density. The typical Eddington ratio at a given mass scales simply as cosmological density (1+z)^3 and the peak is caused by the competition between increased gas density available in more massive hosts, and a decrease due to strong AGN feedback that deprives the black hole of sufficient gas to fuel further rapid growth in the high mass end. In addition to evolution in the mean Eddington ratio, we show that the distribution of Eddington ratio among both mass-selected and luminosity-selected samples is approximately log-normal. We combine these findings into a single log-normal fitting formula for the distribution of Eddington ratios as a function of (M_BH,z). This formula can be used in analytic and semi analytic models for evolving black hole populations, predicting black hole masses of observed quasars, and, in conjunction with the observed distribution of Eddington ratios, can be used to constrain the black hole mass function.
△ Less
Submitted 25 January, 2012;
originally announced January 2012.
-
The Formation of Galaxies Hosting z~6 Quasars
Authors:
Nishikanta Khandai,
Yu Feng,
Colin DeGraf,
Tiziana Di Matteo,
Rupert A. C. Croft
Abstract:
We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simu…
▽ More
We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simulation are compact gas rich systems with high star formations rates of SFR ~ 100-1000 Msun/yr consistent with observed properties of Sloan quasar hosts in the redshift range 5.5 < z < 6.5. We show that the star-forming gas in these galaxies predominantly originates from high density cold streams which are able to penetrate the halo and grow the galaxy at the center. MassiveBlack predicts a deviation from the local Mbh-sigma and Mbh-Mstar relation implying that black holes are relatively more massive for a given stellar host at these redshifts.
△ Less
Submitted 2 November, 2011;
originally announced November 2011.
-
Early Black Holes in Cosmological Simulations: Luminosity Functions and Clustering Behaviour
Authors:
Colin DeGraf,
Tiziana Di Matteo,
Nishikanta Khandai,
Rupert Croft,
Julio Lopez,
Volker Springel
Abstract:
We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a reds…
▽ More
We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a redshift-independent occupation distribution of BHs among dark matter halos (which we provide) such that the evolution of the QLF follows that of the halo mass function. The sole exception is the bright-end at z=6 and 7, where BHs in high-mass halos tend to be unusually bright due to extended periods of Eddington growth caused by high density cold flows into the halo center. We further use these luminosity functions to make predictions for the number density of quasars in upcoming surveys, predicting there should be ~119+-28 (~87+-28) quasars detectable in the F125W band of the WIDE (DEEP) fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) from z=5-6, ~19+-7 (~18+-9) from z=6-7, and ~1.7+-1.5 (~1.5+-1.5) from z=7-8. We also investigate quasar clustering, finding that the correlation length is fully consistent with current constraints for Sloan quasars (r_0~17 h^{-1} Mpc at z=4 for quasars above m_i = 20.2), and grows slowly with redshift up to z=6 (r_0~22 h^{-1} Mpc). Finally, we note that the quasar clustering strength depends weakly on luminosity for low L_BH, but gets stronger at higher L_BH as the BHs are found in higher mass halos.
△ Less
Submitted 6 July, 2011;
originally announced July 2011.
-
Cold flows and the first quasars
Authors:
Tiziana Di Matteo,
Nishikanta Khandai,
Colin DeGraf,
Yu Feng,
Rupert Croft,
Julio Lopez,
Volker Springel
Abstract:
Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing questio…
▽ More
Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing question whether and how such high black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new cosmological hydrodynamic simulation (MassiveBlack) covering a volume (0.75 \Gpc)^3 appropriate for studying the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensurate with the existence of ~10^9 solar mass black holes as early as z~7. We find that under these conditions quasar feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench the accretion until the host galaxy reaches M_halo > 10^{12} solar masses. This cold-flow driven scenario for the formation of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs.
△ Less
Submitted 6 July, 2011;
originally announced July 2011.
-
The Halo Occupation Distribution of Active Galactic Nuclei
Authors:
Suchetana Chatterjee,
Colin DeGraf,
Jonathan Richardson,
Zheng Zheng,
Daisuke Nagai,
Tiziana Di Matteo
Abstract:
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 int…
▽ More
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.
△ Less
Submitted 28 September, 2011; v1 submitted 18 April, 2011;
originally announced April 2011.
-
The Halo Occupation Distribution of Black Holes: Dependence on Mass
Authors:
Colin DeGraf,
Matthew Oborski,
Tiziana Di Matteo,
Suchetana Chatterjee,
Daisuke Nagai,
Zheng Zheng,
Jonathan Richardson
Abstract:
We investigate the halo occupation distribution (HOD) of black holes within a hydrodynamic cosmological simulation that directly follows black hole growth. Similar to the HOD of galaxies/subhalos, we find that the black hole occupation number can be described by the form N_BH proportional to 1+ (M_Host)^alpha where alpha evolves mildly with redshift indicating that a given mass halo (M_Host) at lo…
▽ More
We investigate the halo occupation distribution (HOD) of black holes within a hydrodynamic cosmological simulation that directly follows black hole growth. Similar to the HOD of galaxies/subhalos, we find that the black hole occupation number can be described by the form N_BH proportional to 1+ (M_Host)^alpha where alpha evolves mildly with redshift indicating that a given mass halo (M_Host) at low redshift tends to host fewer BHs than at high redshift (as expected as a result of galaxy and BH mergers). We further divide the occupation number into contributions from black holes residing in central and satellite galaxies within a halo. The distribution of M_BH within halos tends to consist of a single massive BH (distributed about a peak mass strongly correlated with M_Host), and a collection of relatively low-mass secondary BHs, with weaker correlation with M_Host. We also examine the spatial distribution of BHs within their host halos, and find they typically follow a power-law radial distribution (i.e. much more centrally concentrated than the subhalo distribution). Finally, we characterize the host mass for which BH growth is feedback dominated (e.g. star formation quenched). We show that halos with M_Host > 3 * 10^12 M_sun have primary BHs that are feedback dominated by z~3 with lower mass halos becoming increasingly more affected at lower redshift.
△ Less
Submitted 7 February, 2011;
originally announced February 2011.
-
Quasar Clustering in Cosmological Hydrodynamic Simulations: Evidence for mergers
Authors:
Colin DeGraf,
Tiziana Di Matteo,
Volker Springel
Abstract:
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 fro…
▽ More
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.
△ Less
Submitted 10 May, 2010;
originally announced May 2010.
-
Faint-end Quasar Luminosity Functions from Cosmological Hydrodynamic Simulations
Authors:
Colin DeGraf,
Tiziana Di Matteo,
Volker Springel
Abstract:
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 si…
▽ More
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.
△ Less
Submitted 9 October, 2009;
originally announced October 2009.