The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems.

We report on the first simultaneous $NICER$ and $NuSTAR$ observations of the neutron star (NS) low-mass X-ray binary 4U 1735$-$44, obtained in 2018 August. The source was at a luminosity of $\sim1.8~(D/5.6\ \mathrm{kpc})^{2}\times10^{37}$ ergs s$^{-1}$ in the $0.4-30$ keV band. We account for the continuum emission with two different continuum descriptions that have been used to model the source previously. Despite the choice in continuum model, the combined passband reveals a broad Fe K line indicative of reflection in the spectrum. In order to account for the reflection spectrum we utilize a modified version of the reflection model RELXILL that is tailored for thermal emission from accreting NSs. Alternatively, we also use the reflection convolution model of RFXCONV to model the reflected emission that would arise from a Comptonized thermal component for comparison. We determine that the innermost region of the accretion disk extends close to the innermost stable circular orbit ($R_{\mathrm{ISCO}}$) at the 90% confidence level regardless of reflection model. Moreover, the current flux calibration of $NICER$ is within 5% of the $NuSTAR$/FPMA(B).

It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below $\sim 10^{-2}$ of the Eddington limit ($L_{\rm Edd}$). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608-52, to study changes in the reflection spectrum. We find that the broad Fe-K$\alpha$ line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is $\sim 10^{37}$ erg/s ($\sim 0.05$ $L_{\rm Edd}$), disappear during the decay of the outburst when the source luminosity drops to $\sim 4.5 \times 10^{35}$ erg/s ($\sim 0.002$ $L_{\rm Edd}$). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disk radius, disk ionisation, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disk ionisation ($\log \xi \geq 4.1$) or a much decreased reflection fraction. A changing disk truncation alone, however, cannot account for the lack of reprocessed Fe-K$\alpha$ emission. The required increase in ionisation parameter could occur if the inner accretion flow evaporates from a thin disk into a geometrically thicker flow, such as the commonly assumed formation of an radiatively inefficient accretion flow at lower mass accretion rates.

NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

The Neutron Star Interior Composition Explorer (NICER) has extensively monitored the August 2019 outburst of the 401 Hz millisecond X-ray pulsar SAX J1808.4-3658. In this Letter, we report on the detection of a bright helium-fueled Type I X-ray burst. With a bolometric peak flux of $(2.3\pm0.1)\times 10^{-7}$ erg/cm^2/s, this was the brightest X-ray burst among all bursting sources observed with NICER to date. The burst shows a remarkable two-stage evolution in flux, emission lines at $1.0$ keV and $6.7$ keV, and burst oscillations at the known pulsar spin frequency, with $\approx4$\% fractional sinusoidal amplitude. We interpret the burst flux evolution as the detection of the local Eddington limits associated with the hydrogen and helium layers of the neutron star envelope. The emission lines are likely associated with Fe, due to reprocessing of the burst emission in the accretion disk.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

We report for the first time below 1.5 keV, the detection of a secondary peak in an Eddington-limited thermonuclear X-ray burst observed by the Neutron Star Interior Composition Explorer (NICER) from the low-mass X-ray binary 4U 1608-52. Our time-resolved spectroscopy of the burst is consistent with a model consisting of a varying-temperature blackbody, and an evolving persistent flux contribution, likely attributed to the accretion process. The dip in the burst intensity before the secondary peak is also visible in the bolometric flux. Prior to the dip, the blackbody temperature reached a maximum of $\approx3$ keV. Our analysis suggests that the dip and secondary peak are not related to photospheric expansion, varying circumstellar absorption, or scattering. Instead, we discuss the observation in the context of hydrodynamical instabilities, thermonuclear flame spreading models, and re-burning in the cooling tail of the burst.

The Neutron Star Interior Composition Explorer (NICER) has observed seven thermonuclear X-ray bursts from the Low Mass X-ray Binary (LMXB) neutron star 4U 1728-34 from the start of the mission's operations until February of 2019. Three of these bursts show oscillations in their decaying tail with frequencies that are within 1 Hz of the previously detected burst oscillations from this source. Two of these burst oscillations have unusual properties: They have large fractional rms amplitudes of $ 48 \pm 9 \%$ and $ 46 \pm 9 \%$, and they are detected only at photon energies above 6 keV. By contrast, the third detected burst oscillation is compatible with previous observations of this source, with a fractional rms amplitude of $7.7 \pm 1.5\%$ rms in the 0.3 to 6.2 keV energy band. We discuss the implications of these large-amplitude burst oscillations, finding they are difficult to explain with the current theoretical models for X-ray burst tail oscillations.

We report on X-ray and radio observations of the ultra-compact X-ray binary 4U 1543-624 taken in August 2017 during an enhanced accretion episode. We obtained NICER monitoring of the source over a $\sim10$ day period during which target-of-opportunity observations were also conducted with Swift, INTEGRAL, and ATCA. Emission lines were measured in the NICER X-ray spectrum at $\sim0.64$ keV and $\sim6.4$ keV that correspond to O and Fe, respectively. By modeling these line components, we are able to track changes in the accretion disk throughout this period. The innermost accretion flow appears to move inwards from hundreds of gravitational radii ($R_{g}=GM/c^{2}$) at the beginning of the outburst to $<8.7$ $R_{g}$ at peak intensity. We do not detect the source in radio, but are able to place a $3\sigma$ upper limit on the flux density at $27$ $\mu$Jy beam$^{-1}$. Comparing the radio and X-ray luminosities, we find that the source lies significantly away from the range typical of black holes in the ${L}_{{r}}$-${L}_{{x}}$ plane, suggesting a neutron star (NS) primary. This adds to the evidence that NSs do not follow a single track in the ${L}_{{r}}$-${L}_{{x}}$ plane, limiting its use in distinguishing between different classes of NSs based on radio and X-ray observations alone.

We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of narrow emission and absorption lines during photospheric radius expansion (PRE) X-ray bursts from the ultracompact binary 4U 1820-30. NICER observed the source in 2017 August accumulating about 60 ks of exposure. Five thermonuclear X-ray bursts were detected of which four showed clear signs of PRE. We extracted spectra during the PRE phases and fit each to a model that includes a comptonized component to describe the accretion-driven emission, and a black body for the burst thermal radiation. The temperature and spherical emitting radius of the fitted black body are used to assess the strength of PRE in each burst. The two strongest PRE bursts (burst pair 1) had black body temperatures of approximately 0.6 keV and emitting radii of 100 km (at a distance of 8.4 kpc). The other two bursts (burst pair 2) had higher temperatures (~0.67 keV) and smaller radii (75 km). All of the PRE bursts show evidence of narrow line emission near 1 keV. By co-adding the PRE phase spectra of burst pairs 1 and, separately, 2 we find, in both co-added spectra, significant, narrow, spectral features near 1.0 (emission), 1.7 and 3.0 keV (both in absorption). Remarkably, all the fitted line centroids in the co-added spectrum of burst pair 1 appear systematically blue-shifted by a factor of $1.046 \pm 0.006$ compared to the centroids of pair 2, strongly indicative of a gravitational red-shift, a wind-induced blue-shift, or more likely some combination of both effects. The observed shifts are consistent with this scenario in that the stronger PRE bursts in pair 1 reach larger photospheric radii, and thus have weaker gravitational red-shifts, and they generate faster outflows, yielding higher blue-shifts. We discuss possible elemental identifications for the observed features in the context of recent burst-driven wind models.

We report on the Neutron Star Interior Composition Explorer (NICER) monitoring campaign of the 468 Hz accreting millisecond X-ray pulsar IGR J17379-3747. From a detailed spectral and timing analysis of the coherent pulsations we find that they show a strong energy dependence, with soft thermal emission lagging about 640 microseconds behind the hard, Comptonized emission. Additionally, we observe uncommonly large pulse fractions, with measured amplitudes in excess of 20% sinusoidal fractional amplitude across the NICER passband and fluctuations of up to ~70%. Based on a phase-resolved spectral analysis, we suggest that these extreme properties might be explained if the source has an unusually favorable viewing geometry with a large magnetic misalignment angle. Due to these large pulse fractions, we were able to detect pulsations down to quiescent luminosities (~5 x 10^33 erg s^-1). We discuss these low-luminosity pulsations in the context of transitional millisecond pulsars.

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

We present the discovery of a low-frequency $\approx 5.7$ Hz quasi-periodic oscillation (QPO) feature in observations of the black hole X-ray binary MAXI J1535-571 in its soft-intermediate state, obtained in September-October 2017 by the Neutron Star Interior Composition Explorer (NICER). The feature is relatively broad (compared to other low-frequency QPOs; quality factor $Q\approx 2$) and weak (1.9% rms in 3-10 keV), and is accompanied by a weak harmonic and low-amplitude broadband noise. These characteristics identify it as a weak Type A/B QPO, similar to ones previously identified in the soft-intermediate state of the transient black hole X-ray binary XTE J1550-564. The lag-energy spectrum of the QPO shows increasing soft lags towards lower energies, approaching 50 ms at 1 keV (with respect to a 3-10 keV continuum). This large phase shift has similar amplitude but opposite sign to that seen in Rossi X-ray Timing Explorer data for a Type B QPO from the transient black hole X-ray binary GX 339-4. Previous phase-resolved spectroscopy analysis of the Type B QPO in GX 339-4 pointed towards a precessing jet-like corona illuminating the accretion disk as the origin of the QPO signal. We suggest that this QPO in MAXI J1535-571 may have the same origin, with the different lag sign depending on the scale height of the emitting region and the observer inclination angle.

We report on the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum is dominated by Comptonisation of soft thermal seed photons with a temperature of ~0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.

Accretion disks around neutron stars regularly undergo sudden strong irradiation by Type I X-ray bursts powered by unstable thermonuclear burning on the stellar surface. We investigate the impact on the disk during one of the first X-ray burst observations with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. The burst is seen from Aql X-1 during the hard spectral state. In addition to thermal emission from the neutron star, the burst spectrum exhibits an excess of soft X-ray photons below 1 keV, where NICER's sensitivity peaks. We interpret the excess as a combination of reprocessing by the strongly photoionized disk and enhancement of the pre-burst persistent flux, possibly due to Poynting Robertson drag or coronal reprocessing. This is the first such detection for a short sub-Eddington burst. As these bursts are observed frequently, NICER will be able to study how X-ray bursts affect the disk and corona for a range of accreting neutron star systems and disk states.

We present results of recent Neutron Star Interior Composition Explorer observations of the accreting millisecond X-ray pulsar IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38 minute orbital period. NICER observed the source for approximately 26 ksec over a 5.3 day span in 2017 August, and again for 14 and 11 ksec in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hz in Rossi X-ray Timing Explorer data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z^2 method finds a best-fitting circular orbit with a period of 2278.21 s (37.97 min), a projected semi-major axis of 0.00390 lt-sec, and a barycentric pulsar frequency of 163.6561105 Hz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 e-8} solar masses, presently the smallest known for a stellar binary. The minimum donor mass ranges from about 0.005 - 0.007 solar masses, for a neutron star mass from 1.2 - 2 solar masses. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175 - 0.0155 solar masses and 19 deg < i < 27.5 deg, where the lower and upper bounds correspond to 1.4 and 2 solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 +- 0.11 % (0.3 - 3.2 keV).

We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of mHz X-ray brightness oscillations from the "clocked burster." NICER observed the source in the periods 2017 June 20 - 29, July 11 - 13, and September 9 - 15, for a total useful exposure of 34 ks. Two consecutive dwells obtained on 2017 September 9 revealed highly significant oscillations at a frequency of 8 mHz. The fractional, sinusoidal modulation amplitude increases from 0.7 % at 1 keV to approximately 2 % at 6 keV. Similar oscillations were also detected at lower significance in three additional dwells. The oscillation frequency and amplitude are consistent with those of mHz QPOs reported in other accreting neutron star systems. A thermonuclear X-ray burst was also observed on 2017 June 22. The burst properties and X-ray colors are both consistent with the source being in a soft spectral state during these observations, findings that are confirmed by ongoing monitoring with MAXI and SWIFT-BAT. Assuming that the mHz oscillations are associated with black body emission from the neutron star surface, modeling of the phase-resolved spectra shows that the oscillation is consistent with being produced by modulation of the temperature component of this emission. In this interpretation, the black body normalization, proportional to the emitting surface area, is consistent with being constant through the oscillation cycle. We place the observations in the context of the current theory of marginally stable burning and briefly discuss the potential for constraining neutron star properties using mHz oscillations.

We report on the detection of a kilohertz quasi-periodic oscillation (QPO) with the Neutron Star Interior Composition Explorer (NICER). Analyzing approximately 165 ks of NICER exposure on the X-ray burster 4U 0614+09, we detect multiple instances of a single-peak upper kHz QPO, with centroid frequencies that range from 400 Hz to 750 Hz. We resolve the kHz QPO as a function of energy, and measure, for the first time, the QPO amplitude below 2 keV. We find the fractional amplitude at 1 keV is on the order of 2% rms, and discuss the implications for the QPO emission process in the context of Comptonization models.

We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar Swift J1756.9$-$2508 during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combining our NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we also studied the long-term spin and orbital evolution of this source. We find that the binary system is well described by a constant orbital period model, with an upper limit on the orbital period derivative of $|\dot P_b| < 7.4\times 10^{-13}$ s s$^{-1}$. Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times, finding R.A. = $17^{\rm h}56^{\rm m}57.18^{\rm s}\pm0.08^{\rm s}$ and Decl. = $-25^{\circ}06'27.8''\pm3.5''$, while simultaneously measuring the long-term spin frequency derivative as $\dot\nu = -7.3\times 10^{-16}$ Hz s$^{-1}$. We briefly discuss the implications of these measurements in the context of the wider population of accreting millisecond pulsars.

We report on a NICER observation of the Galactic X-ray binary and stellar-mass black hole candidate, MAXI J1535-571. The source was likely observed in an "intermediate" or "very high" state, with important contributions from both an accretion disk and hard X-ray corona. The 2.3-10 keV spectrum shows clear hallmarks of relativistic disk reflection. Fits with a suitable model strongly indicate a near-maximal spin parameter of a = cJ/GM^2 = 0.994(2) and a disk that extends close to the innermost stable circular orbit, r/r_ISCO = 1.08(8) (1-sigma statistical errors). In addition to the relativistic spectrum from the innermost disk, a relatively narrow Fe K emission line is also required. The resolution of NICER reveals that the narrow line may be asymmetric, indicating a specific range of emission radii. Fits with a relativistic line model suggest an inner radius of r = 144 (+140,-60) GM/c^2 for the putative second reflection geometry; full reflection models suggest that radii a few times larger are possible. The origin of the narrow line is uncertain but a warp likely provides the most physically plausible explanation. We discuss our results in terms of the potential for NICER to reveal new features of the inner and intermediate accretion disk around black holes.

We report on a spectral-timing analysis of the neutron star low-mass X-ray binary Aql~X-1 with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. Aql~X-1 was observed with NICER during a dim outburst in 2017 July, collecting approximately $50$ ks of good exposure. The spectral and timing properties of the source correspond to that of an (hard) extreme island state in the atoll classification. We find that the fractional amplitude of the low frequency ($<0.3$ Hz) band-limited noise shows a dramatic turnover as a function of energy: it peaks at 0.5 keV with nearly 5\% rms, drops to $12\%$ rms at 2 keV, and rises to $15\%$ rms at 10 keV. Through the analysis of covariance spectra, we demonstrate that band-limited noise exists in both the soft thermal emission and the power-law emission. Additionally, we measure hard time lags, indicating the thermal emission at $0.5$ keV leads the power-law emission at 10 keV on a timescale of $\sim100$ ms at $0.3$ Hz to $\sim10$ ms at $3$ Hz. Our results demonstrate that the thermal emission in the hard state is intrinsically variable, and driving the modulation of the higher energy power-law. Interpreting the thermal spectrum as disk emission, we find our results are consistent with the disk propagation model proposed for accretion onto black holes.

We present Neutron Star Interior Composition Explorer (NICER) observations of the neutron star low-mass X-ray binary Serpens X-1 during the early mission phase in 2017. With the high spectral sensitivity and low-energy X-ray passband of NICER, we are able to detect the Fe L line complex in addition to the signature broad, asymmetric Fe K line. We confirm the presence of these lines by comparing the NICER data to archival observations with XMM-Newton/RGS and NuSTAR. Both features originate close to the innermost stable circular orbit (ISCO). When modeling the lines with the relativistic line model RELLINE, we find the Fe L blend requires an inner disk radius of $1.4_{-0.1}^{+0.2}$ $R_{\mathrm{ISCO}}$ and Fe K is at $1.03_{-0.03}^{+0.13}$ $R_{\mathrm{ISCO}}$ (errors quoted at 90%). This corresponds to a position of $17.3_{-1.2}^{+2.5}$ km and $12.7_{-0.4}^{+1.6}$ km for a canonical neutron star mass ($M_{\mathrm{NS}}=1.4\ M_{\odot}$) and dimensionless spin value of $a=0$. Additionally, we employ a new version of the RELXILL model tailored for neutron stars and determine that these features arise from a dense disk and supersolar Fe abundance.

We introduce a new method for analysing the aperiodic variability of coherent pulsations in accreting millisecond X-ray pulsars. Our method involves applying a complex frequency correction to the time-domain light curve, allowing for the aperiodic modulation of the pulse amplitude to be robustly extracted in the frequency domain. We discuss the statistical properties of the resulting modulation spectrum and show how it can be correlated with the non-pulsed emission to determine if the periodic and aperiodic variability are coupled processes. Using this method, we study the 598.88 Hz coherent pulsations of the accreting millisecond X-ray pulsar IGR J00291+5934 as observed with the Rossi X-ray Timing Explorer and XMM-Newton. We demonstrate that our method easily confirms the known coupling between the pulsations and a strong 8 mHz QPO in XMM-Newton observations. Applying our method to the RXTE observations, we further show, for the first time, that the much weaker 20 mHz QPO and its harmonic are also coupled the pulsations. We discuss the implications of this coupling and indicate how it may be used to extract new information on the underlying accretion process.

In this work I report on the stochastic X-ray variability of the 340 Hz accreting millisecond pulsar MAXI J0911-655. Analyzing pointed observations of the XMM-Newton and NuSTAR observatories I find that the source shows broad band-limited stochastic variability in the $0.01-10$ Hz range, with a total fractional variability of $\sim24\%$ rms in the $0.4-3$ keV energy band, which increases to $\sim40\%$ rms in the $3-10$ keV band. Additionally a pair of harmonically related quasi-periodic oscillations are discovered. The fundamental frequency of this harmonic pair is observed between frequencies of $62$ mHz and $146$ mHz. Like the band-limited noise, the amplitude of the QPOs show a steep increase as a function of energy, suggesting they share a similar origin, which is likely the inner accretion flow. Based on their energy dependence and their frequency relation with respect to the noise terms, the QPOs are identified as Low-Frequency oscillations, and discussed in terms of Lense-Thirring precession model.

The accreting millisecond X-ray pulsar (AMXP) SAX J1808.4-3658, shows a peculiar orbital evolution that proceeds at a much faster pace than predicted by conservative binary evolution models. It is important to identify the underlying mechanism responsible for this behavior because it can help to understand how this system evolves. It has also been suggested that, when in quiescence, SAX J1808.4-3658 turns on as a radio pulsar, a circumstance that might provide a link between AMXPs and black-widow radio pulsars. In this work we report the results of a deep radio pulsation search at 2 GHz using the Green Bank Telescope in August 2014 and an X-ray monitoring of the 2015 outburst with Chandra, Swift, and INTEGRAL. In particular, we present the X-ray timing analysis of a 30-ks Chandra observation executed during the 2015 outburst. We detect no radio pulsations, and place the strongest limit to date on the pulsed radio flux density of any AMXP. We also find that the orbit of SAX J1808.4-3658 continues evolving at a fast pace and we compare it to the bhevior of other accreting and non-accreting binaries. We discuss two scenarios: either the neutron star has a large moment of inertia (I>1.7x10^45 g cm^2) and is ablating the donor (by using its spin-down power) thus generating mass-loss with an efficiency of 40% or the donor star is undergoing quasi-cyclic variations due to a varying mass-quadrupole induced by either a strong (1 kG) field or by some unidentified mechanism probably linked to irradiation.

We study the 205.9 Hz pulsations of the accreting millisecond X-ray pulsar NGC 6440 X-2 across all outbursts observed with the Rossi X-ray Timing Explorer over a period of 800 days. We find the pulsations are highly sinusoidal with a fundamental amplitude of 5%-15% rms and a second harmonic that is only occasionally detected with amplitudes of <2% rms. By connecting the orbital phase across multiple outbursts, we obtain an accurate orbital ephemeris for this source and constrain its 57 min orbital period to sub-ms precision. We do not detect an orbital period derivative to an upper limit of $ | \dot{P} | \leq 8 \times 10^{-11}$ s/s. We investigate the possibility of coherently connecting the pulse phase across all observed outbursts, but find that due to the poorly constrained systematic uncertainties introduced by a flux-dependent bias in the pulse phase, multiple statistically acceptable phase-connected timing solutions exist.

In this work we have estimated upper and lower limits to the strength of the magnetic dipole moment of all 14 accreting millisecond X-ray pulsars observed with the Rossi X-ray Timing Explorer (RXTE). For each source we searched the archival RXTE data for the highest and lowest flux levels with a significant detection of pulsations. We assume these flux levels to correspond to the closest and farthest location of the inner edge of the accretion disk at which channelled accretion takes place. By estimating the accretion rate from the observed luminosity at these two flux levels, we place upper and lower limits on the magnetic dipole moment of the neutron star, using assumptions from standard magnetospheric accretion theory. Finally, we discuss how our field strength estimates can be further improved as more information on these pulsars is obtained.

We have studied the aperiodic variability of the 401 Hz accreting millisecond X-ray pulsar SAX J1808.4-3658 using the complete data set collected with the Rossi X-ray Timing Explorer over 14 years of observation. The source shows a number of exceptional aperiodic timing phenomena that are observed against a backdrop of timing properties that show consistent trends in all five observed outbursts and closely resemble those of other atoll sources. We performed a detailed study of the enigmatic ~410 Hz QPO, which has only been observed in SAX J1808.4-3658. We find that it appears only when the upper kHz QPO frequency is less than the 401 Hz spin frequency. The difference between the ~410 Hz QPO frequency and the spin frequency follows a similar frequency correlation as the low frequency power spectral components, suggesting that the ~410 Hz QPO is a retrograde beat against the spin frequency of a rotational phenomenon in the 9 Hz range. Comparing this 9 Hz beat feature with the Low-Frequency QPO in SAX J1808.4-3658 and other neutron star sources, we conclude that these two might be part of the same basic phenomenon. We suggest that they might be caused by retrograde precession due to a combination of relativistic, classical and magnetic torques. Additionally we present two new measurements of the lower kHz QPO, allowing us, for the first time, to measure the frequency evolution of the twin kHz QPOs in this source. The twin kHz QPOs are seen to move together over 150 Hz, maintaining a centroid frequency separation of $(0.446 \pm 0.009) \nu_{spin}$.

We study the relation between the 300-700 Hz upper kHz quasi-periodic oscillation (QPO) and the 401 Hz coherent pulsations across all outbursts of the accreting millisecond X-ray pulsar SAX J1808.4-3658 observed with the Rossi X-ray Timing Explorer. We find that the pulse amplitude systematically changes by a factor of ~2 when the upper kHz QPO frequency passes through 401 Hz: it halves when the QPO moves to above the spin frequency and doubles again on the way back. This establishes for the first time the existence of a direct effect of kHz QPOs on the millisecond pulsations and provides a new clue to the origin of the upper kHz QPO. We discuss several scenarios and conclude that while more complex explanations can not formally be excluded, our result strongly suggests that the QPO is produced by azimuthal motion at the inner edge of the accretion disk, most likely orbital motion. Depending on whether this azimuthal motion is faster or slower than the spin, the plasma then interacts differently with the neutron-star magnetic field. The most straightforward interpretation involves magnetospheric centrifugal inhibition of the accretion flow that sets in when the upper kHz QPO becomes slower than the spin.

We report the discovery of a 1--5~Hz X-ray flaring phenomenon observed at $>30$~mCrab near peak luminosity in the 2008 and 2011 outbursts of the accreting millisecond X-ray pulsar SAX J1808.4--3658 in observations with the \textit{Rossi X-ray Timing Explorer}. In each of the two outbursts this high luminosity flaring is seen for $\sim$3 continuous days and switches on and off on a timescale of 1--2~hr. The flaring can be seen directly in the light curve, where it shows sharp spikes of emission at quasi-regular separation. In the power spectrum it produces a broad noise component, which peaks at 1--5~Hz. The total 0.05--10~Hz variability has a fractional rms amplitude of 20\%--45\%, well in excess of the 8\%--12\% rms broad-band noise usually seen in power spectra of SAX J1808.4--3658. We perform a detailed timing analysis of the flaring and study its relation to the 401~Hz pulsations. We find that the pulse amplitude varies proportionally with source flux through all phases of the flaring, indicating that the flaring is likely due to mass density variations created at or outside the magnetospheric boundary. We suggest that this 1--5~Hz flaring is a high mass accretion rate version of the 0.5--2~Hz flaring which is known to occur at low luminosity ($<13$~mCrab), late in the tail of outbursts of SAX J1808.4--3658. We propose the dead-disk instability, previously suggested as the mechanism for the 0.5--2~Hz flaring, as a likely mechanism for the high luminosity flaring reported here.

The accreting millisecond pulsar SAX J1808.4-3658 has shown a peculiar orbital evolution in the past with an orbital expansion much faster than expected from standard binary evolutionary scenarios. Previous limits on the pulsar spin frequency derivative during transient accretion outbursts were smaller than predicted by standard magnetic accretion torque theory, while the spin evolution between outbursts was consistent with magnetic dipole spin-down. In this paper we present the results of a coherent timing analysis of the 2011 outburst observed by the Rossi X-ray Timing Explorer and extend our previous long-term measurements of the orbital and spin evolution over a baseline of thirteen years. We find that the expansion of the 2 hr orbit is accelerating at a rate 1.6E-20 s/s^2 and we interpret this as the effect of short-term angular momentum exchange between the mass donor and the orbit. The gravitational quadrupole coupling due to variations in the oblateness of the companion can be a viable mechanism for explaining the observations. No significant spin frequency derivatives are detected during the 2011 outburst (<4E-13 Hz/s) and the long term spin down remains stable over thirteen years with a rate of approximately -1E-15 Hz/s.