{"id":"813a391b-aabc-4551-9b03-293e575c55f3","arxiv_id":"2607.14955","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A global fit to all 104 HERA exclusive J/ψ data points resolves the proton's small-x gluon hotspots into a 0.105 fm Gaussian core plus a 0.22 fm exponential halo.","lead":"This paper analyzes HERA measurements of exclusive J/ψ production and reports that the proton's small-x gluon cloud has a compact Gaussian core of about 0.105 fm surrounded by a wider exponential halo of about 0.22 fm. The result gives a concrete geometric picture of the proton's gluon distribution that can feed into heavy-ion and future electron-ion collision calculations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed 0.105 fm Gaussian core + 0.22 fm exponential halo is the assumed Gaussian x K0 hotspot profile (Eqs. 21-23), not an unbiased resolution; no alternative shapes are compared.","rationale":"The reader's weakest assumption identifies the key issue: the core-halo structure is imposed by the chosen Gaussian x K0 convolution, and no alternative shapes are tested. This is the load-bearing point because the headline claim is about a resolved two-component geometry. The paper is otherwise internally consistent: the analytic derivation (Eqs. 9-20) is sound, the CM shift is handled correctly, and the fit quality is good. The agreement of lambda with lattice predictions is a point in favor of the model's physical relevance. However, because the form factor is assumed, the extracted core/halo sizes are conditional on that assumption. A test against alternative functional forms would settle whether the data actually resolve the two components. The reader's CONDITIONAL verdict is therefore appropriate; no change is needed.","tokens_in":10400,"tokens_out":12483,"duration_ms":127014,"concrete_test":"Re-run the global fit (104 points, same normalization and fluctuation model) with two alternative hotspot form factors: (i) pure Gaussian \\hat T=exp(-B_eff|t|/2), and (ii) a single power-law \\hat T=(1+B|t|)^{-p} (both with one width parameter, analogous to B_hs). Compare Delta chi^2 and AIC/BIC relative to the Gaussian x K0 form. If the best alternative achieves Delta chi^2 < 9 for the same parameter count (or Delta AIC < 6), the data do not require the exponential halo, and the 'core+halo' claim should be downgraded to a parameterization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim restates the model's assumed functional form. In Eqs. (21)-(22) the hotspot thickness is set to T_hs = Gaussian * K0, giving the factorized form factor \\hat T(t)=exp(-B_hs|t|/2)/(1+lambda^2|t|) (Eq. 23). The fit returns the parameters of this imposed profile. A good chi2 (0.77) demonstrates consistency, but not uniqueness. The small-|t| slope only constrains the combination B_hs/2 + lambda^2, so the separation into a '0.105 fm core' and '0.22 fm halo' is determined by the assumed analytic shape and by the large-|t| tail, which is a particular 1/(lambda^2|t|) Lorentzian. Alternatives such as a pure power-law form factor (1+B|t|)^{-p} or a sum of two Gaussians could plausibly describe the same data with the same number of parameters. No model-selection comparison is shown. Claims of x_IP- and N_hs-independence are likewise tested only within this one-parametric family. The agreement of lambda with the lattice correlation length is suggestive, but not evidence that the data resolve two distinct components.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytic treatment of coherent and incoherent exclusive J/psi photoproduction in the dipole model with a hotspot proton structure. Using the Good-Walker formalism, the authors derive closed-form expressions for the first and second moments of the hotspot thickness function, including a center-of-mass correction and log-normal saturation-scale fluctuations. They fit the resulting model to 104 H1 and ZEUS data points and report chi^2/ndf = 0.77. The main physics claim is that the gluonic hotspots have a resolved transverse structure of a perturbative Gaussian core (radius ~0.105 fm) surrounded by an exponential halo of range ~0.220 fm, independent of x_IP and of the number of hotspots, and consistent with lattice QCD studies of vacuum flux tubes. They also extract a small effective pomeron slope alpha'_eff ~0.046 GeV^-2.","tokens_in":10861,"tokens_out":7222,"duration_ms":74999,"significance":"The analytic derivation is a useful technical advance: it removes the need for Monte Carlo averaging over hotspot configurations and allows a global fit to the full coherent and incoherent t-spectra. The fit quality is good, and the simultaneous description of coherent and incoherent HERA data is a nontrivial success. The connection to lattice-gauge-theory correlation lengths is suggestive. However, the central 'resolved core-halo' claim is not model-independent: the Gaussian-circle?K0 profile is imposed before fitting, and no alternative functional forms are tested. If reframed as a parameter extraction under a specific assumed hotspot profile, the paper is a solid phenomenological contribution; as stated, the abstract overstates the uniqueness of the extracted geometry. The shot-noise interpretation in the abstract is also not derived from data.","major_comments":[{"comment":"The central claim that the data resolve a 0.105 fm Gaussian core and a 0.220 fm exponential halo is not supported as a model-independent measurement. The hotspot profile is assumed to be a convolution of a Gaussian and K0 before fitting (Eqs. 21-22), which fixes the factorized form Eq. (23). At small |t| the coherent cross section constrains essentially the combination B_hs/2 + lambda^2 (together with B_qc), while the separation into a core and a halo is controlled by the assumed analytic form and by the particular 1/(lambda^2|t|) tail. No alternative shapes (e.g., two Gaussians, a power-law form factor, or a single broader Gaussian) are compared. A chi^2/ndf=0.77 shows consistency, not uniqueness. Please either add a model-selection comparison or explicitly present the result as the parameterization of the assumed Gaussian x K0 profile.","section":"Abstract and 'The Model' (Eqs. 21-23)"},{"comment":"The statement that the small-|t| incoherent cross section is dominated by shot-noise fluctuations of order one gluon per hotspot is not a result of the fit. sigma_S is a fitted parameter controlling the variance of a log-normal fluctuation xi_i. The observation that Var[xi_i]/<xi_i>^2 ~ 1 is a reparameterization of this fitted variance, not an independent measurement. The actual number of gluons per hotspot never enters the calculation. This interpretive claim should be removed or clearly labeled as an analogy, not a determination.","section":"'Comparisons to HERA data', sigma_S paragraph and Abstract"},{"comment":"The conclusion that the hotspot shape is 'frozen in x_IP' overstates the sensitivity. Fits 2 and 3 add one parameter each and improve chi^2 by only about 2 and 0.5, respectively. This means the HERA data are consistent with no x_IP dependence, but they do not positively establish that no dependence exists. Please phrase the result as a null test with an upper bound (e.g., beta' and delta_N limits) and state the sensitivity of the extracted B_hs and lambda to plausible x_IP variations.","section":"'Comparisons to HERA data' (Fit 2 and Fit 3) and Conclusions"}],"minor_comments":[{"comment":"The text says six free parameters are fitted, but Table I lists only five parameter rows (the overall normalization is mentioned in the text but not tabulated). Please include the normalization factor in the table or explicitly state why it is omitted.","section":"Table I and text"},{"comment":"The abstract and conclusions mention independence of 'hotspot repulsion', but the only treatment is the short ad-hoc statement in Conclusions about a 30% reduction of coincident hotspots over 0.3 fm with Delta chi^2=1. No details or fit plots are given. Either provide the analysis or remove this claim from the abstract.","section":"Abstract and Conclusions"},{"comment":"The lattice comparison is somewhat confusing: the flux-tube penetration depth quoted in the Introduction is ~0.1 fm, while the extracted lambda=0.220(2) fm is compared to the vacuum correlation length. The abstract says the halo agrees with the 'core-halo structure of flux tubes', but the numerical agreement is with a different quantity. Please clarify which lattice prediction is being compared and what the quoted uncertainties are.","section":"Introduction, lattice comparison"},{"comment":"There are several typos and OCR artifacts (e.g., 'Mandelstamt', 'IPsat' vs 'IPnonsat' spacing, 'Bessel function' consistency). A careful proofreading pass is needed.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid analytic core and a good global fit, but the headline claim is presented as a measurement when it is actually a fit under an assumed profile. The requested alternative-shape comparison and reframing are feasible within the manuscript's scope; I would support publication after these are addressed. The shot-noise interpretation should be removed from the abstract unless a concrete model links sigma_S to gluon occupancy."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the paper in two sentences: it derives analytic expressions for coherent and incoherent exclusive J/psi cross sections in the hotspot model, including a center-of-mass correction, and fits them to all 104 HERA points with chi2/ndf=0.77. The fit yields a hotspot form factor exp(-B_hs|t|/2)/(1+lambda^2|t|), corresponding to a 0.105 fm Gaussian core plus a 0.220 fm exponential halo, stable across N_hs and x_IP.\n\nThe analytic derivation is clean and, as far as I can tell, correct. The center-of-mass correction in eqs. (16)-(19) properly kills the N_hs=1 variance. The fit is genuinely good over three decades in t, which is a real achievement for this class of models. The stability of the extracted parameters against N_hs and against allowing x_IP dependence (fits 2 and 3) is credible evidence that the data constrains these parameters within the model.\n\nThe soft spot is exactly what the stress-test says: the core-halo structure is put in by hand via eqs. (21)-(23), not measured from the data. At small |t| the slope constrains B_hs/2 + lambda^2, and the separation into core and halo rests on the assumed 1/(1+lambda^2|t|) tail. I would have liked a comparison against, say, a power-law form factor (1+B|t|)^-p or a sum of two Gaussians with the same number of parameters. No such comparison appears. As a result, the abstract's \"resolved into\" overstates what the data alone can say; \"consistent with\" would be more honest. The same caveat applies to the \"shot-noise one gluon per hotspot\" claim, which is an interpretation of sigma_S, not a direct counting. None of this kills the paper, but it should be framed as a model-dependent extraction.\n\nAlso minor: the global normalization comes out at 1.3, a bit large but they argue it absorbs known wave-function uncertainties and doesn't affect the t-shapes. No code or data tables are shipped, which makes independent checking harder.\n\nBottom line: this is a serious paper from a technical standpoint. It deserves peer review. The referee should ask for alternative shape comparisons and ideally the fit implementation. I'd also recommend toning down the abstract.","headline":"Solid fit and clean analytics, but the core+halo result is built into the assumed form factor, not measured from HERA data.","tokens_in":11280,"tokens_out":2951,"would_cite":true,"duration_ms":29934,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new analysis of exclusive J/ψ photoproduction at HERA claims that the proton's small-x gluon distribution has a two-component transverse shape: a perturbative Gaussian core of about 0.105 fm and a nonperturbative exponential halo of about","keywords":["gluon cloud","proton structure","exclusive J/psi photoproduction","hotspot model","coherent diffraction","incoherent diffraction","small-x gluons","core-halo structure"],"falsifier":"A fit of the same 104-point data set with a single-Gaussian hotspot profile (form factor e^{-B|t|}) that yields a comparable χ² would falsify the claim that the data resolve a distinct exponential halo; alternatively, a measurement of the incoherent cross section at |t| > 5 GeV² that follows an exponential rather than the predicted power-law tail (1+λ²|t|)^{-2} would invalidate the core-halo extraction.","tokens_in":10317,"feed_emoji":"⚛️","tokens_out":6329,"duration_ms":60000,"temperature":0.7,"pith_summary":"The paper claims that the transverse shape of the proton's small-x gluon distribution can be extracted from exclusive J/ψ photoproduction data at HERA, and that the shape consists of a perturbative Gaussian core of radius about 0.105 fm surrounded by an exponential halo of range about 0.220 fm. The authors derive analytic expressions for coherent and incoherent cross sections in the hotspot model, enabling a global fit to all 104 available data points without Monte Carlo event generation. If true, this would give a precise determination of the gluon cloud's geometry and connect it to QCD vacuum properties and lattice flux-tube structure. The shape is claimed to be independent of the momentum fraction x_IP and of the assumed number of hotspots.","feed_headline":"Proton's gluon cloud has a 0.105 fm core and a 0.220 fm halo","feed_subtitle":"A fit to 104 HERA data points reveals a perturbative Gaussian core and an exponential halo matching lattice QCD predictions.","key_machinery":"The central object is the proton's thickness function T_p(b) in the hotspot model, built from independent hotspots with a profile given by the convolution of a Gaussian of width √(2B_hs) and a modified Bessel function K0(r/λ). This convolution yields the analytic form factor ˆT(t) = e^{-B_hs|t|/2}/(1+λ²|t|), which combines a Gaussian core with a Lorentzian tail and produces the core-halo structure. Analytic Good–Walker expressions for coherent and incoherent cross sections, including a center-of-mass correction and log-normal saturation-scale fluctuations, allow a direct fit to the data without averaging over thousands of hotspot configurations.","core_discovery":"The central discovery is that the gluon distribution inside the proton is not a simple Gaussian but a composite: a perturbative Gaussian core of size √(2B_hs) = 0.105(2) fm plus a nonperturbative exponential halo of range λ = 0.220(15) fm, obtained by fitting analytic leading-twist expressions to all 104 H1 and ZEUS exclusive J/ψ photoproduction data points across three decades of momentum transfer |t|. The fit quality is χ²/ndf = 0.77, and the geometry remains stable when the number of hotspots or the x_IP dependence is varied. The incoherent cross section at small |t| is dominated by shot-noise fluctuations at the level of about one gluon per hotspot.","pith_inferences":["If the assumed hotspot profile is correct, the core and halo sizes are parameters of that functional form, not model-independent measurements; a direct test would be to refit the same data with alternative smooth profiles (e.g., a single Gaussian or a different tail index) and compare the fit quality.","The result suggests that the Gaussian core and exponential halo might reflect two distinct physical regimes: perturbative bremsstrahlung near the hotspot center and nonperturbative vacuum screening at larger distances. One could test whether the halo scale λ tracks the string tension or the lightest glueball mass as the photon virtuality Q² varies.","If the geometry is truly frozen in x_IP, the energy dependence of exclusive production enters primarily through the DGLAP-evolved gluon density and the diffusion of hotspot centres, not through the shape of individual hotspots. This would simplify initial-state modelling for heavy-ion collisions: the fluctuating gluon profile can be generated once and evolved by simple diffusion.","The shot-noise interpretation (about one gluon per hotspot) is an inference from the fitted log-normal fluctuation width; a direct measurement of the variance of the incoherent cross section at even smaller |t|, where the geometric variance vanishes, could confirm or rule out this picture."],"forward_implications":["If the core-halo geometry is real, the proton's small-x gluon distribution has a well-defined transverse shape that is largely frozen in x_IP, with only a slow diffusion of the hotspot centres.","The extracted halo range of about 0.220 fm matches the gluon-field correlation length of the QCD vacuum and the flux-tube core-halo structure seen on the lattice, suggesting a universal nonperturbative scale.","The fitted effective pomeron slope α'_eff ≈ 0.046 GeV⁻² is smaller than previous HERA Regge extractions but consistent with running-coupling BFKL expectations, affecting how the gluon cloud expands with energy.","The interpretation of the small-|t| incoherent cross section as shot-noise from roughly one gluon per hotspot constrains event-by-event fluctuations of the saturation scale.","The analytic expressions provide a fast, Monte-Carlo-free framework for interpreting future Electron-Ion Collider measurements of exclusive vector meson production."],"fun_headline_variants":["Proton's gluon: Gaussian core + exponential halo","Gluon cloud inside proton: 0.105 fm core, 0.220 fm halo","HERA data reveals proton's gluon core-halo structure","Gluon geometry: core 0.105 fm, halo 0.220 fm","Proton gluons: perturbative core, nonperturbative halo"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the hotspot thickness function has the assumed Gaussian-convolved-K0 profile (Eqs. 21–22); the extracted 'core' and 'halo' are parameters of that assumed shape, not independent measurements, and the paper does not compare against alternative functional forms.","fun_headline_variants_meta":{"raw":{"variants":["Proton's gluon: Gaussian core + exponential halo","Gluon cloud inside proton: 0.105 fm core, 0.220 fm halo","HERA data reveals proton's gluon core-halo structure","Gluon geometry: core 0.105 fm, halo 0.220 fm","Proton gluons: perturbative core, nonperturbative halo"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000277,"raw_usage":{"total_tokens":1510,"prompt_tokens":791,"completion_tokens":719,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":618}},"tokens_in":535,"tokens_out":719,"duration_ms":5885,"temperature":1.0,"reasoning_tokens":618,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T00:36:05.523319+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A fit of the same 104-point data set with a single-Gaussian hotspot profile (form factor e^{-B|t|}) that yields a comparable χ² would falsify the claim that the data resolve a distinct exponential halo; alternatively, a measurement of the incoherent cross section at |t| > 5 GeV² that follows an exponential rather than the predicted power-law tail (1+λ²|t|)^{-2} would invalidate the core-halo extraction.","supporting_citations":[],"review_version":1}