Search

Using hard (E>10 keV) X-ray observations with NuSTAR, we are able to differentiate between accretion states, and thus compact object types, of neutron stars and black holes in X-ray binaries (XRBs) in M31, our nearest Milky Way-type neighbor. Using ten moderate-depth (20-50 ks) observations of the disk of M31 covering a total of ~0.45 deg$^{2}$, we detect 20 sources at 2$\sigma$ in the 4-25 keV band pass, 14 of which we consider to be XRB candidates. This complements an existing deeper (100-400 ks) survey covering ~0.2 deg$^{2}$ of the bulge and the northeastern disk. We make tentative classifications of 9 of these sources with the use of diagnostic color-intensity and color-color diagrams, which separate sources into various neutron star and black hole regimes, identifying 3 black holes and 6 neutron stars. In addition, we create X-ray luminosity functions for both the full (4-25 keV) and hard (12-25 keV) band, as well as sub-populations of the full band based on compact object type and association with globular clusters. Our best fit globular cluster XLF is shallower than the field XLF, and preliminary BH and NS XLFs suggest a difference in shape based on compact object type. We find that the cumulative disk XLFs in the full and hard band are best fit by power laws with indices of 1.32 and 1.28 respectively. This is consistent with models of the Milky Way XLF from Grimm et al. (2002), Voss & Ajello (2010), and Doroshenko et al. (2014)., Comment: 19 pages, 6 Figures, Published in ApJ

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., Comment: Accepted for publication in Frontiers in Astronomy and Space Sciences

The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives., Comment: Published in Proceedings of SPIE Optics & Photonics 2023, San Diego

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 $\lambda$1640,4686, [Ne V] $\lambda$3426,14.3$\mu$m, and [O IV] 25.9$\mu$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., Comment: accepted for publication in ApJ; simulation results available through FSPS (see Section 5.4)

In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneous soft X-ray coverage. FORCE aims to be launched in the early 2030s, providing a perfect hard X-ray complement to the ESA flagship mission Athena. FORCE will be the most powerful X-ray probe for discovering obscured/hidden black holes and studying high energy particle acceleration in our Universe and will address how relativistic processes in the universe are realized and how these affect cosmic evolution. FORCE, which will operate over 1--79 keV, is equipped with two identical pairs of supermirrors and wideband X-ray imagers. The mirror and imager are connected by a high mechanical stiffness extensible optical bench with alignment monitor systems with a focal length of 12~m. A light-weight silicon mirror with multi-layer coating realizes a high angular resolution of $<15''$ in half-power diameter in the broad bandpass. The imager is a hybrid of a brand-new SOI-CMOS silicon-pixel detector and a CdTe detector responsible for the softer and harder energy bands, respectively. FORCE will play an essential role in the multi-messenger astronomy in the 2030s with its broadband X-ray sensitivity., Comment: 12 pages, 8 figures. Proceedings of SPIE Astronomical Telescopes and Instrumentation 2022

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., Comment: Accepted for publication in ApJ (25 pages, 9 figures, and 4 tables)

We present an analysis of the ultraluminous X-ray source (ULX) population in 75 Virgo cluster late-type galaxies, including all those with a star formation rate >~ 1 M_{sun}/yr and a representative sample of the less star-forming ones. This study is based on 110 observations obtained over 20 years with the Chandra X-ray Observatory Advanced Camera for Imaging Spectroscopy. As part of a Large Chandra Program, new observations were obtained for 52 of these 75 galaxies. The data are complete to a sensitivity of about 10^{39} erg/s, with a typical detection limit of about 3 x 10^{38} erg/s for the majority of the sources. The catalogue contains about 80 ULXs (0.3-10 keV luminosity >10^{39} erg/s), and provides their location, observed flux, de-absorbed luminosity, and (for the 25 most luminous ones) simple X-ray spectral properties. We discuss the ULX luminosity function in relation to the mass and star formation rate of the sample galaxies. We show that recent models of low-mass plus high-mass X-ray binary populations (scaling with stellar mass and star formation rate, respectively) are mostly consistent with our observational results. We tentatively identify the most luminous X-ray source in the sample (a source in IC 3322A with L_{X} ~ 6 x 10^{40} erg/s) as a recent supernova or its young remnant. The properties of the sample galaxies (morphologies, stellar masses, star formation rates, total X-ray luminosities from their point-source population) are also summarised., Comment: 24 pages, 200 KB, accepted by MNRAS (2022 January 17)

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., Comment: ApJ accepted. 25 pages, 8 figures, 4 tables

Galactic winds are associated with intense star formation and AGNs. Depending on their formation mechanism and velocity they may remove a significant fraction of gas from their host galaxies, thus suppressing star formation, enriching the intergalactic medium, and shaping the circumgalactic gas. However, the long-term evolution of these winds remains mostly unknown. We report the detection of a wind from NGC 3079 to at least 60 kpc from the galaxy. We detect the wind in FUV line emission to 60 kpc (as inferred from the broad FUV filter in GALEX) and in X-rays to at least 30~kpc. The morphology, luminosities, temperatures, and densities indicate that the emission comes from shocked material, and the O/Fe ratio implies that the X-ray emitting gas is enriched by Type II supernovae. If so, the speed inferred from simple shock models is about 500 km/s, which is sufficient to escape the galaxy. However, the inferred kinetic energy in the wind from visible components is substantially smaller than canonical hot superwind models., Comment: 15 pages, 8 figures, 3 tables. Accepted for publication in ApJ

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., Comment: Accepted for publication in ApJS. Data products and catalogs available at https://lehmer.uark.edu/downloads/

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., Comment: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

Stellar feedback -- stars regulating further star formation through the injection of energy and momentum into the interstellar medium -- operates through a complex set of processes that originate in star clusters but shape entire galaxies. A mature theory of stellar feedback is essential to a complete theory of star and galaxy formation, but the energy and momentum injected by hot gas into its surroundings remains unclear. With a next-generation X-ray observatory, breakthrough progress can be made through precision measurements of the temperature, density, velocity, and abundances of hot gas on scales of star clusters to galactic superwinds., Comment: white paper submitted to the Astro2020 Decadal Survey

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain., Comment: Accepted for publication in the Astrophysical Journal

We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combination of a hot thermal plasma model, a power-law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with a $3\sigma$ reflection fraction upper-limits $R<0.3$ and $R<0.18$ for NGC 3998 and NGC 4579, respectively; (2) the 6-7 keV line complex in NGC 4579 could either be fit with a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line, or 3 relatively narrow lines, which includes contribution from Fe XXVI; (3) NGC 4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy $E_{\rm cut}=107_{-18}^{+27}$ keV in NGC 3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) NGC 3998 spectrum is consistent with a Comptonization model with either a sphere ($\tau\approx3\pm1$) or slab ($\tau\approx1.2\pm0.6$) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGN, other LLAGN, and hot accretion flow models., Comment: 14 pages, 11 figures, 4 tables, accepted for publication in ApJ

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

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

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient., Comment: 12 pages, 3 figures, accepted to MNRAS

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], Comment: Accepted for publication in ApJ (13 pages, 10 figures)

The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit., Comment: Accepted for publication by ApJ

Deep Chandra ACIS observations of the region around the putative pulsar, CXOU J061705.3+222127, in the supernova remnant IC443 reveal an ~5$^{\prime\prime}$-radius ring-like structure surrounding the pulsar and a jet-like feature oriented roughly north-south across the ring and through the pulsar's location at 06$^{\rm h}$17$^{\rm m}$5.200$^{\rm s}$ +22$^{\circ}$21$^{\prime}$27.52$^{\prime\prime}$ (J2000.0 coordinates). The observations further confirm that (1) the spectrum and flux of the central object are consistent with a rotation-powered pulsar, (2) the non-thermal spectrum and morphology of the surrounding nebula are consistent with a pulsar wind and, (3) the spectrum at greater distances is consistent with thermal emission from the supernova remnant. The cometary shape of the nebula, suggesting motion towards the southwest, appears to be subsonic: There is no evidence either spectrally or morphologically for a bow shock or contact discontinuity; the nearly circular ring is not distorted by motion through the ambient medium; and the shape near the apex of the nebula is narrow. Comparing this observation with previous observations of the same target, we set a 99% confidence upper limit to the proper motion of CXOU J061705.3+222127 to be less than 44 mas/yr (310 km/s for a distance of 1.5 kpc), with the best-fit (but not statistically significant) projected direction toward the west., Comment: Accepted for publication in the Astrophysical Journal

Prior to the launch of NuSTAR, it was not feasible to spatially resolve the hard (E > 10 keV) emission from galaxies beyond the Local Group. The combined NuSTAR dataset, comprised of three ~165 ks observations, allows spatial characterization of the hard X-ray emission in the galaxy NGC 253 for the first time. As a follow up to our initial study of its nuclear region, we present the first results concerning the full galaxy from simultaneous NuSTAR, Chandra, and VLBA monitoring of the local starburst galaxy NGC 253. Above ~10 keV, nearly all the emission is concentrated within 100" of the galactic center, produced almost exclusively by three nuclear sources, an off-nuclear ultraluminous X-ray source (ULX), and a pulsar candidate that we identify for the first time in these observations. We detect 21 distinct sources in energy bands up to 25 keV, mostly consisting of intermediate state black hole X-ray binaries. The global X-ray emission of the galaxy - dominated by the off-nuclear ULX and nuclear sources, which are also likely ULXs - falls steeply (photon index >~ 3) above 10 keV, consistent with other NuSTAR-observed ULXs, and no significant excess above the background is detected at E > 40 keV. We report upper limits on diffuse inverse Compton emission for a range of spatial models. For the most extended morphologies considered, these hard X-ray constraints disfavor a dominant inverse Compton component to explain the {\gamma}-ray emission detected with Fermi and H.E.S.S. If NGC 253 is typical of starburst galaxies at higher redshift, their contribution to the E > 10 keV cosmic X-ray background is < 1%., Comment: 20 pages, 14 figures, accepted for publication in ApJ

Hinge clumps are luminous knots of star formation near the base of tidal features in some interacting galaxies. We use archival Hubble Space Telescope UV/optical/IR images and Chandra X-ray maps along with GALEX UV, Spitzer IR, and ground-based optical/near-IR images to investigate the star forming properties in a sample of 12 hinge clumps in five interacting galaxies. The most extreme of these hinge clumps have star formation rates of 1 - 9 M(sun)/yr, comparable to or larger than the `overlap' region of intense star formation between the two disks of the colliding galaxy system the Antennae. In the HST images, we have found remarkably large (~70 pc) and luminous (M(I) ~ 12.2 to -16.5) sources at the centers of these hinge clumps, sometimes embedded in a linear ridge of fainter star clusters. We have found strong X-ray emission from several of these hinge clumps. In most cases, this emission is well-resolved with Chandra and has a thermal X-ray spectrum, thus it is likely due to hot gas associated with the star formation. The ratio of the extinction-corrected diffuse X-ray luminosity to the mechanical energy rate (the X-ray production efficiency) for the hinge clumps is similar to that in the Antennae galaxies, but higher than those for regions in the normal spiral galaxy NGC 2403. Two of the hinge clumps have point-like X-ray emission much brighter than expected for hot gas; these sources are likely `ultra-luminous X-ray sources' (ULXs) due to accretion disks around black holes. The most extreme of these sources, in Arp 240, has a hard X-ray spectrum and an absorbed X-ray luminosity of ~2 X 10^41 erg/s; more than expected by single high mass X-ray binaries (HMXBs), thus it may be either a collection of HMXBs or an intermediate-mass black hole (>=80 M(sun)). [ abridged ], Comment: Astronomical Journal, in press

Observational confirmation of hot accretion model predictions has been hindered by the challenge to resolve spatially the Bondi radii of black holes with X-ray telescopes. Here, we use the Megasecond Chandra X-ray Visionary Project (XVP) observation of the NGC~3115 supermassive black hole to place the first direct observational constraints on the spatially and spectroscopically resolved structures of the X-ray emitting gas inside the Bondi radius of a black hole. We measured temperature and density profiles of the hot gas from a fraction out to tens of the Bondi radius (R_B = 2.4-4.8 arcsec = 112-224 pc). The projected temperature jumps significantly from ~0.3 keV beyond 5 arcsec to ~0.7 keV within ~4-5 arcsec, but then abruptly drops back to ~0.3 keV within ~3 arcsec. This is contrary to the expectation that the temperature should rise toward the center for a radiatively inefficient accretion flow. A hotter thermal component of ~1 keV inside 3 arcsec (~150 pc) is revealed using a two component thermal model, with the cooler ~0.3 keV thermal component dominating the spectra. We argue that the softer emission comes from diffuse gas physically located within $\sim 150$~pc from the black hole. The density profile is broadly consistent with rho ~ r^{-1} within the Bondi radius for either the single temperature or the two-temperature model. The X-ray data alone with physical reasoning argue against the absence of a black hole, supporting that we are witnessing the onset of the gravitational influence of the supermassive black hole., Comment: 19 pages including 4 appendix pages, 10 figures, accepted for publication in the Astrophysical Journal

(Abridged) We present a deep Chandra observation of the late-type barred spiral galaxy NGC 2903. The Chandra data reveal soft (kT_e ~ 0.2-0.5keV) diffuse emission in the nuclear starburst region and extending ~5kpc to the north and west of the nucleus. Much of this soft hot gas is likely to be from local active star-forming regions; however, besides the nuclear region, the morphology of hot gas does not strongly correlate with sites of active star formation. The central ~650 pc radius starburst zone exhibits much higher surface brightness diffuse emission than the surrounding regions and a harder spectral component in addition to its soft component. We interpret the hard component as being of thermal origin with kT_e~3.6keV and to be directly associated with a wind fluid produced by supernovae and massive star winds. The inferred terminal velocity for this hard component, ~1100 km/s, exceeds the local galaxy escape velocity suggesting a potential outflow. The softer extended emission does not display an obvious outflow geometry. However, the column density through which the X-rays are transmitted is lower to the west of the nucleus compared to the east and the surface brightness is higher there suggesting some soft hot gas originates from above the disk; viewed directly from the western zone but through the intervening galaxy disk from the eastern zone. There are several point-like sources in the nuclear region with X-ray spectra typical of compact binaries. None of these are coincident with the mass center of the galaxy and we place an upper limit luminosity from any point-like nuclear source to be < 2x10^38 ergs/s in the 0.5-8.0keV band which indicates that NGC 2903 lacks an active galactic nucleus. Heating from the nuclear starburst and a galactic wind may be responsible for preventing cold gas from accreting onto the galactic center., Comment: 17 pages accepted to ApJ

The M81 group member dwarf galaxy IC 2574 hosts a supergiant shell of current and recent star-formation activity surrounding a 1000 x 500 pc hole in the ambient Hi gas distribution. Chandra X-ray Observatory imaging observations reveal a luminous, L_x ~ 6.5 x 10^{38} erg/s in the 0.3 - 8.0 keV band, point-like source within the hole but offset from its center and fainter diffuse emission extending throughout and beyond the hole. The star formation history at the location of the point source indicates a burst of star formation beginning ~25 Myr ago and currently weakening and there is a young nearby star cluster, at least 5 Myr old, bracketing the likely age of the X-ray source at between 5 and ~25 Myr. The source is thus likely a bright high-mass X-ray binary --- either a neutron star or black hole accreting from an early B star undergoing thermal-timescale mass transfer through Roche lobe overflow. The properties of the residual diffuse X-ray emission are consistent with those expected from hot gas associated with the recent star-formation activity in the region., Comment: 5 pages, accepted for publication in the Astrophysical Journal Letters

One hundred seven ultraluminous X-ray (ULX) sources with 0.3-10.0 keV luminosities in excess of 1e39 erg/s are identified in a complete sample of 127 nearby galaxies. The sample includes all galaxies within 14.5 Mpc above the completeness limits of both the Uppsala Galaxy Catalog and the Infrared Astronomical Satellite survey. The galaxy sample spans all Hubble types, a four decade range in mass and in star-formation rate. ULXs are detected in this sample at rates of one per 3.2e10 solar mass, one per 0.5 solar mass/year star-formation rate, and one per 57 cubic Mpc corresponding to a luminosity density of ~2e37 erg/s/Mpc3. At these rates we estimate as many as 19 additional ULXs remain undetected in fainter dwarf galaxies within the survey volume. An estimated 14 or 13%, of the 107 ULX candidates are expected to be background sources. The differential ULX luminosity function shows a power law slope of -1.2 to -2.0 with an exponential cutoff at 2e40 erg/s with precise values depending on the model and on whether the ULX luminosities are estimated from their observed numbers of counts or, for a subset of candidates, from their spectral shapes. Extrapolating the observed luminosity function predicts at most one very luminous ULX, L~1e41 erg/s, within a distance as small as 100 Mpc. The luminosity distribution of ULXs within the local universe cannot account for the recent claims of luminosities in excess of 2e41 erg/s requiring a new population class to explain these extreme objects., Comment: 13 pages accepted to ApJ, revised Figure 6 and Table 2, thanks to Hua Feng

Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are resolving, for the first time, the accretion flow within the Bondi radius of an SMBH. We show that the temperature is rising toward the galaxy center as expected in all accretion models in which the black hole is gravitationally capturing the ambient gas. There is no hard central point source that could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 4 arcsec-5 arcsec (188-235 pc), suggesting an SMBH of 2 x 10^9 M_sun that is consistent with the upper end of the optical results. The density profile within the Bondi radius has a power-law index of 1.03^{+0.23}_{-0.21} which is consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate at the Bondi radius is determined to be {\dot M}_B = 2.2 x 10^{-2} M_sun yr^{-1}. Thus, the accretion luminosity with 10% radiative efficiency at the Bondi radius (10^{44} erg s^{-1}) is about six orders of magnitude higher than the upper limit of the X-ray luminosity of the nucleus., Comment: 5 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters. Version 2 is slightly modified to match the publication

Archival Chandra observations are used to study the X-ray emission associated with star formation in the central region of the nearby SAB(s)cd galaxy NGC 2403. The distribution of X-ray emission is compared to the morphology visible at other wavelengths using complementary Spitzer, GALEX, and ground-based Halpha imagery. In general, the brightest extended X-ray emission is associated with HII regions and to other star-forming structures but is more pervasive; existing also in regions devoid of strong Halpha and UV emission. This X-ray emission has the spectral properties of diffuse hot gas (kT ~ 0.2keV) whose likely origin is in gas shock-heated by stellar winds and supernovae with < 20% coming from faint unresolved X-ray point sources. This hot gas may be slowly-cooling extra-planar remnants of past outflow events, or a disk component that either lingers after local star formation activity has ended or that has vented from active star-forming regions into a porous interstellar medium., Comment: 25 pages, accepted to AJ

We report the discovery of an X-ray source coincident with the nuclear star cluster at the dynamical center of the nearby late-type spiral galaxy NGC 2403. The X-ray luminosity of this source varies from below detection levels, ~1e35 erg/s in the 0.5-8.0 keV band, to 7e38 erg/s on timescales between observations of less than 2 months. The X-ray spectrum is well-fit by an accretion disk model consisting of multiple blackbody components and corresponding physically to a compact object mass of greater than approximately 5 solar masses. No pulsations nor aperiodic behavior is evident in its X-ray light curve on the short timescales of the individual observations. The X-ray properties of the source are more similar to those of the nuclear source X-8 in M33, believed to be a low-mass X-ray binary, then to those of the low-luminosity active galactic nucleus in NGC 4395. The brightness of the nuclear star cluster, M_I ~ -11.8 mag, is typical of clusters in late-type spirals but its effective radius, r_e \~ 12 pc, is several times larger than average indicating a relatively relaxed cluster and a low probability of a central massive object. The cluster has a mass 3e6 solar masses and an age of 1.4 Gyr estimating from its observed colors and brightness., Comment: 7 pages, submitted to ApJ

We present submillimeter continuum observations at 870 microns of the cores of low redshift 3CRR radio galaxies, observed at the Heinrich Hertz Submillimeter Telescope. The cores are nearly flat spectrum between the radio and submillimeter which implies that the submillimeter continuum is likely to be synchrotron emission and not thermal emission from dust. The emitted power from nuclei detected at optical wavelengths and in the X-rays is similar in the submillimeter, optical and X-rays. The submillimeter to optical and X-ray power ratios suggest that most of these sources resemble misdirected BL Lac type objects with synchrotron emission peaking at low energies. However we find three exceptions, the FR I galaxy 3C264 and the FR II galaxies 3C390.3 and 3C338 with high X-ray to submillimeter luminosity ratios. These three objects are candidate high or intermediate energy peaked BL Lac type objects. With additional infrared observations and from archival data, we compile spectral energy distributions (SEDs) for a subset of these objects. The steep dips observed near the optical wavelengths in many of these objects suggest that extinction inhibits the detection and reduces the flux of optical continuum core counterparts. High resolution near or mid-infrared imaging may provide better measurements of the underlying synchrotron emission peak., Comment: accepted for publication in AJ

We measure reddenings to HII regions in NGC 2903, NGC 1512, M51, NGC 4449 and NGC 6946 from Hubble Space Telescope Pa alpha and H alpha images. Extinctions range from A_V ~ 5 - 0 depending upon the galaxy. For the galaxies with HST images in both lines, NGC 2903, NGC 1512 and M51, the Pa alpha and H alpha emission are almost identical in morphology which implies that little emission from bright HII regions is hidden from view by regions of comparatively high extinction. The scatter in the measured extinctions is only +- 0.5 mag. We compare the reddenings we measure in five galaxies using the Pa alpha to H alpha ratios to those measured previously from the Balmer decrement in the LMC and as a function of radius in M101 and M51. We find that luminosity weighted mean extinctions of these ensembles of HI regions are correlated with gas surface density and metallicity. The correlation is consistent with the mean extinction depending on dust density where the dust to gas mass ratio scales with the metallicity. This trend is expected if HII regions tend to be located near the mid-plane of a gas disk and emerge from their parent molecular clouds soon after birth. In environments with gas densities below a few hundred Msol/pc^2 star formation rates estimated from integrated line fluxes and mean extinctions are likely to be fairly accurate., Comment: accepted for publication in AJ

We report the analysis of simultaneous XMM-Newton+Nuclear Spectroscopic Telescope Array (NuSTAR) observations of two low-luminosity active galactic nuclei (LLAGNs), NGC3998 and NGC4579. We do not detect any significant variability in either source over the ∼3 day length of the NuSTAR observations. The broadband 0.5–60 keV spectrum of NGC3998 is best fit with a cutoff power law, while the one for NGC4579 is best fit with a combination of a hot thermal plasma model, a power law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with 3σ reflection fraction upper limits of R < 0.3 and R < 0.18 for NGC3998 and NGC4579, respectively; (2) the 6–7 keV line complex in NGC4579 could be fit with either a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line or with three relatively narrow lines, which include contribution from Fe XXVI; (3) the NGC4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy = - E 107+ cut 18 27 keV in NGC3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) the NGC3998 spectrum is consistent with a Comptonization model with either a sphere (τ ≈ 3 ± 1) or slab (τ ≈ 1.2 ± 0.6) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGNs, other LLAGNs, and hot accretion flow models.

We present the development of the reflective coating by magnetron sputtering deposition onto precisely-fabricated thin X-ray mirrors. Our goal is to remove distortion induced by the coating and then keep their surface pro les. We first addressed the uniform coating to minimize the distortion by introducing a mask to control the spatial distribution of the coating thickness. The uniformity was finally achieved within 1%. We next tried a platinum single-layer coating on a glass substrate with a dimension of 200 mm 125 mm. The distortion caused by the frontside coating with a thickness of 320 A was found to be at most 1 m, smaller than the previous results obtained from the non-uniform coating. We then carried out the platinum coating with the same amount of the thickness on the backside surface of the glass substrate. The surface pro le of the glass substrate was fully recovered, indicating that the residual stress was successfully balanced by the backside coating. Furthermore, we tried to an iridium single-layer coating with a thickness of 150 Aon the silicon mirrors. The frontside coating caused the degradation of the imaging quality by 7:5 arcsec in half-power width. However, the backside coating with the same amount of the thickness reduced this degradation to be 3:4 arcsec. Finally, an additional backside coating with a thickness of 100 A and the annealing to relax the residual stress were found to eliminate the distortion completely; the final degradation of the imaging quality was only 0:4 arcsec.

The recent discovery by Bachetti et al. of a pulsar in M82 that can reach luminosities of up to 10(exp 40) erg s(exp -1), a factor of approximately 100 times the Eddington luminosity for a 1.4 solar mass compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.58 keV X-ray emission from this source from 15 years of Chandra observations. We analyze 19 ACIS observations where the point-spread function (PSF) of the pulsar is not contaminated by nearby sources. We fit the Chandra spectra of the pulsar with a power-law model and a disk blackbody model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where four observations have a pile-up fraction of 10, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (LX greater than 10(exp 39) erg s(exp -1) is equal to gamma 1.33 +/-.0.15. For the disk blackbody model, the average temperature is T(sub in) 3.24 +/- 0.65 keV, the spectral shape being consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cutoff, where gamma is equal to 0.6 +/- 0.3 and E(sub C) is equal to 14(exp +5) (sub -3)) keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10( exp (39) erg s (exp −1), greater than 10 times its Eddington limit.

We present the X-ray luminosity function (XLF) of ultraluminous X-ray (ULX) sources with 0.3-10.0 keV luminosities in excess of 10(sup 39) erg/s in a complete sample of nearby galaxies. The XLF shows a break or cut-off at high luminosities that deviates from its pure power law distribution at lower luminosities. The cut-off is at roughly the Eddington luminosity for a 90-140 solar mass accretor. We examine the effects on the observed XLF of sample biases, of small-number statistics (at the high luminosity end) and of measurement uncertainties. We consider the physical implications of the shape and normalization of the XLF. The XLF is also compared and contrasted to results of other recent surveys.