Transient Neutron Star Research Articles

The Influence of the Magnetic Braking Laws on the Evolution of Persistent and Transient Low-mass X-Ray Binaries

Swift J1858.6−0814 (hereafter J1858) is a transient neutron star (NS) low-mass X-ray binary (LMXB). There is controversy regarding its donor mass derived from observations and theoretical calculations. In this paper, we adopt seven magnetic braking (MB) prescriptions suggested in the literature and different metallicity Z to simulate the evolution of the LMXB. Our results show that, employing the MB model proposed by A. Reiners & S. Mohanty (“rm12”), the convection- and rotation-boosted (“carb”) model, and the intermediate (“inter”) and convection-boosted (“cboost”) models in K. X. Van et al. can match (part of) the observational parameters of J1858 well. We then apply our method to other observed LMXBs and find that the “rm12” and “inter” MB laws are most promising in explaining transient LMXBs. In comparison, the simulations with the “cboost” and “carb” MB laws are more inclined to reproduce persistent LMXBs and ultracompact X-ray binaries, respectively. Our results, though subject to computational and/or observational bias, show that it is challenging to find a unified MB law that applies to the NS LMXB subpopulations simultaneously, indicating our lack of understanding of the true MB law. In addition, we explore the influence of various MB laws on the magnitude of the bifurcation periods in LMXBs.

A Study of Interstellar Medium in the Line of Sight of Transient Neutron Star Low-Mass X-ray Binary, MXB 1659-298, by Timing and Spectral Analysis

This work is dedicated to the study of interstellar medium (ISM) along the line of sight (LOS) of the transient low-mass X-ray binary, MXB 1659-298, capitalizing the high resolving power of XMM-Newton in the soft energy range. We emphasized the analysis of reflection grating spectrometer (RGS) data in the energy range 0.5–2.15 keV, suitable for the study of ISM. The paper includes an explanation of why, in the soft X-ray energy range, only two observations (out of seven) were deemed eligible for analysis. Three absorption lines associated with highly ionized Fe XX (1s22p2-2p2 (3p) 4d), Si XIV (1s2-1s2p), and Mg XI (1s2-1s6p) were identified in the observations, with IDs of 8620701(2001) and 748391601(2015). These new absorption lines and the absorption edge due to the neutral oxygen K edge seen in the spectra validate the multiphase structure of ISM. The predominance of interstellar medium over the ionized absorber is established along the direction of the source. The equivalent hydrogen column density measured is nearly equal to the galactic HI value derived previously. The small value of the ionic column density of Fe, Si, and Mg in the site of the high-temperature region resembles previous findings.

Probing outbursts of the transient neutron star low-mass X-ray binary Aql X-1 with NICER: a study of spectral evolution

ABSTRACT X-ray observations of neutron star (NS) low-mass X-ray binaries (LMXBs) are useful for probing the physical processes close to the NS and for constraining source parameters. Aql X-1 is a transient NS LMXB that frequently undergoes outbursts and provides an excellent opportunity to study source properties and accretion mechanisms in a strong-gravity regime over a wide range of accretion rates. In this work, we systematically investigate the spectral evolution of Aql X-1 usingNICER observations during the source outbursts in 2019 and 2020. The NICER observations cover the complete transition of the source from its canonical hard state to a soft state and back. The spectra extracted from most observations can be explained by a partially Comptonized accretion disc. We find that the system can be described by an accretion disc with an inner temperature of $\sim 0.7$ keV and a Comptonizing medium of thermal electrons at $\sim 2$ keV, while the photon index is strongly degenerate with the covering fraction of the medium. We also find evidence of Fe K$\alpha$ fluorescence emission in the spectra, indicating reprocessing of the Comptonized photons. We observe an absorption column density higher than the Galactic column density for most of the observations, indicating a significant local absorption. For some of the observations in the 2020 outburst, however, the local absorption is negligible.

Constraining the physical parameters of XTE J1701−462 through NuSTAR observations

ABSTRACT The spectral properties of the transient neutron star low-mass X-ray binary XTE J1701−462 were studied using the data obtained from FPMA/B detectors onboard NuSTAR during its second known outburst (2022 September). The physical parameters of the system were derived from the analysis of the data in the 3.0−30.0 keV energy range. The patterns displayed on the hardness–intensity diagram of the three observations closely resembled the banana branch/normal branch, a vertex of horizontal and normal branch of the Z-track and a transition from normal branch to flaring branch. Spectral analysis of the source revealed the presence of Fe K emission complex. The source spectra were fitted with a multitemperature blackbody ($\tt {diskbb}$) component in conjunction with the reflection model ($\tt {relxillNS}$). The values of temperature (kTin) and radius (Rin) of the inner accretion disc obtained from the spectral fitting with the model combination – $\tt {constant} \times \tt {tbabs} (\tt {diskbb}$ + $\tt {relxillNS})$ showed the source to be in its soft spectral state during the observations. The inclination angle (θ) of the source was estimated to be between 19° and 33° and the inner disc radius (Rin) was found to be 17.4 km. Assuming the case of magnetic truncation of accretion disc, the upper limits for the magnetic dipole moment (μ) and the magnetic field strength (B) at the poles of the neutron star in the system were found to be 5.78 × 1026 G cm3 and 8.23 × 108 G, respectively, for kA = 1.

X-ray iron absorption line in Swift J1858.6−0814

Abstract We present the spectral analysis of bright steady states in an outburst of the transient neutron star low-mass X-ray binary Swift J1858.6−0814 observed with the Neutron star Interior Composition Explorer (NICER). We detected an ionized iron K absorption line (H-like Fe) at 6.97 keV in the spectrum. We estimate the photoionization parameter using the ratio of the equivalent widths of the Fe xxvi (H-like) (17 ± 5 eV) and Fe xxv (He-like) (<3 eV) and discuss the origin of the iron absorption line. The irradiated gas producing the absorption line would locate within (3–6) × 109 cm of the X-ray source. We suggest that the observed H-like Fe absorption line originates from the highly ionized gas in the inner accretion disk in Swift J1858.6−0814.

Discovery of strongly variable X-ray polarization in the neutron star low-mass X-ray binary transient XTE J1701−462

Context. After about 16 years since its first outburst, the transient neutron star low-mass X-ray binary XTE J1701−462 turned on again in September 2022, allowing for the first study of its X-ray polarimetric characteristics by a dedicated observing program with the Imaging X-ray Polarimeter Explorer (IXPE). Aims. Polarimetric studies of XTE J1701−462 have been expected to improve our understanding of accreting weakly magnetized neutron stars, in particular, the physics and the geometry of the hot inner regions close to the compact object. Methods. The IXPE data of two triggered observations were analyzed using time-resolved spectroscopic and polarimetric techniques, following the source along its Z-track of the color–color diagram. Results. During the first pointing on 2022 September 29, an average 2–8 keV polarization degree of (4.6 ± 0.4)% was measured, the highest value found up to now for this class of sources. Conversely, only a ∼0.6% average degree was obtained during the second pointing ten days later. Conclusions. The polarimetric signal appears to be strictly related to the higher energy blackbody component associated with the boundary layer (BL) emission and its reflection from the inner accretion disk, and it is as strong as 6.1% and 1.2% (> 95% significant) above 3–4 keV for the two measurements, respectively. The variable polarimetric signal is apparently related to the spectral characteristics of XTE J1701−462, which is the strongest when the source was in the horizontal branch of its Z-track and the weakest in the normal branch. These IXPE results provide new important observational constraints on the physical models and geometry of the Z-sources. Here, we discuss the possible reasons for the presence of strong and variable polarization among these sources.

AstroSat and NuSTAR observations of XTE J1739−285 during the 2019–2020 outburst

ABSTRACT We report results from a study of XTE J1739−285, a transient neutron star low mass X-ray binary observed with AstroSat and NuSTAR during its 2019–2020 outburst. We detected accretion-powered X-ray pulsations at 386 Hz during very short intervals (0.5–1 s) of X-ray flares. These flares were observed during the 2019 observation of XTE J1739−285. During this observation, we also observed a correlation between intensity and hardness ratios, suggesting an increase in hardness with the increase in intensity. Moreover, a thermonuclear X-ray burst detected in our AstroSat observation during the 2020 outburst revealed the presence of coherent burst oscillations at 383 Hz during its decay phase. The frequency drift of 3 Hz during X-ray burst can be explained with r modes. Thus, making XTE J1739−285 belong to a subset of NS-LMXBs which exhibit both nuclear- and accretion-powered pulsations. The power density spectrum created using the AstroSat-laxpc observations in 2020 showed the presence of a quasi-periodic oscillation at ∼0.83 Hz. Our X-ray spectroscopy revealed significant changes in the spectra during the 2019 and 2020 outburst. We found a broad iron line emission feature in the X-ray spectrum during the 2020 observation, while this feature was relatively narrow and has a lower equivalent width in 2019, when the source was accreting at higher rates than 2020. Hard X-ray tail was observed during the 2019 observations, indicating the presence of non-thermal component in the X-ray spectra.

Photospheric Radius Expansion and a Double-peaked Type-I X-Ray Burst from GRS 1741.9–2853

Abstract We present an analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star (NS) low-mass X-ray binary GRS 1741.9–2853 during a period of outburst in 2020 May. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe Kα emission, which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of L bol = 7.03 − 0.05 + 0.04 × 10 36 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 NSs that 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 peak = 2.94 − 0.26 + 0.28 × 10 − 8 erg cm − 2 s − 1 . This yields a new distance estimate of d = 9.0 ± 0.5 kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical NS. Additionally, we performed a detailed timing analysis that failed to find evidence for quasi-periodic 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.

Constraining the properties of dense neutron star cores: the case of the low-mass X-ray binary HETE J1900.1–2455

ABSTRACT Measuring the time evolution of the effective surface temperature of neutron stars can provide invaluable information on the properties of their dense cores. Here, we report on a new Chandra observation of the transient neutron star low-mass X-ray binary HETE J1900.1–2455, which was obtained ≈2.5 yr after the end of its ≈10-yr long accretion outburst. The source is barely detected during the observation, collecting only six net photons, all below 2 keV. Assuming that the spectrum is shaped as a neutron star atmosphere model, we perform a statistical analysis to determine a 1σ confidence upper range for the neutron star temperature of ≈30–39 eV (for an observer at infinity), depending on its mass, radius, and distance. Given the heat injected into the neutron star during the accretion outburst, estimated from data provided by all-sky monitors, the inferred very low temperature suggests that the core either has a very high heat capacity or undergoes very rapid neutrino cooling. While the present data do not allow us to disentangle these two possibilities, both suggest that a significant fraction of the dense core is not superfluid/superconductor. Our modelling of the thermal evolution of the neutron star predicts that it may still cool further, down to a temperature of ≃15 eV. Measuring such a low temperature with a future observation may provide constraints on the fraction of baryons that is paired in the stellar core.

The variable radio counterpart of Swift J1858.6-0814

ABSTRACT Swift J1858.6-0814 is a transient neutron star X-ray binary discovered in 2018 October. Multiwavelength follow-up observations across the electromagnetic spectrum revealed many interesting properties, such as erratic flaring on minute time-scales and evidence for wind outflows at both X-ray and optical wavelengths, strong and variable local absorption, and an anomalously hard X-ray spectrum. Here, we report on a detailed radio observing campaign consisting of one observation at 5.5/9 GHz with the Australia Telescope Compact Array, and nine observations at 4.5/7.5 GHz with the Karl G. Jansky Very Large Array. A radio counterpart with a flat to inverted radio spectrum is detected in all observations, consistent with a compact jet being launched from the system. Swift J1858.6-0814 is highly variable at radio wavelengths in most observations, showing significant variability when imaged on 3-to-5-min time-scales and changing up to factors of 8 within 20 min. The periods of brightest radio emission are not associated with steep radio spectra, implying they do not originate from the launching of discrete ejecta. We find that the radio variability is similarly unlikely to have a geometric origin, be due to scintillation, or be causally related to the observed X-ray flaring. Instead, we find that it is consistent with being driven by variations in the accretion flow propagating down the compact jet. We compare the radio properties of SwiftJ1858.6-0814 with those of Eddington-limited X-ray binaries with similar X-ray and optical characteristics, but fail to find a match in radio variability, spectrum, and luminosity.

Thermal state of transiently accreting neutron stars with additional heating beyond deep crustal heating

Abstract As some neutron star transients require an additional unknown heat source (referred to as “shallow heating”) to explain their high temperatures at the beginning of quiescence, we investigate the effect of shallow heating as well as compressional heating on the thermal state of transiently accreting neutron stars with the use of evolutionary calculations in the present work. Through comparing our theoretical predictions of the equilibrium redshifted luminosities $(L_{\gamma}^{\infty})$ produced by both deep crustal heating and shallow heating/compressional heating for different time-averaged mass-accretion rates $\langle\dot{M}\rangle$ with 35 updated observations of soft X-ray transients, the results show that both shallow heating and compressional heating make significant contributions to the equilibrium redshifted luminosity. The hotter sources (XTE J1701, MAXI J0556, EXO 0748, Aql X-1 etc.) with higher accretion rates are more likely to be explained with the effect of shallow heating or compressional heating. In addition, for a proper shallow heat $q_\mathrm{sh}$ and mass-accretion rate $\dot{M}$, the effect of shallow heating could be simulated by compressional heating.

Long-term temperature evolution of neutron stars undergoing episodic accretion outbursts

Context. Transient neutron star low-mass X-ray binaries undergo episodes of accretion, alternated with quiescent periods. During an accretion outburst, the neutron star heats up due to exothermic accretion-induced processes taking place in the crust. Besides the long-known deep crustal heating of nuclear origin, a likely non-nuclear source of heat, dubbed “shallow heating”, is present at lower densities. Most of the accretion-induced heat slowly diffuses into the core on a timescale of years. Over many outburst cycles, a state of equilibrium is reached when the core temperature is high enough that the heating and cooling (photon and neutrino emission) processes are in balance.Aims. We investigate how stellar characteristics and outburst properties affect the long-term temperature evolution of a transiently accreting neutron star. For the first time the effects of crustal properties are considered, particularly that of shallow heating.Methods. Using our codeNSCool, we tracked the thermal evolution of a neutron star undergoing outbursts over a period of 105yr. The outburst sequence is based on the regular outbursts observed from the neutron star transient Aql X-1. For each model we calculated the timescale over which equilibrium was reached and we present these timescales along with the temperature and luminosity parameters of the equilibrium state.Results. We performed several simulations with scaled outburst accretion rates, to vary the amount of heating over the outburst cycles. The results of these models show that the equilibrium core temperature follows a logarithmic decay function with the equilibrium timescale. Secondly, we find that shallow heating significantly contributes to the equilibrium state. Increasing its strength raises the equilibrium core temperature. We find that if deep crustal heating is replaced by shallow heating alone, the core would still heat up, reaching only a 2% lower equilibrium core temperature. Deep crustal heating may therefore not be vital to the heating of the core. Additionally, shallow heating can increase the quiescent luminosity to values higher than previously expected. The thermal conductivity in the envelope and crust, including the potentially low-conductivity pasta layer at the bottom of the crust, is unable to significantly alter the long-term internal temperature evolution. Stellar compactness and nucleon pairing in the core change the specific heat and the total neutrino emission rate as a function of temperature, with the consequences for the properties of the equilibrium state depending on the exact details of the assumed pairing models. The presence of direct Urca emission leads to the lowest equilibrium core temperature and the shortest equilibrium timescale.

Deep crustal heating by neutrinos from the surface of accreting neutron stars

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.

Searching for γ-ray counterparts to very faint X-ray transient neutron star binaries

Very faint X-ray transients (VFXTs) are a group of X-ray binaries with low luminosities, displaying peak X-ray luminosities during their outbursts of only 1034–1036 erg s−1. Using γ-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope (Fermi), we investigate their possible nature of containing rotation-powered pulsars, or more specifically being transitional millisecond pulsars (MSPs). Among more than 40 known VFXTs, we select 12 neutron star systems. We analyze the LAT data for the fields of 12 VFXTs in the energy range 0.2–300 GeV, but do not find any counterparts likely detected by Fermi. We obtain luminosity upper limits for the 12 sources. While the distances to the sources are largely uncertain, the upper limits are comparable to the luminosities of two transitional systems, PSR J1023−0038 and XSS J12270−4859. From our study, we conclude that no evidence is found at γ-rays for the suggestion that some VFXTs could contain rotation-powered MSPs (or be transitional MSP systems).

Extreme quiescent variability of the transient neutron star low-mass X-ray binary EXO 1745−248 in Terzan 5

EXO 1745-248 is a transient neutron-star low-mass X-ray binary that resides in the globular cluster Terzan 5. We studied the transient during its quiescent state using 18 Chandra observations of the cluster acquired between 2003 and 2016. We found an extremely variable source, with a luminosity variation in the 0.5-10 keV energy range of $\sim3$ orders of magnitude (between $3\times10^{31}$ erg s$^{-1}$ and $2\times10^{34}$ erg s$^{-1}$) on timescales from years down to only a few days. Using an absorbed power-law model to fit its quiescent spectra, we obtained a typical photon index of $\sim1.4$, indicating that the source is even harder than during outburst and much harder than typical quiescent neutron stars if their quiescent X-ray spectra are also described by a single power-law model. This indicates that EXO 1745-248 is very hard throughout the entire observed X-ray luminosity range. At the highest luminosity, the spectrum fits better when an additional (soft) component is added to the model. All these quiescent properties are likely related to strong variability in the low-level accretion rate in the system. However, its extreme variable behavior is strikingly different from the one observed for other neutron star transients that are thought to still accrete in quiescence. We compare our results to these systems. We also discuss similarities and differences between our target and the transitional millisecond pulsar IGR J18245-2452, which also has hard spectra and strong variability during quiescence.

Different Accretion Heating of the Neutron Star Crust during Multiple Outbursts in MAXI J0556–332

Abstract The transient neutron star (NS) low-mass X-ray binary MAXI J0556−332 provides a rare opportunity to study NS crust heating and subsequent cooling for multiple outbursts of the same source. We examine MAXI, Swift, Chandra, and XMM-Newton data of MAXI J0556−332 obtained during and after three accretion outbursts of different durations and brightnesses. We report on new data obtained after outburst III. The source has been tracked up to ∼1800 days after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from ∼333 eV to ∼146 eV was observed during the first ∼1200 days. Outburst III reheated the crust up to ∼167 eV, after which the crust cooled again to ∼131 eV in ∼350 days. We model the thermal evolution of the crust and find that this source required a different strength and depth of shallow heating during each of the three outbursts. The shallow heating released during outburst I was ∼17 MeV nucleon−1 and outburst III required ∼0.3 MeV nucleon−1. These cooling observations could not be explained without shallow heating. The shallow heating for outburst II was not well constrained and could vary from ∼0 to 2.2 MeV nucleon−1, i.e., this outburst could in principle be explained without invoking shallow heating. We discuss the nature of the shallow heating and why it may occur at different strengths and depths during different outbursts.

Discovery of the New X-Ray Transient MAXI J1807+132: A Candidate of a Neutron Star Low-mass X-Ray Binary

Abstract We report on the detection and follow-up multi-wavelength observations of the new X-ray transient MAXI J1807+132 with the MAXI/GSC, Swift, and ground-based optical telescopes. The source was first recognized with the MAXI/GSC on 2017 March 13. About a week later, it reached maximum intensity (∼10 mCrab in 2–10 keV), and then gradually faded in ∼10 days by more than one order of magnitude. Time-averaged Swift/X-ray Telescope spectra in the decaying phase can be described by a blackbody with a relatively low temperature (0.1–0.5 keV), plus a hard power-law component with a photon index of ∼2. These spectral properties are similar to those of neutron star low-mass X-ray binaries (LMXBs) in their dim periods. The blackbody temperature and the radius of the emission region varied in a complex manner as the source became dimmer. The source was detected in the optical wavelength on March 27–31 as well. The optical flux decreased monotonically as the X-ray flux decayed. The correlation between the X-ray and optical fluxes is found to be consistent with those of known neutron star LMXBs, supporting the idea that the source is likely to be a transient neutron star LMXB.

NuSTAR + XMM-Newton monitoring of the neutron star transient AX J1745.6-2901

AX J1745.6-2901 is a high-inclination (eclipsing) transient neutron star (NS) Low Mass X-ray Binary (LMXB) showcasing intense ionised Fe K absorption. We present here the analysis of 11 XMM-Newton and 15 NuSTAR new data-sets (obtained between 2013-2016), therefore tripling the number of observations of AX J1745.6-2901 in outburst. Thanks to simultaneous XMM-Newton and NuSTAR spectra, we greatly improve on the fitting of the X-ray continuum. During the soft state the emission can be described by a disk black body ($kT\sim1.1-1.2$ keV and inner disc radius $r_{DBB}\sim14$ km), plus hot ($kT\sim2.2-3.0$ keV) black body radiation with a small emitting radius ($r_{BB}\sim0.5-0.8$ km) likely associated with the boundary layer or NS surface, plus a faint Comptonisation component. Imprinted on the spectra are clear absorption features created by both neutral and ionised matter. Additionally, positive residuals suggestive of an emission Fe K$\alpha$ disc line and consistent with relativistic ionised reflection are present during the soft state, while such residuals are not significant during the hard state. The hard state spectra are characterised by a hard ($\Gamma\sim1.9-2.1$) power law, showing no evidence for a high energy cut off ($kT_e>60-140$ keV) and implying a small optical depth ($\tau<1.6$). The new observations confirm the previously witnessed trend of exhibiting strong Fe K absorption in the soft state, that significantly weakens during the hard state. Optical (GROND) and radio (GMRT) observations suggest for AX J1745.6-2901 a standard broad band SED as typically observed in accreting neutron stars.

Near-infrared counterparts of three transient very faint neutron star X-ray binaries

We present near-infrared (NIR) imaging observations of three transient neutron star X-ray binaries, SAX J1753.5-2349, SAX J1806.5-2215 and AX J1754.2-2754. All three sources are members of the class of `very faint' X-ray transients which exhibit X-ray luminosities $L_X\lesssim10^{36}$ erg s$^{-1}$. The nature of this class of sources is still poorly understood. We detect NIR counterparts for all three systems and perform multi-band photometry for both SAX J1753.5-2349 and SAX J1806.5-2215, including narrow-band Br$_{\gamma}$ photometry for SAX J1806.5-2215. We find that SAX J1753.5-2349 is significantly redder than the field population, indicating that there may be absorption intrinsic to the system, or perhaps a jet is contributing to the infrared emission. SAX J1806.5-2215 appears to exhibit absorption in Br$_{\gamma}$, providing evidence for hydrogen in the system. Our observations of AX J1754.2--2754 represent the first detection of a NIR counterpart for this system. We find that none of the measured magnitudes are consistent with the expected quiescent magnitudes of these systems. Assuming that the infrared radiation is dominated by either the disc or the companion star, the observed magnitudes argue against an ultracompact nature for all three systems.

Late-time Cooling of Neutron Star Transients and the Physics of the Inner Crust

Abstract An accretion outburst onto a neutron star transient heats the neutron star’s crust out of thermal equilibrium with the core. After the outburst, the crust thermally relaxes toward equilibrium with the neutron star core, and the surface thermal emission powers the quiescent X-ray light curve. Crust cooling models predict that thermal equilibrium of the crust will be established into quiescence. Recent observations of the cooling neutron star transient MXB 1659-29, however, suggest that the crust did not reach thermal equilibrium with the core on the predicted timescale and continued to cool after into quiescence. Because the quiescent light curve reveals successively deeper layers of the crust, the observed late-time cooling of MXB 1659-29 depends on the thermal transport in the inner crust. In particular, the observed late-time cooling is consistent with a low thermal conductivity layer near the depth predicted for nuclear pasta that maintains a temperature gradient between the neutron star’s inner crust and core for thousands of days into quiescence. As a result, the temperature near the crust–core boundary remains above the critical temperature for neutron superfluidity, and a layer of normal neutrons forms in the inner crust. We find that the late-time cooling of MXB 1659-29 is consistent with heat release from a normal neutron layer near the crust–core boundary with a long thermal time. We also investigate the effect of inner crust physics on the predicted cooling curves of the accreting transient KS 1731-260 and the magnetar SGR 1627-41.