skip to main content

Title: The Impact of Neutron Transfer Reactions on the Heating and Cooling of Accreted Neutron Star Crusts

Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and long-term cooling in transiently accreting systems. It was recently suggested that previously ignored neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.

more » « less
Award ID(s):
2011890 1913554
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 205
["Article No. 205"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution codedStarto further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

    more » « less
  2. Abstract Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on the accreting layers of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modeling for observed X-ray light curves. In this study, we investigate theoretical X-ray burst models compared with burst light curves with GS 1826-24 observations. We focus on the impacts of the NS mass and radius and base heating on the NS surface using the MESA code. We find a monotonic correlation between the NS mass and the parameters of the light curve. The higher the mass, the longer the recurrence time and the greater the peak luminosity. While the larger the radius, the longer the recurrence time, the peak luminosity remains nearly constant. In the case of increasing base heating, both the recurrence time and peak luminosity decrease. We also examine the above results with a different numerical code, HERES , based on general relativity and consider the central NS. We find that the burst rate, energy, and strength are almost the same in two X-ray burst codes by adjusting the base heat parameter in MESA (the relative errors ≲5%), while the duration and rise times are significantly different between (the relative error is possibly ∼50%). The peak luminosity and the e-folding time change irregularly between two codes for different accretion rates. 
    more » « less

    Three sudden spin-down events, termed ‘antiglitches’, were recently discovered in the accreting pulsar NGC 300 ULX-1 by the Neutron Star Interior Composition Explorer mission. Unlike previous antiglitches detected in decelerating magnetars, these are the first antiglitches recorded in an accelerating pulsar. One standard theory is that pulsar spin-up glitches are caused by avalanches of collectively unpinning vortices that transfer angular momentum from the superfluid interior to the crust of a neutron star. Here, we test whether vortex avalanches are also consistent with the antiglitches in NGC 300 ULX-1, with the angular momentum transfer reversed. We perform N-body simulations of up to 5 × 103 pinned vortices in two dimensions in secularly accelerating and decelerating containers. Vortex avalanches routinely occur in both scenarios, propagating inwards and outwards, respectively. The implications for observables, such as size and waiting time statistics, are considered briefly.

    more » « less

    We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

    more » « less

    The evolution of accreting X-ray binary systems is closely coupled to the properties of their donor stars. Consequently, we can constrain the evolutionary track a system is by establishing the nature of its donor. Here, we present far-ultraviolet (far-UV) spectroscopy of the transient neutron-star low-mass X-ray binary J1858 in different accretion states (low-hard, high-hard, and soft). All of these spectra exhibit anomalous N v, C iv, Si iv, and He ii lines, suggesting that its donor star has undergone CNO processing. We also determine the donor’s effective temperature, Td ≃ 5700 K, and radius, Rd ≃ 1.7 R⊙, based on photometric observations obtained during quiescence. Lastly, we leverage the transient nature of the system to set an upper limit of $\dot{M}_{\rm acc} \lesssim 10^{-8.5}~{\rm M}_{\odot }~\mathrm{ yr}^{-1}$ on the present-day mass-transfer rate. Combining these with the orbital period of the system, Porb = 21.3 h, we search for viable evolution paths. The initial donor masses in the allowed solutions span the range 1 M⊙ ≲ Md,i ≲ 3.5 M⊙. All but the lowest masses in this range are consistent with the strong CNO-processing signature in the UV line ratios. The present-day donor mass in the permitted tracks are 0.5 M⊙ ≲ Md,obs ≲ 1.3 M⊙, higher than suggested by eclipse modelling. Since Porb is close to the so-called bifurcation period, both converging and diverging binary tracks are permitted. If the former is confirmed, J1858 will end its life as an ultracompact system with a substellar donor.

    more » « less