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.

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams.

X-ray data for quasar 4C 74.26 have previously been modeled with a broad Fe K$\alpha$ emission line and reflection continuum originating in the inner part of the accretion disk around the central supermassive black hole (SMBH), i.e. the strong gravity regime. We modeled broadband X-ray spectra from $Suzaku$ and $NuSTAR$ with MYTORUS, self-consistently accounting for Fe K$\alpha$ line emission, as well as direct and reflected continuum emission, from finite column density matter. A narrow Fe K$\alpha$ emission line originating in an X-ray reprocessor with solar Fe abundance far from the central SMBH is sufficient to produce excellent fits for all spectra. For the first time, we are able to measure the global, out of the line-of-sight column density to be in the range $\sim$$1.5$ to $\sim$$2.9\times10^{24}$ cm$^{-2}$, i.e. in the Compton thick regime, while the line-of-sight column density is Compton thin in all observations. The Fe K$\alpha$ emission line is unresolved in all but one observations. The Compton scattered continuum from distant matter removes the need for relativistic broadening of the Fe K$\alpha$ emission line, which is required for SMBH spin measurements. The resolved line observation can alternatively be modeled with a relativistic model but we do not find evidence for a truncated accretion disk model. We conclude that the X-ray emission in these 4C 74.26 data is unlikely to originate in the inner accretion disk region and thus cannot be used to measure SMBH spin.

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.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

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.

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.

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.

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.

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.

We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard state black holes, pulsars, and non-magnetized neutron stars associated with optical point source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates --- X-ray sources with UV-bright point source optical counterpart candidates --- are found in regions with star formation bursts less than 50 Myr ago, with 3 associated with young stellar ages (<10 Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403.

We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr.

We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age.

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary.

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

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.

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]

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.

Lyman Break Analogs (LBA) are local proxies of high-redshift Lyman Break Galaxies (LBG). Studies of nearby starbursts have shown that Lyman continuum and line emission are absorbed by dust and that the Lyman-alpha is resonantly scattered by neutral hydrogen. A source of feedback is required to prevent scattering and allow the Lyman-alpha emission to escape. There are two X-ray point sources embedded in the Lyman Break Analog (LBA) galaxy Haro 11. Haro 11 X-1 is an extremely luminous (L$_{X} \sim 10^{41}$ ergs s$^{-1}$), spatially compact source with a hard X-ray spectrum. Haro 11 X-1 is similar to the extreme Black Hole Binary (BHB) M82 X-1. The hard X-ray spectrum indicates Haro 11 X-1 may be a Black Hole Binary (BHB) in a low accretion state. The very high X-ray luminosity suggests an intermediate mass black hole that could be the seed for formation of a supermassive black hole. Source Haro 11 X-2 has an X-ray luminosity L$_{X} \sim 5\times10^{40}$ ergs s$^{-1}$ and a soft X-ray spectrum. This strongly suggests that Haro 11 X-2 is an X-ray binary in the ultra luminous state. Haro 11 X-2 is coincident with the star forming knot that is the source of the Lyman-alpha emission, raising the possibility that strong winds from X-ray binaries play an important part in injecting mechanical power into the Interstellar Medium (ISM), thus blowing away neutral material from the starburst region and allowing the Lyman-alpha to escape. We suggest that feedback from X-ray binaries may play a significant role in allowing Lyman-alpha emission to escape from galaxies in the early universe.

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.

We present nearly simultaneous Chandra and NuSTAR observations of two actively star-forming galaxies within 50 Mpc: NGC 3256 and NGC 3310. Both galaxies are detected by both Chandra and NuSTAR, which together provide the first-ever spectra of these two galaxies spanning 0.3-30 keV. The X-ray emission from both galaxies is spatially resolved by Chandra; we find that hot gas dominates the E < 1-3 keV emission while ultraluminous X-ray sources (ULXs) dominate at E > 1-3 keV. The NuSTAR galaxy-wide spectra of both galaxies follow steep power-law distributions with Gamma ~ 2.6 at E > 5-7 keV, similar to the spectra of bright individual ULXs and other galaxies that have been studied by NuSTAR. We find that both NGC 3256 and NGC 3310 have X-ray detected sources coincident with nuclear regions; however, the steep NuSTAR spectra of both galaxies restricts these sources to be either low luminosity AGN or non-AGN in nature (e.g., ULXs or crowded X-ray sources that reach L2-10 keV ~ 10^40 erg/s cannot be ruled out). Combining our constraints on the 0.3-30 keV spectra of NGC 3256 and NGC 3310 with equivalent measurements for nearby star-forming galaxies M83 and NGC 253, we analyze the SFR-normalized spectra of these starburst galaxies. The spectra of all four galaxies show sharply declining power-law slopes above 3-6 keV due to ULX populations. Our observations therefore constrain the average spectra of luminous accreting binaries (i.e., ULXs). This result is similar to the super-Eddington accreting ULXs that have been studied individually in a targeted NuSTAR ULX program. We also find that NGC 3310 exhibits a factor of ~3-10 elevation of X-ray emission over the other star-forming galaxies. We argue that the excess is most likely explained by the relatively low metallicity of the young stellar population in NGC 3310.

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR.

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

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

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 ]

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.

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

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.

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.

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.

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.

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.

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.

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.