Researcher profile

Naveen A. Reddy

Naveen A. Reddy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
3topics
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

14 published item(s)

preprint2023arXiv

Excitation and Ionization Properties of Star-forming Galaxies at z=2.0-9.3 with JWST/NIRSpec

We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at $z=2.0-9.3$ from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the [OIII]/H$β$ vs. [NII]/H$α$, [SII]/H$α$, and [OI]/H$α$ ``BPT&#39;&#39; diagrams at $z>2.7$, and the first analysis of the O32 vs. R23 diagram at $z>4$ with a large sample. We divide the sample into five redshift bins containing 30-40 galaxies each. The subsamples at $z\sim2.3$, $z\sim3.3$, and $z\sim4.5$ are representative of the main-sequence star-forming galaxy population at these redshifts, while the $z\sim5.6$ and $z\sim7.5$ samples are likely biased toward high specific star-formation rate due to selection effects. Using composite spectra, we find that each subsample at $z=2.0-6.5$ falls on the same excitation sequence in the [NII] and [SII] BPT diagrams and the O32-R23 diagram on average, offset from the sequences followed by $z=0$ HII regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between $z\sim2$ and $z\sim6$. Overall, the rest-optical line ratios suggest the $z=2.7-9.3$ CEERS/NIRSpec galaxies at log($M_*/M_{\odot})\sim7.5-10$ have high degrees of ionization and moderately low oxygen abundances ($\sim0.1-0.3~Z_{\odot}$), but are not extremely metal poor ($<0.1~Z_{\odot}$) even at $z>6.5$.

preprint2022arXiv

An uncontaminated measurement of the escaping Lyman continuum at $z\sim3$

Observations of reionization-era analogs at $z\sim3$ are a powerful tool for constraining reionization. Rest-ultraviolet observations are particularly useful, in which both direct and indirect tracers of ionizing-photon production and escape can be observed. We analyse a sample of 124 $z\sim3$ galaxies from the Keck Lyman Continuum Spectroscopic Survey, with sensitive spectroscopic measurements of the Lyman continuum (LyC) region. We present a method of removing foreground contamination from our sample using high-resolution, multiband Hubble Space Telescope imaging. We re-measure the global properties of the cleaned sample of 13 individually-detected LyC sources and 107 individually-undetected sources, including a sample-averaged absolute escape fraction of $f_{esc,abs}=0.06\pm0.01$ and a sample-averaged ratio of ionizing to non-ionizing ultraviolet flux density of $<f_{900}/f_{1500}>_{out}=0.040\pm0.006$, corrected for attenuation from the intergalactic and circumgalactic media. Based on composite spectra, we also recover a strong positive correlation between $<f_{900}/f_{1500}>_{out}$ and Ly$α$ equivalent width (W$_λ$(Ly$α$)) and a negative correlation between $<f_{900}/f_{1500}>_{out}$ and UV luminosity. As in previous work, we interpret the relationship between $<f_{900}/f_{1500}>_{out}$ and W$_λ$(Ly$α$) in terms of the modulation of the escape of ionizing radiation from star-forming galaxies based on the covering fraction of neutral gas. We also use a W$_λ$(Ly$α$)-weighted $<f_{900}/f_{1500}>_{out}$ to estimate an ionizing emissivity from star-forming galaxies at $z\sim3$ as $ε_{\rm LyC}\simeq5.5\times10^{24}$erg s$^{-1}$Hz$^{-1}$Mpc$^{-3}$. This estimate, evaluated using the uncontaminated sample of this work, affirms that the contribution of galaxies to the ionizing background at $z\sim3$ is comparable to that of active galactic nuclei.

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.

preprint2022arXiv

Reconciling the Results of the z~2 MOSDEF and KBSS-MOSFIRE Surveys

The combination of the MOSDEF and KBSS-MOSFIRE surveys represents the largest joint investment of Keck/MOSFIRE time to date, with ~3000 galaxies at 1.4<=z<=3.8, roughly half of which are at z~2. MOSDEF is photometric- and spectroscopic-redshift selected with a rest-optical magnitude limit, while KBSS-MOSFIRE is primarily selected based on rest-UV colors and a rest-UV magnitude limit. Analyzing both surveys in a uniform manner with consistent spectral-energy-distribution (SED) models, we find that the MOSDEF z~2 targeted sample has a higher median M_* and redder rest U-V color than the KBSS-MOSFIRE z~2 targeted sample, and a smaller median SED-based SFR and sSFR (SFR(SED) and sSFR(SED)). Specifically, MOSDEF targeted a larger population of red galaxies with U-V and V-J >=1.25, while KBSS-MOSFIRE contains more young galaxies with intense star formation. Despite these differences in the z~2 targeted samples, the subsets of the surveys with multiple emission lines detected and analyzed in previously published work are much more similar. All median host-galaxy properties with the exception of stellar population age -- i.e., M_*, SFR(SED), sSFR(SED), A_V, and UVJ colors -- agree within the uncertainties. Additionally, when uniform emission-line fitting and stellar Balmer absorption correction techniques are applied, there is no significant offset between the two samples in the [OIII]$λ$5008/H$β$ vs. [NII]$λ$6585/H$α$ diagnostic diagram, in contrast to previously-reported discrepancies. We can now combine the MOSDEF and KBSS-MOSFIRE surveys to form the largest z~2 sample with moderate-resolution rest-optical spectra and construct the fundamental scaling relations of star-forming galaxies during this important epoch.

preprint2022arXiv

Searching for the connection between ionizing-photon escape and the surface density of star formation at z~3

The connection between the escape fraction of ionizing photons ($f_{\rm esc}$) and star-formation rate surface density ($Σ_{\rm SFR}$) is a key input for reionization models, but remains untested at high redshift. We analyse 35 z~3 galaxies from the Keck Lyman Continuum Survey (KLCS) covered by deep, rest far-UV spectra of the Lyman continuum (LyC) and high-resolution HST V$_{606}$ imaging, enabling estimates of both $f_{\rm esc}$ and rest-UV sizes. Using Sérsic profile fits to HST images and spectral-energy distribution fits to multi-band photometry, we measure effective sizes and star-formation rates for the galaxies in our sample, and separate the sample into two bins of $Σ_{\rm SFR}$. Based on composite spectra, we estimate <$f_{\rm esc}$> for both $Σ_{\rm SFR}$ subsamples, finding no significant difference in <$f_{\rm esc}$> between the two. To test the representativeness of the KLCS HST sample and the robustness of this result, we attempt to recover the well-established correlation between $f_{\rm esc}$ and Ly$α$ equivalent width. This correlation is not significant within the KLCS HST sample, indicating that the sample is not sufficient for correlating $f_{\rm esc}$ and galaxy properties such as $Σ_{\rm SFR}$. We perform stacking simulations using the KLCS parent sample to determine the optimal sample size for robust probes of the $f_{\rm esc}$-$Σ_{\rm SFR}$ connection to inform future HST LyC observing programs. For a program with a selection independent of ionizing properties, >= 90 objects are required; for one preferentially observing strongly-leaking LyC sources, >= 58 objects are required. More generally, measuring the connection between $f_{\rm esc}$ and $Σ_{\rm SFR}$ requires a larger, representative sample spanning a wide dynamic range in galaxies properties such as $Σ_{\rm SFR}$.

preprint2022arXiv

The MOSDEF Survey: Implications of the Lack of Evolution in the Dust Attenuation-Mass Relation to z~2

We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H$α$ and H$β$ emission lines, and photometric measurements of the rest-UV stellar continuum. The H$α$/H$β$ ratio is used to determine the magnitude of attenuation at the wavelength of H$α$, $A_{{\rm H}α}$. Rest-UV colors and spectral-energy-distribution fitting are used to estimate $A_{1600}$, the magnitude of attenuation at a rest wavelength of 1600Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and $M_*$ over the redshift range $z\sim0$ to $z\sim2.3$. Folding in the latest estimates of the evolution of $M_{\rm dust}$, $({M_{\rm dust}}/{M_{\rm gas}})$, and gas surface density at fixed $M_*$, we find that the expected $M_{\rm dust}$ and dust mass surface density are both significantly higher at $z\sim2.3$ than at $z\sim0$. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation vs. $M_*$, it is essential to determine the relationship between metallicity and $({M_{\rm dust}}/{M_{\rm gas}})$, the dust mass absorption coefficient, and dust geometry, and the evolution of these relations and quantities from $z\sim0$ to $z\sim2.3$.

preprint2021arXiv

The Effects of Stellar Population and Gas Covering Fraction on the Emergent Lyman Alpha Emission of High-Redshift Galaxies

We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85<z<3.49 star-forming galaxies to deduce key properties of the massive stars, ionized ISM, and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Ly-alpha (Lya) photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Ly-beta through CIII] 1907, 1909; and [OII], [NeIII], H-beta, [OIII], H-alpha, [NII], and [SII]. Spectral and photoionization modeling indicate that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Lya equivalent width, W(Lya), and only a marginal anti-correlation between age and W(Lya). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W(Lya). Rather, the covering fraction of optically-thick HI appears to be the principal factor modulating the escape of Lya, with most of the Lya photons in down-the-barrel observations of galaxies escaping through low-column-density or ionized channels in the ISM. Our analysis shows that a high star-formation-rate surface density, Sigma_SFR, particularly when coupled with a low galaxy potential (i.e., low stellar mass), can aid in reducing the covering fraction and ease the escape of Lya photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.

preprint2021arXiv

The MOSDEF Survey: A Comprehensive Analysis of the Rest-optical Emission-line Properties of $z\sim 2.3$ Star-forming Galaxies

We analyze the rest-optical emission-line spectra of $z\sim2.3$ star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O III]5008/H$β$ vs. [N II]6585/H$α$ (&#34;[N II] BPT&#34;) diagram, we define two populations of $z\sim2.3$ MOSDEF galaxies. These include the &#34;high&#34; population that is offset towards higher [O III]5008/H$β$ and/or [N II]6585/H$α$ with respect to the local SDSS sequence and the &#34;low&#34; population that overlaps the SDSS sequence. These two groups are also segregated within the [O III]5008/H$β$ vs. [S II]6718,6733/H$α$ and the [O III]4960,5008/[O II]3727,3730 (O$_{32}$) vs. ([O III]4960,5008+[O II]3727,3730)/H$β$ (R$_{23}$) diagram, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the &#34;high&#34; sample is on average smaller in size and less massive, but has higher specific star-formation rate and star-formation-rate surface density values and is slightly younger compared to the &#34;low&#34; population. From Cloudy+BPASS photoionization models, we estimate that the &#34;high&#34; population has a lower stellar metallicity (i.e., harder ionizing spectrum) but slightly higher nebular metallicity and higher ionization parameter compared to the &#34;low&#34; population. While the &#34;high&#34; population is more $α$-enhanced (i.e., higher $α$/Fe) than the &#34;low&#34; population, both samples are significantly more $α$-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies -- not only those &#34;offset&#34; from local excitation sequences.

preprint2021arXiv

The MOSDEF Survey: Environmental dependence of the gas-phase metallicity of galaxies at $1.4 \leq z \leq 2.6$

Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the role of local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [NII]$\lambda6584$ and H$α$, we measure gas-phase oxygen abundance of 167 galaxies at $1.37\leq z\leq1.7$ and 303 galaxies at $2.09\leq z\leq2.61$, located in diverse environments. We find that at $z\sim1.5$, the average metallicity of galaxies in overdensities with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.8}M_\odot$ is higher relative to their field counterparts by $0.094\pm0.051$, $0.068\pm0.028$ and $0.052\pm0.043$ dex, respectively. However, this metallicity enhancement does not exist at higher redshift, $z\sim2.3$, where, compared to the field galaxies, we find $0.056\pm0.043$, $0.056\pm0.028$ and $0.096\pm 0.034$ dex lower metallicity for galaxies in overdense environments with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.7}M_\odot$, respectively. Our results suggest that, at $1.37\leq z\leq2.61$, the variation of mass-metallicity relation with local environment is small ($<0.1$dex), and reverses at $z\sim2$. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the Universe evolves to lower redshifts ($z\lesssim2$), shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities.

preprint2020arXiv

The MOSDEF Survey: Calibrating the relationship between H$α$ star-formation rate and radio continuum luminosity at $1.4 < z < 2.6$

The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$α$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$α$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$α$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$.

preprint2020arXiv

The MOSDEF Survey: Kinematic and Structural Evolution of Star-Forming Galaxies at $1.4\leq z\leq 3.8$

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt-Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, $V/σ_{V,0}$. We find that $V/σ_{V,0}$ increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ($Σ_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $Σ_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $Σ_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

preprint2020arXiv

The MOSDEF-LRIS Survey: The Interplay Between Massive Stars and Ionized Gas in High-Redshift Star-Forming Galaxies

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

preprint2020arXiv

The Redshift Evolution of Rest-UV Spectroscopic Properties to z~5

We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$α$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$α$ equivalent width, performing corrections to Ly$α$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$α$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$α$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$α$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$α$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$α$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$α$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$α$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$α$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$α$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($ξ_{\rm ion}$).

preprint2019arXiv

The MOSDEF Survey: The Metallicity Dependence of X-ray Binary Populations at $z\sim2$

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.