We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability.

Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows.

We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $\alpha=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.

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.

The launch of the Nuclear Spectroscopic Telescope Array (NuSTAR) heralded a new era of sensitive high energy X-ray spectroscopy for X-ray binaries (XRBs). In this paper we show how multiple physical parameters can be measured from the accretion disk spectrum when the high-energy side of the disk spectrum can be measured precisely using NuSTAR. This immediately makes two exciting developments possible. If the mass and distance of the source are known, the continuum fitting method can be used to calculate the spin and inner disk inclination independently of the iron line fitting method. If the mass and distance are unknown, the two methods can be combined to constrain these values to a narrow region of parameter space. In this paper we perform extensive simulations to establish the reliability of these techniques. We find that with high quality spectra, spin and inclination can indeed be simultaneously measured using the disk spectrum. These measurements are much more precise at higher spin values, where the relativistic effects are stronger. The inclusion of a soft X-ray snapshot observation alongside the NuSTAR data significantly improves the reliability, particularly for lower temperature disks, as it gives a greatly improved measurement of the disk peak. High signal to noise data are not necessary for this, as measuring the peak temperature is relatively easy. We discuss the impact of systematic effects on this technique, and the implications of our results such as robust measurements of accretion disk warps and XRB mass surveys.

Here we study a 1-200 keV energy spectrum of the black hole binary Cygnus X-1 taken with NuSTAR and Suzaku. This is the first report of a NuSTAR observation of Cyg X-1 in the intermediate state, and the observation was taken during the part of the binary orbit where absorption due to the companion's stellar wind is minimal. The spectrum includes a multi-temperature thermal disk component, a cutoff power-law component, and relativistic and non-relativistic reflection components. Our initial fits with publicly available constant density reflection models (relxill and reflionx) lead to extremely high iron abundances (>9.96 and 10.6(+1.6)(-0.9) times solar, respectively). Although supersolar iron abundances have been reported previously for Cyg X-1, our measurements are much higher and such variability is almost certainly unphysical. Using a new version of reflionx that we modified to make the electron density a free parameter, we obtain better fits to the spectrum even with solar iron abundances. We report on how the higher density (n_e = (3.98(+0.12)(-0.25))E20 cm-3) impacts other parameters such as the inner radius and inclination of the disk.

We have discovered heavy obscuration in the dual active galactic nucleus (AGN) in the Swift/Burst Alert Telescope (BAT) source SWIFT J2028.5+2543 using Nuclear Spectroscopic Telescope Array (NuSTAR). While an early XMM-Newton study suggested the emission was mainly from NGC 6921, the superior spatial resolution of NuSTAR above 10 keV resolves the Swift/BAT emission into two sources associated with the nearby galaxies MCG +04-48-002 and NGC 6921 (z = 0.014) with a projected separation of 25.3 kpc (91"). NuSTAR's sensitivity above 10 keV finds both are heavily obscured to Compton-thick levels (N H=(1-2)x10^24 cm-2) and contribute equally to the BAT detection ({L}10-50 {keV}{{int}}= 6x10^42 erg s-1). The observed luminosity of both sources is severely diminished in the 2-10 keV band, illustrating the importance of >10 keV surveys like those with NuSTAR and Swift/BAT. Compared to archival X-ray data, MCG +04-48-002 shows significant variability (>3) between observations. Despite being bright X-ray AGNs, they are difficult to detect using optical emission-line diagnostics because MCG +04-48-002 is identified as a starburst/composite because of the high rates of star formation from a luminous infrared galaxy while NGC 6921 is only classified as a LINER using line detection limits. SWIFT J2028.5+2543 is the first dual AGN resolved above 10 keV and is the second most heavily obscured dual AGN discovered to date in the X-rays other than NGC 6240.

We present NuSTAR and Swift observations of the neutron star Aquila X-1 during the peak of its July 2014 outburst. The spectrum is soft with strong evidence for a broad Fe K\alpha line. Modeled with a relativistically broadened reflection model, we find that the inner disk is truncated with an inner radius of 15+/-3 R_G. The disk is likely truncated by either the boundary layer and/or a magnetic field. Associating the truncated inner disk with pressure from a magnetic field gives an upper limit of B<5+/-2x10^8G. Although the radius is truncated far from the stellar surface, material is still reaching the neutron star surface as evidenced by the X-ray burst present in the t NuSTAR observation.

The black hole candidate XTE J1908+094 went into outburst for the first time since 2003 in October 2013. We report on an observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) and monitoring observations with Swift during the outburst. NuSTAR caught the source in the soft state: the spectra show a broad relativistic iron line, and the light curves reveal a ~40 ks flare with the count rate peaking about 40% above the non-flare level and with significant spectral variation. A model combining a multi-temperature thermal component, a power-law, and a reflection component with an iron line provides a good description of the NuSTAR spectrum. Although relativistic broadening of the iron line is observed, it is not possible to constrain the black hole spin with these data. The variability of the power-law component, which can also be modeled as a Comptonization component, is responsible for the flux and spectral change during the flare, suggesting that changes in the corona (or possibly continued jet activity) are the likely cause of the flare.

We report on a 10 ks simultaneous Chandra/HETG-NuSTAR observation of the Bursting Pulsar, GRO J1744-28, during its third detected outburst since discovery and after nearly 18 years of quiescence. The source is detected up to 60 keV with an Eddington persistent flux level. Seven bursts, followed by dips, are seen with Chandra, three of which are also detected with NuSTAR. Timing analysis reveals a slight increase in the persistent emission pulsed fraction with energy (from 10% to 15%) up to 10 keV, above which it remains constant. The 0.5-70 keV spectra of the persistent and dip emission are the same within errors, and well described by a blackbody (BB), a power-law with an exponential rolloff, a 10 keV feature, and a 6.7 keV emission feature, all modified by neutral absorption. Assuming that the BB emission originates in an accretion disc, we estimate its inner (magnetospheric) radius to be about 4x10^7 cm, which translates to a surface dipole field B~9x10^10 G. The Chandra/HETG spectrum resolves the 6.7 keV feature into (quasi-)neutral and highly ionized Fe XXV and Fe XXVI emission lines. XSTAR modeling shows these lines to also emanate from a truncated accretion disk. The burst spectra, with a peak flux more than an order of magnitude higher than Eddington, are well fit with a power-law with an exponential rolloff and a 10~keV feature, with similar fit values compared to the persistent and dip spectra. The burst spectra lack a thermal component and any Fe features. Anisotropic (beamed) burst emission would explain both the lack of the BB and any Fe components.

We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $\sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} \times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40\pm0.20 \times 10^{39}$ erg s$^{-1}$ for IC 342 X-2, so that both were observed in a similar, luminous state. Both sources have a high absorbing column in excess of the Galactic value. Neither source has a spectrum consistent with a black hole binary in low/hard state, and both ULXs exhibit strong curvature in their broadband X-ray spectra. This curvature rules out models that invoke a simple reflection-dominated spectrum with a broadened iron line and no cutoff in the illuminating power-law continuum. X-ray spectrum of IC 342 X-1 can be characterized by a soft disk-like black body component at low energies and a cool, optically thick Comptonization continuum at high energies, but unique physical interpretation of the spectral components remains challenging. The broadband spectrum of IC 342 X-2 can be fit by either a hot (3.8 keV) accretion disk, or a Comptonized continuum with no indication of a seed photon population. Although the seed photon component may be masked by soft excess emission unlikely to be associated with the binary system, combined with the high absorption column, it is more plausible that the broadband X-ray emission arises from a simple thin blackbody disk component. Secure identification of the origin of the spectral components in these sources will likely require broadband spectral variability studies.

We present an analysis of a short NuSTAR observation of the stellar-mass black hole and low-mass X-ray binary 4U 1630$-$472. Reflection from the inner accretion disk is clearly detected for the first time in this source, owing to the sensitivity of NuSTAR. With fits to the reflection spectrum, we find evidence for a rapidly spinning black hole, $a_{*}=0.985^{+0.005}_{-0.014}$ ($1\sigma$ statistical errors). However, archival data show that the source has relatively low radio luminosity. Recently claimed relationships between jet power and black hole spin would predict either a lower spin or a higher peak radio luminosity. We also report the clear detection of an absorption feature at $7.03\pm0.03$ keV, likely signaling a disk wind. If this line arises in dense, moderately ionized gas ($\log\xi=3.6^{+0.2}_{-0.3}$) and is dominated by He-like Fe XXV, the wind has a velocity of $v/c=0.043^{+0.002}_{-0.007}$ (12900$^{+600}_{-2100}$ km s$^{-1}$). If the line is instead associated with a more highly ionized gas ($\log\xi=6.1^{+0.7}_{-0.6}$), and is dominated by Fe XXVI, evidence of a blue-shift is only marginal, after taking systematic errors into account. Our analysis suggests the ionized wind may be launched within 200$-$1100 Rg, and may be magnetically driven.