Light nuclei in supernova envelopes: a quasiparticle gas model

We present an equation of state and the composition of low-density supernova matter composed of light nuclei with mass number A \le 13. We work within the quasiparticle gas model, which accounts for bound states with decay time scales larger than the relevant time scale of supernova and protoneutron star evolution. The mean-field contribution is included in terms of Skyrme density functional. Deuterons, tritons, and 3H(e) nuclei appear in matter in concentrations that are substantially higher than those of heavier nuclei. We calculate the critical temperature of deuteron condensation in such matter, and demonstrate that the appearance of clusters substantially lowers the critical temperature.

💡 Research Summary

This paper develops a comprehensive equation of state (EOS) and composition model for low‑density supernova matter that explicitly includes light nuclei with mass numbers up to A = 13. The authors adopt a quasiparticle gas model (QGM), treating each bound nuclear species as a quasiparticle whose lifetime exceeds the characteristic dynamical timescales of a core‑collapse supernova and the subsequent proto‑neutron star (PNS) evolution (milliseconds to seconds). In this framework, the statistical occupation of each species is governed by Boltzmann factors that incorporate the chemical potential, binding energy, and temperature, while the mean‑field interaction among nucleons is described by a Skyrme density functional. The Skyrme parametrization supplies density‑ and temperature‑dependent effective masses and mean‑field potentials, ensuring that the chemical equilibrium conditions among neutrons, protons, and the various clusters are satisfied self‑consistently.

The authors first derive the chemical equilibrium relations for the light clusters (μ_d = μ_n + μ_p, μ_t = 2μ_n + μ_p, μ_³He = μ_n + 2μ_p, etc.) and then solve for the abundances Y_i as functions of baryon density ρ_B (10⁻⁴–10⁻² fm⁻³) and temperature T (1–10 MeV). Their calculations reveal that deuterons (d), tritons (t), and helium‑3 (³He) dominate the composition in this regime, reaching mass fractions that are one to two orders of magnitude larger than those of heavier nuclei (A ≥ 4). For example, at ρ_B ≈ 10⁻³ fm⁻³ and T ≈ 5 MeV, deuterons account for roughly 30 % of the total baryonic mass. This high concentration of light clusters has several astrophysical implications: it enhances neutron–proton recombination rates, modifies the electron fraction Y_e, and can affect neutrino opacity through altered charged‑current reaction rates.

A central focus of the work is the possibility of deuteron condensation (Bose–Einstein condensation of deuterons) in supernova matter. By evaluating the deuteron chemical potential as a function of temperature and density, the authors identify a critical temperature T_c at which μ_d reaches zero, signalling the onset of macroscopic occupation of the deuteron ground state. In a medium that includes only free nucleons, T_c is found to be of order 1 MeV. However, when the full spectrum of light clusters is present, the critical temperature drops dramatically to ≲ 0.5 MeV. The reduction arises because the formation of clusters lowers the effective deuteron binding energy and enhances deuteron–deuteron interactions, thereby requiring less thermal energy to achieve condensation. This result demonstrates that the presence of clusters substantially suppresses the temperature window for deuteron superfluidity in the supernova envelope.

The inclusion of light clusters also softens the EOS. Compared with a traditional nucleonic EOS that neglects clusters, the QGM EOS exhibits lower pressure at a given density and a reduced incompressibility modulus. This softening can influence the dynamics of core bounce, the propagation of the shock wave, and the subsequent PNS cooling phase. Moreover, the altered composition changes the electron capture rates (e⁻ + p → n + ν_e) because the abundance of free protons is reduced in favor of bound clusters, potentially affecting the lepton‑richness of the early ejecta.

The paper emphasizes the flexibility of the QGM approach. While the present study limits the cluster spectrum to A ≤ 13, the formalism can be extended to include α‑particles, ⁶Li, ⁷Be, and even heavier nuclei, as well as more sophisticated nucleon‑nucleon interactions derived from chiral effective field theory. Adjusting the Skyrme parameters to reflect different nuclear matter constraints would allow the model to be calibrated for a variety of astrophysical simulations, ranging from multidimensional supernova explosion models to long‑term PNS cooling calculations.

In summary, the authors provide a quantitatively robust EOS that incorporates light nuclear clusters via a quasiparticle gas model and a Skyrme mean‑field. Their results show that deuterons, tritons, and helium‑3 dominate the composition of low‑density supernova matter, that these clusters dramatically lower the critical temperature for deuteron condensation, and that the overall EOS becomes softer. These findings have direct implications for supernova dynamics, neutrino transport, and the thermal evolution of proto‑neutron stars, and they establish a solid foundation for future work that seeks to integrate a richer nuclear composition into astrophysical simulations.