Bayesian Inference of the Landau Parameter $G'_0$ from Joint Gamow-Teller Measurements

Using Bayesian inference and a self-consistent Skyrme Random Phase Approximation model constrained by joint Gamow-Teller resonance measurements in $^{208}\mathrm{Pb}$, $^{132}\mathrm{Sn}$, and $^{90}\mathrm{Zr}$, this study reports the first quantified determination of the Landau-Migdal parameter $G'_0$ as $0.48\pm0.034$, offering new insights for constructing energy density functionals that accurately describe spin-isospin interactions in dense nuclear matter.

The Gist

Imagine the atomic nucleus not as a solid ball, but as a bustling, chaotic dance floor filled with tiny dancers called protons and neutrons. These dancers are constantly interacting, pushing and pulling on each other.

In the world of nuclear physics, scientists have a specific rulebook to describe how these dancers interact when they spin and swap identities (a proton turning into a neutron, or vice versa). This rulebook is called the Landau Parameter $G'_0$. Think of $G'_0$ as the "tension setting" on the dance floor. If the tension is high, the dancers push each other away strongly when they try to spin; if it's low, they are more relaxed.

Knowing the exact value of this tension is crucial. It helps us understand:

• How unstable atoms decay (beta decay).
• What happens inside exploding stars (supernovae).
• How neutron stars (the densest objects in the universe) behave.

The Problem: Guessing the Tension

For decades, scientists tried to figure out the value of this "tension setting" ($G'_0$) by looking at a specific dance move called the Gamow-Teller Resonance (GTR). This is like a synchronized spin that happens when you hit a nucleus with a particle.

However, previous attempts were like trying to guess the weight of a hidden object by looking at a blurry photo.

1. They used "schematic" models: They made simplified drawings of the dance floor that didn't match reality perfectly.
2. They made assumptions: They guessed how heavy the dancers felt (effective mass) without measuring it directly.
3. The results were messy: Different scientists got different answers, ranging from 1.0 to 1.5, with no clear way to say which one was right or how wrong they might be.

The New Approach: A Bayesian Detective Story

In this paper, the authors (Lin, Colò, Steiner, and Stinson) decided to solve this mystery using a modern detective tool called Bayesian Inference.

The Analogy: The Master Chef and the Taste Test
Imagine you are trying to recreate a famous soup recipe, but you don't have the exact list of ingredients. You only know how the soup tastes (the experimental data).

• Old Method: You guess a few ingredients, taste the soup, and say, "It needs more salt." You guess again. It's a hit-or-miss process.
• The New Bayesian Method: You create a massive digital simulation of every possible soup recipe (thousands of variations of the nuclear model). You then "taste" every single one of them against the real experimental data (the GTR measurements from three specific atoms: Lead-208, Tin-132, and Zirconium-90).

Instead of picking just one recipe, the computer tells you: "Out of all the possible recipes, 95% of the ones that taste right have a specific amount of salt." It gives you a probability map rather than a single guess.

What They Did

1. Built a Self-Consistent Kitchen: They used a sophisticated model (Skyrme Random Phase Approximation) where the "ingredients" (nuclear forces) and the "cooking process" (how the nucleus moves) are perfectly linked. No more guessing parts of the recipe separately.
2. The Joint Taste Test: They didn't just look at one atom. They looked at three different "dance floors" (Lead, Tin, Zirconium) simultaneously. If a recipe works for all three, it's a winner.
3. Accounting for Uncertainty: They admitted, "Our model isn't perfect." They added "fudge factors" (called Intrinsic Scattering parameters) to account for the fact that our current understanding of nuclear physics might have small blind spots.

The Big Discovery

After running millions of simulations, they found the answer.

• The Old Guess: Scientists used to think the tension ($G'_0$) was around 1.0 to 1.5.
• The New Result: The "tension" is actually much lower: 0.48 (with a very small margin of error, $\pm 0.034$).

Why the difference?
The authors suggest that previous models were like using a heavy, stiff dance floor, while the real nucleus is more like a springy, flexible one. The "effective mass" of the nucleons (how heavy they feel inside the nucleus) is lighter than previously thought. Because the dancers are lighter, they don't need as much "tension" to push them apart during the spin.

Why Does This Matter?

This isn't just about numbers on a page.

• Better Star Models: If we know the correct "tension," we can simulate supernova explosions and neutron star mergers much more accurately. This helps us understand how heavy elements (like gold and uranium) are forged in the universe.
• New Recipes: This result acts as a strict guideline for building the next generation of nuclear models. Any new theory must now produce a tension of roughly 0.48 to be considered valid.

In a Nutshell

The authors used a powerful statistical method (Bayesian inference) to combine data from three different atoms and a highly accurate computer model. They discovered that the "pushiness" between protons and neutrons is about half of what scientists previously believed. This corrects our understanding of the atomic dance floor and helps us better predict the violent, beautiful events happening in the deepest corners of the cosmos.

Technical Summary

1. Problem Statement

The Landau-Migdal parameter $G'_0$ characterizes the spin-isospin dependent part of the nucleon-nucleon interaction. It is a critical quantity for understanding:

• The Gamow-Teller Resonance (GTR) in finite nuclei.
• Beta and double-beta decay rates.
• Neutrino-nucleon absorption rates in core-collapse supernovae (CCSNe) and binary neutron star (BNS) mergers.
• The critical density for pion condensation in neutron stars.

Historical Challenges:

• Inconsistency: Previous extractions of $G'_0$ yielded widely varying values (e.g., ranging from $\approx 1.0$ to $1.5$) depending on the nucleus and the model used.
• Lack of Self-Consistency: Traditional methods relied on schematic $\pi$- (or $\pi+\rho$-) exchange models combined with empirical single-particle states (Woods-Saxon) and simplified linear response theory. These approaches lacked a self-consistent link between the ground state and excited states.
• Uncertainty Quantification: Previous studies often failed to provide robust statistical uncertainties or account for model systematic errors.
• Effective Mass Ambiguity: Converting phenomenological parameters ($g'$) to the Landau parameter ($G'_0$) requires knowledge of the nucleon effective mass ($m^*$), which was often estimated roughly, introducing uncontrolled uncertainties.

2. Methodology

The authors employ a Bayesian inference framework to extract $G'_0$ with quantified uncertainties, utilizing a self-consistent theoretical model.

• Theoretical Framework:

  • Skyrme Energy Density Functionals (EDFs): Used to describe the mean field.
  • Random Phase Approximation (RPA): A self-consistent Quasiparticle RPA (QRPA) scheme is used to calculate GTR features based on the Hartree-Fock (HF) ground states.
  • Residual Interaction: The spin-isospin dependent part of the interaction is dominated by the density-dependent Landau parameter $G'_0(\rho)$.

• Bayesian Inference Setup:

  • Priors: Uniform priors were defined for Skyrme parameters ($t_0, t_1, t_2, t_3, x_0, x_1, x_2, x_3, \alpha, b_4, b'_4$) based on the range of 11 widely used Skyrme models.
  • Likelihood & Constraints: The model parameters were constrained by a joint dataset including:
    • GTR Observables: Peak energy ($E_{GT}$) and the percentage of the Ikeda sum rule ($m_{0\%}$) for $^{208}\text{Pb}$, $^{132}\text{Sn}$, and $^{90}\text{Zr}$.
    • Ground State Properties: Binding energies and charge radii for $^{208}\text{Pb}$, $^{132}\text{Sn}$, $^{90}\text{Zr}$, and $^{48}\text{Ca}$.
    • Spin-Orbit Splittings: 9 specific splitting values across the same nuclei.
  • Systematic Uncertainties: The authors introduced Intrinsic Scattering (IS) parameters to account for model deficiencies (systematic uncertainties). These parameters act as additional degrees of freedom to inflate credible regions when the model cannot perfectly describe the data.

• Sampling: Markov Chain Monte Carlo (MCMC) methods were used to sample the posterior distribution of the Skyrme parameters, from which the posterior distribution of $G'_0$ was derived.

3. Key Contributions

• First Self-Consistent Extraction: This is the first study to extract $G'_0$ using a fully self-consistent Skyrme-RPA framework constrained by joint experimental data, eliminating the need for schematic assumptions or empirical single-particle inputs.
• Quantified Uncertainty: The study provides a rigorous statistical uncertainty for $G'_0$ ($0.48 \pm 0.034$), moving beyond point estimates.
• Resolution of Tension: The work highlights and investigates the tension between values extracted from phenomenological pion-exchange models (typically $G'_0 > 1.0$) and those predicted by standard Skyrme EDFs (typically $G'_0 < 0.6$).
• Correlation Analysis: The study systematically analyzes the correlations between $G'_0$, the effective mass ($m^*$), and GTR observables ($E_{GT}$, $m_{0\%}$).

4. Results

• Extracted Value: The Bayesian posterior distribution yields a value of $G'_0 = 0.48 \pm 0.034$ at nuclear saturation density.
• Comparison with Models:
  • This value is consistent with a few existing Skyrme models that specifically consider spin-isospin observables (e.g., SGII, SAMI).
  • It is significantly smaller than values derived from traditional phenomenological pion-exchange models (which often report $G'_0 \approx 1.0 - 1.5$).
• Role of Effective Mass ($m^*$):
  • The study finds a strong correlation between $G'_0$ and the effective mass ratio $m^*/m$.
  • The Bayesian inference suggests a lower effective mass ($m^*/m \approx 0.67$) compared to the values ($0.8 - 1.0$) often assumed in phenomenological models.
  • Explanation of Discrepancy: The authors argue that the higher $G'_0$ values in previous works are likely artifacts of using simplified assumptions (e.g., schematic RPA, constant attenuation factors) and overestimating the effective mass. In self-consistent models, a higher $G'_0$ would lead to an overestimation of the Ikeda sum rule fraction ($m_{0\%}$).
• Predictive Power: The "GTAll" posterior distribution successfully reproduces the experimental GTR peak energies and sum rule fractions for the constrained nuclei while maintaining consistency with ground-state properties.

5. Significance

• Astrophysical Implications: The extracted $G'_0$ and the associated Skyrme models provide a more reliable basis for calculating neutrino-nucleon interaction rates in dense matter, which are crucial for simulating core-collapse supernovae and neutron star mergers.
• Model Construction: The results guide the construction of new Energy Density Functionals (EDFs). Future EDFs aiming to describe dense matter in the spin-isospin channel should target $G'_0 \approx 0.48$ rather than the historically higher values.
• Methodological Advancement: The paper demonstrates the power of Bayesian inference combined with intrinsic scattering parameters to handle model systematic uncertainties in nuclear physics, setting a precedent for future studies of nuclear forces and astrophysical phenomena.
• Future Directions: The authors note that while the current extraction is robust within the RPA framework, future work should explore extensions to second RPA or particle-vibration coupling (PVC) models to further reduce model dependency.