Low-mass Halos Research Articles

The effect of cosmic web filaments on galaxy properties in the RESOLVE and ECO surveys

ABSTRACT Galaxy environment plays an important role in driving the transformation of galaxies from blue and star forming to red and quenched. Recent works have focused on the role of cosmic web filaments in galaxy evolution and have suggested that stellar mass segregation, quenching of star formation, and gas-stripping may occur within filaments. We study the relationship between distance to filament and the stellar mass, colour, and H i gas content of galaxies using data from the REsolved Spectroscopy of a Local VolumE survey and Environmental COntext (ECO) catalogue, two overlapping census-style, volume-complete surveys. We use the Discrete Persistence Structures Extractor to identify cosmic web filaments over the full ECO area. We find that galaxies close to filaments have higher stellar masses, in agreement with previous results. Controlling for stellar mass, we find that galaxies also have redder colours and are more gas poor closer to filaments. When accounting for group membership and halo mass, we find that these trends in colour and gas content are dominated by the increasing prevalence of galaxy group environments close to filaments, particularly for high-halo mass and low-stellar mass galaxies. Filaments have an additional small effect on the gas content of galaxies in low-mass haloes, possibly due to cosmic web stripping.

Probing wave-optics effects and low-mass dark matter halos with lensing of gravitational waves from massive black holes

Probing wave-optics effects and low-mass dark matter halos with lensing of gravitational waves from massive black holes

The mass accretion history of dark matter haloes down to Earth mass

ABSTRACT We take advantage of the unprecedented dynamical range provided by the ‘Cosmic-Zoom’ project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the Lambda cold dark matter paradigm when the dark matter is made of weakly interacting massive particles of mass 100 GeV. In particular, we complement previous studies by exploring the MAHs of haloes with mass from $10^8\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$ down to Earth mass, $10^{-6}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$. The formation redshift of low-mass haloes anticorrelates weakly with mass, peaking at z = 3 for haloes of mass $10^{-4}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$. Even lower masses are affected by the free-streaming cut-off in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press–Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological N-body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, $M\gt 10^{10}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$, the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than 20 per cent at high redshift. If this result can be extrapolated to lower halo masses, we conclude that EPS theory predicts the hierarchical build up of dark matter haloes quite well over the entire halo mass range.

A Self-interacting Dark Matter Solution to the Extreme Diversity of Low-mass Halo Properties

The properties of low-mass dark matter (DM) halos appear to be remarkably diverse relative to cold, collisionless DM predictions, even in the presence of baryons. We show that self-interacting DM (SIDM) can simultaneously explain observations of halo diversity at two opposite extremes—the inner density profile of the dense substructure perturbing the strong lens galaxy SDSSJ0946+1006 and the rotation curves of isolated, gas-rich ultradiffuse galaxies (UDGs). To achieve this, we present the first cosmological zoom-in simulation featuring strong DM self-interactions in a galaxy group environment centered on a 1013 M ⊙ host halo. In our SIDM simulation, most surviving subhalos of the group-mass host are deeply core-collapsed, yielding excellent candidates for the observed dense strong-lens perturber. Self-interactions simultaneously create kiloparsec-scale cores in low-concentration isolated halos, which could host the observed UDGs. Our scenario can be further tested with observations of DM structure and galaxies over a wide mass range.

Bridging the Gap between Cosmic Dawn and Reionization Favors Models Dominated by Faint Galaxies

It has been claimed that traditional models struggle to explain the tentative detection of the 21 cm absorption trough centered at z ∼ 17 measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power law of deep Hubble Space Telescope observations of 4 < z < 9 galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction at z ∼ 6–8, Thomson-scattering optical depth, and ionizing emissivity at z ∼ 5. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile existing high-redshift (cosmic dawn) and low-redshift (reionization) constraints. In particular, we find that high contributions from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with a gradual and longer duration, resulting in a higher optical depth. While the models dominated by faint galaxies successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically motivated semianalytical prescription produces results that are consistent with the (sophisticated) state-of-the-art THESAN radiation-magnetohydrodynamic simulation of the reionization.

Ghostly Galaxies: Accretion-dominated Stellar Systems in Low-mass Dark Matter Halos

Wide-area deep imaging surveys have discovered large numbers of extremely low surface brightness (LSB) dwarf galaxies, which challenge galaxy formation theory and, potentially, offer new constraints on the nature of dark matter. Here we discuss one as-yet-unexplored formation mechanism that may account for a fraction of LSB dwarfs. We call this the “ghost galaxy” scenario. In this scenario, inefficient radiative cooling prevents star formation in the “main branch” of the merger tree of a low-mass dark matter halo, such that almost all its stellar mass is acquired through mergers with less massive (but nevertheless star-forming) progenitors. Present-day systems formed in this way would be “ghostly” isolated stellar halos with no central galaxy. We use merger trees based on the extended Press–Schechter formalism and the Copernicus Complexio cosmological N-body simulation to demonstrate that mass assembly histories of this kind can occur for low-mass halos in ΛCDM, but they are rare. They are most probable in isolated halos of present-day mass ∼4 × 109 M ⊙, occurring for ∼5% of all halos of that mass under standard assumptions about the timing and effect of cosmic reionization. The stellar masses of star-forming progenitors in these systems are highly uncertain; abundance-matching arguments imply a bimodal present-day mass function having a brighter population (median M ⋆ ∼ 3 × 106 M ⊙) consistent with the tail of the observed luminosity function of ultradiffuse galaxies. This suggests that observable analogs of these systems may await discovery. We find that a stronger ionizing background (globally or locally) produces brighter and more extended ghost galaxies.

The first fireworks: A roadmap to Population III stars during the epoch of reionization through pair-instability supernovae

ABSTRACT With the launch of JWST and other scheduled missions aimed at probing the distant universe, we are entering a new promising era for high-z astronomy. One of our main goals is the detection of the first population of stars (Population III or Pop III stars), and models suggest that Pop III star formation is allowed well into the Epoch of Reionization (EoR), rendering this an attainable achievement. In this paper, we focus on our chance of detecting massive Pop IIIs at the moment of their death as Pair-Instability Supernovae (PISNe). We estimate the probability of discovering PISNe during the EoR in galaxies with different stellar masses (7.5 ≤ Log(M⋆/M⊙) ≤ 10.5) from six dustyGadget simulations of 50 h−1 cMpc per side. We further assess the expected number of PISNe in surveys with JWST/NIRCam and Roman/WFI. On average, less than one PISN is expected in all examined JWST fields at z ≃ 8 with Δz = 1, and O(1) PISN may be found in a ∼1 deg2 Roman field in the best-case scenario, although different assumptions on the Pop III IMF and/or Pop III star formation efficiency can decrease this number substantially. Including the contribution from unresolved low-mass haloes holds the potential for increased discoveries. JWST/NIRCam and Roman/WFI allow the detection of massive-progenitor (∼250 M⊙) PISNe throughout all the optimal F200W-F356W, F277W-F444W, and F158-F213 colours. PISNe are also pre-dominantly located at the outskirts of their hosting haloes, facilitating the disentangling of underlying stellar emission thanks to the spatial-resolution capabilities of the instruments.

Tracking the evolution of satellite galaxies: mass stripping and dark-matter deficient galaxies

ABSTRACT Satellite galaxies undergo a variety of physical processes when they are accreted by groups and clusters, often resulting in the loss of baryonic and dark matter (DM) mass. In this work, we evaluate the predictions from the IllustrisTNG hydrodynamical simulation regarding the evolution of the matter content of satellites, focusing on a population that are accreted at z &amp;gt; 1 and retain their identity as satellites down to z = 0. At fixed host halo mass, the amount of DM and stellar mass stripped depends mostly on the pericentric distance, dperi, here normalized by host halo virial radius. The closest encounters result in significant loss of DM, with subhaloes retaining between 20 and a few per cent of their z = 1 mass. At fixed dperi, DM mass stripping seems more severe in lower mass haloes. Conversely, the average satellite in higher mass haloes has its stellar mass growth halted earlier, having lost a higher fraction of stellar mass by z = 0. We also show that mass stripping has a strong impact on the quenched fractions. The IllustrisTNG boxes are qualitatively consistent in these predictions, with quantitative differences mostly originating from the distinct subhalo mass ranges covered by the boxes. Finally, we have identified DM-deficient systems in all TNG boxes. These objects are preferentially found in massive clusters (Mhost ≳ 1013 M⊙), had very close encounters with their central galaxies ($d_{\rm peri}\simeq 0.05\, R_{\rm vir}$), and were accreted at high redshift (zinfall ≳ 1.4), reinforcing the notion that tidal stripping is responsible for their remarkable lack of DM.

The X-ray invisible Universe. A look into the haloes undetected by eROSITA

ABSTRACT The paper presents the analysis of optically selected GAMA groups and clusters in the SRG/eROSITA X-ray map of eROSITA Final Equatorial Depth Survey, in the halo mass range 1013−5 × 1014 M⊙ and at z &amp;lt; 0.2. All X-ray detections have a clear GAMA counterpart, but most of the GAMA groups in the halo mass range 1013−1014 M⊙ remain undetected. We compare the X-ray surface brightness profiles of the eROSITA detected groups with the mean stacked profile of the undetected low-mass haloes at fixed halo mass. Overall, we find that the undetected groups exhibit less concentrated X-ray surface brightness, dark matter, and galaxy distributions with respect to the X-ray-detected haloes. The mean gas mass fraction profiles are consistent in the two samples within 1.5σ, indicating that the gas follows the dark matter profile. The low-mass concentration and the magnitude gap indicate that these systems are young. They reside with a higher probability in filaments while X-ray-detected groups favour the nodes of the Cosmic Web. Because of the lower central emission, the undetected systems tend to be X-ray underluminous at fixed halo mass and to lie below the LX−Mhalo relation. Interestingly, the X-ray-detected systems inhabiting the nodes scatter the less around the relation, while those in filaments tend to lie below it. We do not observe any strong relationship between the system X-ray appearance and the active galactic nucleus (AGN) activity. We cannot exclude the role of the past AGN feedback in affecting the gas distribution over the halo lifetime. However, the data suggests that the observed differences might be related to the halo assembly bias.

Signatures of cosmic ray heating in 21-cm observables

ABSTRACT Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D seminumerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass haloes, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman–Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal.

Thesan-hr: how does reionization impact early galaxy evolution?

ABSTRACT The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir &amp;lt; 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 = − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

Halo scaling relations and hydrostatic mass bias in the simba simulation from realistic mock X-ray catalogues

ABSTRACT We present a new end-to-end pipeline for Mock Observations of X-ray Haloes and Analysis (moxha) for hydrodynamic simulations of massive haloes, and use it to investigate X-ray scaling relations and hydrostatic mass bias in the simba cosmological hydrodynamic simulation for haloes with M500 ∼ 1013−15M⊙. moxha ties together existing yT-based software packages and adds new functionality to provide an end-to-end pipeline for generating mock X-ray halo data from large-scale or zoom simulation boxes. We compare moxha-derived halo properties in simba to their emission-weighted counterparts, and forecast the systematic mass bias in mock Athena observations. Overall, we find inferred hydrostatic masses are biased low compared to true simba values. For simple mass-weighting, we find $b_\text{MW} = 0.15^{+0.15}_{-0.14}$ (16–84 per cent range), while emission-weighting increases this to $b_\text{LW}=0.30^{+0.19}_{-0.10}$. The larger bias versus mass-weighted values we attribute to the spectroscopic and emission-weighted temperatures being biased systematically lower than mass-weighted temperatures. The full moxha pipeline recovers the emission-weighted hydrostatic masses at R500 reasonably well, yielding $b_\text{X}=0.33^{+0.28}_{-0.34}$. moxha-derived halo X-ray scalings are in very good agreement with observed scaling relations, with the inclusion of lower mass groups significantly steepening the LX − M500, M500 − TX, and LX − TX relations. This indicates the strong effect the simba feedback model has on low-mass haloes, which strongly evacuates poor groups but still retains enough gas to reproduce observations. We find similar trends for analogous scaling relations measured at R500, as expected for halo-wide gas evacuation.

The impact of UV variability on the abundance of bright galaxies at z ≥ 9

ABSTRACT JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the Λ cold dark matter (ΛCDM) cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function of galaxies at z ≥ 9. Our study differentiates between the median UV radiation yield and the variability of UV luminosities of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, σUV ≈ 0.75 mag, matches the observational constraints at z ≈ 9. However, increasingly larger σUV is necessary when moving to higher redshifts, reaching $\sigma _{\rm UV} \approx 2.0\, (2.5)\, {\rm mag}$ at z ≈ 12 (16). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and dust clearance during strong feedback cycles. Photometric constraints from JWST at z ≳ 9 therefore can be reconciled with a standard ΛCDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield.

Pushing the limits of detectability: mixed dark matter from strong gravitational lenses.

One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter haloes. There is a large volume of information in strong gravitational lens images - the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter haloes relative to warm dark matter. Since the strong lensing signal receives a contribution from a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter haloes, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with 40 lenses with Bayesian odds of 30:1.

Effect of Neutrinos on Angular Momentum of Dark Matter Halo

Massive neutrinos are expected to affect the large-scale structure formation, including the major component of solid substances, dark matter halos. How halos are influenced by neutrinos is vital and interesting, and angular momentum (AM) as a significant feature provides a statistical perspective for this issue. Exploring halos from TianNu N-body cosmological simulation with the co-evolving neutrino particles, we obtain some concrete conclusions. First, by comparing the same halos with and without neutrinos, in contrast to the neutrino-free case, over 89.71% of halos have smaller halo moduli, over 71.06% have smaller particle-mass-reduced (PMR) AM moduli, and over 95.44% change their orientations of less than 0.°65. Moreover, the relative variation of PMR modulus is more visible for low-mass halos. Second, to explore the PMR moduli of halos in dense or sparse areas, we divide the whole box into big cubes, and search for halos within a small spherical cell in a single cube. From the two-level divisions, we discover that in denser cubes, the variation of PMR moduli with massive neutrinos decreases more significantly. This distinction suggests that neutrinos exert heavier influence on halos’ moduli in compact regions. With massive neutrinos, most halos (86.60%) have lower masses than without neutrinos.

The MillenniumTNG Project: an improved two-halo model for the galaxy–halo connection of red and blue galaxies

ABSTRACT Approximate methods to populate dark-matter haloes with galaxies are of great utility to galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, z = 1 and z = 0, in the large-volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to r ∼ 0.1 h−1 Mpc, while here we focus on modelling galaxy assembly bias in the two-halo regime, r ≳ 1 h−1 Mpc. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales (r ∼ 20 h−1 Mpc), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10–15 per cent discrepancy is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties (‘shear’ being the best-performing one) rather than intrinsic ones (e.g. concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and ‘shear’) is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher order statistics such as nearest neighbour counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters.

JADES Imaging of GN-z11: Revealing the Morphology and Environment of a Luminous Galaxy 430 Myr after the Big Bang

We present JWST NIRCam nine-band near-infrared imaging of the luminous z = 10.6 galaxy GN-z11 from the JWST Advanced Deep Extragalactic Survey of the GOODS-N field. We find a spectral energy distribution (SED) entirely consistent with the expected form of a high-redshift galaxy: a clear blue continuum from 1.5 to 4 μm with a complete dropout in F115W. The core of GN-z11 is extremely compact in JWST imaging. We analyze the image with a two-component model, using a point source and a Sérsic profile that fits to a half-light radius of 200 pc and an index n = 0.9. We find a low-surface-brightness haze about 0.″4 to the northeast of the galaxy, which is most likely a foreground object but might be a more extended component of GN-z11. At a spectroscopic redshift of 10.60 (Bunker et al. ), the comparison of the NIRCam F410M and F444W images spans the Balmer jump. From population-synthesis modeling, here assuming no light from an active galactic nucleus, we reproduce the SED of GN-z11, finding a stellar mass of ∼109 M ⊙, a star formation rate of ∼20 M ⊙ yr−1, and a young stellar age of ∼20 Myr. Since massive galaxies at high redshift are likely to be highly clustered, we search for faint neighbors of GN-z11, finding nine galaxies out to ∼5 comoving Mpc transverse with photometric redshifts consistent with z = 10.6, and a tenth more tentative dropout only 3″ away. This is consistent with GN-z11 being hosted by a massive dark-matter halo (≈8 × 1010 M ⊙), though lower halo masses cannot be ruled out.

Environmental Dependence of the Mass–Metallicity Relation in Cosmological Hydrodynamical Simulations

We investigate the environmental dependence of the gas-phase metallicity for galaxies at z = 0 to z ≳ 2 and the underlying physical mechanisms driving this dependence using state-of-the-art cosmological hydrodynamical simulations. We find that, at fixed stellar mass, central galaxies in massive halos have lower gas-phase metallicity than those in low-mass halos, while satellite galaxies residing in more massive halos are more metal-rich. The combined effect is that massive galaxies are more metal-poor in massive halos and low-mass galaxies are more metal-rich in massive halos. By inspecting the environmental dependence of other galaxy properties, we identify that the accretion of low-metallicity gas is responsible for the environmental dependence of central galaxies at high z, whereas the AGN feedback processes play a crucial role at low z. For satellite galaxies, we find that both the suppression of gas accretion and the stripping of existing gas are responsible for their environmental dependence, with negligible effect from the AGN feedback. Finally, we show that the difference of gas-phase metallicity as a function of stellar mass between protocluster and field galaxies agrees with recent observational results, for example from the MAMMOTH-Grism survey.

Microlensing of gravitational waves by dark matter structures

Gravitational lensing of gravitational waves provides a potential new probe of dark matter structures. In this work, we consider the microlensing effect on gravitational wave signals from black hole binaries induced by low-mass dark matter halos that do not retain enough baryonic matter to hold stars. We clarify systematically when this microlensing effect is relevant and study in detail its detectability by future gravitational wave observatories. We consider lensing by cold dark matter halos and by solitonic cores that reside in fuzzy dark matter halos. Our results show that although the effect can be detectable at relatively large impact parameters, the probability of detecting such lensed events is low. In particular, we find that the expected number of events lensed by cold dark matter halos is 𝒪(0.01) per year for BBO and the expected number of events lensed by solitonic cores inside fuzzy dark matter halos is 𝒪(0.01) per year for ET. In the case that a significant fraction of dark matter consists of 𝒪(100 M ⊙) objects that are relatively compact, R < 𝒪(0.1 pc), we show that the expected number of lensed events per year for ET can be very large, 𝒪(1000).

Starbursts in low-mass haloes at Cosmic Dawn. I. The critical halo mass for star formation

ABSTRACT The first stars, galaxies, star clusters, and direct-collapse black holes are expected to have formed in low-mass (∼105–109 M⊙) haloes at Cosmic Dawn (z ∼ 10–30) under conditions of efficient gas cooling, leading to gas collapse towards the centre of the halo. The halo mass cooling threshold has been analysed by several authors using both analytical models and numerical simulations, with differing results. Since the halo number density is a sensitive function of the halo mass, an accurate model of the cooling threshold is needed for (semi-)analytical models of star formation at Cosmic Dawn. In this paper, the cooling threshold mass is calculated (semi-)analytically, considering the effects of H2-cooling and formation (in the gas phase and on dust grains), cooling by atomic metals, Lyman-α cooling, photodissociation of H2 by Lyman–Werner photons (including self-shielding by H2), photodetachment of H− by infrared photons, photoevaporation by ionization fronts, and the effect of baryon streaming velocities. We compare the calculations to several high-resolution cosmological simulations, showing excellent agreement. We find that in regions of typical baryon streaming velocities, star formation is possible in haloes of mass ≳ 1–2 × 106 M⊙ for z ≳ 20. By z ∼ 8, the expected Lyman–Werner background suppresses star formation in all minihaloes below the atomic cooling threshold (Tvir = 104 K). The halo mass cooling threshold increases by another factor of ∼4 following reionization, although this effect is slightly delayed (z ∼ 4–5) because of effective self-shielding.