On unveiling buried nuclei with JWST: A technique for hunting the most obscured galaxy nuclei from local to high redshift
Abstract:
We analyze JWST NIRSpec+MIRI/MRS observations of the infrared (IR) polycyclic aromatic hydrocarbon (PAH) features in the central regions (a at 6 μm; a 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 astar-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 a3 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 za 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.Inferring the ionizing photon contributions of high-redshift galaxies to reionization with JWST NIRCam photometry
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.Impact of star formation models on the growth of simulated galaxies at high redshifts
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>Evaluating the variance of individual halo properties in constrained cosmological simulations
Abstract:
Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation. We introduce an internal method for quantifying the extent to which the variance of individual halo properties is suppressed by the constraints imposed on the initial conditions. We apply it to the Constrained Simulations in BORG (CSiBORG) suite of 101 high-resolution realizations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, measuring the degree to which the haloes' particles are initially coincident. This addresses the extent to which the imposed large-scale structure constraints reduce the variance of individual halo properties. We find consistent reconstructions of M≳1014M⊙h-1 haloes, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions, and peculiar velocities of clusters to high precision, though we do not assess how well they reproduce observations of the local Universe. The effect of the constraints tapers off towards lower mass, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions and describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG.Dwarf galaxies as a probe of a primordially magnetized Universe
Abstract:
Aims: The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By comparing both we can determine the dominant agent in shaping galaxy evolution.
Methods: We use the magneto-hydrodynamics code RAMSES, to generate multiple new zoom-in simulations for eight different host halos of dwarf galaxies across a wide luminosity range of 103 − 106 L⊙. These halos are selected from a ΛCDM cosmological box, tracking their evolution down to redshift z = 0. We explore a variety of primordial magnetic field (comoving) strengths Bλ ranging from 0.05 to 0.50 nG.
Results: We find magnetic fields in the interstellar medium not only modify star formation in dwarf spheroidal galaxies but also completely prevent the formation of stars in less compact ultra-faints with halo mass and stellar mass below ∼ 2.5 · 109 and 3 · 106 M⊙, respectively. At high redshifts, the impact of PMFs on host halos of dwarf galaxies through the modification of the matter power spectrum is more dominant than the influence of magneto-hydrodynamics in shaping their gaseous structure. Through the amplification of small perturbations ranging in mass from 107 to 109 M⊙ in the ΛCDM+PMFs matter power spectrum, primordial fields expedite the formation of the first dark matter halos, leading to an earlier onset and a higher star formation rate at redshifts z > 12. We investigate the evolution of various energy components and demonstrate that magnetic fields with an initial strength of Bλ ≥ 0.05 nG exhibit a strong growth of magnetic energy, accompanied by a saturation phase, that starts quickly after the growth phase. These trends persist consistently, regardless of the initial conditions, whether it is the classical ΛCDM or modified by PMFs. Lastly, we investigate the impact of PMFs on the present-time observable properties of dwarf galaxies, namely, the half light radius, V-band luminosity, mean metallicity and velocity dispersion profile. We find that PMFs with moderate strengths of Bλ ≤ 0.10 nG show great agreement with the scaling relations of the observed Local group dwarfs. However, stronger fields lead to large sizes and high velocity dispersion.