Resolving the structure of bound states using lattice quantum field theories

This paper presents the first lattice calculation of a two-to-two particle matrix element using pionless effective field theory to demonstrate that while finite-volume formalism is unnecessary for deep-bound states, it is critical for accurately determining the infinite-volume elastic form factors and charge radii of shallow-bound states like the deuteron.

The Gist

Imagine you are trying to understand the shape and structure of a tiny, fragile bubble floating in a room. This bubble represents a deuteron, a simple atomic nucleus made of two particles (a proton and a neutron) stuck together.

In the real world, this bubble floats freely in infinite space. But in the world of Lattice Quantum Field Theory (the computer simulations physicists use to study the universe's building blocks), we can't simulate infinite space. Instead, we have to put our bubble inside a digital box (a finite 3D grid).

Here is the problem: When you put a delicate bubble inside a small box, the walls of the box push on it, squishing it and changing its shape. If you try to measure the bubble's size while it's being squished by the walls, you get a distorted, wrong answer.

This paper is about a new, clever mathematical "magic trick" that allows physicists to look at the bubble inside the box, understand how the walls are squishing it, and mathematically undo that distortion to see what the bubble actually looks like in the real, infinite world.

The Two Scenarios: The Rock vs. The Bubble

The authors tested their method on two very different types of "bubbles" (bound states):

1. The Deep-Bound State (The Rock): Imagine a heavy rock glued to the floor. It's so tightly stuck that the walls of the room don't bother it at all.

   • The Result: When the authors simulated this, they found that the "box walls" didn't matter. The rock looked the same inside the box as it would in the open field. Their fancy new math wasn't strictly necessary here, but it worked perfectly fine.

2. The Shallow-Bound State (The Bubble): Now imagine a soap bubble floating just inches from the wall. The air pressure from the wall is pushing on it, making it wobble and change shape.

   • The Result: Without their new math, the simulation gave a chaotic, nonsensical result. The "bubble" looked like it had multiple sizes at once (a mathematical glitch called being "multi-valued").
   • The Fix: When they applied their new formalism (the "magic trick"), the distortion vanished. The chaotic data smoothed out into a single, clear, sensible shape. This proved that for delicate, shallow bound states, this new method is absolutely critical.

The "Magic Trick": How It Works

Think of the simulation as a photographer taking a picture of an object in a funhouse mirror (the finite box). The mirror distorts the image.

1. The Snapshot: They take a picture of the object inside the box (the "finite-volume matrix element").
2. The Mirror Map: They use a complex set of rules (called Lüscher formalism and its newer extensions) to understand exactly how the mirror bends the light.
3. The Correction: They use these rules to mathematically "straighten out" the image, revealing the true, undistorted object (the infinite-volume form factor).

Why Does This Matter?

For decades, physicists have been great at simulating how particles bounce off each other (scattering). But simulating how they react to a probe (like a photon hitting a nucleus) is much harder, especially when the nucleus is barely holding together.

This paper is a proof of concept. It says:

"We have successfully built a bridge between the distorted world of our computer boxes and the real world. We can now take messy, distorted data from a simulation and extract the true, physical properties of delicate nuclear structures."

The "Charge Radius" Discovery

One of the things they measured was the charge radius (essentially, how big the bubble is).

• They found that as the bubble gets "shallower" (looser), it gets physically larger.
• Their results matched a famous theoretical prediction called the "Anomalous Threshold."
• The Analogy: It's like predicting that as a soap bubble gets closer to popping, it gets bigger and more sensitive to the slightest breeze. Their math confirmed this intuition perfectly.

The Bottom Line

This paper is a major step forward. It shows that we can now reliably use supercomputers to study the internal structure of light nuclei (like the deuteron) even when they are barely holding together.

Previously, trying to study these "shallow" states in a computer box was like trying to weigh a feather on a scale that was shaking violently. This paper provides the shock absorbers for that scale, allowing us to finally get a precise measurement. This opens the door to understanding how the universe's most basic building blocks interact, which is crucial for everything from understanding how stars shine to searching for new physics beyond our current theories.

Technical Summary

1. Problem Statement

The primary goal of nuclear physics is to predict the properties and reactions of nuclei directly from Quantum Chromodynamics (QCD). While Lattice QCD (LQCD) has successfully determined elastic form factors for single nucleons, extending these calculations to light nuclei (specifically two-nucleon systems like the deuteron) remains a significant challenge.

The core difficulty lies in extracting electroweak transition amplitudes (specifically $2 + J \to 2$, where a current $J$ interacts with a two-particle state) from finite-volume lattice simulations.

• Finite-Volume Artifacts: LQCD calculations are performed in a finite spatial volume ($L$). For shallow bound states (like the deuteron near the physical point), the wavefunction extends significantly beyond the lattice size. Consequently, the matrix elements calculated directly in finite volume are multi-valued and violate analyticity, rendering them physically meaningless without correction.
• Formalism Gap: While formalisms exist for single-particle currents ($J \to 2$) or purely hadronic scattering ($2 \to 2$), a rigorous, non-perturbative framework to map finite-volume matrix elements to infinite-volume $2 + J \to 2$ amplitudes (coupling two-particle states via a current insertion) was previously untested in a lattice calculation.

2. Methodology

The authors perform the first lattice calculation of a two-to-two particle matrix element using a leading-order (LO) pionless Effective Field Theory (EFT/$\pi$) as a toy model. This approach allows for exact diagonalization of the Hamiltonian in moderate volumes.

Key Methodological Steps:

1. Model Setup:

   • The system consists of two non-relativistic nucleons (proton and neutron) in a finite 3D cubic volume.
   • The interaction is a contact term tuned to support a bound state. By varying the coupling constant ($g_0$ or $c$), the authors simulate both deeply bound states (large binding energy) and shallow bound states (binding energy near zero).
   • The transfer matrix is constructed and diagonalized to obtain the finite-volume energy spectrum and eigenvectors.

2. Matrix Element Calculation:

   • The authors calculate the finite-volume matrix elements of a conserved local vector current ($J^\mu$) between the ground state and boosted states.
   • The current couples only to the proton (charged), while the neutron is neutral.

3. Application of the Formalism (Refs. [41–43]):

   • Lüscher Quantization: The finite-volume spectrum is used to constrain the infinite-volume hadronic scattering amplitude ($M$) via the Lüscher formalism.
   • Lellouch-Lüscher Factor: The finite-volume matrix elements are corrected using the Lellouch-Lüscher factor ($R$) to account for the normalization difference between finite and infinite volume states.
   • Amplitude Decomposition: The infinite-volume amplitude $W_\mu$ for the process $2 + J \to 2$ is decomposed into:
     • A pole term (single-particle propagation).
     • A divergence-free term ($W_{df}$) containing the dynamical information.
   • Triangle Functions: The formalism utilizes relativistic triangle functions ($G, G_\mu$) arising from the loop integrals, which encode the anomalous threshold singularities.
   • Extraction of Form Factor: The bound-state elastic form factor ($f_B$) is extracted from the residue of the amplitude at the bound-state pole. The relation is:
     $$(P_f + P_i)^\mu f_B(Q^2) = g^2 \left( A^\mu + (P_f + P_i)^\mu G - 2G^\mu \right)$$
     where $A^\mu$ is a short-distance, real function determined by the lattice data.

3. Key Contributions

• First Implementation: This is the first lattice QFT calculation implementing the formalism for $2 + J \to 2$ amplitudes to extract bound-state form factors.
• Validation of Formalism: The study provides a proof-of-principle that the formalism correctly maps finite-volume observables to infinite-volume quantities, specifically handling the complex analytic structure required for bound states.
• Deep vs. Shallow Bound States: The work explicitly demonstrates the necessity of the formalism.
  • For deeply bound states, finite-volume effects are exponentially suppressed ($e^{-\kappa L}$), and the raw matrix elements approximate the form factor well.
  • For shallow bound states, the raw matrix elements are multi-valued and unphysical. The formalism is absolutely critical to remove these artifacts and recover a single-valued, analytic form factor.
• Anomalous Threshold Dominance: The analysis reveals that for shallow bound states, the form factor's momentum dependence is driven almost entirely by the anomalous threshold (the singularity in the triangle diagram), rather than the short-distance physics ($A^\mu$).

4. Results

• Spectrum and Scattering: The authors successfully reproduced the finite-volume spectrum for various couplings and extracted the scattering length and effective range, confirming the presence of bound states.
• Form Factor Extraction:
  • The extracted form factors $f_B(Q^2)$ are smooth, single-valued functions of $Q^2$, satisfying charge conservation ($f_B(0)=1$).
  • The function $A_0$ (short-distance contribution) was found to be nearly constant with respect to $Q^2$, confirming that the $Q^2$ dependence is dominated by the kinematic triangle functions.
• Charge Radius:
  • The charge radius $\langle r_C^2 \rangle$ was calculated from the derivative of the form factor at $Q^2=0$.
  • As the binding energy approaches zero (shallow limit), the calculated radius converges to the universal prediction derived from the asymptotic wavefunction:
    $$\langle r_C^2 \rangle = \frac{1}{8\kappa^2}$$
    where $\kappa$ is the binding momentum. This agreement validates the formalism's ability to capture the correct long-distance physics.
• Comparison with LO EFT: The lattice results show excellent agreement with the analytic leading-order pionless EFT prediction for the form factor, $f_B^{LO}(Q^2) = \frac{4\kappa}{Q} \tan^{-1}(\frac{Q}{4\kappa})$.

5. Significance and Outlook

• Path to Physical Nuclei: This work establishes the necessary theoretical and computational framework to study the structure of the deuteron and other light nuclei directly from LQCD at physical quark masses. Without this formalism, extracting meaningful structural information (like charge radii) for shallow bound states from lattice data is impossible due to finite-volume artifacts.
• Electroweak Reactions: The methodology enables the calculation of crucial electroweak processes relevant to neutrino physics and astrophysics, such as:
  • Coherent neutrino-deuteron scattering ($\nu + d \to \nu + d$).
  • Photodisintegration of the deuteron ($\gamma^* + d \to n + p$).
• Future Work: The authors note that the current formalism assumes spinless particles. The next critical step is generalizing this framework to include the intrinsic spin of nucleons and extending the analysis to virtual bound states (resonances), which will allow for a complete resolution of nuclear structure from first principles.

In summary, this paper successfully bridges the gap between finite-volume lattice data and infinite-volume nuclear structure, proving that the complex formalism required for shallow bound states is not only necessary but highly effective in recovering physical observables.