Search

The discovery of pulsations in (at least) six ultraluminous X-ray sources (ULXs) has shown that neutron stars can accrete at (highly) super-Eddington rates, challenging the standard accretion theories. M51 ULX-7, with a spin signal of $P\simeq2.8$ s, is the pulsating ULX (PULX) with the shortest known orbital period ($P_\mathrm{orb}\simeq2$ d) and has been observed multiple times by XMM-Newton, Chandra, and NuSTAR. We report on the timing and spectral analyses of three XMM-Newton observations of M51 ULX-7 performed between the end of 2021 and the beginning of 2022, together with a timing re-analysis of XMM-Newton, Chandra, and NuSTAR archival observations. We investigated the spin signal by applying accelerated search techniques and studied the power spectrum through the fast Fourier transform, looking for (a)periodic variability in the source flux. We analysed the energy spectra of the 2021-2022 observations and compared them to the older ones. We report the discovery of a recurrent, significant ($>$3$\sigma$) broad complex at mHz frequencies in the power spectra of M51 ULX-7. We did not detect the spin signal, setting a 3$\sigma$ upper limit on the pulsed fraction of $\lesssim10\%$ for the single observation. The complex is significantly detected also in five Chandra observations performed in 2012. M51 ULX-7 represents the second PULX for which we have a significant detection of mHz-QPOs at super-Eddington luminosities. These findings suggest that one should avoid using the observed QPO frequency to infer the mass of the accretor in a ULX. The absence of spin pulsations when the broad complex is detected suggests that the mechanism responsible for the aperiodic modulation also dampens the spin signal's pulsed fraction. If true, this represents an additional obstacle in the detection of new PULXs, suggesting an even larger occurrence of PULXs among ULXs., Comment: 14 pages (12 main text + Appendix), 5 figures. Accepted for publication on A&A

We report on the orbit of the binary system powering the most extreme ultraluminous X-ray pulsar known to date: NGC 5907 ULX-1 (hereafter ULX1). ULX1 has been the target of a substantial multi-instrument campaign, mainly in the X-ray band, but no clear counterparts are known in other bands. Although ULX1 is highly variable and pulsations can be transient (regardless of the source flux), the timing data collected so far allow us to investigate the orbit of this system. We find an orbital period $P_{orb}=5.7^{+0.1}_{-0.6}\text{ d}$ and a projected semi-axis $A_1 =3.1^{+0.8}_{-0.9}\text{ lts}$. The most likely ephemeris is: $P_{orb}=5.6585(6)\text{ d}$, $A_1 = 3.1(4)\text{ lts}$, and the epoch of ascending nodes passage is: $T_{asc} = 57751.37(5)\text{ MJD}$. However, there are 6 similar solutions, acceptable within $3\,\sigma$. We find further indications that ULX1 is a high-mass X-ray binary. This implies that we are observing its orbit face-on, with an inclination $<5\text{ deg}$., Comment: 11 pages, 2 tables, 6 figures

We observed the nearby and relatively understudied ultraluminous X-ray source (ULX) NGC 4190 ULX-1 jointly with NICER and NuSTAR to investigate its broadband spectrum, timing properties, and spectral variation over time. We found NGC 4190 ULX-1 to have a hard spectrum characterized by two thermal components (with temperatures ~0.25keV and ~1.6keV) and a high-energy excess typical of the ULX population, although the spectrum turns over at an unusually low energy. While no pulsations were detected, (with pulsed fraction 3-sigma upper limits of 16% for NICER and 35% for NuSTAR), the source shows significant stochastic variability and the covariance spectrum indicates the presence of a high-energy cut-off power-law component, potentially indicative of an accretion column. Additionally, when fitting archival XMM-Newton data with a similar model, we find that the luminosity-temperature evolution of the hot thermal component follows the behavior of a super-Eddington slim disk though the expected spectral broadening for such a disk is not seen, suggesting that the inner accretion disk may be truncated by a magnetic field. Therefore, despite the lack of detected pulsations, there is tantalizing evidence for NGC 4190 ULX-1 being a candidate neutron star accretor, although further broadband observations will be required to confirm this behavior., Comment: 14 pages, 9 figures, 2 tables, accepted for publication in The Astrophysical Journal

We present radio observations of 24 confirmed and candidate strongly lensed quasars identified by the Gaia Gravitational Lenses (GraL) working group. We detect radio emission from 8 systems in 5.5 and 9 GHz observations with the Australia Telescope Compact Array (ATCA), and 12 systems in 6 GHz observations with the Karl G. Jansky Very Large Array (VLA). The resolution of our ATCA observations is insufficient to resolve the radio emission into multiple lensed images, but we do detect multiple images from 11 VLA targets. We have analysed these systems using our observations in conjunction with existing optical measurements, including measuring offsets between the radio and optical positions, for each image and building updated lens models. These observations significantly expand the existing sample of lensed radio quasars, suggest that most lensed systems are detectable at radio wavelengths with targeted observations, and demonstrate the feasibility of population studies with high resolution radio imaging.

Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion -- hence their progenitors, lifetimes, and future evolution -- is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational wave sources. Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities. HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR. We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies., Comment: 24 pages, 9 figures. Accepted for publication in Frontiers of Astronomy and Space Science. Minor corrections included

Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discovered in Swift/XRT observations. Swift monitoring of these sources have provided well sampled lightcurves, allowing for us to model the lightcurves with the disk instability model of Hameury & Lasota (2020) which implies durations of 60-400 days and that the mass accretion rate through the disk is close to or greater than the Eddington rate. Of the three source regions with prior HST imaging, color magnitude diagrams of the potential stellar counterparts show varying ages of the possible stellar counterparts. Our estimation of the rates of these sources in these three galaxies is 0.4-1.3 year$^{-1}$. We find that while persistent ULXs dominate the high end of galaxy luminosity functions, the number of systems that produce ULX luminosities are likely dominated by transient sources., Comment: Accepted for publication in ApJ

We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs.

We present a high-density relativistic reflection analysis of 21 spectra of six black hole X-ray binaries in the hard state with data from \textit{NuSTAR} and \textit{Swift}. We find that 76\% of the observations in our sample require a disk density higher than the 10$^{15}$~cm$^{-3}$ assumed in the previous reflection analysis. Compared with the measurements from active galactic nuclei, stellar mass black holes have higher disk densities. Our fits indicate that the inner disk radius is close to the innermost stable circular orbit in the hard state. The coronal temperatures are significantly lower than the prediction of a purely thermal plasma, which can be explained with a hybrid plasma model. If the disk density is fixed at 10$^{15}$~cm$^{-3}$, the disk ionization parameter would be overestimated while the inner disk radius is unaffected., Comment: 17 pages, 9 figures, accepted by ApJ

The nuclear transient AT2019cuk/Tick Tock/SDSS J1430+2303 has been suggested to harbor a supermassive black hole (SMBH) binary near coalescence. We report results from high-cadence NICER X-ray monitoring with multiple visits per day from January-August 2022, as well as continued optical monitoring during the same time period. We find no evidence of periodic/quasi-periodic modulation in the X-ray, UV, or optical bands, however we do observe exotic hard X-ray variability that is unusual for a typical AGN. The most striking feature of the NICER light curve is repetitive hard (2-4 keV) X-ray flares that result in distinctly harder X-ray spectra compared to the non-flaring data. In its non-flaring state, AT2019cuk looks like a relatively standard AGN, but it presents the first case of day-long, hard X-ray flares in a changing-look AGN. We consider a few different models for the driving mechanism of these hard X-ray flares, including: (1) corona/jet variability driven by increased magnetic activity, (2) variable obscuration, and (3) self-lensing from the potential secondary SMBH. We prefer the variable corona model, as the obscuration model requires rather contrived timescales and the self-lensing model is difficult to reconcile with a lack of clear periodicity in the flares. These findings illustrate how important high-cadence X-ray monitoring is to our understanding of the rapid variability of the X-ray corona and necessitate further high-cadence, multi-wavelength monitoring of changing-look AGN like AT2019cuk to probe the corona-jet connection., Comment: 23 pages, 13 figures, 1 table, accepted for publication in ApJL

The stability criteria of rapid mass transfer and common envelope evolution are fundamental in binary star evolution. They determine the mass, mass ratio and orbital distribution of many important systems, such as X-ray binaries, Type Ia supernovae and merging gravitational wave sources. We use our adiabatic mass-loss model to systematically survey the intermediate-mass stars' thresholds for dynamical-timescale mass transfer. The impact of metallicity on the stellar responses and critical mass ratios is explored. Both tables ($Z=0.001$) and fitting formula ($Z=0.001$ and $Z=0.02$) of critical mass ratios of intermediate-mass stars are provided. An application of our results to intermediate-mass X-ray binaries (IMXBs) is discussed. We find that the predicted upper limit to mass ratios, as a function of orbital period, is consistent with the observed IMXBs that undergo thermal or nuclear timescale mass transfer. According to the observed peak X-ray luminosity $L_\mathrm{X}$, we predict the range of $L_\mathrm{X}$ for IMXBs as a function of the donor mass and the mass transfer timescale., Comment: 16 pages, 19 figures, 3 tables, accepted for publication in The Astrophysical Journal

Ultra-luminous X-ray sources (ULXs) are the most extreme members of the X-ray binary population, exhibiting X-ray luminosities that can surpass the 10^39 erg/s threshold (by orders of magnitude). They are mainly seen in external galaxies and are most preferentially found in star-forming galaxies with lower metallicities. The vast majority of these systems are now understood to be powered by super-Eddington accretion of matter onto stellar-mass compact objects (black holes and neutron stars). This is driven by the discovery of coherent pulsations, cyclotron lines and powerful winds in members of the ULX population. The latter was possible thanks to high-resolution X-ray spectrometers such as those aboard XMM-Newton. ULX winds carry a huge amount of power owing to their relativistic speeds (0.1-0.3 c) and are likely responsible for the ~100 pc superbubbles observed around many ULXs. The winds also regulate the amount of matter that can reach the central accretor. Their study is, therefore, essential to understanding how quickly compact objects can grow and how strong their feedback onto the surrounding medium can be. This may also be relevant to understand supermassive black hole growth, particularly in the early Universe. Here we provide an overview on ULX phenomenology, highlight some recent exciting results, and show how future missions such as XRISM and ATHENA will drive further significant progress in this field., Comment: Invited review chapter for the book High-Resolution X-Ray Spectroscopy: Instrumentation, Data Analysis, and Science (Eds. C. Bambi and J. Jiang, Springer Singapore, expected in 2023)

4U 1957+11 is a black hole candidate system that has been in a soft X-ray spectral state since its discovery. We present analyses of recent joint NICER and NuSTAR spectra, which are extremely well-described by a highly inclined disk accreting into a near maximally spinning black hole. Owing to the broad X-ray coverage of NuSTAR the fitted spin and inclination are strongly constrained for our hypothesized disk models. The faintest spectra are observed out to 20 keV, even though their hard tail components are almost absent when described with a simple corona. The hard tail increases with luminosity, but shows clear two track behavior with one track having appreciably stronger tails. The disk spectrum color-correction factor is anti-correlated with the strength of the hard tail (e.g., as measured by the Compton $y$ parameter). Although the spin and inclination parameters are strongly constrained for our chosen model, the mass and distance are degenerate parameters. We use our spectral fits, along with a theoretical prior on color-correction, an observational prior on likely fractional Eddington luminosity, and an observational prior on distance obtained from Gaia studies, to present mass and distance contours for this system. The most likely parameters, given our presumed disk model, suggest a 4.6 $\mathrm{M_\odot}$ black hole at 7.8 kpc observed at luminosities ranging from $\approx 1.7\%$--$9\%$ of Eddington. This would place 4U 1957+11 as one of the few actively accreting sources within the `mass gap' of ${\approx} 2$--$5\,\mathrm{M_\odot}$ where there are few known massive neutron stars or low mass black holes. Higher mass and distance, however, remain viable., Comment: 21 pages. To be published in the Astrophysical Journal

We present X-ray, UV, optical, and radio observations of the nearby ($\approx78$ Mpc) tidal disruption event (TDE) AT2021ehb/ZTF21aanxhjv during its first 430 days of evolution. AT2021ehb occurs in the nucleus of a galaxy hosting a $\approx 10^{7}\,M_\odot$ black hole ($M_{\rm BH}$ inferred from host galaxy scaling relations). High-cadence Swift and NICER monitoring reveals a delayed X-ray brightening. The spectrum first undergoes a gradual ${\rm soft }\rightarrow{\rm hard}$ transition and then suddenly turns soft again within 3 days at $\delta t\approx 272$ days during which the X-ray flux drops by a factor of ten. In the joint NICER+NuSTAR observation ($\delta t =264$ days, harder state), we observe a prominent non-thermal component up to 30 keV and an extremely broad emission line in the iron K band. The bolometric luminosity of AT2021ehb reaches a maximum of $6.0^{+10.4}_{-3.8}\% L_{\rm Edd}$ when the X-ray spectrum is the hardest. During the dramatic X-ray evolution, no radio emission is detected, the UV/optical luminosity stays relatively constant, and the optical spectra are featureless. We propose the following interpretations: (i) the ${\rm soft }\rightarrow{\rm hard}$ transition may be caused by the gradual formation of a magnetically dominated corona; (ii) hard X-ray photons escape from the system along solid angles with low scattering optical depth ($\sim\,$a few) whereas the UV/optical emission is likely generated by reprocessing materials with much larger column density -- the system is highly aspherical; (iii) the abrupt X-ray flux drop may be triggered by the thermal-viscous instability in the inner accretion flow leading to a much thinner disk., Comment: 35 pages, 21 figures, accepted by ApJ

We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find strong evidence for a change in period over the time that NGC 925 ULX-3 has been observed, although the source may have been in a much lower flux state when first observed with Chandra in 2005. We do not detect pulsations, and we place an upper limit on the pulsed fraction of ~40% in the XMM-Newton band, consistent with some previous pulsation detections at low energies in other ULXs. The source exhibits a typical ULX spectrum that turns over in the NuSTAR band and can be fitted using two thermal components. These components have a high temperature ratio that may indicate the lack of extreme inner disk truncation by a magnetar-level magnetic field. We examine the implications for a number of different models for superorbital periods in ULXs, finding that a neutron star with a magnetic field of ~10^12 G may be plausible for this source. The future detection of pulsations from this source would allow for the further testing and constraining of such models., Comment: 13 pages, 6 figures, 2 tables. Accepted for publication in ApJ

We report on VLT/FORS2 imaging polarimetry observations in the $R_{\rm special}$ band of WISE J011601.41-050504.0 (W0116-0505), a heavily obscured hyper-luminous quasar at $z=3.173$ classified as a Hot, Dust-Obscured Galaxy (Hot DOG) based on its mid-IR colors. Recently, Assef et al. (2020) identified W0116-0505 as having excess rest-frame optical/UV emission, and concluded this excess emission is most likely scattered light from the heavily obscured AGN. We find that the broad-band rest-frame UV flux is strongly linearly polarized (10.8$\pm$1.9\%, with a polarization angle of 74$\pm$9~deg), confirming this conclusion. We analyze these observations in the context of a simple model based on scattering either by free electrons or by optically thin dust, assuming a classical dust torus with polar openings. Both can replicate the degree of polarization and the luminosity of the scattered component for a range of geometries and column densities, but we argue that optically thin dust in the ISM is the more likely scenario. We also explore the possibility that the scattering medium corresponds to an outflow recently identified for W0116-0505. This is a feasible option if the outflow component is bi-conical with most of the scattering occurring at the base of the receding outflow. In this scenario the quasar would still be obscured even if viewed face on, but might appear as a reddened type 1 quasar once the outflow has expanded. We discuss a possible connection between blue-excess Hot DOGs, extremely red quasars (ERQs), reddened type 1 quasars, and unreddened quasars that depends on a combination of evolution and viewing geometry., Comment: 19 pages, 10 figures, 2 tables. Resubmitted to ApJ after first round of referee comments

We present an analysis of the NuSTAR and Swift spectra of the black hole candidate MAXI J1813-095 in a failed-transition outburst in 2018. The NuSTAR observations show evidence of reflected emission from the inner region of the accretion disc. By modelling the reflection component in the spectra, we find a disc inner radius of $R_{\rm in}<7$ $r_{\rm g}$. This result suggests that either a slightly truncated disc or a non-truncated disc forms at a few percent of the Eddington limit in MAXI J1813-095. Our best-fit reflection models indicate that the geometry of the innermost accretion remains consistent during the period of NuSTAR observations. The spectral variability of MAXI J1813-095 from multi-epoch observations is dominated by the variable photon index of the Comptonisation emission., Comment: 9 pages, 12 figures, accepted for publication in MNRAS

We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3\sigma$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monitoring reveals a highly variable source changing by over an order of magnitude from maximum to minimum. The origin of the absorption line appears likely to be by highly ionized iron with a blue shift of 0.19$c$, indicating an ultrafast outflow (UFO). However, the large equivalent width of the line (EW$=-0.22^{+0.08}_{-0.09}$ keV) paired with the lack of other absorption lines detected are difficult to reconcile with models. An alternative explanation is that the line is due to a cyclotron resonance scattering feature (CRSF), produced by the interaction of X-ray photons with the powerful magnetic field of a neutron star., Comment: Accepted for publication in ApJ

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenarios. In this Paper, we follow the orbit of the neutron star for seven years, measure the decay of the orbit ($\dot{P}_{orb}/{P}_{orb}\approx-8\cdot10^{-6}\mathrm{yr}^{-1}$), and argue that this orbital decay is driven by extreme mass transfer of more than 150 times the mass transfer limit set by the Eddington luminosity. If this is true, the mass available to the accretor is more than enough to justify its luminosity, with no need for beaming. This also strongly favors models where the accretor is a highly-magnetized neutron star., Comment: 18 pages, 6 figures, accepted for publication in ApJ

M51 ULX7 is among a small group of known ultraluminous X-ray pulsars (ULXP). The neutron star powering the source has a spin period of 2.8s, orbits its companion star with a period of 2 days, and a super-orbital period of 38 days is evident in its X-ray lightcurve. Here we present NuSTAR and XMM-Newton data on the source from 2019 obtained when the source was near its peak brightness. We detect the pulsations, having spun up at a rate of 3$\pm0.5\times10^{-10}$ s s$^{-1}$ since they were previously detected in 2018. The data also provide the first high-quality broadband spectrum of the source. We find it to be very similar to that of other ULXPs, with two disk-like components, and a high energy tail. When combined with XMM-Newton data obtained in 2018, we explore the evolution of the spectral components with super-orbital phase, finding that the luminosity of the hotter component drives the super-orbital flux modulation. The inclination the disk components appear to change with phase, which may support the idea that these super-orbital periods are caused by disk precession. We also reexamine the super-orbital period with 3 years of Swift/XRT monitoring, finding that the period is variable, increasing from 38.2$\pm0.5$ days in 2018--2019 to 44.2$\pm0.9$ days in 2020--2021, which rules out alternative explanations for the super-orbital period., Comment: Accepted for publication in ApJ

Ultra-fast outflows (UFOs) have been detected in the high-quality X-ray spectra of a number of active galactic nuclei (AGN) with fairly high accretion rates and are thought to significantly contribute to the AGN feedback. After a decade of dedicated study, their launching mechanisms and structure are still not well understood, but variability techniques may provide useful constraints. In this work, therefore, we perform a flux-resolved X-ray spectroscopy on a highly accreting and variable NLS1 AGN, 1H 0707-495, using all archival XMM-Newton observations to study the structure of the UFO. We find that the wind spectral lines weaken at higher luminosities, most likely due to an increasing ionization parameter as previously found in a few similar sources. Instead, the velocity is anticorrelated with the luminosity, which is opposite to the trend observed in the NLS1 IRAS 13224-3809. Furthermore, the detection of the emission lines, which are not observed in IRAS 13224-3809, indicates a wind with a larger opening angle in 1H 0707-495, presumably due to a higher accretion rate. The emitting gas is found to remain broadly constant with the luminosity. We describe the variability of the wind with a scenario where the strong radiation extends the launch radius outwards and shields the outer emitting gas, similarly to super-Eddington compact objects, although other possible explanations are discussed. Our work provides several hints for a multi-phase outflow in 1H 0707-495., Comment: 20 pages, 9 figures, accepted for publication in MNRAS

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA., Comment: 11 pages with 11 pages of appendices. Submitted to the Astrophysical Journal. Comments welcome

We present results for the first observed outburst from the transient X-ray binary source MAXI J0637-430. This study is based on eight observations from the Nuclear Spectroscopic Telescope Array (NuSTAR) and six observations from the Neil Gehrels Swift Observatory X-Ray Telescope (Swift/XRT) collected from 2019 November 19 to 2020 April 26 as the 3-79 keV source flux declined from 8.2e-10 to 1.4e-12 erg/cm^2/s. We see the source transition from a soft state with a strong disk-blackbody component to a hard state dominated by a power-law or thermal Comptonization component. NuSTAR provides the first reported coverage of MAXI J0637-430 above 10 keV, and these broadband spectra show that a two-component model does not provide an adequate description of the soft state spectrum. As such, we test whether blackbody emission from the plunging region could explain the excess emission. As an alternative, we test a reflection model that includes a physical Comptonization continuum. Finally, we also test a spectral component based on reflection of a blackbody illumination spectrum, which can be interpreted as a simple approximation to the reflection produced by returning disk radiation due to the bending of light by the strong gravity of the black hole. We discuss the physical implications of each scenario and demonstrate the value of constraining the source distance., Comment: 17 pages, 10 figures, accepted to The Astrophysical Journal (ApJ)

We present a detailed spectral analysis of the joint XMM-Newton and NuSTAR observations of the active galactic nuclei (AGN) in the Seyfert 1.5 Galaxy ESO 362-G18. The broadband ($0.3\mbox{--}79$ keV) spectrum shows the presence of a power-law continuum with a soft excess below $2$ keV, iron K$\alpha$ emission ($\sim 6.4$ keV), and a Compton hump (peaking at $\sim 20$ keV). We find that the soft excess can be modeled by two different possible scenarios: a warm ($kT_\mathrm{e}\sim0.2$ keV) and optically thick ($\tau\sim34$) Comptonizing corona; or with relativistically-blurred reflection off a high-density ($\log{[n_\mathrm{e}/\mathrm{cm}^{-3}]}>18.3$) inner disk. These two models cannot be easily distinguished solely from their fit statistics. However, the low temperature ($kT_\mathrm{e}\sim20$ keV) and the thick optical depth ($\tau\sim5$) of the hot corona required by the warm corona scenario are uncommon for AGNs. We also fit a 'hybrid' model, which includes both disk reflection and a warm corona. Unsurprisingly, as this is the most complex of the models considered, this provides the best fit, and more reasonable coronal parameters. In this case, the majority of the soft excess flux arises in the warm corona component. However, based on recent simulations of warm coronae, it is not clear whether such a structure can really exist at the low accretion rates relevant for ESO 362-G18 ($\dot{m}\sim0.015$). This may therefore argue in favor of a scenario in which the soft excess is instead dominated by the relativistic reflection. Based on this model, we find that the data would require a compact hot corona ($h\sim3\,R_\mathrm{Horizon}$) around a highly spinning ($a_\star>0.927$) black hole., Comment: 18 pages, 9 figures, Accepted by ApJ

KUG 1141+371 is a Seyfert 1 galaxy that shows a simultaneous flux increase in the optical and UV bands in the past decade. For instance, the latest Swift observation in 2019 shows that the UVW2 flux of the AGN in KUG 1141+371 has increased by over one order of magnitude since 2009. Meanwhile, the soft X-ray flux of KUG 1141+371 also shows a steady increase by one order of magnitude since 2007. The significant multi-wavelength luminosity change is likely due to a boost of mass accretion rate from approximately 0.6% of the Eddington limit to 3.2%, assuming a black hole mass of $10^{8}M_{\odot}$. In this work, we conduct detailed multi-epoch X-ray spectral analysis focusing on the variability of the X-ray continuum emission and the puzzling soft excess emission. In addition, our SED models also suggest a simultaneous increase of disc temperature and a decreasing inner disc radius along with the increasing accretion rate. Finally, we discuss possible connection between KUG 1141+371 and black hole transients in outburst., Comment: 18 pages, 14 figures, accepted by MNRAS

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission is the first focusing X-ray telescope in the hard X-ray (3-79 keV) band. Among the phenomena that can be studied in this energy band, some require high time resolution and stability: rotation-powered and accreting millisecond pulsars, fast variability from black holes and neutron stars, X-ray bursts, and more. Moreover, a good alignment of the timestamps of X-ray photons to UTC is key for multi-instrument studies of fast astrophysical processes. In this Paper, we describe the timing calibration of the NuSTAR mission. In particular, we present a method to correct the temperature-dependent frequency response of the on-board temperature-compensated crystal oscillator. Together with measurements of the spacecraft clock offsets obtained during downlinks passes, this allows a precise characterization of the behavior of the oscillator. The calibrated NuSTAR event timestamps for a typical observation are shown to be accurate to a precision of ~65 microsec., Comment: 13 pages, 10 figures. Comments welcome

We present results from the Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the new black hole X-ray binary candidate MAXI J1631-479 at two epochs during its 2018-2019 outburst, which caught the source in a disk dominant state and a power-law dominant state. Strong relativistic disk reflection features are clearly detected, displaying significant variations in the shape and strength of the broad iron emission line between the two states. Spectral modeling of the reflection spectra reveals that the inner radius of the optically-thick accretion disk evolves from $<1.9$ $r_{\rm g}$ to $12\pm1$ $r_{\rm g}$ (statistical errors at 90% confidence level) from the disk dominant to the power-law dominant state. Assuming in the former case that the inner disk radius is consistent with being at the ISCO, we estimate a black hole spin of $a^*>0.94$. Given that the bolometric luminosity is similar in the two states, our results indicate that the disk truncation observed in MAXI J1631-479 in the power-law dominant state is unlikely to be driven by a global variation in the accretion rate. We propose that it may instead arise from local instabilities in the inner edge of the accretion disk at high accretion rates. In addition, we find an absorption feature in the spectra centered at $7.33\pm0.03$ keV during the disk dominant state, which is evidence for a rare case that an extremely fast disk wind ($v_{\rm out}=0.067^{+0.001}_{-0.004}~c$) is observed in a low-inclination black hole binary, with the viewing angle of $29\pm1^{\circ}$ as determined by the reflection modeling., Comment: 12 pages, 6 figures, accepted for publication in ApJ

We present results from NuSTAR and XMM-Newton observations of the new black hole X-ray binary MAXI J1820+070 at low accretion rates (below 1% of the Eddington luminosity). We detect a narrow Fe K$\alpha$ emission line, in contrast to the broad and asymmetric Fe K$\alpha$ line profiles commonly present in black hole binaries at high accretion rates. The narrow line, with weak relativistic broadening, indicates that the Fe K$\alpha$ line is produced at a large disk radius. Fitting with disk reflection models assuming standard disk emissivity finds a large disk truncation radius (a few tens to a few hundreds of gravitational radii, depending on the disk inclination). In addition, we detect a quasi-periodic oscillation (QPO) varying in frequency between $11.6\pm0.2$~mHz and $2.8\pm0.1$~mHz. The very low QPO frequencies suggest a large size for the optically-thin Comptonization region according to the Lense-Thirring precession model, supporting that the accretion disk recedes from the ISCO and is replaced by advection-dominated accretion flow at low accretion rates. We also discuss the possibility of an alternative accretion geometry that the narrow Fe K$\alpha$ line is produced by a lamppost corona with a large height illuminating the disk., Comment: 8 pages, 5 figures, accepted for publication in ApJ

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

M82 hosts two well-known ultraluminous X-ray sources (ULXs). X-1, an intermediate-mass black hole (IMBH) candidate, and X-2, an ultraluminous X-ray pulsar (ULXP). Here we present a broadband X-ray spectral analysis of both sources based on ten observations made simultaneously with Chandra and NuSTAR. Chandra provides the high spatial resolution to resolve the crowded field in the 0.5--8 keV band, and NuSTAR provides the sensitive hard X-ray spectral data, extending the bandpass of our study above 10 keV. The observations, taken in the period 2015--2016, cover a period of flaring from X-1, allowing us to study the spectral evolution of this source with luminosity. During four of these observations, X-2 was found to be at a low flux level, allowing an unambiguous view of the emission from X-1. We find that the broadband X-ray emission from X-1 is consistent with that seen in other ULXs observed in detail with NuSTAR, with a spectrum that includes a broadened disk-like component and a high-energy tail. We find that the luminosity of the disk scales with inner disk temperature as L~T^-3/2 contrary to expectations of a standard accretion disk and previous results. These findings rule out a thermal state for sub-Eddington accretion and therefore do not support M82 X-1 as an IMBH candidate. We also find evidence that the neutral column density of the material in the line of sight increases with L$_X$, perhaps due to an increased mass outflow with accretion rate. For X-2, we do not find any significant spectral evolution, but we find the spectral parameters of the phase-averaged broadband emission are consistent with the pulsed emission at the highest X-ray luminosities., Comment: In press at ApJ

X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel path between the direct emission from a hot X-ray corona and the reprocessed emission from the inner edge of the accretion disc. However, there is an inherent degeneracy in the reverberation signal between black hole mass, inner disc radius and height of the illuminating corona above the disc. Here, we use a long X-ray observation of the highly-variable active galaxy, IRAS 13224-3809, to track the reverberation signal as the system evolves on timescales of a day. With the inclusion of all the relativistic effects, modelling reveals that the height of the X-ray corona increases with increasing luminosity, providing a dynamic view of the inner accretion region. This simultaneous modelling allows us to break the inherent degeneracies and obtain an independent timing-based estimate for the mass and spin of the black hole. The uncertainty on black hole mass is comparable to the leading optical reverberation method, making X-ray reverberation a powerful technique, particularly for sources with low optical variability., Comment: Accepted Nat. Ast. version

Swift J1858.6-0814 was discovered by Swift-BAT on October 25, 2018. Here we report on the first follow-up NuSTAR observation of the source, which shows variability spanning two orders of magnitude in count rate on timescales of ~10-100 s. The power-spectrum of the source does not show any quasi-periodic oscillations or periodicity, but has a large fractional rms amplitude of 147%$\pm3$%, exhibiting a number of large flares throughout the observation. The hardness ratio (defined as $R_{10-79 \rm keV}/R_{3-10 \rm keV}$) of the flares tends to be soft, while the source spans a range of hardness ratios during non-flaring periods. The X-ray spectrum of the source shows strong reflection features, which become more narrow and peaked during the non-flaring intervals. We fit an absorbed relativistic reflection model to the source spectra to place physical constraints on the system. Most notably, we find that the source exhibits a large and varying intrinsic absorbing column density ($N_{\rm H}=1.4-4.2\times10^{23}$ cm$^{-2}$). This large intrinsic absorption is further supported by the energy spectra extracted from two flares observed simultaneously by NuSTAR and NICER. We find that the inner accretion disk of the source has a low inclination, $i<29^{\circ}$ ( 3$\sigma$ upper-limit), while the iron abundance in the disk is close to solar, $A_{\rm Fe}=1.0\pm0.3$. We set a 90% confidence upper limit on the inner radius of the accretion disk of $r_{\rm in}<8 \ r_{\rm ISCO}$, and, by fixing $r_{\rm in}$ to be at $r_{\rm ISCO}$, a 90% confidence lower-limit on the spin of the black hole of $a^{*}>0.0$. Lastly, we compare the properties of Swift J1858.6$-$0814 to those of V404 Cygni and V4641 Sgr, which both show rapid flaring and a strong and variable absorption., Comment: Accepted for publication in the Astrophysical Journal

We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$ erg s$^{-1}$. With an additional Chandra non-detection in 2005, this object demonstrates a high dynamic range of flux of factor >26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with $\Gamma=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $\Gamma=1.8^{+0.2}_{-0.1}$ and also well-fitted with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs., Comment: 9 pages, 3 figures, 1 table. Accepted for publication in ApJ; updated to accepted version with minor revisions

We present a detailed analysis of the spectral properties of the black hole transient GRS 1716-249, based on the archival Swift and NuSTAR observations taken during the outburst of this source in 2016-2017. The first six NuSTAR observations show that the source is in a canonical hard state, where the spectrum is dominated by a power-law continuum. The seventh NuSTAR observation is taken during the intermediate state where both a disc thermal component and a power-law continuum are shown. All of our observations show a broad emission line feature in the iron band and a Compton hump above 10 keV. We model the broad band spectra using a high density disc reflection model, where the soft X-ray emission in the hard state is interpreted as part of the disc reflection component. This model enables us to constrain the disc density parameter of GRS 1716-249 in the range of $10^{19}$-$10^{20}$ cm$^{-3}$. We only obtain an upper limit of the inner disc radius using high density disc reflection spectroscopy and the results indicate either a non-truncated disc or a slightly truncated disc with $R_{\rm in}<20r_{\rm g}$., Comment: 10 pages, 7 figures, accepted by MNRAS

We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event., Comment: Submitted to ApJ

We present results from modeling the broadband X-ray spectrum of the Type 1 AGN 2MASX J193013.80+341049.5 using NuSTAR, Swift and archival XMM-Newton observations. We find this source to be highly X-ray obscured, with column densities exceeding 10$^{23}$ cm$^{-2}$ across all epochs of X-ray observations, spanning an 8 year period. However, the source exhibits prominent broad optical emission lines, consistent with an unobscured Type 1 AGN classification. We fit the X-ray spectra with both phenomenological reflection models and physically-motivated torus models to model the X-ray absorption. We examine the spectral energy distribution of this source and investigate some possible scenarios to explain the mismatch between X-ray and optical classifications. We compare the ratio of reddening to X-ray absorbing column density ($E_{B-V}/N_{H}$) and find that 2MASX J193013.80+341049.5 likely has a much lower dust-to-gas ratio relative to the Galactic ISM, suggesting that the Broad Line Region (BLR) itself could provide the source of extra X-ray obscuration, being composed of low-ionization, dust-free gas., Comment: 13 pages, 9 figures. Accepted for publication in ApJ

We present a high density disc reflection spectral analysis of a sample of 17 Seyfert 1 galaxies to study the inner disc densities at different black hole mass scales and accretion rates. All the available XMM-Newton observations in the archive are used. OM observations in the optical/UV band are used to estimate their accretion rates. We find that 65% of sources in our sample show a disc density significantly higher than 1e15 cm^{-3}, which was assumed in previous reflection-based spectral analyses. The best-fit disc densities show an anti-correlation with black hole mass and mass accretion rate. High density disc reflection model can successfully explain the soft excess emission and significantly reduce inferred iron abundances. We also compare our black hole spin and disc inclination angle measurements with previous analyses., Comment: Accepted by MNRAS. The first paper of this series can be found here: arXiv:1901.01739

We report on the spectroscopic analysis of the black hole binary GX 339-4 during its recent 2017-2018 outburst, observed simultaneously by the Swift and NuSTAR observatories. Although during this particular outburst the source failed to make state transitions, and despite Sun constraints during the peak luminosity, we were able to trigger four different observations sampling the evolution of the source in the hard state. We show that even for the lowest luminosity observations the NuSTAR spectra show clear signatures of X-ray reprocessing (reflection) in an accretion disk. Detailed analysis of the highest signal-to-noise spectra with our family of relativistic reflection models RELXILL indicates the presence of both broad and narrow reflection components. We find that a dual-lamppost model provides a superior fit when compared to the standard single lamppost plus distant neutral reflection. In the dual lamppost model two sources at different heights are placed on the rotational axis of the black hole, suggesting that the narrow component of the Fe K emission is likely to originate in regions far away in the disk, but still significantly affected by its rotational motions. Regardless of the geometry assumed, we find that the inner edge of the accretion disk reaches a few gravitational radii in all our fits, consistent with previous determinations at similar luminosity levels. This confirms a very low degree of disk truncation for this source at luminosities above ~1% Eddington. Our estimates of Rin reinforces the suggested behavior for an inner disk that approaches the inner-most regions as the luminosity increases in the hard state., Comment: Accepted for publication in ApJ. 15 pages, 6 figures

We report results from joint NuSTAR, Swift and XMM-Newton observations of the newly discovered black hole X-ray binary candidate Swift J1658.2-4242 in the intermediate state. We observe a peculiar event in this source, with its X-ray flux rapidly decreasing by $\sim$45\% in $\sim$40~s, accompanied by only subtle changes in the shape of the broadband X-ray spectrum. In addition, we find a sudden turn-on of a transient QPO with a frequency of $6-7$~Hz around the time of the flux change, and the total fractional rms amplitude of the power spectrum increases from $\sim$2\% to $\sim$10\%. X-ray spectral and timing analyses indicate that the flux decrease is driven by intrinsic changes in the accretion flow around the black hole, rather than intervening material along the line of sight. In addition, we do not significantly detect any relativistic disk reflection component, indicating it is much weaker than previously observed while the source was in the bright hard state. We propose accretion disk instabilities triggered at a large disk radius as the origin of the fast transition in spectral and timing properties, and discuss possible causes of the unusual properties observed in Swift J1658.2-4242. The prompt flux variation detected along with the emergence of a QPO makes the event an interesting case for investigating QPO mechanisms in black hole X-ray binaries., Comment: accepted for publication in ApJ

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs., Comment: 15 pages, 7 figures, 5 tables. Revised version

Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $\Gamma=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $\Gamma=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $\Gamma=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object., Comment: 19 pages, 11 figures, 6 tables. Accepted for publication by ApJ

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., Comment: 40 pages, 11 figures, accepted for publication in ApJ

We investigated the time-variability and spectral properties of the eclipsing X-ray source Circinus Galaxy X-1 (GG X-1), using Chandra, XMM-Newton and ROSAT. We phase-connected the lightcurves observed over 20 years, and obtained a best-fitting period $P = (25,970.0 \pm 0.1)$ s $\approx$7.2 hr, and a period derivative $\dot{P}/P = ( 10.2\pm4.6) \times 10^{-7}$ yr$^{-1}$. The X-ray lightcurve shows asymmetric eclipses, with sharp ingresses and slow, irregular egresses. The eclipse profile and duration vary substantially from cycle to cycle. We show that the X-ray spectra are consistent with a power-law-like component, absorbed by neutral and ionized Compton-thin material, and by a Compton-thick, partial-covering medium, responsible for the irregular dips. The high X-ray/optical flux ratio rules out the possibility that CG X-1 is a foreground Cataclysmic Variable; in agreement with previous studies, we conclude that it is the first example of a compact ultraluminous X-ray source fed by a Wolf-Rayet star or stripped Helium star. Its unocculted luminosity varies between $\approx$4 $\times 10^{39}$ erg s$^{-1}$ and $\approx$3 $\times 10^{40}$ erg s$^{-1}$. Both the donor star and the super-Eddington compact object drive powerful outflows: we suggest that the occulting clouds are produced in the wind-wind collision region and in the bow shock in front of the compact object. Among the rare sample of Wolf-Rayet X-ray binaries, CG X-1 is an exceptional target for studies of super-critical accretion and close binary evolution; it is also a likely progenitor of gravitational wave events., Comment: 27 pages, 11 figures, accepted by ApJ

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

We present a broad band spectral analysis of the black hole binary GX~339-4 with NuSTAR and Swift using high density reflection model. The observations were taken when the source was in low flux hard states (LF) during the outbursts in 2013 and 2015, and in a very high flux soft state (HF) in 2015. The high density reflection model can explain its LF spectra with no requirement for an additional low temperature thermal component. This model enables us to constrain the density in the disc surface of GX~339-4 in different flux states. The disc density in the LF state is $\log(n_{\rm e}/$ cm$^{-3})\approx21$, 100 times higher than the density in the HF state ($\log(n_{\rm e}/$ cm$^{-3})=18.93^{+0.12}_{-0.16}$). A close-to-solar iron abundance is obtained by modelling the LF and HF broad band spectra with variable density reflection model ($Z_{\rm Fe}=1.50^{+0.12}_{-0.04}Z_{\odot}$ and $Z_{\rm Fe}=1.05^{+0.17}_{-0.15}Z_{\odot}$ respectively)., Comment: Accepted by MNRAS

We present a detailed analysis of the spectral properties of the Seyfert 1 galaxy 1H0419-577, based on the archival XMM-Newton, NuSTAR and simultaneous Swift observations taken between 2002-2015. All the observations show a broad emission line feature at the iron band. We demonstrate that the broad band spectral variability at different levels can be explained by the combination of light-bending effects in the vicinity of the central black hole plus a thin warm absorber. We obtain a black hole spin of a > 0.98 by fitting the multi-epoch spectra with the relativistic disc reflection model. 1H0419-577 is accreting at 40% of its Eddington limit and its X-ray band shows the hardest powerlaw continuum in the highest flux state, which was previously more commonly seen in AGNs with a low accretion rate (e.g. $L_{\rm X} /L_{\rm Edd} < 10^{-2}$). The NuSTAR observation shows a cool coronal temperature of $kT=30^{+22}_{-7}$keV in the high flux state., Comment: Accepted by MNRAS

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

Previous X-ray spectral analysis has revealed an increasing number of AGNs with high accretion rates where an outflow with a mildly relativistic velocity originates from the inner accretion disk. Here we report the detection of a new ultra-fast outflow (UFO) with a velocity of $v_{\rm out}=0.319^{+0.005}_{-0.008}c$ in addition to a relativistic disk reflection component in a poorly studied NLS1 WKK~4438, based on archival \nustar and \suzaku observations. The spectra of both \suzaku and \nustar observations show an Fe~\textsc{xxvi} absorption feature and the \suzaku data also show evidence for an Ar~\textsc{xviii} with the same blueshift. A super-solar argon abundance ($Z^{\prime}_{\rm Ar}>6Z_{\odot}$) and a slight iron over-abundance ($Z^{\prime}_{\rm Fe}=2.6^{+1.9}_{-2.0}Z_{\odot}$) are found in our spectral modelling. Based on Monte-Carlo simulations, the detection of the UFO is estimated to be around at 3$\sigma$ significance. The fast wind most likely arises from a radius of $\geq20r_g$ away from the central black hole. The disk is accreting at a high Eddington ratio ($L_{\rm bol}=0.4-0.7L_{\rm Edd}$). The mass outflow rate of the UFO is comparable with the disk mass inflow rate ($\dot M_{\rm out}>30\%\dot M_{\rm in}$), assuming a maximum covering factor. The kinetic power of the wind might not be high enough to have influence in AGN feedback ($\dot E_{\rm wind}/L_{\rm bol}\approx 3-5\%$) due to a relatively small column density ($12^{+9}_{-4}\times10^{22}$~cm$^{-2}$). However note that both the inferred velocity and the column density could be lower limits owing to the low viewing angle ($i=23^{+3}_{-2}$$^{\circ}$)., Comment: 7 pages, 3 figures, accepted by MNRAS

We present a spectral and timing analysis of the newly reported Galactic X-ray transient Swift J1658.2-4242 observed by NuSTAR and Swift. The broad-band X-ray continuum is typical of a black hole binary in the bright hard state, with a photon index of $\Gamma=1.63\pm0.02$ and a low coronal temperature of $kT_{\rm e}=22\pm1$ keV, corresponding to a low spectral cutoff well constrained by NuSTAR. Spectral modeling of the relativistic disk reflection features, consisting of a broad Fe K$\alpha$ line and the Compton reflection hump, reveals that the black hole is rapidly spinning with the spin parameter of $a^{*}>0.96$, and the inner accretion disk is viewed at a high inclination angle of $i=64^{+2}_{-3}{^\circ}$ (statistical errors, 90% confidence). The high inclination is independently confirmed by dips in the light curves, which can be explained by absorbing material located near the disk plane temporarily obscuring the central region. In addition, we detect an absorption line in the NuSTAR spectra centered at $7.03^{+0.04}_{-0.03}$ keV. If associated with ionized Fe K absorption lines, this provides evidence for the presence of outflowing material in the low/hard state of a black hole binary candidate. A timing analysis shows the presence of a type-C QPO in the power spectrum, with the frequency increasing from $\sim0.14$ Hz to $\sim0.21$ Hz during the single NuSTAR exposure. Our analysis reveals that Swift J1658.2-4242 displays characteristics typical for a black hole binary that is viewed at a high inclination angle, making it a good system for studying the accretion geometry in black hole binaries., Comment: 8 pages, 7 figures, accepted for publication in ApJ

Ultraluminous X-ray sources (ULXs) are bright X-ray sources in nearby galaxies not associated with the central supermassive black hole. Their luminosities imply they are powered by either an extreme accretion rate onto a compact stellar remnant, or an intermediate mass ($\sim100-10^5$ M$_{\odot}$) black hole. The recent detection of coherent pulsations coming from three bright ULXs demonstrates that some of these sources are powered by accretion onto a neutron star, implying accretion rates significantly in excess of the Eddington limit, a high degree of geometric beaming, or both. The physical challenges associated with the high implied accretion rates can be mitigated if the neutron star surface field is high - in the magnetar regime ($10^{14}$ G), since this suppresses the electron scattering cross section, reducing the radiation pressure that chokes off accretion for high luminosities. One of the few ways to determine surface magnetic fields is through the detection of cyclotron resonance scattering features (CRSFs) produced by the transition of charged particles between quantized Landau levels. To date, CRSFs have only been detected in Galactic accreting pulsars. Here we present the detection at 3.8-{\sigma} significance of a strong absorption line at a rest-frame energy of 4.5 keV in the Chandra spectrum of a ULX in M51. We find that this feature is likely to be a CRSF produced by the strong magnetic field of a neutron star. Assuming scattering off electrons, the magnetic field strength is implied to be $\sim10^{11}$ G, however the line is narrower than any electron CRSFs previously observed, and assuming thermal broadening, the implied temperature is significantly cooler than the accretion disk or column. The line shape is, however, consistent with a proton resonance scattering feature, implying that the neutron star has a magnetic field near the surface of B$\sim10^{15}$ G., Comment: Author version of the paper published in Nature Astronomy on 26 February 2018. See free-to-view published version here: http://rdcu.be/HQpR

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., Comment: 14 pages, accepted by ApJ