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 −5 · 10 −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.

We present a broad-band X-ray study of the effect of superorbital periods on X-ray spectra and pulse profiles in the neutron star X-ray binaries LMC X-4 and SMC X-1. These two sources display periodic or quasi-periodic variations in luminosity on the order of tens of days which are known to be superorbital, and are attributed to warped, precessing accretion disks. Using joint NuSTAR and XMM-Newton observations that span a complete superorbital cycle, we examine the broad-band spectra of these sources and find the shape to be well described by an absorbed power law with a soft blackbody component. Changes in spectral shape and pulse profile shape are periodic with superorbital period, as expected from a precessing disk. We perform X-ray tomography using the changes in pulse profiles to model the geometry and kinematics of the inner accretion disk. Our simple beam and inner disk geometric model indicates that the long term changes in soft pulse shape and phase are consistent with reprocessed emission from a precessing inner disk.

We summarize the results of a dedicated effort made between 2012 and 2019 to follow the evolution of the cyclotron line in Her X-1 through repeated NuSTAR observations. The previously observed nearly 20-year-long decay of the cyclotron line energy has ended in 2012: from then on, the pulse-phase-averaged flux-corrected cyclotron line energy has remained stable and constant at an average value of Ecyc = (37.44 ± 0.07) keV (normalized to a flux level of 6.8 RXTE/ASM-cts s−1). The flux dependence of Ecyc discovered in 2007 is now measured with high precision, giving a slope of (0.675 ± 0.075) keV/(ASM-cts s−1), corresponding to an increase of 6.5% of Ecyc for an increase in flux by a factor of two. We also find that all line parameters as well as the continuum parameters show a correlation with X-ray flux. While a correlation between Ecyc and X-ray flux (both positive and negative) is now known for several accreting binaries with various suggestions for the underlying physics, the phenomenon of a long-term decay has so far only been seen in Her X-1 and Vela X-1, with far less convincing explanations.

Two NuSTAR observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse "turn-on" in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

We report on NuSTAR observations of transient pulsations in the neutron star X-ray binary SMC X-1. The transition from non-pulsing to pulsing states was not accompanied by a large change in flux. Instead, both pulsing and non-pulsing states were observed in a single observation during the low-flux super-orbital state. During the high-state, we measure a pulse period of P = 0.70117(9) s at T ref = 56145 MJD. Spectral analysis during non-pulsing and pulsing states reveals that the observations can be consistently modeled by an absorbed power law with a phenomenological cutoff resembling a Fermi-Dirac distribution, or by a partially obscured cutoff power law. The shapes of the underlying continua show little variability between epochs, while the covering fraction and column density vary between super-orbital states. The strength of pulsations also varies, leading us to infer that the absence and reemergence of pulsations are related to changing obscuration, such as by a warped accretion disk. SMC X-1 is accreting near or above its Eddington limit, reaching an unabsorbed X-ray luminosity of L X (2 − 10 keV) ≈ 5 × 10 38 erg s −1 . This suggests that SMC X-1 may be a useful local analog to ultraluminous X-ray pulsars (ULXPs), which likewise exhibit strong variability in their pulsed fractions, as well as flux variability on similar timescales. In particular, the gradual pulse turn-on which has been observed in M82 X-2 is similar to the behavior we observe in SMC X-1. Thus we propose that pulse fraction variability of ULXPs may also be due to variable obscuration.