Radiative Corrections to Elastic Lepton-Proton Scattering with Focus on Two-Photon-Exchange Diagrams

This paper presents a complete next-to-leading-order calculation of QED radiative corrections to elastic electron-proton and muon-proton scattering, with a specific focus on the structure-dependent two-photon-exchange diagrams to address discrepancies like the proton radius puzzle and test lepton universality.

The Gist

Imagine the proton as a tiny, bustling city inside an atom. For decades, scientists have tried to map this city by firing tiny "scouts" (electrons or muons) at it and watching how they bounce off. By studying the bounce, they can figure out how the city's charge and magnetism are distributed.

However, the city isn't just a solid block; it's a complex, fuzzy cloud of particles. When a scout hits the city, the interaction isn't always as simple as a single billiard ball hitting another. Sometimes, the scout and the city exchange two messengers (photons) instead of just one. This is called the Two-Photon Exchange (TPE).

For a long time, scientists used a "one-messenger" rule to calculate these bounces. But as their measuring tools got incredibly precise, they started seeing cracks in the map. Two famous puzzles emerged:

1. The Proton Form Factor Puzzle: Different ways of measuring the city's shape gave conflicting results.
2. The Proton Radius Puzzle: Measuring the city's size with electrons gave a different answer than measuring it with muons (a heavier cousin of the electron).

The authors of this paper, Daniel Crowe, Syed Mehedi Hasan, and Doreen Wackeroth, decided to fix the math behind these measurements. Here is what they did, explained simply:

1. The "Perfect Map" Problem

Think of the old math (called the "Born approximation") as a map that assumes the proton is a perfect, smooth sphere. It works okay for rough estimates, but it misses the details. The authors realized that to get a truly accurate map, they needed to account for the messy reality: the proton is made of quarks, and its "shape" changes depending on how hard you hit it.

They created a complete, high-definition calculation of the "radiative corrections." In everyday terms, this means they calculated all the tiny, invisible "glitches" and "echoes" that happen during the collision. Specifically, they focused on the Two-Photon Exchange, which is the most complex part of the glitch.

2. The "Shape-Shifting" Challenge

The tricky part of their job was that the proton's shape isn't static. It's like a shape-shifting balloon.

• The Old Way: Previous calculations often treated the proton as if its shape was fixed, ignoring how the "messengers" (photons) interacted with the proton's internal structure at different speeds.
• The New Way: The authors built a model where the proton's shape changes dynamically based on the momentum of the messengers. They treated the proton's internal structure as a "loop" that depends on the speed and energy of the particles involved.

To do this, they used two different, powerful mathematical "engines" (Passarino-Veltman reduction and Integration-by-Parts identities). It's like solving a massive jigsaw puzzle using two completely different strategies. When both strategies produced the exact same picture, they knew their map was correct.

3. The Results: Electron vs. Muon

They tested their new map against real-world data from experiments where electrons and muons hit protons.

• The Electron Effect: When electrons hit the proton, the "glitches" (corrections) are huge—sometimes changing the result by 20%. This is because electrons are light and move very fast, making them sensitive to the proton's fuzzy edges.
• The Muon Effect: Muons are much heavier. They act more like a heavy bowling ball hitting a pin, so the "glitches" are much smaller.
• The Two-Photon Surprise: They found that the "two-messenger" exchange (TPE) is significant. It can change the calculated probability of a bounce by up to 15% in certain conditions. This is a big deal because it means the old "one-messenger" maps were missing a major piece of the puzzle.

4. Why This Matters (According to the Paper)

The authors compared their new, detailed map with existing experimental data (from experiments like CLAS and OLYMPUS). They found that their new calculations match the real-world data much better than the old approximations did.

They also compared their results with other theoretical predictions. While there were small differences, they found that these differences often came down to how the proton's shape was described in the math (the "form factor"). Their work shows that to solve the proton puzzles, we need to be very precise about how we describe the proton's internal structure, not just the collision itself.

The Bottom Line

This paper is like a team of cartographers who realized their map of a city was missing the winding alleyways and hidden courtyards. They didn't just draw the main streets; they mapped the entire complex, dynamic structure of the proton's interior.

By doing this, they provided a more accurate "rulebook" for scientists to use when analyzing data from particle accelerators. This helps ensure that when we measure the size or shape of the proton, we aren't being fooled by the messy, invisible "echoes" of the collision. Their work is a foundational step toward finally solving the proton radius and form factor puzzles, ensuring that the map of the atomic world is as accurate as the tools we use to draw it.

Technical Summary

Technical Summary: Radiative Corrections to Elastic Lepton-Proton Scattering with Focus on Two-Photon-Exchange Diagrams

Problem Statement
Elastic lepton-proton ($ep$ and $\mu p$) scattering serves as a primary tool for probing nucleon structure, specifically the spatial distributions of charge and magnetization via electromagnetic form factors. However, recent high-precision experiments have revealed significant discrepancies, notably the "proton form factor puzzle" (divergence between Rosenbluth separation and polarization transfer methods for $G_E/G_M$) and the "proton radius puzzle" (discrepancy between muonic hydrogen spectroscopy and electron scattering/atomic hydrogen results). These tensions necessitate a rigorous theoretical framework for radiative corrections (RC). While the one-photon-exchange (OPE) approximation (Born approximation) is standard, it is insufficient for sub-percent precision. Specifically, the two-photon-exchange (TPE) contribution, which depends on the proton's internal structure, requires a complete calculation beyond soft-photon approximations to accurately interpret experimental data and test lepton universality.

Methodology
The authors present a complete Next-to-Leading Order (NLO) QED calculation for elastic lepton-proton scattering, encompassing both virtual corrections and real soft-photon emission. The calculation explicitly retains finite lepton and proton masses to avoid collinear divergences and regulates infrared (IR) divergences using a fictitious photon mass ($\lambda_\gamma$), with consistency checks performed using dimensional regularization.

Key methodological features include:

1. Loop-Momentum Dependent Form Factors: Unlike many previous treatments that assume point-like interactions or soft-photon approximations for the TPE, this work systematically incorporates proton form factors ($F_1, F_2$) that depend on the loop momentum within the virtual corrections. The form factors are parametrized using monopole ($n=1$) and dipole ($n=2$) functions with a cutoff scale $\Lambda$.
2. Decomposition of Corrections: The virtual matrix element is decomposed by powers of the proton charge $Z$:
   • $Z^0$: Leptonic corrections (lepton self-energy, vertex, vacuum polarization).
   • $Z^2$: Proton-side corrections (proton self-energy, vertex).
   • $Z^1$: Two-Photon Exchange (TPE) diagrams (box and crossed-box).
3. Calculation Techniques: The TPE contribution involves non-standard denominator structures due to the loop-momentum dependent form factors (propagators raised to powers $n$). The authors reduce these to the Passarino-Veltman (PV) scalar basis using two independent, cross-checking methods:
   • Partial Fractioning and Differentiation: Separating form-factor denominators and using differentiation with respect to the cutoff scale $\Lambda$ to handle higher powers of propagators.
   • Integration-by-Parts (IBP): Applying IBP identities to reduce integrals to a set of 23 master integrals (using LiteRed, Package-X, and Collier).
4. Analytic Results: The authors provide complete analytic results for the virtual corrections, including explicit expressions for the form factors of the proton vertex and the TPE box diagrams, suitable for implementation in Monte Carlo event generators.

Key Contributions and Results

• Complete TPE Calculation: The paper provides an independent, complete calculation of the TPE contribution at NLO, retaining the full loop-momentum dependence of the proton form factors. This allows for a model-dependent assessment of the "hard" (IR-finite) part of the TPE correction, distinct from the model-independent IR-divergent parts often subtracted in experimental analyses.
• Numerical Impact:
  • The leptonic ($Z^0$) contribution dominates the total correction, reaching $\sim 20\%$ for $ep$ scattering due to large collinear logarithms, while the $\mu p$ correction is significantly smaller.
  • The proton-side ($Z^2$) correction is small.
  • The structure-dependent TPE ($Z^1$) contribution is found to be significant, reaching up to $-15\%$ for $E_1 = 3$ GeV and $Q^2 = 5$ GeV$^2$.
• Comparison with Existing Models:
  • The results agree with the Maximon-Tjon (MTj) and Mo-Tsai (MoT) soft-photon approximations in the IR-divergent limit.
  • The "hard" TPE contribution ($\delta_{\gamma\gamma}$) is shown to depend on the choice of form-factor parametrization (monopole vs. dipole) and the momentum transfer $Q^2$. Differences between monopole and dipole parametrizations become pronounced at $Q^2 \gtrsim \Lambda^2$.
  • Comparisons with Blunden-Melnitchouk-Tjon (BMT) and other literature results show qualitative agreement, with residual differences attributed to variations in form-factor ansätze.
• Comparison with Data: The calculated ratio of positron to electron scattering cross sections ($R_{2\gamma}$) is compared with CLAS and OLYMPUS data. The results are consistent with experimental uncertainties, showing a TPE effect of at most 3.5% in the dipole case. The authors note that including inelastic intermediate states (e.g., $\Delta(1232)$) could further modify these results.

Significance and Claims
The authors claim that this work provides a necessary theoretical input for the next generation of precision lepton-proton scattering experiments, such as MUSE, PRad, and programs at MAMI and Jefferson Lab. By offering a complete analytic calculation that includes loop-momentum dependent form factors, the paper aims to:

1. Enable more accurate extraction of Sachs form factors and the proton radius by reducing theoretical uncertainties associated with RC.
2. Provide a framework for testing lepton universality by comparing $ep$ and $\mu p$ scattering with consistent theoretical corrections.
3. Supply a stand-alone Fortran code (available upon request) and analytic results to facilitate the implementation of these corrections in experimental event generators.

The paper explicitly limits its scope to the elastic proton intermediate state. It acknowledges that inelastic contributions (e.g., $\pi N$ states, nucleon resonances, and partonic regimes at high $Q^2$) are outside the current scope but represent natural directions for future work. The significance lies in the rigorous treatment of the elastic TPE contribution, which is a prerequisite for isolating and quantifying these higher-order inelastic effects in future analyses.