Strangeness of nucleons from $N_f=2+1+1$ lattice QCD

This paper presents the first continuum-limit calculation of nucleon strange electromagnetic form factors using $N_f=2+1+1$ lattice QCD simulations with physical quark masses, yielding results with statistical and systematic uncertainties an order of magnitude smaller than current experimental determinations.

The Gist

Imagine the nucleus of an atom (specifically, a proton or neutron, collectively called "nucleons") not as a solid marble, but as a bustling, chaotic city.

In this city, there are three main "citizens" that define the building's identity: two "Up" quarks and one "Down" quark. These are the permanent residents, the ones you can see from the outside. But, the city is also filled with a foggy, invisible mist of virtual particles popping in and out of existence. This mist is made of "sea quarks."

Most of this mist is made of the same light particles as the residents (Up and Down), but occasionally, a heavier, stranger particle appears: the Strange quark.

The Mystery of the "Ghost" Quark

The strange quark is like a ghost in the machine. It doesn't live there permanently; it only appears because of the wild energy fluctuations of the vacuum (the empty space inside the atom). Because it's so fleeting, it's incredibly hard to measure.

Scientists have been trying to figure out: How much of this "ghost" mist actually affects the electric charge and magnetic pull of the proton?

For decades, experiments tried to catch a glimpse of this ghost by smashing electrons into protons. They found hints that the ghost was there, but the measurements were like trying to hear a whisper in a hurricane—the signal was there, but the noise (errors) was so loud that they couldn't be sure if the ghost was actually contributing anything or if it was just zero.

The New Approach: Building a Perfect Simulation

This paper presents a breakthrough. Instead of trying to catch the ghost in a noisy real-world experiment, the authors built a perfect digital twin of the proton using a supercomputer. This is called Lattice QCD (Quantum Chromodynamics).

Think of it like this:

• Old Simulations: Previous computer models were like trying to simulate a human by starting with a baby and guessing how they would grow up to be an adult. They had to use "heavier" quarks (easier to calculate) and then mathematically guess what would happen if they were lighter. This guessing game (called "chiral extrapolation") introduced a lot of uncertainty.
• This New Simulation: The authors skipped the guessing. They tuned their digital simulation to use quarks with the exact, real-world weight of the particles found in nature. They didn't have to guess; they just simulated reality directly.

The "City" on Four Different Scales

To make sure their digital city was accurate, they didn't just build one version. They built four versions of the proton simulation, each with a different level of detail (resolution), like taking a photo with a camera at 10 megapixels, then 20, then 40, and finally 80.

By comparing these four versions, they could mathematically remove the "pixelation" (the grid lines of the computer simulation) to see the smooth, perfect picture of reality. This is called taking the continuum limit.

What Did They Find?

Once they had their perfect simulation, they measured the "strange" properties of the proton:

1. The Strange Electric Radius: How far does the ghost's electric charge spread out?
2. The Strange Magnetic Radius: How far does the ghost's magnetic pull spread out?
3. The Strange Magnetic Moment: How strong is the ghost's magnetic personality?

The Result: They found that the strange quark does contribute, but the effect is tiny and negative.

• The Precision: Their measurement is so sharp that their "error bars" (the margin of doubt) are ten times smaller than the best experimental measurements we have. It's like going from a blurry, shaky photo to a crystal-clear 8K image.

Why Does This Matter?

You might ask, "Who cares about a tiny ghost quark?"

1. Testing the Rules of the Universe: The Standard Model is the rulebook of physics. To test if the rules are perfect, we need to know every single detail of how protons behave. If we don't know how the "ghost" quarks behave, we can't be sure if a new discovery is a breakthrough or just a misunderstanding of the proton's internal structure.
2. Future Experiments: There is a new experiment coming up (called MESA) that will try to measure these strange properties directly. This paper gives the experimentalists a "target" to aim for. It's like giving a treasure hunter a precise map instead of a vague clue.

The Bottom Line

This paper is a masterclass in precision. By using supercomputers to simulate the proton with real-world weights and removing all the "guesswork" of previous models, the authors have provided the most accurate map of the "strange" ghost inside the proton we have ever seen. They have turned a blurry whisper into a clear, undeniable fact, helping us understand the hidden, invisible layers of the universe that make up everything we see.

Technical Summary

1. Problem Statement

The internal structure of the nucleon (proton and neutron) is governed by Quantum Chromodynamics (QCD). While valence quarks (up and down) dominate the nucleon's properties, the "sea" of virtual quark-antiquark pairs also contributes. The strange quark, being the lightest non-valence quark, contributes to the nucleon's structure purely through vacuum fluctuations.

• Scientific Challenge: Determining the strange electromagnetic form factors ($G_E^s$ and $G_M^s$) is crucial for understanding non-perturbative QCD vacuum dynamics and is a necessary input for extracting the proton's weak charge from parity-violating electron scattering experiments (e.g., Qweak, MESA).
• Limitations of Previous Work:
  • Experimental: Measurements via parity-violating asymmetries (SAMPLE, A4, HAPPEX, G0) have large uncertainties, often failing to exclude a zero value for the strange magnetic moment and radii.
  • Theoretical (Lattice QCD): Previous lattice calculations either used quark masses heavier than physical values (requiring chiral extrapolations) or mixed physical and unphysical ensembles. These extrapolations introduce significant systematic uncertainties.

2. Methodology

The authors performed a high-precision lattice QCD calculation using a first-principles approach that eliminates the need for chiral extrapolations.

• Lattice Setup:

  • Fermion Action: $N_f = 2+1+1$ (two degenerate light, one strange, and one charm quark) using clover-improved twisted mass fermions.
  • Physical Point: All quark masses were tuned to their physical values.
  • Continuum Limit: Simulations were conducted on four ensembles with different lattice spacings ($a \approx 0.08, 0.068, 0.057, 0.049$ fm) to perform a controlled continuum extrapolation directly at the physical pion mass.
  • Volume: Lattice volumes were chosen to ensure $m_\pi L > 3.6$, minimizing finite-volume effects.

• Computational Techniques:

  • Disconnected Diagrams: The strange quark contribution arises from "disconnected" quark loops (where the current couples to a sea quark loop not connected to the valence quarks). These are computationally expensive and noisy.
  • Noise Reduction: The authors employed Gaussian smearing for nucleon interpolating fields, APE smearing for gauge links, and advanced stochastic techniques for the quark loops, including full dilution (spin/color), hierarchical probing, and low-mode deflation.
  • Correlation Functions: They computed two-point and three-point correlation functions to extract the nucleon matrix elements of the strange vector current.
  • Renormalization: The local vector current was renormalized using non-perturbative renormalization constants ($Z_V$) determined via joint polynomial fits across all lattice spacings.

• Data Analysis & Extrapolation:

  • Plateau Fits: Ground-state matrix elements were extracted using optimized ratios of three-point to two-point functions, fitting for time-independent plateaus.
  • $Q^2$ Dependence: To extract radii and moments, the momentum transfer dependence ($Q^2$) of the form factors was parameterized using three models: Dipole, Galster-like, and $z$-expansion.
  • Continuum & Model Averaging: The lattice spacing dependence was incorporated into the fit parameters ($f(a^2) = f_0 + f_1 a^2$). The final results were obtained by model-averaging the different parameterizations using the Akaike Information Criteria (AIC) to account for systematic uncertainties in the fit model.

3. Key Contributions

• First Physical-Point Continuum Limit: This is the first calculation of strange electromagnetic form factors using only ensembles at physical quark masses to reach the continuum limit. This completely removes the systematic uncertainty associated with chiral extrapolations.
• High Precision: The statistical and systematic errors achieved are an order of magnitude smaller than those of previous experimental determinations.
• Comprehensive Error Budget: The study provides a rigorous separation of statistical errors and systematic uncertainties arising from the $Q^2$-dependence fit and continuum extrapolation.

4. Key Results

The authors extracted the strange electric radius ($\langle r^2_E \rangle_s$), strange magnetic radius ($\langle r^2_M \rangle_s$), and strange magnetic moment ($\mu_s$). The results (including statistical and systematic errors) are:

• Strange Electric Radius:
  $$\langle r^2_E \rangle_s = -0.00545(49)(26) \, \text{fm}^2$$
• Strange Magnetic Radius:
  $$\langle r^2_M \rangle_s = -0.01212(280)(72) \, \text{fm}^2$$
• Strange Magnetic Moment:
  $$\mu_s = -0.01792(195)(18)$$

Note: The negative values indicate a specific distribution of strange charge and magnetization, though the magnitudes are small.

The momentum transfer dependence of the form factors ($G_E^s$ and $G_M^s$) was successfully described using the $z$-expansion, which provided the most reliable constraints compared to dipole or Galster-like fits.

5. Significance

• Theoretical Impact: The results provide a definitive benchmark for non-perturbative QCD, confirming that strange quark contributions are small but non-zero. The precision achieved validates the lattice methodology of working directly at the physical point.
• Experimental Guidance: The high-precision lattice results provide stringent constraints on the strange form factors. As shown in the paper's comparison with experimental 95% confidence contours, the lattice results significantly narrow the allowed parameter space.
• Future Applications: These precise values are critical inputs for:
  • The MESA experiment at Mainz, which aims to measure parity-violating asymmetries at low momentum transfer.
  • Refining the extraction of the proton weak charge, a key test of the Standard Model.
  • Reducing hadronic structure uncertainties in precision electroweak measurements.

In summary, this work represents a major milestone in lattice QCD, demonstrating that high-precision determination of sea-quark effects is achievable without relying on extrapolations, thereby providing the most accurate theoretical determination of the nucleon's strange electromagnetic structure to date.