{"id":"ab67e146-5b6d-421e-995d-89be8ffdcf9e","arxiv_id":"1908.03494","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Exclusive J/psi photoproduction at the LHC favors a gluon distribution with nonlinear saturation corrections over linear BFKL-based predictions, suggesting onset of gluon saturation.","lead":"This paper compares two theoretical descriptions of exclusive J/psi and Upsilon photoproduction at HERA and the LHC, one based on linear QCD evolution and one including nonlinear gluon saturation effects. The authors find that only the nonlinear description matches the energy rise of J/psi data at LHC energies, and they interpret this as evidence for the onset of gluon saturation in the proton.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The saturation claim depends on the LHC-region overshoot of the linear HSS baseline exceeding its renormalization-scale ambiguity; the paper quantifies that ambiguity only for HERA, not for the LHC region where the overshoot is claimed.","rationale":"The paper's logic is a controlled comparison: linear BFKL should fail where nonlinear BK succeeds. For that comparison to constitute evidence, the linear failure must be robust against reasonable theory uncertainties. The paper provides no uncertainty bands on any prediction, and the only quantitative statement about scale sensitivity is in the HERA region, where the claim is not being made. The LHC region is precisely where the overshoot appears after a scale-setting change introduced to stabilize the perturbative expansion; whether that change is representative of the theory uncertainty is not demonstrated. The reader's weakest-assumption pick identifies exactly this gap. I find no stronger objection: the claim is not internally inconsistent, the methods are standard, and the main result is backed by a longer published paper [6]. The weakness is quantitative underdetermination, which supports a conditional verdict rather than rejection. One additional consideration is that the 'linear KS' comparison is a useful internal control, but without a precise definition of the deactivated nonlinear terms, it cannot independently validate the HSS-based conclusion. The proposed test—systematic scale variation with χ²—would settle whether the LHC-region overshoot is real or within theory noise.","tokens_in":6002,"tokens_out":2653,"duration_ms":27735,"concrete_test":"Recompute the HSS and linear-KS J/ψ photoproduction predictions in the LHC region under a systematic renormalization-scale variation, for example μ_R = M/2, M, 2M for the fixed-scale choice and μ_R = μ_r/2, μ_r, 2μ_r for the r-dependent choice, and report χ² values or band coverage for the LHC data points. If the linear predictions' scale envelope overlaps either the data or the nonlinear KS band, the overshoot is not a robust saturation signal; if the envelope lies entirely above the data and above the KS band, the claim survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 establishes the central falsification test: a linear NLO BFKL gluon (HSS) should fail where a nonlinear BK gluon (KS) succeeds. The evidence is the overshoot of the dashed green HSS line relative to LHC J/ψ data after switching from fixed M² = 3.27 GeV² to an r-dependent renormalization scale. The paper explicitly notes that the fixed- versus r-dependent-scale difference is 'of the order of a typical variation of the renormalization scale' for the HERA region, but it never gives the corresponding scale-variation envelope in the LHC region. If that envelope is comparable to the plotted overshoot, the 'breakdown of linear evolution' is not a genuine nonlinear effect but an artifact of one scale choice. The same gap applies to the KS-gluon with non-linearities turned off (the dashed black line), which is presented as the controlled linear limit without a precise definition of which terms are removed and without uncertainty estimates. Without these quantitative baselines, the conclusion in Sec. 3 that non-linear effects are essential, and the 'clear sign' wording, outrun the evidence shown. The concern is not that the fits are wrong; it is that the discriminating quantity—the size of the LHC-region overshoot relative to theory uncertainty—is not reported. This is exactly the condition that must hold for the central claim to be true.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares two unintegrated gluon densities in exclusive J/psi and Upsilon photoproduction: the nonlinear BK-based Kutak-Sapeta (KS) gluon and the linear NLO BFKL-based Hentschinski-Sabio Vera-Salas (HSS) gluon. After adopting an r-dependent renormalization scale for the HSS gluon, the authors report that the linear HSS prediction overshoots LHC-region J/psi data, whereas the nonlinear KS gluon describes the energy dependence; switching off nonlinearities in the KS gluon also leads to an overshoot. The paper interprets this difference as a sign of the onset of gluon saturation.","tokens_in":6298,"tokens_out":4041,"duration_ms":43126,"significance":"If established, this would be a valuable phenomenological indication of nonlinear small-x QCD dynamics in an otherwise standard exclusive process, and the choice of a linear BFKL baseline rather than a DGLAP fit is well motivated. The comparison is not directly circular because the KS and HSS gluon parameters are determined from inclusive HERA data, not from the LHC photoproduction points used for the test. The main strengths are the clear falsification logic and the explicit use of two different evolution frameworks. However, the evidence presented is qualitative: no theory uncertainty bands, no goodness-of-fit measures, and no explicit definition of the 'linear KS' limit. The central claim therefore currently rests on the visual size of the LHC-region overshoot relative to an unquantified renormalization-scale variation.","major_comments":[{"comment":"The central discriminator is whether the LHC-region overshoot of the HSS prediction exceeds the theoretical uncertainty of the linear benchmark. The paper states that the difference between the fixed-scale and r-dependent-scale HSS solutions is 'of the order of a typical variation of the renormalization scale' only for the HERA region; no analogous scale-variation envelope is given in the LHC region, where the overshoot is claimed. Without this envelope, the plotted dashed-green overshoot cannot be distinguished from a scale-setting artifact, and the conclusion that linear evolution fails at LHC energies is not quantitatively supported.","section":"Sec. 3, discussion of Fig. 2"},{"comment":"The 'KS-gluon with non-linearities turned off' is presented as the controlled linear limit that demonstrates the need for nonlinear terms, but the manuscript does not define which terms are removed (for example, whether the BK kernel is set to its linear part or the evolution is stopped after the initial condition) and does not provide uncertainty estimates for this limit. Since this curve is part of the key falsification test, its precise definition and reasonable uncertainty band are required before the comparison can be interpreted.","section":"Sec. 3, dashed black curve in Fig. 1"},{"comment":"No uncertainty bands on any of the theory curves and no chi-squared or other goodness-of-fit values are reported. In the LHC J/psi region the data have appreciable scatter, and the visual difference between the HSS overshoot and the data is comparable in size to that scatter. Reporting the deviation of each theory prediction from the data in units of the experimental uncertainty would turn the claimed failure of HSS and success of KS into a quantitative statement.","section":"Sec. 3 and Fig. 1"},{"comment":"The abstract claims that linear NLO BFKL evolution can describe the highest-energy data only if perturbative corrections are increased to 'unnaturally large values', but the manuscript never quantifies 'unnaturally large'. The text gives the fixed scale M^2 = 3.27 GeV^2 and the r-dependent scale, and Fig. 2 shows the correction becoming large, but no explicit numerical comparison of the correction to the leading term is given for the LHC kinematics. A concrete ratio or scale value would substantiate this claim.","section":"Abstract and Sec. 3"}],"minor_comments":[{"comment":"The phrase 'in low x the proton' should read 'in the low-x proton' or 'in the proton at low x'.","section":"Abstract"},{"comment":"The sentence 'the corrections supersedes the formally leading term' contains a subject-verb disagreement; 'supersede' is the correct form with the plural subject 'corrections'.","section":"Sec. 3, text near Eq. (3.2)"},{"comment":"The figure has no legend and the caption does not identify which line corresponds to which gluon or scale setting; the text description of colors is not usable in monochrome print.","section":"Fig. 1"},{"comment":"The notation for the arguments of the dipole cross-section and the ratio \\bar{\\alpha}_s(M\\cdot Q_0)/\\bar{\\alpha}_s(M^2) should be defined explicitly, since the scale choices are central to the argument.","section":"Eqs. (3.1)-(3.2)"},{"comment":"The statement that collinear PDFs provide no perturbative evolution for J/psi photoproduction would benefit from a quantitative illustration or a reference to a collinear-fit analysis of the same observable, so that the reader can assess the claimed disadvantage relative to the BFKL benchmark.","section":"Sec. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a short proceedings contribution, and the conclusions are stated more strongly than the quantitative evidence shown. The missing scale-variation envelope and the undefined 'linear KS' limit are the key issues; if the authors can provide those, the manuscript could become a useful, if modest, addition. The reference list is heavily self-referential (refs. [3,6,7,14,15]), which is understandable for a methods-focused letter but should be balanced with independent checks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things up front. First, the central result is not new: essentially the same comparison and conclusion appear in the authors' PLB 795, 569 (2019) paper. This LHCP proceedings adds an r-dependent renormalization scale for the HSS gluon and a 'linear KS' variant, but the load-bearing claim—linear BFKL-like evolution fails, BK-based KS succeeds—is unchanged. Second, the new piece is exactly where the argument gets soft. The stress-test note is right: the authors say the fixed and r-dependent scale versions of HSS differ by 'a typical variation of the renormalization scale' in the HERA region, but they never quantify that spread in the LHC region, which is where the claimed overshoot happens. If the two HSS curves diverge by as much as the plotted gap, then 'breakdown of linear evolution' could just be a scale artifact.\n\nOn the positive side, the methodological thinking is sound. They correctly reject DGLAP-based collinear PDFs as a benchmark for this observable, because the x-dependence is effectively a fit at the charm scale. Using a BFKL-based unintegrated gluon, fitted at moderate x and evolved to low x, is a more honest linear baseline. The observation that the NLO BFKL subleading term grows and overtakes the leading term at small x is concrete, and the r-dependent scale is a reasonable attempt to stabilize it. The comparison between the full KS gluon and its linear limit is the right kind of controlled test. The circularity burden is moderate, not high: the fits are to inclusive HERA data, so the LHC exclusive points are a genuine extrapolation. But both models come from the same group, there are no chi-squared values, and no uncertainty bands on any curve. The 'linear KS' definition is vague—'non-linearities turned off'—which matters because the KS initial conditions were fitted with the full nonlinear evolution.\n\nOverall, this is a useful proceedings pointer, but the evidence is qualitative. 'Clear sign' is stronger than what the figures justify. If this goes to a journal, a referee should ask for the LHC-region scale-variation envelope, uncertainty bands, and a precise definition of the linear limit. I would not cite this proceedings instead of the PLB paper unless I specifically wanted the r-dependent scale variant. It could be worth a coffee discussion or a reading group, but it doesn't settle the saturation question.","headline":"A suggestive but non-quantitative proceedings reprint: the central saturation claim is already in the authors' PLB paper, and the one new piece—the r-dependent scale comparison—leaves the LHC-region scale ambiguity unquantified.","tokens_in":6850,"tokens_out":4534,"would_cite":false,"duration_ms":42039,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Exclusive J/psi photoproduction at the LHC shows the onset of gluon saturation.","keywords":["gluon saturation","exclusive vector meson photoproduction","J/psi photoproduction","Balitsky-Kovchegov evolution","BFKL evolution","unintegrated gluon distribution","low-x QCD","ultra-peripheral collisions"],"falsifier":"Compute the stabilized linear prediction for $J/\\Psi$ photoproduction across the full LHC energy range with that internal scale varied over the same range that reproduces the HERA-region variation; if the linear prediction can be made to match the LHC data within that band, the conclusion that nonlinear effects are required is not established.","tokens_in":5757,"feed_emoji":"⚛️","tokens_out":7341,"duration_ms":73840,"temperature":0.7,"pith_summary":"The paper tries to establish that the energy rise of exclusive $J/\\Psi$ photoproduction measured at HERA and the LHC contains a clear signal of gluon saturation in the low-$x$ proton. The strategy is comparative: a linear perturbative baseline built on next-to-leading-order BFKL evolution (the HSS gluon) should describe the data if no saturation were present, while a fit based on nonlinear Balitsky-Kovchegov evolution (the KS gluon) should fail without its nonlinear terms. The paper shows that the stabilized linear baseline overshoots the LHC $J/\\Psi$ data, while the nonlinear KS gluon describes them; the $\\Upsilon$ data, at a larger scale, remain consistent with both. This matters because a failure of the best available linear low-$x$ evolution, together with success of the nonlinear alternative, is the kind of evidence that distinguishes saturation from a generic fit.","feed_headline":"Exclusive J/psi data at the LHC require nonlinear gluon evolution","feed_subtitle":"Linear BFKL growth overshoots the LHC J/psi cross-section unless corrections are inflated; the nonlinear BK fit matches the data.","key_machinery":"The central object is the unintegrated gluon distribution and its dipole cross-section transform, used to compute exclusive vector-meson photoproduction. For the HSS gluon the dipole cross-section factorizes as a formally leading BFKL term plus a next-to-leading-logarithmic correction $\\propto \\bar\\alpha_s^2 \\beta_0 \\chi_0(\\gamma)\\log(1/x)$; that correction is the mechanism that breaks the fixed-scale linear expansion at small $x$ and motivates an $r$-dependent renormalization scale. For the KS gluon the machinery is the nonlinear BK evolution equation, whose damping of the gluon growth at high density tames the cross-section rise. The comparison of these two computations against the $J/\\Psi$ and $\\Upsilon$ photoproduction data is what carries the argument.","core_discovery":"On the paper's own terms, the central discovery is that non-linear corrections to low-$x$ QCD evolution are essential to describe the energy dependence of exclusive $J/\\Psi$ photoproduction in the LHC region, and that this necessity is a sign of the onset of gluon saturation. The authors construct this by taking two fits of the unintegrated gluon distribution: the HSS gluon, from linear NLO BFKL evolution, and the KS gluon, a solution to nonlinear BK evolution fitted to combined HERA data. With the renormalization scale fixed at a heavy-quark scale, the formally sub-leading perturbative correction in the HSS dipole cross-section grows toward small $x$ and eventually dominates, so an $r$-dependent scale is used to stabilize the linear prediction. The stabilized linear HSS prediction agrees with the HERA-region data within typical scale variation but overshoots the LHC $J/\\Psi$ data, and the linearized KS gluon overshoots as well; the full nonlinear KS gluon describes the data. The paper reads this pattern as evidence that a high-density, saturating gluon configuration is being probed.","pith_inferences":["An extension the authors do not make: apply the same two gluons to other exclusive processes at comparable and smaller $x$; the same pattern should appear if the interpretation is right and disappear if the $J/\\Psi$ result is peculiar to that channel.","The comparison's power rests on quantifying the scale-variation band of the linear prediction in the LHC region; an explicit envelope would turn the claimed overshoot from an interpretation into a metrologically testable statement.","A quantitative measure of the onset could be formed from the ratio of nonlinear to linear cross-sections at each LHC energy and compared with saturation-scale estimates from other low-$x$ processes.","If the interpretation is right, the fixed-scale perturbative expansion in $\\log(1/x)$ should break down exactly where the $J/\\Psi$ data deviate, so repeating the analysis at higher resummation order predicts where linear evolution regains validity."],"forward_implications":["A successful nonlinear BK-based gluon fit to inclusive HERA data transfers to exclusive $J/\\Psi$ photoproduction in the LHC region, while the linear BFKL benchmark fails; the low-$x$ gluon growth implied by the data is damped.","The $\\Upsilon$ channel stays in the perturbative regime and should remain describable by linear evolution, providing a scale-dependent cross-check of the saturation interpretation.","Further exclusive $J/\\Psi$ photoproduction measurements at smaller $x$ should show a flatter-than-linear rise if saturation is setting in.","Linear frameworks that want to describe the LHC $J/\\Psi$ data must invoke unnaturally large perturbative corrections, which the paper treats as a symptom rather than a viable alternative.","The conclusion depends on the stabilized linear calculation, so improving the theoretical accuracy of the linear low-$x$ framework directly sharpens the saturation signal."],"supporting_citations":[{"why":"Supplies the HSS unintegrated gluon, the linear NLO BFKL benchmark whose failure drives the argument.","marker":"[3]"},{"why":"Supplies the KS gluon, the nonlinear BK-evolution fit whose success with nonlinear terms is the positive side of the comparison.","marker":"[8]"},{"why":"Contains the detailed derivation of the HSS dipole cross-section decomposition that exposes the breakdown of the fixed-scale linear expansion.","marker":"[6]"},{"why":"Supports the claim that fixed- and r-dependent-scale HSS solutions agree within typical renormalization-scale variation in the HERA region.","marker":"[7]"},{"why":"Defines the BK evolution equation underlying the nonlinear KS gluon.","marker":"[5]"},{"why":"Provides the HERA J/psi photoproduction data used to anchor the low-energy behavior.","marker":"[9]"},{"why":"Provides the LHC J/psi photoproduction data that the linear gluons overshoot and the nonlinear gluon describes.","marker":"[10]"},{"why":"Provides the LHC Upsilon photoproduction data used as a perturbative cross-check at a larger scale.","marker":"[12]"}],"fun_headline_variants":["Nonlinear gluon evolution required to match LHC J/psi data","LHC J/psi data expose need for gluon saturation","J/psi photoproduction reveals gluon saturation onset","Linear BFKL fails LHC J/psi; nonlinear BK succeeds","Exclusive J/psi data signal onset of gluon saturation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the linear benchmark calculation is still trustworthy after its internal energy scale is adjusted; if its overshoot at LHC energies is just a leftover of that adjustment, the claimed saturation signal disappears.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinear gluon evolution required to match LHC J/psi data","LHC J/psi data expose need for gluon saturation","J/psi photoproduction reveals gluon saturation onset","Linear BFKL fails LHC J/psi; nonlinear BK succeeds","Exclusive J/psi data signal onset of gluon saturation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000245,"raw_usage":{"total_tokens":1585,"prompt_tokens":1047,"completion_tokens":538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":450}},"tokens_in":663,"tokens_out":538,"duration_ms":5775,"temperature":1.0,"reasoning_tokens":450,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:11:42.986188+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the stabilized linear prediction for $J/\\Psi$ photoproduction across the full LHC energy range with that internal scale varied over the same range that reproduces the HERA-region variation; if the linear prediction can be made to match the LHC data within that band, the conclusion that nonlinear effects are required is not established.","supporting_citations":[{"cited_title":"The hard to soft Pomeron transition in small x DIS data using optimal renormalization","cited_arxiv_id":"1209.1353","evidence_quote":"Supplies the HSS unintegrated gluon, the linear NLO BFKL benchmark whose failure drives the argument."}],"review_version":1}