Search

We present analysis of multiple Chandra and XMM-Newton spectra, separated by 9-19 years, of four of the youngest central compact objects (CCOs) with ages < 2500 yr: CXOU J232327.9+584842 (Cassiopeia A), CXOU J160103.1-513353 (G330.2+1.0), 1WGA J1713.4-3949 (G347.3-0.5), and XMMU J172054.5-372652 (G350.1-0.3). By fitting these spectra with thermal models, we attempt to constrain each CCO's long-term cooling rate, composition, and magnetic field. For the CCO in Cassiopeia A, 14 measurements over 19 years indicate a decreasing temperature at a ten-year rate of 2.2+/-0.2 or 2.8+/-0.3 percent (1sigma error) for a constant or changing X-ray absorption, respectively. We obtain cooling rate upper limits of 17 percent for CXOU J160103.1-513353 and 6 percent for XMMU J172054.5-372652. For the oldest CCO, 1WGA J1713.4-3949, its temperature seems to have increased by 4+/-2 percent over a ten year period. Assuming each CCO's preferred distance and an emission area that is a large fraction of the total stellar surface, a non-magnetic carbon atmosphere spectrum is a good fit to spectra of all four CCOs. If distances are larger and emission areas are somewhat smaller, then equally good spectral fits are obtained using a hydrogen atmosphere with B <= 7x10^10 G or B >= 10^12 G for CXOU J160103.1-513353, B <= 10^10 G or B >= 10^12 G for XMMU J172054.5-372652, and non-magnetic hydrogen atmosphere for 1WGA J1713.4-3949. In a unified picture of CCO evolution, our results suggest most CCOs, and hence a sizable fraction of young neutron stars, have a surface magnetic field that is low early in their life but builds up over several thousand years., Comment: 15 pages, 9 figures; accepted for publication in MNRAS

We present a new model-independent (applicable for a broad range of equations of state) analysis of the neutrino emissivity due to triplet neutron pairing in neutron star cores. We find that the integrated neutrino luminosity of the Cooper Pair Formation (CPF) process can be written as a product of two factors. The first factor depends on the neutron star mass, radius and maximal critical temperature of neutron pairing in the core, $T_{Cn \mathrm{max}}$, but not on the particular superfluidity model; it can be expressed by an analytical formula valid for many nucleon equations of state. The second factor depends on the shape of the critical temperature profile within the star, the ratio of the temperature $T$ to $T_{Cn \mathrm{max}}$, but not on the maximal critical temperature itself. While this second factor depends on the superfluidity model, it obeys several model-independent constraints. This property allows one to analyse the thermal evolution of neutron stars with superfluid cores without relying on a specific model of their interiors. The constructed expressions allow us to perform a self-consistent analysis of spectral data and neutron star cooling theory. We apply these findings to the cooling neutron star in the Cassiopeia A supernova remnant using 14 sets of observations taken over 19 years. We constrain $T_{Cn\mathrm{max}}$ to the range of $ (5-10)\times 10^8$ K. This value depends weakly on the equation of state and superfluidity model, and will not change much if cooling is slower than the current data suggest. We also constrain the overall efficiency of the CPF neutrino luminosity., Comment: 18 pages, 14 figures, 6 tables; accepted for publication in MNRAS

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies., Comment: 10 pages, 10 figures; Published in MNRAS

The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase., Comment: 7 pages, 4 figures; accepted for publication in MNRAS

The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star., Comment: 7 pages, 7 figures; to appear in the AIP Conference Proceedings of the Xiamen-CUSTIPEN Workshop on the EOS of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy (January 3-7, 2019, Xiamen, China)

A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns y$_H$ > 10$^{7}$ g cm$^{-2}$, while DNB can make helium envelopes more transparent to heat flux for surface temperatures $T_s$ > 2 $\times 10^6$ K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the long-term cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of $\approx$ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope., Comment: 16 pages, 15 figures, 4 tables; accepted for publication in MNRAS

Monitoring the cooling of neutron-star crusts heated during accretion outbursts allows us to infer the physics of the dense matter present in the crust. We examine the crust cooling evolution of the low-mass X-ray binary MXB 1659-29 up to ~505 days after the end of its 2015 outburst (hereafter outburst II) and compare it with what we observed after its previous 1999 outburst (hereafter outburst I) using data obtained from the Swift, XMM-Newton, and Chandra observatories. The observed effective surface temperature of the neutron star in MXB 1659-29 dropped from ~92 eV to ~56 eV from ~12 days to ~505 days after the end of outburst II. The most recently performed observation after outburst II suggests that the crust is close to returning to thermal equilibrium with the core. We model the crust heating and cooling for both its outbursts collectively to understand the effect of parameters that may change for every outburst (e.g., the average accretion rate, the length of outburst, the envelope composition of the neutron star at the end of the outburst) and those which can be assumed to remain the same during these two outbursts (e.g., the neutron star mass, its radius). Our modelling indicates that all parameters were consistent between the two outbursts with no need for any significant changes. In particular, the strength and the depth of the shallow heating mechanism at work (in the crust) were inferred to be the same during both outbursts, contrary to what has been found when modelling the cooling curves after multiple outburst of another source, MAXI J0556-332. This difference in source behaviour is not understood. We discuss our results in the context of our current understanding of cooling of accretion-heated neutron-star crusts, and in particular with respect to the unexplained shallow heating mechanism., Comment: Submitted to A&A. The supplementary video can be found at https://www.youtube.com/watch?v=OpJ053zq9-M

We present a new Chandra observation (performed in July 2016) of the neutron star X-ray transient IGR J17480-2446, located in the globular cluster Terzan 5. We study the continued cooling of the neutron star crust in this system that was heated during the 2010 outburst of the source. This new observation was performed two years after the last observation of IGR J17480-2446, hence, significantly extending the cooling baseline. We reanalysed all available Chandra observations of the source (but excluding observations during which one of the known transients in Terzan 5 was in outburst) and fitted the obtained cooling curve with our cooling code NSCool, which allows for much improved modelling than what was previously performed for the source. The data and our fit models indicate that the crust was still cooling ~5.5 years after the outburst ended. The neutron star crust has likely not reached crust-core thermal equilibrium yet, and further cooling is predicted (which can be confirmed with additional Chandra observations in >5 years). Intriguingly, we find indications that the thermal conductivity might be relatively low in part of the crust compared to what has been inferred for other crust-cooling sources and tentatively suggest that this layer might be located around the neutron drip. The reason for this difference is unclear, but might be related to the fact that IGR J17480-2446 harbours a relatively slowly rotating neutron star (with a spin of 11 Hz) that has a relatively strong inferred surface magnetic field ($10^{9-10}$ Gauss) compared to what is known or typically assumed for other cooling sources., Comment: 17 pages, 10 figures, 4 tables, accepted for publication in MNRAS

In accreting neutron star X-ray transients, the neutron star crust can be substantially heated out of thermal equilibrium with the core during an accretion outburst. The observed subsequent cooling in quiescence (when accretion has halted) offers a unique opportunity to study the structure and thermal properties of the crust. Initially crust cooling modelling studies focussed on transient X-ray binaries with prolonged accretion outbursts (> 1 year) such that the crust would be significantly heated for the cooling to be detectable. Here we present the results of applying a theoretical model to the observed cooling curve after a short accretion outburst of only ~10 weeks. In our study we use the 2010 outburst of the transiently accreting 11 Hz X-ray pulsar in the globular cluster Terzan 5. Observationally it was found that the crust in this source was still hot more than 4 years after the end of its short accretion outburst. From our modelling we found that such a long-lived hot crust implies some unusual crustal properties such as a very low thermal conductivity (> 10 times lower than determined for the other crust cooling sources). In addition, we present our preliminary results of the modelling of the ongoing cooling of the neutron star in MXB 1659-298. This transient X-ray source went back into quiescence in March 2017 after an accretion phase of ~1.8 years. We compare our predictions for the cooling curve after this outburst with the cooling curve of the same source obtained after its previous outburst which ended in 2001., Comment: 4 pages, 1 figure, to appear in the proceedings of "IAUS 337: Pulsar Astrophysics - The Next 50 Years" eds: P. Weltevrede, B.B.P. Perera, L. Levin Preston & S. Sanidas

The dwarf stars in the 26 year period binary alpha Com were predicted to eclipse each other in early 2015. That prediction was based on an orbit model made with over 600 astrometric observations using micrometers, speckle interferometry, and long baseline optical interferometry. Unfortunately, it has been realized recently that the position angle measurements for three of the observations from ~100 years ago were in error by 180 degrees, which skewed the orbital fit. The eclipse was likely 2 months earlier than predicted, at which point the system was low on the horizon at sunrise., Comment: 3 pages, to be submitted to ApJL