BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Gergö Popping

Gergö Popping contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.CO

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

An empirical representation of a physical model for the ISM [CII], CO, and [CI] emission at redshift $1\leq z\leq9$

Sub-millimeter emission lines produced by the interstellar medium (ISM) are strong tracers of star formation and are some of the main targets of line intensity mapping (LIM) surveys. In this work we present an empirical multi-line emission model that simultaneously covers the mean, scatter, and correlations of [CII], CO J=1-0 to J=5-4, and [CI] lines in redshift range $1\leq z\leq9$. We assume the galaxy ISM line emission luminosity versus halo mass relations can be described by double power laws with redshift-dependent log normal scatter. The model parameters are then derived by fitting to the state of the art semi-analytic simulation results that have successfully reproduced multiple sub-millimeter line observations at $0\leq z\lesssim6$. We cross check the line emission statistics predicted by the semi-analytic simulation and our empirical model, finding that at $z\geq1$ our model reproduces the simulated line intensities with fractional error less than about 10%. The fractional difference is less than 25% for the power spectra. Grounded on physically-motivated and self-consistent galaxy simulations, this computationally efficient model will be helpful in forecasting ISM emission line statistics for upcoming LIM surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An upper limit on [OIII] 88 $μ$m and 1.2 mm continuum emission from a JWST $z \approx 12-13$ galaxy candidate with ALMA

A number of new $z>11$ galaxy candidates have recently been identified based on public James Webb Space Telescope (JWST) NIRCam observations. Spectroscopic confirmation of these candidates is necessary to robustly measure their redshift and put them in the context of our understanding of the buildup of galaxies in the early Universe. GLASS-z13 is one of these candidates, with a reported photometric redshift $z>11.9$. I present publicly available Atacama Large Millimeter/submillimeter Array (ALMA) band 6 Director's Discretionary Time observations (project 2021.A.00020.S; PI T. Bakx), taken to acquire a spectroscopic redshift for GLASS-z13 by searching for [OIII] line emission in the redshift range $z=11.9-13.5$. No [OIII] emission is detected in integrated spectra extracted within an aperture around GLASS-z13, nor when using an automated line finding algorithm (applying different uv-weighting strategies for the imaging). 1.2 mm continuum emission associated to GLASS-z13 is not detected either. If GLASS-z13 is at z$\approx$12-13, this implies a 3-$σ$ upper limit on the [OIII] and rest-frame $\sim$90 $μ$m continuum emission of $\sim1\times10^8\,\rm{L}_\odot$ and 10.8 $μ$Jy, respectively. The non-detection of [OIII] and continuum emission does not necessarily imply that GLASS-z13 is not at $z\approx12-13$. It can also be explained by a low metallicity ($\sim 0.2\,\rm{Z}_\odot$ or lower) and/or high-density (at least 100 $\rm{cm}^{-3}$) interstellar medium. This work demonstrates the synergy between ALMA and JWST to study the properties of the first galaxies, although JWST/NIRSpec spectroscopy will be necessary to confirm or reject the high photometric-redshift of GLASS-z13.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CO Emission, Molecular Gas, and Metallicity in Main-Sequence Star-Forming Galaxies at $z\sim2.3$

We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ conversion factor ($α_{\text{CO}}$) and find that $α_{\text{CO}}$ inversely correlates with metallicity at $z\sim2$. We derive molecular gas masses ($M_{\text{mol}}$) and characterize the relations among $M_*$, SFR, $M_{\text{mol}}$, and metallicity. At $z\sim2$, $M_{\text{mol}}$ increases and molecular gas fraction ($M_{\text{mol}}$/$M_*$) decrease with increasing $M_*$, with a significant secondary dependence on SFR. Galaxies at $z\sim2$ lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR-$M_*$, O/H-$M_*$, and $M_{\text{mol}}$-$M_*$ relations is correlated such that, at fixed $M_*$, $z\sim2$ galaxies with larger $M_{\text{mol}}$ have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at $z\sim2$ where O/H is inversely correlated with both SFR and $M_{\text{mol}}$ at fixed $M_*$. These results suggest that the scatter of the $z\sim2$ star-forming main sequence, mass-metallicity relation, and $M_{\text{mol}}$-$M_*$ relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive $z\sim2$ star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

IR SED and Dust Masses of Sub-solar Metallicity Galaxies at z~2.3

We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $χ^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $χ^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observed cosmic evolution of galaxy dust properties with metallicity and tensions with models

The dust abundance of the interstellar medium plays an important role in galaxy physics, the chemical evolution of matter and the absorption and re-emission of stellar light. The last years have seen a surge in observational and theoretical studies constraining the dust-abundance of galaxies up to $z\sim5$. In this work we gather the latest observational measurements (with a focus on absorption studies covering metallicities in the range $6.8 < 12 + \log{(O/H)}<9$) and theoretical predictions (from six different galaxy formation models) for the dust-to-gas (DTG) and dust-to-metal (DTM) ratio of galaxies. The observed trend between DTG and DTM and gas-phase metallicity can be described by a linear relation and shows no evolution from $0<z<5$. Importantly, the fit to the DTG-metallicity relation provides a refined tool for robust dust-based gas mass estimates inferred from millimeter dust-continuum observations. The lack of evolution in the observed relations are indicative of a quickly reached balance (already when the Universe was 1.2 Gyr old) between the formation and destruction of dust and a constant timescale for star-formation at fixed metallicities over cosmic time. None of the models is able to reproduce the observed trends over the entire range in metallicity and redshift probed. The comparison between models and simulations furthermore rules out some of the current implementations for the growth and destruction of dust in galaxy formation models and places tight constraints on the predicted timescale for star-formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Semi-analytic forecasts for JWST -- VI. Simulated lightcones and galaxy clustering predictions

In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new set of 8 ultra-deep lightcones with 132 sq. arcmin footprints, containing galaxies up to z ~ 12 down to the magnitudes expected to be reached in the deepest JWST surveys. These mock lightcones are extracted from dissipationless N-body simulations and populated with galaxies using the well-established, computationally efficient Santa Cruz semi-analytic model for galaxy formation. We provide a wide range of predicted physical properties, and simulated photometry from NIRCam and many other instruments. We explore the predicted counts and luminosity functions and angular two-point correlation functions for galaxies in these simulated lightcones. We also explore the predicted field-to-field variance using multiple lightcone realizations. We find that these lightcones reproduce the available measurements of observed clustering from 0.2 < z < 7.5 very well. We provide predictions for galaxy clustering at high redshift that may be obtained from future JWST observations. All of the lightcones presented here are made available through a web-based, interactive data release portal https://flathub.flatironinstitute.org/group/sam-forecasts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High-Redshift

The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$β$ and [N II]/H$α$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The UV 2175Å Attenuation Bump and its Correlation with PAH Emission at z~2

The UV bump is a broad absorption feature centered at 2175Å that is seen in the attenuation/extinction curve of some galaxies, but its origin is not well known. Here, we use a sample of 86 star-forming galaxies at z=1.7-2.7 with deep rest-frame UV spectroscopy from the MUSE HUDF Survey to study the connection between the strength of the observed UV 2175Å bump and the Spitzer/MIPS 24 micron photometry, which at the redshift range of our sample probes mid-IR polycyclic aromatic hydrocarbon (PAH) emission at ~6-8 micron. The sample has robust spectroscopic redshifts and consists of typical main-sequence galaxies with a wide range in stellar mass (log(Mstar/Msun) ~ 8.5-10.7) and star formation rates (SFRs; SFR ~ 1-100 Msun/yr). Galaxies with MIPS detections have strong UV bumps, except for those with mass-weighted ages younger than ~150 Myr. We find that the UV bump amplitude does not change with SFR at fixed stellar mass but increases with mass at fixed SFR. The UV bump amplitude and the PAH strength (defined as mid-IR emission normalized by SFR) are highly correlated and both also correlate strongly with stellar mass. We interpret these correlations as the result of the mass-metallicity relationship, such that at low metallicities PAH emission is weak due to a lower abundance of PAH molecules. The weak or complete absence of the 2175Å bump feature on top of the underlying smooth attenuation curve at low mass/metallicities is then expected if the PAH carriers are the main source of the additional UV absorption.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tracing Molecular Gas in z $\simeq$ 6 Galaxies with [C${\rm \scriptsize II}$]

We investigate the fine-structure [C${\rm \scriptsize II}$] line at $158\,μ$m as a molecular gas tracer by analyzing the relationship between molecular gas mass ($M_{\rm mol}$) and [C${\rm \scriptsize II}$] line luminosity ($L_{\rm [CII]}$) in 11,125 $z\simeq 6$ star-forming, main sequence galaxies from the SIMBA simulations, with line emission modeled by SÍGAME. Though most ($\sim 50-100\,\%$) of the gas mass in our simulations is ionized, the bulk ($> 50\,\%$) of the [C${\rm \scriptsize II}$] emission comes from the molecular phase. We find a sub-linear (slope $0.78\pm 0.01$) $\log L_{\rm [CII]}-\log M_{\rm mol}$ relation, in contrast with the linear relation derived from observational samples of more massive, metal-rich galaxies at $z \lesssim 6$. We derive a median [C${\rm \scriptsize II}$]-to-$M_{\rm mol}$ conversion factor of $α_{\rm [CII]} \simeq 18\,{\rm M_{\rm \odot}/L_{\rm \odot}}$. This is lower than the average value of $\simeq 30\,{\rm M_{\rm \odot}/L_{\rm \odot}}$ derived from observations, which we attribute to lower gas-phase metallicities in our simulations. Thus, a lower, luminosity-dependent, conversion factor must be applied when inferring molecular gas masses from [C${\rm \scriptsize II}$] observations of low-mass galaxies. For our simulations, [C${\rm \scriptsize II}$] is a better tracer of the molecular gas than CO $J=1-0$, especially at the lowest metallicities, where much of the gas is &#39;CO-dark&#39;. We find that $L_{\rm [CII]}$ is more tightly correlated with $M_{\rm mol}$ than with star-formation rate (${\rm SFR}$), and both the $\log L_{\rm [CII]}-\log M_{\rm mol}$ and $\log L_{\rm [CII]}-\log {\rm SFR}$ relations arise from the Kennicutt-Schmidt relation. Our findings suggest that $L_{\rm [CII]}$ is a promising tracer of the molecular gas at the earliest cosmic epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources&#39; spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe&#39;s first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe&#39;s first dusty galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On estimating the cosmic molecular gas density from CO Line Intensity Mapping observations

The Millimeter-wave Intensity Mapping Experiment (mmIME) recently reported a detection of excess spatial fluctuations at a wavelength of 3 mm, which can be attributed to unresolved emission of several CO rotational transitions between $z\sim1-5$. We study the implications of this data for the high-redshift interstellar medium using a suite of state-of-the-art semianalytic simulations which have successfully reproduced many other sub-millimeter line observations across the relevant redshift range. We find that the semianalytic predictions are mildly in tension with the mmIME result, with a predicted CO power $\sim3.5σ$ below what was observed. We explore some simple modifications to the models which could resolve this tension. Increasing the molecular gas abundance at the relevant redshifts to $\sim10^8\ M_\odot\ \rm{Mpc}^{-3}$, a value well above that obtained from directly imaged sources, would resolve the discrepancy, as would assuming a CO-$H_2$ conversion factor $α_{\rm{CO}}$ of $\sim1.5\ M_{\odot}$ K$^{-1}$ $(\rm{km}/\rm{s})^{-1}$ pc$^{2}$, a value somewhat lower than is commonly assumed. We go on to demonstrate that these conclusions are quite sensitive to the detailed assumptions of our simulations, highlighting the need for more careful modeling efforts as more intensity mapping data become available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A redshift-dependent IRX-$β$ dust attenuation relation for TNG50 galaxies

We study the relation between the UV-slope, $β$, and the ratio between the infrared- and UV-luminosities (IRX) of galaxies from TNG50, the latest installment of the IllustrisTNG galaxy formation simulations. We select 7280 star-forming main-sequence (SFMS) galaxies with stellar mass $10^9M_\odot$ at redshifts $0 \leq z \leq 4$ and perform radiative transfer with SKIRT to model effects of interstellar medium dust on the emitted stellar light. Assuming a Milky Way (MW) dust type and a dust-to-metal ratio of 0.3, we find that TNG50 SFMS galaxies generally agree with observationally-derived IRX-$β$ relations at $z \lesssim 1$. However, we find a redshift-dependent systematic offset with respect to empirically-derived local relations, with the TNG50 IRX-$β$ relation shifting towards lower $β$ and steppening at higher redshifts. This is partially driven by variations in the dust-uncorrected UV-slope of galaxies, due to different star-formation histories of galaxies selected at different cosmic epochs; we suggest the remainder of the effect is caused by differences in the effective dust attenuation curves (EDACs) of galaxies as a function of redshift. We find a typical galaxy-to-galaxy variation of 0.3 dex in IRX at fixed $β$, correlated with intrinsic galaxy properties: galaxies with higher star-formation rates, star-formation efficiencies, gas metallicities and stellar masses exhibit larger IRX values. We demonstrate a degeneracy between stellar age, dust geometry and dust composition: $z=4$ galaxies with a Small Magellanic Cloud dust type follow the same IRX-$β$ relation as low-redshift galaxies with MW dust. We provide a redshift-dependent fitting function for the IRX-$β$ relation for MW dust based on our models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-analytic forecasts for JWST -- III. Intrinsic production efficiency of Lyman-continuum radiation

The James Webb Space Telescope (JWST) is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the Epoch of Reionization (EoR). A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, $ξ_{\rm ion}$, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of $ξ_{\rm ion}$ on galaxy properties and cosmic epoch from $z=4-10$. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_{\rm ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of two more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of $ξ_{\rm ion}$ about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from $z\sim 12-4$ at fixed rest-UV luminosity. We provide predictions of $\dot{N}_{\rm ion}$ as a function of $M_{\rm h}$ and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-analytic forecasts for JWST -- IV. Implications for cosmic reionization and LyC escape fraction

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level&#39; galaxy formation physics to `top-level&#39; cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: A model to explain observed 1.1 and 0.85 millimeter dust continuum number counts

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $μ$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $μ$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $μ$m number counts is almost exhausted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rmμJy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_ν\sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$. At all redshifts, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The ALMA Spectroscopic Survey in the HUDF: CO emission lines and 3 mm continuum sources

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $β=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The ALMA Spectroscopic Survey in the HUDF: Constraining cumulative CO emission at $1 \lesssim z \lesssim 4$ with power spectrum analysis of ASPECS LP data from 84 to 115 GHz

We present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their second moments in the three-dimensional (3D) auto-power spectrum, finding a total CO shot noise power $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}}) \leq 1.9\times10^2$ $μ$K$^2$ (Mpc h$^{-1}$)$^3$. This upper limit ($3σ$) is consistent with the observed ASPECS CO LFs in Decarli et al. 2019, but rules out a large space in the range of $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}})$ inferred from these LFs, which we attribute primarily to large uncertainties in the normalization $Φ_*$ and knee $L_*$ of the Schechter-form CO LFs at $z > 2$. Also, through power spectrum analyses of ASPECS LP data with 415 positions from galaxies with available optical spectroscopic redshifts, we find that contributions to the observed mean CO intensity and shot noise power of MUSE galaxies are largely accounted for by ASPECS blind detections, though there are $\sim20$% contributions to the CO(2-1) mean intensity due to sources previously undetected in the blind line search. Finally, we sum the fluxes from individual blind CO detections to yield a lower limit on the mean CO surface brightness at 99 GHz of $\langle T_{\textrm{CO}} \rangle = 0.55\pm0.02$ $μ$K, which we estimate represents $68$-$80$% of the total CO surface brightness at this frequency.

0 citations0 reviewsNo code linkedOpen access