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× 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 longterm 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 ≈ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.

Transient non-Gaussian noise in gravitational wave detectors, commonly referred to as glitches, pose challenges for detection and inference of the astrophysical properties of detected signals when the two are coincident in time. Current analyses aim toward modeling and subtracting the glitches from the data using a flexible, morphology-independent model in terms of sine-Gaussian wavelets before the signal source properties are inferred using templates for the compact binary signal. We present a new analysis of gravitational wave data that contain both a signal and glitches by simultaneously modeling the compact binary signal in terms of templates and the instrumental glitches using sine-Gaussian wavelets. The model for the glitches is generic and can thus be applied to a wide range of glitch morphologies without any special tuning. The simultaneous modeling of the astrophysical signal with templates allows us to efficiently separate the signal from the glitches, as we demonstrate using simulated signals injected around real O2 glitches in the two LIGO detectors. We show that our new proposed analysis can separate overlapping glitches and signals, estimate the compact binary parameters, and provide ready-to-use glitchsubtracted data for downstream inference analyses.

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.

No abstract

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, TCnmax, 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 TCnmax, 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 TCnmax to the range of (5 − 10) × 108 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.

No abstract

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 <2500yr: 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 (1σ 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 ≤ 7 × 1010G or B ≥ 1012G for CXOU J160103.1−513353 and B ≤ 1010G or B ≥ 1012G 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.

The increasing sensitivity of gravitational-wave detectors has brought about an increase in the rate of astrophysical signal detections as well as the rate of "glitches"; transient and non-Gaussian detector noise. Temporal overlap of signals and glitches in the detector presents a challenge for inference analyses that typically assume the presence of only Gaussian detector noise. In this study we perform an extensive exploration of the efficacy of a recently proposed method that models the glitch with sine-Gaussian wavelets while simultaneously modeling the signal with compact-binary waveform templates. We explore a wide range of glitch families and signal morphologies and demonstrate that the joint modeling of glitches and signals (with wavelets and templates respectively) can reliably separate the two. We find that the glitches that most affect parameter estimation are also the glitches that are well modeled by such wavelets due to their compact time-frequency signature. As a further test, we investigate the robustness of this analysis against waveform systematics like those arising from the exclusion of higher-order modes and spin-precession effects. Our analysis provides an estimate of the signal parameters; the glitch waveform to be subtracted from the data and an assessment of whether some detected excess power consists of a glitch, signal, or both. We analyze the low-significance triggers (191225_215715 and 200114_020818) and find that they are both consistent with glitches overlapping high-mass signals.