⚛️ nuclear theory

Analyzing momentum distributions in light nuclei with the operator product expansion

This paper utilizes the operator product expansion within Pionless effective field theory to analyze high-momentum nucleon distributions in light nuclei, successfully extracting local-operator matrix elements from variational Monte Carlo data to accurately reconstruct single-nucleon distributions and establish approximate relations between single- and two-nucleon contributions in potentials lacking explicit non-nucleon degrees of freedom.

Original authors: Jiexin Yu, Bingwei Long

Published 2026-09-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Jiexin Yu, Bingwei Long

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Inside the heart of every atom lies a bustling crowd of protons and neutrons, bound together by forces so strong they defy our everyday intuition. While we often picture these particles as sitting quietly in neat shells, the reality at the smallest scales is far more chaotic. Occasionally, two particles crash into each other with such ferocity that they momentarily break free from the average flow, shooting off with speeds far higher than the rest of the crowd. These fleeting, high-speed encounters are known as short-range correlations. Understanding how often they happen and how much energy they carry is crucial for physicists, as these rare events hold the key to the fundamental rules governing nuclear matter. To study them, scientists look at the "momentum distribution" of a nucleus—a map showing how many particles are moving at any given speed. The challenge has always been separating the slow, steady motion of the average particle from the wild, high-speed bursts of these rare collisions.

A team of researchers has now developed a new way to untangle this complexity, offering a clearer picture of these high-speed tails in several light nuclei. By combining two powerful theoretical tools, they managed to isolate the specific signatures of these violent encounters. One tool, known as the operator product expansion, acts like a mathematical sieve, separating the short-range, high-energy physics from the long-range, low-energy structure of the nucleus. The second tool, called pionless effective field theory, simplifies the problem by ignoring the exchange of certain heavy particles (pions) that usually complicate nuclear calculations, allowing the researchers to focus purely on the direct interactions between nucleons. Together, these methods allowed the team to analyze the momentum distributions of helium-3, helium-4, lithium-6, and carbon-12 with unprecedented precision.

The researchers began by examining data generated from sophisticated computer simulations, which modeled the behavior of these nuclei using well-established nuclear forces. A major hurdle in previous attempts to understand these distributions was that the signals for the two most common types of high-speed collisions looked almost identical. It was like trying to hear two different instruments playing the same note at the same volume; the individual contributions were impossible to distinguish. To solve this, the team turned their attention to the two-particle distributions. Instead of looking at the speed of a single particle, they analyzed the relative speeds of pairs of particles. By sorting these pairs based on their spin and isospin—quantum properties that describe their internal orientation and type—they could naturally separate the two types of collisions. One type involved particles with opposite spins, while the other involved particles with aligned spins. This separation allowed them to extract the specific "weights" or contributions of each collision type with high accuracy.

Once they had isolated these individual contributions, the team used them to reconstruct the full picture of how single particles move within the nucleus. The results were strikingly accurate. In the high-momentum region, where particles move at speeds exceeding 400 million meters per second, the reconstructed curves matched the computer simulation data almost perfectly. This agreement served as a strong validation of their theoretical framework. The analysis revealed that for all the nuclei they studied, the collisions involving aligned spins were the dominant source of high-momentum particles. This finding aligns with experimental observations that neutron-proton pairs are the primary drivers of these short-range correlations.

The study also uncovered a surprising simplicity in how these distributions behave. In the specific models used for the simulations, which do not include complex internal particle exchanges, the researchers found that the distribution of single particles is directly proportional to the distribution of particle pairs. This means that if you know how the pairs are moving, you can predict how the individual particles are moving with a simple scaling factor. This relationship, which holds true for the high-speed tails of the distribution, suggests that the complex dance of the nucleus can be understood through a few fundamental building blocks. While the method worked exceptionally well for the lighter nuclei, the researchers noted that the agreement began to fade slightly for the heavier carbon-12, hinting that the density of the nucleus plays a role in how these correlations form.

Looking ahead, the authors plan to refine their approach by calculating these contributions directly from first principles, rather than extracting them from existing simulations. They also intend to expand their framework to include more complex, three-particle interactions, which are currently being hunted for in experimental laboratories. By testing their theory against real-world data from high-precision electron-scattering experiments, they hope to confirm whether these universal patterns hold true across the entire periodic table. For now, this work provides a robust, model-independent way to describe the high-energy behavior of atomic nuclei, turning a previously blurry picture of nuclear chaos into a sharp, understandable image.

Technical Summary: Analyzing Momentum Distributions in Light Nuclei with the Operator Product Expansion

Problem Statement
The paper addresses the challenge of disentangling the short-range correlations (SRCs) from the long-distance structure of nucleon momentum distributions in atomic nuclei. While effective field theories (EFTs) and the operator product expansion (OPE) offer a systematic framework for this, a specific difficulty arises in the single-nucleon momentum distribution, $\rho(k)$ . In the large-momentum limit, the leading contributions from the spin-singlet ( $^1S_0$ ) and spin-triplet ( $^3S_1$ ) channels exhibit nearly degenerate momentum dependence in their Wilson coefficients. This degeneracy makes it difficult to independently extract the corresponding local-operator matrix elements (LOMEs) when fitting only the single-nucleon distribution data.

Methodology
The authors employ the OPE in conjunction with Pionless Effective Field Theory (EFT) to analyze high-momentum tails of nucleon distributions. The core methodology involves:

Operator Definition: Defining nonlocal operators for single-nucleon ( $\Omega(k)$ ) and two-nucleon ( $\Pi(k)$ ) momentum distributions within a quantum field theory framework.
OPE Formulation: Expanding these nonlocal operators into a series of local composite operators multiplied by nucleus-independent Wilson coefficients. The leading local operators correspond to the $^1S_0$ ( $O_s$ ) and $^3S_1$ ( $O_t$ ) short-range correlations.
Matching: Determining the Wilson coefficients by matching the OPE-EFT calculations to the AV18 two-nucleon potential and Urbana-X three-nucleon potential in the two-body scattering sector. This is done by equating matrix elements between NN scattering states in the AV18 model and the OPE-EFT.
Decoupling Contributions: To resolve the degeneracy between $^1S_0$ and $^3S_1$ contributions in the single-nucleon distribution, the authors construct the OPE for the two-nucleon momentum distribution. By utilizing spin-isospin dependent pair distributions ( $\rho^{2N}_{ST}$ ) from Variational Monte Carlo (VMC) calculations, they separate the $^1S_0$ and $^3S_1$ channels into distinct spin-isospin components ( $\rho^{2N}_{01}$ and $\rho^{2N}_{10}$ ).
Extraction and Reconstruction: The LOMEs ( $F_s(A)$ and $F_t(A)$ ) are extracted from the VMC two-nucleon distributions by observing the plateau of the ratios $\rho^{2N}_{ST}(k)/\tilde{W}_{s,t}(k)$ at large momenta. These extracted LOMEs are then used to reconstruct the single-nucleon distribution $\rho(k)$ .

Key Contributions

Resolution of Degeneracy: The paper introduces the OPE of the two-nucleon momentum distribution as a necessary tool to separate the $^1S_0$ and $^3S_1$ contributions, which are indistinguishable in fits to the single-nucleon distribution alone due to the near-degeneracy of their Wilson coefficients.
Extraction of LOMEs: The authors successfully extract the LOMEs for $^3\text{He}$ , $^4\text{He}$ , $^6\text{Li}$ , and $^{12}\text{C}$ directly from VMC data rather than calculating them via ab initio Pionless EFT, treating them as free parameters to test the framework.
Approximate Relations in Potential Models: For nuclear potentials lacking explicit non-nucleon degrees of freedom (like AV18), the authors demonstrate that the Wilson coefficients for the single- and two-nucleon OPEs are proportional. This leads to approximate relations between the single-nucleon and two-nucleon distributions (e.g., $\rho(k) \approx 2\rho^{2N}(k)$ ), a result previously derived using the Generalized Contact Formalism (GCF).

Results

Reconstruction Accuracy: Using the LOMEs extracted from two-nucleon distributions, the reconstructed single-nucleon momentum distributions agree accurately with the original VMC data in the high-momentum region ( $k \gtrsim 400$ MeV).
Dominance of $^3S_1$ : The analysis confirms that the $^3S_1$ component dominates the high-momentum tails for all investigated nuclei, consistent with the experimental observation of neutron-proton dominance in SRCs.
Convergence Behavior: The OPE expansion converges faster for $^3\text{He}$ and $^4\text{He}$ compared to $^6\text{Li}$ and $^{12}\text{C}$ . The authors note that while $^4\text{He}$ is dense, the convergence is surprisingly robust, potentially suggesting a larger domain of applicability for contact-only interactions than expected.
Wilson Coefficient Stability: The matching procedure yields Wilson coefficients that are relatively stable (variations $\sim 10\%$ ) with respect to the matching momentum, comparable to the truncation error of Pionless EFT.

Significance and Claims
The paper claims to provide a model-independent description of the high-momentum tails of nucleon momentum distributions for light nuclei. By combining OPE with Pionless EFT, the framework successfully separates short-distance dynamics (encoded in Wilson coefficients) from long-distance nuclear structure (encoded in LOMEs).

The authors emphasize that their work serves as a "proof-of-principle" application. They do not claim to have solved the ab initio calculation of LOMEs within Pionless EFT; rather, they demonstrate that the OPE-EFT framework is consistent with high-precision VMC data when LOMEs are extracted from that data. The success of the reconstruction in the high-momentum regime provides a nontrivial cross-check of the OPE-EFT formalism. The paper concludes that this approach offers a systematic way to study SRC physics, though it acknowledges that future work is needed to calculate LOMEs directly from ab initio methods and to confront the framework with experimental inclusive electron-scattering data.

More like this