The strange and antistrange quark distributions of the proton, s(x)s(x) and sˉ(x)\bar{s}(x), need not be equal: their difference is a purely nonperturbative signature of the proton's light-front Fock-state structure, it is constrained by the lattice QCD strange Dirac form factor F1s(Q2)F_1^s(Q^2), and it enters neutrino and parity-violating deep-inelastic observables. A 2018 analysis combining lattice QCD, light-front holographic QCD (LFHQCD), and the meson–baryon fluctuation model found a first moment ⟨S−⟩=0.0011(4)\langle S_-\rangle=0.0011(4) — far below the ∼0.005\sim 0.005 once invoked to explain the NuTeV weak-mixing-angle anomaly — and showed that the two frameworks make qualitatively different small-xx predictions. This problem asks for a full quantitative determination of the nonperturbative strange–antistrange asymmetry and the intrinsic strange probability: the magnitude and sign pattern of s(x)−sˉ(x)s(x)-\bar{s}(x) over all xx, the intrinsic strange probability IsI_s beyond its positivity lower bound, the separation of ss from sˉ\bar{s} beyond the current one-parameter degeneracy, the small-xx behavior that distinguishes the Regge-like holographic prediction from the vanishing meson-cloud one, and the factorization-scale matching to perturbative QCD evolution. The determination is still open: global fits remain essentially insensitive to the asymmetry at current precision, lattice input reaches only intermediate xx, and the small-xx region where the predictions diverge is unconstrained by any existing data.

Background

A parton distribution function (PDF) f(x)f(x) gives the probability density for finding a parton of flavor ff carrying a fraction xx of a hadron's longitudinal momentum in a fast frame. In the proton, the strange quark and antiquark distributions, s(x)s(x) and sˉ(x)\bar{s}(x), are sea distributions: they are not required by the valence uuduud content. Perturbative QCD generates sea pairs radiatively through gluon splitting (the extrinsic sea); in contrast, the intrinsic sea is the nonperturbative component already present in the proton's light-front Fock-state expansion — the quantum-mechanical superposition of states such as ∣uud⟩|uud\rangle and ∣uud ssˉ⟩|uud\,s\bar{s}\rangle, where the light-front wave function describes the bound state at fixed light-front time. Because the proton carries quark number, the strange quark in the ∣uud ssˉ⟩|uud\,s\bar{s}\rangle component interacts differently with the valence spectators than the antistrange quark does, so s(x)s(x) and sˉ(x)\bar{s}(x) need not coincide. Their difference, the strange–antistrange asymmetry s(x)−sˉ(x)s(x)-\bar{s}(x), is therefore a purely nonperturbative effect (the perturbative contribution first arises at three loops and is tiny).

Two nonperturbative frameworks make predictions. In the meson–baryon fluctuation model, the proton fluctuates part of the time into a virtual kaon plus a hyperon (KΛK\Lambda or KΣK\Sigma); the kaon's valence strange quark and the hyperon's antistrange spectator then acquire different momentum distributions. In light-front holographic QCD (LFHQCD) — a semiclassical approximation to QCD built by mapping light-front dynamics onto a gravitation-like theory in a higher-dimensional curved space — the strange form factor and the asymmetry follow from the Regge trajectory of the vector meson coupled to the strange current, where a Regge trajectory is the relation between a hadron's spin and its squared mass that governs its exchange behavior at high energy; its intercept controls how singular a PDF is as x→0x\to 0. A third, first-principles input is lattice QCD: the strange Dirac form factor F1s(Q2)F_1^s(Q^2), the strange-quark contribution to the proton's Dirac electromagnetic form factor at squared momentum transfer Q2Q^2, is nonzero at Q2>0Q^2>0 when the strange and antistrange distributions differ, s(x)−sˉ(x)≠0s(x)-\bar{s}(x)\neq 0, and its Q2Q^2 behavior constrains the first moment ⟨S−⟩=∫01dx x [s(x)−sˉ(x)]\langle S_-\rangle=\int_0^1 dx\, x\,[s(x)-\bar{s}(x)] through model fits to the lattice data. Lattice computations of F1sF_1^s are difficult because the strange contributions enter through disconnected quark loops with a poor signal-to-noise ratio, and the reachable Q2Q^2 is limited.

The original experimental motivation was the NuTeV anomaly: the NuTeV collaboration extracted the weak mixing angle sin⁡2θW\sin^2\theta_W from the Paschos–Wolfenstein ratio of neutral-current to charged-current neutrino and antineutrino deep-inelastic scattering cross sections on an iron target and found a value about 3σ3\sigma above the Standard Model prediction. Since ⟨S−⟩\langle S_-\rangle enters this ratio, a positive strangeness asymmetry of order ⟨S−⟩∼0.005\langle S_-\rangle\sim 0.005 was one candidate explanation. In the 2018 analysis that defines this problem, the meson–baryon fluctuation model and LFHQCD, normalized to the lattice strange form factor, both give ⟨S−⟩=0.0011(4)\langle S_- \rangle = 0.0011(4) — an order of magnitude too small — and they disagree qualitatively at small xx: LFHQCD predicts Regge behavior s(x)−sˉ(x)≃−0.044 x−0.01s(x)-\bar{s}(x)\simeq -0.044\,x^{-0.01} as x→0x\to 0, while the fluctuation model predicts the asymmetry vanishes. The holographic positivity bound gives a lower limit Is≥0.92%I_s \ge 0.92\% on the intrinsic strange probability Is=∫01dx s(x)=∫01dx sˉ(x)I_s=\int_0^1 dx\, s(x)=\int_0^1 dx\, \bar{s}(x), and the separation of s(x)s(x) from sˉ(x)\bar{s}(x) is fixed only up to a free parameter α\alpha constrained by the positivity of both distributions but not uniquely determined. Finally, comparing these hadronic-scale model distributions with the perturbatively evolved distributions used in global PDF fits requires choosing a factorization scale that neither framework specifies (the 2018 comparison assumed μ=1\mu = 1 GeV).

Problem Statement

Determine the nonperturbative strange–antistrange asymmetry of the nucleon quantitatively: obtain s(x)s(x) and sˉ(x)\bar{s}(x) — equivalently the asymmetry s(x)−sˉ(x)s(x)-\bar{s}(x) together with the intrinsic strange probability Is=∫01dx s(x)=∫01dx sˉ(x)I_s=\int_0^1 dx\, s(x)=\int_0^1 dx\, \bar{s}(x) — for the intrinsic (nonperturbative) strange sea of the proton, with uncertainties, including:

  1. The magnitude and sign of the first moment ⟨S−⟩=∫01dx x [s(x)−sˉ(x)]\langle S_-\rangle=\int_0^1 dx\, x\,[s(x)-\bar{s}(x)] and the sign pattern of s(x)−sˉ(x)s(x)-\bar{s}(x) across the full xx range (the lattice-implied pattern is negative at small xx and positive at large xx).
  2. The value of the intrinsic strange probability IsI_s as a determination rather than the positivity lower bound Is≥0.92%I_s\ge 0.92\%, including resolution of the one-parameter (α\alpha) degeneracy that currently fixes only the separation between s(x)s(x) and sˉ(x)\bar{s}(x).
  3. The small-xx behavior of s(x)−sˉ(x)s(x)-\bar{s}(x) — the Regge-like s(x)−sˉ(x)∼x−0.01s(x)-\bar{s}(x)\sim x^{-0.01} predicted by light-front holographic QCD versus the vanishing behavior of the meson–baryon fluctuation model.
  4. The factorization (hadronic) scale at which the nonperturbative model distributions are matched to perturbative QCD evolution, so that the comparison with global-fit PDFs is unambiguous.

These components jointly serve one overarching objective: a single, constraint-consistent determination of the nonperturbative strange sea, validated against lattice QCD input for the strange form factor F1s(Q2)F_1^s(Q^2) (with controlled continuum and physical-pion-mass systematics, extended to larger Q2Q^2 than presently available) and against experimental constraints such as neutrino dimuon production, W/ZW/Z production at hadron colliders, and parity-violating deep-inelastic scattering (scattering asymmetries generated by the interference of the electromagnetic and weak neutral currents, which are sensitive to individual sea-quark flavors). Resolving whether ⟨S−⟩\langle S_-\rangle can account for the NuTeV weak-mixing-angle discrepancy is part of the lineage of this problem but is not by itself a solution: the determination itself — in particular at small xx, where current data and lattice calculations are blind — is the target.

References

  1. R. S. Sufian, T. Liu, G. F. de Téramond, H. G. Dosch, S. J. Brodsky, A. Deur, M. T. Islam, B.-Q. Ma, 'Nonperturbative strange-quark sea from lattice QCD, light-front holography, and meson-baryon fluctuation models', Phys. Rev. D 98, 114004 (2018). DOI: 10.1103/PhysRevD.98.114004, https://arxiv.org/abs/1809.04975
  2. T.-J. Hou, H.-W. Lin, M. Yan, C.-P. Yuan, 'Impact of Lattice Strangeness Asymmetry Data in the CTEQ-TEA Global Analysis', Phys. Rev. D 107, 076018 (2023). DOI: 10.1103/PhysRevD.107.076018, https://arxiv.org/abs/2211.11064
  3. R. Zhang, H.-W. Lin, B. Yoon, 'Probing nucleon strange and charm distributions with lattice QCD', Phys. Rev. D 104, 094511 (2021). DOI: 10.1103/PhysRevD.104.094511, https://arxiv.org/abs/2005.01124
  4. D. Mason et al. (NuTeV Collaboration), 'Measurement of the Nucleon Strange-Antistrange Asymmetry at Next-to-Leading Order in QCD from NuTeV Dimuon Data', Phys. Rev. Lett. 99, 192001 (2007). DOI: 10.1103/PhysRevLett.99.192001, https://doi.org/10.1103/PhysRevLett.99.192001
  5. G. P. Zeller et al. (NuTeV Collaboration), 'Precise Determination of Electroweak Parameters in Neutrino-Nucleon Scattering', Phys. Rev. Lett. 88, 091802 (2002) [Erratum: Phys. Rev. Lett. 90, 239902 (2003)]. DOI: 10.1103/PhysRevLett.88.091802, https://arxiv.org/abs/hep-ex/0110059
  6. W. Bentz, I. C. Cloet, J. T. Londergan, A. W. Thomas, 'Reassessment of the NuTeV determination of the weak mixing angle', Phys. Lett. B 693, 462 (2010). DOI: 10.1016/j.physletb.2010.09.001, https://arxiv.org/abs/0908.3198
  7. C. Alexandrou, S. Bacchio, M. Constantinou, J. Finkenrath, K. Hadjiyiannakou, K. Jansen, G. Koutsou, 'Nucleon strange electromagnetic form factors', Phys. Rev. D 101, 031501 (2020). DOI: 10.1103/PhysRevD.101.031501, https://arxiv.org/abs/1909.10744
  8. C. Alexandrou, S. Bacchio, M. Bode, J. Finkenrath, A. Herten, C. Iona, G. Koutsou, F. Pittler, B. Prasad, G. Spanoudes, 'Nucleon strange electromagnetic form factors using Nf=2+1+1N_f=2+1+1 twisted-mass fermions at the physical point', arXiv:2605.02049 (2026). https://arxiv.org/abs/2605.02049
  9. R. S. Sufian, T. Liu, A. Alexandru, S. J. Brodsky, G. F. de Téramond, H. G. Dosch, T. Draper, K.-F. Liu, Y.-B. Yang, 'Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD', Phys. Lett. B 808, 135633 (2020). DOI: 10.1016/j.physletb.2020.135633, https://arxiv.org/abs/2003.01078
  10. X. G. Wang, A. W. Thomas, 'Challenges in the extraction of physics beyond the Standard Model from electron scattering', J. Phys. G 52, 015006 (2025). DOI: 10.1088/1361-6471/ad91d1, https://arxiv.org/abs/2403.07327
  11. R. M. Whitehill, M. M. Dalton, T. Liu, W. Melnitchouk, N. Sato, 'Impact of parity-violating deep-inelastic scattering on the weak mixing angle and high-xx parton distributions', arXiv:2606.26056 (2026). https://arxiv.org/abs/2606.26056
  12. J. M. Cruz-Martinez, M. Fieg, T. Giani, P. Krack, T. Mäkelä, T. Rabemananjara, J. Rojo, 'The LHC as a Neutrino-Ion Collider', arXiv:2309.09581 (2023). https://arxiv.org/abs/2309.09581
  13. Y. Duan, S. Xu, S. Cheng, X. Zhao, Y. Li, J. P. Vary, 'Flavor asymmetry from the nonperturbative nucleon sea', Phys. Rev. C 110, 065201 (2024). DOI: 10.1103/PhysRevC.110.065201, https://arxiv.org/abs/2404.07755