Isospin dependence of nuclear EMC effect from global QCD analysis

This paper presents a global QCD analysis incorporating recent MARATHON experiment data to simultaneously determine parton distribution functions and nucleon off-shell corrections, revealing that both isoscalar and isovector nuclear effects are required to describe light nuclei ($A \leq 3$) and that their EMC ratios differ significantly from naive extrapolations of heavy nuclei.

The Gist

Imagine you have a set of building blocks. In the world of physics, these blocks are protons and neutrons (collectively called nucleons), which make up the nucleus of every atom. Inside these nucleons, there are even smaller, faster-moving particles called quarks.

For decades, physicists have known that when you pack these nucleons tightly together inside a heavy nucleus (like gold or lead), the behavior of the quarks inside changes. It's as if the quarks get "squished" or "stretched" by their neighbors. This phenomenon is called the EMC effect (named after the European Muon Collaboration who discovered it).

However, there's a big mystery: Why does this happen? And does it happen the same way for every type of nucleus?

This paper, written by a team of physicists from the US and Australia, tackles this mystery by looking at the lightest, simplest nuclei (like Deuterium, Helium-3, and Tritium) instead of the heavy ones. They used a massive amount of data from the MARATHON experiment at Jefferson Lab to solve a puzzle that previous models couldn't quite fit.

Here is the story of their discovery, broken down into simple concepts:

1. The "Recipe" Problem

Think of trying to understand a cake by looking at the ingredients.

• The Old Way: Previous scientists tried to understand the quarks in a nucleus by assuming the "recipe" (the physics rules) was the same for a tiny cake (Deuterium) as it was for a giant wedding cake (Lead). They used a specific mathematical model (called the KP model) to guess how the ingredients should behave in the tiny cakes because they hadn't measured them directly yet.
• The New Way: The authors of this paper said, "Let's stop guessing." They took all the available data from the world's best experiments and let the data speak for itself. They used a super-computer method (Bayesian Monte Carlo) to simultaneously figure out the "recipe" for the quarks and how the nucleus changes that recipe.

2. The "Off-Shell" Glitch

In physics, a "free" proton is like a dancer spinning perfectly on a stage. But inside a nucleus, the dancer is bumping into others, getting pushed, and losing a bit of energy. This state is called being "off-shell."

• The Analogy: Imagine a group of friends playing tag in a small room. If you are running freely in a big park (a free proton), you move one way. But in a crowded room (a nucleus), you have to dodge people. Your movement changes not just because of the room, but because you are physically interacting with the walls and other people.
• The Discovery: The paper proves that to understand the data from the light nuclei (Helium-3 and Tritium), you must account for this "crowded room" effect. If you ignore it, your predictions are wildly wrong.

3. The "Isospin" Twist (The Left-Handed vs. Right-Handed Effect)

This is the most exciting part. Physicists have two types of quarks in protons and neutrons: Up and Down.

• A Proton has two Ups and one Down.
• A Neutron has one Up and two Downs.

The big question was: Does the nucleus treat the "Up" quarks differently than the "Down" quarks?

• The Old Assumption: Many thought the effect was the same for both (like a generic "crowded room" effect).
• The New Finding: The authors found that the nucleus treats them differently. It's like the crowded room pushes the "Left-handed" dancers one way and the "Right-handed" dancers another way.
  • In Helium-3 (which has two protons and one neutron), the effect is strong and specific.
  • In Tritium (one proton and two neutrons), the effect is different.
  • This proves there is an "Isospin Dependence." The nuclear environment cares about which type of quark it is.

4. Why the Previous Models Were Wrong

The paper points out a major flaw in how the MARATHON experiment's data was previously analyzed.

• The Bias: The previous analysis forced the data to fit the "KP Model" (the giant wedding cake recipe). They adjusted the numbers (normalization) so the tiny cakes looked like they fit the big cake's rules.
• The Result: This was like forcing a square peg into a round hole. It made the data look consistent with the old theory, but it hid the truth.
• The Correction: When this team removed that "forced fit" and let the data breathe, they found the numbers didn't match the old model at all. The "crowded room" effect in light nuclei is actually quite different from what we thought.

The Bottom Line

This paper is a "reset button" for our understanding of atomic nuclei.

1. We need new rules: The way quarks behave in light nuclei (like Helium and Tritium) is not just a smaller version of heavy nuclei. It has its own unique personality.
2. It's not just "squishing": The nucleus changes the behavior of "Up" quarks differently than "Down" quarks.
3. Stop guessing: By using a global analysis that doesn't assume the answer beforehand, they found that the "off-shell" corrections (the crowded room effect) are essential and have a specific, complex shape.

In everyday terms: Imagine you thought all cars drove the same way on a highway. You assumed a tiny go-kart (light nucleus) would react to traffic exactly like a massive semi-truck (heavy nucleus). This paper says, "No! The go-kart reacts differently because it's lighter and has different steering." To understand the physics of the universe, we have to stop treating the tiny and the massive as if they are identical twins. They are cousins, and they have different personalities.

Technical Summary

1. Problem Statement

The paper addresses the long-standing challenge of understanding the nuclear EMC effect—the modification of parton distribution functions (PDFs) when nucleons are bound inside a nucleus compared to free nucleons. While the effect is well-established for heavy nuclei, its origin and specific dependence on isospin (the difference between proton and neutron modifications) remain debated.

Key specific problems addressed include:

• Light Nuclei Constraints: Previous analyses of light nuclei ($A \le 3$, specifically Deuterium, $^3$He, and $^3$H) often relied on the Kulagin-Petti (KP) model, which assumes a universal shape for nuclear corrections across all nuclei (from Pb down to Deuterium). This assumption forces nuclear effects to vanish at a specific Bjorken-$x$ ($x_B \approx 0.31$) and assumes the same correction magnitude for light and heavy nuclei, despite vastly different binding energies.
• Bias in MARATHON Data: Recent high-precision data from the MARATHON experiment at Jefferson Lab (measuring $^3$He/D and $^3$H/D cross-section ratios) were previously normalized to match the KP model. The authors argue this introduces a significant theoretical bias, potentially obscuring the true isospin dependence of the nuclear effects.
• Off-Shell Ambiguity: It is unclear whether nucleon off-shell corrections (modifications due to the nucleon being off its mass shell inside the nucleus) are purely isoscalar (affecting protons and neutrons equally) or if they contain an isovector component (differing between protons and neutrons).

2. Methodology

The authors perform a global QCD analysis using the JAM (Jefferson Lab Angular Momentum) Bayesian Monte Carlo framework.

• Data Set: The analysis incorporates world high-energy scattering data, including:
  • Proton and Deuteron DIS data.
  • New MARATHON data: Individual $^3$He/D and $^3$H/D structure function ratios (replacing the older $^3$He/$^3$H ratio).
  • Other fixed-target, HERA, Drell-Yan, and collider data.
• Theoretical Framework:
  • Convolution Formalism: Nuclear structure functions are modeled as convolutions of free nucleon structure functions with smearing functions (accounting for Fermi motion and binding) and nucleon off-shell corrections.
  • Off-Shell Decomposition: The off-shell corrections ($\delta q_{N/A}$) are decomposed into isoscalar ($\delta q^{(0)}$) and isovector ($\delta q^{(1)}$) components.
    • Isoscalar: Affects $u$ and $d$ quarks similarly in all nuclei.
    • Isovector: Depends on the isospin of the spectator system (e.g., different for $^3$He vs. $^3$H).
  • Parameterization:
    • Free nucleon PDFs are parametrized at the input scale $\mu_0^2 = m_c^2$.
    • Off-shell functions are parametrized independently for isoscalar ($\delta q_0$) and isovector ($\delta u_1, \delta d_1$) components.
    • Higher-twist (HT) corrections are treated using an additive parametrization (unlike the multiplicative form in some previous works) to better separate $1/Q^2$ power corrections from leading-twist logarithmic evolution.
• Fitting Strategy:
  • A simultaneous fit of free nucleon PDFs and off-shell parameters is performed.
  • Crucial Difference: Unlike previous analyses, no model-dependent normalization (like the KP model) is imposed on the MARATHON data. Instead, normalization factors are fitted as free parameters with Gaussian penalties based on experimental uncertainties.
  • Two scenarios are tested: (i) only isoscalar corrections, and (ii) symmetric isoscalar + asymmetric isovector corrections.

3. Key Contributions

1. Removal of Model Bias: The study is the first to analyze MARATHON $^3$He/D and $^3$H/D data without forcing them to align with the KP model. This allows the data to dictate the normalization and the nature of nuclear corrections.
2. Evidence for Isospin Dependence: By allowing both isoscalar and isovector off-shell terms, the analysis provides the first robust global QCD evidence for an isovector component in the nuclear EMC effect for light nuclei.
3. Refutation of Universal KP Assumptions: The analysis demonstrates that the assumption of a universal nuclear correction shape (vanishing at $x_B=0.31$) is inconsistent with the data when treated without bias.
4. Quantification of Off-Shell Effects: The study quantifies the magnitude of off-shell corrections for $A=2$ and $A=3$ systems, showing they are distinct from heavy nuclei.

4. Key Results

• Necessity of Off-Shell Corrections:
  • For the $^3$He/D ratio, a fit without off-shell corrections yields a poor reduced $\chi^2$ ($\chi^2_{red} \approx 5.10$).
  • Including off-shell corrections drastically improves the fit ($\chi^2_{red} \approx 1.04$), providing strong evidence that off-shell effects are essential to describe DIS on light nuclei.
• Isoscalar vs. Isovector:
  • The isoscalar off-shell correction is found to be non-zero, positive, and peaks around $x \approx 0.4$. It is negative in its contribution to the structure function ratio (due to the convolution with negative smearing functions).
  • The isovector correction is small but non-zero. Specifically, the difference $\delta u_1 - \delta d_1$ shows a negative trend at large $x$, indicating a flavor-dependent modification.
  • Scenario Preference: The data slightly favor the scenario with both isoscalar and asymmetric isovector corrections over the purely isoscalar scenario, particularly in describing the high-$x$ region of the $^3$H/D data.
• Normalization Discrepancies:
  • The fitted normalization for $^3$He/D is $0.996(6)$, which is consistent with unity but differs by 3.2$\sigma$ from the KP-model-imposed normalization of $1.021(5)$.
  • This confirms that the previous upward shift applied to MARATHON data to match the KP model was unjustified and introduced a bias.
• EMC Ratios:
  • The extracted EMC ratios ($R_D, R_{^3He}, R_{^3H}$) do not cross unity at $x_B = 0.31$ as assumed by the KP model.
  • At $x_B = 0.31$, the JAM results differ from the KP predictions by $2.3\sigma$ to $3.0\sigma$.
  • The $^3$He/$^3$H super-ratio tends toward values below 1 at large $x$, whereas the KP model predicts values above 1.
• Impact on PDFs:
  • The inclusion of MARATHON data has a minimal impact on the central values of free nucleon PDFs (which are well-constrained by other data) but significantly reduces uncertainties in the $d$-quark PDF at large $x$ ($x \gtrsim 0.6$).

5. Significance

• Dynamical Origin of EMC Effect: The results support dynamical mechanisms involving short-range correlations and mean fields that generate both isoscalar and isovector nuclear effects, challenging models that assume a universal, flavor-independent modification.
• Methodological Shift: The paper establishes a new standard for analyzing light-nucleus data, demonstrating that relying on heavy-nucleus extrapolations (like the KP model) for light nuclei leads to biased conclusions.
• Future Directions: The findings suggest that future experiments (e.g., BONuS12 at Jefferson Lab) focusing on spectator tagging and parity-violating DIS are critical for further constraining the isospin dependence of the EMC effect and the $d/u$ ratio at large $x$.

In summary, this work provides a rigorous, model-independent global QCD analysis that confirms the necessity of nucleon off-shell corrections in light nuclei and reveals a distinct isospin dependence in the nuclear EMC effect, fundamentally challenging previous assumptions derived from the Kulagin-Petti model.