Perturbative calculations of light nuclei up to N$^3$LO in chiral effective field theory

This paper demonstrates that a renormalization-group-invariant power counting scheme in chiral effective field theory, treating subleading two-nucleon interactions perturbatively and calibrating with the triton binding energy, enables robust N$^3$LO predictions for the ground-state energies of $^3$H, $^4$He, and $^6$Li.

The Gist

Imagine you are trying to build a perfect model of a tiny, complex machine using only a few basic Lego bricks. In the world of physics, this machine is an atomic nucleus (the core of an atom), and the bricks are protons and neutrons.

For decades, physicists have been trying to figure out exactly how these bricks snap together. The "blueprint" they use is called Chiral Effective Field Theory (χEFT). Think of this blueprint as a recipe that tells you how to mix the bricks to get the right shape and weight.

However, there's a catch. The recipe has different levels of detail:

• Level 1 (LO): The basic instructions. "Put two bricks together."
• Level 2 (NLO): A bit more detail. "Make sure the bricks are facing the right way."
• Level 3 & 4 (N2LO, N3LO): Extremely fine details. "Add a tiny bit of glue here, adjust the angle by a fraction of a degree."

The problem is that as you add more detail, the math gets incredibly messy and unstable. It's like trying to balance a house of cards in a hurricane. If you try to calculate everything at once, the whole thing collapses.

The Big Breakthrough: "Perturbative" Cooking

This paper, by Oliver Thim and his team at Chalmers University, introduces a clever new way to cook this recipe. Instead of trying to bake the whole cake at once (which causes it to burn), they use a step-by-step approach called perturbation theory.

Here is the analogy:

1. The Base Cake (LO): They bake a simple, stable cake using only the most basic ingredients. This is their "Leading Order" calculation. It's not perfect, but it's solid.
2. The Frosting (Higher Orders): Instead of mixing the frosting into the batter, they calculate how much the frosting would change the weight of the cake if they added it. They do this by looking at the difference between the plain cake and a slightly tweaked version.

This method is called perturbative. It allows them to add the complex, high-level details (up to N3LO, which is the 4th level of precision) without breaking the math.

The "Taste Test" Problem

The team tested this method on three specific "machines":

• Tritium (3H): A nucleus with 3 particles.
• Helium-4 (4He): A nucleus with 4 particles.
• Lithium-6 (6Li): A nucleus with 6 particles.

They found something surprising. When they tried to predict the weight of Helium and Lithium, their results were all over the place. It was like trying to predict the weight of a car by only measuring the engine, but ignoring the wheels.

The Secret Ingredient:
They realized they had to "calibrate" their recipe using the weight of Tritium (3H) first.

• Without Tritium: Their predictions for Helium and Lithium were wild guesses, changing drastically depending on how they adjusted the math.
• With Tritium: Once they forced their recipe to perfectly match the weight of Tritium, the predictions for Helium and Lithium suddenly became incredibly accurate and stable.

It's as if they realized, "Oh, we need to tune our scale using a 3-pound weight before we can trust it for a 4-pound or 6-pound weight."

The "Magic Calculator" Trick

Calculating these tiny changes usually requires knowing every single possible state of the nucleus, which is like trying to count every grain of sand on a beach. For larger nuclei (like Helium and Lithium), this is impossible with current computers.

The authors used a clever mathematical trick called Finite Differences.

• Imagine you want to know how steep a hill is, but you can't measure the slope directly.
• Instead, you take a tiny step forward, measure the height, take a tiny step back, measure the height, and calculate the difference.
• By doing this with their computer code, they could figure out the "slope" (the correction) without needing to map the entire mountain. This allowed them to solve the problem for Helium and Lithium using standard, powerful computers.

Why Does This Matter?

This paper is a huge step forward for two reasons:

1. It works: They proved that this "step-by-step" method can predict the properties of light nuclei with high precision, getting closer to the "gold standard" of physics (Quantum Chromodynamics, or QCD).
2. It's a roadmap: It shows that we don't need super-complex, unstable math to understand the nucleus. We can build it up layer by layer, fixing the recipe as we go.

In a nutshell: The team figured out how to build a stable, high-precision model of atomic nuclei by baking a simple base cake first, tuning their scale with a small 3-piece model, and then using a clever math trick to add the fancy frosting. This brings us one step closer to understanding the fundamental forces that hold our universe together.

Technical Summary

1. Problem Statement

The central challenge in nuclear effective field theory (EFT) is constructing nuclear interactions that are consistent with the symmetries of Quantum Chromodynamics (QCD) and possess predictive power for atomic nuclei.

• Limitations of Standard Power Counting: The standard Weinberg Power Counting (WPC) has been successful but violates Renormalization Group (RG) invariance at the nucleon-nucleon (NN) level. This is due to the singular and attractive nature of one-pion exchange (OPE) potentials, which generate divergences when treated non-perturbatively without sufficient counterterms.
• The Perturbative Alternative: Modified power counting (PC) schemes, such as the one proposed by Long and Yang, treat subleading orders perturbatively to satisfy RG invariance. While these schemes have been tested for NN scattering and light nuclei up to Next-to-Leading Order (NLO), their ability to provide accurate, predictive descriptions of nuclear structure up to Next-to-Next-to-Next-to-Leading Order (N3LO) had not been fully established.
• Computational Hurdle: A major barrier to applying perturbative PC to many-body systems (beyond the NN system) is the computational difficulty of calculating higher-order corrections, which traditionally require knowledge of the full energy spectrum (Rayleigh-Schrödinger perturbation theory), a task often impossible for complex nuclei using current many-body methods.

2. Methodology

The authors address these challenges by combining a modified RG-invariant power counting scheme with advanced computational techniques within the No-Core Shell Model (NCSM).

A. Theoretical Framework

• Power Counting: They employ the Long and Yang PC scheme up to N3LO.
  • Leading Order (LO): Treated non-perturbatively. Includes OPE for specific channels ($^1S_0$, $^3P_0$, $^1P_1$, $^3P_1$, $^3S_1$-$^3D_1$, $^3P_2$-$^3F_2$) and contact interactions.
  • Subleading Orders (NLO, N2LO, N3LO): Treated perturbatively. Includes higher-order contacts and corrections to OPE.
  • Cutoff: A relatively low momentum space cutoff ($\Lambda \approx 500$ MeV) is used to avoid "exceptional cutoff" issues and ensure convergence.
• Calibration: Low-Energy Constants (LECs) are calibrated to neutron-proton phase shifts and binding energies of light nuclei ($^2$H, $^3$H).

B. Computational Techniques

• Rayleigh-Schrödinger Perturbation Theory (RSPT): Used for $^3$H (Triton) using the Jacobi-coordinate NCSM code (py-ncsm). This allows for the explicit calculation of energy corrections ($E^{(1)}, E^{(2)}, E^{(3)}$) using the full LO spectrum.
• Finite-Difference (FD) Method: To overcome the inability to compute full spectra for heavier nuclei ($^4$He, $^6$Li) using the M-scheme NCSM code (pANTOINE), the authors utilize numerical derivatives.
  • They define a Hamiltonian $H(x) = H^{(0)} + \sum x_\nu V^{(\nu)}$.
  • Perturbative corrections are extracted as derivatives of the ground-state energy $E_0(x)$ with respect to the perturbation parameters $x$ at $x=0$.
  • This avoids the need for the full spectrum, requiring only iterative solutions (Lanczos diagonalization) for specific states.
  • The method uses finite-difference stencils (up to 4th order) to achieve high numerical accuracy, validated against exact RSPT results in $^3$H.

3. Key Contributions

1. Extension to N3LO: The first application of RG-invariant perturbative chiral EFT to light nuclei ($^3$H, $^4$He, $^6$Li) up to N3LO.
2. Development of the FD Method: Demonstration that numerical derivatives of non-perturbative eigenvalues can accurately reproduce perturbative corrections for many-body systems, bypassing the need for full spectral resolution.
3. Calibration Sensitivity Analysis: A systematic study of how different calibration strategies for Low-Energy Constants (LECs) affect the predictive power for $A > 2$ nuclei.
4. Validation of RG-Invariant PC: Evidence that this specific power counting scheme can yield quantitatively accurate predictions for bulk nuclear observables.

4. Key Results

• Validation of the FD Method: In $^3$H, the Finite-Difference method matched exact Rayleigh-Schrödinger results with high precision (relative error $\approx 10^{-4}\%$), confirming its reliability for $^4$He and $^6$Li.
• Impact of Calibration on Predictions:
  • Interaction (A): Calibrated solely on NN phase shifts. Failed to reproduce $^4$He and $^6$Li binding energies accurately, showing large subleading corrections and strong cutoff dependence.
  • Interaction (B): Calibrated to include the $^3$H binding energy at N2LO/N3LO and the deuteron binding energy at N3LO. This resulted in significant improvements:
    • Accurate reproduction of $^4$He and $^6$Li ground-state energies.
    • Marked reduction in cutoff dependence.
    • Better order-by-order convergence.
  • Interaction (C): Included the deuteron binding energy at LO. While it improved $^3$H, it led to an overbinding of $^4$He at N3LO, suggesting that forcing strong attraction at LO is not always advantageous.
• Convergence:
  • For $^4$He, the ground-state energy is well-converged up to N3LO.
  • For $^6$Li, convergence is slower but achieves a precision of $\approx 1$ MeV at N3LO.
  • The estimated EFT truncation error at N3LO is approximately 5%.
• Caveats:
  • Including $^3$H in the calibration reduced cutoff dependence in energies but increased it in the $^3S_1$-$^3D_1$ mixing angle. This is likely due to the absence of three-nucleon forces (3NF), which are expected to enter at this order.
  • Some interactions predicted $^6$Li to be unbound (above the $^4$He + $^2$H threshold), a known issue in LO WPC that persists here.

5. Significance

• Bridging QCD and Nuclear Structure: The work demonstrates that nuclear interactions constructed with RG-invariant power counting can successfully predict properties of light nuclei, bringing ab initio nuclear structure closer to a foundation in QCD.
• Computational Scalability: The successful application of the Finite-Difference method to $^4$He and $^6$Li proves that perturbative corrections can be computed for complex nuclei without the prohibitive cost of full spectral calculations. This opens the door for applying these methods to heavier nuclei and other many-body frameworks.
• Predictive Power: The study establishes that including few-body binding energies (specifically $^3$H) in the calibration of NN interactions is essential for robust predictions of larger nuclei.
• Future Directions: The authors identify the inclusion of three-nucleon forces (3NF), Coulomb interactions, and isospin-breaking as necessary next steps to resolve remaining cutoff dependencies and improve accuracy. They also highlight the potential of automatic differentiation and emulators to further systematize these calculations.