From QCD-Based Descriptions to Direct Fits: A Unified Study of Nucleon Electromagnetic Form Factors

This paper presents a unified framework that combines GPD-based contributions and vector-meson exchange to fit nucleon electromagnetic form factors, subsequently deriving stable, global Padé-based analytic parametrizations that accurately describe nucleon structure across a wide range of momentum transfers.

The Gist

Imagine the proton and neutron (the building blocks of atoms) not as solid, tiny billiard balls, but as fuzzy, buzzing clouds of energy. Inside these clouds, tiny particles called quarks and gluons are constantly zipping around, interacting, and changing shape.

Physicists want to take a "snapshot" of these clouds to understand their internal structure. They do this by firing high-speed electrons at them and seeing how the electrons bounce off. The way they bounce is described by something called Electromagnetic Form Factors. Think of these form factors as a "fingerprint" or a "blueprint" that tells us how the charge and magnetism are distributed inside the proton or neutron.

The problem is that taking these fingerprints is tricky. The rules of the game change depending on how hard you hit the target (the energy level).

The Three "Experts"

In this paper, the authors act like a team of detectives trying to solve a mystery using three different experts, each good at a specific type of clue:

1. The Low-Energy Expert (The "VMD" Model):

   • Analogy: Imagine a messenger delivering a package. At low speeds, the package is delivered by a slow, heavy truck (a "vector meson"). This expert is great at describing what happens when the electron hits the proton gently. It sees the proton as a collection of heavy, stable particles.
   • Limitation: If you hit the proton really hard, this truck gets left behind. It can't explain the fast, chaotic motion of the tiny particles inside.

2. The High-Energy Expert #1 (The "ER" Model):

   • Analogy: Now imagine the package is being delivered by a swarm of high-speed drones. This expert looks at the "partons" (the tiny quarks and gluons) and how they behave when the proton is smashed at high speeds. It uses complex math to track how these particles scatter.
   • Limitation: It's great for high speeds but might miss the smooth, gentle transitions that happen at lower speeds.

3. The High-Energy Expert #2 (The "VS24" Model):

   • Analogy: This is another drone expert, but with a slightly different flight pattern. It uses a different set of maps (called "Parton Distribution Functions") to predict how the particles are arranged. It's particularly good at filling in the gaps between the gentle truck and the fast drones.

The Big Idea: The "Master Chef" Recipe

For a long time, physicists tried to use just one of these experts to describe the whole proton. But it was like trying to describe a whole meal using only a recipe for soup, or only a recipe for steak. Neither worked perfectly for the entire menu.

The authors' solution? They decided to hire all three experts and let them work together.

• They created a Combined Model.
• Think of it like a smoothie. You have the heavy truck (low energy), the fast drones (high energy), and the special flight patterns (intermediate energy).
• They didn't just mix them randomly. They added "weights" (like ingredients in a recipe).
  • At low speeds, the "truck" ingredient is dominant.
  • At medium speeds, the "special drone" ingredient takes over.
  • At high speeds, the "fast drone" ingredient becomes the main flavor.

By adjusting these weights and fitting their recipe to real experimental data (the actual photos of the proton), they found the perfect balance. This allowed them to describe the proton's shape accurately from the gentlest touch to the hardest smash.

The "Magic Map" (Padé Approximants)

Once they had their perfect recipe, they wanted to make it easy for other scientists to use. The math was getting very messy and complicated.

So, they created a "Magic Map" (called a Padé Approximant).

• Analogy: Imagine you have a very complicated, winding road with many twists and turns. Instead of giving someone a list of 1,000 specific GPS coordinates to follow, you draw a single, smooth, curved line that perfectly traces the road.
• This "Magic Map" is a simple mathematical formula that captures the complex behavior of the proton. It's smooth, stable, and doesn't break down when you look at extreme speeds.

Why Does This Matter?

This study is a big deal because:

1. It Unifies the View: It bridges the gap between the "slow, heavy" world of particles and the "fast, chaotic" world of quarks.
2. It's Accurate: It fits the real-world data better than previous methods.
3. It's Practical: By creating these smooth "Magic Maps," they give other scientists a reliable tool to use in future experiments, like those at the Large Hadron Collider or electron accelerators.

In short: The authors took three different ways of looking at the proton, mixed them together in the right proportions to match reality, and then drew a simple, smooth map of the result so everyone can understand the shape of the building blocks of our universe.

Technical Summary

1. Problem Statement

Understanding the internal structure of nucleons (protons and neutrons) is a central goal in hadronic physics. Electromagnetic form factors (FFs) encode the spatial and spin distributions of charge and magnetization within nucleons. However, no single theoretical framework currently provides a unified, high-precision description of these form factors across the entire range of momentum transfers ($Q^2$):

• Vector Meson Dominance (VMD) models effectively describe low-$Q^2$ behavior (non-perturbative regime) but fail at higher $Q^2$ where partonic degrees of freedom dominate.
• Generalized Parton Distributions (GPDs) describe the partonic structure and evolution but often struggle to simultaneously reproduce all form factors (Dirac $F_1$, Pauli $F_2$, and Sachs $G_E, G_M$) with high accuracy across all kinematic regions when used in isolation.

The authors aim to bridge this gap by creating a unified framework that combines the strengths of different QCD-inspired approaches and validates them against experimental data.

2. Methodology

The study employs a three-pronged approach, combining three distinct theoretical components into a single weighted model, followed by global analytic parametrization.

A. The Three-Component Combined Model

The authors construct a combined form factor expression for both protons and neutrons by linearly superimposing three models, weighted by parameters $W_1, W_2, W_3$ (where $\sum W_i = 1$):

1. VMD Component ($W_1$): Based on the large-$N_c$ limit of QCD. It models form factors as sums over vector-meson poles (isoscalar $\omega$ and isovector $\rho$ families). This captures low-energy vector meson exchange physics.
2. GPD Component A ($W_2$): Uses the VS24 ansatz combined with KKA10 parton distribution functions (PDFs) at N$^3$LO. This model incorporates flexible exponential tails and is sensitive to PDF evolution.
3. GPD Component B ($W_3$): Uses the ER ansatz (extended Regge-inspired) combined with MRST2002 PDFs at N$^2$LO. This model is designed to capture the sharp fall-off behavior at high $Q^2$ related to the $x \to 1$ limit of GPDs.

The combined expression for Dirac and Pauli form factors is:
$$F_{1,2}^{(p,n)}|_{\text{combined}}(t) = W_1 F_{1,2}^{\text{VMD}}(t) + W_2 F_{1,2}^{\text{VS24}}(t) + W_3 F_{1,2}^{\text{ER}}(t)$$
(Note: The paper defines $W_1$ for VMD, $W_2$ for VS24, and $W_3$ for ER in the text, though Table V assigns dominance differently based on $Q^2$ regions. The fitting procedure determines the optimal weights.)

B. Fitting Strategy

The authors performed global fits to experimental data (specifically Ref. [30], Qattan & Arrington) across four distinct groups of form factors:

• Group 1: Proton $tF_1^p, tF_2^p, G_E^p, G_M^p$.
• Group 2: Proton $G_E^p, G_M^p$ and Neutron $G_M^n$.
• Group 3: Neutron $tF_2^n$ and $G_M^n$.
• Group 4: Neutron $tF_1^n$ and $G_E^n$.

For Group 4 (the most complex neutron channel), a new PDF parametrization was introduced to handle the scarcity of data and large uncertainties. The fitting utilized $\chi^2$ minimization with uncertainties derived from the Minuit covariance matrix.

C. Global Parametrization (Padé Approximants)

To provide stable, analytic representations suitable for phenomenological applications, the authors constructed Padé approximants ($P[N,M]$) based on local Taylor expansions of the fitted data. This step ensures smooth behavior at low $t$ and correct power-law fall-off at high $t$, consistent with perturbative QCD counting rules.

3. Key Contributions

• Unified Framework: Successfully demonstrated that a weighted combination of VMD and two distinct GPD-based models can simultaneously describe proton and neutron form factors across a wide $Q^2$ range ($0 < Q^2 < 5 \text{ GeV}^2$) with high precision.
• Kinematic Dominance Analysis: Identified specific $Q^2$ regions where each physical mechanism dominates:
  • $Q^2 < 1 \text{ GeV}^2$: Dominated by VMD (strong vector meson exchange).
  • $1 < Q^2 < 3 \text{ GeV}^2$: Dominated by VS24 + KKA10 (flexible exponential behavior and PDF sensitivity).
  • $Q^2 > 3 \text{ GeV}^2$: Dominated by ER + MRST2002 (Regge behavior and $x \to 1$ GPD constraints).
• Neutron Form Factor Resolution: Addressed the difficulty of extracting neutron form factors by introducing a two-step fitting procedure and a new PDF parametrization, achieving a $\chi^2/ndf$ of 0.25 for the most difficult channel (indicating excellent agreement within experimental uncertainties).
• Analytic Tools: Provided a set of stable Padé approximant coefficients (Tables VI–IX) for four distinct groups of form factors, offering a reliable tool for future phenomenological studies that avoids extrapolation instabilities.

4. Results

• Fit Quality: The combined model yields excellent agreement with experimental data.
  • Group 1 (Proton): $\chi^2/n = 1.67$.
  • Group 2 (Proton/Neutron M): $\chi^2/n = 1.16$.
  • Group 3 (Neutron): $\chi^2/ndf = 1.31$.
  • Group 4 (Neutron): $\chi^2/ndf = 0.25$.
• Parameter Extraction: The optimal weights ($W_1, W_2, W_3$) and shape parameters (e.g., $\alpha'', \beta, \gamma$ for the VS24 ansatz) were extracted. For instance, in Group 1, the VMD weight ($W_1$) is small (0.035), while the ER component ($W_3$) is dominant (0.925), reflecting the data's sensitivity to the high-$Q^2$ tail in that specific fit.
• Phenomenological Fits:
  • The scaled neutron Dirac form factor ($tF_1^n$) was fitted with a Regge-inspired linear trajectory, yielding a slope $\alpha' = -0.085 \text{ GeV}^{-2}$, significantly flatter than typical baryonic Regge slopes.
  • The neutron magnetic form factor ($G_M^n$) was successfully interpolated using a Padé-like rational function.
• Padé Parametrizations: The constructed Padé functions reproduce the experimental data points and exhibit the correct analytic properties (smoothness at $t \to 0$ and power-law fall-off at $t \to \infty$). Sensitivity bands were generated via Monte Carlo variations of coefficients to estimate model dependence.

5. Significance

This work represents a significant step forward in the phenomenological description of nucleon structure by:

1. Bridging Scales: It effectively unifies the low-energy meson-exchange picture with the high-energy partonic picture, showing that they are not mutually exclusive but rather complementary components of the nucleon wavefunction.
2. Providing a Standard: The extracted Padé coefficients and the combined model parameters serve as a robust, model-dependent reference for future theoretical calculations and experimental analyses.
3. Validating QCD Constraints: The resulting parametrizations satisfy fundamental QCD constraints, including analyticity, Regge behavior, and perturbative counting rules, ensuring they are physically motivated rather than just empirical curve fits.
4. Addressing Neutron Data Scarcity: By developing a specialized fitting strategy for neutron form factors, the paper provides a more reliable extraction of neutron structure information, which is notoriously difficult to obtain due to the lack of free neutron targets.

In conclusion, the paper successfully transitions from disparate QCD-based descriptions to a unified, data-driven framework, offering a precise and analytically stable tool for exploring nucleon electromagnetic structure.