Compositionally-driven convection in the oceans of accreting neutron stars

We discuss the effect of chemical separation as matter freezes at the base of the ocean of an accreting neutron star, and argue that the retention of light elements in the liquid acts as a source of buoyancy that drives a slow but continual mixing of the ocean, enriching it substantially in light elements, and leading to a relatively uniform composition with depth. We first consider the timescales associated with different processes that can redistribute elements in the ocean, including convection, sedimentation, crystallization, and diffusion. We then calculate the steady state structure of the ocean of a neutron star for an illustrative model in which the accreted hydrogen and helium burns to produce a mixture of O and Se. Even though the H/He burning produces only 2% oxygen by mass, the steady state ocean has an oxygen abundance more than ten times larger, almost 40% by mass. Furthermore, we show that the convective motions transport heat inwards, with a flux of ~ 0.2 MeV per nucleon for an O-Se ocean, heating the ocean and steepening the outwards temperature gradient. The enrichment of light elements and heating of the ocean due to compositionally-driven convection likely have important implications for carbon ignition models of superbursts.

💡 Research Summary

The paper investigates a previously under‑appreciated mixing process in the liquid “ocean” that overlies the solid crust of an accreting neutron star. When the ocean material freezes at its base, the solid phase preferentially incorporates the heavier nuclei, leaving the lighter species (e.g., oxygen) in the residual liquid. This compositional segregation creates a buoyancy source that drives a slow, continuous convection distinct from ordinary thermal convection. The authors first quantify the characteristic timescales of four competing transport mechanisms: convection, sedimentation, crystallization, and diffusion. Using standard Boussinesq approximations, they estimate convective velocities and show that, for the relevant density and gravity, compositional buoyancy can generate overturn on timescales far shorter than sedimentation or diffusion, making convection the dominant redistributive process.

To illustrate the effect, they construct a simple one‑dimensional model in which the accreted hydrogen/helium burns to a mixture of oxygen (O) and selenium (Se). Nuclear burning produces only about 2 % O by mass, the rest being heavy Se. As the ocean cools and the base solidifies, Se is incorporated into the crystal lattice while O remains in the liquid, establishing a strong compositional gradient. The resulting buoyancy triggers convection that mixes the entire ocean column. Solving for a steady‑state structure, they find that the oxygen mass fraction rises dramatically—from the initial 2 % up to roughly 40 %—a factor of twenty enrichment. The convective motions also transport heat inward. By evaluating the convective heat flux they obtain an inward energy transport of ≈0.2 MeV per nucleon for the O‑Se ocean, which steepens the outward temperature gradient and raises the temperature of the deep ocean layers.

The authors discuss the astrophysical implications of this enriched, hotter ocean. Superbursts—rare, energetic X‑ray bursts thought to be triggered by carbon ignition in the ocean—are highly sensitive to both the local composition and temperature. An ocean that is both richer in light nuclei and hotter due to composition‑driven convection will alter the ignition depth, the required carbon mass fraction, and the recurrence time of superbursts. Moreover, the enhanced light‑element fraction modifies the electron conductivity and the overall thermal conductivity of the ocean, feeding back into the star’s cooling curve and the observable quiescent emission.

Finally, the paper outlines observational avenues for testing the theory. Changes in burst light curves, the cooling tail after a superburst, and spectral signatures of surface composition could all carry imprints of the deep compositional convection. High‑precision X‑ray timing and spectroscopy with current (e.g., NICER, XMM‑Newton) and upcoming missions (e.g., Athena) may be able to detect the predicted temperature and compositional signatures.

In summary, the work identifies a robust mechanism—chemical‑separation‑induced buoyancy—that drives convection in neutron‑star oceans, leading to substantial enrichment of light elements and additional internal heating. These effects are likely to be crucial for realistic models of superburst ignition, crust‑core thermal coupling, and the interpretation of X‑ray observations of accreting neutron stars.