91̽��

Abstract:

During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, while glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further our understanding of galaxy formation. To tackle this fundamental issue, we investigate the impact of cosmic rays (CRs) and magnetic fields on the interstellar medium and the launching of outflows in a cosmological zoom-in simulation of a Milky Way-like galaxy. We find that including CRs decreases the stellar mass of the galaxy by a factor of 10 at high redshift and ∼4 at cosmic noon, leading to a stellar mass to halo mass ratio in good agreement with abundance matching models. Such decrease is caused by two effects: (i) a reduction of cold, high-density, star-forming gas, and (ii) a larger fraction of supernova (SN) events exploding at lower densities, where they have a higher impact. SN-injected CRs produce enhanced, multiphase galactic outflows, which are accelerated by CR pressure gradients in the circumgalactic medium of the galaxy. While the mass budget of these outflows is dominated by the warm ionized gas, warm neutral and cold gas phases contribute significantly at high redshifts. Importantly, our work shows that future JWST observations of galaxies and their multiphase outflows across cosmic time have the ability to constrain the role of CRs in regulating star formation.

Abstract:

Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction (f_esc) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ‘analogue’ galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14 120 mock spectra of star-forming galaxies with redshifts 4.64 ≤ z ≤ 10 from the SPHINX²⁰ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high f_esc. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must (i) track a high-specific star formation rate, (ii) be sensitive to stellar population age in the range 3.5–10 Myr representing the times when supernova first explode to when LyC production significantly drops, and (iii) include a proxy for neutral gas content and gas density in the interstellar medium. O₃₂, Σ_SFR, M_UV, and H β equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope (β), and E(B − V) match two or more of our criteria, rendering them good f_esc diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting f_esc based on direct observables. When applied to bright z > 6 Ly α emitters observed with JWST, we find that the majority of them have 𝑓esc≲10 per cent⁠.

Abstract:

Thick disks are a prevalent feature observed in numerous disk galaxies, including our own Milky Way. Their significance has been reported to vary widely, ranging from a few percent to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation, which has high spatial and stellar mass resolutions, to investigate the issue of the thick-disk mass fraction. We also use the NewHorizon2 simulation, which was run on the same initial conditions, but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M_* > 10¹⁰M_⊙ in NewHorizon, the contribution of the thick disk was found to be 20% ± 11% in r-band luminosity or 35% ± 15% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick-disk stars based on a single [α/Fe] cut is contaminated by stars of different kinematic properties, while missing the bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic midplane with high circularity and get heated with time via fluctuations in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties.

Abstract:

We study how better resolving the cooling length of galactic outflows affect their energetics. We perform radiativehydrodynamical galaxy formation simulations of an isolated dwarf galaxy (M = 108 M) with the RAMSES-RTZ code, accounting for non-equilibrium cooling and chemistry coupled to radiative transfer. Our simulations reach a spatial resolution of 18 pc in the interstellar medium (ISM) using a traditional quasi-Lagrangian scheme. We further implement a new adaptive mesh refinement strategy to resolve the local gas cooling length, allowing us to gradually increase the resolution in the stellar-feedback-powered outflows, from ≥ 200 pc to 18 pc. The propagation of outflows into the inner circumgalactic medium is significantly modified by this additional resolution, but the ISM, star formation, and feedback remain by and large the same. With increasing resolution in the diffuse gas, the hot outflowing phase (T > 8 × 104 K) systematically reaches overall higher temperatures and stays hotter for longer as it propagates outwards. This leads to two-fold increases in the time-averaged mass and metal outflow loading factors away from the galaxy (r = 5 kpc), a five-fold increase in the average energy loading factor, and a ≈50 per cent increase in the number of sightlines with NO VI ≥ 1013 cm−2. Such a significant boost to the energetics of outflows without new feedback mechanisms or channels strongly motivates future studies quantifying the efficiency with which better-resolved multiphase outflows regulate galactic star formation in a cosmological context.