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.

We present the results of the photometric analysis of a large part of the main body of the Small Magellanic Cloud. Using the 6.5m Magellan Telescope at the Las Campanas Observatory in Chile, we have acquired deep B and I images in four fields (0.44 degree each in diameter), yielding accurate photometry for 1,068,893 stars down to 24th magnitude, with a spatial resolution of 0.20 arcsec per pixel. Colour-magnitude diagrams and (completeness corrected) luminosity functions have been constructed, yielding significant new results that indicate at least two discrete star formation events over a period from 2.7 to 4 Gyr ago. Also, we have derived star formation rates as a function of look back time and have found enhancements of SF between 4-6 Gyr and at younger ages.

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.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Flares are powerful events ignited by a sudden release of magnetic energy. With the aim of studying flares in the 125-Myr-old stars in the Pleiades observed simultaneously in optical and X-ray light, we obtained new XMM-Newton observations of this cluster during the observations of Kepler K2 Campaign 4. Our objective is to characterize the most powerful flares observed in both bands and to constrain the energy released in the optical and X-ray, the geometry of the loops, and their time evolution. We aim to compare our results to existing studies of flares occurring in the Sun and stars at different ages. We selected bright X-ray/optical flares occurred in 12 known members of the Pleiades from their K2 and XMM-Newton light curves. The sample includes ten K-M stars, one F9 star, and one G8 star. Flare average properties were obtained from integrated analysis of the light curves during the flares. The time evolution of the plasma in the magnetic loops is constrained with time-resolved X-ray spectral analysis. Most of the flares studied in this work emitted more energy in optical than in X-rays, as in most solar flares, even if the Pleiades flares output a larger fraction of their total energy in X-rays than typical solar flares do. Additionally, the energy budget in the two bands is weakly correlated. We also found comparable flare duration in optical and X-rays and observed that rapidly rotating stars (e.g., with rotation period shorter than 0.5 days) preferentially host short flares. We estimated the slope of the cooling path of the flares in the log(EM)-versus-log(T) plane. The values we obtained are affected by large uncertainties, but their nominal values suggest that the flares analyzed in this paper are mainly due to single loops with no sustained heating occurring during the cooling phase. We also observed and analyzed oscillations with a period of 500 s during one of the flares.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

Nearby star-forming galaxies offer a unique environment to study the populations of young ($<$100 Myr) accreting binaries. These systems are tracers of past populations of massive stars that heavily affect their immediate environment and parent galaxies. Using a Chandra X-ray Visionary program, we investigate the young neutron-star binary population in the low metallicity of the Small Magellanic Cloud (SMC) by reaching quiescent X-ray luminosity levels ($\sim$few times $10^{32}$ erg/s). We present the first measurement of the formation efficiency of high-mass X-ray binaries (HMXBs) as a function of the age of their parent stellar populations by using 3 indicators: the number ratio of HMXBs to OB stars, to the SFR, and to the stellar mass produced during the specific star-formation burst they are associated with. In all cases, we find that the HMXB formation efficiency increases as a function of time up to $\sim$40-60 Myr, and then gradually decreases.

Almost all confirmed optical counterparts of HMXBs in the SMC are OB stars with equatorial decretion disks (OBe). These sources emit strongly in Balmer lines and standout when imaged through narrow-band H{\alpha} imaging. The lack of secure counterparts for a significant fraction of the HMXBs motivated us to search for more. Using the catalogs for OB/OBe stars (Maravelias et al. 2017) and for HMXBs (Haberl & Sturm 2016) we detect 70 optical counterparts (out of 104 covered by our survey). We provide the first identification of the optical counterpart to the source XTEJ0050-731. We verify that 17 previously uncertain optical counterparts are indeed the proper matches. Regarding 52 confirmed HMXBs (known optical counterparts with H{\alpha} emission), we detect 39 as OBe and another 13 as OB stars. This allows a direct estimation of the fraction of active OBe stars in HMXBs that show H{\alpha} emission at a given epoch to be at least $\sim75\%$ of their total HMXB population.

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.

The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates.

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 report the timing analysis results of X-ray pulsars from a recent deep Chandra survey of the Small Magellanic Cloud (SMC). We have analyzed a total exposure of 1.4 Ms from 31 observations over a 1.2 deg$^2$ region in the SMC under a Chandra X-ray Visionary Program. Using the Lomb-Scargle and epoch folding techniques, we have detected periodic modulations from 20 pulsars and a new candidate pulsar. The survey also covers 11 other pulsars with no clear sign of periodic modulation. The 0.5-8 keV X-ray luminosity ($L_X$) of the pulsars ranges from $10^{34}$ to $10^{37}$ erg s$^{-1}$ at 60 kpc. All the Chandra sources with $L_X$ $\gtrsim 4 \times 10^{35}$ erg s$^{-1}$ exhibit X-ray pulsations. The X-ray spectra of the SMC pulsars (and high mass X-ray binaries) are in general harder than those of the SMC field population. All but SXP~8.02 can be fitted by an absorbed power-law model with a photon index of $\Gamma$ $\lesssim$ 1.5. The X-ray spectrum of the known magnetar SXP~8.02 is better fitted with a two-temperature blackbody model. Newly measured pulsation periods of SXP~51.0, SXP~214 and SXP~701 are significantly different from the previous XMM-Newton and RXTE measurements. This survey provides a rich data set for energy-dependent pulse profile modeling. Six pulsars show an almost eclipse-like dip in the pulse profile. Phase-resolved spectral analysis reveals diverse spectral variation during pulsation cycle: e.g., for an absorbed power-law model, some exhibit an (anti)-correlation between absorption and X-ray flux, while others show more intrinsic spectral variation.

The Small Magellanic Cloud (SMC) hosts a large number of high-mass X-ray binaries, and in particular of Be/X-ray Binaries (BeXRBs; neutron stars orbiting OBe-type stars), offering a unique laboratory to address the effect of metalicity. One key property of their optical companion is H{\alpha} in emission, which makes them bright sources when observed through a narrow-band H{\alpha} filter. We performed a survey of the SMC Bar and Wing regions using wide-field cameras (WFI@MPG/ESO and MOSAIC@CTIO/Blanco) in order to identify the counterparts of the sources detected in our XMM-Newton survey of the same area. We obtained broad-band R and narrow-band H{\alpha} photometry, and identified ~10000 H{\alpha} emission sources down to a sensitivity limit of 18.7 mag (equivalent to ~B8 type Main Sequence stars). We find the fraction of OBe/OB stars to be 13% down to this limit, and by investigating this fraction as a function of the brightness of the stars we deduce that H{\alpha} excess peaks at the O9-B2 spectral range. Using the most up-to-date numbers of SMC BeXRBs we find their fraction over their parent population to be ~0.002-0.025 BeXRBs/OBe, a direct measurement of their formation rate.

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.

Using hundreds of XMM-Newton and Chandra archival observations and nearly a thousand RXTE observations, we have generated a comprehensive library of the known pulsars in the Small and Large Magellanic Clouds (SMC, LMC). The pulsars are detected multiple times across the full parameter spaces of X-ray luminosity ($L_X= 10^{31-38}$~erg/s) and spin period ( P$<$1s -- P$>$1000s) and the library enables time-domain studies at a range of energy scales. The high time-resolution and sensitivity of the EPIC cameras are complemented by the angular resolution of Chandra and the regular monitoring of RXTE. Our processing %$\sim$15 year pipeline uses the latest calibration files and software to generate a suite of useful products for each pulsar detection: event lists, high time-resolution light curves, periodograms, spectra, and complete histories of $\dot{P}$, the pulsed fraction, etc., in the broad (0.2-12 keV), soft (0.2-2 keV), and hard (2-12 keV) energy bands. After combining the observations from these telescopes, we found that 28 pulsars show long-term spin up and 25 long-term spin down. We also used the faintest and brightest sources to map out the lower and upper boundaries of accretion-powered X-ray emission: the propeller line and the Eddington line, respectively. We are in the process of comparing the observed pulse profiles to geometric models of X-ray emission in order to constrain the physical parameters of the pulsars. Finally we are preparing a public release of the library so that it can be used by others in the astronomical community.

In the present work we investigate the link between high-mass X-ray binaries (HMXBs) and star formation in the Large Magellanic Cloud (LMC), our nearest star-forming galaxy. Using optical photometric data, we identify the most likely counterpart of 44 X-ray sources. Among the 40 HMXBs classified in this work, we find 33 Be/X-ray binaries, and 4 supergiant XRBs. Using this census and the published spatially resolved star-formation history map of the LMC, we find that the HMXBs (and as expected the X-ray pulsars) are present in regions with star-formation bursts $\sim$6-25 Myr ago, in contrast to the Small Magellanic Cloud (SMC), for which this population peaks at later ages ($\sim$25-60 Myr ago). We also estimate the HMXB production rate to be equal to 1 system per $\sim43.5\times10^{-3}$ Mo/yr, or 1 system per $\sim$143 Mo of stars formed during the associated star-formation episode. Therefore, the formation efficiency of HMXBs in the LMC is $\sim$17 times lower than that in the SMC. We attribute this difference primarily in the different ages and metallicity of the HMXB populations in the two galaxies. We also set limits on the kicks imparted on the neutron star during the supernova explosion. We find that the time elapsed since the supernova kick is $\sim$3 times shorter in the LMC than the SMC. This in combination with the average offsets of the HMXBs from their nearest star clusters results in $\sim$4 times faster transverse velocities for HMXBs in the LMC than in the SMC.

We have endeavoured to understand the formation and evolution of the black hole (BH) X-ray binary LMC X-3. We estimate the properties of the system at 4 evolutionary stages: 1) at the Zero Age Main Sequence (ZAMS), 2) just prior to the supernova (SN) explosion of the primary, 3) just after the SN, and 4) at the moment of RLO onset.We use a hybrid approach, combining detailed stellar structure and binary evolution calculations with approximate population synthesis models. This allows us to estimate potential natal kicks and the evolution of the BH spin. In the whole analysis we incorporate as model constraints the most up-to-date observational information, encompassing the binary's orbital properties, the companion star mass, effective temperature, surface gravity and radius, as well as the black hole's mass and spin. We find that LMC X-3 began as a ZAMS system with the mass of the primary star in the range $M_{\rm{1,ZAMS}}$ = 22-31 $\rm{M_{\odot}}$ and a secondary star of $M_{\rm{2,ZAMS}} = 5.0-8.3M_{\odot}$, in a wide ($P_{ZAMS} \gtrsim 2.000\, \rm days$) and eccentric ($e_{\rm{ZAMS}} \gtrsim 0.23$) orbit. Just prior to the SN, the primary has a mass of $M_{\mathrm{1,preSN}} = 11.1-18.0\,\rm M_{\odot}$, with the secondary star largely unaffected. The orbital period decreases to $0.6-1.7\, \rm days$, and is still eccentric $0 \leq e_{\rm{preSN}} \leq 0.44$. We find that a symmetric SN explosion with no or small natal kicks (a few tens of $\rm km\, s^{-1}$) imparted on the BH cannot be formally excluded, however, large natal kicks in excess of $\gtrsim 120 \,\rm km\, s^{-1}$ increase the estimated formation rate by an order of magnitude. Following the SN, the system has a BH $M_{\mathrm{BH,postSN}} = 6.4-8.2\,\rm M_{\odot}$ and is put into an eccentric orbit. At the RLO onset the orbit is circularised and it has an orbital period of $P_{\rm{RLO}} = 0.8-1.4\,\rm days$.

Located in the Small Magellanic Cloud (SMC), SXP214 is an X-ray pulsar in a high mass X-ray binary system with a Be-star companion. A recent survey of the SMC under a Chandra X-ray Visionary program found the source in a transition when the X-ray flux was on a steady rise. The Lomb-Scargle periodogram revealed a pulse period of 211.49 +/- 0.42 s, which is significantly (>5sigma) shorter than the previous measurements with XMM-Newton and RXTE. This implies that the system has gone through sudden spin-up episodes recently. The pulse profile shows a sharp eclipse-like feature with a modulation amplitude of >95%. The linear rise of the observed X-ray luminosity from <~2x to 7x10^35 erg s^-1 is correlated with steady softening of the X-ray spectrum, which can be described by the changes in the local absorption from N_H ~ 10^24 to <~10^20 cm^-2 for an absorbed power-law model. The soft X-ray emission below 2 keV was absent in the early part of the observation when only the pulsating hard X-ray component was observed, whereas at later times both soft and hard X-ray components were observed pulsating. A likely explanation is that the neutron star was initially hidden in the circumstellar disk of the companion, and later came out of the disk with the accreted material that continued fueling the observed pulsation.

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

The Magellanic clouds are uniquely placed to study the stellar contribution to dust emission. Individual stars can be resolved in these systems even in the mid-infrared, and they are close enough to allow detection of infrared excess caused by dust.We have searched the Spitzer Space Telescope data archive for all Infrared Spectrograph (IRS) staring-mode observations of the Small Magellanic Cloud (SMC) and found that 209 Infrared Array Camera (IRAC) point sources within the footprint of the Surveying the Agents of Galaxy Evolution in the Small Magellanic Cloud (SAGE-SMC) Spitzer Legacy programme were targeted, within a total of 311 staring mode observations. We classify these point sources using a decision tree method of object classification, based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information. We find 58 asymptotic giant branch (AGB) stars, 51 young stellar objects (YSOs), 4 post-AGB objects, 22 Red Supergiants (RSGs), 27 stars (of which 23 are dusty OB stars), 24 planetary nebulae (PNe), 10Wolf-Rayet (WR) stars, 3 Hii regions, 3 R Coronae Borealis (R CrB) stars, 1 Blue Supergiant and 6 other objects, including 2 foreground AGB stars. We use these classifications to evaluate the success of photometric classification methods reported in the literature.

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 the results of a pair of 100 ksec Chandra observations in the Small Magellanic Cloud to survey high mass X-ray binaries (HMXBs), stars and LMXBs/CVs down to Lx = 4.3 x 10^32 erg/s The two SMC Deep Fields are located in the most active star forming region of the bar, with Deep Field-1 positioned at the most pulsar-rich location identified from previous surveys. Two new pulsars were discovered in outburst: CXOU J004929.7-731058 (P=892s), CXOU J005252.2-721715 (P=326s), and 3 new HMXB candidates were identified. Of 15 Be-pulsars now known in the field, 13 were detected, with pulsations seen in 9 of them. Ephemerides demonstrate that 6 of the 10 pulsars known to exhibit regular outbursts were seen outside their periastron phase, and quiescent X-ray emission at Lx=10^33 - 10^34 is shown to be common. Comparison with ROSAT, ASCA, XMM-Newton catalogs resulted in positive identification of several previously ambiguous sources. Bright optical counterparts exist for 40 of the X-ray sources, of which 33 are consistent with early-type stars Mv<-2, B-V<0.2), and are the subject of a companion paper. The results point to an underlying HMXB population-density up to double that of active systems. The full catalog of 394 point-sources is presented along with detailed analyses of timing and spectral properties.

We have discovered a persistent, but highly variable X-ray source in the nearby starburst galaxy NGC 253. The source varies at the level of a factor of about 5 in count rate on timescales of a few hours. Two long observations of the source with Chandra and XMM-Newton show suggestive evidence for the source having a period of about 14-15 hours, but the time sampling in existing data is insufficient to allow a firm determination that the source is periodic. Given the amplitude of variation and the location in a nuclear starburst, the source is likely to be a Wolf-Rayet X-ray binary, with the tentative period being the orbital period of the system. In light of the fact that we have demonstrated that careful examination of the variability of moderately bright X-ray sources in nearby galaxies can turn up candidate Wolf-Rayet X-ray binaries, we discuss the implications of Wolf-Rayet X-ray binaries for predictions of the gravitational wave source event rate, and, potentially, interpretations of the events.

The Small Magellanic Cloud (SMC) is well known to harbor a large number of High-Mass X-ray Binaries (HMXBs). The identification of their optical counterparts provides information on the nature of the donor stars and can help to constrain the parameters of these systems and their evolution. We obtained optical spectra for a number of HMXBs identified in previous Chandra and XMM-Newton surveys of the SMC using the AAOmega/2dF fiber-fed spectrograph at the Anglo-Australian Telescope. We find 5 new Be/X-ray binaries (BeXRBs; including a tentative one), by identifying the spectral type of their optical counterparts, and we confirm the spectral classification of an additional 15 known BeXRBs. We compared the spectral types, orbital periods, and eccentricities of the BeXRB populations in the SMC and the Milky Way and we find marginal evidence for difference between the spectral type distributions, but no statistically significant differences for the orbital periods and the eccentricities. Moreover, our search revealed that the well known supergiant B[e] star LHA 115-S 18 (or AzV 154) is associated with the weak X-ray source CXOU J005409.57-724143.5. We provide evidence that the supergiant star LHA 115-S 18 is the optical counterpart of the X-ray source, and we discuss different possibilities of the origin of its low X-ray luminosity (Lx ~ 4 x 10^33 erg/s).

We present results from three nearly simultaneous NuSTAR and Chandra monitoring observations between 2012 Sep 2 and 2012 Nov 16 of local star-forming galaxy NGC 253. The 3-40 keV NuSTAR intensity of the inner 20 arcsec (~400 pc) nuclear region varied by a factor of ~2 across the three monitoring observations. The Chandra data reveal that the nuclear region contains three bright X-ray sources, including a luminous (L2-10 keV ~ few x 10^39 erg/s) point source ~1 arcsec from the dynamical center of the galaxy (within the 3sigma positional uncertainty of the dynamical center); this source drives the overall variability of the nuclear region at energies >3 keV. We make use of the variability to measure the spectra of this single hard X-ray source when it was in bright states. The spectra are well described by an absorbed (NH ~ 1.6 x 10^23 cm^-2) broken power-law model with spectral slopes and break energies that are typical of ultraluminous X-ray sources (ULXs), but not AGN. A previous Chandra observation in 2003 showed a hard X-ray point source of similar luminosity to the 2012 source that was also near the dynamical center (~0.4 arcsec); however, this source was offset from the 2012 source position by ~1 arcsec. We show that the probability of the 2003 and 2012 hard X-ray sources being unrelated is >>99.99% based on the Chandra spatial localizations. Interestingly, the Chandra spectrum of the 2003 source (3-8 keV) is shallower in slope than that of the 2012 hard X-ray source. Its proximity to the dynamical center and harder Chandra spectrum indicate that the 2003 source is a better AGN candidate than any of the sources detected in our 2012 campaign; however, we were unable to rule out a ULX nature for this source. Future NuSTAR and Chandra monitoring would be well equipped to break the degeneracy between the AGN and ULX nature of the 2003 source, if again caught in a high state.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

Using Chandra, XMM-Newton, and optical photometric catalogs we study the young X-ray binary (XRB) populations of the Small Magellanic Cloud. We find that the Be/X-ray binaries (Be-XRBs) are observed in regions with star formation rate bursts ~25-60 Myr ago. The similarity of this age with the age of maximum occurrence of the Be phenomenon (~40 Myr) indicates that the presence of a circumstellar decretion disk plays a significant role in the number of observed XRBs in the 10-100 Myr age range. We also find that regions with strong but more recent star formation (e.g., the Wing) are deficient in Be-XRBs. By correlating the number of observed Be-XRBs with the formation rate of their parent populations, we measure a Be-XRB production rate of ~1 system per 3 x 10^(-3) M$_{\odot}$/yr. Finally, we use the strong localization of the Be-XRB systems in order to set limits on the kicks imparted on the neutron star during the supernova explosion.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

We present a large sample (20 in total) of optical spectra of Small Magellanic Cloud (SMC) High-Mass X-ray Binaries obtained with the 2dF spectrograph at the Anglo-Australian Telescope. All of these sources are found to be Be/X-ray binaries (Be-XRBs), while for 5 sources we present original classifications. Several statistical tests on this expanded sample support previous findings for similar spectral-type distributions of Be-XRBs and Be field stars in the SMC, and of Be-XRBs in the Large Magellanic Cloud and the Milky Way, although this could be the result of small samples. On the other hand, we find that Be-XRBs follow a different distribution than Be stars in the Galaxy, also in agreement with previous studies. In addition, we find similar Be spectral type distributions between the Magellanic Clouds samples. These results reinforce the relation between the orbital period and the equivalent width of the Halpha line that holds for Be-XRBs. SMC Be stars have larger Halpha equivalent widths when compared to Be-XRBs, supporting the notion of circumstellar disk truncation by the compact object.

We present the most likely optical counterparts of 113 X-ray sources detected in our Chandra survey of the central region of the Small Magellanic Cloud (SMC) based on the OGLE-II and MCPS catalogs. We estimate that the foreground contamination and chance coincidence probability are minimal for the bright optical counterparts (corresponding to OB type stars; 35 in total). We propose here for the first time 13 High-Mass X-ray Binaries (HMXBs), of which 4 are Be/X-ray binaries (Be-XRBs), and we confirm the previous classification of 18 Be-XRBs. We estimate that the new candidate Be-XRBs have an age of 15-85 Myr, consistent with the age of Be stars. We also examine the "overabundance" of Be-XRBs in the SMC fields covered by Chandra, in comparison with the Galaxy. In luminosities down to about 10^{34} erg/s, we find that SMC Be-XRBs are 1.5 times more common when compared to the Milky Way even after taking into account the difference in the formation rates of OB stars. This residual excess can be attributed to the lower metallicity of the SMC. Finally, we find that the mixing of Be-XRBs with other than their natal stellar population is not an issue in our comparisons of Be-XRBs and stellar populations in the SMC. Instead, we find indication for variation of the SMC XRB populations on kiloparsec scales, related to local variations of the formation rate of OB stars and slight variation of their age, which results in different relative numbers of Be stars and therefore XRBs.