skip to main content


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

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
NSF-PAR ID:
10362680
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
925
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 205
Size(s):
["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
  3. ABSTRACT

    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
  4. ABSTRACT

    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
  5. ABSTRACT

    Theoretically predicted yields of elements created by the rapid neutron capture (r-)process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentally unknown neutron-rich nuclei. This includes two frameworks for calculating the radiative neutron capture rates and 14 different models for nuclear masses, β-decay rates, and fission properties. Our r-process nuclear network calculations are based on detailed hydrodynamical simulations of dynamically ejected material from NS–NS or NS–BH binary mergers plus the secular ejecta from BH–torus systems. The impact of nuclear uncertainties on the r-process abundance distribution and the early radioactive heating rate is found to be modest (within a factor of ∼20 for individual A > 90 abundances and a factor of 2 for the heating rate). However, the impact on the late-time heating rate is more significant and depends strongly on the contribution from fission. We witness significantly higher sensitivity to the nuclear physics input if only a single trajectory is used compared to considering ensembles with a much larger number of trajectories (ranging between 150 and 300), and the quantitative effects of the nuclear uncertainties strongly depend on the adopted conditions for the individual trajectory. We use the predicted Th/U ratio to estimate the cosmochronometric age of six metal-poor stars and find the impact of the nuclear uncertainties to be up to 2 Gyr.

     
    more » « less