Multi-neutron correlations in light nuclei via ab-initio lattice simulations

Using ab initio nuclear lattice effective field theory with Bayesian uncertainty quantification, this study reveals that the ground states of ${}^6$H and ${}^7$H are characterized by valence neutrons forming compact dineutron pairs that predominantly organize into symmetric dineutron-dineutron configurations, thereby providing crucial structural insights into multi-neutron correlations and the nature of four-neutron clusters.

The Gist

Imagine the atomic nucleus as a tiny, bustling city. Usually, this city is a balanced mix of "proton citizens" and "neutron residents." But in some very rare, unstable atoms (like the ones studied in this paper), the city is overflowing with neutrons, creating a chaotic, neutron-rich neighborhood.

For decades, physicists have been trying to answer a big question: Do these extra neutrons ever hang out together in tight little groups? Specifically, can four neutrons stick together to form a temporary "super-neutron" cluster (called a tetraneutron), or do they just float around loosely?

This paper is like a high-tech, super-powerful microscope that allows scientists to zoom in on these neutron-rich atoms and watch exactly how the neutrons behave. Here is the breakdown of their findings using simple analogies:

1. The Problem: A Noisy Crowd

Scientists have been trying to spot these "four-neutron groups" in experiments for a long time. It's like trying to find a specific group of four friends holding hands in a massive, noisy stadium. The data is messy, and different experiments give different answers. Some say the groups exist; others say they don't.

2. The Tool: A Digital Time Machine

The authors used a supercomputer method called Ab Initio Lattice Effective Field Theory.

• The Analogy: Imagine you want to know how a crowd of people will move in a room, but you can't watch them in real life because they move too fast. Instead, you build a perfect digital simulation of the room and the people. You run the simulation forward in "digital time" until the crowd settles into its most natural, comfortable arrangement.
• The Innovation: They didn't just guess the rules of how neutrons interact; they used a massive library of 282 different "rulebooks" (mathematical models based on real-world data) and ran the simulation for all of them. Then, they used a statistical method (Bayesian analysis) to see which rulebooks fit the real-world data best. This gave them a very precise, uncertainty-checked answer.

3. The Discovery: Two Types of "Neutron Parties"

When they looked at the atoms Hydrogen-7 (a proton core with 6 neutrons) and Helium-8 (a core with 4 neutrons), they found the neutrons weren't just floating randomly. They formed two distinct types of "parties":

Party A: The "Back-to-Back" Duet (The Dominant Style)

• What it is: The four extra neutrons form two tight pairs (called dineutrons). Imagine two couples dancing.
• The Arrangement: These two couples stand on opposite sides of the atom's core, facing away from each other. It's like two pairs of dancers standing back-to-back, holding hands within their own pair, but separated from the other pair.
• The Frequency: This is the most common arrangement, happening about 95% of the time.
• Why it matters: This explains why experiments have been confused. When scientists look for a "four-neutron blob," they are often actually seeing these two separate pairs standing apart. It looks like a four-neutron system, but it's really two two-neutron systems playing nice together.

Party B: The "Tight Huddle" (The Rare Style)

• What it is: This is the true "tetraneutron" everyone was looking for. All four neutrons huddle together in one tight, compact ball on one side of the atom.
• The Arrangement: They twist and turn around each other in a complex, non-flat shape (like a twisted knot).
• The Frequency: This happens only about 5% of the time.
• The Twist: Even though it's rare, it proves that a true four-neutron cluster can exist, but it's very hard to catch because it's so fleeting and small compared to the "back-to-back" style.

4. The Big Picture: Why This Matters

• The "Ghost" in the Machine: The paper suggests that when experimentalists think they see a "four-neutron particle," they might actually be seeing the "Back-to-Back" style (the two pairs). This makes it very hard to isolate the "Tight Huddle" (the real four-neutron cluster).
• The Energy Check: They calculated the energy of these atoms and found that the Hydrogen-7 atom is just barely holding together. It's like a house of cards that is about to fall. Because it's so unstable, it prefers to spit out neutrons in groups rather than one by one.
• The Structure: Inside these atoms, the neutrons aren't just a gas; they are organized. They form little "sub-clubs" (dineutrons) that then arrange themselves in specific geometric patterns.

Summary in One Sentence

Using a super-precise digital simulation, scientists discovered that in these neutron-rich atoms, the neutrons mostly form two separate pairs standing on opposite sides of the core, while a true, tight cluster of four neutrons exists only rarely, explaining why spotting them in real experiments has been so difficult.

The Takeaway: Nature loves symmetry. Even in the chaotic world of extra neutrons, they prefer to pair up and stand apart, rather than huddling all together in a single tight ball.

Technical Summary

1. Problem Statement

The nature of multi-neutron systems, particularly candidate three-neutron ($3n$) and four-neutron ($4n$) clusters (tetraneutrons), has been a subject of controversy for over half a century. While recent experimental efforts using quasi-free knockout reactions on neutron-rich light nuclei (like $^8$He and $^7$H) have hinted at preformed multi-neutron configurations, the ground-state energies of the extreme hydrogen isotopes $^6$H and $^7$H remain poorly constrained. There are significant discrepancies between experimental analyses and theoretical predictions. Furthermore, directly calculating four-body (and higher) correlation functions in nuclear systems is computationally prohibitive due to the combinatorial explosion of dimensionality in reduced density matrices. The core challenge is to determine the internal structure of these systems: do they contain compact tetraneutron clusters, or are the correlations dominated by other geometries (e.g., dineutron pairs)?

2. Methodology

The authors employ Ab Initio Nuclear Lattice Effective Field Theory (NLEFT) combined with Auxiliary-Field Monte Carlo (AFMC) sampling to address these challenges.

• Chiral Interactions & Uncertainty Quantification:

  • The study utilizes an ensemble of 282 non-implausible chiral interactions up to Next-to-Next-to-Next-to-Leading Order (N3LO).
  • These interactions were generated via history matching, systematically incorporating uncertainties from nucleon-nucleon (NN) scattering data, chiral hierarchy truncation, and binding energies of selected nuclei.
  • A Bayesian analysis is performed using Posterior Predictive Distributions (PPD) to quantify uncertainties in observables, moving beyond single "best-fit" interactions to provide statistically robust error bars.

• Lattice Setup:

  • Simulations are performed on a periodic cubic lattice with side length $L=12$ and spatial spacing $a \approx 1.32$ fm (momentum cutoff $\Lambda \approx 471$ MeV).
  • Euclidean time evolution is used with a temporal spacing $a_t = 0.001$ MeV$^{-1}$.
  • Finite-volume effects were assessed and found to be subdominant to statistical uncertainties.

• Observables and Correlation Functions:

  • Pinhole Method: To probe intrinsic structure, the authors use the pinhole method to sample the spatial coordinates, spins, and isospins of all nucleons according to the quantum amplitudes of the $A$-body density.
  • Intrinsic Density Alignment: Candidate clusters (triton in $^7$H, $\alpha$-particle in $^8$He) are identified based on spin-isospin quantum numbers and spatial compactness. The density profiles are realigned relative to these cluster centers to reveal intrinsic structures.
  • Correlation Functions:
    • Two-body: Defined as $\rho_2(r, r', \theta)$ to map spatial distributions and opening angles.
    • Reduced Four-body: A novel reduced four-body correlation function $\rho_4(\theta_1, \phi_1, \theta_2, \phi_2; \Theta, \zeta)$ is constructed. This projects the four-body density onto configurations defined by two disjoint nucleon pairs, characterized by intra-pair angles ($\theta, \phi$), inter-pair opening angle ($\Theta$), and torsion angle ($\zeta$).

3. Key Contributions

1. Bayesian Uncertainty Quantification: The first application of a rigorous Bayesian framework with a large ensemble of chiral forces to the ground-state energies of $^6$H and $^7$H, resolving previous theoretical ambiguities.
2. Direct Four-Body Correlation Calculation: Overcoming the "curse of dimensionality" by using AFMC sampling to compute reduced four-body correlation functions without explicitly constructing the full $A$-body density matrix.
3. Geometric Characterization of Multi-Neutron Systems: Providing a microscopic, geometric description of the four-neutron subsystems in $^7$H and $^8$He, distinguishing between extended symmetric configurations and compact tetraneutron-like structures.

4. Key Results

A. Ground-State Energies and Decay Channels

• Energies: The calculated ground-state energies are $E(^6\text{H}) = -4.94^{+0.14}_{-0.14}$ MeV and $E(^7\text{H}) = -5.29^{+0.16}_{-0.16}$ MeV. These values fall within experimental ranges but clarify the energy hierarchy.
• Separation Energy: The single-neutron separation energy for $^7$H is found to be $S_n(^7\text{H}) = 0.35^{+0.32}_{-0.32}$ MeV.
• Implication: The positive (though small) $S_n$ kinematically disfavors sequential decay via $^6$H + $n$. This suggests that multi-neutron emission channels (direct emission of multiple neutrons) are comparatively more relevant, supporting the existence of correlated multi-neutron dynamics.

B. Intrinsic Structure and Clustering

• Core-Valence Picture: Intrinsic densities reveal a compact triton cluster in $^7$H and an $\alpha$-like cluster in $^8$He, surrounded by four diffuse valence neutrons.
• Dineutron Formation: Two-neutron correlation functions show that valence neutrons form compact dineutron pairs in the surface region.

C. Four-Neutron Geometries

The analysis of the reduced four-body correlation function reveals two dominant geometric configurations for the four valence neutrons:

1. Extended Symmetric 2n-2n Configuration (Dominant):

   • Probability: Accounts for approximately 95% of the sampled configurations in $^7$H (and a similar hierarchy in $^8$He).
   • Geometry: Two dineutron pairs are located on opposite sides of the core. The inter-pair opening angle $\Theta$ peaks around $141^\circ$ (near $180^\circ$, i.e., back-to-back).
   • Structure: The pairs exhibit finite torsion ($\zeta \approx 48^\circ$) and dihedral angles, indicating a twisted, non-coplanar arrangement that is approximately spatially symmetric about the center of mass.

2. Compact Tetraneutron-like Configuration (Minor):

   • Probability: Accounts for approximately 5% of the configurations.
   • Geometry: All four neutrons reside on the same side of the core. The inter-pair opening angle $\Theta$ is small ($60^\circ$–$90^\circ$).
   • Structure: Characterized by large torsion and dihedral angles ($60^\circ$–$90^\circ$), forming a strongly non-coplanar, interlaced, and compact geometry. This represents a structural precursor to a genuine tetraneutron cluster with enhanced four-body correlations.

• Comparison: The extended 2n-2n geometry is significantly more prominent in $^7$H than in $^8$He, where the symmetric configuration is somewhat suppressed, likely due to core-induced medium effects from the different proton backgrounds.

5. Significance

• Resolving the Tetraneutron Debate: The results suggest that experimental signatures of "tetraneutrons" in light nuclei may be heavily contaminated by the dominant symmetric 2n-2n modes. Isolating a genuinely compact, four-body tetraneutron signal is difficult because the compact component constitutes only a small fraction (~5%) of the total wavefunction.
• Theoretical Advancement: The work demonstrates that modern ab initio methods can successfully access complex many-body correlations (beyond two-body) in light nuclei, providing a microscopic foundation for understanding emergent clustering in neutron-rich matter.
• Experimental Guidance: The findings imply that future experimental searches for tetraneutron signatures must account for the strong contribution of extended, symmetric dineutron-dineutron configurations to avoid misinterpretation of data. The positive separation energy of $^7$H also guides the interpretation of decay channels in quasi-free knockout experiments.

In summary, this paper provides a definitive ab initio picture of multi-neutron correlations in $^7$H and $^8$He, revealing that while compact tetraneutron-like substructures exist, the system is overwhelmingly dominated by a twisted, extended configuration of two correlated dineutron pairs.