skip to main content

Title: Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars
ABSTRACT

We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z = 6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified ‘delayed cooling’ feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n ∼ 103–105 cm−3). The sites where star formation ignites are characterized by maximum densities of the order of 105 cm−3 and gravitational pressures Pgrav/k >107 K cm−3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z = 6.14 is 3.4$\times 10^7~\rm M_{\odot }$, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of $\sim 10^6~\rm M_{\odot }$ and half-mass sizes of ∼100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one orders of magnitude lower than the observed one. In the size–mass plane, our more » clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

« less
Authors:
; ; ; ; ; ; ; ; ;
Publication Date:
NSF-PAR ID:
10373183
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
4
Page Range or eLocation-ID:
p. 5914-5934
ISSN:
0035-8711
Publisher:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factorsmore »∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.« less
  2. 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).more »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).

    « less
  3. ABSTRACT

    We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including themore »stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.

    « less
  4. ABSTRACT

    Current observations favour that the massive ultraviolet-bright clumps with a median stellar mass of $\sim 10^7\, {\rm M}_{\odot }$, ubiquitously observed in z ∼ 1–3 galaxies, are star-forming regions formed in situ in galaxies. It has been proposed that they result from gas fragmentation due to gravitational instability of gas-rich, turbulent, and high-redshift discs. We bring support to this scenario by reporting the new discovery of giant molecular clouds (GMCs) in the strongly lensed, clumpy, main-sequence galaxy, A521-sys1, at z = 1.043. Its CO(4–3) emission was mapped with the Atacama Large Millimetre/submillimetre Array (ALMA) at an angular resolution of 0.19 × 0.16 arcsec2, reading down to 30 pc, thanks to gravitational lensing. We identified 14 GMCs, most being virialized, with $10^{5.9}-10^{7.9}\, {\rm M}_{\odot }$ masses and a median $800\, {\rm M}_{\odot }~\mathrm{pc}^{-2}$ molecular gas mass surface density, that are, respectively, 100 and 10 times higher than for nearby GMCs. They are also characterized by 10 times higher supersonic turbulence with a median Mach number of 60. They end up to fall above the Larson scaling relations, similarly to the GMCs in another clumpy z ≃ 1 galaxy, the Cosmic Snake, although differences between the two sets of high-redshift GMCs exist. Altogether theymore »support that GMCs form with properties that adjust to the ambient interstellar medium conditions prevalent in the host galaxy whatever its redshift. The detected A521-sys1 GMCs are massive enough to be the parent gas clouds of stellar clumps, with a relatively high star formation efficiency per free-fall time of ∼11 per cent.

    « less
  5. ABSTRACT The galaxy size–stellar mass and central surface density–stellar mass relationships are fundamental observational constraints on galaxy formation models. However, inferring the physical size of a galaxy from observed stellar emission is non-trivial due to various observational effects, such as the mass-to-light ratio variations that can be caused by non-uniform stellar ages, metallicities, and dust attenuation. Consequently, forward-modelling light-based sizes from simulations is desirable. In this work, we use the skirt  dust radiative transfer code to generate synthetic observations of massive galaxies ($M_{*}\sim 10^{11}\, \rm {M_{\odot }}$ at z = 2, hosted by haloes of mass $M_{\rm {halo}}\sim 10^{12.5}\, \rm {M_{\odot }}$) from high-resolution cosmological zoom-in simulations that form part of the Feedback In Realistic Environments project. The simulations used in this paper include explicit stellar feedback but no active galactic nucleus (AGN) feedback. From each mock observation, we infer the effective radius (Re), as well as the stellar mass surface density within this radius and within $1\, \rm {kpc}$ (Σe and Σ1, respectively). We first investigate how well the intrinsic half-mass radius and stellar mass surface density can be inferred from observables. The majority of predicted sizes and surface densities are within a factor of 2 of the intrinsic values.more »We then compare our predictions to the observed size–mass relationship and the Σ1−M⋆ and Σe−M⋆ relationships. At z ≳ 2, the simulated massive galaxies are in general agreement with observational scaling relations. At z ≲ 2, they evolve to become too compact but still star forming, in the stellar mass and redshift regime where many of them should be quenched. Our results suggest that some additional source of feedback, such as AGN-driven outflows, is necessary in order to decrease the central densities of the simulated massive galaxies to bring them into agreement with observations at z ≲ 2.« less