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 (AGNs). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGNs. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include adaptive Markov Chain Monte Carlo (MCMC), affine-invariant MCMC, and Levenberg–Marquardt gradient-descent (
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 constraints 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 IDL and is available at
https://github.com/rafaeleufrasio/lightning
.
Binary population synthesis combined with cosmological models suggests 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 ultraluminous X-ray sources (ULXs) with luminosities above 1040 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, as well as the importance of X-ray emission from star-forming galaxies in the high-redshift universe.
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 ∼10) compared ...with predictions from scaling relations of galaxy-wide high-mass X-ray binary (HMXB) luminosity with star formation rate. In the nearby (≈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 (≲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
IR
within these regions on timescales ≲5 Myr, shorter than the timescale required for the formation of HMXBs.
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 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 yr by including archival X-ray observations of these sources. We compare the detection and nondetection 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.
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 ultraluminous X-ray sources 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 (∼0.1
Z
⊙
), 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 we discuss how the ionization signatures compare to models for sources such as intermediate-mass black holes.
Abstract
We present multiwavelength 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 catalog, which covers the inner disk of M33 at near-IR, optical, and near-UV wavelengths. We perform spectral energy distribution fitting on multiband 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 X-ray binaries (Be-XRBs), as is seen in other Local Group galaxies. We use spatially resolved recent star formation history 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 ∼10 and ∼40 Myr, which match theoretical formation timescales for the most massive HMXBs and Be-XRBs, respectively. We measure an HMXB production rate of 107–136 HMXBs/(
M
⊙
yr
−1
) over the last 50 Myr and 150–199 HMXBs/(
M
⊙
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 black hole HMXBs and main-sequence companion stars in neutron star HMXBs.
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 ∼800 galaxies and integrated emission from ∼6,000 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.
Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks, 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 interstellar medium 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. This white paper is part of a series commissioned for the AXIS Probe mission concept; additional AXIS white papers can be found at the AXIS website.
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 ≈ 200–450 Mpc that have high specific star formation rates of 3–9 Gyr^(−1) and metallicities near Z ≈ 0.3Z_⊙, 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^(HMXB)_(0.5-8keV) / SFR = 40.19 ± 0.06 and log L^(gas)_(0.5-2keV) / SFR = 39.58 ^(+0.17)_(-0.28), significantly elevated compared to local relations. The HMXB scaling is also somewhat higher than L^(HMXB)_(0.5-8keV) -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^(gas)_(0.5-2keV) -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.
In this work, 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 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, showing steadily increasing numbers of luminous and ultraluminous X-ray sources ( (erg s−1) = 38-40.5) with declining metallicity. However, we find that the low-luminosity ( (erg s−1) = 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 , and a new characterization of its 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.