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.
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 code
- PAR ID:
- 10463625
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 955
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 51
- Size(s):
- Article No. 51
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
ABSTRACT Evidence has accumulated for an as-yet unaccounted for source of heat located at shallow depths within the accreted neutron star crust. However, the nature of this heat source is unknown. I demonstrate that the inferred depth of carbon ignition in X-ray superbursts can be used as an additional constraint for the magnitude and depth of shallow heating. The inferred shallow heating properties are relatively insensitive to the assumed crust composition and carbon fusion reaction rate. For low-accretion rates, the results are weakly dependent on the duration of the accretion outburst, so long as accretion has ensued for enough time to replace the ocean down to the superburst ignition depth. For accretion rates at the Eddington rate, results show a stronger dependence on the outburst duration. Consistent with earlier work, it is shown that urca cooling does not impact the calculated superburst ignition depth unless there is some proximity in depth between the heating and cooling sources.
-
more » « less
Abstract At the low temperature and high density conditions of a neutron star crust neutrons are degenerate. In this work, we study the effect of this degeneracy on the capture rates of neutrons on neutron rich nuclei in accreted crusts. We use a statistical Hauser–Feshbach model to calculate neutron capture rates and find that neutron degeneracy can increase rates significantly. Changes increase from a factor of a few to many orders of magnitude near the neutron drip line. We also quantify uncertainties due to model inputs for masses,
γ -strength functions, and level densities. We find that uncertainties increase dramatically away from stability and that degeneracy tends to increase these uncertainties further, except for cases near the neutron drip line where degeneracy leads to more robustness. As in the case of capture of classically distributed neutrons, variations in the mass model have the strongest impact. Corresponding variations in the reaction rates can be as high as 3–4 orders of magnitude, and be more than 5 times larger than under classical conditions. To ease the incorporation of neutron degeneracy in nucleosynthesis networks, we provide tabulated results of capture rates as well as analytical expressions as function of temperature and neutron chemical potential, for proton numbers between 3 ≤Z ≤ 85, derived from fits to our numerical results. Fits are based on a new parametrization that complements previously employed power law approximations with additional Lorentzian terms that account for low energy resonances, significantly improving accuracy. -
The formation of continental crust in magmatic arcs involves cooling of hot magmas to a relatively colder crust enhanced by exhumation and hydrothermal circulation in the upper crust. To quantify the influence of these processes on the thermal and rheological states of the crust, we developed a one-dimensional thermal evolution model, which invokes conductive cooling, advection of crust by erosion-driven exhumation, and cooling by hydrothermal circulation. We parameterized hydrothermal cooling by adopting depth-dependent effective thermal conductivity, which is determined by the crustal permeability structure and the prescribed Nusselt number at the surface. Different combinations of erosion rate and Nusselt number were tested to study the evolution of crustal geotherms, surface heat flux, and cooling rate. Simulations and scaling analyses quantify how erosion and hydrothermal circulation promote cooling via increasing total surface heat flux compared to pure conductive cooling. Hydrothermal circulation imposes intense short-term and persistent long-term cooling effects. Thinner, warmer, fast exhuming crust, with higher permeability and more vigorous hydrothermal circulation, leads to higher steady-state total surface heat flux. Hydrothermal cooling at steady state is more effective when the Péclet number is small. Hydrothermal cooling also changes crustal rheological state and thickens the brittle crust. This in turn promotes the initiation of brittle deformation in the upper crust in magmatic arcs or in regions undergoing exhumation. Interpretation of low-temperature thermochronological data could overestimate average cooling rates if hydrothermal cooling is not considered.more » « less
-
more » « less
ABSTRACT We present an improved version of the 3D Monte Carlo radiative transfer code possis to model kilonovae from neutron star mergers, wherein nuclear heating rates, thermalization efficiencies, and wavelength-dependent opacities depend on local properties of the ejecta and time. Using an axially symmetric two-component ejecta model, we explore how simplistic assumptions on heating rates, thermalization efficiencies, and opacities often found in the literature affect kilonova spectra and light curves. Specifically, we compute five models: one (FIDUCIAL) with an appropriate treatment of these three quantities, one (SIMPLE-HEAT) with uniform heating rates throughout the ejecta, one (SIMPLE-THERM) with a constant and uniform thermalization efficiency, one (SIMPLE-OPAC) with grey opacities, and one (SIMPLE-ALL) with all these three simplistic assumptions combined. We find that deviations from the FIDUCIAL model are of several (∼1–10) magnitudes and are generally larger for the SIMPLE-OPAC and SIMPLE-ALL compared to the SIMPLE-THERM and SIMPLE-HEAT models. The discrepancies generally increase from a face-on to an edge-on view of the system, from early to late epochs and from infrared to ultraviolet/optical wavelengths. This work indicates that kilonova studies using either of these simplistic assumptions ought to be treated with caution and that appropriate systematic uncertainties ought to be added to kilonova light curves when performing inference on ejecta parameters.
