Julien Devriendt

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.

Abstract:

The scaling relations between the gas content and star formation rate of galaxies provide useful insights into the processes governing their formation and evolution. We investigated the emergence and the physical drivers of the global Kennicutt-Schmidt (KS) relation at 0:25 ≤ z ≤ 4 in the cosmological hydrodynamic simulation NewHorizon, capturing the evolution of a few hundred galaxies with a resolution down to 34 pc. The details of this relation vary strongly with the stellar mass of galaxies and the redshift. A power-law relation ΣSFR / Σa gas with a ≈ 1:4, like that found empirically, emerges at z ≈ 2..3 for the more massive half of the galaxy population. However, no such convergence is found in the lower-mass galaxies, for which the relation gets shallower with decreasing redshift. At galactic scales, the star formation activity correlates with the level of turbulence of the interstellar medium, quantified by the Mach number, rather than with the gas fraction (neutral or molecular), confirming the conclusions found in previous works. With decreasing redshift, the number of outliers with short depletion times diminishes, reducing the scatter of the KS relation, while the overall population of galaxies shifts toward low densities. Our results, from parsec-scale star formation models calibrated with local Universe physics, demonstrate that the cosmological evolution of the environmental (e.g., mergers) and internal conditions (e.g., gas fractions) conspire to shape the KS relation. This is an illustration of how the interplay of global and local processes leaves a detectable imprint on galactic-scale observables and scaling relations.