Showing 1–50 of 119 results for author: Lehmer, B

Fast SMBH growth in the SPT2349--56 protocluster at $z=4.3$

Authors: Fabio Vito, William N. Brandt, Andrea Comastri, Roberto Gilli, Rob J. Ivison, Giorgio Lanzuisi, Bret D. Lehmer, Ivan E. Lopez, Paolo Tozzi, Cristian Vignali

Abstract: Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via… ▽ More Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via Chandra (200 ks) observations, and comparing the results with the field environment. We detected two protocluster members, corresponding to an AGN fraction among DSFGs of $\approx10\%$. This value is consistent with other protoclusters at $z=2-4$, but higher than the AGN incidence among DSFGs in the field environment. Both AGN are heavily obscured sources and hosted in star-forming galaxies with $\approx3\times10^{10}\,\mathrm{M_\odot}$ stellar masses. We estimate that the ISM in the host galaxies can contribute significantly to the nuclear obscuration. One of the two AGN is highly luminous ($L_X=2\times10^{45}\,\mathrm{erg\,s^{-1}}$) and Compton-thick ($N_H=2\times10^{24}\,\mathrm{cm^{-2}}$), and likely powered by a $M_{BH}>6\times10^8\,\mathrm{M_\odot}$ SMBH. Its high accretion rate suggests that it is in the phase of efficient growth required to explain the presence of extremely massive SMBHs in the centers of local galaxy clusters. Considering SPT2349-56 and DRC, a similar protocuster at $z=4$, we find that gas-rich protocluster cores at $z\approx4$ enhance the triggering of luminous (log$\frac{L_X}{\mathrm{erg\,s^{-1}}}=45-46$) AGN by 3-5 orders of magnitude with respect to the field environment. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures, and thus represent key science targets to obtain a complete understanding of the relation between environment and galaxy evolution. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Submitted to A&A, referee comments included

The HMXB Luminosity Functions of Dwarf Galaxies

Authors: Robel Geda, Andy D. Goulding, Bret D. Lehmer, Jenny E. Greene, Anish Kulkarni

Abstract: Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M* < 5 x 10^9 M_sun), expanding the specific star-formation rat… ▽ More Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M* < 5 x 10^9 M_sun), expanding the specific star-formation rates (sSFR) and gas-phase metallicities probed in previous investigations. Our analysis of the observed X-ray luminosity function indicates a shallower power-law slope for the dwarf galaxy HMXB population. In our study, we focus on dwarf galaxies that are more representative in terms of sSFR compared to prior work. In this regime, the HMXB luminosity function exhibits significant stochastic sampling at high luminosities. This likely accounts for the pronounced scatter observed in the galaxy-integrated HMXB population's Lx/SFR versus metallicity for our galaxy sample. Our calibration is necessary to understand the AGN content of low mass galaxies identified in current and future X-ray survey fields and has implications for binary population synthesis models, as well as X-ray driven cosmic heating in the early universe. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 19 pages, 7 figures

Monitoring the X-ray Variability of Bright X-ray Sources in M33

Authors: Rebecca Kyer, Shelby Albrecht, Benjamin F. Williams, Kyros Hinton, Breanna Binder, Margaret Lazzarini, Kristen Garofali, Bret Lehmer, Michael Eracleous, Paul P. Plucinsky, Vallia Antoniou

Abstract: We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected i… ▽ More We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected in previous X-ray catalogs of M33 and discuss their possible nature. The final catalog includes 26 known HMXB candidates identified in the literature. We extend the baseline of the X-ray light curves up to 21 years by including archival X-ray observations of these sources. We compare the detection and non-detection epochs of the sources to suites of simulated source duty cycles and infer that most of our detected sources have duty cycles > 30%. We find only four sources whose detection patterns are consistent with having duty cycles below 30%. This large fraction of sources with high duty cycles is unexpected for a population of HMXBs, thus more frequent X-ray monitoring will likely reveal many more low duty cycle HMXBs in M33. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: Accepted for publication in ApJ. 31 pages, 14 figures

The evolution of galaxies and clusters at high spatial resolution with AXIS

Authors: H. R. Russell, L. A. Lopez, S. W. Allen, G. Chartas, P. P. Choudhury, R. A. Dupke, A. C. Fabian, A. M. Flores, K. Garofali, E. Hodges-Kluck, M. J. Koss, L. Lanz, B. D. Lehmer, J. -T. Li, W. P. Maksym, A. B. Mantz, M. McDonald, E. D. Miller, R. F. Mushotzky, Y. Qiu, C. S. Reynolds, F. Tombesi, P. Tozzi, A. Trindade-Falcao, S. A. Walker , et al. (3 additional authors not shown)

Abstract: Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks and producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down sta… ▽ More Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks and producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down star formation and ultimately curtailing the growth of structure after the peak at redshift 2-3. To understand the complex interplay between gravity and feedback, we must resolve both the key physics within galaxies and map the impact of these processes over large scales, out into the cosmic web. The Advanced X-ray Imaging Satellite (AXIS) is a proposed X-ray probe mission for the 2030s with arcsecond spatial resolution, large effective area, and low background. AXIS will untangle the interactions of winds, radiation, jets, and supernovae with the surrounding ISM across the wide range of mass scales and large volumes driving galaxy evolution and trace the establishment of feedback back to the main event at cosmic noon. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 29 pages, 18 figures; this white paper is part of a series commissioned for the AXIS Probe mission concept

The High Energy X-ray Probe (HEX-P): resolving the nature of Sgr A* flares, compact object binaries and diffuse X-ray emission in the Galactic Center and beyond

Authors: Kaya Mori, Gabriele Ponti, Matteo Bachetti, Arash Bodaghee, Jonathan Grindlay, Jaesub Hong, Roman Krivonos, Ekaterina Kuznetsova, Shifra Mandel, Antonio Rodriguez, Giovanni Stel, Shuo Zhang, Tong Bao, Franz Bauer, Maica Clavel, Benjamin Coughenour, Javier A. Garcia, Julian Gerber, Brian Grefenstette, Amruta Jaodand, Bret Lehmer, Kristin Madsen, Melania Nynka, Peter Predehl, Ciro Salcedo , et al. (2 additional authors not shown)

Abstract: HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities' (including XMM-Newton and NuSTAR). These capabilities will enable revolutionary new insights into a variety of important astrophysical problems. We present scientific objectives and simulatio… ▽ More HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities' (including XMM-Newton and NuSTAR). These capabilities will enable revolutionary new insights into a variety of important astrophysical problems. We present scientific objectives and simulations of HEX-P observations of the Galactic Center (GC) and Bulge. We demonstrate the unique and powerful capabilities of the HEX-P observatory for studying both X-ray point sources and diffuse X-ray emission. HEX-P will be uniquely equipped to explore a variety of major topics in Galactic astrophysics, allowing us to (1) investigate broad-band properties of X-ray flares emitted from the supermassive black hole (BH) at Sgr A* and probe the associated particle acceleration and emission mechanisms; (2) identify hard X-ray sources detected by NuSTAR and determine X-ray point source populations in different regions and luminosity ranges; (3) determine the distribution of compact object binaries in the nuclear star cluster and the composition of the Galactic Ridge X-ray emission; (4) identify X-ray transients and measure fundamental parameters such as BH spin; (5) find hidden pulsars in the GC; (6) search for BH-OB binaries and hard X-ray flares from young stellar objects in young massive clusters; (7) measure white dwarf (WD) masses of magnetic CVs to deepen our understanding of CV evolution and the origin of WD magnetic fields; (8) explore primary particle accelerators in the GC in synergy with future TeV and neutrino observatories; (9) map out cosmic-ray distributions by observing non-thermal X-ray filaments; (10) explore past X-ray outbursts from Sgr A* through X-ray reflection components from giant molecular clouds. △ Less

Submitted 8 November, 2023; originally announced November 2023.

The High Energy X-ray Probe: Resolved X-ray Populations in Extragalactic Environments

Authors: Bret D. Lehmer, Kristen Garofali, Breanna A. Binder, Francesca Fornasini, Neven Vulic, Andreas Zezas, Ann Hornschemeier, Margaret Lazzarini, Hannah Moon, Toni Venters, Daniel Wik, Mihoko Yukita, Matteo Bachetti, Javier A. García, Brian Grefenstette, Kristin Madsen, Kaya Mori, Daniel Stern

Abstract: We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging ($<$20 arcsec FWHM), broad spectral coverage (0.2-80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will en… ▽ More We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging ($<$20 arcsec FWHM), broad spectral coverage (0.2-80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will enable HEX-P to detect hard (4-25 keV) X-ray emission from resolved point-source populations within $\sim$800 galaxies and integrated emission from $\sim$6000 galaxies out to 100 Mpc. These galaxies cover wide ranges of galaxy types (e.g., normal, starburst, and passive galaxies) and properties (e.g., metallicities and star-formation histories). In such galaxies, HEX-P will: (1) provide unique information about X-ray binary populations, including accretor demographics (black hole and neutron stars), distributions of accretion states and state transition cadences; (2) place order-of-magnitude more stringent constraints on inverse Compton emission associated with particle acceleration in starburst environments; and (3) put into clear context the contributions from X-ray emitting populations to both ionizing the surrounding interstellar medium in low-metallicity galaxies and heating the intergalactic medium in the $z > 8$ Universe. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): Studying Extreme Accretion with Ultraluminous X-ray Sources

Authors: Matteo Bachetti, Matthew J. Middleton, Ciro Pinto, Andrés Gúrpide, Dominic J. Walton, Murray Brightman, Bret Lehmer, Timothy P. Roberts, Georgios Vasilopoulos, Jason Alford, Roberta Amato, Elena Ambrosi, Lixin Dai, Hannah P. Earnshaw, Hamza El Byad, Javier A. García, Gian Luca Israel, Amruta Jaodand, Kristin Madsen, Chandreyee Maitra, Shifra Mandel, Kaya Mori, Fabio Pintore, Ken Ohsuga, Maura Pilia , et al. (3 additional authors not shown)

Abstract: Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the acc… ▽ More Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion -- hence their progenitors, lifetimes, and future evolution -- is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational wave sources. Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities. HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR. We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies. △ Less

Submitted 10 November, 2023; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 24 pages, 9 figures. Accepted for publication in Frontiers of Astronomy and Space Science. Minor corrections included

Modeling the High-Energy Ionizing Output from Simple Stellar and X-ray Binary Populations

Authors: Kristen Garofali, Antara R. Basu-Zych, Benjamin D. Johnson, Panayiotis Tzanavaris, Anne Jaskot, Chris T. Richardson, Bret D. Lehmer, Mihoko Yukita, Edmund Hodges-Kluck, Ann Hornschemeier, Andrew Ptak, Neven Vulic

Abstract: We present a methodology for modeling the joint ionizing impact due to a "simple X-ray population" (SXP) and its corresponding simple stellar population (SSP), where "simple" refers to a single age and metallicity population. We construct composite spectral energy distributions (SEDs) including contributions from ultra-luminous X-ray sources (ULXs) and stars, with physically meaningful and consist… ▽ More We present a methodology for modeling the joint ionizing impact due to a "simple X-ray population" (SXP) and its corresponding simple stellar population (SSP), where "simple" refers to a single age and metallicity population. We construct composite spectral energy distributions (SEDs) including contributions from ultra-luminous X-ray sources (ULXs) and stars, with physically meaningful and consistent consideration of the relative contributions of each component as a function of instantaneous burst age and stellar metallicity. These composite SEDs are used as input for photoionization modeling with Cloudy, from which we produce a grid for the time- and metallicity-dependent nebular emission from these composite populations. We make the results from the photoionization simulations publicly available. We find that the addition of the SXP prolongs the high-energy ionizing output from the population, and correspondingly increases the intensity of nebular lines such as He II $λ$1640,4686, [Ne V] $λ$3426,14.3$μ$m, and [O IV] 25.9$μ$m by factors of at least two relative to models without an SXP spectral component. This effect is most pronounced for instantaneous bursts of star formation on timescales $>$ 10 Myr and at low metallicities ($\sim$ 0.1 $Z_{\odot}$), due to the imposed time- and metallicity-dependent behavior of the SXP relative to the SSP. We propose nebular emission line diagnostics accessible with JWST suitable for inferring the presence of a composite SXP + SSP, and discuss how the ionization signatures compare to models for sources such as intermediate mass black holes. △ Less

Submitted 7 November, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: accepted for publication in ApJ; simulation results available through FSPS (see Section 5.4)

The Large Deficit of HMXB Emission from Luminous Infrared Galaxies: the Case of the Circumnuclear Starburst Ring in NGC 7552

Authors: Lacey West, Kristen Garofali, Bret D. Lehmer, Andrea Prestwich, Rafael Eufrasio, Wasutep Luangtip, Timothy P. Roberts, Andreas Zezas

Abstract: Luminous infrared galaxies (LIRGs), the most extreme star-forming galaxies in the nearby (D$<$30 Mpc) Universe, show a notable X-ray emission deficiency (up to a factor of $\sim$10) compared with predictions from scaling relations of galaxy-wide high mass X-ray binary (HMXB) luminosity with star-formation rate. In the nearby ($\approx$20 Mpc) LIRG NGC 7552, the majority of the IR emission originat… ▽ More Luminous infrared galaxies (LIRGs), the most extreme star-forming galaxies in the nearby (D$<$30 Mpc) Universe, show a notable X-ray emission deficiency (up to a factor of $\sim$10) compared with predictions from scaling relations of galaxy-wide high mass X-ray binary (HMXB) luminosity with star-formation rate. In the nearby ($\approx$20 Mpc) LIRG NGC 7552, the majority of the IR emission originates in a circumnuclear starburst ring, which has been resolved into several discrete knots of star formation. We present results from recent Chandra observations of NGC 7552, which reveal significant deficits in the 2-7 keV X-ray luminosities from two of the most powerful star-forming knots. We hypothesize that the expected luminous HMXB populations in these knots are either (1) obscured by very large column densities or (2) suppressed due to the knots having relatively high metallicity and/or very young ages ($\lesssim$ 5 Myr). We distinguish between these possibilities using data from recent NuSTAR observations, whose sensitivity above 10 keV is capable of uncovering heavily obscured HMXB populations, since emission at these energies is more immune to absorption effects. We find no evidence of a heavily obscured HMXB population in the central region of NGC 7552, suggesting suppressed HMXB formation. We further show that metallicity-dependent scaling relations cannot fully account for the observed deficit from the most powerful star-forming knots or the central region as a whole. Thus, we suggest that recent bursts in local star formation activity likely drive the high $L_{\rm{IR}}$ within these regions on timescales $\lesssim$ 5 Myr, shorter than the timescale required for the formation of HMXBs. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ (17 pages, 3 tables, 8 figures)

Revisiting the Properties of X-ray AGN in the SSA22 Protocluster: Normal SMBH and Host-Galaxy Growth for AGN in a $z=3.09$ Overdensity

Authors: Erik B. Monson, Keith Doore, Rafael T. Eufrasio, Bret D. Lehmer, David M. Alexander, Chris M. Harrison, Mariko Kubo, Cristian Saez, Hideki Umehata

Abstract: We analyze the physical properties of 8 X-ray selected active galactic nuclei (AGN) and one candidate protoquasar system (ADF22A1) in the $z = 3.09$ SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED fitting code, Lightning. We recover star formation histories (SFH) for 7 of these systems which are well-fit by composite stellar population plus AGN mo… ▽ More We analyze the physical properties of 8 X-ray selected active galactic nuclei (AGN) and one candidate protoquasar system (ADF22A1) in the $z = 3.09$ SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED fitting code, Lightning. We recover star formation histories (SFH) for 7 of these systems which are well-fit by composite stellar population plus AGN models. We find indications that 4/9 of the SSA22 AGN systems we study have host galaxies below the main sequence, with $\rm SFR/SFR_{MS} \leq -0.4$. The remaining SSA22 systems, including ADF22A1, are consistent with obscured supermassive black hole (SMBH) growth in star forming galaxies. We estimate the SMBH accretion rates and masses, and compare the properties and SFH of the 9 protocluster AGN systems with X-ray detected AGN candidates in the Chandra Deep Fields (CDF), finding that the distributions of SMBH growth rates, star formation rates, SMBH masses, and stellar masses for the protocluster AGN are consistent with field AGN. We constrain the ratio between the sample-averaged SSA22 SMBH mass and CDF SMBH mass to $<1.41$. While the AGN are located near the density peaks of the protocluster, we find no statistically significant trends between the AGN or host galaxy properties and their location in the protocluster. We interpret the similarity of the protocluster and field AGN populations together with existing results as suggesting that the protocluster and field AGN co-evolve with their hosts in the same ways, while AGN-triggering events are more likely in the protocluster. △ Less

Submitted 11 July, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: 27 pages, 14 figures, 8 tables. Updated with published version

Journal ref: ApJ 951 (2023), 15

X-ray Binaries in External Galaxies

Authors: Marat Gilfanov, Giuseppina Fabbiano, Bret Lehmer, Andreas Zezas

Abstract: X-ray appearance of normal galaxies is mainly determined by X-ray binaries powered by accretion onto a neutron star or a stellar mass black hole. Their populations scale with the star-formation rate and stellar mass of the host galaxy and their X-ray luminosity distributions show a significant split between star-forming and passive galaxies, both facts being consequences of the dichotomy between h… ▽ More X-ray appearance of normal galaxies is mainly determined by X-ray binaries powered by accretion onto a neutron star or a stellar mass black hole. Their populations scale with the star-formation rate and stellar mass of the host galaxy and their X-ray luminosity distributions show a significant split between star-forming and passive galaxies, both facts being consequences of the dichotomy between high- and low-mass X-ray binaries. Metallicity, IMF and stellar age dependencies, and dynamical formation channels add complexity to this picture. The numbers of high-mass X-ray binaries observed in star-forming galaxies indicate quite high probability for a massive star to become an accretion powered X-ray source once upon its lifetime. This explains the unexpectedly high contribution of X-ray binaries to the Cosmic X-ray Background, of the order of $\sim 10\%$, mostly via X-ray emission of faint star-forming galaxies located at moderate redshifts which may account for the unresolved part of the CXB. Cosmological evolution of the $L_X-{\rm SFR}$ relation can make high-mass X-ray binaries a potentially significant factor in (pre)heating of intergalactic medium in the early Universe. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Invited chapter for the Handbook of X-ray and Gamma-ray Astrophysics. Editors: Cosimo Bambi, Andrea Santangelo. Publisher: Springer Singapore, 2023

Lightning: An X-ray to Submillimeter Galaxy SED-Fitting Code With Physically-Motivated Stellar, Dust, and AGN Models

Authors: Keith Doore, Erik B. Monson, Rafael T. Eufrasio, Bret D. Lehmer, Kristen Garofali, Antara Basu-Zych

Abstract: We present an updated version of Lightning, a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGN). X-ray emission, when utilized, can be modeled as originating from stellar compact binar… ▽ More We present an updated version of Lightning, a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGN). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGN. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include an adaptive Markov-Chain Monte-Carlo (MCMC), affine-invariant MCMC, and Levenberg-Marquardt gradient decent (MPFIT) algorithms. To demonstrate some of the capabilities of Lightning, we present several examples using a variety of observational data. These examples include (1) deriving the spatially resolved stellar properties of the nearby galaxy M81, (2) demonstrating how X-ray emission can provide constrains on the properties of the supermassive black hole of a distant AGN, (3) exploring how to rectify the attenuation effects of inclination on the derived the star formation rate of the edge-on galaxy NGC 4631, (4) comparing the performance of Lightning to similar Bayesian SED fitting codes when deriving physical properties of the star-forming galaxy NGC 628, and (5) comparing the derived X-ray and UV-to-IR AGN properties from Lightning and CIGALE for a distant AGN. Lightning is an open-source application developed in the Interactive Data Language (IDL) and is available at https://github.com/rafaeleufrasio/lightning. △ Less

Submitted 16 July, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 34 pages, 17 figures. Updated with published version

Journal ref: ApJS 266 (2023) 39

X-ray luminosity function of high-mass X-ray binaries: Studying the signatures of different physical processes using detailed binary evolution calculations

Authors: Devina Misra, Konstantinos Kovlakas, Tassos Fragos, Margaret Lazzarini, Simone S. Bavera, Bret D. Lehmer, Andreas Zezas, Emmanouil Zapartas, Zepei Xing, Jeff J. Andrews, Aaron Dotter, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun

Abstract: The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of variou… ▽ More The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of various, poorly-constrained assumptions regarding physical processes such as the common-envelope phase, the core-collapse, and wind-fed accretion. We use the new binary population synthesis code POSYDON, which employs extensive pre-computed grids of detailed stellar structure and binary evolution models, to simulate the evolution of binaries. We generate 96 synthetic XRB populations corresponding to different combinations of model assumptions. The generated HMXB XLFs are feature-rich, deviating from the commonly assumed single-power law. We find a break in our synthetic XLF at luminosity $\sim 10^{38}$ erg s$^{-1}$, similar to observed XLFs. However, we find also a general overabundance of XRBs (up to a factor of $\sim$10 for certain model parameter combinations) driven primarily by XRBs with black hole accretors. Assumptions about the transient behavior of Be-XRBs, asymmetric supernova kicks, and common-envelope physics can significantly affect the shape and normalization of our synthetic XLFs. We find that less well-studied assumptions regarding the circularization of the orbit at the onset of Roche-lobe overflow and criteria for the formation of an X-ray emitting accretion disk around wind-accreting black holes can also impact our synthetic XLFs. Our study reveals the importance of large-scale parameter studies, highlighting the power of XRBs in constraining binary evolution theory. △ Less

Submitted 14 March, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 31 pages, 32 figures, Accepted by A&A. Fixed typos and updated references. Referee's comments were addressed

Journal ref: A&A 672, A99 (2023)

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

Authors: Keith Doore, Rafael T. Eufrasio, Bret D. Lehmer, Erik B. Monson, Antara Basu-Zych, Kristen Garofali

Abstract: 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 d… ▽ More 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. △ Less

Submitted 31 May, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 24 pages, 15 figures. 23 page body and 1 page appendix. Updated with published version

Journal ref: ApJ 931 (2022) 53

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

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Antara Basu-Zych, Kristen Garofali, Woodrow Gilbertson, Andrei Mesinger, Mihoko Yukita

Abstract: 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-format… ▽ More 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. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ (25 pages, 9 figures, and 4 tables)

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

Authors: Guang Yang, Médéric Boquien, William N. Brandt, Véronique Buat, Denis Burgarella, Laure Ciesla, Bret D. Lehmer, Katarzyna E. Małek, George Mountrichas, Casey Papovich, Estelle Pons, Marko Stalevski, Patrice Theulé, Shifu Zhu

Abstract: 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). Re… ▽ More 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 △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: ApJ accepted. The code is available at https://cigale.lam.fr/news/

The Stellar Age Dependence of X-ray Emission from Normal Star-Forming Galaxies in the GOODS Fields

Authors: Woodrow Gilbertson, Bret Lehmer, Keith Doore, Rafael Eufrasio, Antara Basu-Zych, William Brandt, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Bin Luo, Paolo Tozzi, Fabio Vito, Benjamin F. Williams, Yongquan Xue

Abstract: The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve wit… ▽ More The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve with redshift across z = 0-2 following L_HMXB/SFR ~ 1 + z and L_LMXB/M_star ~ (1 + z)^{2-3}. However, these measurements alone do not directly reveal the physical mechanisms behind the redshift evolution of X-ray binaries (XRBs). We derive star-formation histories for a sample of 344 normal galaxies in the CDFs, using spectral energy distribution (SED) fitting of FUV-to-FIR photometric data, and construct a self-consistent, age-dependent model of the X-ray emission from the galaxies. Our model quantifies how X-ray emission from hot gas and XRB populations vary as functions of host stellar-population age. We find that (1) the ratio L_X/M_star declines by a factor of ~1000 from 0-10 Gyr and (2) the X-ray SED becomes harder with increasing age, consistent with a scenario in which the hot gas contribution to the X-ray SED declines quickly for ages above 10 Myr. When dividing our sample into subsets based on metallicity, we find some indication that L_X/M_star is elevated for low-metallicity galaxies, consistent with recent studies of X-ray scaling relations. However, additional statistical constraints are required to quantify both the age and metallicity dependence of X-ray emission from star-forming galaxies. △ Less

Submitted 6 December, 2021; originally announced December 2021.

On the Impact of Inclination-Dependent Attenuation on Derived Star Formation Histories: Results from Disk Galaxies in the Great Observatories Origins Deep Survey Fields

Authors: Keith Doore, Rafael T. Eufrasio, Bret D. Lehmer, Erik B. Monson, Antara Basu-Zych, Kristen Garofali, Andrew Ptak

Abstract: We develop and implement an inclination-dependent attenuation prescription for spectral energy distribution (SED) fitting and study its impact on derived star-formation histories. We apply our prescription within the SED fitting code Lightning to a clean sample of 82, z=0.21-1.35 disk-dominated galaxies in the Great Observatories Origins Deep Survey North and South fields. To compare our inclinati… ▽ More We develop and implement an inclination-dependent attenuation prescription for spectral energy distribution (SED) fitting and study its impact on derived star-formation histories. We apply our prescription within the SED fitting code Lightning to a clean sample of 82, z=0.21-1.35 disk-dominated galaxies in the Great Observatories Origins Deep Survey North and South fields. To compare our inclination-dependent attenuation prescription with more traditional fitting prescriptions, we also fit the SEDs with the inclination-independent Calzetti et al. (2000) attenuation curve. From this comparison, we find that fits to a subset of 58, z < 0.7 galaxies in our sample, utilizing the Calzetti et al. (2000) prescription, recover similar trends with inclination as the inclination-dependent fits for the far-UV-band attenuation and recent star-formation rates. However, we find a difference between prescriptions in the optical attenuation (AV) that is strongly correlated with inclination (p-value < 10^-11). For more face-on galaxies, with i < 50 deg, (edge-on, i = 90 deg), the average derived AV is 0.31 +\- 0.11 magnitudes lower (0.56 +\- 0.16 magnitudes higher) for the inclination-dependent model compared to traditional methods. Further, the ratio of stellar masses between prescriptions also has a significant (p-value < 10^-2) trend with inclination. For i = 0-65 deg, stellar masses are systematically consistent between fits, with log(Mstar_inc/Mstar_Calzetti) = -0.05 +/- 0.03 dex and scatter of 0.11 dex. However, for i = 80-90 deg, derived stellar masses are lower for the Calzetti et al. (2000) fits by an average factor of 0.17 +\- 0.03 dex and scatter of 0.13 dex. Therefore, these results suggest that SED fitting assuming the Calzetti et al. (2000) attenuation law potentially underestimates stellar masses in highly inclined disk-dominated galaxies. △ Less

Submitted 8 December, 2021; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 38 pages, 22 figures. 30 page body and 8 page appendix. Updated with published version

Journal ref: ApJ 923 (2021) 26

On the Nature of AGN and Star Formation Enhancement in the $z = 3.1$ SSA22 Protocluster: The HST WFC3 IR View

Authors: Erik B. Monson, Bret D. Lehmer, Keith Doore, Rafael T. Eufrasio, Brett Bonine, David M. Alexander, Chris M. Harrison, Mariko Kubo, Kameswara B. Mantha, Cristian Saez, Amber Straughn, Hideki Umehata

Abstract: We examine possible environmental sources of the enhanced star formation and active galactic nucleus (AGN) activity in the $z = 3.09$ SSA22 protocluster using Hubble WFC3 F160W ($\sim1.6\ \rm μm$) observations of the SSA22 field, including new observations centered on eight X-ray selected protocluster AGN. To investigate the role of mergers in the observed AGN and star formation enhancement, we ap… ▽ More We examine possible environmental sources of the enhanced star formation and active galactic nucleus (AGN) activity in the $z = 3.09$ SSA22 protocluster using Hubble WFC3 F160W ($\sim1.6\ \rm μm$) observations of the SSA22 field, including new observations centered on eight X-ray selected protocluster AGN. To investigate the role of mergers in the observed AGN and star formation enhancement, we apply both quantitative (Sérsic-fit and Gini-$M_{20}$) and visual morphological classifications to F160W images of protocluster Lyman break galaxies (LBGs) in the fields of the X-ray AGN and $z \sim 3$ field LBGs in SSA22 and GOODS-N. We find no statistically significant differences between the morphologies and merger fractions of protocluster and field LBGs, though we are limited by small number statistics in the protocluster. We also fit the UV-to-near-IR spectral energy distributions (SED) of F160W-detected protocluster and field LBGs to characterize their stellar masses and star formation histories (SFH). We find that the mean protocluster LBG is by a factor of $\sim2$ times more massive and more attenuated than the mean $z \sim 3$ field LBG. We take our results to suggest that ongoing mergers are not more common among protocluster LBGs than field LBGs, though protocluster LBGs appear to be more massive. We speculate that the larger mass of the protocluster LBGs contributes to the enhancement of SMBH mass and accretion rate in the protocluster, which in turn drives the observed protocluster AGN enhancement. △ Less

Submitted 23 September, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 43 pages, 18 figures. 32 page body and 11 page appendix. Updated with published version

Journal ref: ApJ 919 (2021), 51

The Heraklion Extragalactic Catalogue (HECATE): a value-added galaxy catalogue for multi-messenger astrophysics

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Konstantinos Kouroumpatzakis, Bret Lehmer, Andrew Ptak

Abstract: We present the Heraklion Extragalactic Catalogue, or HECATE, an all-sky value-added galaxy catalogue, aiming to facilitate present and future multi-wavelength and multi-messenger studies in the local Universe. It contains 204,733 galaxies up to a redshift of 0.047 (D<200 Mpc), and it is >50% complete in terms of the B-band luminosity density at distances in the 0-170 Mpc range. By incorporating an… ▽ More We present the Heraklion Extragalactic Catalogue, or HECATE, an all-sky value-added galaxy catalogue, aiming to facilitate present and future multi-wavelength and multi-messenger studies in the local Universe. It contains 204,733 galaxies up to a redshift of 0.047 (D<200 Mpc), and it is >50% complete in terms of the B-band luminosity density at distances in the 0-170 Mpc range. By incorporating and homogenising data from astronomical databases and multi-wavelength surveys, the catalogue offers positions, sizes, distances, morphological classifications, star-formation rates, stellar masses, metallicities, and nuclear activity classifications. This wealth of information can enable a wide-range of applications, such as: (i) demographic studies of extragalactic sources, (ii) initial characterisation of transient events, and (iii) searches for electromagnetic counterparts of gravitational-wave events. The catalogue is publicly available to the community at a dedicated portal, which will also host future extensions in terms of the covered volume, and data products. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 21 pages. 11 figures, 5 tables. This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society, published by Oxford University Press on behalf of the Royal Astronomical Society

The XMM-SERVS survey: XMM-Newton point-source catalogs for the W-CDF-S and ELAIS-S1 fields

Authors: Q. Ni, W. N. Brandt, C. -T. Chen, B. Luo, K. Nyland, G. Yang, F. Zou, J. Aird, D. M. Alexander, F. E. Bauer, M. Lacy, B. D. Lehmer, L. Mallick, M. Salvato, D. P. Schneider, P. Tozzi, I. Traulsen, M. Vaccari, C. Vignali, F. Vito, Y. Xue, M. Banerji, K. Chow, A. Comastri, A. Del Moro , et al. (8 additional authors not shown)

Abstract: We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively;… ▽ More We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively; 1.8 Ms and 0.9 Ms exposures remain after flare filtering. The survey in W-CDF-S has a flux limit of 1.0 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 4053 sources are detected in total. The survey in ELAIS-S1 has a flux limit of 1.3 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 2630 sources are detected in total. Reliable optical-to-IR multiwavelength counterpart candidates are identified for $\approx$ 89% of the sources in W-CDF-S and $\approx$ 87% of the sources in ELAIS-S1. 3186 sources in W-CDF-S and 1985 sources in ELAIS-S1 are classified as AGNs. We also provide photometric redshifts for X-ray sources; $\approx$ 84% of the 3319/2001 sources in W-CDF-S/ELAIS-S1 with optical-to-NIR forced photometry available have either spectroscopic redshifts or high-quality photometric redshifts. The completion of the XMM-Newton observations in the W-CDF-S and ELAIS-S1 fields marks the end of the XMM-SERVS survey data gathering. The $\approx$ 12,000 point-like X-ray sources detected in the whole $\approx$ 13 deg$^2$ XMM-SERVS survey will benefit future large-sample AGN studies. △ Less

Submitted 19 June, 2021; originally announced June 2021.

Comments: 40 pages, 35 figures. Accepted for publication in ApJS. Data products available at: https://personal.psu.edu/wnb3/xmmservs/xmmservs.html

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

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Antara Basu-Zych, Keith Doore, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Benjamin J. Williams, Andreas Zezas, Luidhy Santana-Silva

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ

On the X-ray Spectral Energy Distributions of Star-Forming Galaxies: the 0.3-30 keV Spectrum of the Low-Metallicity Starburst Galaxy VV 114

Authors: Kristen Garofali, Bret D. Lehmer, Antara Basu-Zych, Lacey A. West, Daniel Wik, Mihoko Yukita, Neven Vulic, Andrew Ptak, Ann Hornschemeier

Abstract: Binary population synthesis combined with cosmological models suggest that X-ray emission from star-forming galaxies, consisting primarily of emission from X-ray binaries (XRBs) and the hot interstellar medium (ISM), could be an important, and perhaps dominant, source of heating of the intergalactic medium prior to the epoch of reionization. However, such models rely on empirical constraints for t… ▽ More Binary population synthesis combined with cosmological models suggest that X-ray emission from star-forming galaxies, consisting primarily of emission from X-ray binaries (XRBs) and the hot interstellar medium (ISM), could be an important, and perhaps dominant, source of heating of the intergalactic medium prior to the epoch of reionization. However, such models rely on empirical constraints for the X-ray spectral energy distributions (SEDs) of star-forming galaxies, which are currently lacking for low-metallicity galaxies. Using a combination of Chandra, XMM-Newton, and NuSTAR observations, we present new constraints on the 0.3-30 keV SED of the low-metallicity starburst galaxy VV 114, which is known to host several ultra-luminous X-ray sources (ULXs) with luminosities above 10$^{40}$ erg s$^{-1}$. We use an archival Chandra observation of VV 114 to constrain the contributions to the X-ray SED from the major X-ray emitting components of the galaxy, and newly acquired, nearly simultaneous XMM-Newton and NuSTAR observations to extend the spectral model derived from Chandra to cover the 0.3-30 keV range. Using our best-fit galaxy-wide spectral model, we derive the 0.3-30 keV SED of VV 114, which we find is dominated by emission from the XRB population, and in particular ULXs, at energies > 1.5 keV, and which we find to have an elevated galaxy-integrated X-ray luminosity per unit star formation rate relative to higher-metallicity star-forming galaxies. We discuss our results in terms of the effect of metallicity on XRB populations and the hot ISM, and the importance of X-ray emission from star-forming galaxies in the high redshift Universe. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: ApJ accepted. 25 pages, 8 figures, 4 tables

Chandra reveals a luminous Compton-thick QSO powering a $Lyα$ blob in a $z=4$ starbursting protocluster

Authors: Fabio Vito, William Nielsen Brandt, Bret Darby Lehmer, Cristian Vignali, Fan Zou, Franz Erik Bauer, Malcolm Bremer, Roberto Gilli, Rob J. Ivison, Cristiana Spingola

Abstract: Galaxy clusters in the local universe descend from high-redshift overdense regions known as protoclusters. The large gas reservoirs and high rate of galaxy interaction in protoclusters are expected to trigger star-formation activity and luminous SMBH accretion in the host galaxies. We investigated the AGN content of a gas-rich and starbursting protocluster at $z=4$, known as the Distant Red Core (… ▽ More Galaxy clusters in the local universe descend from high-redshift overdense regions known as protoclusters. The large gas reservoirs and high rate of galaxy interaction in protoclusters are expected to trigger star-formation activity and luminous SMBH accretion in the host galaxies. We investigated the AGN content of a gas-rich and starbursting protocluster at $z=4$, known as the Distant Red Core (DRC). We observed with Chandra (139 ks) the 13 identified members of the structure, and searched for luminous and possibly obscured AGN among them. We also tested whether a hidden AGN can power the $Lyα$ blob (LAB) detected with VLT/MUSE in the DRC. We detected obscured X-ray emission from the two most gas-rich members of the DRC, named DRC-1 and DRC-2. Both of them are resolved into multiple interacting clumps in high-resolution ALMA and HST observations. In particular, DRC-2 is found to host a luminous ($L_{2-10\,\mathrm{keV}}\approx3\times10^{45}\,\mathrm{erg\,s^{-1}}$) Compton-thick ($N_H\gtrsim10^{24}\,\mathrm{cm^{-2}}$) QSO, comparable to the most luminous QSOs known at all cosmic times. The AGN fraction among DRC members is consistent with results found for lower redshift protoclusters. However, X-ray stacking analysis reveals that SMBH accretion is likely also taking place in other DRC galaxies that are not detected individually by Chandra. Our results point toward the presence of a strong link between large gas reservoirs, galaxy interactions, and luminous and obscured nuclear activity in protocluster members. The powerful and obscured QSO detected in DRC-2 is likely powering the nearby LAB detected with VLT/MUSE, possibly through photoionization; however, we propose that the diffuse $Lyα$ emission may be due to gas shocked by a massive outflow launched by DRC-2 over a $\approx10$ kpc scale. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 13 pages, accepted for publication in A&A

Journal ref: A&A 642, A149 (2020)

A census of ultraluminous X-ray sources in the local Universe

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Bret Lehmer, Andrew Ptak

Abstract: Using the Chandra Source Catalog 2.0 and a newly compiled catalogue of galaxies in the local Universe, we deliver a census of ultraluminous X-ray source (ULX) populations in nearby galaxies. We find 629 ULX candidates in 309 galaxies with distance smaller than 40\,Mpc. The foreground/background contamination is ${\sim}20\%$. The ULX populations in bona-fide star-forming galaxies scale on average w… ▽ More Using the Chandra Source Catalog 2.0 and a newly compiled catalogue of galaxies in the local Universe, we deliver a census of ultraluminous X-ray source (ULX) populations in nearby galaxies. We find 629 ULX candidates in 309 galaxies with distance smaller than 40\,Mpc. The foreground/background contamination is ${\sim}20\%$. The ULX populations in bona-fide star-forming galaxies scale on average with star-formation rate (SFR) and stellar mass ($M_\star$) such that the number of ULXs per galaxy is $0.45^{+0.06}_{-0.09}\times\frac{\rm SFR}{\rm M_\odot\,yr^{-1}}{+}3.3^{+3.8}_{-3.2}\times\frac{M_\star}{\rm M_\odot}$. The scaling depends strongly on the morphological type. This analysis shows that early spiral galaxies contain an additional population of ULXs that scales with $M_\star$. We also confirm the strong anti-correlation of the ULX rate with the host galaxy's metallicity. In the case of early-type galaxies we find that there is a non-linear dependence of the number of ULXs with $M_\star$, which is interpreted as the result of star-formation history differences. Taking into account age and metallicity effects, we find that the predictions from X-ray binary population synthesis models are consistent with the observed ULX rates in early-type galaxies, as well as, spiral/irregular galaxies. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 21 pages, 19 figures, 6 tables. Accepted for publication in MNRAS

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

Authors: Bret D. Lehmer, Andrew P. Ferrell, Keith Doore, Rafael T. Eufrasio, Erik B. Monson, David M. Alexander, Antara Basu-Zych, William N. Brandt, Greg Sivakoff, Panayiotis Tzanavaris, Mihoko Yukita, Tassos Fragos, Andrew Ptak

Abstract: 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-ra… ▽ More 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. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: Accepted for publication in ApJS. Data products and catalogs available at https://lehmer.uark.edu/downloads/

Neutron Stars and Black Holes in the Small Magellanic Cloud: The SMC NuSTAR Legacy Survey

Authors: M. Lazzarini, B. F. Williams, A. E. Hornschemeier, V. Antoniou, G. Vasilopoulos, F. Haberl, N. Vulic, M. Yukita, A. Zezas, A. Bodaghee, B. D. Lehmer, T. J. Maccarone, A. Ptak, D. Wik, F. M. Fornasini, Jaesub Hong, J. A. Kennea, J. A. Tomsick, T. Venters, A. Udalski, A. Cassity

Abstract: We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$σ$ significance level (4… ▽ More We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$σ$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$σ$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$σ$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 8 tables, 18 figures

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Panayiotis Tzanavaris, Antara Basu-Zych, Tassos Fragos, Andrea Prestwich, Mihoko Yukita, Andreas Zezas, Ann E. Hornschemeier, Andrew Ptak

Abstract: We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs)… ▽ More We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs) and stellar masses (M*) to be measured on subgalactic scales. We divided the 2478 X-ray detected sources into 21 subsamples in bins of specific-SFR (sSFR = SFR/M*) and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs (HMXBs), low-mass XRBs (LMXBs), and background sources from the cosmic X-ray background (CXB) that respectively scale with SFR, M*, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above ~10^38 erg/s and elliptical galaxies that have relatively rich populations of globular clusters (GCs) show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

A magnetar-powered X-ray transient as the aftermath of a binary neutron-star merger

Authors: Y. Q. Xue, X. C. Zheng, Y. Li, W. N. Brandt, B. Zhang, B. Luo, B. B. Zhang, F. E. Bauer, H. Sun, B. D. Lehmer, X. F. Wu, G. Yang, X. Kong, J. Y. Li, M. Y. Sun, J. -X. Wang, F. Vito

Abstract: Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good f… ▽ More Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 29 pages, 4 figures, 3 tables, published in Nature on 11 April 2019

Journal ref: Nature, Vol. 568, 198-201 (2019)

Astro 2020 Science White Paper: Time Domain Studies of Neutron Star and Black Hole Populations: X-ray Identification of Compact Object Types

Authors: N. Vulic, A. E. Hornschemeier, V. Antoniou, A. R. Basu-Zych, B. Binder, F. M. Fornasini, F. Furst, F. Haberl, M. Heida, B. D. Lehmer, T. J. Maccarone, A. F. Ptak, G. R. Sivakoff, P. Tzanavaris, D. R. Wik, B. F. Williams, J. Wilms, M. Yukita, A. Zezas

Abstract: What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, a… ▽ More What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates. △ Less

Submitted 23 March, 2019; originally announced March 2019.

Comments: 9 pages, 2 figures. Submitted to the Astro2020 Decadal Survey

Electromagnetic Window into the Dawn of Black Holes

Authors: Zoltán Haiman, William N. Brandt, Alexey Vikhlinin, Jillian Bellovary, Elena Gallo, Jenny Greene, Kohei Inayoshi, Joseph Lazio, Bret Lehmer, Bin Luo, Piero Madau, Priya Natarajan, Feryal Özel, Fabio Pacucci, Alberto Sesana, Daniel Stern, Christian Vignali, Eli Visbal, Fabio Vito, Marta Volonteri, Joan Wrobel, Emanuele Berti, Volker Bromm, Greg Bryan, Nico Cappelluti , et al. (6 additional authors not shown)

Abstract: Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanyi… ▽ More Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

On the Nature of the X-ray Emission from the Ultraluminous X-ray Source, M33 X-8: New Constraints from NuSTAR and XMM-Newton

Authors: Lacey A. West, Bret D. Lehmer, Daniel Wik, Jun Yang, Dominic J. Walton, Vallia Antoniou, Frank Haberl, Ann Hornschemeier, Thomas J. Maccarone, Paul P. Plucinsky, Andrew Ptak, Benjamin F. Williams, Neven Vulic, Mihoko Yukita, Andreas Zezas

Abstract: We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 k… ▽ More We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 keV spectrum of X-8 can be characterized using an advection-dominated accretion disk model. We find that when fitting to our NuSTAR and XMM-Newton observations, an additional high-energy (>10 keV) Comptonization component is required, which allows us to rule out single advection-dominated disk and classical sub-Eddington models. With our new constraints, we analyze XMM-Newton data taken over the last 17 years to show that small (~30%) variations in the 0.3-10 keV flux of M33 X-8 result in spectral changes similar to those observed for other ULXs. The two most likely phenomenological scenarios suggested by the data are degenerate in terms of constraining the nature of the accreting compact object (i.e., black hole versus neutron star). We further present a search for pulsations using our suite of data; however, no clear pulsations are detected. Future observations designed to observe M33 X-8 at different flux levels across the full 0.3-30 keV range would significantly improve our constraints on the nature of this important source. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted for publication in ApJ (15 pages, 4 tables, 6 figures)

Variability-selected low-luminosity active galactic nuclei candidates in the 7 Ms Chandra Deep Field-South

Authors: N. Ding, B. Luo, W. N. Brandt, M. Paolillo, G. Yang, B. D. Lehmer, O. Shemmer, D. P. Schneider, P. Tozzi, Y. Q. Xue, X. C. Zheng, Q. S. Gu, A. M. Koekemoer, C. Vignali, F. Vito, J. X. Wang

Abstract: In deep X-ray surveys, active galactic nuclei (AGNs) with a broad range of luminosities have been identified. However, cosmologically distant low-luminosity AGN (LLAGN, $L_{\mathrm{X}} \lesssim 10^{42}$ erg s$^{-1}$) identification still poses a challenge due to significant contamination from host galaxies. Based on the 7 Ms Chandra Deep Field-South (CDF-S) survey, the longest timescale (… ▽ More In deep X-ray surveys, active galactic nuclei (AGNs) with a broad range of luminosities have been identified. However, cosmologically distant low-luminosity AGN (LLAGN, $L_{\mathrm{X}} \lesssim 10^{42}$ erg s$^{-1}$) identification still poses a challenge due to significant contamination from host galaxies. Based on the 7 Ms Chandra Deep Field-South (CDF-S) survey, the longest timescale ($\sim 17$ years) deep X-ray survey to date, we utilize an X-ray variability selection technique to search for LLAGNs that remain unidentified among the CDF-S X-ray sources. We find 13 variable sources from 110 unclassified CDF-S X-ray sources. Except for one source which could be an ultraluminous X-ray source, the variability of the remaining 12 sources is most likely due to accreting supermassive black holes. These 12 AGN candidates have low intrinsic X-ray luminosities, with a median value of $7 \times10^{40}$ erg s$^{-1}$. They are generally not heavily obscured, with an average effective power-law photon index of 1.8. The fraction of variable AGNs in the CDF-S is independent of X-ray luminosity and is only restricted by the total number of observed net counts, confirming previous findings that X-ray variability is a near-ubiquitous property of AGNs over a wide range of luminosities. There is an anti-correlation between X-ray luminosity and variability amplitude for high-luminosity AGNs, but as the luminosity drops to $\lesssim 10^{42}$ erg s$^{-1}$, the variability amplitude no longer appears dependent on the luminosity. The entire observed luminosity-variability trend can be roughly reproduced by an empirical AGN variability model based on a broken power-law power spectral density function. △ Less

Submitted 22 October, 2018; originally announced October 2018.

Comments: 18 pages, 11 figures, accepted for publication in ApJ

Black Holes and Neutron Stars in Nearby Galaxies: Insights from NuSTAR

Authors: Neven Vulic, Ann E. Hornschemeier, Daniel R. Wik, Mihoko Yukita, Andreas Zezas, Andrew F. Ptak, Bret D. Lehmer, Valsamo Antoniou, Thomas J. Maccarone, Benjamin F. Williams, Francesca M. Fornasini

Abstract: Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected samp… ▽ More Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc having stellar masses ($M_{\star}$) $10^{7-11}$ $M_{\odot}$ and star formation rates (SFR) $\approx0.01-15$ $M_{\odot}$ yr$^{-1}$. We detect 128 NuSTAR sources to a sensitivity of $\approx10^{38}$ erg s$^{-1}$. Using NuSTAR color-intensity and color-color diagrams we classify 43 of these sources as candidate NS and 47 as candidate BH. We further subdivide BH by accretion states (soft, intermediate, and hard) and NS by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using 8 normal (Milky Way-type) galaxies in the sample, we confirm the relation between SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrain galaxy X-ray point source luminosity using the relation $L_{\rm{X}}=αM_{\star}+β\text{SFR}$, finding agreement with previous work. The XLF of all sources in the 4-25 and 12-25 keV energy bands matches with the $α=1.6$ slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 $M_{\odot}$ NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25keV energy bands. We calculate the overall ratio of BH to NS to be $\approx1$ for 4-25 keV and $\approx2$ for 12-25 keV. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: 38 pages, 12 figures, 8 tables. ApJ, in press

A high dust emissivity index beta for a CO-faint galaxy in a filamentary Lyman-alpha nebula at z=3.1

Authors: Yuta Kato, Yuichi Matsuda, Daisuke Iono, Bunyo Hatsukade, Hideki Umehata, Kotaro Kohno, David M. Alexander, Yiping Ao, Scott C. Chapman, Matthew Hayes, Mariko Kubo, Bret D. Lehmer, Matthew A. Malkan, Tomonari Michiyama, Tohru Nagao, Tomoki Saito, Ichi Tanaka, Yoshiaki Taniguchi

Abstract: We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Lyα nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 σ from one of the ALMA continuum sources associated with the Lyα filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)… ▽ More We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Lyα nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 σ from one of the ALMA continuum sources associated with the Lyα filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of τ(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index β=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high β can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through τ(dep) and β, and thus we can investigate dust chemical compositions even in the early Universe. △ Less

Submitted 23 July, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: 6 pagers, 3 figures

Young Accreting Compact Objects in M31: The Combined Power of NuSTAR, Chandra, and Hubble

Authors: M. Lazzarini, A. E. Hornschemeier, B. F. Williams, D. Wik, N. Vulic, M. Yukita, A. Zezas, A. R. Lewis, M. Durbin, A. Ptak, A. Bodaghee, B. D. Lehmer, V. Antoniou, T. Maccarone

Abstract: We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object… ▽ More We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard state black holes, pulsars, and non-magnetized neutron stars associated with optical point source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates --- X-ray sources with UV-bright point source optical counterpart candidates --- are found in regions with star formation bursts less than 50 Myr ago, with 3 associated with young stellar ages (<10 Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403. △ Less

Submitted 8 June, 2018; originally announced June 2018.

Comments: 16 pages, 9 figures, 4 tables

A long hard-X-ray look at the dual active galactic nuclei of M51 with NuSTAR

Authors: M. Brightman, M. Baloković, M. Koss, D. M. Alexander, A. Annuar, H. Earnshaw, P. Gandhi, F. A. Harrison, A. E. Hornschemeier, B. Lehmer, M. C. Powell, A. Ptak, B. Rangelov, T. P. Roberts, D. Stern, D. J. Walton, A. Zezas

Abstract: We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and bot… ▽ More We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($λ_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources. △ Less

Submitted 13 September, 2018; v1 submitted 30 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJ

The XMM-SERVS survey: new XMM-Newton point-source catalog for the XMM-LSS field

Authors: C. -T. J. Chen, W. N. Brandt, B. Luo, P. Ranalli, G. Yang, D. M. Alexander, F. E. Bauer, D. D. Kelson, M. Lacy, K. Nyland, P. Tozzi, F. Vito, M. Cirasuolo, R. Gilli, M. J. Jarvis, B. D. Lehmer, M. Paolillo, D. P. Schneider, O. Shemmer, I. Smail, M. Sun, M. Tanaka, M. Vaccari, C. Vignali, Y. Q. Xue , et al. (6 additional authors not shown)

Abstract: We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling… ▽ More We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts. △ Less

Submitted 20 April, 2018; originally announced April 2018.

Comments: MNRAS, accepted. 34 pages, 25 figures, and 8 tables. The data products are available at this http url: http://personal.psu.edu/wnb3/xmmservs/xmmservs.html

On the Spatially Resolved Star-Formation History in M51 II: X-ray Binary Population Evolution

Authors: B. D. Lehmer, R. T Eufrasio, L. Markwardt, A. Zezas, A. Basu-Zych, T. Fragos, A. E. Hornschemeier, A. Ptak, P. Tzanavaris, M. Yukita

Abstract: We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next… ▽ More We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ (19 pages, 15 figures, and 2 tables)

On the Spatially Resolved Star Formation History in M51 I: Hybrid UV+IR Star Formation Laws and IR Emission from Dust Heated by Old Stars

Authors: R. T. Eufrasio, B. D. Lehmer, A. Zezas, E. Dwek, R. G. Arendt, A. Basu-Zych, T. Wiklind, M. Yukita, T. Fragos, A. E. Hornschemeier, L. Markwardt, A. Ptak, P. Tzanavaris

Abstract: We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determin… ▽ More We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

High-redshift AGN in the Chandra Deep Fields: the obscured fraction and space density of the sub-$L_*$ population

Authors: F. Vito, W. N. Brandt, G. Yang, R. Gilli, B. Luo, C. Vignali, Y. Q. Xue, A. Comastri, A. M. Koekemoer, B. D. Lehmer, T. Liu, M. Paolillo, P. Ranalli, D. P. Schneider, O. Shemmer, M. Volonteri, J. Wang

Abstract: We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower da… ▽ More We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Comments: Accepted 2017 September 15 (MNRAS). 31 Pages, 25 figures

Galaxy-scale Bars in Late-type Sloan Digital Sky Survey Galaxies Do Not Influence the Average Accretion Rates of Supermassive Black Holes

Authors: Andy D. Goulding, Eve Matthaey, Jenny E. Greene, Ryan C. Hickox, David M. Alexander, William R. Forman, Christine Jones, Bret D. Lehmer, Samuel Griffis, Svilen Kanev, Mehdi Oulmakki

Abstract: Galaxy-scale bars are expected to provide an effective means for driving material towards the central region in spiral galaxies, and possibly feeding supermassive black holes (BHs). Here we present a statistically-complete study of the effect of bars on average BH accretion. From a well-selected sample of 50,794 spiral galaxies (with M* ~ 0.2-30 x 10^10 Msun) extracted from the Sloan Digital Sky S… ▽ More Galaxy-scale bars are expected to provide an effective means for driving material towards the central region in spiral galaxies, and possibly feeding supermassive black holes (BHs). Here we present a statistically-complete study of the effect of bars on average BH accretion. From a well-selected sample of 50,794 spiral galaxies (with M* ~ 0.2-30 x 10^10 Msun) extracted from the Sloan Digital Sky Survey Galaxy Zoo 2 project, we separate those sources considered to contain galaxy-scale bars from those that do not. Using archival data taken by the Chandra X-ray Observatory, we identify X-ray luminous (L_X >~ 10^41 erg/s) active galactic nuclei (AGN) and perform an X-ray stacking analysis on the remaining X-ray undetected sources. Through X-ray stacking, we derive a time-averaged look at accretion for galaxies at fixed stellar mass and star formation rate, finding that the average nuclear accretion rates of galaxies with bar structures are fully consistent with those lacking bars (Mdot_acc ~ 3 x 10^-5 Msun/yr). Hence, we robustly conclude that large-scale bars have little or no effect on the average growth of BHs in nearby (z < 0.15) galaxies over gigayear timescales. △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 10 pages, 5 figures, accepted to ApJ

Further constraints on variations in the IMF from LMXB populations

Authors: Mark B. Peacock, Stephen E. Zepf, Arunav Kundu, Thomas J. Maccarone, Bret D. Lehmer, Claudia Maraston, Anthony H. Gonzalez, Rafael T. Eufrasio, David A. Coulter

Abstract: We present constraints on variations in the initial mass function (IMF) of nine local early-type galaxies based on their low mass X-ray binary (LMXB) populations. Comprised of accreting black holes and neutron stars, these LMXBs can be used to constrain the important high mass end of the IMF. We consider the LMXB populations beyond the cores of the galaxies ($>0.2R_{e}$; covering $75-90\%$ of thei… ▽ More We present constraints on variations in the initial mass function (IMF) of nine local early-type galaxies based on their low mass X-ray binary (LMXB) populations. Comprised of accreting black holes and neutron stars, these LMXBs can be used to constrain the important high mass end of the IMF. We consider the LMXB populations beyond the cores of the galaxies ($>0.2R_{e}$; covering $75-90\%$ of their stellar light) and find no evidence for systematic variations of the IMF with velocity dispersion ($σ$). We reject IMFs which become increasingly bottom heavy with $σ$, up to steep power-laws (exponent, $α>2.8$) in massive galaxies ($σ>300$km/s), for galactocentric radii $>1/4\ R_{e}$. Previously proposed IMFs that become increasingly bottom heavy with $σ$ are consistent with these data if only the number of low mass stars $(<0.5M_{\odot}$) varies. We note that our results are consistent with some recent work which proposes that extreme IMFs are only present in the central regions of these galaxies. We also consider IMFs that become increasingly top-heavy with $σ$, resulting in significantly more LMXBs. Such a model is consistent with these observations, but additional data are required to significantly distinguish between this and an invariant IMF. For six of these galaxies, we directly compare with published IMF mismatch parameters from the Atlas3D survey, $α_{dyn}$. We find good agreement with the LMXB population if galaxies with higher $α_{dyn}$ have more top-heavy IMFs -- although we caution that our sample is quite small. Future LMXB observations can provide further insights into the origin of $α_{dyn}$ variations. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

Black-Hole Growth is Mainly Linked to Host-Galaxy Stellar Mass rather than Star Formation Rate

Authors: G. Yang, C. -T. J. Chen, F. Vito, W. N. Brandt, D. M. Alexander, B. Luo, M. Y. Sun, Y. Q. Xue, F. E. Bauer, A. M. Koekemoer, B. D. Lehmer, T. Liu, D. P. Schneider, O. Shemmer, J. R. Trump, C. Vignali, J. -X. Wang

Abstract: We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or… ▽ More We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs. △ Less

Submitted 25 May, 2017; v1 submitted 21 April, 2017; originally announced April 2017.

Comments: 18 pages, 11 figures, 3 tables; accepted for publication in ApJ

Identification of the Hard X-ray Source Dominating the E > 25 keV Emission of the Nearby Galaxy M31

Authors: M. Yukita, A. Ptak, A. E. Hornschemeier, D. Wik, T. J. Maccarone, K. Pottschmidt, A. Zezas, V. Antoniou, R. Ballhausen, B. D. Lehmer, A. Lien, B. Williams, F. Baganoff, P. T. Boyd, T. Enoto, J. Kennea, K. L. Page, Y. Choi

Abstract: We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observa… ▽ More We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail. △ Less

Submitted 21 March, 2017; originally announced March 2017.

Comments: 10 pages, 6 figures, accepted for publication in ApJ

Deep Chandra observations of NGC~7457, the X-ray point source populations of a low mass early-type galaxy

Authors: Mark B. Peacock, Stephen E. Zepf, Arunav Kundu, Thomas J. Maccarone, Bret D. Lehmer, Anthony H. Gonzalez, Claudia Maraston

Abstract: We present the X-ray point source population of NGC 7457 based on 124 ks of Chandra observations. Previous deep Chandra observations of low mass X-ray binaries (LMXBs) in early-type galaxies have typically targeted the large populations of massive galaxies. NGC 7457 is a nearby, early-type galaxy with a stellar luminosity of $1.7\times10^{10} L_{K\odot}$, allowing us to investigate the populations… ▽ More We present the X-ray point source population of NGC 7457 based on 124 ks of Chandra observations. Previous deep Chandra observations of low mass X-ray binaries (LMXBs) in early-type galaxies have typically targeted the large populations of massive galaxies. NGC 7457 is a nearby, early-type galaxy with a stellar luminosity of $1.7\times10^{10} L_{K\odot}$, allowing us to investigate the populations in a relatively low mass galaxy. We classify the detected X-ray sources into field LMXBs, globular cluster LMXBs, and background AGN based on identifying optical counterparts in new HST/ACS images. We detect 10 field LMXBs within the $r_{ext}$ ellipse of NGC 7457 (with semi-major axis $\sim$ 9.1 kpc, ellipticity = 0.55). The corresponding number of LMXBs with $L_{x}>2\times10^{37}erg/s$ per stellar luminosity is consistent with that observed in more massive galaxies, $\sim 7$ per $10^{10} L_{K\odot}$. We detect a small globular cluster population in these HST data and show that its colour distribution is likely bimodal and that its specific frequency is similar to that of other early type galaxies. However, no X-ray emission is detected from any of these clusters. Using published data for other galaxies, we show that this non-detection is consistent with the small stellar mass of these clusters. We estimate that 0.11 (and 0.03) LMXBs are expected per $10^{6}M_{\odot}$ in metal-rich (and metal-poor) globular clusters. This corresponds to 1100 (and 330) LMXBs per $10^{10} L_{K\odot}$, highlighting the enhanced formation efficiency of LMXBs in globular clusters. A nuclear X-ray source is detected with $L_{x}$ varying from $2.8-6.8\times10^{38}erg/s$. Combining this $L_{x}$ with a published dynamical mass estimate for the central SMBH in NGC 7457, we find that $L_{x}/L_{Edd}$ varies from $0.5-1.3\times10^{-6}$. △ Less

Submitted 2 January, 2017; originally announced January 2017.

Comments: 9 pages, 7 figures, 1 table, accepted for publication in MNRAS

Testing the Universality of the Stellar IMF with Chandra and HST

Authors: D. A. Coulter, B. D. Lehmer, R. T. Eufrasio, A. Kundu, T. Maccarone, M. Peacock, A. E. Hornschemeier, A. Basu-Zych, A. H. Gonzalez, C. Maraston, S. E. Zepf

Abstract: The stellar initial mass function (IMF), which is often assumed to be universal across unresolved stellar populations, has recently been suggested to be "bottom-heavy" for massive ellipticals. In these galaxies, the prevalence of gravity-sensitive absorption lines (e.g. Na I and Ca II) in their near-IR spectra implies an excess of low-mass ($m <= 0.5$ $M_\odot$) stars over that expected from a can… ▽ More The stellar initial mass function (IMF), which is often assumed to be universal across unresolved stellar populations, has recently been suggested to be "bottom-heavy" for massive ellipticals. In these galaxies, the prevalence of gravity-sensitive absorption lines (e.g. Na I and Ca II) in their near-IR spectra implies an excess of low-mass ($m <= 0.5$ $M_\odot$) stars over that expected from a canonical IMF observed in low-mass ellipticals. A direct extrapolation of such a bottom-heavy IMF to high stellar masses ($m >= 8$ $M_\odot$) would lead to a corresponding deficit of neutron stars and black holes, and therefore of low-mass X-ray binaries (LMXBs), per unit near-IR luminosity in these galaxies. Peacock et al. (2014) searched for evidence of this trend and found that the observed number of LMXBs per unit $K$-band luminosity ($N/L_K$) was nearly constant. We extend this work using new and archival Chandra X-ray Observatory (Chandra) and Hubble Space Telescope (HST) observations of seven low-mass ellipticals where $N/L_K$ is expected to be the largest and compare these data with a variety of IMF models to test which are consistent with the observed $N/L_K$. We reproduce the result of Peacock et al. (2014), strengthening the constraint that the slope of the IMF at $m >= 8$ $M_\odot$ must be consistent with a Kroupa-like IMF. We construct an IMF model that is a linear combination of a Milky Way-like IMF and a broken power-law IMF, with a steep slope ($α_1=$ $3.84$) for stars < 0.5 $M_\odot$ (as suggested by near-IR indices), and that flattens out ($α_2=$ $2.14$) for stars > 0.5 $M_\odot$, and discuss its wider ramifications and limitations. △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: Accepted for publication in ApJ; 7 pages, 2 figures, 1 table

The Chandra Deep Field-South Survey: 7 Ms Source Catalogs

Authors: B. Luo, W. N. Brandt, Y. Q. Xue, B. Lehmer, D. M. Alexander, F. E. Bauer, F. Vito, G. Yang, A. R. Basu-Zych, A. Comastri, R. Gilli, Q. -S. Gu, A. E. Hornschemeier, A. Koekemoer, T. Liu, V. Mainieri, M. Paolillo, P. Ranalli, P. Rosati, D. P. Schneider, O. Shemmer, I. Smail, M. Sun, P. Tozzi, C. Vignali , et al. (1 additional authors not shown)

Abstract: We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV… ▽ More We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$). △ Less

Submitted 18 December, 2016; v1 submitted 10 November, 2016; originally announced November 2016.

Comments: 31 pages, 31 figures, 8 tables, ApJS in press, minor text updates; full-resolution PDF version and data products available at http://www.astro.psu.edu/users/niel/cdfs/cdfs-chandra.html

The growth of typical star-forming galaxies and their super massive black holes across cosmic time since z~2

Authors: João Calhau, David Sobral, Andra Stroe, Philip Best, Ian Smail, Bret Lehmer, Chris Harrison, Alasdair Thomson

Abstract: Understanding galaxy formation and evolution requires studying the interplay between the growth of galaxies and the growth of their black holes across cosmic time. Here we explore a sample of Ha-selected star-forming galaxies from the HiZELS survey and use the wealth of multi-wavelength data in the COSMOS field (X-rays, far-infrared and radio) to study the relative growth rates between typical gal… ▽ More Understanding galaxy formation and evolution requires studying the interplay between the growth of galaxies and the growth of their black holes across cosmic time. Here we explore a sample of Ha-selected star-forming galaxies from the HiZELS survey and use the wealth of multi-wavelength data in the COSMOS field (X-rays, far-infrared and radio) to study the relative growth rates between typical galaxies and their central supermassive black holes, from z=2.23 to z=0. Typical star-forming galaxies at z~1-2 have black hole accretion rates (BHARs) of 0.001-0.01 Msun/yr and star formation rates (SFRs) of ~10-40 Msun/yr, and thus grow their stellar mass much quicker than their black hole mass (~3.3 orders of magnitude faster). However, ~3% of the sample (the sources detected directly in the X-rays) show a significantly quicker growth of the black hole mass (up to 1.5 orders of magnitude quicker growth than the typical sources). BHARs fall from z=2.23 to z=0, with the decline resembling that of star formation rate density or the typical SFR. We find that the average black hole to galaxy growth (BHAR/SFR) is approximately constant for star-forming galaxies in the last 11 Gyrs. The relatively constant BHAR/SFR suggests that these two quantities evolve equivalently through cosmic time and with practically no delay between the two. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 10 pages, 8 figures, accepted to MNRAS

A resolved map of the infrared excess in a Lyman Break Galaxy at z=3

Authors: M. P. Koprowski, K. E. K. Coppin, J. E. Geach, N. K. Hine, M. Bremer, S. C. Chapman, L. J. M. Davies, T. Hayashino, K. K. Knudsen, M. Kubo, B. D. Lehmer, Y. Matsuda, D. J. B. Smith, P. P. van der Werf, G. Violino, T. Yamada

Abstract: We have observed the dust continuum of ten z=3.1 Lyman Break Galaxies with the Atacama Large Millimeter/Submillimeter Array at ~450 mas resolution in Band 7. We detect and resolve the 870um emission in one of the targets with an integrated flux density of S(870)=(192+/-57) uJy, and measure a stacked 3-sigma signal of S(870)=(67+/-23) uJy for the remaining nine. The total infrared luminosities esti… ▽ More We have observed the dust continuum of ten z=3.1 Lyman Break Galaxies with the Atacama Large Millimeter/Submillimeter Array at ~450 mas resolution in Band 7. We detect and resolve the 870um emission in one of the targets with an integrated flux density of S(870)=(192+/-57) uJy, and measure a stacked 3-sigma signal of S(870)=(67+/-23) uJy for the remaining nine. The total infrared luminosities estimated from full spectral energy distribution fits are L(8-1000um)=(8.4+/-2.3)x10^10 Lsun for the detection and L(8-1000um)=(2.9+/-0.9)x10^10 Lsun for the stack. With HST ACS I-band imaging we map the rest-frame UV emission on the same scale as the dust, effectively resolving the 'infrared excess' (IRX=L_FIR/L_UV) in a normal galaxy at z=3. Integrated over the galaxy we measure IRX=0.56+/-0.15, and the galaxy-averaged UV slope is beta=-1.25+/-0.03. This puts the galaxy a factor of ~10 below the IRX-beta relation for local starburst nuclei of Meurer et al. (1999). However, IRX varies by more than a factor of 3 across the galaxy, and we conclude that the complex relative morphology of the dust relative to UV emission is largely responsible for the scatter in the IRX-beta relation at high-z. A naive application of a Meurer-like dust correction based on the UV slope would dramatically over-estimate the total star formation rate, and our results support growing evidence that when integrated over the galaxy, the typical conditions in high-z star-forming galaxies are not analogous to those in the local starburst nuclei used to establish the Meurer relation. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: Accepted in The Astrophysical Journal Letters