Cosmological Hydrodynamical Simulations Research Articles

The SRG/eROSITA All-Sky Survey

Context. Galaxy groups lying between galaxies and galaxy clusters in the mass spectrum of dark matter halos play a crucial role in the evolution and formation of the large-scale structure. Their shallower potential wells compared to clusters of galaxies make them excellent sources to constrain non-gravitational processes such as feedback from the central active galactic nuclei (AGN). Aims. We investigate the impact of feedback, particularly from AGN, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA’s first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms and help guide future AGN feedback implementations in numerical simulations. Methods. We analyzed the deeper eROSITA observations of 1178 galaxy groups detected in the eRASS1. We divided the sample into 271 subsamples based on their physical and statistical properties and extracted average thermodynamic properties, including the electron number density, temperature, and entropy, at three characteristic radii from cores to outskirts along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. Results. We present the tightest constraints with unprecedented statistical precision on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium (IGrM) temperatures, while for the low-mass groups with cooler (T < 1.44 keV) IGrM temperatures, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. Additionally, comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL) after applying the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Conclusions. These measurements will pave the way for achieving more realistic AGN feedback implementations in numerical simulations. The future eROSITA Surveys will enable the extension of the entropy measurements in even cooler IGrM temperatures below 0. 5 keV, allowing for the testing of the AGN feedback models in this regime.

SAGAbg. II. The Low-mass Star-forming Sequence Evolves Significantly between 0.05 < z < 0.21

The redshift-dependent relation between galaxy stellar mass and star formation rate (SFR), known as the star-forming sequence (SFS), is a key observational yardstick for galaxy assembly. We use the SAGAbg-A sample of background galaxies from the Satellites Around Galactic Analogs (SAGA) Survey to model the low-redshift evolution of the low-mass SFS. The sample is comprised of 23,258 galaxies with Hα-based SFRs spanning 6<log10(M⋆/[M⊙])<10 and z < 0.21 (t < 2.5 Gyr). Although it is common to bin or stack galaxies at z ≲ 0.2 for galaxy population studies, the difference in lookback time between z = 0 and z = 0.21 is comparable to the time between z = 1 and z = 2. We develop a model to account for both the physical evolution of low-mass SFS and the selection function of the SAGA Survey, allowing us to disentangle redshift evolution from redshift-dependent selection effects across the SAGAbg-A redshift range. Our findings indicate significant evolution in the SFS over the last ∼2.5 Gyr, with a rising normalization: 〈SFR(M⋆=108.5M⊙)〉(z)=1.24−0.23+0.25z−1.47−0.03+0.03 . We also identify the redshift limit at which a static SFS is ruled out at the 95% confidence level, which is z = 0.05 based on the precision of the SAGAbg-A sample. Comparison with cosmological hydrodynamic simulations reveals that some contemporary simulations underpredict the recent evolution of the low-mass SFS. This demonstrates that the recent evolution of the low-mass SFS can provide new constraints on the assembly of the low-mass Universe and highlights the need for improved models in this regime.

Sequential formation of supermassive stars and heavy seed BHs through the interplay of cosmological cold accretion and stellar radiative feedback

ABSTRACT Supermassive stars (SMSs) and heavy seed black holes, as their remnants, are promising candidates for supermassive black hole (SMBH) progenitors, especially for ones observed in the early universe $z\simeq 8.5-10$ by recent JWST observations. Expected cradles of SMSs are the atomic cooling haloes ($M_{\rm halo}\simeq 10^7\ \mathrm{M}_{\odot }$), where ‘cold accretion’ emerges and possibly forms SMSs. We perform a suit of cosmological radiation hydrodynamics simulations and investigate star formation after the emergence of cold accretion, solving radiative feedback from stars inside the halo. We follow the mass growth of the protostars for $\sim 3\ \mathrm{Myr}$, resolving the gas inflow down to $\sim 0.1\ \mathrm{pc}$ scales. We discover that, after cold accretion emerges, multiple SMSs of $m_{\star }\gtrsim 10^5\ \mathrm{M}_{\odot }$ form at the halo centre with the accretion rates maintained at $\dot{m}_{\star }\simeq 0.04\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ for $\lesssim 3\ \mathrm{Myr}$. Cold accretion supplies gas at a rate of $\dot{M}_{\rm gas}\gtrsim 0.01-0.1\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ from outside the halo virial radius to the central gas disc. Gravitational torques from spiral arms transport gas further inwards, which feeds the SMSs. Radiative feedback from stars suppresses $\mathrm{H}_2$ cooling and disc fragmentation, while photoevaporation is prevented by a dense envelope, which attenuates ionizing radiation. Our results suggest that cold accretion can bring efficient BH mass growth after seed formation in the later universe. Moreover, cold accretion and gas migration inside the central disc increase the mass concentration and provide a promising formation site for the extremely compact stellar clusters observed by JWST.

Slick: Modeling a Universe of Molecular Line Luminosities in Hydrodynamical Simulations

We present slick (the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [C i], and [C ii] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer code despotic, slick computes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their physical properties and individual environments. We validate our results by applying slick to the high-resolution run of the Simba simulations, testing the derived luminosities against empirical and theoretical/analytic relations. To simulate the line emission from a universe of emitting clouds, we have incorporated random forest machine learning (ML) methods into our approach, allowing us to predict cosmologically evolving properties of CO, [C i], and [C ii] emission from galaxies such as luminosity functions. We tested this model in 100,000 gas particles, and 2500 galaxies, reaching an average accuracy of ∼99.8% for all lines. Finally, we present the first model light cones created with realistic and ML-predicted CO, [C i], and [C ii] luminosities in cosmological hydrodynamical simulations, from z = 0 to z = 10.

Galaxy shapes in Magneticum

Context. Despite being one of the most fundamental properties of galaxies that dictate the form of the potential, the 3D shapes are intrinsically difficult to determine from observations. The improving quality of triaxial modeling methods in recent years has made it possible to measure these shapes more accurately. Aims. This study provides a comprehensive understanding of the stellar and dark matter (DM) shapes of galaxies and the connection between them. As these shapes are the result of the formation history of a galaxy, we investigate which galaxy properties they are correlated with, which will be especially useful for interpreting the results from dynamical modeling. Methods. Using the hydrodynamical cosmological simulation Magneticum Pathfinder Box4 (uhr), we computed the stellar and DM intrinsic shapes of 690 simulated galaxies with stellar masses above 2 × 1010 M⊙ at three different radii with an iterative unweighted method. We also determined their morphologies, their projected morphological and kinematic parameters, and their fractions of in situ formed stars. Results. The DM follows the stellar component in shape and orientation at three half-mass radii, indicating that DM is heavily influenced by the baryonic potential in the inner parts of the halo. The outer DM halo is independent of the inner properties such as the DM shape or galaxy morphology, however, and is more closely related to the large-scale anisotropy of the gas inflow. Overall, DM halo shapes are prolate, consistent with previous literature. The stellar shapes of galaxies are correlated with their morphology, with early-type galaxies featuring more spherical and prolate shapes than late-type galaxies out to 3 R1/2. Galaxies with more rotational support are flatter, and the stellar shapes are connected to the mass distribution, though not to the mass itself. In particular, more extended early-type galaxies have larger triaxialities at a given mass. Finally, the shapes can be used to better constrain the in situ fraction of stars when combined with the stellar mass. Conclusions. The relations between shape, mass distribution, and in situ formed star fraction of galaxies show that the shapes depend on the details of the accretion history through which the galaxies are formed. The similarities between DM and stellar shapes in the inner regions of galaxy halos signal the importance of baryonic matter for the behavior of DM in galaxies and will be of use for improving the underlying assumptions of dynamical models for galaxies in the future. However, at large radii the shapes of the DM are completely decoupled from the central galaxy, and their shapes and spin are coupled more to the large scale inflow than to the galaxy in the center.

Tracing the History of Obscured Star Formation with the SIMBA Cosmological Galaxy Evolution Simulation

We explore the cosmic evolution of the fraction of dust-obscured star formation predicted by the simba cosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up to z = 3, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, at z > 3 we find strong evolution at fixed galaxy stellar mass toward greater amounts of obscured star formation, in tension with high-redshift observations. We explain the trend of increasing obscuration at higher redshifts by evolving star-dust geometry, as the dust-to-stellar mass ratios remain relatively constant across cosmic time. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust-to-stellar mass ratios at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust-obscured star formation history peaks around z ∼ 2−3, and becomes subdominant at z ≳ 5. The dominance of obscured star formation at redshifts z ≲ 4 is consistent with our results for the evolution of the obscured star formation fraction at fixed stellar mass to higher values at higher redshift because there exist fewer massive, heavily obscured galaxies at high redshift.

A calibrated model for N-body dynamical friction acting on supermassive black holes

Abstract Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is responsible for driving massive black holes towards the centres of galaxies, cannot be accurately modelled with softened N-body interactions. A number of subgrid models have been proposed to mimic dynamical friction or directly include its full effects in simulations. Each of these methods has its individual benefits and shortcomings, while all suffer from a common issue of being unable to represent black holes with masses below a few times the simulated dark matter particle mass. In this paper, we propose a correction for unresolved dynamical friction, which has been calibrated on simulations run with the code KETJU, in which gravitational interactions of black holes are not softened. We demonstrate that our correction is able to sink black holes with masses greater than the dark matter particle mass at the correct rate. We show that the impact of stochasticity is significant for low-mass black holes (MBH ≤ 5MDM) and propose a correction for stochastic heating. Combined, this approach is applicable to next generation cosmological hydrodynamics simulations that jointly track galaxy and black hole growth with realistic black hole orbits.

Testing Lyα Emitters and Lyman-break Galaxies as Tracers of Large-scale Structures at High Redshifts

We test whether Lyα emitters (LAEs) and Lyman-break galaxies (LBGs) can be good tracers of high-z large-scale structures, using the Horizon Run 5 cosmological hydrodynamical simulation. We identify LAEs using the Lyα emission line luminosity and its equivalent width, and LBGs using the broadband magnitudes at z ∼ 2.4, 3.1, and 4.5. We first compare the spatial distributions of LAEs, LBGs, all galaxies, and dark matter around the filamentary structures defined by dark matter. The comparison shows that both LAEs and LBGs are more concentrated toward the dark matter filaments than dark matter. We also find an empirical fitting formula for the vertical density profile of filaments as a binomial power-law relation of the distance to the filaments. We then compare the spatial distributions of the samples around the filaments defined by themselves. LAEs and LBGs are again more concentrated toward their filaments than dark matter. We also find the overall consistency between filamentary structures defined by LAEs, LBGs, and dark matter, with the median spatial offsets that are smaller than the mean separation of the sample. These results support the idea that the LAEs and LBGs could be good tracers of large-scale structures of dark matter at high redshifts.

CTA and SWGO can discover Higgsino dark matter annihilation

Thermal Higgsino dark matter (DM), with a mass near 1.1 TeV, is one of the most well-motivated and untested DM candidates. Leveraging recent hydrodynamic cosmological simulations that give DM density profiles in Milky Way analog galaxies we show that the linelike gamma-ray signal predicted from Higgsino annihilation in the Galactic Center could be detected at high significance with the upcoming Cherenkov Telescope Array (CTA) and Southern Wide-field Gamma-ray Observatory (SWGO) for all but the most pessimistic DM profiles. We perform the most sensitive search to-date for the linelike signal using 15 years of data from the Fermi Large Area Telescope, coming within an order one factor of the necessary sensitivity to detect the Higgsino for some Milky Way analog DM density profiles. We show that H.E.S.S. has subleading sensitivity relative to Fermi for the Higgsino at present. In contrast, we analyze H.E.S.S. inner Galaxy data for the thermal wino model with a mass near 2.8 TeV; we find no evidence for a DM signal and exclude the wino by over a factor of two in cross section for all DM profiles considered. In the process, we identify and attempt to correct what appears to be an inconsistency in previous H.E.S.S. inner Galaxy analyses for DM annihilation related to the analysis effective area, which may weaken the DM cross-section sensitivity claimed in those works by around an order of magnitude. Published by the American Physical Society 2024

The Negative Baryon Acoustic Oscillation Shift in the Lyα Forest from Cosmological Simulations

We present the first measurement of the Lyα forest baryon acoustic oscillations (BAO) shift parameter from cosmological simulations. In particular, we generate a suite of 1000 accurate effective field-level bias-based Lyα forest simulations of volume V=(1h−1Gpc)3 at z = 2, both in real and redshift space, calibrated upon two fixed-and-paired cosmological hydrodynamic simulations. To measure the BAO, we stack the 3D power spectra of the 1000 different realizations, compute the average, and use a model accounting for a proper smooth-peak component decomposition of the power spectrum, to fit it via an efficient Markov Chain Monte Carlo scheme estimating the covariance matrices directly from the simulations. We report the BAO shift parameters to be α=0.9969−0.0014+0.0014 and α=0.9905−0.0027+0.0027 in real and redshift space, respectively. We also measure the bias b lya and the BAO broadening parameter Σnl, finding blya=−0.1786−0.0001+0.0001 and Σnl=3.87−0.20+0.20 in real space, and blya=−0.073−0.004+0.005 and Σnl=6.55−0.22+0.23 in redshift space. Moreover, we measure the linear Kaiser factor βlya=1.39−0.18+0.24 from the isotropic redshift space fit. Overall, we find evidence for a negative shift of the BAO peak at the ∼2.2σ and ∼3.5σ levels in real and redshift space, respectively. This work sets new important theoretical constraints on the Lyα forest BAO scale and offers a potential solution to the tension emerging from previous observational analysis, in light of ongoing and upcoming Lyα forest spectroscopic surveys, such as DESI, the Prime Focus Spectrograph Survey, and WEAVE-QSO.

Can We Constrain Warm Dark Matter Masses with Individual Galaxies?

We study the impact of warm dark matter (WDM) mass on the internal properties of individual galaxies using a large suite of 1024 state-of-the-art cosmological hydrodynamic simulations from the DREAMS project. We take individual galaxies’ properties from the simulations, which have different cosmologies, astrophysics, and WDM masses, and train normalizing flows to learn the posterior of the parameters. We find that our models cannot infer the value of the WDM mass, even when the values of the cosmological and astrophysical parameters are given explicitly. This result holds for galaxies with stellar mass larger than 2 × 108 M ⊙ h −1 at both low and high redshifts. We calculate the mutual information and find no significant dependence between the WDM mass and galaxy properties. On the other hand, our models can infer the value of Ωm with a ∼10% accuracy from the properties of individual galaxies while marginalizing astrophysics and WDM masses.

A warm dark matter cosmogony may yield more low-mass galaxy detections in 21-cm surveys than a cold dark matter one

ABSTRACT The 21-cm spectral line widths, $w_{50}$, of galaxies are an approximate tracer of their dynamical masses, such that the dark matter halo mass function is imprinted in the number density of galaxies as a function of $w_{50}$. Correcting observed number counts for survey incompleteness at the level of accuracy needed to place competitive constraints on warm dark matter (WDM) cosmological models is very challenging, but forward-modelling the results of cosmological hydrodynamical galaxy formation simulations into observational data space is more straightforward. We take this approach to make predictions for an ALFALFA-like survey from simulations using the EAGLE galaxy formation model in both cold (CDM) and WDM cosmogonies. We find that for WDM cosmogonies more galaxies are detected at the low-$w_{50}$ end of the 21-cm velocity width function than in the CDM cosmogony, contrary to what might naïvely be expected from the suppression of power on small scales in such models. This is because low-mass galaxies form later and retain more gas in WDM cosmogonies (with EAGLE). While some shortcomings in the treatment of cold gas in the EAGLE model preclude placing definitive constraints on WDM scenarios, our analysis illustrates that near-future simulations with more accurate modelling of cold gas will likely make strong constraints possible, especially in conjunction with new 21-cm surveys such as WALLABY.

Why does the Milky Way have a metallicity floor?

ABSTRACT The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim 10^{-3.5}$ Z$_{\odot }$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihaloes present during the Galaxy’s initial formation. To arrive at this conclusion, we analyse a cosmological radiation hydrodynamics simulation that follows the concurrent evolution of multiple Population III star-forming minihaloes. The main driver for the central gas within minihaloes is the steady increase in hydrostatic pressure as the haloes grow. We incorporate this insight into a hybrid one-zone model that switches between pressure-confined and modified free-fall modes to evolve the gas density with time according to the ratio of the free-fall and sound-crossing time-scales. This model is able to accurately reproduce the density and chemo-thermal evolution of the gas in each of the simulated minihaloes up to the point of runaway collapse. We then use this model to investigate how the gas responds to the absence of H$_{2}$. Without metals, the central gas becomes increasingly stable against collapse as it grows to the atomic cooling limit. When metals are present in the halo at a level of $\sim 10^{-3.7}$ Z$_{\odot }$, however, the gas is able to achieve gravitational instability while still in the minihalo regime. Thus, we conclude that the Galaxy’s metallicity floor is set by the balance within minihaloes of gas-phase metal cooling and the radiation background associated with its early formation environment.

Correction to: The role of environment and AGN feedback in quenching local galaxies: comparing cosmological hydrodynamical simulations to the SDSS

Correction to: The role of environment and AGN feedback in quenching local galaxies: comparing cosmological hydrodynamical simulations to the SDSS

Star formation properties of z ∼ 1 galaxy clusters and groups from Horizon Run 5

ABSTRACT Quiescent galaxies are predominantly observed in local galaxy clusters. However, the fraction of quiescent galaxies in high-redshift clusters significantly varies among different clusters. In this study, we present the results of an analysis of the star formation (SF) properties of $z \sim 0.87$ clusters and groups from a cosmological hydrodynamical simulation Horizon Run 5. We investigate the correlation between the quiescent galaxy fraction (QF) of these model clusters/groups and their various internal or external properties. We find that halo mass is one of the most important characteristics as higher mass clusters and groups have higher QFs. We also find that other properties such as stellar-mass ratio and Friends-of-Friends fraction, which measures the proportion of the area around a cluster occupied by dense structures, may mildly affect the QFs of clusters and groups. This may indicate that the evolutionary history as well as the large-scale environment of clusters and groups also play a certain role in determining the SF status of high-redshift galaxy clusters and groups.

The baryon cycle in modern cosmological hydrodynamical simulations

ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx 0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2\!-\!3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $R_{\rm 200c}$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx 10^{13.5}\, \mathrm{M}_{\odot }$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2\!-\!3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform subgrid feedback models in the next generation of simulations.

The great escape: understanding the connection between Ly α emission and LyC escape in simulated JWST analogues

ABSTRACT Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly$\, \alpha$ emission line correlate with the inferred LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) of low-redshift galaxies. Using a data set of 9600 mock Ly$\, \alpha$ spectra of star-forming galaxies at $4.64 \le z \le 6$ from the SPHINX$^{20}$ cosmological radiation hydrodynamical simulation, we study the physics controlling the escape of Ly$\, \alpha$ and LyC photons. We find that our mock Ly$\, \alpha$ observations are representative of high-redshift observations and that typical observational methods tend to overpredict the Ly$\, \alpha$ escape fraction ($f_{\rm esc}^{\rm Ly\, \alpha }$) by as much as 2 dex. We investigate the correlations between $f_{\rm esc}^{\rm LyC}$ and $f_{\rm esc}^{\rm Ly\, \alpha }$, Ly$\, \alpha$ equivalent width ($W_{\lambda }({\rm Ly\, \alpha })$), peak separation ($v_{\rm sep}$), central escape fraction ($f_{\rm cen}$), and red peak asymmetry ($A_f^{\rm red}$). We find that $f_{\rm esc}^{\rm Ly\, \alpha }$ and $f_{\rm cen}$ are good diagnostics for LyC leakage, selecting for galaxies with lower neutral gas densities and less UV attenuation that have recently experienced supernova feedback. In contrast, $W_{\lambda }({\rm Ly\, \alpha })$ and $v_{\rm sep}$ are found to be necessary but insufficient diagnostics, while $A_f^{\rm red}$ carries little information. Finally, we use stacks of Ly$\, \alpha$, H$\, \alpha$, and F150W mock surface brightness profiles to find that galaxies with high $f_{\rm esc}^{\rm LyC}$ tend to have less extended Ly$\, \alpha$ and F150W haloes but larger H$\, \alpha$ haloes than their non-leaking counterparts. This confirms that Ly$\, \alpha$ spectral profiles and surface brightness morphology can be used to better understand the escape of LyC photons from galaxies during the epoch of reionization.

Black Hole–Halo Mass Relation from UNIONS Weak Lensing

This Letter presents, for the first time, direct constraints on the black hole–halo mass relation using weak gravitational-lensing measurements. We construct type I and type II active galactic nucleus (AGN) samples from the Sloan Digital Sky Survey, with a mean redshift of 0.4 (0.1) for type I (type II) AGNs. This sample is cross correlated with weak-lensing shear from the Ultraviolet Near Infrared Optical Northern Survey. We compute the excess surface mass density of the halos associated with 36,181 AGNs from 94,308,561 lensed galaxies and fit the halo mass in bins of black hole mass. We find that more massive AGNs reside in more massive halos. The relation between halo mass and black hole mass is well described by a power law of slope 0.6 for both type I and type II samples, in agreement with models that link black hole growth to baryon feedback. We see no dependence on AGN type or redshift in the black hole–halo mass relation below a black hole mass of 108.5 M ⊙. Above that mass, we find more massive halos for the low-z type II sample compared to the high-z type I sample, but this difference may be interpreted as systematic error in the black hole mass measurements. Our results are consistent with previous measurements for non-AGN galaxies. At a fixed black hole mass, our weak-lensing halo masses are consistent with galaxy rotation curves but significantly lower than galaxy-clustering measurements. Finally, our results are broadly consistent with state-of-the-art hydrodynamical cosmological simulations, providing a new constraint for black hole masses in simulations.

Violent Starbursts and Quiescence Induced by Far-ultraviolet Radiation Feedback in Metal-poor Galaxies at High Redshift

JWST observations of galaxies at z ≳ 8 suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one galaxy at subparsec resolution with different radiative feedback models. We find that the suppression of cooling by far-UV (FUV) radiation (i.e., H2 dissociating radiation) from Population II stars is the main physical process triggering the formation of compact and massive star clusters and is responsible for the bursty star formation observed in metal-poor galaxies at z ≳ 10. Indeed, artificially suppressing FUV radiation leads to a less intense continuous mode of star formation distributed into numerous but low-mass open star clusters. Due to the intense FUV field, low-metallicity clouds remain warm (∼104 K) until they reach a relatively high density (≳103 cm−3), before becoming self-shielded and transitioning to a colder (∼100 K), partially molecular phase. As a result, star formation is delayed until the clouds accumulate enough mass to become gravitationally unstable. At this point, the clouds undergo rapid star formation, converting gas into stars with high efficiency. We therefore observe exceptionally bright galaxies (10 times brighter than for continuous star formation) and subsequent quenched “dead” galaxies that did not form stars for tens of Myr.

What drives the corpulence of galaxies?

Nearby dwarf galaxies display a variety of effective radii (sizes) at a given stellar mass, suggesting different evolution scenarios according to their final “stellar” size. The TNG hydrodynamical simulations present a bimodality in the z = 0 size–mass relation (SMRz0) of dwarf galaxies, at r1/2, ⋆ ∼ 450 pc. Using the TNG50 simulation, we explored the evolution of the most massive progenitors of dwarf galaxies (z = 0 log(M⋆/M⊙) between 8.4 and 9.2) that end up as central galaxies of their groups. We split these dwarfs into three classes of the SMRz0: “Normals” from the central spine of the main branch, and “Compacts” from the secondary branch as well as the lower envelope of the main branch. Both classes of Compacts see their stellar sizes decrease from z ∼ 1 onwards in contrast to Normals, while the sizes of the gas and dark matter (DM) components continue to increase (as for Normals). A detailed analysis reveals that Compacts live in poorer environments, and thus suffer fewer major mergers from z = 0.8 onwards, which otherwise would pump angular momentum into the gas, allowing strong gas inflows, producing inner star formation, and thus leading to the buildup of a stellar core. Compacts are predicted to be rounder and to have bluer cores. Compact dwarfs of similar sizes are observed in the GAMA survey, but the bimodality in size is less evident and the most compact dwarfs tend to be passive rather than star forming, as in TNG50. Our conclusions should therefore be confirmed with future cosmological hydrodynamical simulations.