51�������

Abstract:

We use the cosmological hydrodynamical simulation NewHorizon to study the effects of the baryonic component on the inner mass profile of dark matter haloes of isolated galaxies (). Dark matter deficits (‘cores’) develop only in galaxies in a narrow range of stellar mass, . The lower stellar mass limit arises because a minimum amount of star formation is required to drive the baryonic outflows that redistribute dark matter and create a core. The upper limit roughly coincides with the total amount of dark matter initially contained within the innermost 2 kpc (), which roughly coincides with the stellar half-mass radius of these dwarfs. This enclosed mass is quite insensitive to the total virial mass of the system. The same upper limit applies to other simulations, like NIHAO and EAGLE-CHT10, despite their rather different galaxy formation efficiencies. This suggests that it is the galaxy total stellar mass that determines when a core is formed, and not the galaxy-to-dark halo mass ratio, as argued in earlier work. This is consistent with a back-of-the-envelope estimate for a SN-induced rate of orbital diffusion. Although NewHorizon dwarfs reproduce the observed diversity of rotation curves better than other simulations, there are significant differences in the gravitational importance of baryons in the inner regions of dwarfs compared to observations. These differences prevent us from concluding that cosmological simulations are currently fully able to account for the observed diversity of rotation curve shapes.

Abstract:

We analyze JWST NIRSpec+MIRI/MRS observations of the infrared (IR) polycyclic aromatic hydrocarbon (PAH) features in the central regions (∼0.26″ at 6 μm; ∼50–440 pc depending on the source) of local luminous IR galaxies. In this work, we examine the effect of nuclear obscuration on the PAH features of deeply obscured nuclei, predominantly found in local luminous IR galaxies, and we compare these nuclei with “normal” star-forming regions. We extend previous work to include shorter wavelength PAH ratios now available with the NIRSpec+MIRI/MRS spectral range. We introduce a new diagnostic diagram for selecting deeply obscured nuclei based on the 3.3 and 6.2 μm PAH features and/or mid-IR continuum ratios at ∼3 and 5 μm. We find that the PAH equivalent width ratio of the brightest PAH features at shorter wavelengths (at 3.3 and 6.2 μm) is impacted by nuclear obscuration. Although the sample of luminous IR galaxies used in this analysis is relatively small, we find that sources exhibiting a high silicate absorption feature cluster tightly in a specific region of the diagram, whereas star-forming regions experiencing lower extinction levels occupy a different area in the diagram. This demonstrates the potential of this technique to identify buried nuclei. To leverage the excellent sensitivity of the MIRI imager on board JWST, we extend our method of identifying deeply obscured nuclei at higher redshifts using a selection of MIRI filters. Specifically, the combination of various MIRI JWST filters enables the identification of buried sources beyond the local Universe and up to z ∼ 3, where other commonly used obscuration tracers such as the 9.7 μm silicate band, are out of the spectral range of MRS. Our results pave the way for identifying distant deeply obscured nuclei with JWST.

Abstract:

JWST observations are providing unprecedented constraints on the history of reionization owing to the ability to detect faint galaxies at z ≫ 6. Modelling this history requires understanding both the ionizing photon production rate (ξion) and the fraction of those photons that escape into the intergalactic medium (fesc). Observational estimates of these quantities generally rely on spectroscopy for which large samples with well-defined selection functions remain limited. To overcome this challenge, we present and release a novel implicit likelihood inference pipeline, PHOTONIOn, trained on mock photometry to predict the escaped ionizing luminosity of individual galaxies (N ion) based on photometric magnitudes and redshifts. We show that PHOTONIOn is able to reliably infer N ion from photometry. This is in contrast to traditional spectral energy distribution-fitting approaches which rely on fesc prescriptions that often overpredict N ion for Lyman Continuum (LyC)-dim galaxies, even when given access to spectroscopic data. We have deployed PHOTONIOn on a sample of 4559 high-redshift galaxies from the JWST Advanced Deep Extragalactic Survey (JADES), finding gentle redshift evolutions of log10(N ion) = (0.08 ± 0.01)z + (51.60 ± 0.06) and log10(fescξion) = (0.07 ± 0.01)z + (24.12 ± 0.07). Late-time values for the ionizing photon production rate density are consistent with both theoretical models and observations. Finally, we measure the evolution of the intergalactic medium ionized fraction to find that observed populations of star-forming galaxies are capable of driving reionization in this field to completion by z ∼ 5.3 without the need for active galactic nucleus or other exotic sources, consistent with other studies of the same field. The 20 per cent of UV-brightest galaxies (MUV < −18.5) reionize roughly 35 per cent of the survey volume, demonstrating that UV faint LyC emitters are crucial for reionization.

Abstract:

<jats:p>Star formation is a key process that governs the baryon cycle within galaxies, however, the question of how it controls their growth remains elusive due to modeling uncertainties. To understand the impact of star formation models on galaxy evolution, we performed cosmological zoom-in radiation-hydrodynamic simulations of a dwarf dark matter halo, with a virial mass of <jats:italic>M</jats:italic><jats:sub>vir</jats:sub> ∼ 10<jats:sup>9</jats:sup> <jats:italic>M</jats:italic><jats:sub>⊙</jats:sub> at <jats:italic>z</jats:italic> = 6. We compared two different star formation models: a multi-freefall model combined with a local gravo-thermo-turbulent condition and a more self-consistent model based on a sink particle algorithm, where gas accretion and star formation are directly controlled by the gas kinematics. As the first study in this series, we used cosmological zoom-in simulations with different spatial resolutions and found that star formation is more bursty in the runs with the sink algorithm, generating stronger outflows than in the runs with the gravo-thermo-turbulent model. The main reason for the increased burstiness is that the gas accretion rates on the sinks are high enough to form stars on very short timescales, leading to more clustered star formation. As a result, the star-forming clumps are disrupted more quickly in the sink run due to more coherent radiation and supernova feedback. The difference in burstiness between the two star formation models becomes even more pronounced when the supernova explosion energy is artificially increased. Our results suggest that improving the modeling of star formation on small, sub-molecular cloud scales can significantly impact the global properties of simulated galaxies.</jats:p>