Pion Parton Distribution Functions in the Light-Cone Quark Model and Experimental Constraints

This paper investigates pion valence quark parton distribution functions using a light-cone quark model, demonstrating that the evolved distributions and predicted structure functions align well with experimental data from DESY-HERA and theoretical extractions while also providing predictions for future electron-ion collider kinematics.

The Gist

Imagine a pion not as a tiny, boring dot, but as a bustling, chaotic city. Inside this city live the "citizens" of the subatomic world: quarks (the main residents) and gluons (the glue holding them together).

For decades, physicists have been trying to map this city. They know the map for the "nucleon" cities (like protons and neutrons) is quite detailed. But the pion city? It's been a bit of a mystery. We know it exists, but we don't have a clear blueprint of how its citizens share the city's energy and momentum.

This paper is like a team of cartographers (Hari Govind P and his colleagues) who have built a new, highly detailed 3D map of the pion city using a special tool called the Light-Cone Quark Model.

Here is the breakdown of their journey, explained simply:

1. The Starting Point: The "Snapshot" Model

Imagine taking a high-speed photo of the pion city at a very specific, low-energy moment. In this photo, the city is very simple. It only has two main citizens: a quark and an antiquark (like a husband and wife). They are holding hands, and there are no extra guests (gluons or sea-quarks) yet.

The authors used a mathematical "camera" (the Light-Cone Quark Model) to calculate exactly how these two citizens share the city's total energy. They figured out that if you ask, "What percentage of the energy does the quark carry?" they found a specific answer that depends on how fast the quark is moving relative to the pion.

2. The Time Machine: Zooming into the Future (Evolution)

Here is the tricky part: The "snapshot" they took was at a low energy level. But real-world experiments (like those at the Large Hadron Collider or the upcoming Electron-Ion Collider) happen at huge energy levels.

If you zoom in on a city with a powerful microscope, you see more details. You see the main residents, but you also see the "glue" (gluons) vibrating and "sea-quarks" (temporary guests popping in and out) appearing.

To bridge the gap between their simple snapshot and the complex reality of high-energy experiments, the authors used a Time Machine called the DGLAP equations.

• The Analogy: Think of the pion as a balloon. At low pressure (low energy), it's small and simple. As you pump more air into it (increase the energy), it expands. As it expands, the "glue" and "sea-quarks" start to appear and take up space.
• The authors used this math to "inflate" their simple model, predicting how the quark, gluon, and sea-quark distributions change as the energy gets higher.

3. The Reality Check: Does the Map Match the Territory?

A map is useless if it doesn't match the real world. The authors took their "inflated" map and compared it against real data collected from massive experiments over the last 40 years (like FNAL, CERN, and HERA).

• The Result: Their map matched the real-world data surprisingly well!
• The Structure Function ($F_2$): They also predicted a specific "fingerprint" of the pion called the $F_2$ structure function. This is like a unique ID card that tells us how the pion interacts with light. They compared their predicted ID card with data from the HERA lab in Germany, and it fit perfectly.

4. The Big Test: The Drell-Yan Collision

To really prove their map was good, they simulated a crash test. They imagined smashing a pion into a heavy nucleus (like a tungsten target) to create a pair of muons (heavy electrons). This is called the Drell-Yan process.

They ran their simulation using their new pion map and compared the results to actual crash data from the COMPASS and FNAL experiments.

• The Verdict: The simulation matched the real crash data almost perfectly. This confirms that their "Light-Cone Quark Model" is a reliable way to understand the pion.

Why Does This Matter?

You might ask, "Why do we care about a pion?"

1. The Missing Piece: Pions are the lightest particles made of quarks, and they are the "glue" that holds atomic nuclei together. Understanding them is key to understanding why matter has mass.
2. Future Colliders: We are building new, massive machines called Electron-Ion Colliders (EIC). These machines will smash particles together with incredible force to see inside them. Before we turn these machines on, we need accurate maps (PDFs) to know what we are looking for. This paper provides a high-quality map for the pion, helping scientists interpret the data from these future experiments.
3. The "Sea" Surprise: They found that at very high energies, the "main residents" (valence quarks) only carry about 40% of the momentum. The rest is carried by the "glue" (gluons) and the "sea guests." This helps us understand how energy is distributed in the subatomic world.

In a Nutshell

The authors built a theoretical model of a pion, used math to predict how it changes at high speeds, and then proved their model works by showing it matches decades of real-world experimental data. It's like building a model of a hurricane, simulating how it grows, and then showing that your simulation perfectly predicts the wind speeds recorded in actual storms.

This work gives physicists a solid foundation to explore the deepest secrets of matter in the coming years.

Technical Summary

1. Problem Statement

Understanding the internal structure of hadrons, particularly mesons like the pion, remains a significant challenge in Quantum Chromodynamics (QCD). While Parton Distribution Functions (PDFs) for nucleons (baryons) are well-constrained by extensive experimental data, the partonic structure of mesons is less understood due to:

• Lack of Stable Targets: Pions are unstable, making direct experimental measurement of pion structure functions difficult.
• Theoretical Gaps: Direct calculations from the QCD Lagrangian are hindered by non-perturbative effects (confinement and chiral symmetry breaking).
• Data Scarcity: Existing experimental constraints on pion PDFs come from limited fixed-target Drell-Yan experiments (e.g., FNAL, CERN) and leading-neutron electroproduction at HERA, leaving significant uncertainties, especially at high energy scales.

The paper aims to address these gaps by calculating pion valence quark PDFs using a non-perturbative model, evolving them to experimental scales, and validating them against a wide range of existing and upcoming experimental data.

2. Methodology

The authors employ a multi-stage theoretical framework combining non-perturbative modeling with perturbative QCD evolution:

• Light-Cone Quark Model (LCQM):

  • The pion is treated as a bound state of a quark-antiquark pair ($q\bar{q}$) within the Light-Front (LF) framework.
  • The calculation focuses on the leading-order Fock state ($|q\bar{q}\rangle$), neglecting explicit gluon and sea-quark components at the initial scale.
  • Wave Functions: The total wave function is constructed as the product of a spin wave function ($S$) and a momentum-space radial wave function ($\phi$).
    • The radial wave function follows the Brodsky-Huang-Lepage (BHL) prescription.
    • The spin wave function is derived from the quark-meson vertex using proper Dirac spinors.
  • Parameters: The model uses constituent quark masses ($m_q = m_{\bar{q}} = 0.20$ GeV) and a harmonic scale parameter ($\beta = 0.410$ GeV), fitted to reproduce the pion mass and decay constant ($f_\pi \approx 122.9$ MeV).

• Derivation of Initial PDFs:

  • The unpolarized valence quark PDF, $f(x)$, is derived by solving the quark-quark correlation function using the overlap of the LCQM wave functions.
  • At the initial scale ($\mu_0$), the PDFs satisfy momentum sum rules where the valence quarks carry 100% of the momentum.

• QCD Evolution (DGLAP):

  • The initial valence PDFs are evolved from the low non-perturbative scale ($\mu_0$) to higher experimental scales ($\mu^2 = 5, 16, \dots, 1000$ GeV$^2$) using the Dokshitzer–Gribov–Lipatov–Altarelli–Parisi (DGLAP) equations.
  • Evolution is performed at Leading Order (LO), Next-to-Leading Order (NLO), and Next-to-Next-to-Leading Order (NNLO) using the HOPPET toolkit.
  • This evolution generates gluon and sea-quark distributions dynamically from the initial valence state.

• Phenomenological Applications:

  • Structure Function ($F_2^\pi$): Calculated at NLO accuracy and compared with HERA (ZEUS and H1) leading-neutron electroproduction data.
  • Drell-Yan Cross-Sections: Calculated for pion-nucleus collisions ($\pi A \to \ell^+\ell^- X$) using the evolved pion PDFs and nuclear PDFs (nCTEQ15 and nNNPDF20) from the LHAPDF library.

3. Key Contributions

• First NLO Prediction for $F_2^\pi$: The authors present, for the first time, a prediction for the pion structure function $F_2$ at Next-to-Leading Order (NLO) accuracy within the LCQM framework.
• Comprehensive Evolution: They provide a complete set of evolved valence, gluon, and sea-quark PDFs for the pion across a wide range of energy scales, including Mellin moments up to $n=4$.
• Extensive Experimental Validation: The model is rigorously tested against a vast array of historical and modern data:
  • FNAL-E-0615 and modified FNAL-E-0615 data.
  • CERN-WA-011, WA-039, NA-010, NA-003.
  • HERA (ZEUS/H1) leading-neutron data.
  • COMPASS-II (190 GeV) Drell-Yan data.
• Future Projections: The study includes predictions for the scale evolution of the pion structure function at the kinematics of the upcoming Electron-Ion Collider (EIC).

4. Key Results

• Initial Scale Determination: By fitting the evolved PDFs to modified FNAL-E-0615 data at $\mu^2 = 16$ GeV$^2$, the initial scale was determined to be $\mu_0 = 0.6 \pm 0.1$ GeV. The NNLO fit yielded the best $\chi^2$/d.o.f. (1.51).
• PDF Behavior:
  • The valence quark PDF $xf(x)$ peaks around $x > 0.5$ at the initial scale.
  • Upon evolution, the valence momentum fraction decreases as energy increases (e.g., dropping to ~39% at $\mu^2 = 49$ GeV$^2$), while gluon and sea-quark fractions increase, eventually dominating the momentum budget at high scales.
  • The calculated Mellin moments ($\langle x^n \rangle$) show excellent agreement with lattice QCD simulations and global fits (JAM, xFitter, MAP).
• Structure Function ($F_2^\pi$):
  • The NLO calculated $F_2^\pi$ agrees well with HERA data, particularly with the Effective One-Pion-Exchange (EF) flux model at high scales ($\mu^2 > 60$ GeV$^2$).
  • Sea-quark contributions are found to be significant, surpassing valence and gluon contributions in the low-$\beta$ region at high energies.
• Drell-Yan Cross-Sections:
  • The model successfully reproduces Drell-Yan cross-sections for various targets (Tungsten, Platinum, Carbon, Beryllium) and energies (39.5 GeV to 252 GeV).
  • The nCTEQ15 nuclear PDF set provided a superior description of the data compared to nNNPDF20, particularly for Tungsten targets.
  • The results show good agreement with recent COMPASS-II data, confirming the universality of the extracted pion PDFs.

5. Significance

• Validation of LCQM: The study demonstrates that the Light-Cone Quark Model, despite starting with a simple valence-only Fock state, can accurately describe complex hadronic structures when combined with perturbative QCD evolution.
• Constraint on Pion Structure: The work provides robust theoretical constraints on pion PDFs, filling a critical gap in the Standard Model's description of mesons.
• Preparation for Future Experiments: The predictions for the EIC and the detailed comparison with COMPASS-II data are vital for interpreting results from ongoing and future high-energy experiments.
• Foundation for Future Work: The authors identify that while the current model is successful, future refinements will require the inclusion of higher Fock states (explicit gluons and sea quarks) to further reduce discrepancies at very high $x$ and improve the description of gluon distributions.

In conclusion, this paper establishes a consistent and predictive framework for pion parton distributions, bridging the gap between non-perturbative quark models and high-energy experimental phenomenology.