An official website of the United States government
ABSTRACT While the first “seeds” of supermassive black holes (BH) can range from $$\sim 10^2-10^6 \rm ~{\rm M}_{\odot }$$, the lowest mass seeds ($$\lesssim 10^3~\rm {\rm M}_{\odot }$$) are inaccessible to most cosmological simulations due to resolution limitations. We present our new BRAHMA simulations that use a novel flexible seeding approach to predict the $$z\ge 7$$ BH populations for low-mass seeds. We ran two types of boxes that model $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds using two distinct but mutually consistent seeding prescriptions at different simulation resolutions. First, we have the highest resolution $$[9~\mathrm{Mpc}]^3$$ (BRAHMA-9-D3) boxes that directly resolve $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds and place them within haloes with dense, metal-poor gas. Second, we have lower resolution, larger volume $$[18~\mathrm{Mpc}]^3$$ (BRAHMA-18-E4), and $$\sim [36~\mathrm{Mpc}]^3$$ (BRAHMA-36-E5) boxes that seed their smallest resolvable $$\sim 10^4~\&~10^5~\mathrm{{\rm M}_{\odot }}$$ BH descendants using new stochastic seeding prescriptions calibrated using BRAHMA-9-D3. The three boxes together probe key BH observables between $$\sim 10^3\,\mathrm{ and}\,10^7~\rm {\rm M}_{\odot }$$. The active galactic nuclei (AGN) luminosity function variations are small (factors of $$\sim 2-3$$) at the anticipated detection limits of potential future X-ray facilities ($$\sim 10^{43}~ \mathrm{ergs~s^{-1}}$$ at $$z\sim 7$$). Our simulations predict BHs $$\sim 10-100$$ times heavier than the local $$M_*$$ versus $$M_{\mathrm{ bh}}$$ relations, consistent with several JWST-detected AGN. For different seed models, our simulations merge binaries at $$\sim 1-15~\mathrm{kpc}$$, with rates of $$\sim 200-2000$$ yr−1 for $$\gtrsim 10^3~\rm {\rm M}_{\odot }$$ BHs, $$\sim 6-60$$ yr−1 for $$\gtrsim 10^4~\rm {\rm M}_{\odot }$$ BHs, and up to $$\sim 10$$ yr−1 amongst $$\gtrsim 10^5~\rm {\rm M}_{\odot }$$ BHs. These results suggest that Laser Interferometer Space Antenna mission has promising prospects for constraining seed models.
more »
« less
The mass-loss rates of star clusters with stellar-mass black holes: implications for the globular cluster mass function
ABSTRACT Stellar-mass black holes (BHs) can be retained in globular clusters (GCs) until the present. Simulations of GC evolution find that the relaxation driven mass-loss rate is elevated if BHs are present, especially near dissolution. We capture this behaviour in a parametrized mass-loss rate, bench marked by results from N-body simulations, and use it to evolve an initial GC mass function (GCMF), similar to that of young massive clusters in the Local Universe, to an age of 12 Gyr. Low-metallicity GCs ([Fe/H] ≲ −1.5) have the highest mass-loss rates, because of their relatively high BH masses, which combined with their more radial orbits and stronger tidal field in the past explains the high turnover mass of the GCMF ($$\sim 10^5\, {\rm M}_\odot$$ ) at large Galactic radii ($$\gtrsim 10\, {\rm kpc}$$ ). The turnover mass at smaller Galactic radii is similar because of the upper mass truncation of the initial GCMF and the lower mass-loss rate due to the higher metallicities. The density profile in the Galaxy of mass lost from massive GCs ($$\gtrsim 10^{5}\, {\rm M}_\odot$$ ) resembles that of nitrogen-rich stars in the halo, confirming that these stars originated from GCs. We conclude that two-body relaxation is the dominant effect in shaping the GCMF from a universal initial GCMF, because including the effect of BHs reduces the need for additional disruption mechanisms.
more »
« less
- Award ID(s):
- 1909063
- PAR ID:
- 10413883
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 522
- Issue:
- 4
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 5340-5357
- Size(s):
- p. 5340-5357
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $$M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($$M_{*}\lesssim 10^7\, \rm {M}_{\odot }$$) are mostly quiescent and higher mass satellites ($$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$) are mostly star forming, with intermediate-mass satellites ($$M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($$M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $$M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$, they do not generally reproduce SAGA’s turnover at lower masses.more » « less
-
ABSTRACT The origin of the ‘seeds’ of supermassive black holes (BHs) continues to be a puzzle, as it is currently unclear if the imprints of early seed formation could survive to today. We examine the signatures of seeding in the local Universe using five $$[18~\mathrm{Mpc}]^3$$BRAHMA simulation boxes run to $z=0$. They initialize $$1.5\times 10^5~\rm {M}_{\odot }$$ BHs using different seeding models. The first four boxes initialize BHs as heavy seeds using criteria that depend on dense and metal-poor gas, Lyman–Werner radiation, gas spin, and environmental richness. The fifth box initializes BHs as descendants of lower mass seeds ($$\sim 10^3~\rm {M}_{\odot }$$) using a new stochastic seed model built in our previous work. In our simulations, we find that the abundances and properties of $$\sim 10^5-10^6~\rm {M}_{\odot }$$ local BHs hosted in $$M_*\lesssim 10^{9}~\rm {M}_{\odot }$$ dwarf galaxies, are sensitive to the assumed seeding criteria. This is for two reasons: (1) there is a substantial population of local $$\sim 10^5~\rm {M}_{\odot }$$ BHs that are ungrown relics of early seeds from $$z\sim 5-10$$; (2) BH growth up to $$\sim 10^6~\rm {M}_{\odot }$$ is dominated by mergers in our simulations all the way down to $$z\sim 0$$. As the contribution from gas accretion increases, the signatures of seeding start to weaken in more massive $$\gtrsim 10^6~\rm {M}_{\odot }$$ BHs, and they are erased for $$\gtrsim 10^7~\rm {M}_{\odot }$$ BHs. The different seed models explored here predict abundances of local $$\sim 10^6~\rm {M}_{\odot }$$ BHs ranging from $$\sim 0.01-0.05~\mathrm{Mpc}^{-3}$$ with occupation fractions of $$\sim 20-100~{{\ \rm per\ cent}}$$ for $$M_*\sim 10^{9}~\rm {M}_{\odot }$$ galaxies. These results highlight the potential for placing constraints on seeding models using local $$\sim 10^5-10^6~\rm {M}_{\odot }$$ BHs hosted in dwarf galaxies. Since merger dynamics and accretion physics impact the persistence of seeding signatures, and both high and low mass seed models can produce similar local BH populations, disentangling their roles will require combining high and low redshift constraints.more » « less
-
ABSTRACT We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($$M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$$), run to z = 0 with $$30\, \mathrm{M}_{\odot }$$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $$M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $$M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($$M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core $${\gt}200\rm \, pc$$, while lower mass UFDs exhibit cusps down to $${\lesssim}100\rm \, pc$$, as expected from energetic arguments. Our dwarfs with $$M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).more » « less
-
ABSTRACT The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMSs, M>100 $$\rm {\rm M}_{\odot }$$) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50–500 $$\rm {\rm M}_{\odot }$$, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne–Na and Mg–Al cycles. VMS with enhanced winds eject 5–10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anticorrelations, such as C–N and Na–O, in globular clusters. We find that for VMS 95 per cent of the total wind yields is produced on the main sequence, while only ∼ 5 per cent is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf–Rayet winds had been suggested to be responsible for galactic 26Al enrichment. Finally, 200 $$\rm {\rm M}_{\odot }$$ stars eject 100 times more of each heavy element in their winds than 50 $$\rm {\rm M}_{\odot }$$ stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50 $$\rm {\rm M}_{\odot }$$ stars.more » « less
