KICC papers

arXiv:2409.11469v2 Announce Type: replace Abstract: We study the implications of relaxing the requirement for ultralight axions to account for all dark matter in the Universe by examining mixed dark matter (MDM) cosmologies with axion fractions $f \leq 0.3$ within the fuzzy dark matter (FDM) window $10^{-25}$ eV $\lesssim m \lesssim 10^{-23}$ eV. Our simulations, using a new MDM gravity solver implemented in AxiREPO, capture wave dynamics across various scales with high accuracy down to redshifts $z\approx 1$. We identify halos with Rockstar using the CDM component and find good agreement of inferred halo mass functions (HMFs) and concentration-mass relations with theoretical models across redshifts $z=1-10$. This justifies our halo finder approach a posteriori as well as the assumptions underlying the MDM halo model AxionHMcode. Using the inferred axion halo mass-cold halo mass relation $M_{\text{a}}(M_{\text{c}})$ and calibrating a generalised smoothing parameter $\alpha$ to our MDM simulations, we present a new version of AxionHMcode. The code exhibits excellent agreement with simulations on scales $k< 20 \ h$ cMpc$^{-1}$ at redshifts $z=1-3.5$ for $f\leq 0.1$ around the fiducial axion mass $m = 10^{-24.5}$ eV $ = 3.16\times 10^{-25}$ eV, with maximum deviations remaining below 10%. For axion fractions $f\leq 0.3$, the model maintains accuracy with deviations under 20% at redshifts $z\approx 1$ and scales $k< 10 \ h$ cMpc$^{-1}$, though deviations can reach up to 30% for higher redshifts when $f=0.3$. Reducing the run-time for a single evaluation of AxionHMcode to below $1$ minute, these results highlight the potential of AxionHMcode to provide a robust framework for parameter sampling across MDM cosmologies in Bayesian constraint and forecast analyses.

arXiv:2310.03067v3 Announce Type: replace Abstract: Merger events can trigger gas accretion onto supermassive black holes (SMBHs) located at the centre of galaxies, and form close pairs of AGN. The fraction of AGN in pairs offers critical insights into the dynamics of galaxy interactions, SMBH growth, and their co-evolution with host galaxies. However, the identification of dual AGN is difficult, as it requires high-quality spatial and spectral data; hence, only very few pairs have been found in the distant Universe so far. This study aims to provide a first observational estimate of the fraction of dual AGN at 2

arXiv:2503.04893v1 Announce Type: new Abstract: We present measurements of the spatially averaged HeII photo-ionization rate ($\langle \Gamma_{\rm HeII} \rangle$), mean free path of HeII ionizing photons ($\lambda_{\rm mfp, HeII}$), and HeII fraction ($f_{\rm HeII}$) across seven redshift bins within the redshift range $2

arXiv:2503.04894v1 Announce Type: new Abstract: We explore the evolution of galaxy sizes at high redshift ($36$ show rising star-formation histories, and therefore are in a compaction phase. We demonstrate that H$\alpha$ emission is systematically extended relative to the UV continuum by a median factor of 1.7, consistent with recent JWST studies. However, in contrast to previous interpretations that link extended H$\alpha$ sizes to inside-out growth, we find that Lyman-continuum (LyC) emission is spatially disconnected from H$\alpha$. Instead, a simple Str\"{o}mgren sphere argument reproduces observed trends, suggesting that extreme LyC production during central starbursts is the primary driver of extended nebular emission.

arXiv:2503.04243v1 Announce Type: new Abstract: Understanding the physical mechanisms that drive star formation is crucial for advancing our knowledge of galaxy evolution. We explore the interrelationships between key galaxy properties associated with star formation, with a particular focus on the impact of dark matter halos. Given the sensitivity of atomic hydrogen (HI) to external processes, we concentrate exclusively on central spiral galaxies. We find that the molecular-to-atomic gas mass ratio ($M_{\rm H_2}/M_{\rm HI}$) strongly depends on stellar mass and specific star formation rate (sSFR). In the star formation efficiency (SFE)-sSFR plane, most galaxies fall below the H$_2$ fundamental formation relation (FFR), with SFE$_{\rm HI}$ being consistently lower than SFE$_{\rm H_2}$. Using the improved halo masses derived by Zhao et al. (2025), for star-forming galaxies, both SFE$_{\rm HI}$ and $M_{\rm H_2}/M_{\rm HI}$ increase rapidly and monotonically with halo mass, indicating a higher efficiency in converting HI to H$_2$ in more massive halos. This trend ultimately leads to the unsustainable state where SFE$_{\rm HI}$ exceeds SFE$_{\rm H_2}$ at halo mass around $10^{12} \hbox{$M_{\odot}$}$. For halos with masses exceeding $10^{12} \hbox{$M_{\odot}$}$, galaxies predominantly experience quenching. We propose a plausible evolutionary scenario in which the growth of halo mass regulates the conversion of HI to H$_2$, star formation, and the eventual quenching of galaxies. The disk size, primarily regulated by the mass, spin and concentration of the dark matter halo, also significantly influences HI to H$_2$ conversion and star formation. These findings underscore the critical role of dark matter halos as a global regulator of galaxy-wide star formation, a key factor that has been largely underappreciated in previous studies.

arXiv:2503.03806v1 Announce Type: new Abstract: Population III (Pop III) stars are the first stars in the Universe, forming from pristine, metal-free gas and marking the end of the cosmic dark ages. Their formation rate is expected to sharply decline after redshift $z \approx 15$ due to metal enrichment from previous generations of stars. In this paper, we analyze 14 zoom-in simulations from the THESAN-ZOOM project, which evolves different haloes from the THESAN-1 cosmological box down to redshift $z=3$. The high mass resolution of up to $142 M_\odot$ per cell in the gas phase combined with a multiphase model of the interstellar medium (ISM), radiative transfer including Lyman-Werner radiation, dust physics, and a non-equilibrium chemistry network that tracks molecular hydrogen, allows for a realistic but still approximate description of Pop III star formation in pristine gas. Our results show that Pop III stars continue to form in low-mass haloes ranging from $10^6 M_\odot$ to $10^9 M_\odot$ until the end of reionization at around $z=5$. At this stage, photoevaporation suppresses further star formation in these minihaloes, which subsequently merge into larger central haloes. Hence, the remnants of Pop III stars primarily reside in the satellite galaxies of larger haloes at lower redshifts. While direct detection of Pop III stars remains elusive, these results hint that lingering primordial star formation could leave observable imprints or indirectly affect the properties of high-redshift galaxies. Explicit Pop III feedback and specialized initial mass function modelling within the THESAN-ZOOM framework would further help interpreting emerging constraints from the James Webb Space Telescope.

arXiv:2503.03964v1 Announce Type: new Abstract: We present a new pipeline designed for the robust inference of cosmological parameters using both second- and third-order shear statistics. We build a theoretical model for rapid evaluation of three-point correlations using our fastnc code and integrate it into the CosmoSIS framework. We measure the two-point functions $\xi_{\pm}$ and the full configuration-dependent three-point shear correlation functions across all auto- and cross-redshift bins. We compress the three-point functions into the mass aperture statistic $\langle M_{\rm ap}^3\rangle$ for a set of 796 simulated shear maps designed to model the Dark Energy Survey (DES) Year 3 data. We estimate from it the full covariance matrix and model the effects of intrinsic alignments, shear calibration biases and photometric redshift uncertainties. We apply scale cuts to minimize the contamination from the baryonic signal as modeled through hydrodynamical simulations. We find a significant improvement of $83\%$ on the Figure of Merit in the $\Omega_{\rm m}$-$S_8$ plane when we add the $\langle M_{\rm ap}^3\rangle$ data to the $\xi_{\pm}$ information. We present our findings for all relevant cosmological and systematic uncertainty parameters and discuss the complementarity of third-order and second-order statistics.

arXiv:2503.03829v1 Announce Type: new Abstract: We measure H$\alpha$ luminosity functions (LFs) at redshifts $z \sim 4.5$ and 6.3 using the JWST MAGNIF (Medium-band Astrophysics with the Grism of NIRCam In Frontier fields) survey. MAGNIF obtained NIRCam grism spectra with the F360M and F480M filters in four Frontier Fields. We identify 248 H$\alpha$ emitters based on the grism spectra and photometric redshifts from combined HST and JWST imaging data. The numbers of the H$\alpha$ emitters show a large field-to-field variation, highlighting the necessity of multiple fields to mitigate cosmic variance. We calculate both observed and dust-corrected H$\alpha$ LFs in the two redshift bins. Thanks to the gravitational lensing, the measured H$\alpha$ LFs span three orders of magnitude in luminosity, and the faint-end luminosity reaches $L_{\mathrm{H}\alpha} \sim 10^{40.3} \mathrm{erg} \mathrm{s}^{-1}$ at $z \sim 4.5$ and $10^{41.5} \mathrm{erg} \mathrm{s}^{-1}$ at $z \sim 6.3$, corresponding to star-formation rates (SFRs) of $\sim$ 0.1 and 1.7 $\mathrm{M}_\odot \mathrm{yr}^{-1}$. We conclude no or weak redshift evolution of the faint-end slope of H$\alpha$ LF across $z\simeq0.4-6.3$, and the comparison with the faint-end slopes of UV LF indicates stochastic star formation history among low-mass H$\alpha$ emitters. The derived cosmic SFR densities are $0.058^{+0.008}_{-0.006}\ \ M_\odot\ \mathrm{yr}^{-1}\ \mathrm{Mpc}^{-3}$ at $z \sim 4.5$ and $0.025^{+0.009}_{-0.007}\ \ M_\odot\ \mathrm{yr}^{-1}\ \mathrm{Mpc}^{-3}$ at $z \sim 6.3$. These are approximately 2.2 times higher than previous estimates based on dust-corrected UV LFs, but consistent with recent measurements from infrared surveys. We discuss uncertainties in the H$\alpha$ LF measurements, including those propagate from the lens models, cosmic variance, and AGN contribution.

arXiv:2409.20533v2 Announce Type: replace Abstract: JADES-GS-z14-0 is the most distant spectroscopically confirmed galaxy yet, at $z>14$. With a UV magnitude of -20.81, it is one of the most luminous galaxies at cosmic dawn and its half-light radius of 260 pc means that stars dominate the observed UV emission. We report ALMA detection of [OIII]88$\mu$m line emission with a significance of 6.67$\sigma$ and at a frequency of 223.524~GHz, corresponding to a redshift of $14.1796\pm0.0007$, which is consistent with the candidate CIII] line detected in the NIRSpec spectrum. At this spectroscopic redshift, the Lyman-$\alpha$ break identified with NIRSpec requires a damped Lyman-$\alpha$ absorber with a column density of $\log(N_{\rm HI}/\mathrm{cm}^{-2})=21.96$. The total [O\,{\sc iii}]88$\mu$m luminosity (log$(L_{\rm [OIII]}/L_\odot) = 8.3\pm 0.1$) is fully consistent with the local $L_{\rm [OIII]}-SFR$ relation and indicating a gas-phase metallicity $>0.1~{\rm Z_{\rm \odot}}$. Using \texttt{prospector} SED modeling and combining the ALMA data with JWST observations, we find $Z=0.17~{\rm Z_{\rm \odot}}$ and a non-zero escape fraction of ionizing photons ($\sim11\%$), which is necessary by the code to reproduce the UV spectrum. We measure an ${\rm [O III]}5007$\r{A}/[O III]88$\mu$m line flux ratio between 1 and 20, resulting in an upper limit to the electron density of roughly 700 cm$^{-3}$ assuming a single-cloud photoionization model. The [OIIII]88$\mu$m emission line is spectrally resolved, with a FWHM of 100 km/s, resulting in a dynamical mass of log($M_{\rm dyn}/M_\odot$) = 9.0$\pm0.2$. When compared to the stellar mass, this value represents a conservative upper limit on the gas mass fraction, which ranges from 50\% to 80\%, depending on the assumed star formation history. Past radiation-driven outflows may have cleared the galaxy from the gas, reducing the gas fraction and thus increasing the escape fraction of ionizing photons.

arXiv:2503.02927v1 Announce Type: new Abstract: We investigate the impact of ionizing external ultraviolet (UV) radiation on low-mass haloes ($M_{h}<10^{10}M_\odot$) at high redshift using $1140M_\odot$ baryonic resolution zoom-in simulations of seven regions from the THESAN-ZOOM project. We compare three simulation sets that differ in the treatment of external UV radiation: one employing a uniform UV background initiated at z=10.6 in addition to radiation transport for local sources, another with the same background starting at z=5.5, and the default configuration in which the large-scale radiation field from the parent THESAN-1 simulation box acts as a boundary condition. The multi-phase interstellar medium (ISM) model, combined with its high mass resolution, allows us to resolve all star-forming haloes and capture the back-reaction of ionizing radiation on galaxy properties during the epoch of reionization. When present, external UV radiation efficiently unbinds gas in haloes with masses below $10^9M_\odot$ and suppresses subsequent star formation. As a result, in simulations with early reionization, minihaloes fail to form stars from pristine gas, leading to reduced metal enrichment of gas later accreted by more massive haloes. Consequently, haloes with masses below $10^{10}M_\odot$ at all simulated epochs (z>3) exhibit lower metallicities and altered metallicity distributions. The more accurate and realistic shielding from external UV radiation, achieved through self-consistent radiative transfer, permits the existence of a cold but low-density gas phase down to z=3. These findings highlight the importance of capturing a patchy reionization history in high-resolution simulations targeting high-redshift galaxy formation. We conclude that at minimum, a semi-numerical model that incorporates spatially inhomogeneous reionization and a non-uniform metallicity floor is necessary to accurately emulate metal enrichment in minihaloes.

arXiv:2503.01949v1 Announce Type: new Abstract: Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($\epsilon^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $\epsilon^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $\epsilon^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $\epsilon^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

arXiv:2503.00919v1 Announce Type: new Abstract: The Largest Cluster Statistics\,(LCS) analysis of the redshifted 21\,cm maps has been demonstrated to be an efficient and robust method for following the time evolution of the largest ionized regions\,(LIRs) during the Epoch of Reionization\,(EoR). The LCS can, in principle, constrain the reionization model and history by quantifying the morphology of neutral hydrogen\,(\HI) distribution during the different stages of the EoR. Specifically, the percolation transition of ionized regions, quantified and constrained via LCS, provides a crucial insight about the underlying reionization model. The previous LCS analysis of EoR 21\,cm maps demonstrates that the convolution of the synthesized beam of the radio interferometric arrays, e.g. SKA1-Low with the target signal, shifts the apparent percolation transition of ionized regions towards the lower redshifts. In this study, we present an optimal thresholding strategy to reduce this bias in the recovered percolation transition. We assess the robustness of LCS analysis of the 21\,cm maps in the presence of antenna-based gain calibration errors and instrumental noise for SKA1-Low. This analysis is performed using synthetic observations simulated by the \textsc{21cmE2E} pipeline, considering SKA1-Low AA4 configuration within a radius of 2\,km from the array centre. Our findings suggest that a minimum of $1500$\,hours of observation (SNR $\gtrapprox 3$) are required for the LCS analysis to credibly suppress the confusion introduced by thermal noise. Further, we also demonstrate that for a maximum antenna-based calibration error tolerance of $\sim 0.05\%$ (post calibration), the reionization history can be recovered in a robust and relatively unbiased manner using the LCS.

arXiv:2503.00636v1 Announce Type: new Abstract: We describe updated scientific goals for the wide-field, millimeter-wave survey that will be produced by the Simons Observatory (SO). Significant upgrades to the 6-meter SO Large Aperture Telescope (LAT) are expected to be complete by 2028, and will include a doubled mapping speed with 30,000 new detectors and an automated data reduction pipeline. In addition, a new photovoltaic array will supply most of the observatory's power. The LAT survey will cover about 60% of the sky at a regular observing cadence, with five times the angular resolution and ten times the map depth of Planck. The science goals are to: (1) determine the physical conditions in the early universe and constrain the existence of new light particles; (2) measure the integrated distribution of mass, electron pressure, and electron momentum in the late-time universe, and, in combination with optical surveys, determine the neutrino mass and the effects of dark energy via tomographic measurements of the growth of structure at $z < 3$; (3) measure the distribution of electron density and pressure around galaxy groups and clusters, and calibrate the effects of energy input from galaxy formation on the surrounding environment; (4) produce a sample of more than 30,000 galaxy clusters, and more than 100,000 extragalactic millimeter sources, including regularly sampled AGN light-curves, to study these sources and their emission physics; (5) measure the polarized emission from magnetically aligned dust grains in our Galaxy, to study the properties of dust and the role of magnetic fields in star formation; (6) constrain asteroid regoliths, search for Trans-Neptunian Objects, and either detect or eliminate large portions of the phase space in the search for Planet 9; and (7) provide a powerful new window into the transient universe on time scales of minutes to years, concurrent with observations from Rubin of overlapping sky.

arXiv:2503.00106v1 Announce Type: new Abstract: Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is typically inferred from observations. We can reproduce the flatter observed evolution by filtering out weakly star-forming galaxies, implying that current observational fits are biased due to a missing population of lulling galaxies or overestimated star-formation rates. We also explore star-formation variability using the scatter of galaxies around the SFMS ($\sigma_{\mathrm{MS}}$). At the population level, the scatter around the SFMS increases with cosmic time, driven by the increased importance of long-term environmental effects in regulating star formation at later times. To study short-term star-formation variability, or ''burstiness'', we isolate the scatter on timescales shorter than 50 Myr. The short-term scatter is larger at higher redshift, indicating that star formation is indeed more bursty in the early Universe. We identify two starburst modes: (i) externally driven, where rapid large-scale inflows trigger and fuel prolonged, extreme star formation episodes, and (ii) internally driven, where cyclical ejection and re-accretion of the interstellar medium in low-mass galaxies drive bursts, even under relatively steady large-scale inflow. Both modes occur at all redshifts, but the increased burstiness of galaxies at higher redshift is due to the increasing prevalence of the more extreme external mode of star formation.

arXiv:2502.20741v1 Announce Type: new Abstract: We present a follow-up study to the method recently proposed by Namikawa & Sherwin (2023) to probe gravitational waves using cross-correlations between two CMB $B$-modes and a large-scale structure tracer. We first improve on the previous forecast by including the impact of CMB component separation and find that we can achieve $\sigma_r\simeq3.5\times10^{-3}$ by combining upcoming experiments, i.e., LiteBIRD, CMB-S4 and the Advanced Simons Observatory. With a more futuristic experiment, we can achieve even tighter constraints on $r$ if improved delensing can be realized. Using a simulated analysis pipeline, we also explore possible biases from higher-order terms in the lensing potential, which were previously not examined in detail. We find that these bias terms are negligible compared to a detectable signal from inflationary gravitational waves. Our simulated results confirm that this new method is capable of obtaining powerful constraints on $r$. The method is immune to Gaussian Galactic foregrounds and has a different response to non-Gaussian Galactic foregrounds than the $B$-mode power spectrum, offering an independent cross-check of $r$ constraints from the standard power spectrum analysis.

arXiv:2408.15615v2 Announce Type: replace Abstract: In this work, we spectroscopically select narrow-line AGN (NLAGN) among the $\sim 300$ publicly available medium-resolution spectra of the CEERS Survey. Using both traditional and newly identified emission line NLAGN diagnostics diagrams, we identified 52 NLAGN at $2\lesssim z\lesssim 9$ on which we performed a detailed multiwavelength analysis. We also identified 4 new $z\lesssim 2$ broad-line AGN (BLAGN), in addition to the 8 previously reported high-$z$ BLAGN. We found that the traditional BPT diagnostic diagrams are not suited to identify high-$z$ AGN, while most of the high-$z$ NLAGN are selected using the recently proposed AGN diagnostic diagrams based on the [OIII]$\lambda$4363 auroral line or high-ionization emission lines. We compared the emission line velocity dispersion and the obscuration of the sample of NLAGN with those of the parent sample without finding significant differences between the two distributions, suggesting a population of AGN heavily buried and not significantly impacting the host galaxies' physical properties, as further confirmed by SED-fitting. The bolometric luminosities of the high-$z$ NLAGNs selected in this work are well below those sampled by surveys before JWST, potentially explaining the weak impact of these AGN. Finally, we investigate the X-ray properties of the selected NLAGN and of the sample of high-$z$ BLAGN. We find that all but 4 NLAGN are undetected in the deep X-ray image of the field, as well as all the high-$z$ BLAGN. We do not obtain a detection even by stacking the undetected sources, resulting in an X-ray weakness of $\sim 1-2$ dex from what is expected based on their bolometric luminosities. To discriminate between a heavily obscured AGN scenario or an intrinsic X-ray weakness of these sources, we performed a radio stacking analysis, which did not reveal any detection leaving open the questions about the origin of the X-ray weakness.

arXiv:2406.10348v2 Announce Type: replace Abstract: We present JWST NIRSpec observations in IFS mode of the galaxy GS5001 at redshift z=3.47, the central member of a candidate protocluster in the GOODS-S field. The data cover a field of view (FoV) of 4''$\times$4'' (~$30\times30$~kpc$^2$) and were obtained as part of the GA-NIFS GTO program. The observations include both high (R~2700) and low (R~100) spectral resolution data, spanning the rest-frame wavelength ranges 3700-6780A and 1300-11850A, respectively. We analyse the spatially resolved ionised gas kinematics and interstellar medium properties, including obscuration, gas metallicity, excitation, ionisation parameter, and electron density. In addition to the main galaxy (GS5001), the NIRSpec FoV covers three components in the south, with velocities blue-shifted by -150 km/s with respect to the main galaxy, and another source in the north redshifted by ~200 km/s. Optical line ratio diagnostics indicate star formation ionisation and electron densities of ~500 cm$^{-3}$ across all sources in the FoV. The gas-phase metallicity in the main galaxy is 12+log(O/H) $= 8.45\pm0.04$, and slightly lower in the companions (12+log(O/H)$ = 8.34-8.42$), consistent with the mass-metallicity relation at $z\sim3$. We find peculiar line ratios (high log [NII]/H$\alpha$, low log [OIII]/H$\beta$) in the northern part of GS5001. These could be attributed to either higher metallicity, or to shocks resulting from the interaction of the main galaxy with the northern source. We identify a spatially resolved outflow in the main galaxy, with an extension of about 3 kpc. We find maximum outflow velocities of ~400 km/s, an outflow mass of $(1.7\pm0.4)\times 10^8$ M$_{\odot}$, a mass outflow rate of $23\pm5$ M$_{\odot}$ yr$^{-1}$ and a mass loading factor of 0.23. These properties are compatible with star formation being the driver of the outflow.

arXiv:2502.20437v1 Announce Type: new Abstract: We introduce the THESAN-ZOOM project, a comprehensive suite of high-resolution zoom-in simulations of $14$ high-redshift ($z>3$) galaxies selected from the THESAN simulation volume. This sample encompasses a diverse range of halo masses, with $M_\mathrm{halo} \approx 10^8 - 10^{13}~\mathrm{M}_\odot$ at $z=3$. At the highest-resolution, the simulations achieve a baryonic mass of $142~\mathrm{M}_\odot$ and a gravitational softening length of $17~\mathrm{cpc}$. We employ a state-of-the-art multi-phase interstellar medium (ISM) model that self-consistently includes stellar feedback, radiation fields, dust physics, and low-temperature cooling through a non-equilibrium thermochemical network. Our unique framework incorporates the impact of patchy reionization by adopting the large-scale radiation field topology from the parent THESAN simulation box rather than assuming a spatially uniform UV background. In total, THESAN-ZOOM comprises $60$ simulations, including both fiducial runs and complementary variations designed to investigate the impact of numerical and physical parameters on galaxy properties. The fiducial simulation set reproduces a wealth of high-redshift observational data such as the stellar-to-halo-mass relation, the star-forming main sequence, the Kennicutt-Schmidt relation, and the mass-metallicity relation. While our simulations slightly overestimate the abundance of low-mass and low-luminosity galaxies they agree well with observed stellar and UV luminosity functions at the higher mass end. Moreover, the star-formation rate density closely matches the observational estimates from $z=3-14$. These results indicate that the simulations effectively reproduce many of the essential characteristics of high-redshift galaxies, providing a realistic framework to interpret the exciting new observations from JWST.

arXiv:2502.19484v1 Announce Type: new Abstract: Using data from the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), we tentatively detect two molecular hydrogen (H$_2$) fluorescent emission features in high-redshift galaxies at 2.3 and 3.1$\sigma$. These features consist of many blended emission lines that result from the de-excitation cascade of H$_2$ molecules that have absorbed Lyman-Werner band photons. Our study targets galaxies at redshifts $z \geq 7$ galaxies of the early Universe, as they host some of the most extreme conditions in terms of star formation, molecular gas content, and the possible presence of outflows driven by starbursts and active galactic nuclei. To enhance the signal-to-noise ratio of H$_2$ emission features in the rest-frame wavelength range of 155-191 nm, we stack JWST/NIRSpec spectra from $z=7.0-13.4$. These stacked spectra also exhibit atomic emission features, such as CIV emission with a P-Cygni profile, as well as CIII], OIII], and HeII. The presence of these features and the slightly blue-shifted fluorescent H$_2$ lines suggest active multiphase atomic and molecular outflows may be common at these redshifts. Our results suggest that it is possible to search for FUV fluorescent H$_2$ lines in high redshift galaxies using JWST, which would enable the characterization of interstellar radiation fields, densities, and temperatures in extreme photodissociation environments.

arXiv:2409.11464v2 Announce Type: replace Abstract: For the first time, we present a systematic search for galaxies with extended emission line and potential outflow features using \textit{JWST} medium-band images in the GOODS-S field. This is done by comparing the morphology in medium-band images to adjacent continuum and UV bands. We look for galaxies that have a maximum extent 50\% larger, an excess area 30\% greater, or an axis ratio difference of more than 0.3 in the medium band compared to the reference bands. After visual inspection, we find 326 candidate galaxies at $1.4 < z < 8.4$, with a peak in the population near cosmic noon, benefiting from the good coverage of the medium-band filters. By fitting their SEDs, we find that the candidate galaxies are at least 20\% more bursty in their star-forming activity and have 50\% more young stellar populations compared to a control sample selected based on the continuum band flux. Additionally, these candidates exhibit a significantly higher production rate of ionizing photons. We further find that candidates hosting known AGN produce extended emission that is more anisotropic compared to non-AGN candidates. A few of our candidates have been spectroscopically confirmed to have prominent outflow signatures through NIRSpec observations, showcasing the robustness of the photometric selection. Future spectroscopic follow-up will better help verify and characterize the kinematics and chemical properties of these systems.