We analyze seven NICER and NuSTAR epochs of the black hole X-ray binary GX 339-4 in the hard state during its two most recent hard-only outbursts in 2017 and 2019. These observations cover the 1-100 keV unabsorbed luminosities between 0.3% and 2.1% of the Eddington limit. With NICER's negligible pile-up, high count rate and unprecedented time resolution, we perform a spectral-timing analysis and spectral modeling using relativistic and distant reflection models. Our spectral fitting shows that as the inner disk radius moves inwards, the thermal disk emission increases in flux and temperature, the disk becomes more highly ionized and the reflection fraction increases. This coincides with the inner disk increasing its radiative efficiency around ~1% Eddington. We see a hint of hysteresis effect at ~0.3% of Eddington: the inner radius is significantly truncated during the rise ($>49R_{g}$), while only a mild truncation ($\sim5R_g$) is found during the decay. At higher frequencies ($2-7$~Hz) in the highest luminosity epoch, a soft lag is present, whose energy dependence reveals a thermal reverberation lag, with an amplitude similar to previous findings for this source. We also discuss the plausibility of the hysteresis effect and the debate of the disk truncation problem in the hard state.

The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS.

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

The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.

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.

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3\sigma$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.

We report the quasi-simultaneous INTEGRAL, SWIFT, and NuSTAR observations showing spectral state transitions in the neutron star low mass X-ray binary 1RXS J180408.9-342058 during its 2015 outburst. We present results of the analysis of high-quality broad energy band (0.8-200 keV) data in three different spectral states: high/soft, low/very-hard, and transitional state. The broad band spectra can be described in general as the sum of thermal Comptonization and reflection due to illumination of an optically-thick accretion disc. During the high/soft state, blackbody emission is generated from the accretion disc and the surface of the neutron star. This emission, measured at a temperature of kTbb ~1.2 keV, is then Comptonized by a thick corona with an electron temperature of ~2.5 keV. For the transitional and low/very-hard state, the spectra are successfully explained with emission from a double Comptonizing corona. The first component is described by thermal Comptonization of seed disc/neutron-star photons (kTbb ~1.2 keV) by a cold corona cloud with kT e ~8-10 keV, while the second one originates from lower temperature blackbody photons (kTbb~0.1 keV) Comptonized by a hot corona (kTe~35 keV). Finally, from NuSTAR observations, there is evidence that the source is a new clocked burster. The average time between two successive X-ray bursts corresponds to ~7.9 ks and ~4.0 ks when the persistent emission decreases by a factor ~2, moving from very hard to transitional state. The accretion rate and the decay time of the X-ray bursts longer than ~30 s suggest that the thermonuclear emission is due to mixed H/He burning triggered by thermally unstable He ignition.

We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.

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.

We present an analysis of the time-averaged spectrum of the Seyfert-2 active galaxy NGC 4388, obtained by NICER. The intrinsic strength of the reflection spectrum in NGC 4388, the large collecting area and favorable pass band of NICER, and a net exposure of 105.6 ks yielded an exceptionally sensitive spectrum. Using two independent families of models, the intrinsic spectrum from the central engine is found to be highly obscured but not Compton-thick. Enforcing physical self-consistency within each model, the independent treatments give formally consistent results: N_H = 2.67 (-0.03,+0.02) E+23 cm^-2 or N_H = 2.64 (-0.03, +0.03) E+23 cm^-2. Past measurements made with Suzaku and XMM-Newton are in broad agreement with these column density values. A more recent measurement with NuSTAR (in late 2013) recorded a column density about twice as large; the robustness of this variability is reinforced by the use of consistent models and procedures. The neutral Fe K-alpha line in the NICER spectrum is nominally resolved and consistent with an origin in the optical broad line region (BLR). The data also require ionized absorption in the Fe K band, similar to the "warm absorbers" detected in Seyfert-1 active galactic nuclei (AGN). The low-energy spectrum is consistent with a set of ionized plasma components. We discuss these findings and note that the geometric inferences that derive from this analysis can be tested with XRISM and Athena.

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.

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.

Atoll sources are accreting neutron star (NS) low-mass X-ray binaries. We present a spectral analysis of four persistent atoll sources (GX 3+1, 4U 1702$-$429, 4U 0614+091, and 4U 1746$-$371) observed for $\sim20$ ks each with NuSTAR to determine the extent of the inner accretion disk. These sources range from an apparent luminosity of $0.006-0.11$ of the Eddington limit (assuming the empirical limit of $3.8\times10^{38}$ ergs s$^{-1}$). Broad Fe emission features shaped by Doppler and relativistic effects close to the NS were firmly detected in three of these sources. The position of the disk appears to be close to the innermost stable circular orbit (ISCO) in each case. For GX 3+1, we determine $R_{in}=1.8^{+0.2}_{-0.6}\ R_{\mathrm{ISCO}}$ (90% confidence level) and an inclination of $27^{\circ}-31^{\circ}$. For 4U 1702$-$429, we find a $R_{in}=1.5_{-0.4}^{+1.6}\ R_{\mathrm{ISCO}}$ and inclination of $53^{\circ}-64^{\circ}$. For 4U 0614+091, the disk has a position of $R_{in}=1.3_{-0.2}^{+5.4}\ R_{\mathrm{ISCO}}$ and inclination of $50^{\circ}-62^{\circ}$. If the disk does not extend to the innermost stable circular orbit, we can place conservative limits on the magnetic field strength in these systems in the event that the disk is truncated at the Alfv\'{e}n radius. This provides the limit at the poles of $B\leq6.7\times10^{8}$ G, $3.3\times10^{8}$ G, and $14.5\times10^{8}$ G for GX 3+1, 4U 1702$-$429, and 4U 0614+091, respectively. For 4U 1746$-$371, we argue that the most plausible explanation for the lack of reflection features is a combination of source geometry and strong Comptonization. We place these sources among the larger sample of NSs that have been observed with NuSTAR.

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.

Although the most luminous class of neutron star low mass X-ray binaries, known as Z sources, have been well studied, their behavior is not fully understood. In particular, what causes these sources to trace out the characteristic Z-shaped pattern on color-color or hardness-intensity diagrams is not well known. By studying the physical properties of the different spectral states of these sources, we may better understand such variability. With that goal in mind, we present a recent NuSTAR observation of the Z source GX 349+2, which spans approximately 2 days, and covers all its spectral states. By creating a hardness-intensity diagram we were able to extract five spectra and trace the change in spectral parameters throughout the Z-track. GX 349+2 shows a strong, broad Fe K$\alpha$ line in all states, regardless of the continuum model used. Through modeling of the reflection spectrum and Fe K$\alpha$ line we find that in most states the inner disk radius is consistent with remaining unchanged at an average radius of 17.5 $R_g$ or 36.4 km for a canonical 1.4 $M_\odot$ neutron star. During the brightest flaring branch, however, the inner disk radius from reflection is not well constrained.

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 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 present NuSTAR observations of the luminous neutron star low-mass X-ray binary (NS LMXB) and Z source GX 5-1. During our three observations made with separations of roughly two days, the source traced out an almost complete Z track. We extract spectra from the various branches and fit them with a continuum model that has been successfully applied to other Z sources. Surprisingly, and unlike most of the (luminous) NS-LMXBs observed with NuSTAR, we do not find evidence for reflection features in any of the spectra of GX 5-1. We discuss several possible explanations for the absence of reflection features. Based on a comparison with other accreting neutron star systems and given the high luminosity of GX 5-1 (~1.6-2.3 times the Eddington luminosity, for a distance of 9 kpc), we consider a highly ionized disk the most likely explanation for the absence of reflection features in GX 5-1.

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.

The bright, erratic black hole X-ray binary GRS 1915+105 has long been a target for studies of disk instabilities, radio/infrared jets, and accretion disk winds, with implications that often apply to sources that do not exhibit its exotic X-ray variability. With the launch of NICER, we have a new opportunity to study the disk wind in GRS 1915+105 and its variability on short and long timescales. Here we present our analysis of 39 NICER observations of GRS 1915+105 collected during five months of the mission data validation and verification phase, focusing on Fe XXV and Fe XXVI absorption. We report the detection of strong Fe XXVI in 32 (>80%) of these observations, with another four marginal detections; Fe XXV is less common, but both likely arise in the well-known disk wind. We explore how the properties of this wind depends on broad characteristics of the X-ray lightcurve: mean count rate, hardness ratio, and fractional RMS variability. The trends with count rate and RMS are consistent with an average wind column density that is fairly steady between observations but varies rapidly with the source on timescales of seconds. The line dependence on spectral hardness echoes known behavior of disk winds in outbursts of Galactic black holes; these results clearly indicate that NICER is a powerful tool for studying black hole winds.

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 perform a reflection study on a new observation of the neutron star low-mass X-ray binary Aquila X-1 taken with NuSTAR during the August 2016 outburst and compare with the July 2014 outburst. The source was captured at $\sim32\%\ L_{\mathrm{Edd}}$, which is over four times more luminous than the previous observation during the 2014 outburst. Both observations exhibit a broadened Fe line profile. Through reflection modeling, we determine that the inner disk is truncated $R_{in,\ 2016}=11_{-1}^{+2}\ R_{g}$ (where $R_{g}=GM/c^{2}$) and $R_{in,\ 2014}=14\pm2\ R_{g}$ (errors quoted at the 90% confidence level). Fiducial neutron star parameters (M$_{NS}=1.4$ M$_{\odot}$, $R_{NS}=10$ km) give a stellar radius of $R_{NS}=4.85\ R_{g}$; our measurements rule out a disk extending to that radius at more than the $6\sigma$ level of confidence. We are able to place an upper limit on the magnetic field strength of $B\leq3.0-4.5\times10^{9}$ G at the magnetic poles, assuming that the disk is truncated at the magnetospheric radius in each case. This is consistent with previous estimates of the magnetic field strength for Aquila X-1. However, if the magnetosphere is not responsible for truncating the disk prior to the neutron star surface, we estimate a boundary layer with a maximum extent of $R_{BL,\ 2016}\sim10\ R_{g}$ and $R_{BL,\ 2014}\sim6\ R_{g}$. Additionally, we compare the magnetic field strength inferred from the Fe line profile of Aquila X-1 and other neutron star low-mass X-ray binaries to known accreting millisecond X-ray pulsars.

We perform the first reflection study of the soft X-ray transient and Type 1 burst source XTE J1709-267 using NuSTAR observations during its 2016 June outburst. There was an increase in flux near the end of the observations, which corresponds to an increase from $\sim$0.04 L$_{\mathrm{Edd}}$ to $\sim$0.06 L$_{\mathrm{Edd}}$ assuming a distance of 8.5 kpc. We have separately examined spectra from the low and high flux intervals, which were soft and show evidence of a broad Fe K line. Fits to these intervals with relativistic disk reflection models have revealed an inner disk radius of $13.8_{-1.8}^{+3.0}\ R_{g}$ (where $R_{g} = GM/c^{2}$) for the low flux spectrum and $23.4_{-5.4}^{+15.6}\ R_{g}$ for the high flux spectrum at the 90\% confidence level. The disk is likely truncated by a boundary layer surrounding the neutron star or the magnetosphere. Based on the measured luminosity and using the accretion efficiency for a disk around a neutron star, we estimate that the theoretically expected size for the boundary layer would be $\sim0.9-1.1 \ R_{g}$ from the neutron star's surface, which can be increased by spin or viscosity effects. Another plausible scenario is that the disk could be truncated by the magnetosphere. We place a conservative upper limit on the strength of the magnetic field at the poles, assuming $a_{*}=0$ and $M_{NS}=1.4\ M_{\odot}$, of $B\leq0.75-3.70\times10^{9}$ G, though X-ray pulsations have not been detected from this source.

We present $\emph{NuSTAR}$ observations of neutron star (NS) low-mass X-ray binaries: 4U 1636-53, GX 17+2, and 4U 1705-44. We observed 4U 1636-53 in the hard state, with an Eddington fraction, $F_{\mathrm{Edd}}$, of 0.01; GX 17+2 and 4U 1705-44 were in the soft state with fractions of 0.57 and 0.10, respectively. Each spectrum shows evidence for a relativistically broadened Fe K$_{\alpha}$ line. Through accretion disk reflection modeling, we constrain the radius of the inner disk in 4U 1636-53 to be $R_{in}=1.03\pm0.03$ ISCO (innermost stable circular orbit) assuming a dimensionless spin parameter $a_{*}=cJ/GM^{2}=0.0$, and $R_{in}=1.08\pm0.06$ ISCO for $a_{*}=0.3$ (errors quoted at 1 $\sigma$). This value proves to be model independent. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, for example, $1.08\pm0.06$ ISCO translates to a physical radius of $R=10.8\pm0.6$ km, and the neutron star would have to be smaller than this radius (other outcomes are possible for allowed spin parameters and masses). For GX 17+2, $R_{in}=1.00-1.04$ ISCO for $a_{*}=0.0$ and $R_{in}=1.03-1.30$ ISCO for $a_{*}=0.3$. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, $R_{in}=1.03-1.30$ ISCO translates to $R=10.3-13.0$ km. The inner accretion disk in 4U 1705-44 may be truncated just above the stellar surface, perhaps by a boundary layer or magnetosphere; reflection models give a radius of 1.46-1.64 ISCO for $a_{*}=0.0$ and 1.69-1.93 ISCO for $a_{*}=0.3$. We discuss the implications that our results may have on the equation of state of ultradense, cold matter and our understanding of the innermost accretion flow onto neutron stars with low surface magnetic fields, and systematic errors related to the reflection models and spacetime metric around less idealized neutron stars.

1RXS J180408.9-342058 is a transient neutron star low-mass X-ray binary that exhibited a bright accretion outburst in 2015. We present Nustar, Swift, and Chandra observations obtained around the peak of this outburst. The source was in a soft X-ray spectral state and displayed an X-ray luminosity of Lx~(2-3)E37 (D/5.8 kpc)^2 erg cm-2 (0.5-10 keV). The Nustar data reveal a broad Fe-K emission line that we model as relativistically broadened reflection to constrain the accretion geometry. We found that the accretion disk is viewed at an inclination of i~27-35 degrees and extended close to the neutron star, down to Rin~5-7.5 gravitational radii (~11-17 km). This inner disk radius suggests that the neutron star magnetic field strength is B<2E8 G. We find a narrow absorption line in the Chandra/HEG data at an energy of ~7.64 keV with a significance of ~4.8 sigma. This feature could correspond to blue-shifted Fe xxvi and arise from an accretion disk wind, which would imply an outflow velocity of v~0.086c (~25800 km s-1). However, this would be extreme for an X-ray binary and it is unclear if a disk wind should be visible at the low inclination angle that we infer from our reflection analysis. Finally, we discuss how the X-ray and optical properties of 1RXS J180408.9-342058 are consistent with a relatively small (Porb<3 hr) binary orbit.

We report on observations of the neutron star (NS) residing in the low-mass X-ray binary 1RXS J180408.9-34205 taken 2015 March by $\emph{NuSTAR}$ and $\emph{XMM-Newton}$ while the source was in the hard spectral state. We find multiple reflection features (Fe K$_{\alpha}$ detected with $\emph{NuSTAR}$; N VII, O VII, and O VIII detected in the RGS) from different ionization zones. Through joint fits using the self consistent relativistic reflection model {\sc relxill}, we determine the inner radius to be $\leq 11.1\ R_{g}$. For a 1.4 M$_{\odot}$ NS with a spin of $a_{*}=0$, this is an inner disk radius of $\leq22.2$ km. We find the inclination of the system to be between $18^{\circ}$-$29^{\circ}$. If the disk is truncated at a radius greater than the neutron star radius, it could be truncated by a boundary layer on the neutron star surface. It is also possible that the disk is truncated at the magnetospheric radius; conservative estimates would then imply $B\leq(0.3 -1.0)\times10^{9}$ G at the magnetic poles, though coherent pulsations have not been detected and the source is not identified as a pulsar.

We systematically reanalyze two previous observations of the black hole (BH) GX 339-4 in the very high and intermediate state taken with $\emph{XMM-Newton}$ and $\emph{Suzaku}$. We utilize up-to-date data reduction procedures and implement the recently developed, self-consistent model for X-ray reflection and relativistic ray tracing, {\sc relxill}. In the very high and intermediate state, the rate of accretion is high and thus the disk remains close to the innermost stable circular orbit (ISCO). We require a common spin parameter and inclination when fitting the two observations since these parameters should remain constant across all states. This allows for the most accurate determination of the spin parameter of this galactic black hole binary from fitting the Fe K$\alpha$ emission line and provides a chance to test previous estimates. We find GX 339-4 to be consistent with a near maximally spinning black hole with a spin parameter $a_{*}$ $>0.97$ with an inclination of $36 \pm 4$ degrees. This spin value is consistent with previous high estimates for this object. Further, if the inner disk is aligned with the binary inclination, this modest inclination returns a high black hole mass, but they need not be aligned. Additionally, we explore how the spin is correlated with the power of the jet emitted but find no correlation between the two.

We study the X-ray properties of a sample of 14 optically-selected low-mass AGN whose masses lie within the range 1E5 -2E6 M(solar) with XMM-Newton. Only six of these low-mass AGN have previously been studied with sufficient quality X-ray data, thus, we more than double the number of low-mass AGN observed by XMM-Newton with the addition of our sample. We analyze their X-ray spectral properties and variability and compare the results to their more massive counterparts. The presence of a soft X-ray excess is detectable in all five objects which were not background dominated at 2-3 keV. Combined with previous studies, this gives a total of 8 low-mass AGN with a soft excess. The low-mass AGN exhibit rapid, short-term variability (hundreds to thousands of seconds) as well as long-term variability (months to years). There is a well-known anti-correlation between black hole mass and variability amplitude (normalized excess variance). Comparing our sample of low-mass AGN with this relation we find that all of our sample lie below an extrapolation of the linear relation. Such a flattening of the relation at low masses (below about 1E6 M(solar)) is expected if the variability in all AGN follows the same shape power spectrum with a break frequency that is dependent on mass. Finally, we also found two objects that show significant absorption in their X-ray spectrum, indicative of type 2 objects, although they are classified as type 1 AGN based on optical spectra.