Pith. sign in

REVIEW 2 major objections 5 minor 17 references

Nuclear suppression in diffractive vector meson production within the color glass condensate framework

T0 review · 2 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Within the color glass condensate, coherent J/psi photoproduction on heavy nuclei is suppressed far below naive A^(4/3) scaling, down to about a factor of 0.15 at high energy, with the suppression growing with both nuclear size and…

desk verdict A credible, honest model study, but the suppression numbers are not out-of-sample predictions because the K-factor and the data used to fit it shape the A-dependence. read the letter →

arxiv 2508.21562 v1 pith:5PBTDHAM submitted 2025-08-29 hep-ph nucl-th

Heikki Mäntysaari , Hendrik Roch , Farid Salazar , Björn Schenke , Chun Shen , Wenbin Zhao This is my paper · ORCID
classification hep-phnucl-th
keywords colorglasscondensatediffractivevectormesonproductioncoherentJ/psiphotoproductionnuclearsuppressiongluonsaturationBayesianparameterestimationultraperipheralheavy-ioncollisionssmall-xphysics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper seeks to establish that, within a Color Glass Condensate framework, coherent $\mathrm{J}/\psi$ photoproduction on heavy nuclei is strongly suppressed by gluon saturation once proton and nuclear data are fitted together. It reports that at photon–nucleus energies $W=31.5$ GeV and $W=813$ GeV, the cross section ratio for heavy nuclei relative to the naive no-saturation $A^{4/3}$ scaling falls to about $0.3$ and $0.15$, respectively. A global $K$-factor of about $0.33$ allows the same framework to describe both $\gamma+p$ and $\gamma+\mathrm{Pb}$ data, and because this rescaling is compensated by a larger color charge density, it strengthens the predicted nuclear suppression. If these predictions hold, exclusive $\mathrm{J}/\psi$ measurements at a future electron–ion collider would provide a direct probe of the onset of gluon saturation.

What carries the argument

The machinery is the coherent diffractive cross section formula $d\sigma_{\gamma+A\to \mathrm{J}/\psi+A}/dt = \frac{K}{16\pi}|\langle\mathcal{A}_{\gamma^*+p\to V+p}\rangle_\Omega|^2$, averaged over target configurations. The amplitude uses the photon–vector-meson wave function overlap and the dipole–target amplitude derived from Wilson lines in the McLerran–Venugopalan model, with the Wilson lines evolved event-by-event using the JIMWLK equation. The load-bearing parameter is the global $K$-factor, about $0.33$, which absorbs wave-function and higher-order uncertainties; the fit compensates $K<1$ with a larger color charge density, and that compensation is what drives the strong $A$-dependent suppression.

What would settle it

Measure the ratio of the coherent $\mathrm{J}/\psi$ cross section on a mid-size and a heavy nucleus to that on the proton at $W\approx 813$ GeV: the paper predicts a monotonic fall of this ratio to the $A^{4/3}$ baseline, reaching about $0.15$ for the heaviest nuclei, and a weaker fall at $W=31.5$ GeV. An experimental ratio consistent with $1$ at both energies would falsify the central claim.

Watch

Extended reading notes

Core claim

The central claim is that a CGC calculation whose parameters come from a Bayesian fit to HERA and LHC data produces a clear departure from the $A^{4/3}$ scaling expected without saturation. The coherent $\mathrm{J}/\psi$ cross section on nuclei, normalized to the proton case, is suppressed to roughly $0.3$ at $W=31.5$ GeV and to roughly $0.15$ at $W=813$ GeV for the heaviest nuclei, with the suppression growing monotonically with nuclear mass number $A$ and with energy. The authors attribute the improved simultaneous description of proton and lead data to an overall $K$-factor of about $0.33$, whose effect is offset in the fit by a larger color charge density, making the nucleons denser and enhancing the nuclear suppression.

Load-bearing premise

The prediction rests on the unsuppressed baseline being $A^{4/3}$ and on the fitted $K$-factor being the same for protons and all nuclei at all collision energies; if nuclear form factors shift the baseline or $K$ depends on $A$ or $W$, the reported suppression values change.

Editorial extensions

If this is right

  • Coherent $\mathrm{J}/\psi$ production on heavy nuclei should visibly violate the $A^{4/3}$ scaling at both $W=31.5$ and $813$ GeV, with about a factor of two more suppression at the higher energy for large $A$.
  • The global $K$-factor of about $0.33$ means the leading-order CGC cross section needs a substantial overall rescaling, so any parameter inference from this framework depends on the universality of $K$.
  • Future electron–ion collider measurements of exclusive $\mathrm{J}/\psi$ production would give a direct, $Q^2$-resolved test of where gluon saturation sets in, complementing the two-energy predictions made here.
  • Extending the same analysis to other quarkonium states or higher $Q^2$ would show whether the suppression is a generic saturation feature or particular to the $\mathrm{J}/\psi$ wave function.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the $K$-factor turns out to depend on $A$ or $W$, the quoted suppression factors would likely shift, and the authors' own caveat that nuclear form factors must be included implies the present numbers are best read as an upper bound on the saturation-driven suppression.
  • A measurement spanning several nuclei at one fixed energy in the same detector would separate the trivial geometric $A^{4/3}$ factor from genuine saturation effects more cleanly than the two-point comparison in this paper.
  • The large deviation of $K$ from $1$ hints at missing next-to-leading-order contributions or wave-function normalization; if those were computed, the fitted color charge density would probably decrease and with it the predicted suppression.
  • The framework implies that raising photon virtuality $Q^2$ at the EIC should weaken the suppression, which would help distinguish saturation from nuclear shadowing.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. This manuscript studies the A-dependence of coherent diffractive J/psi photoproduction in a Color Glass Condensate framework with JIMWLK evolution. The authors use the posterior distribution from a companion global Bayesian analysis (Ref. [13]) that fits model parameters, including an overall K-factor of about 0.33, to HERA gamma-p and LHC gamma-Pb data. From 25 posterior samples they compute the ratio of the nuclear to proton cross section as a function of A at W = 31.5 GeV (x = 0.01) and W = 813 GeV (x = 1.5e-5). They report a strong suppression relative to the A^{4/3} scaling expected without saturation, with the ratio reduced to about 0.3 at W = 31.5 GeV and about 0.15 at W = 813 GeV for the heaviest nuclei, and they interpret this as evidence for gluon saturation. The paper concludes by noting that a more quantitative assessment would require including nuclear form factors.

Significance. If the quantitative suppression factors were robust, this would be a useful demonstration of the sensitivity of coherent J/psi production to nuclear gluon saturation and would provide a target for EIC measurements. The manuscript benefits from using a modern CGC framework with a Bayesian uncertainty estimate, and it explicitly propagates parameter uncertainties from a global fit. However, the significance is substantially tempered by two issues. First, the parameters, including the K-factor and color charge density, are fitted to the very gamma-Pb data that the suppression is meant to characterize, so the quoted suppression is a posterior statement rather than an independent prediction. Second, the A^{4/3} baseline is an approximation that neglects nuclear form factors, a limitation the authors themselves acknowledge. As a result, the central numerical claims are not isolated from normalization degeneracies and finite-nuclear-size effects.

major comments (2)
  1. [Section 3, Eq. (1) and surrounding text] The central suppression factors (~0.3 at W = 31.5 GeV and ~0.15 at W = 813 GeV) are not independent predictions: the model parameters, including the K-factor and the color charge density, are fitted to the same gamma+Pb data that the comparison is intended to explain. The paper explicitly states that 'A value K < 1 is compensated in the fits by a larger color charge density, which corresponds to denser nucleons and, in turn, stronger nuclear suppression.' This compensation means the A-dependent result is shaped by the fit's normalization degeneracy, not purely by CGC dynamics. The authors should show a comparison with a fit where K is fixed to unity, or with the Pb data removed from the fit, to demonstrate that the reported A-dependence is robust. Without such a control, the claim that the results show a 'clear departure from A^{4/3} scaling' is overstated.
  2. [Section 4, Fig. 1] The A^{4/3} no-saturation baseline is an approximate scaling derived for a uniform nuclear density and does not account for nuclear form factors, finite nuclear size, or the t-integration. The authors themselves state in Section 4 that 'a more quantitative assessment of the actual saturation effect will require taking into account the nuclear form factors.' Therefore the quoted suppression factors conflate genuine saturation effects with nuclear-geometry effects. The no-saturation reference should be computed within the same model, for example by switching off the JIMWLK evolution or linearizing the dipole amplitude, rather than using a phenomenological A^{4/3} scaling. Until that is done, the quantitative suppression numbers cannot be attributed uniquely to saturation.
minor comments (5)
  1. [Section 2] The manuscript uses only 25 posterior samples from Ref. [13] to represent the full posterior and to define the 2-sigma error bands in Fig. 1. Given the likely degeneracy between K and the color charge density, 25 samples is small; the authors should report the effective sample size or demonstrate that the quoted central values and error bands are stable with respect to the number of samples.
  2. [Section 2] The values W = 31.5 and 813 GeV are said to correspond to x = 0.01 and 1.5e-5, but the scale Q^2 at which x is evaluated is not specified. Please state the relevant scale or the kinematics used to relate W and x.
  3. [Section 2] There is a typo in 'McLerran-Venguopalan'; the correct name is McLerran-Venugopalan.
  4. [Figure 1 caption] The caption does not explain the content of the two panels shown in the figure (upper panel with sigma/sigma_p and the A^{4/3} dashed line, and lower panel with the ratio to A^{4/3}). Please describe both panels explicitly so the reader can interpret the figure without the main text.
  5. [Section 2] The statement that adding a K-factor is 'favored' by the Bayesian analysis would be more informative if the manuscript quoted the Bayesian evidence or the effective chi-square improvement from Ref. [13].

Circularity Check

2 steps flagged · score 6.0 of 10

The quoted nuclear suppression factors are not independent CGC predictions: they inherit the K-factor and color-charge-density posterior from Ref. [13], which was fit to the same gamma+Pb data that the suppression is compared with.

  1. fitted input called prediction [Sec. 2 (Model) and Sec. 3 (Results), Eq. (1)]
    "We make use of this minimally extended model and the generated posterior distribution from Ref. [13] to make predictions using 25 parameter samples from the posterior distribution... A value K < 1 is compensated in the fits by a larger color charge density, which corresponds to denser nucleons and, in turn, stronger nuclear suppression."

    The central prediction, namely the A-dependence of coherent J/psi production and the quoted suppression factors (0.3 at W=31.5 GeV and 0.15 at W=813 GeV for heavy nuclei), is computed directly from the posterior distribution obtained in Ref. [13]. That posterior was obtained by fitting the K-factor and the color charge density to the same gamma+Pb data that the present paper's suppression comparison is meant to describe. The paper itself explains that K<1 is compensated by a larger color charge density, which directly increases nuclear suppression. Therefore the A-dependence is not a parameter-free CGC prediction; it is a refit-informed statement whose saturation-scale input is constrained by the very Pb data whose suppression is being reported.

  2. self citation load bearing [Sec. 2 (Model)]
    "In Ref. [13], it was explored whether additional parameters... It turns out that the combined data from a Bayesian fit to the gamma+p and gamma+Pb data disfavors most extensions of the model with additional parameters, but we find that the addition of a K-factor – scaling all cross sections by a constant – is favored."

    The justification for using the K-factor-extended model and for the specific posterior values is taken entirely from Ref. [13], a companion paper by the same authors. The present paper does not independently test or derive this model extension; it imports the posterior as input and then presents the resulting A-dependence as a prediction. Thus the load-bearing model choice and parameters rest on a self-citation that is not independently established within the present paper.

full rationale

The paper is transparent that it uses the posterior distribution from Ref. [13], and the K-factor is admittedly a global rescaling determined from the combined proton and lead data. This transparency does not remove the circularity in calling the resulting A-dependence a prediction: the parameters that set the nuclear saturation scale were constrained by the same gamma+Pb measurements whose suppression is being quantified. A global K-factor alone would cancel in the ratio sigma_A/(A^{4/3} sigma_p), but the paper explicitly states that K<1 is compensated by a larger color charge density, which changes the saturation scale and therefore directly shapes the A-dependence. Consequently the numerical suppression factors are partly determined by the Pb data used for the fit, not purely by CGC dynamics. The self-citation to Ref. [13] is load-bearing because the model extension and the parameter posterior are taken from that same-author work without independent verification here. The paper also notes in Sec. 4 that a more quantitative saturation assessment would require nuclear form factors, which further indicates that the quoted suppression values are not a standalone first-principles result. I therefore score partial circularity at 6 rather than a higher score, because the ratio is still a nontrivial model output and the K-factor does cancel in the ratio for fixed parameters; it is the correlation between K and Q_s, fitted to Pb data, that makes the prediction statistically dependent on the input data.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The model predictions rest on the CGC/JIMWLK framework, on parameters fitted to the same gamma+p and gamma+Pb data in the companion paper, and on a simplified no-saturation baseline. No new entities are introduced. The global K-factor is the only parameter newly emphasized here, but it is fitted, not derived, and it interacts with the color charge density in setting the A-dependence.

free parameters (2)
  • Global K-factor = ~0.33 (from Ref [13])
    Overall normalization multiplying all cross sections in Eq. (1); fitted to gamma+p and gamma+Pb data in the companion global Bayesian analysis. It absorbs vector meson wave function and higher-order uncertainties.
  • Color charge density and related MV/JIMWLK parameters = posterior distribution, 25 samples (values not listed in this paper)
    Parameters of the dipole-target amplitude and initial conditions fitted in Ref [13]; the text says K<1 is compensated by larger color charge density, which directly affects the A-dependence of the suppression.
assumptions (5)
  • domain assumption JIMWLK equation describes the small-x evolution of Wilson lines for both protons and nuclei.
    Used to evolve the dipole-target amplitude from initial conditions; standard but unproven in this paper.
  • domain assumption McLerran-Venugopalan model provides the initial Wilson line distribution.
    Initial condition for the evolution; standard CGC modeling assumption.
  • ad hoc to paper The K-factor is a single global constant, independent of target A and energy W.
    The paper introduces K to improve the simultaneous fit; no A- or W-dependence is modeled or tested, yet the suppression ratio and its energy dependence depend on this assumption.
  • ad hoc to paper A^(4/3) is the correct no-saturation baseline for the coherent cross section ratio.
    Used to define the suppression; the conclusion says nuclear form factors are needed for a quantitative assessment, so this baseline is a simplifying assumption.
  • ad hoc to paper 25 posterior samples from Ref [13] adequately represent the parameter posterior.
    The paper propagates uncertainty using 25 samples; the adequacy of this sample size is not demonstrated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nuclear suppression in diffractive vector meson production within the color glass condensate framework." pith.science (2026). https://pith.science/paper/5PBTDHAM

@misc{pith2026250821562,
  author       = {Pith},
  title        = {Pith review of: Nuclear suppression in diffractive vector meson production within the color glass condensate framework},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5PBTDHAM}},
  note         = {Machine review of arXiv:2508.21562}
}
abstract

We perform a global Bayesian analysis of diffractive $\mathrm{J}/\psi$ production in $\gamma+p$ and $\gamma+\mathrm{Pb}$ collisions within a Color Glass Condensate based framework. Using data from HERA and the LHC, we find that a simultaneous description of $\gamma+p$ and $\gamma+\mathrm{Pb}$ observables is challenging. Introducing a global $K$-factor to account for theoretical uncertainties improves the agreement with data and enhances the framework's predictive power. We present predictions for integrated $\mathrm{J}/\psi$ cross sections at different photon-nucleus energies and study their $A$-dependence relative to a no-saturation baseline, quantifying nuclear suppression and providing insights into the onset of saturation effects.

Figures

Figures reproduced from arXiv: 2508.21562 by the authors.

Figure 1
Figure 1. Ratio of the integrated coherent J/Ψ cross section to the γ + p case as a function of the nuclear mass number A for W = 31.5 and 813 GeV. The dashed line shows the expected A 4/3 scaling in the absence of saturation effects. Error bars indicate the 2σ interval from 25 posterior samples. Using this framework, we presented predictions for the A-dependence of the coherent J/ψ cross section at two photon-nucleus center-… view at source ↗

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 12 linked inside Pith

  1. [13]

    Mäntysaari, H

    H. Mäntysaari, H. Roch, F. Salazar, B. Schenke, C. Shen, W. Zhao, Global Bayesian Analysis of J/ψ Photoproduction on Proton and Lead Targets (2025), 2507.14087

    arXiv 2025

  2. [1]

    Gribov, E.M

    L.V . Gribov, E.M. Levin, M.G. Ryskin, Semihard Processes in QCD, Phys. Rept. 100, 1 (1983). 10.1016/0370-1573(83)90022-4

    doi:10.1016/0370-1573(83)90022-4 1983

  3. [2]

    McLerran, R

    L.D. McLerran, R. Venugopalan, Computing quark and gluon distribution functions for very large nuclei, Phys. Rev. D 49, 2233 (1994), hep-ph/9309289. 10.1103/Phys- RevD.49.2233

    arXiv 1994

  4. [3]

    McLerran, R

    L.D. McLerran, R. Venugopalan, Gluon distribution functions for very large nuclei at small transverse momentum, Phys. Rev. D 49, 3352 (1994), hep-ph/9311205. 10.1103/PhysRevD.49.3352

    arXiv 1994

  5. [4]

    Ryskin, Di ffractive J/ψ electroproduction in LLA QCD, Z

    M.G. Ryskin, Di ffractive J/ψ electroproduction in LLA QCD, Z. Phys. C57, 89 (1993). 10.1007/BF01555742

    doi:10.1007/bf01555742 1993

  6. [5]

    Mäntysaari, Review of proton and nuclear shape fluctuations at high energy, Rept

    H. Mäntysaari, Review of proton and nuclear shape fluctuations at high energy, Rept. Prog. Phys. 83, 082201 (2020), 2001.10705. 10.1088/1361-6633/aba347

    arXiv 2020

  7. [6]

    Eskola, C.A

    K.J. Eskola, C.A. Flett, V . Guzey, T. Löytäinen, H. Paukkunen, Exclusive J /ψ photo- production in ultraperipheral Pb+Pb collisions at the CERN Large Hadron Collider cal- culated at next-to-leading order perturbative QCD, Phys. Rev. C 106, 035202 (2022), 2203.11613. 10.1103/PhysRevC.106.035202

    arXiv 2022

  8. [7]

    Lappi, H

    T. Lappi, H. Mantysaari, J/ψ production in ultraperipheral Pb+Pb and p+Pb collisions at energies available at the CERN Large Hadron Collider, Phys. Rev. C 87, 032201 (2013), 1301.4095. 10.1103/PhysRevC.87.032201

    arXiv 2013

Show all 17 references
  1. [8]

    Mäntysaari, F

    H. Mäntysaari, F. Salazar, B. Schenke, Nuclear geometry at high energy from ex- clusive vector meson production, Phys. Rev. D 106, 074019 (2022), 2207.03712. 10.1103/PhysRevD.106.074019

    2022 arXiv

  2. [9]

    Mäntysaari, B

    H. Mäntysaari, B. Schenke, Evidence of strong proton shape fluctuations from inco- herent diffraction, Phys. Rev. Lett. 117, 052301 (2016), 1603.04349. 10.1103/Phys- RevLett.117.052301

    2016 arXiv

  3. [10]

    Mäntysaari, B

    H. Mäntysaari, B. Schenke, Confronting impact parameter dependent JIMWLK evolu- tion with HERA data, Phys. Rev. D 98, 034013 (2018), 1806.06783. 10.1103/Phys- RevD.98.034013

    2018 arXiv

  4. [11]

    Acharya et al

    S. Acharya et al. (ALICE), Energy dependence of coherent photonuclear production of J/ψ mesons in ultra-peripheral Pb-Pb collisions at √sNN = 5.02 TeV, JHEP 10, 119 (2023), 2305.19060. 10.1007/JHEP10(2023)119

    2023 arXiv

  5. [12]

    Tumasyan et al

    A. Tumasyan et al. (CMS), Probing Small Bjorken-x Nuclear Gluonic Structure via Coherent J/ψ Photoproduction in Ultraperipheral Pb-Pb Collisions at sNN =5.02 TeV, Phys. Rev. Lett. 131, 262301 (2023), 2303.16984. 10.1103/PhysRevLett.131.262301

    2023 arXiv

  6. [14]

    Mueller, A Simple derivation of the JIMWLK equation, Phys

    A.H. Mueller, A Simple derivation of the JIMWLK equation, Phys. Lett. B 523, 243 (2001), hep-ph/0110169. 10.1016/S0370-2693(01)01343-0

    2001 arXiv

  7. [15]

    Mäntysaari, R

    H. Mäntysaari, R. Venugopalan, Systematics of strong nuclear amplification of gluon saturation from exclusive vector meson production in high energy electron–nucleus col- lisions, Phys. Lett. B 781, 664 (2018), 1712.02508. 10.1016/j.physletb.2018.04.044

    2018 arXiv

  8. [16]

    Pordes et al., The Open Science Grid, J

    R. Pordes et al., The Open Science Grid, J. Phys. Conf. Ser. 78, 012057 (2007). 10.1088/1742-6596/78/1/012057

    2007 doi

  9. [17]

    Sfiligoi, D.C

    I. Sfiligoi, D.C. Bradley, B. Holzman, P. Mhashilkar, S. Padhi, F. Wurthwrin, The pi- lot way to Grid resources using glideinWMS, WRI World Congress 2, 428 (2009). 10.1109/CSIE.2009.950

    2009 doi

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.