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 »

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
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. 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 »

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 »

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 »

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 »