An analysis of nuclear parton distribution function based on Kullback-Leibler divergence

This paper proposes a Kullback-Leibler divergence-based methodology to quantify differences between nuclear and free nucleon parton distribution functions, successfully modeling the EMC effect for quarks and revealing that the EPPS21 gluon nPDFs align more closely with the "minimum relative entropy" hypothesis than nNNPDF3.0, thereby offering a novel tool for future global fits in data-scarce regimes.

The Gist

Here is an explanation of the paper using simple language, everyday analogies, and creative metaphors.

The Big Picture: The "Crowded Room" Problem

Imagine a proton (a tiny particle inside an atom) as a solo musician playing a guitar. When the musician is alone in a quiet room, we know exactly how they play: the notes they hit, the rhythm, and the volume. In physics, we call this the "Parton Distribution Function" (PDF). It's a map showing where the energy and momentum are inside the particle.

Now, imagine putting that musician into a packed, noisy concert hall full of other musicians (this is an atomic nucleus). The music changes! The soloist has to shout over the crowd, adjust their rhythm to avoid bumping into others, and maybe even change the song entirely.

In physics, when protons are stuck inside a nucleus, their "music" (their internal structure) changes. We call this the EMC effect. Scientists have been trying to map exactly how the music changes for decades. They usually do this by measuring the concert hall with microphones (experiments) and then trying to guess the rules of the noise.

The New Tool: The "Information Thermometer"

The authors of this paper decided to try a completely new way to solve this puzzle. Instead of just looking at the raw data, they used a concept from Quantum Information Theory called Kullback-Leibler (KL) Divergence.

The Analogy:
Think of KL Divergence as an "Information Thermometer" or a "Confusion Meter."

• It measures how much "surprise" or "extra information" you need to explain one thing (the noisy concert) when you are already expecting another thing (the solo performance).
• If the soloist and the concert hall sound exactly the same, the "Confusion Meter" reads zero.
• If the concert hall sounds totally different, the meter reads high.

The authors used this meter to measure the difference between a free proton and a proton inside a nucleus.

The "Minimum Energy" Hypothesis

Here is the clever part. The authors proposed a rule they call the "Minimum Relative Entropy" hypothesis.

The Analogy:
Imagine you are trying to figure out the shape of a muddy footprint in the snow. You know the heel and the toe, but the middle is covered in snow.

• Old way: Guess randomly based on what other footprints look like.
• This paper's way: Nature is lazy. It always takes the path of least resistance. The authors hypothesize that when a proton gets squeezed into a nucleus, it rearranges its internal parts in the most efficient, "lazy" way possible to minimize the "Confusion Meter" (KL Divergence).

It's like a ball rolling down a hill; it will naturally stop at the very bottom (the minimum point). The authors assume the proton's structure does the same thing: it settles into the shape that creates the least amount of "informational friction" with its surroundings.

What They Found

1. For Quarks (The "Strings" of the Guitar):
   When they applied this "lazy ball" rule to the quarks (the main parts of the proton), the shape they calculated matched perfectly with the best existing maps made by other scientists.

   • Translation: Their new thermometer works! It confirms that nature really does seem to follow this "minimum confusion" rule.

2. For Gluons (The "Glue" holding it together):
   This is where it gets exciting. Gluons are the particles that stick quarks together. Scientists are very confused about how gluons behave in a nucleus because the data is messy and different groups of scientists disagree on the map.

   • The authors used their "Minimum Entropy" rule to predict what the gluon map should look like.
   • They compared their prediction to two famous existing maps: EPPS21 and nNNPDF3.0.
   • The Result: Their prediction matched the EPPS21 map much better than the nNNPDF3.0 map.

Why This Matters

Think of this paper as a new quality control test for nuclear physics.

• Before, if two scientists had different maps of the nucleus, there was no easy way to say which one was "more right" without doing more expensive experiments.
• Now, the authors say: "Let's check the 'Confusion Meter.' The map that requires the least amount of extra information to explain the change from a free proton to a nuclear proton is likely the correct one."

The Takeaway

This paper suggests that the universe has a hidden "efficiency principle" for how particles behave when they are crowded together. By using a tool from information theory (KL Divergence), the authors found a new way to predict the structure of atomic nuclei.

It's like realizing that even in a chaotic, noisy concert hall, the musicians are actually following a very simple, efficient rule to stay in tune. This discovery gives physicists a powerful new compass to navigate the messy world of nuclear gluons, potentially leading to better maps of the universe's building blocks.

Technical Summary

Here is a detailed technical summary of the paper "An analysis of nuclear parton distribution function based on Kullback-Leibler divergence."

1. Problem Statement

The paper addresses the challenge of understanding Nuclear Parton Distribution Functions (nPDFs), which describe how partons (quarks and gluons) are distributed within nucleons bound inside an atomic nucleus.

• The EMC Effect: A primary motivation is the "EMC effect," where quark distributions in bound nucleons deviate from free nucleons in the intermediate-$x$ region ($0.25 \lesssim x \lesssim 0.65$).
• Theoretical Limitations: Non-perturbative Quantum Chromodynamics (QCD) cannot analytically solve for nPDFs, and Lattice QCD is currently computationally infeasible for complex nuclei. Consequently, current nPDFs rely entirely on global fits to experimental data.
• Gluon Uncertainty: While quark nPDFs are relatively well-constrained by data, gluon nPDFs remain highly uncertain. Different global fitting groups (e.g., EPPS21, nNNPDF3.0) produce significantly divergent predictions for gluons, particularly in the intermediate-$x$ region.
• Goal: The authors propose a new, model-independent methodology to quantify nuclear modifications and determine the shape of nPDFs using principles from Quantum Information Theory, specifically the Kullback-Leibler (KL) divergence.

2. Methodology

The authors introduce an information-theoretic framework to treat the transformation of a free nucleon's Parton Distribution Function (PDF) into a nuclear nPDF as an optimization problem.

• KL Divergence (Relative Entropy):
  The core metric is the KL divergence, defined as:
  $$D_{KL}(p \parallel q) = \int p(x) \ln \frac{p(x)}{q(x)} dx$$
  Here, $q(x)$ represents the reference distribution (the PDF of a free nucleon, approximated by the deuteron using CT18A free PDFs), and $p(x)$ represents the true distribution (the nPDF of a nucleus $A$).

• The "Minimum Relative Entropy" Hypothesis:
  The paper hypothesizes that when a free nucleon is bound in a nucleus, its PDF transforms into an nPDF in a way that minimizes the relative entropy (KL divergence) between the two, subject to physical constraints (specifically, fixed boundary conditions at the endpoints of the intermediate-$x$ region).

  • This is analogous to the "principle of least action" or the "brachistochrone problem" in classical mechanics, where a system evolves to minimize a specific functional.

• Implementation:

  1. Normalization: Since numerical integration is performed over a restricted range (e.g., $x \in [0.25, 0.65]$ for quarks), the distributions are renormalized within this range to ensure the KL divergence remains non-negative and physically meaningful.
  2. Parametrization: To find the function $\bar{p}_A(x)$ that minimizes the KL divergence, the authors test two functional forms:
     • Polynomial: $\bar{p}_A(x) = (a_1 x^3 + b_1 x^2 + c_1 x + d_1) q_d(x)$
     • Canonical: $\bar{p}_A(x) = (a_2 x^{b_2} (1-x)^{c_2} e^{d_2 x}) q_d(x)$
  3. Constraints: The coefficients are determined by:
     • Fixing the values at the endpoints ($x_{min}$ and $x_{max}$) based on existing global fit data (EPPS21 or nNNPDF3.0).
     • Minimizing the KL divergence with respect to the remaining free parameters.

3. Key Contributions

• Novel Metric for Nuclear Structure: The paper successfully applies KL divergence as a quantitative measure of nuclear modification, moving beyond simple ratio plots to an information-theoretic distance metric.
• The Minimum Relative Entropy Principle: It proposes a new physical principle for determining nPDF shapes in data-scarce regions, framing the nuclear modification as an entropy-minimization process.
• Benchmarking Tool: The methodology provides a rigorous, independent benchmark to evaluate the quality of global fitting groups. If a global fit's result is close to the "minimum entropy" solution, it suggests the fit is physically consistent with the hypothesis.
• Robustness Check: The authors demonstrate that the results are robust across different parametrization forms (polynomial vs. canonical), indicating the solution is driven by the entropy principle rather than the specific mathematical ansatz.

4. Results

A. Quark nPDFs (Validation)

• The authors applied the method to quark nPDFs in the intermediate-$x$ region ($x \in [0.25, 0.65]$) for nuclei including $^4$He, $^{12}$C, $^{56}$Fe, and $^{208}$Pb.
• Agreement with Data: The structure functions derived from the "minimum relative entropy" hypothesis showed excellent agreement with the latest global fits (EPPS21).
• Quantitative Metrics: The $L_2$ norm (distance) between the entropy-minimized curves and the EPPS21 data was very small ($10^{-3}$to$10^{-4}$), whereas reference curves with non-minimal entropy had much larger deviations.
• Conclusion: This validates the hypothesis for quarks, suggesting that the EMC effect can be effectively described by minimizing relative entropy.

B. Gluon nPDFs (Prediction and Comparison)

The method was then applied to gluon nPDFs, where experimental constraints are weak. The authors compared the entropy-minimized results against two major global fitting groups:

1. Comparison with EPPS21:

   • The gluon nPDFs derived from the minimum relative entropy hypothesis closely aligned with the EPPS21 results.
   • The $L_2$ norms were small, and the KL divergence values were minimized.
   • Implication: EPPS21 appears to be more consistent with the "minimum relative entropy" physical principle.

2. Comparison with nNNPDF3.0:

   • The entropy-minimized results showed significant deviations from the nNNPDF3.0 predictions, particularly for heavier nuclei ($^{56}$Fe, $^{208}$Pb).
   • The $L_2$ norms were approximately four times larger than those observed for EPPS21.
   • Implication: The nNNPDF3.0 gluon parametrization may not satisfy the minimum relative entropy condition as well as EPPS21, suggesting potential areas for refinement in the nNNPDF3.0 global fit.

5. Significance and Conclusion

• New Perspective on QCD: The work bridges high-energy physics and quantum information theory, suggesting that the complex, non-perturbative dynamics of nuclear binding may be governed by information-theoretic principles (entropy minimization).
• Guidance for Future Fits: The proposed methodology offers a powerful tool for future global fitting efforts. It suggests that fitting groups should aim for solutions that minimize the KL divergence relative to free nucleons, providing a theoretical constraint where experimental data is lacking.
• Resolution of Gluon Uncertainty: By identifying that EPPS21 aligns better with the entropy hypothesis than nNNPDF3.0, the paper provides a theoretical criterion to resolve discrepancies in gluon nPDFs, which are crucial for precision physics at the LHC and future colliders.
• Future Work: The authors note that extending this method to the full $x$-range and incorporating more complex constraints are necessary next steps.

In summary, the paper establishes that the Kullback-Leibler divergence and the minimum relative entropy hypothesis provide a robust, model-independent framework for quantifying nuclear modifications and offer a new, theoretically grounded benchmark for determining the elusive gluon nPDFs.