We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z$ ≳ 6 quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in the IllustrisTNG simulations. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim 10^{12}\, {\rm M}_{\odot }\,h^{-1}$ haloes by $z$ = 7) using a (9 Mpc h−1)3 simulated volume. BH growth is maximal within haloes that are compact and have a low tidal field. For these haloes, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4\!-\!10^6\, {\rm M}_{\odot }\,h^{-1}$ seeds are inserted in haloes above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of Mseed). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z$ ∼ 6 quasar with $\sim 10^9\, {\rm M}_{\odot }$ mass and ∼1047 erg s−1 luminosity. BH mergers play a crucial role at $z$ ≳ 9, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. With more stringent seeding conditions (e.g. $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z$ ≳ 6 quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors ≳10) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results demonstrate that progenitors of $z$ ∼ 6 quasars have distinct BH merger histories for different seeding models, which will be distinguishable with Laser Interferometer Space Antenna observations.
We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z ∼ 4–7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ∼1″ to $\lesssim \!\!{0{^{\prime \prime }_{.}}1}$, which is comparable to that of the Hubble Space Telescope, allowing us to identify z ∼ 4–7 bright major mergers at a high completeness value of $\gtrsim \!\!90\%$. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z ∼ 4, 5, 6, and 7, respectively, in a UV luminosity range of LUV ∼ 3–$15\, L_{\rm UV}^*$ corresponding to −24 ≲ MUV ≲ −22. The major merger fractions are estimated to be $f_{\rm merger}\sim 10\%$–$20\%$ at z ∼ 4 and $\sim 50\%$–$70\%$ at z ∼ 5–7, which shows no fmerger difference compared to those of a control faint galaxy sample. Based on the fmerger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z ∼ 4–7 galaxy UV LFs. While these two effects partly explain the DPL shape at LUV ∼ 3–$10\, L_{\rm UV}^*$, the DPL shape cannot be explained at the very bright end of $L_{\rm UV}\gtrsim 10\, L_{\rm UV}^*$, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs.
more » « less- NSF-PAR ID:
- 10363666
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Publications of the Astronomical Society of Japan
- Volume:
- 74
- Issue:
- 1
- ISSN:
- 0004-6264
- Format(s):
- Medium: X Size: p. 73-91
- Size(s):
- ["p. 73-91"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
null (Ed.)ABSTRACT We study the projected spatial offset between the ultraviolet continuum and Ly α emission for 65 lensed and unlensed galaxies in the Epoch of Reionization (5 ≤ z ≤ 7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly α properties of galaxies observed in slit spectroscopy. While we find that ∼40 per cent of galaxies in our sample show significant projected spatial offsets ($|\Delta _{\rm {Ly}\alpha -\rm {UV}}|$), we find a relatively modest average projected offset of $|\widetilde{\Delta }_{\rm {Ly}\alpha -\rm {UV}}|$ = 0.61 ± 0.08 proper kpc for the entire sample. A small fraction of our sample, ∼10 per cent, exhibit offsets in excess of 2 proper kpc, with offsets seen up to ∼4 proper kpc, sizes that are considerably larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta _{\rm {Ly}\alpha -\rm {UV}}|$ with redshift. In our sample, ultraviolet (UV)-bright galaxies ($\widetilde{L_{\mathrm{ UV}}}/L^{\ast }_{\mathrm{ UV}}=0.67$) showed offsets a factor of three greater than their fainter counterparts ($\widetilde{L_{\mathrm{ UV}}}/L^{\ast }_{\mathrm{ UV}}=0.10$), 0.89 ± 0.18 versus 0.27 ± 0.05 proper kpc, respectively. The presence of companion galaxies and early stage merging activity appeared to be unlikely causes of these offsets. Rather, these offsets appear consistent with a scenario in which internal anisotropic processes resulting from stellar feedback, which is stronger in UV-brighter galaxies, facilitate Ly α fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly α flux due to offsets was quantified. It was found that the differential loss of Ly α photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4 arcsec). However, for the largest offsets, if they are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of ≤1 arcsec were employed, such as those planned for James Webb Space Telescope/NIRSpec observations.more » « less
-
more » « less
ABSTRACT Direct collapse black holes (BHs) are promising candidates for producing massive z ≳ 6 quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: (1) low gas angular momentum (spin), and (2) a minimum incident Lyman–Werner (LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4\!-\!10^6\, \mathrm{M}_{\odot }\, h^{-1})$ in haloes with a total mass >3000 × Mseed and a dense, metal-poor gas mass >5 × Mseed. Within this framework, we find that the seed-forming haloes have a prior history of star formation and metal enrichment, but they also contain pockets of dense, metal-poor gas. When seeding is further restricted to haloes with low gas spins, the number of seeds formed is suppressed by factors of ∼6 compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal-poor gas is required to have a critical LW flux (Jcrit). Even for Jcrit values as low as 50J21, no $8\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$ seeds are formed. While lower mass ($1.25\times 10^{4},1\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$) seeds do form, they are strongly suppressed (by factors of ∼10–100) compared to the baseline model at gas mass resolutions of $\sim 10^4~\mathrm{M}_{\odot }\, h^{-1}$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, none of the seeds are able to grow to the supermassive regime ($\gtrsim 10^6~\mathrm{M}_{\odot }\, h^{-1}$) by z = 7. Our results hint that producing the bulk of the z ≳ 6 supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities.
-
more » « less
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 an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the 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{Myr}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $\lesssim 20{{\rm per \,cent}}$ 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 $\gt 10^9\, {\rm 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 ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $\gtrsim 60\,{\rm{per\,cent}}$ seed-seed mergers, $\sim 30\,{\rm{per\,cent}}$ involving only one seed-mass MBH, and $\sim 10\,{\rm{per\,cent}}$ mergers of non-seed MBHs.
-
ABSTRACT In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simulations, at any redshift.more » « less
