-
Inferring the Ionizing Photon Contributions of High-Redshift Galaxies to Reionization with JWST NIRCam Photometry
Authors:
Nicholas Choustikov,
Richard Stiskalek,
Aayush Saxena,
Harley Katz,
Julien Devrient,
Adrianne Slyz
Abstract:
JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergala…
▽ More
JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergalactic medium. Observational estimates of these values generally rely on spectroscopy for which large samples with well-defined selection functions remain limited. To overcome this challenge, we present an implicit likelihood inference (ILI) pipeline trained on mock photometry to predict the escaped ionizing luminosity of individual galaxies ($\dot{n}_{\rm ion}$) based on photometric magnitudes and redshifts. Compared to traditional SED-fitting methods, the new ILI pipeline is consistently more accurate and significantly faster. We deploy the method on a sample of 4,559 high-redshift galaxies from the JADES Deep survey, finding a gentle redshift evolution of $\log_{10}(\dot{n}_{\rm ion}) = (0.08\pm0.01)z + (51.60\pm0.06)$, with late-time values for $\dot{N}_{\rm ion}$ consistent with theoretical models and observations. We measure the evolution of the volume-averaged ionized fraction and optical depth to find that observed populations of star-forming galaxies are capable of driving reionization to completion at $z\sim 5.3$ without the need for AGN or other exotic sources. The $20\%$ of UV-brightest galaxies ($M_{\rm UV}<-18.5$) in our sample can reionize only $\sim30\%$ of the survey volume, demonstrating that faint LyC emitters are crucial for reionization.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
The sizes of bright Lyman-break galaxies at $z\simeq3-5$ with JWST PRIMER
Authors:
R. G. Varadaraj,
R. A. A. Bowler,
M. J. Jarvis,
N. J. Adams,
N. Choustikov,
A. M. Koekemoer,
A. C. Carnall,
D. J. McLeod,
J. S. Dunlop,
C. T. Donnan,
N. A. Grogin
Abstract:
We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distribution…
▽ More
We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distributions of luminous sources. We fit Sérsic profiles to eight NIRCam bands and also measure a non-parametric half-light radius. We find that the size distributions with both measurements are well-fit by a log-normal distribution at all redshifts, consistent with disk formation models where size is governed by host dark-matter halo angular momentum. We find a size-redshift evolution of $R_{e} = 3.51(1+z)^{-0.60\pm0.22}$ kpc, in agreement with JWST studies. When considering the typical (modal) size over $z=3-5$, we find little evolution with bright LBGs remaining compact at $R_{e}\simeq0.7-0.9$ kpc. Simultaneously, we find evidence for a build-up of large ($R_{e} > 2$ kpc) galaxies by $z=3$. We find some evidence for a negatively sloped size-mass relation at $z=5$ when Sérsic profiles are used to fit the data in F200W. The intrinsic scatter in our size-mass relations increases at higher redshifts. Additionally, measurements probing the rest-UV (F200W) show larger intrinsic scatter than those probing the rest-optical (F356W). Finally, we leverage rest-UV and rest-optical photometry to show that disky galaxies are well established by $z=5$, but are beginning to undergo dissipative processes, such as mergers, by $z=3$. The agreement of our size-mass and size-luminosity relations with simulations provides tentative evidence for centrally concentrated star formation at high-redshift.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
The Great Escape: Understanding the Connection Between Ly$α$ Emission and LyC Escape in Simulated JWST Analogues
Authors:
Nicholas Choustikov,
Harley Katz,
Aayush Saxena,
Thibault Garel,
Julien Devriendt,
Adrianne Slyz,
Taysun Kimm,
Jeremy Blaizot,
Joki Rosdahl
Abstract:
Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly$α$ emission line correlate with the inferred LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) of low-redshift galaxies. Using a data-set of 9,600 mock Ly$α$ spectra of star-forming galaxies…
▽ More
Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly$α$ emission line correlate with the inferred LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) of low-redshift galaxies. Using a data-set of 9,600 mock Ly$α$ spectra of star-forming galaxies at $4.64 \leq z \leq 6$ from the SPHINX$^{20}$ cosmological radiation hydrodynamical simulation, we study the escape of Ly$α$ and LyC photons. We find that our mock Ly$α$ observations are representative of high-redshift observations and that typical observational methods tend to over-predict the Ly$α$ escape fraction ($f_{\rm esc}^{\rm Lyα}$) by as much as two dex. We investigate the correlations between $f_{\rm esc}^{\rm LyC}$ and $f_{\rm esc}^{\rm Lyα}$, Ly$α$ equivalent width ($W_λ({\rm Lyα})$), peak separation ($v_{\rm sep}$), central escape fraction ($f_{\rm cen}$), and red peak asymmetry ($A_f^{\rm red}$). We find that $f_{\rm esc}^{\rm Lyα}$ and $f_{\rm cen}$ are good diagnostics for LyC leakage, selecting for galaxies with lower neutral gas densities and less UV attenuation that have recently experienced supernova feedback. In contrast, $W_λ({\rm Lyα})$ and $v_{\rm sep}$ are found to be necessary but insufficient diagnostics, while $A_f^{\rm red}$ carries little information. Finally, we use stacks of Ly$α$, H$α$, and F150W mock surface brightness profiles to find that galaxies with high $f_{\rm esc}^{\rm LyC}$ have less extended Ly$α$ and F150W but larger H$α$ haloes than their non-leaking counterparts. This confirms that Ly$α$ spectral profiles and surface brightness morphology can be used to better understand the escape of LyC photons from galaxies during the Epoch of Reionization.
△ Less
Submitted 24 June, 2024; v1 submitted 17 January, 2024;
originally announced January 2024.
-
The Sphinx Public Data Release: Forward Modelling High-Redshift JWST Observations with Cosmological Radiation Hydrodynamics Simulations
Authors:
Harley Katz,
Joki Rosdahl,
Taysun Kimm,
Jeremy Blaizot,
Nicholas Choustikov,
Marion Farcy,
Thibault Garel,
Martin G. Haehnelt,
Leo Michel-Dansac,
Pierre Ocvirk
Abstract:
The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medi…
▽ More
The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medium (ISM) in high-redshift galaxies. To this end, we present Version 1 of the Sphinx$^{20}$ public data release. Sphinx$^{20}$ is a full box cosmological radiation hydrodynamics simulation that simultaneously models the large-scale process of cosmic reionization and the detailed physics of a multiphase ISM, providing a statistical sample of galaxies akin to those currently being observed by JWST. The data set contains $\sim14,000$ mock images and spectra of the stellar continuum, nebular continuum, and 52 nebular emission lines, including Ly$α$, for each galaxy in Sphinx$^{20}$ with a star formation rate $\geq0.3\ {\rm M_{\odot}\ yr^{-1}}$. All galaxy emission has been processed with dust radiative transfer and/or resonant line radiative transfer, and data is provided for ten viewing angles for each galaxy. Additionally, we provide a comprehensive set of intrinsic galaxy properties, including halo masses, stellar masses, star formation histories, and ISM characteristics (e.g., metallicity, ISM gas densities, LyC escape fractions). This paper outlines the data generation methods, presents a comparative analysis with JWST ERS and Cycle 1 observations, and addresses data set limitations. The Sphinx$^{20}$ data release can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data
△ Less
Submitted 3 December, 2023; v1 submitted 6 September, 2023;
originally announced September 2023.
-
Optimizing the Evolution of Perturbations in the $Λ$CDM Universe
Authors:
Nicholas Choustikov,
Zvonimir Vlah,
Anthony Challinor
Abstract:
Perturbation theory is a powerful tool for studying large-scale structure formation in the universe and calculating observables such as the power spectrum or bispectrum. However, beyond linear order, typically this is done by assuming a simplification in the time-dependence of gravitational-coupling kernels between the matter and velocity fluctuations. Though the true dependencies are known for La…
▽ More
Perturbation theory is a powerful tool for studying large-scale structure formation in the universe and calculating observables such as the power spectrum or bispectrum. However, beyond linear order, typically this is done by assuming a simplification in the time-dependence of gravitational-coupling kernels between the matter and velocity fluctuations. Though the true dependencies are known for Lambda cold dark matter cosmologies, they are ignored due to the computational costs associated with considering them in full and, instead, are replaced by simpler dependencies valid for an Einstein--de-Sitter cosmology. Here we develop, implement and demonstrate the effectiveness of a new numerical method for finding the full dynamical evolution of these kernels to all perturbative orders based upon spectral methods using Chebyshev polynomials. This method is found to be orders of magnitude more efficient than direct numerical solvers while still producing highly accurate and reliable results. A code implementation of the Chebyshev spectral method is then presented and characterised. The code has been made publicly available alongside this paper. We expect our method to be of use for interpretation of upcoming galaxy clustering measurements.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
The Physics of Indirect Estimators of Lyman Continuum Escape and their Application to High-Redshift JWST Galaxies
Authors:
Nicholas Choustikov,
Harley Katz,
Aayush Saxena,
Alex J. Cameron,
Julien Devriendt,
Adrianne Slyz,
Joki Rosdahl,
Jeremy Blaizot,
Leo Michel-Dansac
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_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ''analog'' galaxies, it remains unclear whether these are applicable at h…
▽ More
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_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ''analog'' 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 \leq z \leq 10$ from the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high $f_{\rm 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. ${\rm O}_{32}$, $Σ_{\rm SFR}$, M$_{\rm 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 ${\rm E(B-V)}$ match two or more of our criteria, rendering them good $f_{\rm esc}$ diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting $f_{\rm esc}$ based on direct observables. When applied to bright $z > 6$ Ly$α$ emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10\%$.
△ Less
Submitted 13 March, 2024; v1 submitted 17 April, 2023;
originally announced April 2023.
-
The Einstein Toolkit: A Student's Guide
Authors:
Nicholas Choustikov
Abstract:
The Einstein Toolkit represents a unique opportunity for students to explore the world of numerical relativity, without the need for high-level computing power or knowledge of the mathematics behind the simulations themselves. This document aims to act as a resource for students aiming to learn about and get involved with Numerical Relativity, focusing specifically on the modelling of Binary Black…
▽ More
The Einstein Toolkit represents a unique opportunity for students to explore the world of numerical relativity, without the need for high-level computing power or knowledge of the mathematics behind the simulations themselves. This document aims to act as a resource for students aiming to learn about and get involved with Numerical Relativity, focusing specifically on the modelling of Binary Black Hole Mergers in order to produce sample Gravitational Waveforms as would be measured by Gravitational Wave detectors such as LIGO.
The feasibility of using the Einstein Toolkit on Amazon's Web Services (AWS) has also been explored for the first time, confirming a potential - yet expensive - method for students or citizen scientists to gain access to high-level computing power sufficient to play with and run their own numerical relativity simulations.
△ Less
Submitted 26 November, 2020;
originally announced November 2020.
