Source author record

Bret D. Lehmer

Bret D. Lehmer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

14 published item(s)

preprint2022arXiv

Elevated Hot Gas and High-Mass X-ray Binary Emission in Low Metallicity Galaxies: Implications for Nebular Ionization and Intergalactic Medium Heating in the Early Universe

High-energy emission associated with star formation has been proposed as a significant source of interstellar medium (ISM) ionization in low-metallicity starbursts and an important contributor to the heating of the intergalactic medium (IGM) in the high-redshift ($z > 8$) Universe. Using Chandra observations of a sample of 30 galaxies at $D \approx$~200--450 Mpc that have high specific star-formation rates of 3--9 Gyr$^{-1}$ and metallicities near $Z \approx 0.3 Z_\odot$, we provide new measurements of the average 0.5--8 keV spectral shape and normalization per unit star-formation rate (SFR). We model the sample-combined X-ray spectrum as a combination of hot gas and high-mass X-ray binary (HMXB) populations and constrain their relative contributions. We derive scaling relations of $\log L_{\rm 0.5-8 keV}^{\rm HMXB}$/SFR $= 40.19 \pm 0.06$ and $\log L_{\rm 0.5-2 keV}^{\rm gas}$/SFR $= 39.58^{+0.17}_{-0.28}$; significantly elevated compared to local relations. The HMXB scaling is also somewhat higher than $L_{\rm 0.5-8 keV}^{\rm HMXB}$-SFR-$Z$ relations presented in the literature, potentially due to our galaxies having relatively low HMXB obscuration and young and X-ray luminous stellar populations. The elevation of the hot gas scaling relation is at the level expected for diminished attenuation due to a reduction of metals; however, we cannot conclude that an $L_{\rm 0.5-2 keV}^{\rm gas}$-SFR-$Z$ relation is driven solely by changes in ISM metal content. Finally, we present SFR-scaled spectral models (both emergent and intrinsic) that span the X-ray--to--IR band, providing new benchmarks for studies of the impact of ISM ionization and IGM heating in the early Universe.

preprint2022arXiv

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr

preprint2022arXiv

The Impact of Inclination-dependent Attenuation on Ultraviolet Star Formation Rate Tracers

We examine and quantify how hybrid (e.g., UV+IR) star formation rate (SFR) estimators and the $A_{\rm FUV}$-$β$ relation depend on inclination for disk-dominated galaxies using spectral energy distribution modeling that utilizes the inclination-dependent attenuation curves described in Doore et al. We perform this analysis on a sample of 133 disk-dominated galaxies from the CANDELS fields and 18 disk galaxies from the Spitzer Infrared Nearby Galaxies Survey and Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel samples. We find that both the hybrid SFR estimators and the $A_{\rm FUV}$-$β$ relation present clear dependencies on inclination. To quantify this dependence in the hybrid SFR estimators, we derive an inclination and a far-UV and near-IR color-dependent parametric relation for converting observed UV and IR luminosities into SFRs. For the $A_{\rm FUV}$-$β$ relation, we introduce an inclination-dependent component that accounts for the majority of the inclination dependence with the scatter of the relation increasing with inclination. We then compare both of these inclination-dependent relations to similar inclination-independent relations found in the literature. From this comparison, we find that the UV+IR correction factor and $A_{\rm FUV}$ for our hybrid and $A_{\rm FUV}$-$β$ relations, respectively, result in a reduction in the residual scatter of our sample by approximately a factor of 2. Therefore, we demonstrate that inclination must be considered in hybrid SFR estimators and the $A_{\rm FUV}$-$β$ relation to produce more accurate SFR estimates in disk-dominated galaxies.

preprint2020arXiv

The Metallicity Dependence of the High-Mass X-ray Binary Luminosity Function

We present detailed constraints on the metallicity dependence of the high mass X-ray binary (HMXB) X-ray luminosity function (XLF). We analyze ~5 Ms of Chandra data for 55 actively star-forming galaxies at D < 30 Mpc with gas-phase metallicities spanning 12 + log(O/H) = 7-9.2. Within the galactic footprints, our sample contains a total of 1311 X-ray point sources, of which ~49% are expected to be HMXBs, with the remaining sources likely to be low-mass X-ray binaries (LMXBs; ~22%) and unrelated background sources (~29%). We construct a model that successfully characterizes the average HMXB XLF over the full metallicity range. We demonstrate that the SFR-normalized HMXB XLF shows clear trends with metallicity, with steadily increasing numbers of luminous and ultraluminous X-ray sources (logL(erg/s) = 38-40.5) with declining metallicity. However, we find that the low-luminosity (logL(erg/s) = 36-38) HMXB XLF appears to show a nearly constant SFR scaling and slope with metallicity. Our model provides a revised scaling relation of integrated LX/SFR versus 12 + log(O/H) and a new characterization of its the SFR-dependent stochastic scatter. The general trend of this relation is broadly consistent with past studies based on integrated galaxy emission; however, our model suggests that this relation is driven primarily by the high-luminosity end of the HMXB XLF. Our results have implications for binary population synthesis models, the nature of super-Eddington accreting objects (e.g., ultraluminous X-ray sources), recent efforts to identify active galactic nucleus candidates in dwarf galaxies, and the X-ray radiation fields in the early Universe during the epoch of cosmic heating at z > 10.

preprint2020arXiv

X-ray Binary Luminosity Function Scaling Relations in Elliptical Galaxies: Evidence for Globular Cluster Seeding of Low-Mass X-ray Binaries in Galactic Fields

We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ray-detected sources as being either coincident with unrelated background/foreground objects, GCs, or sources that are within the fields of the galaxy targets. We model the GC and field LMXB population XLFs for all galaxies separately, and then construct global models characterizing how the LMXB XLFs vary with galaxy stellar mass and SN. We find that our field LMXB XLF models require a component that scales with SN, and has a shape consistent with that found for the GC LMXB XLF. We take this to indicate that GCs are "seeding" the galactic field LMXB population, through the ejection of GC-LMXBs and/or the diffusion of the GCs in the galactic fields themselves. However, we also find that an important LMXB XLF component is required for all galaxies that scales with stellar mass, implying that a substantial population of LMXBs are formed "in situ," which dominates the LMXB population emission for galaxies with SN < 2. For the first time, we provide a framework quantifying how directly-associated GC LMXBs, GC-seeded LMXBs, and in-situ LMXBs contribute to LMXB XLFs in the broader early-type galaxy population.

preprint2016arXiv

Demonstrating the likely neutron star nature of five M31 globular cluster sources with Swift-NuSTAR spectroscopy

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient.

preprint2016arXiv

The deepest X-ray view of high-redshift galaxies: constraints on low-rate black-hole accretion

We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7σ$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7σ$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.

preprint2014arXiv

A new candidate Wolf-Rayet X-ray binary in NGC 253

We have discovered a persistent, but highly variable X-ray source in the nearby starburst galaxy NGC 253. The source varies at the level of a factor of about 5 in count rate on timescales of a few hours. Two long observations of the source with Chandra and XMM-Newton show suggestive evidence for the source having a period of about 14-15 hours, but the time sampling in existing data is insufficient to allow a firm determination that the source is periodic. Given the amplitude of variation and the location in a nuclear starburst, the source is likely to be a Wolf-Rayet X-ray binary, with the tentative period being the orbital period of the system. In light of the fact that we have demonstrated that careful examination of the variability of moderately bright X-ray sources in nearby galaxies can turn up candidate Wolf-Rayet X-ray binaries, we discuss the implications of Wolf-Rayet X-ray binaries for predictions of the gravitational wave source event rate, and, potentially, interpretations of the events.

preprint2014arXiv

Evidence for a constant IMF in early-type galaxies based on their X-ray binary populations

A number of recent studies have proposed that the stellar initial mass function (IMF) of early type galaxies varies systematically as a function of galaxy mass, with higher mass galaxies having bottom heavy IMFs. These bottom heavy IMFs have more low-mass stars relative to the number of high mass stars, and therefore naturally result in proportionally fewer neutron stars and black holes. In this paper, we specifically predict the variation in the number of black holes and neutron stars based on the power-law IMF variation required to reproduce the observed mass-to-light ratio trends with galaxy mass. We then test whether such variations are observed by studying the field low-mass X-ray binary populations (LMXBs) of nearby early-type galaxies. In these binaries, a neutron star or black hole accretes matter from a low-mass donor star. Their number is therefore expected to scale with the number of black holes and neutron stars present in a galaxy. We find that the number of LMXBs per K-band light is similar among the galaxies in our sample. These data therefore demonstrate the uniformity of the slope of the IMF from massive stars down to those now dominating the K-band light, and are consistent with an invariant IMF. Our results are inconsistent with an IMF which varies from a Kroupa/Chabrier like IMF for low mass galaxies to a steep power-law IMF (with slope $x$=2.8) for high mass galaxies. We discuss how these observations constrain the possible forms of the IMF variations and how future Chandra observations can enable sharper tests of the IMF.

preprint2014arXiv

Spatially Resolving a Starburst Galaxy at Hard X-ray Energies: NuSTAR, Chandra, AND VLBA Observations of NGC 253

Prior to the launch of NuSTAR, it was not feasible to spatially resolve the hard (E > 10 keV) emission from galaxies beyond the Local Group. The combined NuSTAR dataset, comprised of three ~165 ks observations, allows spatial characterization of the hard X-ray emission in the galaxy NGC 253 for the first time. As a follow up to our initial study of its nuclear region, we present the first results concerning the full galaxy from simultaneous NuSTAR, Chandra, and VLBA monitoring of the local starburst galaxy NGC 253. Above ~10 keV, nearly all the emission is concentrated within 100" of the galactic center, produced almost exclusively by three nuclear sources, an off-nuclear ultraluminous X-ray source (ULX), and a pulsar candidate that we identify for the first time in these observations. We detect 21 distinct sources in energy bands up to 25 keV, mostly consisting of intermediate state black hole X-ray binaries. The global X-ray emission of the galaxy - dominated by the off-nuclear ULX and nuclear sources, which are also likely ULXs - falls steeply (photon index >~ 3) above 10 keV, consistent with other NuSTAR-observed ULXs, and no significant excess above the background is detected at E > 40 keV. We report upper limits on diffuse inverse Compton emission for a range of spatial models. For the most extended morphologies considered, these hard X-ray constraints disfavor a dominant inverse Compton component to explain the γ-ray emission detected with Fermi and H.E.S.S. If NGC 253 is typical of starburst galaxies at higher redshift, their contribution to the E > 10 keV cosmic X-ray background is < 1%.

preprint2013arXiv

Energy Feedback from X-ray Binaries in the Early Universe

X-ray photons, because of their long mean-free paths, can easily escape the galactic environments where they are produced, and interact at long distances with the inter-galactic medium, potentially having a significant contribution to the heating and reionization of the early Universe. The two most important sources of X-ray photons in the Universe are active galactic nuclei (AGN) and X-ray binaries (XRBs). In this Letter we use results from detailed, large scale population synthesis simulations to study the energy feedback of XRBs, from the first galaxies (z~ 20) until today. We estimate that X-ray emission from XRBs dominates over AGN at z>6-8. The shape of the spectral energy distribution of the emission from XRBs shows little change with redshift, in contrast to its normalization which evolves by ~4 orders of magnitude, primarily due to the evolution of the cosmic star-formation rate. However, the metallicity and the mean stellar age of a given XRB population affect significantly its X-ray output. Specifically, the X-ray luminosity from high-mass XRBs per unit of star-formation rate varies an order of magnitude going from solar metallicity to less than 10% solar, and the X-ray luminosity from low-mass XRBs per unit of stellar mass peaks at an age of ~300 Myr and then decreases gradually at later times, showing little variation for mean stellar ages > 3 Gyr. Finally, we provide analytical and tabulated prescriptions for the energy output of XRBs, that can be directly incorporated in cosmological simulations.

preprint2013arXiv

Evidence for Elevated X-ray Emission in Local Lyman Break Galaxy Analogs

In this paper, we study the relationship between the 2-10 keV X-ray luminosity (L_X), assumed to originate from X-ray binaries (XRBs), and star formation rate (SFR) in UV-selected z<0.1 Lyman break analogs (LBAs). We present Chandra observations for four new GALEX-selected LBAs. Including previously studied LBAs, Haro 11 and VV 114, we find that LBAs demonstrate L_X/SFR ratios that are elevated by ~1.5sigma compared to local galaxies, similar to the ratios found for stacked LBGs in the early Universe (z>2). We show that these LBAs are unlikely to harbor AGN, based on their optical and X-ray spectra and the spatial distribution of the X-rays in three spatially extended cases. We expect that high-mass X-ray binaries (HMXBs) dominate the X-ray emission in these galaxies, based on their high specific SFRs (sSFRs=SFR/M* > 10^{-9}/yr), which suggest the prevalence of young stellar populations. Since both LBAs and LBGs have lower dust attenuations and metallicities compared to similar samples of more typical local galaxies, we investigate the effects of dust extinction and metallicity on the L_X/SFR for the broader population of galaxies with high sSFRs (>10^{-10}/yr). The estimated dust extinctions (corresponding to column densities of N_H<10^{22}/cm^2) are expected to have insignificant effects on observed L_X/SFR ratio for the majority of galaxy samples. We find that the observed relationship between L_X/SFR and metallicity appears consistent with theoretical expectations from X-ray binary population synthesis models. Therefore, we conclude that lower metallicities, related to more luminous HMXBs such as ultraluminous X-ray sources (ULXs), drive the elevated L_X/SFR observed in our sample of z<0.1 LBAs. The relatively metal-poor, active mode of star formation in LBAs and distant z>2 LBGs may yield higher total HMXB luminosity than found in typical galaxies in the local Universe.

preprint2012arXiv

The X-ray Star Formation Story as Told by Lyman Break Galaxies in the 4 Ms CDF-S

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

preprint2010arXiv

Towards a Complete Census of AGNs in Nearby Galaxies: The Incidence of Growing Black Holes

We investigate the local supermassive black hole (SMBH) density function and relative mass accretion rates of all active galactic nuclei (AGNs) identified in a volume-limited sample of infrared (IR) bright galaxies (L_IR > 3 x 10^9 L_sun) to D<15 Mpc (Goulding & Alexander 2009). A database of accurate SMBH mass (M_BH) estimates is compiled from literature sources using physically motivated AGN modeling techniques (reverberation mapping, maser mapping and gas kinematics) and well-established indirect M_BH estimation methods (the M-sigma and M_BH-L_(K,bul) relations). For the three sources without previously published M_BH estimates, we use 2MASS K-band imaging and GALFIT to constrain the bulge luminosities, and hence SMBH masses. In general, we find the AGNs in the sample host SMBHs which are spread over a wide mass range (M_BH ~ (0.1-30) x 10^7 M_sun), but with the majority in the poorly studied M_BH ~ 10^6-10^7 M_sun region. Using sensitive hard X-ray (2-10 keV) and mid-IR constraints we calculate the bolometric luminosities of the AGNs (L_(Bol,AGN)) and use them to estimate relative mass accretion rates. We use these data to calculate the volume-average SMBH growth rate of galaxies in the local Universe and find that the AGNs hosting SMBHs in the mass range M_BH ~ 10^6-10^7 M_sun are dominated by optically unidentified AGNs. These relatively small SMBHs are acquiring a significant proportion of their mass in the present-day, and are amongst the most rapidly growing in the local Universe (SMBH mass doubling times of ~6 Gyrs). Additionally, we find tentative evidence for an increasing volume-weighted AGN fraction with decreasing SMBH mass in the M_BH ~ 10^6-10^8 M_sun range. Overall, we conclude that significant mass accretion onto small SMBHs may be missed in even the most sensitive optical surveys due to absent or weak optical AGN signatures.