Nucleon strange electromagnetic form factors using $N_f=2+1+1$ twisted-mass fermions at the physical point

This study presents a lattice QCD calculation of the nucleon's strange electromagnetic form factors, including electric and magnetic radii and the magnetic moment, using $N_f=2+1+1$ twisted-mass fermions at the physical point and multiple lattice spacings to achieve a continuum limit with high-precision stochastic noise mitigation for disconnected contributions.

The Gist

Imagine the nucleus of an atom (the proton or neutron) not as a solid marble, but as a bustling, chaotic city. Inside this city, there are three main "citizens" that define its identity: two up-quarks and one down-quark. These are the valence quarks. They are the permanent residents who give the city its name and basic structure.

However, the city is also filled with a swirling, invisible fog of "sea quarks"—quarks and anti-quarks that pop in and out of existence constantly. Among this fog, there is a specific type of citizen called the strange quark. It's the lightest of the "non-resident" quarks. Even though they aren't permanent residents, they still carry an electric charge and a magnetic personality. The question physicists have been asking for decades is: How much does this invisible strange fog actually contribute to the proton's overall electric and magnetic personality?

This paper is the report from a team of scientists who built a digital simulation of this city to answer that question with unprecedented precision.

The Digital City: Lattice QCD

To study these invisible particles, the scientists used a method called Lattice QCD (Quantum Chromodynamics). Think of this as building a giant, 4D digital grid (a lattice) that acts like a pixelated universe. They populated this grid with the rules of physics to simulate how quarks and gluons interact.

Usually, these simulations are like taking a blurry photo: you have to guess what the final picture looks like by taking pictures at different resolutions and trying to smooth them out. This team, however, did something special. They ran their simulation on four different grid sizes (ranging from coarse to very fine) and, crucially, they tuned the "mass" of the particles in the simulation to match the exact, real-world values found in nature.

• The Analogy: Imagine trying to measure the height of a tree. Most people might measure it on a small, low-resolution map and guess the real height. This team measured it on four different maps, all calibrated to the exact real-world scale, and then combined them to get a crystal-clear, "continuum" (perfectly smooth) image without any pixelation.

The Challenge: The "Ghost" Signal

The tricky part of this experiment is that the strange quarks don't stick to the main proton; they float in the "sea." In the simulation, this creates a "disconnected" signal. It's like trying to hear a whisper in a stadium full of cheering fans. The signal from the strange quarks is incredibly faint and gets lost in the "noise" of the simulation.

To fix this, the team used advanced "noise-canceling" techniques:

• Spin-Color Dilution: Imagine trying to listen to a specific instrument in an orchestra by asking the musicians to play one by one in a specific order, rather than all at once. This helps isolate the specific sound.
• Hierarchical Probing: This is like using a high-tech flashlight that scans the stadium in layers, ensuring no dark corner is missed, allowing them to find the faint whisper of the strange quark.

The Findings: What the Strange Quarks Do

Once they cleaned up the noise, they measured two main things:

1. The Strange Electric Radius: How "spread out" the strange quark's electric charge is inside the proton.
2. The Strange Magnetic Moment: How much the strange quark contributes to the proton's magnetism.

The Results:

• The Magnetic Moment: They found that the strange quark does have a magnetic personality, but it is very small. It's like a tiny, barely noticeable tug on the proton's overall magnetism. Their result is consistent with previous studies but is much more precise because they didn't have to guess or "extrapolate" from heavier, unrealistic simulations.
• The Electric Radius: They calculated how far the strange charge extends. Their data suggests a small but measurable spread.
• The Big Picture: When they compared their results to other experiments (which use particle beams to measure these properties indirectly), their numbers fit perfectly within the "confidence zones" of those experiments.

Why This Matters (According to the Paper)

The paper claims this is the first time these specific measurements have been made using a simulation that is:

1. At the physical point (using real-world particle masses, not heavier "fake" ones).
2. In the continuum limit (removing the digital grid artifacts to get a smooth, real-world answer).

By doing this, they provided a very strict "ruler" for experimentalists. If future experiments measure the proton's properties and find a value that doesn't match this simulation, it might mean our understanding of the "sea" of quarks is incomplete. For now, however, the simulation and the experiments agree, giving us a clearer picture of the invisible, strange fog that swirls inside every proton in the universe.

In short: The scientists built a perfect digital model of a proton, filtered out the static noise to hear the faint voice of the strange quark, and confirmed that while this quark is a minor player in the proton's magnetic and electric life, its contribution is now measured with the highest precision ever achieved in a computer simulation.

Technical Summary

1. Problem Statement

The electromagnetic form factors of the nucleon (proton and neutron) are fundamental observables describing their internal structure. While the contributions from valence quarks ($u$ and $d$) are well understood, the contribution from the strange quark sea remains elusive.

• Physics Motivation: Strange quarks constitute the lightest non-valence quarks in the nucleon. Their contribution to electromagnetic form factors ($G_E^s$ and $G_M^s$) provides a direct probe of non-perturbative sea quark dynamics in Quantum Chromodynamics (QCD).
• Experimental Challenges: Experimentally, strange form factors are measured indirectly via parity-violating electron scattering (interference between electromagnetic and weak interactions). Current experimental data (e.g., from SAMPLE, A4, HAPPEX, G0) cannot exclude a zero value for the strange magnetic moment ($\mu_s$) or the electric/magnetic radii, leaving significant uncertainty.
• Lattice QCD Challenges: Previous lattice QCD calculations of strange form factors faced two major limitations:
  1. Disconnected Diagrams: Strange contributions arise from "disconnected" quark loops, which are computationally expensive and suffer from severe statistical noise.
  2. Systematic Errors: Many previous results were calculated at unphysical pion masses and required chiral extrapolation, or were not extrapolated to the continuum limit (zero lattice spacing), introducing significant systematic uncertainties.

2. Methodology

The authors performed a high-precision Lattice QCD calculation addressing the above limitations.

A. Simulation Setup

• Fermion Action: Used $N_f=2+1+1$ twisted-mass fermions with clover improvement. This includes two degenerate light quarks, one strange, and one charm quark.
• Physical Point: Quark masses were tuned to their physical values, ensuring the pion mass is physical ($m_\pi \approx 135-140$ MeV).
• Ensembles: Four distinct ensembles were analyzed with lattice spacings ($a$) of 0.080 fm, 0.068 fm, 0.057 fm, and 0.049 fm. All ensembles maintained similar physical volumes ($m_\pi L \approx 3.6 - 3.9$) to control finite-volume effects.
• Statistics: High-statistics two-point functions were generated. For the disconnected strange loops, stochastic noise mitigation techniques were employed, including:
  • Spin-color dilution.
  • Hierarchical probing (4-dimensional coloring with distance 8).
  • The "generalized one-end trick" specific to twisted-mass fermions.

B. Extraction of Matrix Elements

• Correlation Functions: Calculated nucleon two-point and three-point correlation functions (Fig. 1).
• Ratio Method: Used an optimized ratio of three-point to two-point functions to isolate the ground-state matrix element, canceling overlap with excited states.
• Source-Sink Separation: Analyzed multiple source-sink separations ($t_s$) to ensure ground-state dominance. Plateau fits were performed, showing mild dependence on $t_s$, allowing for a weighted average of results.
• Momentum Transfer: Utilized boosted frames (non-zero sink momenta $\vec{p}'$) in addition to the lab frame. This significantly increased the number of accessible momentum transfer squared ($Q^2$) values without additional computational cost, improving the resolution of the form factor shape.
• Renormalization: The local vector current used for strange contributions was renormalized using constants ($Z_V$) determined in previous works.

C. Continuum Extrapolation and Fitting

• One-Step Approach: The authors employed a global fit to the data that simultaneously parameterizes the $Q^2$ dependence and the lattice spacing ($a^2$) dependence. This allows for a direct extraction of the continuum limit ($a \to 0$) without a separate two-step extrapolation.
• Ansätze: Three functional forms were tested for the form factors:
  1. Dipole Ansatz: For the magnetic form factor.
  2. Galster-like: For the electric form factor.
  3. $z$-expansion: A model-independent expansion used for both, incorporating cut-off effects via $a^2$ terms.
• Model Averaging: Results were averaged using the Akaike Information Criterion (AIC) to account for model dependence.

3. Key Contributions

1. First Continuum Limit at Physical Point: This is the first calculation of nucleon strange electromagnetic form factors that achieves the continuum limit using ensembles strictly at the physical pion mass with $N_f=2+1+1$ flavor dynamical quarks.
2. Advanced Noise Reduction: Demonstrated the efficacy of combining high statistics with hierarchical probing and dilution to control the noise in disconnected strange loops.
3. Boosted Frame Analysis: Showed that utilizing boosted frames significantly enhances the kinematic coverage ($Q^2$ range), leading to more robust fits for radii and moments.
4. Systematic Control: By using four lattice spacings and a one-step global fit, the study rigorously controls discretization errors and chiral extrapolation uncertainties.

4. Results

The study provides precise values for the strange electric radius ($\langle r_E^2 \rangle_s$), magnetic radius ($\langle r_M^2 \rangle_s$), and magnetic moment ($\mu_s$) in the continuum limit.

• Form Factors: The extracted $G_E^s(Q^2)$ and $G_M^s(Q^2)$ show a smooth $Q^2$ dependence consistent with previous lattice data but with significantly reduced errors.
• Strange Magnetic Moment ($\mu_s$): The result is consistent with zero within uncertainties but has a specific sign and magnitude. The authors note that while the sign is consistent with indirect experimental/lattice combinations, the magnitude is smaller than some previous indirect determinations.
• Radii: Precise bounds were placed on the strange electric and magnetic radii.
• Comparison:
  • Lattice: The results agree well with previous lattice calculations (LHPC, $\chi$QCD, Mainz) but offer superior precision due to the continuum/physical-point setup.
  • Experiment: The results provide a tight constraint on the $Q^2=0.1$ GeV$^2$ region, consistent with the 95% confidence ellipses derived from experimental parity-violation data.

5. Significance

• Fundamental QCD: This work provides a rigorous, first-principles determination of the strange sea's contribution to nucleon structure, validating the non-perturbative methods of Lattice QCD.
• Experimental Constraints: The precise lattice results serve as a critical constraint for interpreting parity-violating electron scattering experiments, helping to refine the determination of the proton's weak charge ($Q_{weak}$).
• Methodological Benchmark: The successful application of the one-step continuum extrapolation with physical masses sets a new standard for future lattice calculations of disconnected observables.
• Future Physics: These results are essential for understanding the "proton spin puzzle" and the distribution of momentum and charge within the nucleon, which are crucial for interpreting data from current and future facilities like the Electron-Ion Collider (EIC).

In summary, this paper represents a milestone in lattice QCD by delivering the first high-precision, continuum-limit determination of nucleon strange form factors at the physical point, effectively bridging the gap between theoretical predictions and experimental measurements of sea quark dynamics.