Fast Radio Bursts (FRBs) are millisecond-long radio pulses of extragalactic origin with peak luminosities far exceeding any Milky Way sources. The prevalent invocation for the FRB origin involves magnetars: young, magnetically powered neutron stars with the strongest magnetic fields in the Universe. A magnetar-defining signature is the emission of bright, hard X-ray bursts of sub-second duration. These occur in isolation or during a burst storm, when several hundred are observed within minutes to hours. On April 27th 2020, the Galactic magnetar SGR J1935+2154 entered an active period, emitting hundreds of X-ray bursts in a few hours. Remarkably, only one of these temporally coincided with an FRB-like radio burst. Here we report on the spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with NASA NICER and Fermi/GBM missions in their combined energy range. We demonstrate that the FRB-associated X-ray burst is very similar temporally, albeit strikingly different spectrally, from the 24 NICER/GBM bursts. If the FRB-associated burst were drawn from this magnetar burst population, its occurrence rate would be at most around 1 in 7000. This rarity combined with the unusual X-ray burst spectrum is perhaps indicative of an uncommon locale for the origin of the FRB-associated burst. We suggest that this unique event originated in quasi-polar open or closed magnetic field lines extending to high altitudes where radio emission can be generated, possibly from a collimated plasma flow.

We present a recent Chandra observation of the quiescent low-mass X-ray binary containing a neutron star, located in the globular cluster M30. We fit the thermal emission from the neutron star to extract its mass and radius. We find no evidence of flux variability between the two observations taken in 2001 and 2017, nor between individual 2017 observations, so we analyse them together to increase the signal to noise. We perform simultaneous spectral fits using standard light-element composition atmosphere models (hydrogen or helium), including absorption by the interstellar medium, correction for pile-up of X-ray photons on the detector, and a power-law for count excesses at high photon energy. Using a Markov-chain Monte Carlo approach, we extract mass and radius credible intervals for both chemical compositions of the atmosphere: $R_{\textrm{NS}}=7.94^{+0.76}_{-1.21}$ km and $M_{\textrm{NS}}<1.19$ M$_{\odot}$ assuming pure hydrogen, and $R_{\textrm{NS}}=10.50^{+2.88}_{-2.03}$ km and $M_{\textrm{NS}}<1.78$ M$_{\odot}$ for helium, where the uncertainties represent the 90% credible regions. For H, the small radius is difficult to reconcile with most current nuclear physics models (especially for nucleonic equations of state) and with other measurements of neutron star radii, with recent preferred values generally in the 11-14 km range. Whereas for He, the measured radius is consistent with this range. We discuss possible sources of systematic uncertainty that may result in an underestimation of the radius, identifying the presence of surface temperature inhomogeneities as the most relevant bias. According to this, we conclude that either the atmosphere is composed of He, or it is a H atmosphere with a significant contribution of hot spots to the observed radiation.

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

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

We aimed to identify the variable X-ray source 3XMM J000511.8+634018, which was serendipitously discovered through routine inspections while the 3XMM catalogue was compiled. We analysed the archival XMM-Newton observation of the source, obtained BUSCA photometry in three colours, and performed optical spectroscopy with the LBT. These data were supplemented by archival observations from the Zwicky Transient Facility. Based on its optical and X-ray properties, 3XMM J000511.8+634018 is classified as a magnetic cataclysmic variable, or polar. The flux is modulated with a period of 2.22 hr (8009.1 +- 0.2 s), which we identify with the orbital period. The bright phases are highly variable in X-ray luminosity from one cycle to the next. The source shows a thermal plasma spectrum typical of polars without evidence of a luminous soft blackbody-like component. It is non-eclipsing and displays one-pole accretion. The X-ray and BUSCA light curves show a stream absorption dip, which suggests an inclination 50degr < i < 75degr. The phasing of this feature, which occurs at the end of the bright phase, requires a somewhat special accretion geometry with a stream running far around the white dwarf before it is magnetically channelled. The period of this polar falls within the period gap of the cataclysmic variables (2.15-3.18 h), but appears to fall just below the minimum period when only polars are considered.

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution.

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.

We investigate the long-term variability of the K$\alpha$ line of iron in the spectra of two Ultra Compact X-ray Sources (UCXBs) with C/O-rich donors. We revisit archival observations by five different X-ray telescopes, over a ~twenty year period. Adopting physically motivated models for the spectral continuum, we probe the long-term evolution of the source emission in a self-consistent manner enabling physical interpretation of potential variability in the primary X-ray emission continuum and/or any emission lines from reflection off the accretion disk. We find that the spectral shape and flux of the source emission (for both objects) has remained almost constant throughout all the observations, displaying only minor variability in some spectral parameters and the source flux (largest variation is a ~25% drop in the flux of Swift J1756.9-2508). We note a striking variability of the Fe K$\alpha$ line which fluctuates from a notable equivalent width of ~66-100 eV in 4U 1543-624 and ~170 eV in Swift J1756.9-2508, to non-detections with upper limits of 2-8 eV. We argue that the disappearance of the iron line is due to the screening of the Fe K$\alpha$ line by the overabundant oxygen in the C/O-rich UCXBs. This effect is cancelled when oxygen becomes fully ionized in the inner disk region, resulting in the variability of the Fe K$\alpha$ line in an otherwise unaltered spectral shape. This finding supports earlier predictions on the consequences of H-poor, C/O-rich accretion disk on reflection induced fluorescent lines in the spectra of UCXBs.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

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

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.

Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

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.

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 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 analyzed the thermal emission from the entire surface of the millisecond pulsar PSR J0437-4715 observed in the ultraviolet and soft X-ray bands. For this, we calculated non-magnetized, partially ionized atmosphere models of hydrogen, helium, and iron compositions and included plasma frequency effects that may affect the emergent spectrum. This is particularly true for the coldest atmospheres composed of iron (up to a few percent changes in the soft X-ray flux). Employing a Markov chain Monte Carlo method, we found that the spectral fits favour a hydrogen atmosphere, disfavour a helium composition and rule out iron atmosphere and blackbody models. By using a Gaussian prior on the dust extinction, based on the latest 3D map of Galactic dust, and accounting for the presence of hot polar caps found in previous work, we found that the hydrogen atmosphere model results in a well-constrained neutron star radius $R_\mathrm{NS} = 13.6^{+0.9}_{-0.8}~\mathrm{km}$ and bulk surface temperature $T^{\infty}_\mathrm{eff}=\left(2.3\pm0.1\right)\times10^5~\mathrm{K}$. This relatively large radius favours a stiff equation of state and disfavours a strange quark composition inside neutron stars.

This paper presents a new analysis of the thermal emission from the neutron star surface to constrain the dense matter equation of state. It is based on the use of a Markov-Chain Monte Carlo algorithm combined with an empirical parametrization of the equation of state, as well as the consistent treatment of seven neutron star quiescent low-mass X-ray binaries in globular clusters with well-measured distances. Previous analyses have indicated that the thermal emission of these neutron stars tends to prefer low neutron star radii, questioning basic knowledge from nuclear physics. We show that it is possible to reconcile the thermal emission analyses with nuclear physics knowledge, with or without including a prior on the slope of the symmetry energy $L_{\rm sym}$. We obtain radii of the order of about 12~km without worsening the fit statistic. With an empirical parametrization of the equation of state, we obtain the following values for the slope of the symmetry energy, its curvature $K_{\rm sym}$, and the isoscalar skewness parameter $Q_{\rm sat}$: $L_{\rm sym}=37.2^{+9.2}_{-8.9}$ MeV, $K_{\rm sym}=-85^{+82}_{-70}$ MeV, and $Q_{\rm sat}=318^{+673}_{-366}$ MeV. For the first time, we measure the values of the empirical parameters $K_{\rm sym}$ and $Q_{\rm sat}$. These values are only weakly impacted by our assumptions, such as the distances or the number of free empirical parameters, provided they are taken within a reasonable range. We also study the weak sensitivity of our results to the set of sources analyzed, and we identify a group of sources that dominates the constraints. The resulting masses and radii obtained are also discussed in the context of the independent constraints from GW 170817 and its electromagnetic counterpart, AT 2017gfo.

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

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

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

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

Pulsating thermal X-ray emission from millisecond pulsars can be used to obtain constraints on the neutron star equation of state, but to date only five such sources have been identified. Of these five millisecond pulsars, only two have well constrained neutron star masses, which improve the determination of the radius via modelling of the X-ray waveform. We aim to find other millisecond pulsars that already have well constrained mass and distance measurements that show pulsed thermal X-ray emission in order to obtain tight constraints on the neutron star equation of state. The millisecond pulsar PSR~J1909--3744 has an accurately determined mass, M = 1.54$\pm$0.03 M$_\odot$ (1 $\sigma$ error) and distance, D = 1.07$\pm$0.04 kpc. We analysed {\em XMM-Newton} data of this 2.95 ms pulsar to identify the nature of the X-ray emission. We show that the X-ray emission from PSR~J1909--3744 appears to be dominated by thermal emission from the polar cap. Only a single component model is required to fit the data. The black-body temperature of this emission is kT=0.26\ud{0.03}{0.02} keV and we find a 0.2--10 keV un-absorbed flux of 1.1 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ or an un-absorbed luminosity of 1.5 $\times$ 10$^{30}$ erg s$^{-1}$. Thanks to the previously determined mass and distance constraints of the neutron star PSR~J1909--3744, and its predominantly thermal emission, deep observations of this object with future X-ray facilities should provide useful constraints on the neutron star equation of state.

Thermally emitting X-ray isolated neutron stars represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer physical parameters of the neutron star. Among the seven known members of this class, RX J1605.3+3249 is the only one that still lacks confirmation of its spin period. Here we analyze NICER and XMM-Newton observations of RX J1605.3+3249, in order to address its timing and spectral behavior. Contrary to a previous tentative detection, but in agreement with the recent work by Pires et al. (2019), we find no significant pulsation with pulsed fraction higher than 1.3% (3{\sigma}) for periods above 150 ms. We also find a limit of 2.6% for periods above 2 ms, despite searches in different energy bands. The X-ray spectrum can be fit by either a double-blackbody model or by a single-temperature magnetized atmosphere model, both modified by a Gaussian absorption line at ~0.44 keV. The origin of the absorption feature as a proton cyclotron line or as an atomic transition in the neutron star atmosphere is discussed. The predictions of the best-fit X-ray models extended to IR, optical and UV bands are compared with archival data. Our results are interpreted in the framework of a fallback disk scenario.

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

We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231$-$1411 using NICER. The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s$^{-1}$ in the 0.35-1.5 keV band. We analyzed Fermi LAT data to update the pulse timing model covering 10 years of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and $\gamma$ rays provides insight into the geometry of the system.

We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse.

We present the earliest available soft X-ray observations of XTE J1810-187, the prototypical transient magnetar, obtained 75--84 days after its 2018 outburst with the Neutron Star Interior Composition Explorer (NICER). Using a series of observations covering eight days we find that its decreasing X-ray flux is well-described by either a blackbody plus power-law or a two-blackbody spectral model. The 2-10 keV flux of the source varied from (1.206+/-0.007)x10^{-10} to (1.125+/-0.004)x10^{-10} erg s^{-1} cm^{-2}, a decrease of about 7% within our observations and 44% from that measured 7-14 days after the outburst with NuSTAR. We confirm that the pulsed fraction and spin pulse phase of the neutron star are energy dependent up to at least 8 keV. Phase resolved spectroscopy of the pulsar suggests magnetospheric variations relative to the line of sight.

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

Studying the physical processes occurring in the region just above the magnetic poles of strongly magnetized, accreting binary neutron stars is essential to our understanding of stellar and binary system evolution. Perhaps more importantly, it provides us with a natural laboratory for studying the physics of high temperature and high density plasmas exposed to extreme radiation, gravitational, and magnetic fields. Observations over the past decade have shed new light on the manner in which plasma falling at velocities near the speed of light onto a neutron star surface is halted. Recent advances in modeling these processes have resulted in direct measurement of the magnetic fields and plasma properties. On the other hand, numerous physical processes have been identified that challenge our current picture of how the accretion process onto neutron stars works. Observation and theory are our essential tools in this regime because the extreme conditions cannot be duplicated on Earth. This white paper gives an overview of the current theory, the outstanding theoretical and observational challenges, and the importance of addressing them in contemporary astrophysics research.

We report non-detections of the $\sim 3\times 10^8$ yr old, slow, isolated, rotation-powered pulsar PSR J2144$-$3933 in observations with the Hubble Space Telescope in one optical band (F475X) and two far-ultraviolet bands (F125LP and F140LP), yielding upper bounds $F_{\rm F475X}< 22.7$ nJy, $F_{\rm F125LP}< 5.9$ nJy, $F_{\rm F140LP}< 19.5$ nJy, at the pivot wavelengths 4940 \AA, 1438 \AA\ and 1528 \AA, respectively. Assuming a blackbody spectrum, we deduce a conservative upper bound on the surface (unredshifted) temperature of the pulsar of $T<42,000$ K. This makes PSR~J2144--3933 the coldest known neutron star, allowing us to study thermal evolution models of old neutron stars. This temperature is consistent with models with either direct or modified Urca reactions including rotochemical heating, and, considering frictional heating from the motion of neutron vortex lines, it puts an upper bound on the excess angular momentum in the neutron superfluid, $J<10^{44}\,\mathrm{erg\,s}$.

The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.

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

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a power law. Our refined timing measurements enable deeper searches for pulsations at other wavelengths and sensitive targeted searches by LIGO/Virgo for continuous gravitational waves from these neutron stars. Using the sensitivity of LIGO's first observing run, we estimate constraints that a gravitational wave search of these pulsars would obtain on the size of their mass deformation and r-mode fluid oscillation.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

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

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

The Neutron Star Interior Composition Explorer (NICER) on the International Space Station (ISS) observed strong photospheric expansion of the neutron star in 4U 1820-30 during a Type I X-ray burst. A thermonuclear helium flash in the star's envelope powered a burst that reached the Eddington limit. Radiation pressure pushed the photosphere out to ~200 km, while the blackbody temperature dropped to 0.45 keV. Previous observations of similar bursts were performed with instruments that are sensitive only above 3 keV, and the burst signal was weak at low temperatures. NICER's 0.2-12 keV passband enables the first complete detailed observation of strong expansion bursts. The strong expansion lasted only 0.6 s, and was followed by moderate expansion with a 20 km apparent radius, before the photosphere finally settled back down at 3 s after the burst onset. In addition to thermal emission from the neutron star, the NICER spectra reveal a second component that is well fit by optically thick Comptonization. During the strong expansion, this component is six times brighter than prior to the burst, and it accounts for 71% of the flux. In the moderate expansion phase, the Comptonization flux drops, while the thermal component brightens, and the total flux remains constant at the Eddington limit. We speculate that the thermal emission is reprocessed in the accretion environment to form the Comptonization component, and that changes in the covering fraction of the star explain the evolution of the relative contributions to the total flux.

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

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

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

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

Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G.

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 report on a Hubble Space Telescope detection of the nearby, old pulsar B0950+08 ($d\simeq 262$ pc, spin-down age 17.5 Myr) in two far-ultraviolet (FUV) bands. We measured the mean flux densities $\bar{f}_\nu = 109\pm 6$ nJy and $83\pm 14$ nJy in the F125LP and F140LP filters (pivot wavelengths 1438 and 1528 \AA). Using the FUV data together with previously obtained optical-UV data, we conclude that the optical-FUV spectrum consists of two components -- a nonthermal (presumably magnetospheric) power-law spectrum ($f_\nu\propto \nu^\alpha$) with slope $\alpha\sim -1.2$ and a thermal spectrum emitted from the bulk of the neutron star surface with a temperature in the range of $(1-3)\times 10^5$ K, depending on interstellar extinction and neutron star radius. These temperatures are much higher than predicted by neutron star cooling models for such an old pulsar, which means that some heating mechanisms operate in neutron stars. A plausible mechanism responsible for the high temperature of PSR B0950+08 is the interaction of vortex lines of the faster rotating neutron superfluid with the slower rotating normal matter in the inner neutron star crust (vortex creep heating).

We observed a nearby millisecond pulsar J2124-3358 with the Hubble Space Telescope in broad far-UV (FUV) and optical filters. The pulsar is detected in both bands with fluxes F(1250-2000 A)= (2.5+/-0.3)x10^-16 erg/s/cm^2 and F(3800-6000 A)=(6.4+/-0.4)x10^-17 erg/s/cm^2, which correspond to luminosities of ~5.8x10^27 and 1.4x10^27 erg/s, for d=410 pc and E(B-V)=0.03. The optical-FUV spectrum can be described by a power-law model, f_nu~nu^alpha, with slope alpha=0.18-0.48 for a conservative range of color excess, E(B-V)=0.01-0.08. Since a spectral flux rising with frequency is unusual for pulsar magnetospheric emission in this frequency range, it is possible that the spectrum is predominantly magnetospheric (power law with alpha<0) in the optical while it is dominated by thermal emission from the neutron star surface in the FUV. For a neutron star radius of 12 km, the surface temperature would be between 0.5x10^5 and 2.1x10^5 K, for alpha ranging from -1 to 0, E(B-V)=0.01-0.08, and d=340-500 pc. In addition to the pulsar, the FUV images reveal extended emission spatially coincident with the known Halpha bow shock, making PSR J2124-3358 the second pulsar (after PSR J0437-4715) with a bow shock detected in FUV.

We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $\Gamma= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase-resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsar's phase-folded lightcurve with the pulsar's well-defined mass and distance from radio timing observations.

We report on simultaneous NuSTAR and XMM-Newton observations of the magnetar 1E 1048.1$-$5937 obtained in July 2013, along with Rossi X-ray Timing Explorer (RXTE) data for the same source obtained from December 2002 through March 2012. The NuSTAR data provide a clear detection of this magnetar's persistent emission up to 20 keV. We detect a previously unreported small secondary peak in the average pulse profile in the 7-10 keV band, which grows to an amplitude comparable to that of the main peak in the 10--20 keV band. We show using RXTE data that this secondary peak is likely transient. We find that the pulsed fraction increases with energy from a value of $\sim$0.55 at $\sim$2 keV to a value of $\sim$0.75 near 8 keV but shows evidence for decreasing at higher energies. After filtering out multiple bright X-ray bursts during the observation, we find that the phase-averaged spectrum from combined NuSTAR and XMM data is well described by an absorbed double blackbody plus power-law model, with no evidence for the spectral turn-up near $\sim$10 keV as has been seen in some other magnetars. Our data allow us to rule out a spectral turn-up similar to those seen in magnetars 4U 0142+61 and 1E 2259+586 of $\Delta\Gamma >= 2$, where $\Delta\Gamma$ is the difference between the soft-band and hard-band photon indexes. The absence of a significant spectral turn-up is consistent with what has been observed from a particularly active subset of magnetars having high spin-inferred magnetic fields, and with previously reported trends suggesting the degree of spectral turn-up is correlated with spin-down rate and/or spin-inferred magnetic field strength.

We present a phase-coherent timing solution for PSR J1640-4631, a young 206 ms pulsar using X-ray timing observations taken with NuSTAR. Over this timing campaign, we have measured the braking index of PSR J1640-4631 to be n = 3.15+/-0.03. Using a series of simulations, we argue that this unusually high braking index is not due to timing noise, but is intrinsic to the pulsar's spin-down. We cannot, however, rule out contamination due to an unseen glitch recovery, although the recovery timescale would have to be longer than most yet observed. If this braking index is eventually proven to be stable, it demonstrates that pulsar braking indices greater than 3 are allowed in nature, hence other physical mechanisms such as mass or magnetic quadrupoles are important in pulsar spin-down. We also present a 3-sigma upper limit on the pulsed flux at 1.4 GHz of 0.018 mJy.

We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of ~2 Msun neutron stars, we place strong and quantitative constraints on the properties of the equation of state between approximately 2-8 times the nuclear saturation density. We find that around M=1.5 Msun, the preferred equation of state predicts radii between 10.1 - 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential.

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6x10^12 G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless, the modelling allows us to place constraints on the system geometry (i.e. the angles $\psi$ and $\xi$ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte-Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to $65\deg < \psi+\xi < 125\deg$ or $235\deg < \psi+\xi <295\deg$; (b) the temperature contrast between the poles and the equator must be at least a factor of $\sim6$, and (c) the size of the hottest region ranges between 0.2-0.7 km (including uncertainties on the source distance). Last, we interpret our findings within the context of internal and external heating models.

We present temporal and spectral analysis of simultaneous 0.5-79 keV Swift-XRT and NuSTAR observations of the magnetar 4U 0142+61. The pulse profile changes significantly with photon energy between 3 and 35 keV. The pulse fraction increases with energy, reaching a value of ~20%, similar to that observed in 1E 1841-045 and much lower than the ~80% pulse fraction observed in 1E 2259+586. We do not detect the 55-ks phase modulation reported in previous Suzaku-HXD observations. The phase-averaged spectrum of 4U 0142+61 above 20 keV is dominated by a hard power law with a photon index, $\Gamma$ ~ 0.65, and the spectrum below 20 keV can be described by two blackbodies, a blackbody plus a soft power law, or by a Comptonized blackbody model. We study the full phase-resolved spectra using the electron-positron outflow model of Beloborodov (2013). Our results are consistent with the parameters of the active j-bundle derived from INTEGRAL data by Hascoet et al. (2014). We find that a significant degeneracy appears in the inferred parameters if the footprint of the j-bundle is allowed to be a thin ring instead of a polar cap. The degeneracy is reduced when the footprint is required to be the hot spot inferred from the soft X-ray data.

Neutrons stars are unique laboratories to discriminate between the various proposed equations of state of matter at and above nuclear density. One sub-class of neutron stars - those inside quiescent low-mass X-ray binaries (qLMXBs) - produce a thermal surface emission from which the neutron star radius (R_NS) can be measured, using the widely accepted observational scenario for qLMXBs, assuming unmagnetized H atmospheres. In a combined spectral analysis, this work first reproduces a previously published measurement of the \rns, assumed to be the same for all neutron stars, using a slightly expanded data set. The radius measured is R_NS = 9.4 +/-1.2 km. On the basis of spectral analysis alone, this measured value is not affected by imposing an assumption of causality in the core. However, the assumptions underlying this R_NS measurement would be falsified by the observation of any neutron star with a mass >2.6 Msun, since radii <11 km would be rejected if causality is assumed, which would exclude most of the R_NS parameter space obtained in this analysis. Finally, this work directly tests a selection of dense matter equations of states: WFF1, AP4, MPA1, PAL1, MS0, and three versions of equations of state produced through chiral effective theory. Two of those, MSO and PAL1, are rejected at the 99% confidence level, accounting for all quantifiable uncertainties, while the other cannot be excluded at >99% certainty.

This paper presents the measurement of the neutron star (NS) radius using the thermal spectra from quiescent low-mass X-ray binaries (qLMXBs) inside globular clusters (GCs). Recent observations of NSs have presented evidence that cold ultra dense matter -- present in the core of NSs -- is best described by "normal matter" equations of state (EoSs). Such EoSs predict that the radii of NSs, Rns, are quasi-constant (within measurement errors, of ~10%) for astrophysically relevant masses (Mns > 0.5 Msun). The present work adopts this theoretical prediction as an assumption, and uses it to constrain a single Rns value from five qLMXB targets with available high signal-to-noise X-ray spectroscopic data. Employing a Markov-Chain Monte-Carlo approach, we produce the marginalized posterior distribution for Rns, constrained to be the same value for all five NSs in the sample. An effort was made to include all quantifiable sources of uncertainty into the uncertainty of the quoted radius measurement. These include the uncertainties in the distances to the GCs, the uncertainties due to the Galactic absorption in the direction of the GCs, and the possibility of a hard power-law spectral component for count excesses at high photon energy, which are observed in some qLMXBs in the Galactic plane. Using conservative assumptions,we found that the radius, common to the five qLMXBs and constant for a wide range of masses, lies in the low range of possible NS radii, Rns=9.1(+1.3)(-1.5) km (90%-confidence). Such a value is consistent with low-res equations of state. We compare this result with previous radius measurements of NSs from various analyses of different types of systems. In addition, we compare the spectral analyses of individual qLMXBs to previous works.

This paper reports the search for quiescent low-mass X-ray binaries (qLMXBs) in the globular cluster (GC) NGC 6553 using an XMM-Newton observation designed specifically for that purpose. We spectrally identify one candidate qLMXB in the core of the cluster, based on the consistency of the spectrum with a neutron star H-atmosphere model at the distance of NGC 6553. Specifically, the best-fit radius found using the three XMM European Photon Imaging Camera spectra is R_NS=6.3(+2.3)(-0.8) km (for M_NS=1.4 Msun) and the best-fit temperature is kTeff=136 (+21)(-34) eV. Both physical parameters are in accordance with typical values of previously identified qLMXBs in GC and in the field, i.e., R_NS~5-20 km and kTeff~50-150 eV. A power-law (PL) component with a photon index Gamma=2.1(+0.5)(-0.8) is also required for the spectral fit and contributes to ~33% of the total flux of the X-ray source. A detailed analysis supports the hypothesis that the PL component originates from nearby sources in the core, unresolved with XMM. The analysis of an archived Chandra observation provides marginal additional support to the stated hypothesis. Finally, a catalog of all the sources detected within the XMM field of view is presented here.

This paper reports the spectral and timing analyses of the quiescent low-mass X-ray binary U24 observed during five archived Chandra-ACIS exposures of the nearby globular cluster NGC 6397, for a total of 350 ksec. We find that the X-ray flux and the parameters of the hydrogen atmosphere spectral model are consistent with those previously published for this source. On short timescales, we find no evidence of aperiodic intensity variability, with 90% confidence upper limits during five observations ranging between <8.6% rms and <19% rms, in the 0.0001-0.1 Hz frequency range (0.5-8.0 keV); and no evidence of periodic variability, with maximum observed powers in this frequency range having a chance probability of occurrence from a Poisson-deviated light curve in excess of 10%. We also report the improved neutron star physical radius measurements, with statistical accuracy of the order of ~10%: R_ns = 8.9(+0.9)(-0.6) km for M_ns = 1.4 Msun. Alternatively, we provide the confidence regions in mass-radius space as well as the best-fit projected radius R_infinity= 11.9(+1.0)(-0.8)km, as seen by an observer at infinity. The best-fit effective temperature, kTeff = 80(+4)(-5) eV, is used to estimate the neutron star core temperature which falls in the range T_core = (3.0 - 9.8) x10 7 K, depending on the atmosphere model considered. This makes U24 the third most precisely measured neutron star radius among qLMXBs, after those in OmCen and in M13.

This paper presents the analysis of candidate quiescent low mass xray binarie (qLMXBs) observed during a short Chandra/ACIS observation of the globular cluster (GC) NGC 6304. Two out of the three candidate qLMXBs of this cluster, XMMU 171433-292747 and XMMU 171421-292917, lie within the field of view. This permits comparison with the discovery observation of these sources. The one in the GC core -- XMMU 171433-292747 -- is spatially resolved into two separate X-ray sources, one of which is consistent with a pure H-atmosphere qLMXB, and the other is an X-ray power-law spectrum source. These two spectral components separately account for those observed from XMMU 171433-292747 in its discovery observation. We find that the observed flux and spectral parameters of the H-atmosphere spectral components are consistent with the previous observation, as expected from a qLMXB powered by deep crustal heating. XMMU 171421-292917 also has neutron star atmosphere spectral parameters consistent with those in the XMM-Newton observation and the observed flux has decreased by a factor 0.54^{+0.30}_{-0.24}.

We report the search for low-mass X-ray binaries in quiescence (qLMXBs) in the globular cluster NGC 6304 using XMM observations. We present the spectral analysis leading to the identification of three candidate qLMXBs within the field of this globular cluster (GC), each consistent with the X-ray spectral properties of previously identified qLMXBs in the field and in other globular clusters -- specifically, with a hydrogen atmosphere neutron star with radius between 5--20\km. One (source 4, with R=11.7^{+8.3}_{-0.4} (D/5.97 kpc) km and kT_eff=117^{+59}_{-44} eV) is located within one core radius (r_c) of the centre of NGC 6304. This candidate also presents a spectral power-law component contributing 49 per cent of the 0.5-10 keV flux. A second one (source 9 with R=15.3^{+11.2}_{-6.5} (D/5.97 kpc) km and kT_eff=100^{+24}_{-19} eV) is found well outside the optical core (at 32 r_c) but still within the tidal radius. From spatial coincidence, we identify a bright 2MASS infrared counterpart which, at the distance of NGC 6304, seems to be a post-asymptotic giant branch star. The third qLMXB (source 5 with R=23^{+38}_{-14} (D/5.97 kpc) km and kT_eff=70^{+28}_{-20} eV) is a low signal-to-noise candidate for which we also identify from spatial coincidence a bright 2MASS infrared counterpart, with 99.916 per cent confidence. Three qLMXBs from this GC is marginally consistent with that expected from the encounter rate of NGC 6304. We also report a low signal-to-noise source with an unusually hard photon index (\alpha=-2.0^{+1.2}_{-2.2}). Finally, we present an updated catalogue of the X-ray sources lying in the field of NGC 6304, and compare this with the previous catalogue compiled from ROSAT observations.