Nuclear Pasta and Crustal Quasi-Periodic Oscillations in Neutron Star

We investigate the impact of nuclear pasta on crustal structure and torsional oscillations using a Bayesian ensemble of unified neutron-star equations of state based on relativistic mean-field models constrained by nuclear experiments, empirical saturation properties, chiral effective field theory, and multimessenger observations. For each posterior sample, we compute the pasta sequence within a compressible liquid-drop model and quantify the onset density, thickness, and mass fraction of the pasta layers. We show that the appearance and extent of nuclear pasta are primarily controlled by the symmetry-energy slope parameter $L$. While spherical and rod-like pasta configurations are present for all equations of state, only a small fraction of the posterior supports slab, tube, or bubble geometries. The transition from spherical nuclei to rods is tightly constrained to occur at a density of $ρ_{\rm sr} = 0.0588^{+0.0045}{-0.0065},\mathrm{fm^{-3}}$. We further predict that the nuclear pasta layer occupies a relative radial thickness of $ΔR{\rm pasta}/ΔR_{\rm c} = 0.140^{+0.025}{-0.036}$ and contributes a relative mass fraction of $ΔM{\rm pasta}/ΔM_{\rm c} = 0.475^{+0.071}{-0.113}$. Using the resulting crust models, we present the first quasi-periodic oscillations (QPOs) analysis based on a Bayesian posterior ensemble of neutron-star equations of state and systematically assess their compatibility with observed low-frequency quasi-periodic oscillations. We find that the predicted QPO frequencies are strongly correlated with the curvature of the symmetry energy evaluated at sub-saturation density, $K{\rm sym}(ρ_0/2)$, and that uncertainties in the equation of state translate into a range of angular indices $\ell$ consistent with the observed frequencies.

💡 Research Summary

This paper presents a comprehensive Bayesian investigation of nuclear pasta in neutron‑star crusts and its impact on crustal torsional quasi‑periodic oscillations (QPOs). The authors start from a posterior ensemble of unified neutron‑star equations of state (EOS) generated with relativistic mean‑field (RMF) models that are constrained by a wide range of nuclear‑experimental data, empirical saturation properties, chiral effective‑field theory at sub‑saturation densities, and multimessenger astrophysical observations. Approximately 40 000 EOS samples are drawn from this posterior.

For each EOS, the inner‑crust composition is calculated using a compressible liquid‑drop model (CLDM) within the Wigner‑Seitz approximation. Five canonical geometries—spherical nuclei, cylindrical rods, planar slabs, cylindrical tubes, and spherical bubbles—are considered. At each baryon density the total energy (bulk, surface, curvature, Coulomb, and electron contributions) is evaluated for all geometries, and the energetically favored configuration is selected. This procedure yields, for every EOS, the full pasta sequence, the transition densities between successive shapes, and the overall thickness and mass fraction of the pasta layer.

The statistical analysis shows that the onset of non‑spherical pasta (the transition from spherical nuclei to rods) is tightly constrained to a density ρ_sr = 0.0588 +0.0045/‑0.0065 fm⁻³. The relative radial thickness of the pasta region is ΔR_pasta/ΔR_c = 0.140 +0.025/‑0.036, and its mass fraction is ΔM_pasta/ΔM_c = 0.475 +0.071/‑0.113. The presence and extent of pasta are found to be primarily governed by the symmetry‑energy slope parameter L; larger L values produce earlier onset and a more extensive pasta layer. While spherical and rod phases appear for virtually all EOS, slab, tube, and bubble phases are supported only by a minority of the posterior (≈30 %).

To assess the observable consequences, the authors compute torsional shear‑mode frequencies of the crust using a plane‑parallel slab approximation. The shear modulus for spherical nuclei follows the standard BCC lattice expression, μ = 0.1194 (1‑0.010 Z^{2/3}) n_i (Z e)² a. In the pasta region the shear modulus is smoothly suppressed via a quadratic function μ̄ = c₁(ρ‑ρ_c)²(ρ‑c₂), ensuring continuity at the pasta onset density and vanishing at the crust‑core transition. This phenomenological softening captures the expected reduction of rigidity in disordered pasta structures. Solving the resulting one‑dimensional eigenvalue problem yields the fundamental and overtone torsional frequencies for each EOS realization.

A key result is the strong correlation between the predicted QPO frequencies and the curvature of the symmetry energy at half saturation, K_sym(ρ₀/2). The EOS uncertainties translate into a spread of angular quantum numbers ℓ that can reproduce the observed low‑frequency QPOs (≈30–150 Hz) from magnetar giant flares. In other words, rather than a single ℓ value, a range ℓ ≈ 2–10 becomes compatible once pasta‑induced softening and EOS variance are accounted for. This demonstrates that nuclear pasta can significantly lower crustal torsional frequencies and modify mode spacing, providing a plausible mechanism for the observed QPO spectrum.

Overall, the paper delivers three major contributions: (1) a statistically robust quantification of pasta onset density, thickness, and mass fraction grounded in a Bayesian EOS posterior; (2) a physically motivated model for the reduction of shear rigidity in the pasta layer; and (3) a direct link between these microphysical properties and observable QPO frequencies, offering an indirect but testable probe of nuclear pasta in neutron stars. The methodology sets a new standard for propagating nuclear‑physics uncertainties through to astrophysical observables.