Showing 1–22 of 22 results for author: Garofali, K

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

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

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

Submitted 5 December, 2023; originally announced December 2023.

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

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

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

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

Submitted 13 November, 2023; originally announced November 2023.

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

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

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

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

Submitted 8 November, 2023; originally announced November 2023.

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

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

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

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

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

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

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

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

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

Submitted 30 May, 2023; originally announced May 2023.

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

Multiwavelength Characterization of the High Mass X-ray Binary Population of M33

Authors: Margaret Lazzarini, Kyros Hinton, Cheyanne Shariat, Benjamin F. Williams, Kristen Garofali, Julianne J. Dalcanton, Meredith Durbin, Vallia Antoniou, Breanna Binder, Michael Eracleous, Neven Vulic, Jun Yang, Daniel Wik, Aria Gasca, Quetzalcoatl Kuauhtzin

Abstract: We present multi-wavelength characterization of 65 high mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalog, which covers the inner disk of M33 a… ▽ More We present multi-wavelength characterization of 65 high mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalog, which covers the inner disk of M33 at near infrared, optical, and near ultraviolet wavelengths. We perform spectral energy distribution (SED) fitting on multi-band photometry for each point source optical counterpart to measure its physical properties including mass, temperature, luminosity, and radius. We find that the majority of the HMXB companion star candidates are likely B-type main sequence stars, suggesting that the HMXB population of M33 is dominated by Be-XRBs, as is seen in other Local Group galaxies. We use spatially-resolved recent star formation history (SFH) maps of M33 to measure the age distribution of the HMXB candidate sample and the HMXB production rate for M33. We find a bimodal distribution for the HMXB production rate over the last 80 Myr, with a peak at $\sim$10 Myr and $\sim$40 Myr, which match theoretical formation timescales for the most massive HMXBs and Be X-ray binaries (Be-XRBs), respectively. We measure an HMXB production rate of 107$-$136 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 50 Myr and 150$-$199 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 80 Myr. For sources with compact object classifications from overlapping NuSTAR observations, we find a preference for giant/supergiant companion stars in BH-HMXBs and main sequence companion stars in neutron star HMXBs (NS-HMXBs). △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ, 37 pages, 6 figures, 8 tables

The Spatial Correlation of High Mass X-ray Binaries and Young Star Clusters in Nearby Star-Forming Galaxies

Authors: Breanna A. Binder, Ashley K. Anderson, Kristen Garofali, Margaret Lazzarini, Benjamin F. Williams

Abstract: We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we… ▽ More We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we find that the peak value of the spatial correlation function correlates strongly with the specific star formation rate of the host galaxy. We additionally perform an X-ray stacking analysis of 211 non-X-ray detected YSCs in M31 and 463 YSCs in M33. We do not detect excess X-ray emission at the stacked cluster locations down to 3$σ$ upper limits of $\sim10^{33}$ erg s$^{-1}$ (0.35-8 keV) in both galaxies, which strongly suggests that dynamical formation within YSCs is not a major HMXB formation channel. We interpret our results in the context of (1) the recent star formation histories of the galaxies, which may produce differences in the demographics of compact objects powering the HMXBs, and (2) the differences in natal kicks experienced by compact objects during formation, which can eject newly-formed HMXB from their birth clusters. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Accepted to MNRAS; 11 pages, 7 figures

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

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 Stellar Age Dependence of X-ray Emission from Normal Star-Forming Galaxies in the GOODS Fields

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

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

Submitted 6 December, 2021; originally announced December 2021.

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

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

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

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

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

Journal ref: ApJ 923 (2021) 26

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

Multiwavelength Study of the X-Ray Bright Supernova Remnant N300-S26 in NGC 300

Authors: Jacob Gross, Benjamin F. Williams, Thomas G. Pannuti, Breanna Binder, Kristen Garofali, Zachary G. Hanvey

Abstract: We present a multiwavelength examination of the supernova remnant (SNR) S26 in the nearby galaxy NGC 300 using data from Chandra X-ray Observatory, XMM-Newton X-ray Observatory, Hubble Space Telescope (HST), the Very Large Array, and the Australia Telescope Compact Array. We simultaneously fit all of the available X-ray data with a thermal plasma model and find a temperature of $0.77 \pm 0.13$ keV… ▽ More We present a multiwavelength examination of the supernova remnant (SNR) S26 in the nearby galaxy NGC 300 using data from Chandra X-ray Observatory, XMM-Newton X-ray Observatory, Hubble Space Telescope (HST), the Very Large Array, and the Australia Telescope Compact Array. We simultaneously fit all of the available X-ray data with a thermal plasma model and find a temperature of $0.77 \pm 0.13$ keV with a hydrogen column density of ($9.7^{+6.4}_{-4.8}$)$\times 10^{20}$ cm$^{-2}$. HST imaging allows us to measure a semimajor axis of $0.78 \pm 0.10$ arcsec ($7.5 \pm 1.0$ pc) and a semiminor axis of $0.69^{+0.14}_{-0.12}$ arcsec ($6.7^{+1.2}_{-1.4}$ pc). This precise size helps to constrain the age and velocity of the shock to be ($3.3^{+0.7}_{-0.6}$)$\times 10^{3}$ yr and $411^{+275}_{-122}$ km s$^{-1}$. We also fit photometry of the surrounding stars to infer the age and mass of the progenitor star to be $8 \pm 1$ Myr and $25^{+1}_{-5}$ M$_{\odot}$. Based on measured radio properties of the source and assuming equipartition, the estimated radio luminosity of $\sim 1.7 \times 10^{34}$ erg s$^{-1}$ over the $10^{8}-10^{11}$ Hz frequency range results in a minimum magnetic field associated with this SNR of $0.067$ mG and the minimum energy needed to power the observed synchrotron emission of $1.5 \times 10^{49}$ erg. The size and temperature of N300-S26 appear to be similar to the Galactic SNR G311.5-0.3 except that G311.5-0.3 has a significantly lower X-ray luminosity, is older, and has a slower shock velocity. △ Less

Submitted 11 June, 2019; originally announced June 2019.

Comments: 37 pages, 6 figures, 3 tables

Journal ref: ApJ, 877, 15 (2019)

The First Candidate Colliding-Wind Binary in M33

Authors: Kristen Garofali, Emily M. Levesque, Philip Massey, Benjamin F. Williams

Abstract: We present the detection of the first candidate colliding-wind binary (CWB) in M33, located in the giant H II region NGC 604. The source was first identified in archival {\it Chandra} imaging as a relatively soft X-ray point source, with the likely primary star determined from precise astrometric alignment between archival {\it Hubble Space Telescope} and {\it Chandra} imaging. The candidate prima… ▽ More We present the detection of the first candidate colliding-wind binary (CWB) in M33, located in the giant H II region NGC 604. The source was first identified in archival {\it Chandra} imaging as a relatively soft X-ray point source, with the likely primary star determined from precise astrometric alignment between archival {\it Hubble Space Telescope} and {\it Chandra} imaging. The candidate primary star in the CWB is classified for the first time in this work as a carbon-rich Wolf-Rayet star with a likely O star companion based on spectroscopy obtained from Gemini-North. We model the X-ray spectrum using {\it Chandra} and {\it XMM-Newton} observations, and find the CWB is well-fit as a $\sim$ 1 keV thermal plasma with a median unabsorbed luminosity in the 0.5--2.0 keV band of $L_{\rm X}$ $\sim$ 3 $\times$ 10$^{35}$ erg s$^{-1}$, making this source among the brightest of CWBs observed to date. We present a long term light curve for the candidate CWB from archival {\it Chandra} and {\it XMM-Newton} observations, and discuss the constraints placed on the binary by this light curve, as well as the X-ray luminosity at maximum. Finally, we compare this candidate CWB in M33 to other well-studied, bright CWBs in the Galaxy and Magellanic Clouds, such as $η$ Car. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: 21 pages, 8 figures, 4 tables, Accepted for publication in ApJ

Supernova Remnants in High Definition

Authors: Laura A. Lopez, Brian J. Williams, Samar Safi-Harb, Sangwook Park, Paul P. Plucinsky, David Pooley, Tea Temim, Katie Auchettl, Aya Bamba, Carles Badenes, Daniel Castro, Kristen Garofali, Denis Leahy, Patrick Slane, Jacco Vink, Benjamin F. Williams, J. Craig Wheeler

Abstract: Supernova remnants (SNRs) offer the means to study SN explosions, dynamics, and shocks at sub-parsec scales. X-ray observations probe the hot metals synthesized in the explosion and the TeV electrons accelerated by the shocks, and thus they are key to test recent, high-fidelity three-dimensional SN simulations. In this white paper, we discuss the major advances possible with X-ray spectro-imaging… ▽ More Supernova remnants (SNRs) offer the means to study SN explosions, dynamics, and shocks at sub-parsec scales. X-ray observations probe the hot metals synthesized in the explosion and the TeV electrons accelerated by the shocks, and thus they are key to test recent, high-fidelity three-dimensional SN simulations. In this white paper, we discuss the major advances possible with X-ray spectro-imaging at arcsecond scales, with a few eV spectral resolution and a large effective area. These capabilities would revolutionize SN science, offering a three-dimensional view of metals synthesized in explosions and enabling population studies of SNRs in Local Group galaxies. Moreover, this future X-ray mission could detect faint, narrow synchrotron filaments and shock precursors that will constrain the diffusive shock acceleration process. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures, Astro2020 Science White Paper

Formation Timescales for High-Mass X-ray Binaries in M33

Authors: Kristen Garofali, Benjamin F. Williams, Tristan Hillis, Karoline M. Gilbert, Andrew E. Dolphin, Michael Eracleous, Breanna Binder

Abstract: We have identified 55 candidate high-mass X-ray binaries (HMXBs) in M33 using available archival {\it HST} and {\it Chandra} imaging to find blue stars associated with X-ray positions. We use the {\it HST} photometric data to model the color-magnitude diagrams in the vicinity of each candidate HMXB to measure a resolved recent star formation history (SFH), and thus a formation timescale, or age fo… ▽ More We have identified 55 candidate high-mass X-ray binaries (HMXBs) in M33 using available archival {\it HST} and {\it Chandra} imaging to find blue stars associated with X-ray positions. We use the {\it HST} photometric data to model the color-magnitude diagrams in the vicinity of each candidate HMXB to measure a resolved recent star formation history (SFH), and thus a formation timescale, or age for the source. Taken together, the SFHs for all candidate HMXBs in M33 yield an age distribution that suggests preferred formation timescales for HMXBs in M33 of $<$ 5 Myr and $\sim$ 40 Myr after the initial star formation episode. The population at 40 Myr is seen in other Local Group galaxies, and can be attributed to a peak in formation efficiency of HMXBs with neutron stars as compact objects and B star secondary companions. This timescale is preferred as neutron stars should form in abundance from $\sim$ 8 M$_{\odot}$ core-collapse progenitors on these timescales, and B stars are shown observationally to be most actively losing mass around this time. The young population at $<$ 5 Myr has not be observed in other Local Group HMXB population studies, but may be attributed to a population of very massive progenitors forming black holes very early on. We discuss these results in the context of massive binary evolution, and the implications for compact object binaries and gravitational wave sources. △ Less

Submitted 18 June, 2018; originally announced June 2018.

Comments: 20 pages, 9 figures, 4 tables, Accepted for publication in MNRAS

Pre-MAP Search for Transiting Objects Orbiting White Dwarfs

Authors: Aislynn Wallach, Brett M. Morris, Doug Branton, Teagan O'Reilly, Brittany Platt, Ada Beale, Andrew Yetter, Katie Reil, Kristen Garofali, Eric Agol

Abstract: Metal pollution in white dwarf atmospheres may be the accreted remnants of planetary objects. After the discovery of disintegrating planetary objects transiting WD 1145+017 (Vanderburg et al. 2015), undergraduates in the University of Washington's Pre- Major in Astronomy Program (Pre-MAP) were inspired to collect photometry of the brightest white dwarfs to hunt for similar transiting objects aroun… ▽ More Metal pollution in white dwarf atmospheres may be the accreted remnants of planetary objects. After the discovery of disintegrating planetary objects transiting WD 1145+017 (Vanderburg et al. 2015), undergraduates in the University of Washington's Pre- Major in Astronomy Program (Pre-MAP) were inspired to collect photometry of the brightest white dwarfs to hunt for similar transiting objects around other metal-polluted white dwarfs. Prior surveys have yet to make a detection of a transiting planet orbiting a white dwarf, yet white dwarfs are still an attractive target for searches of small, rocky planetary material. Since a typical white dwarf is Earth-sized, transits of Earth-sized planets could have depths >50%, so even low S/N photometry has a chance at discovering transiting material. We identified bright, northern, metal-polluted white dwarfs in the SDSS DR10 white dwarf catalog, and observed five targets with the Astrophysical Research Consortium Small Aperture Telescope (ARCSAT) 0.5-meter telescope at Apache Point Observatory. The ARCSAT photometry had sufficient precision to detect Moon-sized objects or larger at short orbital periods, though no such planets were detected for these targets. We look forward to surveys which may find planets orbiting white dwarfs, such as NASA's TESS, ESA's PLATO, and the Evryscope. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: Accepted to RNAAS

Supernova Remnants in M33: X-ray Properties as Observed by XMM-Newton

Authors: Kristen Garofali, Benjamin F. Williams, Paul P. Plucinsky, Terrance J. Gaetz, Brian Wold, Frank Haberl, Knox S. Long, William P. Blair, Thomas G. Pannuti, P. Frank Winkler, Jacob Gross

Abstract: We have carried out a study of the X-ray properties of the supernova remnant (SNR) population in M33 with XMM-Newton, comprising deep observations of 8 fields in M33 covering all of the area within the D$_{25}$ contours, and with a typical luminosity of 7.1$\times$10$^{34}$ erg s$^{-1}$ (0.2-2.0 keV) . Here we report our work to characterize the X-ray properties of the previously identified SNRs i… ▽ More We have carried out a study of the X-ray properties of the supernova remnant (SNR) population in M33 with XMM-Newton, comprising deep observations of 8 fields in M33 covering all of the area within the D$_{25}$ contours, and with a typical luminosity of 7.1$\times$10$^{34}$ erg s$^{-1}$ (0.2-2.0 keV) . Here we report our work to characterize the X-ray properties of the previously identified SNRs in M33, as well as our search for new X-ray detected SNRs. With our deep observations and large field of view we have detected 105 SNRs at the 3$σ$ level, of which 54 SNRs are newly detected in X-rays, and three are newly discovered SNRs. Combining XMM-Newton data with deep Chandra survey data allows detailed spectral fitting of 15 SNRs, for which we have measured temperatures, ionization timescales, and individual abundances. This large sample of SNRs allows us to construct an X-ray luminosity function, and compare its shape to luminosity functions from host galaxies of differing metallicities and star formation rates to look for environmental effects on SNR properties. We conclude that while metallicity may play a role in SNR population characteristics, differing star formation histories on short timescales, and small-scale environmental effects appear to cause more significant differences between X-ray luminosity distributions. In addition, we analyze the X-ray detectability of SNRs, and find that in M33 SNRs with higher [SII]/H$α$ ratios, as well as those with smaller galactocentric distances, are more detectable in X-rays. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: 29 pages, 14 figures, 5 tables, accepted for publication in MNRAS

A Deep XMM-Newton Survey of M33: Point Source Catalog, Source Detection and Characterization of Overlapping Fields

Authors: Benjamin F. Williams, Brian Wold, Frank Haberl, Kristen Garofali, William P. Blair, Terrance J. Gaetz, K. D. Kuntz, Knox S. Long, Thomas G. Pannuti, Wolfgang Pietsch, Paul P. Plucinsky, P. Frank Winkler

Abstract: We have obtained a deep 8-field XMM-Newton mosaic of M33 covering the galaxy out to the D$_{25}$ isophote and beyond to a limiting 0.2--4.5 keV unabsorbed flux of 5$\times$10$^{-16}$ erg cm$^{-2}$ s$^{-1}$ (L${>}$4$\times$10$^{34}$ erg s$^{-1}$ at the distance of M33). These data allow complete coverage of the galaxy with high sensitivity to soft sources such as diffuse hot gas and supernova remna… ▽ More We have obtained a deep 8-field XMM-Newton mosaic of M33 covering the galaxy out to the D$_{25}$ isophote and beyond to a limiting 0.2--4.5 keV unabsorbed flux of 5$\times$10$^{-16}$ erg cm$^{-2}$ s$^{-1}$ (L${>}$4$\times$10$^{34}$ erg s$^{-1}$ at the distance of M33). These data allow complete coverage of the galaxy with high sensitivity to soft sources such as diffuse hot gas and supernova remnants. Here we describe the methods we used to identify and characterize 1296 point sources in the 8 fields. We compare our resulting source catalog to the literature, note variable sources, construct hardness ratios, classify soft sources, analyze the source density profile, and measure the X-ray luminosity function. As a result of the large effective area of XMM-Newton below 1 keV, the survey contains many new soft X-ray sources. The radial source density profile and X-ray luminosity function for the sources suggests that only $\sim$15% of the 391 bright sources with L${>}$3.6$\times$10$^{35}$ erg s$^{-1}$ are likely to be associated with M33, and more than a third of these are known supernova remnants. The log(N)--log(S) distribution, when corrected for background contamination, is a relatively flat power-law with a differential index of 1.5, which suggests many of the other M33 sources may be high-mass X-ray binaries. Finally, we note the discovery of an interesting new transient X-ray source, which we are unable to classify. △ Less

Submitted 17 March, 2015; originally announced March 2015.

Comments: 26 pages, 6 tables, 13 figures, accepted for publication in ApJS

Discovery of Chromium, Manganese, Nickel, and Copper Isotopes

Authors: K. Garofali, R. Robinson, M. Thoennessen

Abstract: Twenty-seven chromium, twenty-five manganese, thirty-one nickel and twenty-six copper isotopes have so far been observed and the discovery of these isotopes is discussed. For each isotope a brief synopsis of the first refereed publication, including the production and identification method, is presented. Twenty-seven chromium, twenty-five manganese, thirty-one nickel and twenty-six copper isotopes have so far been observed and the discovery of these isotopes is discussed. For each isotope a brief synopsis of the first refereed publication, including the production and identification method, is presented. △ Less

Submitted 7 December, 2010; originally announced December 2010.

Comments: to be published in Atomic Data Nuclear Data Tables