Cutoff effects in Hartree-Fock calculations at leading order of chiral effective field theory

The study demonstrates that using leading-order chiral effective field theory potentials with regularization cutoffs exceeding the nucleon mass prevents the Hartree-Fock approximation from yielding self-consistent, bound mean-field solutions for nuclei like $^{16}$O due to spurious deeply bound states, thereby necessitating that corrections for these artifacts be incorporated into residual interactions for beyond-mean-field treatments.

The Gist

Imagine the atomic nucleus as a tiny, bustling city made of protons and neutrons (nucleons) living together in a very small space. To understand how this city holds together, physicists use a set of mathematical rules called "Chiral Effective Field Theory." Think of this theory as a rulebook for how these tiny citizens interact.

However, there's a catch: when you try to do the math for these interactions, the numbers can blow up to infinity, making the calculations impossible. To fix this, scientists use a "cutoff." Imagine this cutoff as a speed limit or a blur filter on a camera. It says, "We will ignore any interactions that happen at speeds or distances smaller than this specific limit." This keeps the math manageable.

The paper by Sánchez Sánchez, Duc, and Bonneau investigates what happens when you set this "speed limit" (the cutoff) very high—higher than the mass of the nucleons themselves. They specifically looked at the Oxygen-16 nucleus (a city with 16 citizens) using a method called the Hartree-Fock approximation.

The Problem: Ghosts in the Machine

When the scientists set this speed limit very high, they ran into a strange problem. The math started creating "spurious deeply bound states."

Think of these as ghosts. In the real world, these specific combinations of particles shouldn't exist; they are mathematical errors caused by the high speed limit. But in the calculation, the math says, "Hey, look! There's a super-tightly bound particle here!" These ghosts are so attractive that they mess up the entire calculation, making the nucleus look like it's collapsing or behaving in ways that don't match reality.

The Experiment: Cleaning Up the City

The researchers tried to fix this by "kicking out" these ghosts. They added a mathematical tool (a projector) to push these fake, ghostly states away, similar to how you might use a magnet to push away a piece of metal that doesn't belong.

They tested this in two ways:

1. Low Speed Limit (500 MeV): The math was calm. No ghosts appeared. The calculation worked perfectly, and the nucleus looked stable.
2. High Speed Limit (1500 MeV): Ghosts appeared. When the scientists tried to kick them out, they found a major issue.

The Big Discovery: The "Goldilocks" Failure

Here is the core finding, explained simply:

When they tried to remove the ghosts in the high-speed limit scenario, the nucleus refused to stay together.

• The Analogy: Imagine trying to build a house of cards. If you blow too hard (high cutoff), the cards fly apart. You try to tape them down (remove the ghosts), but the tape is so strong that it actually pushes the cards apart even more.
• The Result: The "Hartree-Fock" method, which is supposed to be the foundation for understanding the nucleus, broke down. It couldn't produce a stable nucleus that was free of these ghost errors.

The paper concludes that if you use this specific high-speed rulebook (the Nogga–Timmermans–van Kolck power counting) with a high cutoff, the Hartree-Fock method cannot give you a clean, stable starting point. It's like trying to build a skyscraper on a foundation that keeps sinking whenever you try to fix the cracks.

What Can Be Done?

The authors suggest a compromise. You can't fix everything at the basic level (the Hartree-Fock level).

• You can fix some of the ghost problems to get a stable nucleus, but you have to leave the rest of the mess for more advanced, complex calculations to handle later.
• Specifically, they found that while they could fix the ghosts in some channels (like the 3D2 channel) without breaking the nucleus, trying to fix the ghosts in the most important channel (3S1) completely destroyed the nucleus's stability.

Summary

In short, this paper is a warning label for nuclear physicists. It says: "If you use a very high speed limit in your calculations, the simple 'average' method (Hartree-Fock) will fail to give you a stable nucleus because the mathematical 'ghosts' are too strong to remove without breaking the whole system."

To get accurate results with these high limits, you have to accept that the simple method isn't enough; you need to use much more complex methods to clean up the remaining mess.

Technical Summary

1. Problem Statement

The paper investigates the behavior of nuclear bulk properties when using chiral Effective Field Theory (EFT) potentials within the Nogga–Timmermans–van Kolck (NTvK) power counting scheme, specifically at Leading Order (LO).

• Context: Standard ab initio calculations often use Weinberg power counting. However, the NTvK scheme advocates for non-perturbative renormalization of attractive, singular partial waves (where the one-pion exchange tensor force is attractive) to ensure renormalization-group invariance. This requires large regularization cutoffs ($\Lambda$), often exceeding the nucleon mass and the chiral symmetry-breaking scale ($\sim 1$ GeV).
• The Issue: Using large cutoffs in the NTvK scheme leads to the appearance of spurious deeply bound states (unphysical bound states) in the two-nucleon system. While these can be corrected in few-body systems (like $A=3$) or diagonalization-based many-body methods (like No-Core Shell Model), their impact on Hartree-Fock (HF) mean-field approximations is poorly understood.
• Objective: The authors aim to determine if a self-consistent HF mean field can be constructed using LO NTvK potentials with large cutoffs ($\Lambda > m_N$) that is free of spurious bound-state effects and suitable as a reference state for beyond-mean-field methods (e.g., Coupled Cluster, In-Medium SRG).

2. Methodology

The authors perform Hartree-Fock calculations for the $^{16}$O nucleus using a specific framework designed to handle the complexities of the NTvK potential.

• Potential Model:
  • The LO potential consists of the One-Pion Exchange (OPE) and a contact term.
  • Power Counting: They adopt the NTvK scheme, promoting counter-terms in singular attractive channels ($^3S_1$, $^3P_0$, $^3P_2$, $^3D_2$) to LO for non-perturbative renormalization.
  • Spurious State Correction: To remove unphysical deeply bound states arising from large $\Lambda$, they apply an energy-shift projector method: $\hat{V}' = \hat{V} + \sum E^{(i)}_{\text{shift}} |\chi_i\rangle\langle\chi_i|$, where $|\chi_i\rangle$ are the spurious bound states and $E_{\text{shift}}$ are large positive energy shifts.
• Numerical Framework:
  • Basis: A confined plane-wave basis in a cubic box of edge length $L$. This avoids the need for center-of-mass to laboratory frame transformations required in harmonic oscillator bases.
  • Symmetry: The basis respects the octahedral symmetry of the cube ($O_{2h}^D$), allowing for the description of triaxial deformations and time-reversal symmetry breaking.
  • Convergence: They systematically vary the box size $L$ (infrared cutoff) and the single-particle momentum cutoff $k_{\text{max}}$ (ultraviolet cutoff) to ensure convergence.
• Cutoff Values: Calculations are performed for regularization cutoffs $\Lambda = 500, 1000,$ and $1500$ MeV.

3. Key Contributions and Findings

A. Cutoff Dependence and Spurious States

• $\Lambda = 500$ MeV: No spurious bound states appear. The HF solution converges well with respect to UV and IR cutoffs, yielding a bound ground state with a charge radius close to experimental values.
• $\Lambda = 1000$ MeV: A spurious bound state appears in the $^3P_0$ channel. The uncorrected HF solution yields a positive (unbound) ground-state energy and an oscillating, unphysically large charge radius. This is attributed to a strongly repulsive contact potential in the $^3S_1$ channel.
• $\Lambda = 1500$ MeV: Spurious bound states appear in the coupled $(^3S_1, ^3D_1)$, $^3P_0$, and $^3D_2$ channels.

B. The Failure of Full Correction in Mean-Field

The core finding is the extreme sensitivity of the HF solution to the energy shifts ($E_{\text{shift}}$) used to remove spurious states:

• Repulsion vs. Binding: Removing spurious states requires adding a repulsive projector term. In the HF approximation, this repulsion is distributed over the mean field.
• The $^3S_1$ and $^3P_0$ Channels: Even modest energy shifts in these channels (which dominate the interaction) cause the HF ground state to become unbound. The repulsive mean field pushes single-particle states above the Fermi level, causing the iterative solution to diverge.
• The $^3D_2$ Channel: In contrast, the $^3D_2$ channel is less sensitive. It is possible to apply "full" corrections (using recommended shifts of 10–15 GeV) in this channel while maintaining a bound HF solution.

C. Implications for Beyond-Mean-Field Methods

• No Valid Reference State: The authors conclude that it is impossible to generate a self-consistent HF mean field free of spurious bound-state effects using LO NTvK potentials with large cutoffs.
• Breakdown of Variational Principle: Unlike diagonalization methods (e.g., NCSM) where the energy becomes independent of large $E_{\text{shift}}$ once spurious states are decoupled, the HF Slater determinant breaks translation symmetry and lacks the correlations necessary to absorb the strong repulsive shifts.
• Proposed Solution: A "partial" correction strategy is proposed. One can incorporate a partial correction for spurious states (e.g., in $^3D_2$) into the mean field to achieve UV convergence, but the remaining correction must be treated as part of the residual interaction in a beyond-mean-field calculation.

4. Significance

• Limitation of Mean-Field Approaches: The study demonstrates a fundamental limitation of the Hartree-Fock approximation when applied to chiral EFT potentials renormalized non-perturbatively at large cutoffs. It cannot serve as a standalone reference state for these specific power counting schemes.
• Guidance for Future Calculations: The results suggest that for NTvK power counting with large cutoffs, the treatment of spurious bound states cannot be fully relegated to the mean-field level. Future many-body calculations (Coupled Cluster, IM-SRG) must handle the residual corrections to the potential explicitly.
• Validation of Power Counting Debates: The work highlights the practical difficulties of the NTvK scheme in many-body systems, reinforcing the debate regarding the necessity of large cutoffs versus the stability of the resulting mean fields.

Conclusion

The paper establishes that while the NTvK power counting ensures renormalization-group invariance in the two-nucleon sector, its application in the Hartree-Fock approximation for nuclei (specifically $^{16}$O) fails when using large cutoffs ($\Lambda \ge 1000$ MeV). The necessary corrections for spurious deeply bound states introduce repulsion that destroys the binding of the mean-field solution, particularly due to the sensitivity of the $^3S_1$ channel. Consequently, a "full" correction cannot be achieved at the mean-field level; a hybrid approach where partial corrections are applied to the mean field and the remainder is treated in the residual interaction is required.