Search

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., Comment: Accepted for publication in MNRAS

We present analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star low-mass X-ray binary GRS 1741.9-2853 during a period of outburst in May 2020. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe K$\alpha$ emission which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of $L_{\mathrm{bol}}=7.03^{+0.04}_{-0.05}\times10^{36}\,\mathrm{erg\,s^{-1}}$, assuming a distance of 7 kpc, corresponding to an Eddington ratio of $4.5\%$. This persistent luminosity combined with light curve analysis leads us to infer that the bursts were the result of pure He burning rather than mixed H/He burning. Time-resolved spectroscopy reveals that the bolometric flux of the first burst exhibits a double-peaked structure, placing the source within a small population of accreting neutron stars which exhibit multiple-peaked type-I X-ray bursts. We find that the second, brighter burst shows evidence for photospheric radius expansion (PRE) and that at its peak, this PRE event had an unabsorbed bolometric flux of $F_{\mathrm{peak}}=2.94^{+0.28}_{-0.26}\times10^{-8}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$. This yields a new distance estimate of $d=9.0\pm0.5$ kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical neutron star. Additionally, we performed a detailed timing analysis which failed to find evidence for quasiperiodic oscillations or burst oscillations, and we place an upper limit of $16\%$ on the rms variability around 589 Hz, the frequency at which oscillations have previously been reported., Comment: 15 pages, 4 figures, 4 tables, accepted for publication in ApJ

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740$+$6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument (NICER XTI) event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint NANOGrav and CHIME/Pulsar wideband radio timing measurements of arXiv:2104.00880. We use XMM European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740$+$6620 to be $12.39_{-0.98}^{+1.30}$ km and $2.072_{-0.066}^{+0.067}$ M$_{\odot}$ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully-ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is $\log_{10}(T$ [K]$)=5.99_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples., Comment: 40 pages, 16 figures (3 of which are figure sets), 1 animation, 2 tables. ApJL accepted version. Software: https://github.com/ThomasEdwardRiley/xpsi. Zenodo: https://doi.org/10.5281/zenodo.4697624. (v1 = submitted version; v2 = accepted version; v3 = metadata edits)

The hard X-ray spectrum of magnetic cataclysmic variables can be modelled to provide a measurement of white dwarf mass. This method is complementary to radial velocity measurements, which depend on the (typically rather uncertain) binary inclination. Here we present results from a Legacy Survey of 19 magnetic cataclysmic variables with NuSTAR. We fit accretion column models to their 20-78 keV spectra and derive the white dwarf masses, finding a weighted average $\bar{M}_{\rm WD}=0.77\pm0.02$ $M_{\odot}$, with a standard deviation $\sigma=0.10$ $M_{\odot}$, when we include the masses derived from previous NuSTAR observations of seven additional magnetic cataclysmic variables. We find that the mass distribution of accreting magnetic white dwarfs is consistent with that of white dwarfs in non-magnetic cataclysmic variables. Both peak at a higher mass than the distributions of isolated white dwarfs and post-common-envelope binaries. We speculate as to why this might be the case, proposing that consequential angular momentum losses may play a role in accreting magnetic white dwarfs and/or that our knowledge of how the white dwarf mass changes over accretion-nova cycles may also be incomplete., Comment: 14 pages. 5 figures in main paper, 1 figure in appendix. Accepted for publication in MNRAS

We present a detailed X-ray spectral and variability study of the full 2018 outburst of MAXI J1727-203 using NICER observations. The outburst lasted approximately four months. Spectral modelling in the 0.3-10 keV band shows the presence of both a soft thermal and a hard Comptonised component. The analysis of these components shows that MAXI J1727-203 evolved through the soft, intermediate and hard spectral states during the outburst. We find that the soft (disc) component was detected throughout almost the entire outburst, with temperatures ranging from ~0.4 keV, at the moment of maximum luminosity, to ~0.1 keV near the end of the outburst. The power spectrum in the hard and intermediate states shows broadband noise up to 20 Hz, with no evidence of quasi-periodic oscillations. We also study the rms spectra of the broadband noise at 0.3-10 keV of this source. We find that the fractional rms increases with energy in most of the outburst except during the hard state, where the fractional rms remains approximately constant with energy. We also find that, below 3 keV, the fractional rms follows the same trend generally observed at energies >3 keV, a behaviour known from previous studies of black holes and neutron stars. The spectral and timing evolution of MAXI J1727-203, as parametrised by the hardness-intensity, hardness-rms, and rms-intensity diagrams, suggest that the system hosts a black hole, although we could not rule out a neutron star., 15 pages, 9 figures, 1 table, Accepted for publication in MNRAS main journal

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., 18 pages, 8 figures. Accepted for publication in ApJ Letters

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)., Comment: Accepted for publication in ApJ, 8 pages, 1 table, 5 figures

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

We present the first joint NuSTAR and NICER observations of the ultra-compact X-ray binary (UCXB) 4U 1543$-$624 obtained in 2020 April. The source was at a luminosity of $L_{0.5-50\ \mathrm{keV}} = 4.9 (D/7\ \mathrm{kpc})^{2}\times10^{36}$ ergs s$^{-1}$ and showed evidence of reflected emission in the form of an O VIII line, Fe K line, and Compton hump within the spectrum. We used a full reflection model, known as xillverCO, that is tailored for the atypical abundances found in UCXBs, to account for the reflected emission. We tested the emission radii of the O and Fe line components and conclude that they originate from a common disk radius in the innermost region of the accretion disk ($R_{\rm in} \leq1.07\ R_{\mathrm{ISCO}}$). Assuming that the compact accretor is a neutron star (NS) and the position of the inner disk is the Alfv\'{e}n radius, we placed an upper limit on the magnetic field strength to be $B\leq0.7(D/7\ \mathrm {kpc})\times10^{8}$ G at the poles. Given the lack of pulsations detected and position of $R_{\rm in}$, it was likely that a boundary layer region had formed between the NS surface and inner edge of the accretion disk with an extent of 1.2 km. This implies a maximum radius of the neutron star accretor of $R_{\mathrm{NS}}\leq 12.1$ km when assuming a canonical NS mass of 1.4 $M_{\odot}$., Comment: Accepted for publication in ApJ, 10 pages, 6 figures, 1 table

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., Comment: Submitted to ApJ: June 1, 2017. Accepted for publication in ApJ: October 11, 2019

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 $, Comment: 13 pages, 9 figures, 2 tables, accepted for publication in ApJ

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., Accepted for publication in ApJ

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., Comment: 24 pages, 9 figures, accepted for publication in The Astrophysical Journal Letters

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., 28 pages, 9 figures; accepted for publication in the Astrophysical Journal Letters

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

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., Comment: 10 pages, 5 figures, submitted to ApJ Letters

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., Comment: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State; 76 pages, 24 figures, 7 tables, 8 figure sets (available in the online journal or from the authors)

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., 17 pages, 9 figures, accepted for publication in ApJ

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., Comment: 49 pages, 16 figures, part of The Astrophysical Journal Letters focus issue on the Neutron Star Interior Composition Explorer

We present StrayCats: a catalog of NuSTAR stray light observations of X-ray sources. Stray light observations arise for sources 1--4$^{\circ}$ away from the telescope pointing direction. At this off-axis angle, X-rays pass through a gap between optics and aperture stop and so do not interact with the X-ray optics but, instead, directly illuminate the NuSTAR focal plane. We have systematically identified and examined over 1400 potential observations resulting in a catalog of 436 telescope fields and 78 stray light sources that have been identified. The sources identified include historically known persistently bright X-ray sources, X-ray binaries in outburst, pulsars, and Type I X-ray bursters. In this paper we present an overview of the catalog and how we identified the StrayCats sources and the analysis techniques required to produce high level science products. Finally, we present a few brief examples of the science quality of these unique data., Comment: 17 pages, 14 figures, Accepted in ApJ

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., Accepted for publication in ApJ Letters

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., 6 pages, 5 figures, published by ApJ Letters

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., Accepted for publication in ApJL. Comments welcome

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., Comment: proof correction applied. 7 pages, 3 figures, 1 table, accepted for publication in ApJ Letters

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., Comment: 6 pages, 4 figures, 1 table, accepted for publication in ApJL

We study the X-ray properties of a sample of 14 optically-selected low-mass AGN whose masses lie within the range 10 5 2 × 10 6 M⊙ 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 � 10 6 M⊙) 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.

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., Comment: Accepted for publication in ApJ, 12 pages, 5 Figures, 6 Tables

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., Comment: Updated to match the published Journal version (2018, ApJ, 853, 157)

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., Comment: Accepted for publication in ApJ, 7 pages, 2 Tables, 5 Figures

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., 8 pages, 6 figures, Accepted for publication in ApJ

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, 10 pages, 3 tables, 5 figures, accepted for publication in MNRAS

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

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., Accepted for publication in ApJ Letters

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., Comment: 7 pages, 7 figures, 3 tables, Accepted to ApJ

We systematically reanalyze two previous observations of the black hole (BH) GX 339-4 in the very high and intermediate state taken with XMM-Newton and 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, 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. 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 BH binary from fitting the Fe Kα emission line and provides a chance to test previous estimates. We find GX 339-4 to be consistent with a near maximally spinning BH with a spin parameter with an inclination of 36° ± 4°. 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 BH 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.