-
LtU-ILI: An All-in-One Framework for Implicit Inference in Astrophysics and Cosmology
Authors:
Matthew Ho,
Deaglan J. Bartlett,
Nicolas Chartier,
Carolina Cuesta-Lazaro,
Simon Ding,
Axel Lapel,
Pablo Lemos,
Christopher C. Lovell,
T. Lucas Makinen,
Chirag Modi,
Viraj Pandya,
Shivam Pandey,
Lucia A. Perez,
Benjamin Wandelt,
Greg L. Bryan
Abstract:
This paper presents the Learning the Universe Implicit Likelihood Inference (LtU-ILI) pipeline, a codebase for rapid, user-friendly, and cutting-edge machine learning (ML) inference in astrophysics and cosmology. The pipeline includes software for implementing various neural architectures, training schemata, priors, and density estimators in a manner easily adaptable to any research workflow. It i…
▽ More
This paper presents the Learning the Universe Implicit Likelihood Inference (LtU-ILI) pipeline, a codebase for rapid, user-friendly, and cutting-edge machine learning (ML) inference in astrophysics and cosmology. The pipeline includes software for implementing various neural architectures, training schemata, priors, and density estimators in a manner easily adaptable to any research workflow. It includes comprehensive validation metrics to assess posterior estimate coverage, enhancing the reliability of inferred results. Additionally, the pipeline is easily parallelizable and is designed for efficient exploration of modeling hyperparameters. To demonstrate its capabilities, we present real applications across a range of astrophysics and cosmology problems, such as: estimating galaxy cluster masses from X-ray photometry; inferring cosmology from matter power spectra and halo point clouds; characterizing progenitors in gravitational wave signals; capturing physical dust parameters from galaxy colors and luminosities; and establishing properties of semi-analytic models of galaxy formation. We also include exhaustive benchmarking and comparisons of all implemented methods as well as discussions about the challenges and pitfalls of ML inference in astronomical sciences. All code and examples are made publicly available at https://github.com/maho3/ltu-ili.
△ Less
Submitted 2 July, 2024; v1 submitted 6 February, 2024;
originally announced February 2024.
-
Field-level simulation-based inference with galaxy catalogs: the impact of systematic effects
Authors:
Natalí S. M. de Santi,
Francisco Villaescusa-Navarro,
L. Raul Abramo,
Helen Shao,
Lucia A. Perez,
Tiago Castro,
Yueying Ni,
Christopher C. Lovell,
Elena Hernandez-Martinez,
Federico Marinacci,
David N. Spergel,
Klaus Dolag,
Lars Hernquist,
Mark Vogelsberger
Abstract:
It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocit…
▽ More
It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocities of galaxies that are robust to uncertainties in astrophysics and subgrid models. However, observations are affected by many effects, including 1) masking, 2) uncertainties in peculiar velocities and radial distances, and 3) different galaxy selections. Moreover, observations only allow us to measure redshift, intertwining galaxies' radial positions and velocities. In this paper we train and test our models on galaxy catalogs, created from thousands of state-of-the-art hydrodynamic simulations run with different codes from the CAMELS project, that incorporate these observational effects. We find that, although the presence of these effects degrades the precision and accuracy of the models, and increases the fraction of catalogs where the model breaks down, the fraction of galaxy catalogs where the model performs well is over 90 %, demonstrating the potential of these models to constrain cosmological parameters even when applied to real data.
△ Less
Submitted 9 May, 2024; v1 submitted 23 October, 2023;
originally announced October 2023.
-
Constraints on the Epoch of Reionization with Roman Space Telescope and the Void Probability Function of Lyman-Alpha Emitters
Authors:
Lucia A. Perez,
Sangeeta Malhotra,
James E. Rhoads,
Isak G. B. Wold
Abstract:
We use large simulations of Lyman-Alpha Emitters with different fractions of ionized intergalactic medium to quantify the clustering of Ly$α$ emitters as measured by the Void Probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace R…
▽ More
We use large simulations of Lyman-Alpha Emitters with different fractions of ionized intergalactic medium to quantify the clustering of Ly$α$ emitters as measured by the Void Probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace Roman Space Telescope. Since Roman will be able to carry out blind spectroscopic surveys of Ly$α$ emitters continuously between $7<z<12$ to sensitivities of at least $10^{-17}$ erg sec$^{-1}$ over a wide field of view, it can measure the epoch of reionization as well as the pace of ionization of the intergalactic medium (IGM). We compare deep Roman surveys covering roughly 1, 4, and 16 deg$^2$, and quantify what constraints on reionization the VPF may find for these surveys. A survey of 1 deg$^2$ would distinguish between very late reionization and early reionization to 3$σ$ near $z=7.7$ with the VPF. The VPF of a 4 deg$^2$ survey can distinguish between slow vs.\ fast, and early vs.\ late, reionization at $> 3-4σ$ at several redshifts between $7<z<9$. However, a survey of 13-16 deg$^2$ would allow the VPF to give several robust constraints ($>5-8σ$) across the epoch of reionization, and would yield a detailed history of the reionization of the IGM and its effect on Lyman-$α$ Emitter clustering.
△ Less
Submitted 4 April, 2023;
originally announced April 2023.
-
Probing Patchy Reionization with the Void Probability Function of Lyman-$α$ Emitters
Authors:
Lucia A. Perez,
Sangeeta Malhotra,
James E. Rhoads,
Peter Laursen,
Isak G. B. Wold
Abstract:
We probe what constraints for the global ionized hydrogen fraction the Void Probability Function (VPF) clustering can give for the Lyman-Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyman-alpha emission, and measuring the drop in the luminosity function of Lyman-$α$ emitters (LAEs) has been used to constrain the i…
▽ More
We probe what constraints for the global ionized hydrogen fraction the Void Probability Function (VPF) clustering can give for the Lyman-Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyman-alpha emission, and measuring the drop in the luminosity function of Lyman-$α$ emitters (LAEs) has been used to constrain the ionization fraction in narrowband surveys. However, the clustering of LAEs is independent from the luminosity function's inherent evolution, and can offer additional constraints for reionization under different models. The VPF measures how likely a given circle is to be empty. It is a volume-averaged clustering statistic that traces the behavior of higher order correlations, and its simplicity offers helpful frameworks for planning surveys.
Using the \citet{Jensen2014} simulations of LAEs within various amount of ionized intergalactic medium, we predict the behavior of the VPF in one (301x150.5x30 Mpc$^3$), four (5.44$\times 10^6$ Mpc$^3$), or eight (1.1$\times 10^7$ Mpc$^3$) fields of LAGER imaging. We examine the VPF at 5 and 13 arcminutes, corresponding to the minimum scale implied by the LAE density and the separation of the 2D VPF from random, and the maximum scale from the 8-field 15.5 deg$^2$ LAGER area. We find that even a single DECam field of LAGER (2-3 deg$^2$) could discriminate between mostly neutral vs. ionized. Additionally, we find four fields allows the distinction between 30, 50, and 95 percent ionized; and that eight fields could even distinguish between 30, 50, 73, and 95 percent ionized.
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
Constraining cosmology with machine learning and galaxy clustering: the CAMELS-SAM suite
Authors:
Lucia A. Perez,
Shy Genel,
Francisco Villaescusa-Navarro,
Rachel S. Somerville,
Austen Gabrielpillai,
Daniel Anglés-Alcázar,
Benjamin D. Wandelt,
L. Y. Aaron Yung
Abstract:
As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophy…
▽ More
As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project: CAMELS-SAM, encompassing one thousand dark-matter only simulations of (100 $h^{-1}$ cMpc)$^3$ with different cosmological parameters ($Ω_m$ and $σ_8$) and run through the Santa Cruz semi-analytic model for galaxy formation over a broad range of astrophysical parameters. As a proof-of-concept for the power of this vast suite of simulated galaxies in a large volume and broad parameter space, we probe the power of simple clustering summary statistics to marginalize over astrophysics and constrain cosmology using neural networks. We use the two-point correlation function, count-in-cells, and the Void Probability Function, and probe non-linear and linear scales across $0.68<$ R $<27\ h^{-1}$ cMpc. Our cosmological constraints cluster around 3-8$\%$ error on $Ω_{\text{M}}$ and $σ_8$, and we explore the effect of various galaxy selections, galaxy sampling, and choice of clustering statistics on these constraints. We additionally explore how these clustering statistics constrain and inform key stellar and galactic feedback parameters in the Santa Cruz SAM. CAMELS-SAM has been publicly released alongside the rest of CAMELS, and offers great potential to many applications of machine learning in astrophysics: https://camels-sam.readthedocs.io.
△ Less
Submitted 22 May, 2023; v1 submitted 5 April, 2022;
originally announced April 2022.
-
The CAMELS project: public data release
Authors:
Francisco Villaescusa-Navarro,
Shy Genel,
Daniel Anglés-Alcázar,
Lucia A. Perez,
Pablo Villanueva-Domingo,
Digvijay Wadekar,
Helen Shao,
Faizan G. Mohammad,
Sultan Hassan,
Emily Moser,
Erwin T. Lau,
Luis Fernando Machado Poletti Valle,
Andrina Nicola,
Leander Thiele,
Yongseok Jo,
Oliver H. E. Philcox,
Benjamin D. Oppenheimer,
Megan Tillman,
ChangHoon Hahn,
Neerav Kaushal,
Alice Pisani,
Matthew Gebhardt,
Ana Maria Delgado,
Joyce Caliendo,
Christina Kreisch
, et al. (22 additional authors not shown)
Abstract:
The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present…
▽ More
The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$α$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}.
△ Less
Submitted 4 January, 2022;
originally announced January 2022.
-
New spectroscopic confirmations of Lyman-$α$ emitters at z $\sim$ 7 from the LAGER survey
Authors:
Santosh Harish,
Isak G. B. Wold,
Sangeeta Malhotra,
James Rhoads,
Weida Hu,
Junxian Wang,
Zhen-ya Zheng,
L. Felipe Barrientos,
Jorge González-López,
Lucia A. Perez,
Ali Ahmad Khostovan,
Leopoldo Infante,
Chunyan Jiang,
Cristóbal Moya-Sierralta,
John Pharo,
Francisco Valdes,
Huan Yang
Abstract:
We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sen…
▽ More
We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sensitive probes of cosmic reionization and narrow-band imaging is a robust and effective method for selecting a large number of LAEs. In this work, we present results from the spectroscopic follow-up of LAE candidates in two LAGER fields, COSMOS and WIDE-12, using observations from Keck/LRIS. We report the successful detection of Ly$α$ emission in 15 candidates (11 in COSMOS and 4 in WIDE-12 fields). Three of these in COSMOS have matching confirmations from a previous LAGER spectroscopic follow-up and are part of the overdense region, LAGER-$z7$OD1. Additionally, two candidates that were not detected in the LRIS observations have prior spectroscopic confirmations from Magellan. Including these, we obtain a spectroscopic confirmation success rate of $\sim$$80$% for LAGER LAE candidates. Apart from Ly$α$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$σ$ upper limit for the ratio of NV/Ly$α$, $f_{NV}/f_{Lyα} \lesssim 0.27$, which implies that ionizing emission from these sources is mostly dominated by star formation. Including confirmations from this work, a total of 33 LAE sources from LAGER are now spectroscopically confirmed. LAGER has more than doubled the sample of spectroscopically confirmed LAE sources at $z\sim7$.
△ Less
Submitted 1 November, 2021;
originally announced November 2021.
-
The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence
Authors:
Francisco Villaescusa-Navarro,
Shy Genel,
Daniel Angles-Alcazar,
Leander Thiele,
Romeel Dave,
Desika Narayanan,
Andrina Nicola,
Yin Li,
Pablo Villanueva-Domingo,
Benjamin Wandelt,
David N. Spergel,
Rachel S. Somerville,
Jose Manuel Zorrilla Matilla,
Faizan G. Mohammad,
Sultan Hassan,
Helen Shao,
Digvijay Wadekar,
Michael Eickenberg,
Kaze W. K. Wong,
Gabriella Contardo,
Yongseok Jo,
Emily Moser,
Erwin T. Lau,
Luis Fernando Machado Poletti Valle,
Lucia A. Perez
, et al. (3 additional authors not shown)
Abstract:
We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light year…
▽ More
We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
LAGER Ly$α$ Luminosity Function at $z\sim7$, Implications for Reionization
Authors:
Isak G. B. Wold,
Sangeeta Malhotra,
James Rhoads,
Junxian Wang,
Weida Hu,
Lucia A. Perez,
Zhen-Ya Zheng,
Ali Ahmad Khostovan,
Alistair R. Walker,
L. Felipe Barrientos,
Jorge González-López,
Santosh Harish,
Leopoldo Infante,
Chunyan Jiang,
John Pharo,
Cristóbal Moya-Sierralta,
Franz E. Bauer,
Gaspar Galaz,
Francisco Valdes,
Huan Yang
Abstract:
We present a new measurement of the Ly$α$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest…
▽ More
We present a new measurement of the Ly$α$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest $z\sim 7$ survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify $N=95$ new $z=6.9$ Ly$α$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$α$ luminosity functions. The best-fit Schechter luminosity function parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Ly$α$ luminosity density declines at the same rate as the UV continuum LF from $z=5.7$ to $z=6.9$. This is consistent with an intergalactic medium that was fully ionized as early as redshift $z\sim 7$, or with a volume-averaged neutral hydrogen fraction of $x_{HI} < 0.33$ at $1σ$.
△ Less
Submitted 20 January, 2022; v1 submitted 25 May, 2021;
originally announced May 2021.
-
Correlations between H$α$ Equivalent Width and Galaxy Properties at $z = 0.47$: Physical or Selection-driven?
Authors:
Ali Ahmad Khostovan,
Sangeeta Malhotra,
James E. Rhoads,
Santosh Harish,
Chunyan Jiang,
Junxian Wang,
Isak Wold,
Zhen-Ya Zheng,
L. Felipe Barrientos,
Alicia Coughlin,
Weida Hu,
Leopoldo Infante,
Lucia A. Perez,
John Pharo,
Francisco Valdes,
Alistair R. Walker
Abstract:
The H$α$ equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic star-formation activity. Previous assessments show that EW strongly anti-correlates with stellar mass as $M^{-0.25}$ similar to the sSFR -- stellar mass relation. However, such a correlation may be driven/formed by selection effects. In this study, we investigate how H$α$ EWs c…
▽ More
The H$α$ equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic star-formation activity. Previous assessments show that EW strongly anti-correlates with stellar mass as $M^{-0.25}$ similar to the sSFR -- stellar mass relation. However, such a correlation may be driven/formed by selection effects. In this study, we investigate how H$α$ EWs correlate with galaxy properties and how selection biases could alter such correlations using a $z = 0.47$ narrowband-selected sample of 1572 H$α$ emitters from the Ly$α$ Galaxies in the Epoch of Reionization (LAGER) survey. The sample covers 3 deg$^2$ of COSMOS and $1.1\times10^5$ cMpc$^3$. We assume an intrinsic EW distribution to form mock samples of H$α$ emitters (HAEs) and propagate the selection criteria to match observations, giving us control on how selection biases can affect the underlying results. We find EW intrinsically correlates with stellar mass as $W_0 \propto M^{-0.16\pm0.03}$ and decreases by a factor of $\sim 3$ from $10^{7}$ to $10^{10}$ M$_\odot$. We find low-mass HAEs to be $\sim 320$ times more likely to have rest-frame EW$ > 200$Åcompared to high-mass HAEs. Combining the intrinsic EW -- stellar mass correlation with an observed SMF correctly reproduces the observed H$α$ LF, while not correcting for selection effects underestimates the number of bright HAEs. This suggests that the intrinsic EW -- stellar mass correlation is physically significant and reproduces three statistical distributions of galaxy populations (LF, SMF, EW distribution). At lower masses, we find there are more high-EW outliers compared to high masses, even after taking into account selection effects. Our results suggest that high sSFR outliers indicative of bursty SF activity are intrinsically more prevalent in low-mass HAEs and not a byproduct of selection effects.
△ Less
Submitted 19 March, 2021;
originally announced March 2021.
-
A Lyman-α protocluster at redshift 6.9
Authors:
Weida Hu,
Junxian Wang,
Leopoldo Infante,
James E. Rhoads,
Zhen-Ya Zheng,
Huan Yang,
Sangeeta Malhotra,
L. Felipe Barrientos,
Chunyan Jiang,
Jorge González-López,
Gonzalo Prieto,
Lucia A. Perez,
Pascale Hibon,
Gaspar Galaz,
Alicia Coughlin,
Santosh Harish,
Xu Kong,
Wenyong Kang,
Ali Ahmad Khostovan,
John Pharo,
Francisco Valdes,
Isak Wold,
Alistair R. Walker,
XianZhong Zheng
Abstract:
Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million ye…
▽ More
Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million years old and could be experiencing rapid evolution of the neutral hydrogen fraction in the intergalactic medium. The protocluster is identified by an overdensity of 6 times the average galaxy density, and with 21 narrowband selected Lyman-$α$ galaxies, among which 16 have been spectroscopically confirmed. At redshifts similar to or above this record, smaller protogroups with fewer members have been reported. LAGER-z7OD1 shows an elongated shape and consists of two subprotoclusters, which would have merged into one massive cluster with a present-day mass of $3.7 \times 10^{15}$ solar masses. The total volume of the ionized bubbles generated by its member galaxies is found to be comparable to the volume of the protocluster itself, indicating that we are witnessing the merging of the individual bubbles and that the intergalactic medium within the protocluster is almost fully ionized. LAGER-z7OD1 thus provides a unique natural laboratory to investigate the reionization process.
△ Less
Submitted 25 January, 2021;
originally announced January 2021.
-
Void Probability Function of Simulated Surveys of high-redshift Lyman-Alpha Emitters
Authors:
Lucia A. Perez,
Sangeeta Malhotra,
James E. Rhoads,
Vithal Tilvi
Abstract:
We calculate the void probability function (VPF) in simulations of Lyman-$α$ emitters (LAEs) across a wide redshift range ($z=3.1,\ 4.5,\ 5.7,\ 6.6$). The VPF measures the zero-point correlation function (i.e. places devoid of galaxies) and naturally connects to higher order correlation functions while being computationally simple to calculate. We explore the Poissonian and systematic errors on th…
▽ More
We calculate the void probability function (VPF) in simulations of Lyman-$α$ emitters (LAEs) across a wide redshift range ($z=3.1,\ 4.5,\ 5.7,\ 6.6$). The VPF measures the zero-point correlation function (i.e. places devoid of galaxies) and naturally connects to higher order correlation functions while being computationally simple to calculate. We explore the Poissonian and systematic errors on the VPF, specify its accuracy as a function of average source density and the volume probed, and provide the appropriate size scales to measure the VPF. At small radii the accuracy of the VPF is limited by galaxy density, while at large radii the VPF is limited by the number of independent volumes probed. We also offer guidelines for understanding and quantifying the error in the VPF. We approximate the error in the VPF by using independent sub-volumes of the catalogs, after finding that jackknife statistics underestimate the uncertainty. We use the VPF to probe the strength of higher order correlation functions by measuring and examining the hierarchical scaling between the correlation functions using count-in-cells. The negative binomial model (NBM) has been shown to best describe the scaling between the two point correlation function and VPF for low-redshift galaxy observations. We further test the fit of the NBM by directly deriving the volume averaged two-point correlation function from the VPF and vice versa. We find the NBM best describes the $z=3.1, 4.5, 5.7$ simulated LAEs, with a 1$σ$ deviation from the model in the $z=6.6$ catalog. This suggests that LAEs show higher order clustering terms similar to those of normal low redshift galaxies.
△ Less
Submitted 6 November, 2020;
originally announced November 2020.
-
A comprehensive study of H$α$ emitters at $z \sim$ 0.62 in the DAWN survey: the need for deep and wide regions
Authors:
Santosh Harish,
Alicia Coughlin,
James E. Rhoads,
Sangeeta Malhotra,
Steven L. Finkelstein,
Matthew Stevans,
Vithal S. Tilvi,
Ali Ahmad Khostovan,
Sylvain Veilleux,
Junxian Wang,
Pascale Hibon,
Johnnes Zabl,
Bhavin Joshi,
John Pharo,
Isak Wold,
Lucia A. Perez,
Zhen-Ya Zheng,
Ronald Probst,
Rob Swaters,
Bahram Mobasher,
Tianxing Jiang,
Huan Yang
Abstract:
We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$α$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$α$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total…
▽ More
We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$α$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$α$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$α$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}α}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$α$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $φ^*=10^{-2.85}$ Mpc$^{-3}$, $α= -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $φ^*=10^{-2.8}$ Mpc$^{-3}$, $α=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $φ^*=10^{-2.91}$ Mpc$^{-3}$, $α= -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$α$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $ρ_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $ρ_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $ρ_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$α$ surveys.
△ Less
Submitted 28 January, 2020;
originally announced January 2020.
-
A large, deep 3 deg$^2$ survey of H$α$, [OIII], and [OII] emitters from LAGER: constraining luminosity functions
Authors:
Ali Ahmad Khostovan,
Sangeeta Malhotra,
James E. Rhoads,
Chunyan Jiang,
Junxian Wang,
Isak Wold,
Zhen-Ya Zheng,
L. Felipe Barrientos,
Alicia Coughlin,
Santosh Harish,
Weida Hu,
Leopoldo Infante,
Lucia A. Perez,
John Pharo,
Francisco Valdes,
Alistair R. Walker,
Huan Yang
Abstract:
We present our measurements of the H$α$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$α$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$σ$ depths of…
▽ More
We present our measurements of the H$α$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$α$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$σ$ depths of $8.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and comoving volumes of $(1-7)\times10^{5}$ Mpc$^3$ making our survey one of the deepest narrowband surveys. We measure the observed luminosity functions and find best-fits of $φ^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$α$, $φ^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc$^{-3}$ and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s$^{-1}$ for [OIII], and $φ^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $α$ fixed to $-1.75$, $-1.6$, and $-1.3$, respectively. An excess of bright $> 10^{42}$ erg s$^{-1}$ [OIII] emitters is observed and may be due to AGN contamination. Dust corrections are applied assuming $A_{\rm{H}α} = 1$ mag. We also design our own empirical rest-frame $g - r$ calibration using SDSS DR12 data, test it against our $z = 0.47$ H$α$ emitters with $z$COSMOS $1$D spectra, and calibrate it for $(g - r)$ between $-0.8$ and $1.3$ mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as $\log_{10}ρ_{\rm{SFR}}/(\rm{M}_\odot\ \rm{yr}^{-1}\ \rm{Mpc}^{-3}) = -1.63\pm0.04$, $-1.07\pm0.06$, and $-0.90\pm0.10$ for H$α$, [OIII], and [OII], respectively. We find our [OIII] and [OII] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multi-wavelength SFRDs, while the H$α$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD.
△ Less
Submitted 14 January, 2020;
originally announced January 2020.
-
The Ly$α$ Luminosity Function and Cosmic Reionization at $z \sim$ 7.0: a Tale of Two LAGER Fields
Authors:
Weida Hu,
Junxian Wang,
Zhen-Ya Zheng,
Sangeeta Malhotra,
James E. Rhoads,
Leopoldo Infante,
L. Felipe Barrientos,
Huan Yang,
Chunyan Jiang,
Wenyong Kang,
Lucia A. Perez,
Isak Wold,
Pascale Hibon,
Linhua Jiang,
Ali Ahmad Khostovan,
Francisco Valdes,
Alistair R. Walker,
Gaspar Galaz,
Alicia Coughlin,
Santosh Harish,
Xu Kong,
John Pharo,
XianZhong Zheng
Abstract:
We present the largest-ever sample of 79 Ly$α$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two l…
▽ More
We present the largest-ever sample of 79 Ly$α$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$α$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Lyα}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$α$ luminosity function.
△ Less
Submitted 2 December, 2019; v1 submitted 21 March, 2019;
originally announced March 2019.
-
Design for the First Narrowband Filter for the Dark Energy Camera: Optimizing the LAGER Survey for z ~ 7 Galaxies
Authors:
Zhen-Ya Zheng,
James E. Rhoads,
Junxian Wang,
Sangeeta Malhotra,
Alistair Walker,
Thomas Mooney,
Chunyan Jiang,
Weida Hu,
Pascale Hibon,
Linhua Jiang,
Leopoldo Infante,
L. Felipe Barrientos,
Gaspar Galaz,
Francisco Valdes,
William Wester,
Huan Yang,
Alicia Coughlin,
Santosh Harish,
Wenyong Kang,
Ali Ahmad Khostovan,
Xu Kong,
Lucia A. Perez,
John Pharo,
Isak Wold,
XianZhong Zheng
Abstract:
We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in t…
▽ More
We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.
△ Less
Submitted 21 March, 2019;
originally announced March 2019.
-
Calibrating the relation of low-frequency radio continuum to star formation rate at 1 kpc scale with LOFAR
Authors:
V. Heesen,
E. Buie II,
C. J. Huff,
L. A. Perez,
J. G. Woolsey,
D. A. Rafferty,
A. Basu,
R. Beck,
E. Brinks,
C. Horellou,
E. Scannapieco,
M. Brüggen,
R. -J. Dettmar,
K. Sendlinger,
B. Nikiel-Wroczyński,
K. T. Chyży,
P. N. Best,
G. H. Heald,
R. Paladino
Abstract:
Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($Σ_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC…
▽ More
Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($Σ_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $Σ_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $Σ_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $μ\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $Σ_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(Σ_{\rm SFR})_{\rm RC} \propto [(Σ_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $σ= 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($α> -0.65$; $I_ν\propto ν^α$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $σ= 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $Σ_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.