Showing 1–50 of 54 results for author: Basu-Zych, A

The first all-sky survey of star-forming galaxies with eROSITA: Scaling relations and a population of X-ray luminous starbursts

Authors: E. Kyritsis, A. Zezas, F. Haberl, P. Weber, A. Basu-Zych, N. Vulic, C. Maitra, S. Hämmerich, J. Wilms, M. Sasaki, A. Hornschemeier, A. Ptak, A. Merloni, J. Comparat

Abstract: We present a study of X-ray normal galaxies using data from the first all-sky scan of the eROSITA X-ray survey. eRASS1 provides the first unbiased X-ray census of normal galaxies allowing us to study the X-ray emission from XRBs and the hot ISM in the full range of stellar population parameters present in the local Universe. By combining the HECATE value-added galaxy catalogue with the eRASS1, we… ▽ More We present a study of X-ray normal galaxies using data from the first all-sky scan of the eROSITA X-ray survey. eRASS1 provides the first unbiased X-ray census of normal galaxies allowing us to study the X-ray emission from XRBs and the hot ISM in the full range of stellar population parameters present in the local Universe. By combining the HECATE value-added galaxy catalogue with the eRASS1, we study the X-ray emission from normal galaxies as a function of their SFR, M$_{*}$, Metallicity, and stellar population age. After applying optical and mid-IR activity classification criteria, we constructed a sample of 18790 star-forming galaxies with measurements of their L$_{X}$. By stacking the X-ray data in SFR-M$_{*}$-distance bins we study the correlation between the average L$_{X}$ and stellar population parameters. We also present updated L$_{\rm{X}}$-SFR and L$_{\rm{X}}$/SFR-Metallicity scaling relations accounting for the scatter dependence on the SFR. We find that the integrated L$_{X}$ of the HEC-eR1 star-forming galaxies is significantly elevated with respect to that expected from the current scaling relations. The observed scatter is also significantly larger. This excess persists even when we measure the average L$_{X}$ of galaxies in SFR-M$_{*}$-distance and metallicity bins and it is stronger in lower SFRs. The excess is not the result of hot gas, LMXBs, background AGN, LLAGN (including TDEs), or stochastic sampling of the XRB XLF. We find that while the excess correlates with lower metallicity, its primary driver is the age of the stellar populations. Our analysis reveals a sub-population of X-ray luminous starbursts with high sSFRs, low metallicities, and young stellar populations. This population drives upwards the X-ray scaling relations for star-forming galaxies, and has important implications for understanding the population of XRBs in the local and high-z Universe. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: Submitted to A&A

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)

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

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. H. Buckley, A. Brill, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A Farrell, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, C. Giuri, O. Gueta , et al. (64 additional authors not shown)

Abstract: The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four… ▽ More The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hours of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star forming galaxies, dwarf-spheroidal galaxies, and highly-variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. △ Less

Submitted 13 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: Version with corrected author list

Journal ref: Research Notes of the AAS, Volume 7, Number 1, 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)

The 21-cm signal from the Cosmic Dawn: metallicity dependence of high mass X-ray binaries

Authors: Harman Deep Kaur, Yuxiang Qin, Andrei Mesinger, Andrea Pallottini, Tassos Fragos, Antara Basu-Zych

Abstract: X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift… ▽ More X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs ($L_{\rm X}/{\rm SFR}$) should increase in metal-poor environments, typical of early galaxies. Here we study the impact of the metallicity ($Z$) dependence of $L_{\rm X}/{\rm SFR}$ during the EoH. For our fiducial models, galaxies with star formation rates of order $10^{-3} - 10^{-1}$ $M_\odot$ yr$^{-1}$ and metallicities of order $10^{-3} - 10^{-2}$ $Z_\odot$ are the dominant contributors to the X-ray background (XRB) during this period. Different $L_{\rm X}/{\rm SFR}$-$Z$ relations result in factors of $\sim$ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000h integration with the upcoming Square Kilometer Array (SKA), for two different $L_{\rm X}/{\rm SFR}$-$Z$ relations. We perform inference on these mock observations using the common simplification of a constant $L_{\rm X}/{\rm SFR}$, finding that constant $L_{\rm X}/{\rm SFR}$ models can recover the IGM evolution of the more complicated $L_{\rm X}/{\rm SFR}$-$Z$ simulations only during the EoH. At $z<10$, where the typical galaxies are more polluted, constant $L_{\rm X}/{\rm SFR}$ models over-predict the XRB and its relative contribution to the early stages of the reionization. △ Less

Submitted 3 May, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

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.

Quenching, bursting and galaxy shapes: colour transformation as a function of morphology

Authors: Camila de Sá-Freitas, Thiago Signorini Gonçalves, Reinaldo R. de Carvalho, Karín Menéndez-Delmestre, Paulo H. Barchi, Vitor M. Sampaio, Antara Basu-Zych, Behnam Darvish, Christopher Martin

Abstract: Different mechanisms for quenching star formation in galaxies are commonly invoked in the literature, but the relative impact of each one at different cosmic epochs is still unknown. In particular, the relation between these processes and morphological transformation remains poorly understood. In this work, we measure the effectiveness of changes in star formation rates by analysing a new paramete… ▽ More Different mechanisms for quenching star formation in galaxies are commonly invoked in the literature, but the relative impact of each one at different cosmic epochs is still unknown. In particular, the relation between these processes and morphological transformation remains poorly understood. In this work, we measure the effectiveness of changes in star formation rates by analysing a new parameter, the Star Formation Acceleration (SFA), as a function of galaxy morphology. This methodology is capable of identifying both bursting and quenching episodes that occurred in the preceding 300 Myrs. We use morphological classification catalogs based on Deep learning techniques. Our final sample has $\sim$14,200 spirals and $\sim$2,500 ellipticals. We find that elliptical galaxies in the transition region have median shorter quenching timescales ( $τ$ < 1 Gyr) than spirals ($τ\geq 1$ Gyr). This result conforms to the scenario in which major mergers and other violent processes play a fundamental role in galaxy evolution for most ellipticals, not only quenching star formation more rapidly but also playing a role in morphological transformation. We also find that $\sim$two thirds of galaxies bursting in the green valley in our sample are massive spirals ($M_\star \geq 10^{11.0}M_\odot$) with signs of disturbance. This is in accordance with the scenario where low mass galaxies are losing their gas in a interaction with a massive galaxy: while the former is quenching, the last is being refueled and going through a burst, showing signs of recent interaction. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 16 pages, 12 figures, accepted by MNRAS

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

The eROSITA Final Equatorial-Depth Survey (eFEDS): Presenting The Demographics of X-ray Emission From Normal Galaxies

Authors: N. Vulic, A. E. Hornschemeier, F. Haberl, A. R. Basu-Zych, E. Kyritsis, A. Zezas, M. Salvato, A. Ptak, A. Bogdan, K. Kovlakas, J. Wilms, M. Sasaki, T. Liu, A. Merloni, T. Dwelly, H. Brunner, G. Lamer, C. Maitra, K. Nandra, A. Santangelo

Abstract: The $\it{eROSITA}$ Final Equatorial Depth Survey (eFEDS), completed during the calibration and performance verification phase of the $\it{eROSITA}$ instrument on $\it{Spectrum\, Roentgen\, Gamma}$, delivers data at and beyond the final depth of the four-year $\it{eROSITA}$ all-sky survey (eRASS:8), $f_{0.5-2\,\text{ keV}}$ = $1.1\times10^{-14}$ erg s$^{-1}$ cm$^{2}$, over 140 deg$^{2}$. It provide… ▽ More The $\it{eROSITA}$ Final Equatorial Depth Survey (eFEDS), completed during the calibration and performance verification phase of the $\it{eROSITA}$ instrument on $\it{Spectrum\, Roentgen\, Gamma}$, delivers data at and beyond the final depth of the four-year $\it{eROSITA}$ all-sky survey (eRASS:8), $f_{0.5-2\,\text{ keV}}$ = $1.1\times10^{-14}$ erg s$^{-1}$ cm$^{2}$, over 140 deg$^{2}$. It provides the first view of normal galaxy X-ray emission from X-ray binaries (XRBs) and the hot interstellar medium at the full depth of eRASS:8. We use the Heraklion Extragalactic Catalogue (HECATE) of galaxies to correlate with eFEDS X-ray sources and identify 94 X-ray detected normal galaxies. We classify galaxies as star-forming, early-type, composite, and AGN using SDSS and 6dF optical spectroscopy. The eFEDS field harbours 37 normal galaxies: 36 late-type (star-forming) galaxies and 1 early-type galaxy. There are 1.9 times as many normal galaxies as predicted by scaling relations via SIXTE simulations, with an overabundance of late-type galaxies and a dearth of early-type galaxies. Dwarf galaxies with high specific star formation rate (SFR) have elevated L$_{\text{X}}$/SFR when compared with specific SFR and metallicity, indicating an increase in XRB emission due to low-metallicity. We expect that eRASS:8 will detect 12,500 normal galaxies, the majority of which will be star-forming, with the caveat that there are unclassified sources in eFEDS and galaxy catalogue incompleteness issues that could increase the actual number of detected galaxies over these current estimates. eFEDS observations detected a rare population of galaxies -- the metal-poor dwarf starbursts -- that do not follow known scaling relations. eRASS is expected to discover significant numbers of these high-redshift analogues, which are important for studying the heating of the intergalactic medium at high-redshift. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 17 pages, 10 figures, 3 tables. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Journal ref: A&A 661, A16 (2022)

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 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

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

The Next Generation X-ray Galaxy Survey with eROSITA

Authors: Antara R. Basu-Zych, Ann E. Hornschemeier, Frank Haberl, Neven Vulic, Jorn Wilms, Andreas Zezas, Konstantinos Kovlakas, Andrew Ptak, Thomas Dauser

Abstract: eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientifi… ▽ More eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e., not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50-200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect $\gtrsim$ 15 000 galaxies (3 $σ$ significance) at 50-200 Mpc, which is ~100X more normal galaxies than detected in any X-ray survey to date. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: 18 pages, 15 figures. 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/

The Environment of Lyman Break Analogues (ELBA) survey: Star-forming galaxies in small groups

Authors: L. Santana-Silva, T. S. Gonçalves, A. Basu-Zych, M. Soares-Santos, K. Menéndez-Delmestre, A. Drlica-Wagner, L. Riguccini, N. P. Kuropatkin, B. Yanny, R. T. Eufrasio

Abstract: The Environment of Lyman Break Analogues (ELBA) survey is an imaging survey of 33 $deg^{2}$ of the southern sky. The survey was observed in {\it u}, {\it g}, {\it r}, and {\it i} bands with the Dark Energy Camera (DECam) on the Blanco telescope. The main goal of this project is to investigate the environment of Lyman break analogues (LBAs), low-redshift (z $\sim $0.2) galaxies that are remarkably… ▽ More The Environment of Lyman Break Analogues (ELBA) survey is an imaging survey of 33 $deg^{2}$ of the southern sky. The survey was observed in {\it u}, {\it g}, {\it r}, and {\it i} bands with the Dark Energy Camera (DECam) on the Blanco telescope. The main goal of this project is to investigate the environment of Lyman break analogues (LBAs), low-redshift (z $\sim $0.2) galaxies that are remarkably similar to typical star-forming galaxies at z $\sim$ 3. We explore whether the environment has any influence on the observed properties of these galaxies, providing valuable insight on the formation and evolution of galaxies over cosmic time. Using the Nearest Neighbour method, we measure the local density of each object ranging from small to large scales (clusters of galaxies). Comparing the environment around LBAs with that of the general galaxy population in the field, we conclude that LBAs, on average, populate denser regions at small scales ($\sim$ $1.5Mpc$), but are located in similar environment to other star-forming galaxies at larger scales ($\sim$ $3.0 Mpc$). This offers evidence that nearby encounters such as mergers may influence the star formation activity in LBAs, before infall onto larger galaxy clusters. We interpret this an indication of galaxy preprocessing, in agreement with theoretical expectations for galaxies at z $\sim$ 2 -3 where the gravitational interactions are more intense in early formation processes of this objects △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 12 pages, 12 figures, submitted to MNRAS

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/)

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

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

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 deepest X-ray view of high-redshift galaxies: constraints on low-rate black-hole accretion

Authors: Fabio Vito, Roberto Gilli, Cristian Vignali, William N. Brandt, Andrea Comastri, Guang Yang, Bret D. Lehmer, Bin Luo, Antara Basu-Zych, Franz E. Bauer, Nico Cappelluti, Anton Koekemoer, Vincenzo Mainieri, Maurizio Paolillo, Piero Ranalli, Ohad Shemmer, Jonathan Trump, Junxian Wang, Yongquan Xue

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

Submitted 8 August, 2016; originally announced August 2016.

Comments: Submitted to MNRAS. Accepted 2016 August 5. Received 2016 August 5; in original form 2016 April 18

The Evolution of Normal Galaxy X-ray Emission Through Cosmic History: Constraints from the 6 Ms Chandra Deep Field-South

Authors: B. D. Lehmer, A. R. Basu-Zych, S. Mineo, W. N. Brandt, R. T. Eufrasio, T. Fragos, A. E. Hornschemeier, B. Luo, Y. Q. Xue, F. E. Bauer, M. Gilfanov, P. Ranalli, D. P. Schneider, O. Shemmer, P. Tozzi, J. R. Trump, C. Vignali, J. -X. Wang, M. Yukita, A. Zezas

Abstract: We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small c… ▽ More We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Comments: Accepted for publication in ApJ (23 pages, 13 figures, 5 tables)

Exploring the Overabundance of ULXs in Metal- and Dust-poor Local Lyman Break Analogs

Authors: Antara R. Basu-Zych, Bret Lehmer, Tassos Fragos, Ann Hornschemeier, Mihoko Yukita, Andreas Zezas, Andy Ptak

Abstract: We have studied high mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z>2) Lyman break galaxies, and within the larger sample of Lyman break analogs (LBAs) are sufficiently nearby (<87 Mpc) to be spatially-resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the… ▽ More We have studied high mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z>2) Lyman break galaxies, and within the larger sample of Lyman break analogs (LBAs) are sufficiently nearby (<87 Mpc) to be spatially-resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12+log[O/H]= 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. We have performed an in-depth study of the only two LBAs that have spatially-resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L$_X>10^{39}$ erg/s; ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on 8 detected ULXs. Comparing with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al. (2012), Haro 11 and VV 114 host ~4 times more L$_X>10^{40}$ erg/s sources than expected given their SFRs. We simulate the effects of source blending from crowded lower luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF shapes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright end slope than the standard XLF. If we conservatively assume that the brightest X-ray source from each galaxy is powered by an AGN rather than a HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF. [Abridged] △ Less

Submitted 4 January, 2016; originally announced January 2016.

Comments: Accepted for publication in ApJ (13 pages, 10 figures)

Molecular gas properties of UV-luminous star-forming galaxies at low redshift

Authors: T. S. Gonçalves, A. Basu-Zych, R. A. Overzier, L. Pérez, D. C. Martin

Abstract: Lyman break analogues (LBAs) are a population of star-forming galaxies at low redshift (z ~ 0.2) selected in the ultraviolet (UV). These objects present higher star formation rates and lower dust extinction than other galaxies with similar masses and luminosities in the local universe. In this work we present results from a survey with the Combined Array for Research in Millimetre-wave Astronomy (… ▽ More Lyman break analogues (LBAs) are a population of star-forming galaxies at low redshift (z ~ 0.2) selected in the ultraviolet (UV). These objects present higher star formation rates and lower dust extinction than other galaxies with similar masses and luminosities in the local universe. In this work we present results from a survey with the Combined Array for Research in Millimetre-wave Astronomy (CARMA) to detect CO(1-0) emission in LBAs, in order to analyse the properties of the molecular gas in these galaxies. Our results show that LBAs follow the same Schmidt-Kennicutt law as local galaxies. On the other hand, they have higher gas fractions (up to 66%) and faster gas depletion time-scales (below 1 Gyr). These characteristics render these objects more akin to high-redshift star-forming galaxies. We conclude that LBAs are a great nearby laboratory for studying the cold interstellar medium in low-metallicity, UV-luminous compact star-forming galaxies. △ Less

Submitted 4 June, 2014; originally announced June 2014.

Comments: 13 pages, 8 figures; accepted for publication in MNRAS

The X-ray Luminosity Functions of Field Low Mass X-ray Binaries in Early-Type Galaxies: Evidence for a Stellar Age Dependence

Authors: B. D Lehmer, M. Berkeley, A. Zezas, D. M. Alexander, A. Basu-Zych, F. E. Bauer, W. N. Brandt, T. Fragos, A. E. Hornschemeier, V. Kalogera, A. Ptak, G. R. Sivakoff, P. Tzanavaris, M. Yukita

Abstract: We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 3379, and 3384. The luminosity-weighted stellar ages of our sample span 3-1… ▽ More We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 3379, and 3384. The luminosity-weighted stellar ages of our sample span 3-10 Gyr. X-ray binary population synthesis models predict that the field LMXBs associated with younger stellar populations should be more numerous and luminous per unit stellar mass than older populations due to the evolution of LMXB donor star masses. Crucially, the combination of deep Chandra and HST observations allows us to test directly this prediction by identifying and removing counterparts to X-ray point sources that are unrelated to the field LMXB populations, including LMXBs that are formed dynamically in globular clusters, Galactic stars, and background AGN/galaxies. We find that the "young" early-type galaxy NGC 3384 (~2-5 Gyr) has an excess of luminous field LMXBs (L_X > (5-10) x 10^37 erg/s) per unit K-band luminosity (L_K; a proxy for stellar mass) than the "old" early-type galaxies NGC 3115 and 3379 (~8-10 Gyr), which results in a factor of ~2-3 excess of LX/LK for NGC 3384. This result is consistent with the X-ray binary population synthesis model predictions; however, our small galaxy sample size does not allow us to draw definitive conclusions on the evolution field LMXBs in general. We discuss how future surveys of larger galaxy samples that combine deep Chandra and HST data could provide a powerful new benchmark for calibrating X-ray binary population synthesis models. △ Less

Submitted 8 May, 2014; originally announced May 2014.

Comments: Accepted for publication in ApJ (16 pages, 9 figures, 6 tables)

Modeling X-ray binary evolution in normal galaxies: Insights from SINGS

Authors: P. Tzanavaris, T. Fragos, M. Tremmel, L. Jenkins, A. Zezas, B. D. Lehmer, A. Hornschemeier, V. Kalogera, A. Ptak, A. R. Basu-Zych

Abstract: We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthet… ▽ More We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration ~=0.1, and a 50% uniform -- 50% "twins" initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L_x - star formation rate and L_x - stellar mass relations of Lehmer et al 2010. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift. △ Less

Submitted 21 June, 2013; originally announced June 2013.

Comments: Accepted by ApJ

Journal ref: ApJ, 2013, 774, 136

Energy Feedback from X-ray Binaries in the Early Universe

Authors: Tassos Fragos, Bret D. Lehmer, Smadar Naoz, Andreas Zezas, Antara R. Basu-Zych

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

Submitted 19 September, 2013; v1 submitted 6 June, 2013; originally announced June 2013.

Comments: Accepted for publication to ApJ Letters, 6 pages, 2 figures, 2 table. Significant changes to figure 2b

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

Authors: Antara R. Basu-Zych, Bret D. Lehmer, Ann E. Hornschemeier, Thiago S. Goncalves, Tassos Fragos, Tim Heckman, Roderik A. Overzier, Andrew F. Ptak, David Schiminovich

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

Submitted 11 June, 2013; v1 submitted 4 June, 2013; originally announced June 2013.

Comments: 11 pages, 7 figures, Submitted to ApJ (references updated in v2)

Modeling the Redshift Evolution of the Normal Galaxy X-ray Luminosity Function

Authors: M. Tremmel, T. Fragos, B. D. Lehmer, P. Tzanavaris, K. Belczynski, V. Kalogera, A. R. Basu-Zych, W. M. Farr, A. Hornschemeier, L. Jenkins, A. Ptak, A. Zezas

Abstract: Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and evolution of X-ray binaries on very long timescales. In this paper, we present results from large-scale population synthesis models of binary populations in galaxies from z = 0 to 20. We use as input into our model… ▽ More Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and evolution of X-ray binaries on very long timescales. In this paper, we present results from large-scale population synthesis models of binary populations in galaxies from z = 0 to 20. We use as input into our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history (SFH) and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the parameters known from previous studies to affect the evolution of XRBs. We use our models and observationally derived prescriptions for hot gas emission to create theoretical galaxy X-ray luminosity functions (XLFs) for several redshift bins. Models with low CE efficiencies, a 50% twins mass ratio distribution, a steeper IMF exponent, and high stellar wind mass loss rates best match observational results from Tzanavaris & Georgantopoulos (2008), though they significantly underproduce bright early-type and very bright (Lx > 10d41) late-type galaxies. These discrepancies are likely caused by uncertainties in hot gas emission and SFHs, AGN contamination, and a lack of dynamically formed Low-mass XRBs. In our highest likelihood models, we find that hot gas emission dominates the emission for most bright galaxies. We also find that the evolution of the normal galaxy X-ray luminosity density out to z = 4 is driven largely by XRBs in galaxies with X-ray luminosities between 10d40 and 10d41 erg/s. △ Less

Submitted 1 February, 2013; v1 submitted 26 October, 2012; originally announced October 2012.

Comments: Accepted into ApJ, 17 pages, 3 tables, 7 figures. Text updated to address referee's comments

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

Authors: Antara R. Basu-Zych, Bret D. Lehmer, Ann E. Hornschemeier, Rychard J. Bouwens, Tassos Fragos, Pascal A. Oesch, Krzysztof Belczynski, W. N. Brandt, Vassiliki Kalogera, Bin Luo, Neal Miller, James R. Mullaney, Panayiotis Tzanavaris, Yongquan Xue, Andreas Zezas

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

Submitted 11 October, 2012; originally announced October 2012.

Comments: 16 pages, 9 figures. Accepted for publication in ApJ

X-ray Binary Evolution Across Cosmic Time

Authors: Tassos Fragos, Bret Lehmer, Michael Tremmel, Panayiotis Tzanavaris, Antara Basu-Zych, Krzysztof Belczynski, Ann Hornschemeier, Leigh Jenkins, Vassiliki Kalogera, Andrew Ptak, Andreas Zezas

Abstract: High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the… ▽ More High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray binary populations with properties such as star-formation rate and stellar mass, and the evolution of these relations with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) from X-ray binaries in our Universe today is dominated by low-mass X-ray binaries, and it is only at z>2.5 that high-mass X-ray binaries become dominant. We also find that there is a delay of ~1.1 Gyr between the peak of X-ray emissivity from low-mass Xray binaries (at z~2.1) and the peak of star-formation rate density (at z~3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z~3.9), happens ~0.8 Gyr before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe. △ Less

Submitted 20 November, 2012; v1 submitted 11 June, 2012; originally announced June 2012.

Comments: 14 pages, 8 figures, Accepted for publication in ApJ. Main changes from last revision: updated Figure 2, added section 4.3

A search for AGN in the most extreme UV-selected starbursts using the European VLBI Network

Authors: R. Alexandroff, R. A. Overzier, Z. Paragi, A. Basu-Zych, T. Heckman, G. Kauffmann, S. Bourke, A. Lobanov, A. Ptak, D. Schiminovich

Abstract: We have used the European VLBI Network (EVN) to observe a sample of Lyman Break Analogs (LBAs), nearby (z < 0.3) galaxies with properties similar to the more distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed-back relationship between black holes (BHs) and their host galaxies. Previous VLA observations have shown that the kpc-scale radio emission from LBAs is dominated… ▽ More We have used the European VLBI Network (EVN) to observe a sample of Lyman Break Analogs (LBAs), nearby (z < 0.3) galaxies with properties similar to the more distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed-back relationship between black holes (BHs) and their host galaxies. Previous VLA observations have shown that the kpc-scale radio emission from LBAs is dominated by starbursts. The main targets of this VLBI experiment were selected because they possessed emission-line properties between starbursts and Type 2 (obscured) AGN. Eight targets (three star-forming LBAs, four composite LBAs, and one Type 1 AGN) were observed at 5 GHz, four of which were also observed at 1.7 GHz. One star-forming LBA and one composite LBA were detected above 5 σat 1.7 GHz (only), while the AGN was detected at 5 GHz. In both LBAs, the radio luminosity (LR) exceeds that expected from supernovae (remnants) based on a comparison with Arp220, Arp229A and Mrk273, by factors of 2 - 8. The composite LBA exhibits a compact core emitting around 10% of the VLA flux density. The high Tb of 3.5E7 K and excess core L_R with respect to the L_R/L_X relation of radio-quiet AGN indicate that this LBA possesses an obscured AGN (MBH ~ E5-E7 M_sun). While weak AGN may co-exist with the starbursts as shown in at least one of the LBAs, their contribution to the total radio flux is fairly minimal. Our results show that the detection of such weak AGN presents a challenge at radio, X-ray and optical emission-line wavelengths at z ~ 0.2, indicating the great difficulties that need to be overcome in order to study similar processes at high redshift when these types of galaxies were common. △ Less

Submitted 21 March, 2012; v1 submitted 20 March, 2012; originally announced March 2012.

Comments: 10 pages, 4 figures, accepted for publication in MNRAS

A Search for Lyman Break Galaxies in the CDF-S Using Swift UVOT

Authors: Antara R. Basu-Zych, Ann E. Hornschemeier, Erik A. Hoversten, Bret Lehmer, Caryl Gronwall

Abstract: (Abridged) While the Swift satellite is primarily designed to study gamma-ray bursts, its ultraviolet and optical imaging and spectroscopy capabilities are also being used for a variety of scientific programs. In this study, we use the UV/Optical Telescope (UVOT) instrument aboard Swift to discover 0.5<z<2 Lyman break galaxies (LBGs). UVOT has covered ~266 arcmin^2 at >60ks exposure time, achievin… ▽ More (Abridged) While the Swift satellite is primarily designed to study gamma-ray bursts, its ultraviolet and optical imaging and spectroscopy capabilities are also being used for a variety of scientific programs. In this study, we use the UV/Optical Telescope (UVOT) instrument aboard Swift to discover 0.5<z<2 Lyman break galaxies (LBGs). UVOT has covered ~266 arcmin^2 at >60ks exposure time, achieving a limiting magnitude of u<24.5, in the Chandra Deep Field South (CDF-S). Applying UVOT near-ultraviolet color selection, we select 50 UV-dropouts from this UVOT CDF-S data. We match the selected sources with available multiwavelength data from GOODS-South, MUSYC, and COMBO-17 to characterize the spectral energy distributions for these galaxies and determine stellar masses, star formation rates (SFRs), and dust attenuations. We compare these properties for LBGs selected in this paper versus z~3 LBGs and other CDF-S galaxies in the same redshift range (0.5<z<2), identified using photometric redshift techniques. The z~1 LBGs have slightly lower stellar masses compared to z~3 LBGs and slightly higher stellar masses compared to the z~1 CDF-S galaxies. Similarly, our sample of z~1 LBGs has SFRs (derived using both ultraviolet and infrared data, where available) nearly an order of magnitude lower than z~3 LBGs but slightly higher than the comparison z~1 sample of CDF-S galaxies. We find that our z~1 UV-dropouts have A_FUV higher than z~3 LBGs, but is similar to the distribution of dust attenuations in the other CDF-S galaxies. Using the GOODS-South multiwavelength catalog of galaxies, we simulate a larger and fainter sample of LBGs to compare their properties with those of the UVOT-selected LBG sample. We conclude that UVOT can be useful for finding and studying the bright end of 0.5<z<2.0 LBGs. △ Less

Submitted 19 July, 2011; originally announced July 2011.

Comments: 18 pages, 12 figures. Accepted by ApJ

Swift UV/Optical Telescope Imaging of Star Forming Regions in M81 and Holmberg IX

Authors: E. A. Hoversten, C. Gronwall, D. E. Vanden Berk, A. R. Basu-Zych, A. A. Breeveld, P. J. Brown, N. P. M. Kuin, M. J. Page, P. W. A. Roming, M. H. Siegel

Abstract: We present Swift UV/Optical Telescope (UVOT) imaging of the galaxies M81 and Holmberg IX. We combine UVOT imaging in three near ultraviolet (NUV) filters (uvw2: 1928 Å, uvm2: 2246 Å, and uvw1: 2600 Å) with ground based optical imaging from the Sloan Digital Sky Survey to constrain the stellar populations of both galaxies. Our analysis consists of three different methods. First we use the NUV imagi… ▽ More We present Swift UV/Optical Telescope (UVOT) imaging of the galaxies M81 and Holmberg IX. We combine UVOT imaging in three near ultraviolet (NUV) filters (uvw2: 1928 Å, uvm2: 2246 Å, and uvw1: 2600 Å) with ground based optical imaging from the Sloan Digital Sky Survey to constrain the stellar populations of both galaxies. Our analysis consists of three different methods. First we use the NUV imaging to identify UV star forming knots and then perform SED modeling on the UV/optical photometry of these sources. Second, we measure surface brightness profiles of the disk of M81 in the NUV and optical. Last we use SED fitting of individual pixels to map the properties of the two galaxies. In agreement with earlier studies we find evidence for a burst in star formation in both galaxies starting ~200 Myr ago coincident with the suggested time of an M81-M82 interaction. In line with theories of its origin as a tidal dwarf we find that the luminosity weighted age of Holmberg IX is a few hundred million years. Both galaxies are best fit by a Milky Way dust extinction law with a prominent 2175 Å bump. In addition, we describe a stacked median filter technique for modeling the diffuse background light within a galaxy, and a Markov chain method for cleaning segment maps generated by SExtractor. △ Less

Submitted 8 April, 2011; originally announced April 2011.

Comments: 21 pages, 16 figures, accepted to AJ

Lyman Break Analogs: Constraints on the Formation of Extreme Starbursts at Low and High Redshift

Authors: Thiago S. Gonçalves, Roderik Overzier, Antara Basu-Zych, D. Christopher Martin

Abstract: Lyman Break Analogs (LBAs), characterized by high far-UV luminosities and surface brightnesses as detected by GALEX, are intensely star-forming galaxies in the low-redshift universe ($z\sim 0.2$), with star formation rates reaching up to 50 times that of the Milky Way. These objects present metallicities, morphologies and other physical properties similar to higher redshift Lyman Break Galaxies (L… ▽ More Lyman Break Analogs (LBAs), characterized by high far-UV luminosities and surface brightnesses as detected by GALEX, are intensely star-forming galaxies in the low-redshift universe ($z\sim 0.2$), with star formation rates reaching up to 50 times that of the Milky Way. These objects present metallicities, morphologies and other physical properties similar to higher redshift Lyman Break Galaxies (LBGs), motivating the detailed study of LBAs as local laboratories of this high-redshift galaxy population. We present results from our recent integral-field spectroscopy survey of LBAs with Keck/OSIRIS, which shows that these galaxies have the same nebular gas kinematic properties as high-redshift LBGs. We argue that such kinematic studies alone are not an appropriate diagnostic to rule out merger events as the trigger for the observed starburst. Comparison between the kinematic analysis and morphological indices from HST imaging illustrates the difficulties of properly identifying (minor or major) merger events, with no clear correlation between the results using either of the two methods. Artificial redshifting of our data indicates that this problem becomes even worse at high redshift due to surface brightness dimming and resolution loss. Whether mergers could generate the observed kinematic properties is strongly dependent on gas fractions in these galaxies. We present preliminary results of a CARMA survey for LBAs and discuss the implications of the inferred molecular gas masses for formation models. △ Less

Submitted 13 February, 2011; v1 submitted 8 February, 2011; originally announced February 2011.

Comments: To appear in the proceedings of IAU Symposium 277, "Tracing the Ancestry of Galaxies on the Land of our Ancestors", eds. C. Carignan, K.C. Freeman, and F. Combes

Extreme Feedback and the Epoch of Reionization: Clues in the Local Universe

Authors: Timothy M. Heckman, Sanchayeeta Borthakur, Roderik Overzier, Guinevere Kauffmann, Antara Basu-Zych, Claus Leitherer, Ken Sembach, D. Chris Martin, R. Michael Rich, David Schiminovich, Mark Seibert

Abstract: The source responsible for reionizing the universe at z > 6 remains uncertain. While an energetically adequate population of star-forming galaxies may be in place, it is unknown whether a large enough fraction of their ionizing radiation can escape into the intergalactic medium. Attempts to measure this escape-fraction in intensely star-forming galaxies at lower redshifts have largely yielded uppe… ▽ More The source responsible for reionizing the universe at z > 6 remains uncertain. While an energetically adequate population of star-forming galaxies may be in place, it is unknown whether a large enough fraction of their ionizing radiation can escape into the intergalactic medium. Attempts to measure this escape-fraction in intensely star-forming galaxies at lower redshifts have largely yielded upper limits. In this paper we present new HST COS and archival FUSE far-UV spectroscopy of a sample of eleven Lyman Break Analogs (LBAs), a rare population of local galaxies that strongly resemble the high-z Lyman Break galaxies. We combine these data with SDSS optical spectra and Spitzer photometry. We also analyze archival FUSE observations of fifteen typical UV-bright local starbursts. We find evidence of small covering factors for optically-thick neutral gas in 3 cases. This is based on two independent pieces of evidence: a significant residual intensity in the cores of the strongest interstellar absorption-lines tracing neutral gas and a small ratio of extinction-corrected H-alpha to UV plus far-IR luminosities. These objects represent three of the four LBAs that contain a young, very compact (~100pc), and highly massive (~10^9 Mo) dominant central object (DCO). These three objects also differ from the other galaxies in showing a significant amount of blueshifted Ly-alpha emission, which may be related to the low covering factor of neutral gas. All four LBAs with DCOs in our sample show extremely high velocity outflows of interstellar gas, with line centroids blueshifted by about 700km/s and maximum outflow velocities reaching at least 1500km/s. We show that these properties are consistent with an outflow driven by a powerful starburst that is exceptionally compact. We speculate that such extreme feedback may be required to enable the escape of ionizing radiation from star forming galaxies. △ Less

Submitted 30 January, 2011; v1 submitted 21 January, 2011; originally announced January 2011.

Comments: Accepted for publication in ApJ. 14 pages, 8 figures, 2 tables

Dust Attenuation in UV-selected Starbursts at High Redshift and their Local Counterparts: Implications for the Cosmic Star Formation Rate Density

Authors: Roderik Overzier, Tim Heckman, Jing Wang, Lee Armus, Veronique Buat, Justin Howell, Gerhardt Meurer, Mark Seibert, Brian Siana, Antara Basu-Zych, Stéphane Charlot, Thiago S. Gonçalves, D. Christopher Martin, James D. Neill, R. Michael Rich, Samir Salim, David Schiminovich

Abstract: We present a new analysis of the dust obscuration in starburst galaxies at low and high redshift. This study is motivated by our unique sample of the most extreme UV-selected starburst galaxies in the nearby universe (z<0.3), found to be good analogs of high-redshift Lyman Break Galaxies (LBGs) in most of their physical properties. We find that the dust properties of the Lyman Break Analogs (LBAs)… ▽ More We present a new analysis of the dust obscuration in starburst galaxies at low and high redshift. This study is motivated by our unique sample of the most extreme UV-selected starburst galaxies in the nearby universe (z<0.3), found to be good analogs of high-redshift Lyman Break Galaxies (LBGs) in most of their physical properties. We find that the dust properties of the Lyman Break Analogs (LBAs) are consistent with the relation derived previously by Meurer et al. (M99) that is commonly used to dust-correct star formation rate measurements at a very wide range of redshifts. We directly compare our results with high redshift samples (LBGs, BzK, and sub-mm galaxies at z=2-3) having IR data either from Spitzer or Herschel. The attenuation in typical LBGs at z=2-3 and LBAs is very similar. Because LBAs are much better analogs to LBGs compared to previous local star-forming samples, including M99, the practice of dust-correcting the SFRs of high redshift galaxies based on the local calibration is now placed on a much more solid ground. We illustrate the importance of this result by showing how the locally calibrated relation between UV measurements and extinction is used to estimate the integrated, dust-corrected star formation rate density at z=2-6. △ Less

Submitted 28 November, 2010; originally announced November 2010.

Comments: Accepted for publication in ApJ Letters (6 pages, 4 figures, 1 table)

The kinematics of ionized gas in Lyman-Break Analogs at z ~ 0.2

Authors: Thiago S. Gonçalves, Antara Basu-Zych, Roderik Overzier, D. Christopher Martin, David R. Law, David Schiminovich, Ted K. Wyder, Ryan Mallery, R. Michael Rich, Timothy H. Heckman

Abstract: We present results for 19 "Lyman Break Analogs" (LBAs) observed with Keck/OSIRIS with an AO-assisted spatial resolution of less than 200 pc. We detect satellites/companions, diffuse emission and velocity shear, all with high signal-to-noise ratios. These galaxies present remarkably high velocity dispersion along the line of sight(- 70 km s-1), much higher than standard star-forming spirals in the… ▽ More We present results for 19 "Lyman Break Analogs" (LBAs) observed with Keck/OSIRIS with an AO-assisted spatial resolution of less than 200 pc. We detect satellites/companions, diffuse emission and velocity shear, all with high signal-to-noise ratios. These galaxies present remarkably high velocity dispersion along the line of sight(- 70 km s-1), much higher than standard star-forming spirals in the low-redshift universe. We artificially redshift our data to z - 2.2 to allow for a direct comparison with observations of high-z LBGs and find striking similarities between both samples. This suggests that either similar physical processes are responsible for their observed properties, or, alternatively, that it is very difficult to distinguish between different mechanisms operating in the low versus high redshift starburst galaxies based on the available data. The comparison between morphologies in the UV/optical continuum and our kinemetry analysis often shows that neither is by itself sufficient to confirm or completely rule out the contribution from recent merger events. We find a correlation between the kinematic properties and stellar mass, in that more massive galaxies show stronger evidence for a disk-like structure. This suggests a co-evolutionary process between the stellar mass build-up and the formation of morphological and dynamical sub-structure within the galaxy. △ Less

Submitted 24 September, 2010; originally announced September 2010.

Comments: 18 pages; accepted for publication in ApJ

The GALEX Arecibo SDSS Survey II: The Star Formation Efficiency of Massive Galaxies

Authors: David Schiminovich, Barbara Catinella, Guinevere Kauffmann, Silvia Fabello, Jing Wang, Cameron Hummels, Jenna Lemonias, Sean M. Moran, Ronin Wu, Riccardo Giovanelli, Martha P. Haynes, Timothy M. Heckman, Antara R. Basu-Zych, Michael R. Blanton, Jarle Brinchmann, Tamas Budavari, Thiago Goncalves, Benjamin D. Johnson, Robert C. Kennicutt, Barry F. Madore, Christopher D. Martin, Michael R. Rich, Linda J. Tacconi, David A. Thilker, Vivienne Wild , et al. (1 additional authors not shown)

Abstract: We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star f… ▽ More We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged] △ Less

Submitted 28 June, 2010; originally announced June 2010.

Comments: 18 pages, 11 figures. Accepted for publication in MNRAS. GASS publications and released data can be found at http://www.mpa-garching.mpg.de/GASS/index.php

The GALEX Arecibo SDSS Survey. I. Gas Fraction Scaling Relations of Massive Galaxies and First Data Release

Authors: Barbara Catinella, David Schiminovich, Guinevere Kauffmann, Silvia Fabello, Jing Wang, Cameron Hummels, Jenna Lemonias, Sean M. Moran, Ronin Wu, Riccardo Giovanelli, Martha P. Haynes, Timothy M. Heckman, Antara R. Basu-Zych, Michael R. Blanton, Jarle Brinchmann, Tamás Budavári, Thiago Gonçalves, Benjamin D. Johnson, Robert C. Kennicutt, Barry F. Madore, Christopher D. Martin, Michael R. Rich, Linda J. Tacconi, David A. Thilker, Vivienne Wild , et al. (1 additional authors not shown)

Abstract: We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mas… ▽ More We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged] △ Less

Submitted 9 December, 2009; originally announced December 2009.

Comments: 25 pages, 12 figures. Accepted for publication in MNRAS. Version with high resolution figures available at http://www.mpa-garching.mpg.de/GASS/pubs.php

Morphologies of local Lyman break galaxy analogs II: A Comparison with galaxies at z=2-4 in ACS and WFC3 images of the Hubble Ultra Deep Field

Authors: R. A. Overzier, T. M. Heckman, D. Schiminovich, A. Basu-Zych, T. Goncalves, D. C. Martin, R. M. Rich

Abstract: Lyman break galaxies (LBGs) display a range in structures (from single/compact to clumpy/extended) that is different from typical local star-forming galaxies. Recently, we have introduced a sample of rare, nearby (z<0.3) starbursts that appear to be good analogs of LBGs. These "Lyman Break Analogs" (LBAs) provide an excellent training set for understanding starbursts at different redshifts. We p… ▽ More Lyman break galaxies (LBGs) display a range in structures (from single/compact to clumpy/extended) that is different from typical local star-forming galaxies. Recently, we have introduced a sample of rare, nearby (z<0.3) starbursts that appear to be good analogs of LBGs. These "Lyman Break Analogs" (LBAs) provide an excellent training set for understanding starbursts at different redshifts. We present an application of this by comparing the rest-frame UV/optical morphologies of 30 LBAs with those of sBzK galaxies at z~2, and LBGs at z~3-4 in the HUDF. The UV/optical colors and sizes of LBAs and LBGs are very similar, while the BzK galaxies are somewhat redder and larger. There is significant overlap between the morphologies (G, C, A and M_20) of the local and high-z samples, although the latter are somewhat less concentrated and clumpier. We find that in the majority of LBAs the starbursts appear to be triggered by interactions/mergers. When the images of the LBAs are degraded to the same sensitivity and resolution as the images of LBGs and BzK galaxies, these relatively faint asymmetric features are no longer detectable. This effect is particularly severe in the rest-frame UV. It has been suggested that high-z galaxies experience intense bursts unlike anything seen locally, possibly due to cold flows and instabilities. In part, this is based on the fact that the majority (~70%) of LBGs do not show morphological signatures of mergers. Our results suggest that this evidence is insufficient, since a large fraction of such signatures would likely have been missed in current observations of z>2 galaxies. This leaves open the possibility that clumpy accretion and mergers remain important in driving the evolution of these starbursts, together with rapid gas accretion through other means. △ Less

Submitted 24 December, 2009; v1 submitted 6 November, 2009; originally announced November 2009.

Comments: ApJ, In Press (14 pages, 7 figures; minor changes since v1). For background material, see http://www.mpa-garching.mpg.de/~overzier/index.html

Journal ref: Astrophys.J.710:979-991,2010

Local Lyman Break Galaxy Analogs: The Impact of Massive Star-forming Clumps on the Interstellar Medium and the Global Structure of Young, Forming Galaxies

Authors: R. A. Overzier, T. M. Heckman, C. Tremonti, L. Armus, A. Basu-Zych, T. Goncalves, R. M. Rich, D. C. Martin, A. Ptak, D. Schiminovich, H. C. Ford, B. Madore, M. Seibert

Abstract: We present HST UV/optical imaging, Spitzer mid-IR photometry, and optical spectroscopy of a sample of 30 low-redshift (z=0.1-0.3) galaxies chosen from SDSS/GALEX to be accurate local analogs of the high-z Lyman Break Galaxies. The Lyman Break Analogs (LBAs) are similar in mass, metallicity, dust, SFR, size and gas velocity dispersion, thus enabling a detailed investigation of processes that are… ▽ More We present HST UV/optical imaging, Spitzer mid-IR photometry, and optical spectroscopy of a sample of 30 low-redshift (z=0.1-0.3) galaxies chosen from SDSS/GALEX to be accurate local analogs of the high-z Lyman Break Galaxies. The Lyman Break Analogs (LBAs) are similar in mass, metallicity, dust, SFR, size and gas velocity dispersion, thus enabling a detailed investigation of processes that are important at high-z. The optical emission line properties of LBAs are also similar to those of LBGs, indicating comparable conditions in their ISM. In the UV, LBAs are characterized by complexes of massive star-forming "clumps", while in the optical they most often show evidence for (post-)mergers/interactions. In 6 cases, we find an extremely massive (>10^9 Msun) compact (R~100 pc) dominant central object (DCO). The DCOs are preferentially found in LBAs with the highest mid-IR luminosities and correspondingly high SFRs (15-100 Msun/yr). We show that the massive SF clumps (including the DCOs) have masses much larger than the nuclear super star clusters seen in normal late type galaxies. However, the DCOs have masses, sizes, and densities similar to the excess-light/central-cusps seen in typical elliptical galaxies with masses similar to the LBA galaxies. We suggest that the DCOs form in present-day examples of the dissipative mergers at high redshift that are believed to have produced the central-cusps in local ellipticals. More generally, the properties of the LBAs are consistent with the idea that instabilities in a gas-rich disk lead to very massive star-forming clumps that eventually coalesce to form a spheroid. We speculate that the DCOs are too young at present to be growing a supermassive black hole because they are still in a supernova-dominated outflow phase. △ Less

Submitted 7 October, 2009; originally announced October 2009.

Comments: The Astrophysical Journal, In Press (22 pages, 16 figures). For the full version with high-resolution colour figures, see: http://www.mpa-garching.mpg.de/~overzier/Overzier_LBApaper09.pdf

Journal ref: Astrophys.J.706:203-222,2009

An OSIRIS study of the gas kinematics in a sample of UV-selected galaxies: Evidence of "Hot and Bothered" starbursts in the local Universe

Authors: Antara R. Basu-Zych, Thiago S. Goncalves, Roderik Overzier, David R. Law, David Schiminovich, Tim Heckman, Chris Martin, Ted Wyder, Matt O'Dowd

Abstract: We present data from Integral Field Spectroscopy for 3 supercompact UV-Luminous Galaxies (ScUVLGs). As nearby (z~0.2), compact (R_50~1-2 kpc), bright Paschen-alpha sources, with unusually high star formation rates (SFR=3-100 M_sun/yr), ScUVLGs are an ideal population for studying detailed kinematics and dynamics in actively star-forming galaxies. In addition, ScUVLGs appear to be excellent analo… ▽ More We present data from Integral Field Spectroscopy for 3 supercompact UV-Luminous Galaxies (ScUVLGs). As nearby (z~0.2), compact (R_50~1-2 kpc), bright Paschen-alpha sources, with unusually high star formation rates (SFR=3-100 M_sun/yr), ScUVLGs are an ideal population for studying detailed kinematics and dynamics in actively star-forming galaxies. In addition, ScUVLGs appear to be excellent analogs to high redshift Lyman Break Galaxies (LBGs) and our results may offer additional insight into the dynamics of LBGs. Previous work by our team has shown that the morphologies of these galaxies exhibit tidal features and companions, and in this study we find that the dynamics of ScUVLGs are dominated by disturbed kinematics of the emission line gas-- suggestive that these galaxies have undergone recent feedback, interactions or mergers. While 2 of the 3 galaxies do display rotation, v/sigma < 1 -- suggesting dispersion dominated kinematics rather than smooth rotation. We also simulate how these observations would appear at z~2. Lower resolution and loss of low surface brightness features causes some apparent discrepancies between the low-z (observed) and high-z (simulated) interpretations and quantitatively gives different values for v/sigma, yet simulations of these low-z analogs manage to detect the brightest regions well and resemble actual high-z observations of LBGs. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: 4 pages, 4 figures (bitmapped), accepted for publication in ApJL

Journal ref: Astrophys.J.699:L118-L124,2009

Studying Large and Small Scale Environments of Ultraviolet Luminous Galaxies

Authors: Antara R. Basu-Zych, David Schiminovich, Sebastien Heinis, Roderik Overzier, Tim Heckman, Michel Zamojski, Olivier Ilbert, Anton M. Koekemoer, Tom A. Barlow, Luciana Bianchi, Tim Conrow, Jose Donas, Karl G. Forster, Peter G. Friedman, Young-Wook Lee, Barry F. Madore, D. Christopher Martin, Bruno Milliard, Patrick Morrissey, Susan G. Neff, R. Michael Rich Samir Salim, Mark Seibert, Todd A. Small, Alex S. Szalay, Ted K. Wyder , et al. (1 additional authors not shown)

Abstract: Studying the environments of 0.4<z<1.2 UV-selected galaxies, as examples of extreme star-forming galaxies (with star formation rates in the range of 3-30 M_sol/yr), we explore the relationship between high rates of star-formation, host halo mass and pair fractions. We study the large-scale and small-scale environments of local Ultraviolet Luminous Galaxies (UVLGs) by measuring angular correlatio… ▽ More Studying the environments of 0.4<z<1.2 UV-selected galaxies, as examples of extreme star-forming galaxies (with star formation rates in the range of 3-30 M_sol/yr), we explore the relationship between high rates of star-formation, host halo mass and pair fractions. We study the large-scale and small-scale environments of local Ultraviolet Luminous Galaxies (UVLGs) by measuring angular correlation functions. We cross-correlate these systems with other galaxy samples: a volume-limited sample (ALL), a Blue Luminous Galaxy sample (BLG) and a Luminous Red Galaxy sample (LRG). We determine the UVLG comoving correlation length to be r_0=4.8(+11.6/-2.4) h^-1 Mpc at <z> =1.0, which is unable to constrain the halo mass for this sample. However, we find that UVLGs form close (separation < 30 kpc) pairs with the ALL sample, but do not frequently form pairs with LRGs. A rare subset of UVLGs, those with the highest FUV surface brightnesses, are believed to be local analogs of high redshift Lyman Break Galaxies (LBGs) and are called Lyman Break Analogs (LBAs). LBGs and LBAs share similar characteristics (i.e., color, size, surface brightness, specific star formation rates, metallicities, and dust content). Recent HST images of z~0.2 LBAs show disturbed morphologies, signs of mergers and interactions. UVLGs may be influenced by interactions with other galaxies and we discuss this result in terms of other high star-forming, merging systems. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: 35 pages, 10 Figures, accepted for publication in ApJ

Journal ref: Astrophys.J.699:1307-1320,2009

The Chandra Deep Protocluster Survey: Ly-alpha Blobs are powered by heating, not cooling

Authors: J. E. Geach, D. M. Alexander, B. D. Lehmer, Ian Smail, Y. Matsuda, S. C. Chapman, C. A. Scharf, R. J. Ivison, M. Volonteri, T. Yamada, A. W. Blain, R. G. Bower, F. E. Bauer, A. Basu-Zych

Abstract: We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices imply… ▽ More We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices implying obscuring column densities of N_H > 10^23 cm^-2. The AGN fraction should be considered a lower limit, since several more LABs not detected with Chandra show AGN signatures in their mid-infrared emission. We show that the UV luminosities of the AGN are easily capable of powering the extended Ly-alpha emission via photo-ionization alone. When combined with the UV flux from a starburst component, and energy deposited by mechanical feedback, we demonstrate that 'heating' by a central source, rather than gravitational cooling is the most likely power source of LABs. We argue that all LABs could be powered in this manner, but that the luminous host galaxies are often just below the sensitivity limits of current instrumentation, or are heavily obscured. No individual LABs show evidence for extended X-ray emission, and a stack equivalent to a >9Ms exposure of an average LAB also yields no statistical detection of a diffuse X-ray component. The resulting diffuse X-ray/Ly-alpha luminosity limit implies there is no hot (T>10^7 K) gas component in these halos, and also rules out inverse Compton scattering of cosmic microwave background photons, or local far-infrared photons, as a viable power source for LABs. △ Less

Submitted 23 June, 2009; v1 submitted 3 April, 2009; originally announced April 2009.

Comments: 9 pages, 4 figures. Accepted for publication in the Astrophysical Journal. This version updated to accepted version

Journal ref: Astrophys.J.700:1-9,2009

The Chandra Deep Protocluster Survey: Evidence for an Enhancement in the AGN Activity in the SSA22 Protocluster at z = 3.09

Authors: B. D. Lehmer, D. M. Alexander, J. E. Geach, Ian Smail, A. Basu-Zych, F. E. Bauer, S. C. Chapman, Y. Matsuda, C. A. Scharf, M. Volonteri, T. Yamada

Abstract: We present results from a new ultra-deep 400 ks Chandra observation of the SSA22 protocluster at z = 3.09. We have studied the X-ray properties of 234 z ~ 3 Lyman break galaxies (LBGs; protocluster and field) and 158 z = 3.09 Ly-alpha emitters (LAEs) in SSA22 to measure the influence of the high-density protocluster environment on the accretion activity of supermassive black holes (SMBHs) in the… ▽ More We present results from a new ultra-deep 400 ks Chandra observation of the SSA22 protocluster at z = 3.09. We have studied the X-ray properties of 234 z ~ 3 Lyman break galaxies (LBGs; protocluster and field) and 158 z = 3.09 Ly-alpha emitters (LAEs) in SSA22 to measure the influence of the high-density protocluster environment on the accretion activity of supermassive black holes (SMBHs) in these UV-selected star forming populations. We detect individually X-ray emission from active galactic nuclei (AGNs) in six LBGs and five LAEs; due to small overlap between the LBG and LAE source population, ten of these sources are unique. At least six and potentially eight of these sources are members of the protocluster. These sources have rest-frame 8-32 keV luminosities in the range of L_8-32 keV = (3-50) X 10^{43} ergs/s and an average observed-frame 2-8 keV to 0.5-2 keV band-ratio of ~0.8 (mean effective photon index of Gamma_eff = 1.1), suggesting significant absorption columns of N_H > 10^{22}-10^{24} cm^{-2}. We find that the fraction of LBGs and LAEs in the z = 3.09 protocluster harboring an AGN with L_8-32 keV > 3 X 10^{43} ergs/s is 9.5^{+12.7}_{-6.1}% and 5.1^{+6.8}_{-3.3}%, respectively. These AGN fractions are somewhat larger (by a mean factor of 6.1^{+10.3}_{-3.6}; significant at the 95% confidence level) than z ~ 3 sources found in lower-density "field" environments. Theoretical models imply that these results may be due to the presence of more actively growing and/or massive SMBHs in LBGs and LAEs within the protocluster compared to the field. Such a result is expected in a scenario where enhanced merger activity in the protocluster drives accelerated galaxy and SMBH growth at z > 2-3. (abridged) △ Less

Submitted 29 September, 2008; originally announced September 2008.

Comments: Accepted for publication in the Astrophysical Journal; 9 pages, 4 figures, and 2 tables

Journal ref: Astrophys.J.691:687-695,2009