Neutron star crust and outer core equation of state from chiral effective field theory with quantified uncertainties

This paper employs a Bayesian framework with a two-dimensional Gaussian process to quantify chiral effective field theory uncertainties in the equation of state for asymmetric nuclear matter up to twice saturation density, thereby constructing consistent models for the inner crust and outer core of neutron stars.

The Gist

Imagine a neutron star as the ultimate cosmic pressure cooker. Inside, matter is squeezed so tightly that atoms collapse, leaving behind a dense soup of neutrons and a few protons. To understand how this soup behaves—how hard it pushes back, how heavy it is, and how it changes under pressure—scientists need a "recipe" called an Equation of State (EOS).

This paper is about writing a new, highly precise recipe for the inner layers of a neutron star, while also attaching a "confidence meter" to every ingredient.

Here is the breakdown of what the scientists did, using everyday analogies:

1. The Problem: Guessing the Recipe

For decades, scientists have tried to calculate how nuclear matter behaves using a theory called Chiral Effective Field Theory (EFT). Think of EFT like a recipe that gets better the more steps you add.

• Step 1 (Low Order): You get a rough idea of the taste.
• Step 10 (High Order): You get a gourmet meal.
• The Catch: Every time you add a step, you introduce a tiny bit of uncertainty because you have to guess what the next step would look like.

The authors wanted to know: How much can we trust our recipe at different densities and different mixtures of neutrons and protons?

2. The Solution: A "Smart Predictor" (The Gaussian Process)

Instead of just calculating one single number for every condition, the team built a two-dimensional Gaussian Process (GP).

• The Analogy: Imagine you are trying to predict the weather. You have data for "Temperature" and "Humidity." A simple map might just show dots. A Gaussian Process is like a smart, flexible rubber sheet stretched over those dots. It doesn't just connect the dots; it learns the shape of the terrain between them.
• Why it's special here: Most previous maps only looked at "pure neutron weather" (like a desert) or "perfectly balanced weather" (like a tropical island). This new map covers the whole spectrum, from pure neutrons to mixtures with protons.
• The "Uncertainty" Feature: This rubber sheet doesn't just give you a temperature; it tells you, "I'm 95% sure the temperature is between 20°C and 22°C." This is crucial because it quantifies the "truncation error"—the guesswork involved in stopping the recipe at a certain step.

3. Refining the Ingredients (The Reference Energy)

When they first tried to stretch their rubber sheet, they noticed it was getting a bit wobbly, especially when there were more protons in the mix.

• The Fix: They realized their "baseline" (the reference energy) wasn't quite right. It was like trying to measure the height of a building using a ruler that shrinks when the sun gets hot.
• The Adjustment: They tweaked the baseline to account for the fact that three-neutron interactions (a specific type of nuclear force) become very important at high pressures. By adjusting this baseline, the "rubber sheet" smoothed out, and the predictions became much more consistent across all types of neutron-rich matter.

4. The Result: A Map of the Neutron Star Crust

With their new, smoothed-out, confidence-metered map, they looked at the inner crust of the neutron star. This is the layer just below the surface where the matter starts to get weird.

• The "Drip" Phenomenon: In the crust, neutrons are usually trapped inside atomic nuclei (like marbles in a jar). But as you go deeper, the pressure gets so high that neutrons start to "drip" out of the jar and float freely.
• The Surprise: The paper also found evidence for "Proton Drip." Just like neutrons, protons can eventually get squeezed out of their jars and float freely too.
• The Uncertainty: The team calculated exactly where this happens. They found that depending on how "stiff" or "soft" the nuclear matter is (based on their uncertainty bands), proton drip might happen in a very narrow range of depths, or it might almost disappear entirely.

5. Why This Matters (According to the Paper)

The authors didn't predict black holes or new medicines. Instead, they provided a consistent, mathematically rigorous framework.

• They showed that you can calculate the properties of neutron star matter (pressure, energy, chemical balance) without having to make up "fake" numbers to fill in the gaps between pure neutron matter and normal matter.
• They proved that even with the uncertainties in our current nuclear theories, we can still pinpoint specific phases of matter, like the "proton drip" phase, with a known level of confidence.

In short: The authors built a high-tech, flexible map of the inside of a neutron star. They didn't just draw the lines; they added a "fuzziness" indicator to show exactly where the map is blurry and where it is sharp, specifically focusing on the moment when protons and neutrons start to escape their atomic cages.

Technical Summary

Technical Summary: Neutron Star Crust and Outer Core Equation of State from Chiral Effective Field Theory with Quantified Uncertainties

Problem Statement
The nuclear equation of state (EOS) is fundamental to understanding the interior of neutron stars, core-collapse supernovae, and neutron star mergers. While recent advances in chiral effective field theory (EFT) have enabled systematic calculations of asymmetric nuclear matter up to approximately twice nuclear saturation density ($2n_0$), a robust framework for quantifying theoretical uncertainties across arbitrary densities and proton fractions has been lacking. Specifically, existing methods often treated pure neutron matter (PNM) and symmetric nuclear matter (SNM) separately, requiring phenomenological interpolations for intermediate proton fractions. Furthermore, constructing consistent EOSs for the neutron star inner crust, which involves non-uniform matter and finite-size effects, requires a uniform matter baseline with well-defined error bands to ensure physical consistency.

Methodology
The authors develop a Bayesian framework to quantify EFT truncation uncertainties for asymmetric nuclear matter. The core of this approach is a two-dimensional Gaussian Process (2D GP) emulator trained on many-body perturbation theory results derived from chiral two- and three-nucleon interactions up to next-to-next-to-next-to-leading order (N3LO).

Key methodological steps include:

1. Order-by-Order Expansion: The energy per particle $E(n, x)$ is expressed as an EFT expansion around a reference energy $E_{\text{ref}}$, with expansion coefficients $c_i(n, x)$ and an expansion parameter $Q(n, x)$ based on the Fermi momentum and the chiral breakdown scale ($\Lambda_b = 600$ MeV).
2. Reference Energy Optimization: To ensure the expansion coefficients exhibit stationarity and naturalness across different proton fractions ($x$), the authors introduce a proton-fraction-dependent reference energy, $E_{\text{ref}}^{BQ+3B}$. This includes a term motivated by three-nucleon (3N) forces, which become significant at higher densities and proton fractions, thereby reducing the magnitude of expansion coefficients and establishing a common scale.
3. 2D Gaussian Process: A 2D GP is constructed where the expansion coefficients are modeled as joint Gaussian distributions dependent on density ($n$) and proton fraction ($x$). The covariance kernel uses a radial basis function (RBF) to account for correlations in both variables. Hyperparameters (marginal variance and correlation lengths) are optimized using an inverse-$\chi^2$ prior for the variance and uniform priors for correlation lengths, trained on data up to N3LO.
4. Model Validation: The GP emulator is validated using model-checking diagnostics, including credibility interval diagnostics and Mahalanobis distance analysis, confirming that the constructed uncertainty bands reasonably capture the order-by-order convergence behavior.
5. Crust Construction: Using the 2D GP emulator for the uniform matter bulk energy, the authors construct the inner crust EOS using a compressible liquid drop model within the Wigner-Seitz approximation. This includes surface and Coulomb corrections. The equilibrium conditions (pressure and chemical potential matching between phases) are solved for the mean EOS and for the upper and lower boundaries of the 68% credibility intervals.

Key Contributions and Results

• Unified Uncertainty Quantification: The paper provides the first 2D GP framework that quantifies EFT truncation uncertainties for asymmetric nuclear matter as a continuous function of density and proton fraction, eliminating the need for phenomenological interpolation between PNM and SNM.
• Thermodynamic Properties: The framework yields the energy per particle, pressure, and chemical potentials for matter in $\beta$-equilibrium up to $2n_0$. At N3LO, the proton fraction in $\beta$-equilibrium is found to increase monotonously with density, reaching approximately 7.5% at $2n_0$.
• Crust-Core Transition: The study constructs consistent EOSs for the neutron star inner crust. The inclusion of surface and Coulomb corrections, combined with the EFT uncertainties from the uniform phase, allows for the determination of the crust-core transition density.
  • The calculated crust-core transition densities range between 0.062 and 0.088 fm$^{-3}$, consistent with other studies.
  • The analysis confirms the existence of a proton drip phase (where protons drip out of nuclei into the uniform phase) in all three considered cases (mean, $+1\sigma$, and $-1\sigma$), though the density range for this phase varies. For the stiffest EOS ($+1\sigma$), the proton drip region is significantly reduced, suggesting it might be skipped in certain scenarios.
• Uncertainty Bands: The 68% credibility bands for the crust EOS are shown to match smoothly to the BPS crust EOS at low densities. The study demonstrates that the EFT uncertainties in the uniform phase propagate consistently to the crust regime, providing a statistically rigorous error estimate for the inner crust.

Significance
The paper claims that its primary significance lies in providing a consistent and uncertainty-quantified EOS for both the uniform outer core and the non-uniform inner crust of neutron stars, derived directly from chiral EFT. By extending the Bayesian GP framework to two dimensions (density and proton fraction), the work enables the calculation of thermodynamic derivatives (such as pressure and chemical potentials) analytically and robustly. The authors emphasize that this approach avoids phenomenological assumptions for the proton fraction dependence and provides a solid theoretical foundation for future astrophysical Bayesian inference, particularly when combining these results with observational data (e.g., from NICER). The work also establishes that proton drip is a robust feature of the inner crust within the current theoretical uncertainties, even when accounting for finite-size effects.