An official website of the United States government
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
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
- PAR ID:
- 10362680
- 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
-
-
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
-
The effect of isospin-dependent nuclear forces on the inner crust of neutron stars is modeled within the framework of Quantum Molecular Dynamics (QMD). To successfully control the density dependence of the symmetry energy of neutron-star matter below nuclear saturation density, a mixed vector-isovector potential is introduced. This approach is inspired by the baryon density and isospin density-dependent repulsive Skyrme force of asymmetric nuclear matter. In isospin-asymmetric nuclear matter, the system shows nucleation, as nucleons are arranged into shapes resembling nuclear pasta. The dependence of clusterization in the system on the isospin properties is also explored by calculating two-point correlation functions. We show that, as compared to previous results that did not involve such mixed interaction terms, the energy symmetry slope L is successfully controlled by varying the corresponding coupling strength. Nevertheless, the effect of changing the slope of the nuclear symmetry energy L on the crust-core transition density does not seem significant. To the knowledge of the authors, this is the first implementation of such a coupling in a QMD model for isospin asymmetric matter, which is relevant to the inner crust of neutron and proto-neutron stars.more » « less
-
Abstract An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm−3and 0.1 fm−3, matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.more » « less
-
We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2 M_Sun neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above 1.8 M_Sun.more » « less
