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.

We can learn about the formation and evolution of compact objects, such as neutron stars and black holes (BHs), by studying the X-ray emission from accreting systems in nearby star-forming galaxies. The hard (E > 10 keV) X-ray emission in particular allows strong discrimination among the accretion states and compact object types. We conducted a NuSTAR survey (∼600 ks) of the Local Group spiral galaxy M33 to study the distribution of X-ray binary (XRB) accretors in an actively star-forming environment. We constructed color–intensity and color–color diagrams to infer XRB accretion states. Using these diagrams, we have classified 28 X-ray sources in M33 by comparing their hard X-ray colors to those of known systems. Four sources lie in the parameter space occupied by X-ray pulsars, while 8, 10, and 4 sources lie in the parameter space occupied by BHs in the hard, intermediate, and soft states, respectively. The known ultraluminous X-ray source M33 X-8 is also found to be consistent with that source type. Some sources overlap within the Z/Atoll sources due to the overlap of the two categories of BHs and Z/Atoll sources. In contrast to a similar NuSTAR survey of M31 (with a low-mass XRB-dominant population), the source population in M33 is dominated by high-mass XRBs (HMXBs), allowing the study of a very different population with similar sensitivity due to the galaxy's similar distance. This characterization of a population of HMXB accretion states will provide valuable constraints for theoretical XRB population synthesis studies to their formation and evolution.

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

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.