Compositional Freeze-Out of Neutron Star Crusts

We have investigated the crustal properties of neutron stars without fallback accretion. We have calculated the chemical evolution of the neutron star crust in three different cases (a modified Urca process without the thermal influence of a crust, a thick crust, and a direct Urca process with a thin crust) in order to determine the detailed composition of the envelope and atmosphere as the nuclear reactions freeze out. Using a nuclear reaction network up to technetium, we calculate the distribution of nuclei at various depths of the neutron star. The nuclear reactions quench when the cooling timescale is shorter than the inverse of the reaction rate. Trace light elements among the calculated isotopes may have enough time to float to the surface before the layer crystallizes and form the atmosphere or envelope of the neutron star. The composition of the neutron-star envelope determines the total photon flux from the surface, and the composition of the atmosphere determines the emergent spectrum. Our calculations using each of the three cooling models indicate that without accretion of fallback the neutron star atmospheres are dependent on the assumed cooling process of the neutron star. Each of the cooling methods have different elements composing the atmosphere: for the modified Urca process the atmosphere is $^{28}$Si, the thick crust has an atmosphere of $^{50}$Cr, and the thin crust has an atmosphere of $^{40}$Ca. In all three cases the atmospheres are composed of elements which are lighter then iron.

💡 Research Summary

The paper investigates the composition of neutron‑star crusts and atmospheres in the absence of fallback accretion by modeling three distinct cooling scenarios: a modified Urca process (with negligible crustal heating), a thick‑crust model, and a direct Urca process coupled with a thin crust. Using a nuclear reaction network that follows isotopes up to technetium (Z = 43), the authors compute the temperature‑density evolution of the star and track the abundances of hundreds of nuclei at various depths. The key physical criterion is the “freeze‑out” condition: when the cooling timescale becomes shorter than the inverse of a given nuclear reaction rate, the reaction effectively ceases and the existing isotopic distribution is frozen in place.

After freeze‑out, light nuclei can buoy upward through the liquid layer. The authors estimate buoyancy velocities based on the balance of gravitational acceleration, ionic charge‑to‑mass ratio, and the viscosity of the surrounding plasma. If the buoyancy time is shorter than the crystallization time of the solid crust, the element reaches the surface and becomes part of the outer envelope or atmosphere. This competition determines which species dominate the observable layers.

The three cooling models yield markedly different atmospheric compositions. In the modified Urca case, silicon‑28 (^28Si) has sufficient time to float before the layer solidifies, making it the primary constituent of the atmosphere. In the thick‑crust scenario, chromium‑50 (^50Cr) dominates, while the thin‑crust/direct‑Urca model produces an atmosphere rich in calcium‑40 (^40Ca). All three atmospheres consist of elements lighter than iron, contrary to the common assumption that neutron‑star surfaces are iron‑dominated.

These compositional differences have direct implications for the emergent photon spectrum. Light elements possess distinct bound‑free edges (e.g., Si K‑edge at ~1.84 keV, Ca K‑edge at ~4.0 keV, Cr K‑edge at ~5.4 keV) and different opacities, which modify the shape of the X‑ray continuum and imprint characteristic absorption features. Consequently, the thermal photon flux and spectral shape observed from isolated neutron stars can be used to infer the underlying cooling mechanism and crustal physics.

The study also highlights the sensitivity of the final composition to the thermal history. The modified Urca process cools relatively slowly, allowing extended nuclear processing that favors silicon production. The thick crust retains heat longer, enabling more extensive alpha‑capture chains that culminate in chromium. The thin crust with direct Urca cools rapidly, truncating the reaction flow and leaving calcium as the most abundant light nucleus.

While the reaction network stops at technetium, the authors note that inclusion of heavier nuclei (e.g., silver, gold) could reveal additional abundance peaks, especially in scenarios with residual fallback material. Future work is suggested to couple the present compositional models with detailed radiative‑transfer calculations and to compare the predicted spectral signatures with high‑resolution X‑ray observations (e.g., from NICER or Athena). Moreover, exploring mixed scenarios where a small amount of fallback accretion coexists with the cooling‑driven freeze‑out would clarify the robustness of the light‑element atmospheres.

In summary, the paper demonstrates that the atmospheric composition of a neutron star without fallback accretion is not universal but depends critically on the cooling pathway. Modified Urca, thick‑crust, and thin‑crust/direct‑Urca models each produce distinct light‑element atmospheres (Si, Cr, Ca respectively), which in turn affect the observable photon flux and spectral features. This work provides a framework for using surface spectroscopy as a diagnostic of neutron‑star interior physics and cooling processes.