{"id":"662f504d-8a85-4f10-9d74-2a016aebd86f","arxiv_id":"2502.01712","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"After summing large Pomeron loops, the dipole-nucleus amplitude has the same energy dependence as dipole-dipole scattering, limiting the BK equation to z' below roughly 2 sqrt(kappa c) A^{1/6}.","lead":"A theory paper calculates dipole-nucleus scattering at very high energy by summing large Pomeron loop corrections that the standard BK equation misses. It finds the same energy dependence as dipole-dipole scattering and concludes the BK equation is reliable only below an energy bound that grows as the sixth root of the atomic number.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central 1/(4κ) exponent inherits an unproven ansatz for the projectile dipole densities in Eq. (17); if those factorial moments are not the true QCD densities, the main claim does not follow.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing point: the projectile dipole densities of Eq. (17) are an unproven ansatz imported from previous work, and the loop summation inherits this input. I agree with the reader's conditional verdict: the calculation is internally coherent and the qualitative claim is plausible, but the central energy-dependence exponent is not independently secured. The paper's own Conclusions list unproven assumptions, which further supports the reader's caution. I do not see a stronger objection that would require rejection: the derivation is analytic, the j=1 case reproduces the known rare-fluctuation result, and the final S-matrix is positive. The quantitative A^{1/6} bound is additionally sensitive to the unproven C^j(z)→const assumption and to the nuclear profile model, but those weaken only the numerical bound, not the qualitative claim. The sign discrepancy in Eq. (45c) appears to be a typo and does not affect the final summed result, so it is not load-bearing for the central claim. Therefore the read should remain CONDITIONAL, with the request that Eq. (17) be derived or checked against the full evolution equations before the strong version of the claim is accepted.","tokens_in":15788,"tokens_out":11320,"duration_ms":104288,"concrete_test":"Independently derive Eq. (17) from the evolution equation for the generating functional Z in Eqs. (2)-(4), or at minimum verify that the n=2 and n=3 factorial moments obtained from Eq. (17) satisfy the QCD dipole evolution and recurrence relations of Refs. [20,26,55]. If the n=2 moment does not match the BFKL 1→2 dipole splitting from a single initial dipole, Eq. (17) is not the true projectile density and the z'^2/(4κ) exponent in Eq. (49) is an artifact of the ansatz.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The entire loop resummation is controlled by the projectile dipole densities ρ_n^P in Eq. (17), imported from the author's Ref. [58]. These densities are an ansatz: they were chosen so that one particular projection, the S-matrix in Eq. (23a), reproduces the known BK solution (5). The same factorial moments Γ(ω+n)/Γ(ω) then enter the target density Eq. (27), and the saddle-point evaluation leading to Eq. (49) and finally to S(z') ~ exp(−z'^2/(4κ)) in Eq. (52) inherits this construction. Eq. (23a) constrains only the weighted sum with γ^n, so it does not validate the individual n-dipole factorial moments that control the multi-nucleon configurations. If the true QCD dipole densities beyond BK have different correlations, the claimed equality of the dipole-nucleus and dipole-dipole energy dependence would not follow. The paper does not rederive Eq. (17) here, and the Conclusions explicitly concede that the estimates rely on assumptions 'which we cannot prove', including the separate assumption that C^j(z) tends to a constant at large z. The quantitative BK-validity bound z' ≤ 2√(κ c) A^{1/6} is therefore not established beyond the plausibility of this ansatz.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper addresses the summation of large BFKL Pomeron loops in dipole-nucleus scattering in the saturation region. Starting from the t-channel unitarity representation of the scattering amplitude and using dipole densities introduced in the author's previous work (Ref. [58]), the paper computes the S-matrix for configurations with j nucleons at the target, performs a saddle-point evaluation for general j, and sums over j to obtain S(z') = C S_A e^{-C S_A} e^{-z'^2/(4\\kappa)}. Comparing with the BK asymptotic S_BK(z') = \\pi R_A^2 C e^{-z'^2/(2\\kappa)}, it concludes that the loop-summed amplitude has the same energy dependence as dipole-dipole scattering and that BK applies only for z' <= 2\\sqrt{\\kappa c} A^{1/6}. The analytic manipulations for j=1,2,3 and the general saddle point are presented in detail.","tokens_in":16079,"tokens_out":5636,"duration_ms":50842,"significance":"If the result were established, it would be an important step in understanding Pomeron-loop effects in high-energy QCD: it would show that the BK equation misses the dominant loop contribution and that dipole-nucleus and dipole-dipole amplitudes share the same energy dependence in the saturation region, with a quantitative bound on the BK validity range. The paper is transparent about its limitations, and the analytic derivations are sufficiently detailed to be checked. Its main limitation is that the entire construction inherits the multi-dipole factorial moments of Eq. (17) from a previous paper without a derivation here, and the final sum over nucleons relies on an additional unproven smoothness/constancy assumption; until these are justified, the claimed exponent and crossover are conditional.","major_comments":[{"comment":"The projectile dipole densities in Eq. (17) are imported from Ref. [58] and are not derived in this manuscript. The check in Eq. (23a) only shows that the combination with a single \\gamma^n insertion reproduces the known BK solution (5); it does not validate the individual factorial moments \\Gamma(\\omega+n)/\\Gamma(\\omega) that determine the multi-Pomeron correlations entering the target densities (27) and, through the saddle point (48)-(49a), the final exponent z'^2/(4\\kappa). Since the central claim that the dipole-nucleus amplitude has the same energy dependence as dipole-dipole scattering depends entirely on this ansatz, the paper should either derive Eq. (17) from QCD evolution equations for dipole densities or provide an independent test of the higher factorial moments.","section":"II.B, Eq. (17); III.E, Eq. (50)"},{"comment":"The closed-form result S = C S_A e^{-C S_A} e^{-z'^2/(4\\kappa)} is obtained by summing Eq. (49b) over j, which the author explicitly justifies by two unproven assumptions stated in the Conclusions: that the smooth function is C^j(z)=C^j and that C^j(z) tends to a constant at large z. These assumptions are load-bearing: without them the sum over j cannot be performed and the nuclear factor e^{-C S_A}, and therefore the crossover estimate z' \\le 2\\sqrt{\\kappa c} A^{1/6}, do not follow. The manuscript should state these as explicit hypotheses and provide either a derivation or a controlled estimate of the error incurred by replacing C^j(z) by C^j.","section":"IV, Conclusions; Eq. (50)"},{"comment":"The large-z saddle-point evaluation is performed for fixed j, and Eq. (49b) retains only the j'=1 term. The subsequent summation over j requires uniformity of this asymptotic approximation in j, because the alternating series has terms growing like (C S_A)^j/(j-1)! before cancellation. No uniformity argument is provided, so the exponential e^{-C S_A} in Eq. (50) is not rigorously controlled. The author should specify the domain of validity in (j,z') or justify interchanging the large-z and large-j limits.","section":"III.E, Eqs. (48)-(50)"}],"minor_comments":[{"comment":"The constant C in Eq. (17) is treated as a normalization of dipole densities, but Eq. (51) states that C has dimension of cross section; please clarify the dimensions and relations between the various C and C' constants.","section":"II.B and III.E, notation"},{"comment":"The summation limits in Eq. (26) are difficult to parse; please define the ranges of k_1,...,k_j more transparently and state explicitly that each k_i runs at least from 1 to n-j+1 with the constraint \\sum k_i = n.","section":"II.C, Eq. (26)"},{"comment":"Eq. (40) contains \\Gamma(\\omega_3) although only \\omega_1 and \\omega_2 are introduced for the j=2 case; this appears to be a typo and should be corrected.","section":"III.C, Eq. (40)"},{"comment":"The transition from the cylindrical-model expression in Eq. (51) to the simplified form in Eq. (52) should define c in terms of A and R_A explicitly, since the notation c is otherwise new.","section":"III.E, Eqs. (51)-(52)"},{"comment":"The abstract uses z \\le 2\\sqrt{\\kappa C} A^{1/6}, while the body of the paper uses z' and the coefficient c; please align the notation between the abstract and the main text.","section":"Abstract and Eq. (34)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is essentially an application of the author's own previous dipole-density construction, and the lack of an independent derivation makes the main claim fragile. The Conclusions' candid admission that the estimates are 'not very reliable' suggests the paper would benefit either from a derivation of Eq. (17) or from a clear statement that the central result is conditional on that ansatz. As it stands, the framing 'It turns out that...' overstates the status of a result that depends on unproven input assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper is a genuine extension, not a repackaging: it takes the dipole densities from Levin's dipole-dipole paper (Ref. [58]), builds the nuclear target density as a sum over j nucleons each running its own cascade (Eq. 27), and shows the loop-summed dipole-nucleus S-matrix factorizes into a Glauber factor e^{-CSA} times the dipole-dipole result exp(-z'^2/4κ). The new A^{1/6} bound on BK validity is a real prediction. Second, the load-bearing input is Eq. (17), an ansatz for the projectile dipole densities imported from the author's own earlier work and not derived in this paper. Everything — the 1/(4κ) exponent included — inherits it.\n\nCredit where due: the calculation is transparent and checkable. I went through the j=1, 2, and 3 cases; the saddle-point steps are coherent, and the general-j formula (49a), with j'=1 dominating at large z, follows. The author also states plainly in the Conclusions which assumptions he cannot prove and says the quantitative estimates are not very reliable. The result he claims survives those assumptions — the equality with dipole-dipole energy dependence — is the part most robust within the framework.\n\nSoft spots, in proportion. Eq. (17) is the big one. The paper shows the ansatz reproduces the BK solution via Eq. (23a), but that constrains the generating function along one kinematic slice; it does not pin down the individual multi-dipole factorial moments that control the multi-nucleon configurations. The reader's worry that the 1/(4κ) exponent would not follow if the true QCD densities differ is legitimate, and the stress-test concern lands on the same spot the author himself concedes. None of this makes the paper unserious; the weakness is in the input, not the execution. Second, the quantitative BK-validity bound depends on C^j(z) tending to a constant and on the cylindrical nuclear model — the author flags both. Treat the A^{1/6} number as illustrative; the qualitative statement is the contribution. Third, the abstract states the bound as a headline result without the Conclusions' caveats, and it has a C/c notation slip. That mismatch should be fixed.\n\nAudience: people tracking the Pomeron-loop problem in small-x/saturation physics. The paper does not settle whether Eq. (17) is the true QCD density, but it works out the nuclear consequences of that choice and gives a concrete, falsifiable bound. No code or data, so reproducibility means analytic re-derivation. The self-citation is heavy but appropriate for a single-author program. I'd send it to review and let a referee push on Eq. (17) and the abstract alignment. Expected outcome: revised version, not rejection.","headline":"A checkable nuclear extension of Levin's dipole-dipole loop sum that yields a 1/(4κ) energy dependence and an A^{1/6} BK-validity bound — but the input densities in Eq. (17) are an imported, unproven ansatz, and the abstract oversells what the Conclusions concede.","tokens_in":16640,"tokens_out":14645,"would_cite":true,"duration_ms":118914,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.60.Hb","12.38.Cy"],"model":"deepseek-v4-flash","headline":"Summing large Pomeron loops replaces the BK exponential decay of the dipole-nucleus S-matrix, $\\exp(-z'^2/2\\kappa)$, by the slower $\\exp(-z'^2/4\\kappa)$, so the BK equation is valid only for $z' \\le 2\\sqrt{\\kappa c}A^{1/6}$.","keywords":["Pomeron loops","BK equation","dipole-nucleus scattering","saturation region","BFKL Pomeron","t-channel unitarity","geometric scaling","enhanced diagrams"],"falsifier":"Compute the projectile dipole densities $\\rho^P_n$ of Eq. (17) directly from a numerical solution of the full evolution hierarchy without assuming the specific $\\Gamma(\\omega+n)/\\Gamma(\\omega)$ form, and check whether the summed $S$-matrix still switches from the BK-like $\\exp(-z'^2/(2\\kappa))$ at intermediate energies to $\\exp(-z'^2/(4\\kappa))$ at large $z'$; if the factorial moments of the density deviate from the assumed form, the predicted exponent and crossover shift.","tokens_in":15514,"feed_emoji":"⚛️","tokens_out":9412,"duration_ms":71099,"temperature":0.7,"pith_summary":"The paper confronts the Balitsky-Kovchegov (BK) equation, the standard nonlinear description of dipole-nucleus scattering in the saturation regime, with the contribution of large BFKL Pomeron loops that BK leaves out. Its central claim is that after summing these enhanced loop diagrams through $t$-channel unitarity, the $S$-matrix for dipole-nucleus scattering decays with rapidity as $\\exp(-z'^2/(4\\kappa))$ — the same energy dependence as the dipole-dipole amplitude — rather than the BK form $\\exp(-z'^2/(2\\kappa))$. Because the loop-summed amplitude falls more slowly, it overtakes the BK amplitude at $z'\\approx 2\\sqrt{\\kappa c}A^{1/6}$; beyond that point the BK equation no longer captures the scattering. The conclusion matters because it delimits where the otherwise successful BK framework can be trusted, and it makes the asymptotic dipole-nucleus amplitude as universal as the dipole-dipole one.","feed_headline":"Pomeron loops beat BK at ultrahigh energies","feed_subtitle":"Summing large loops changes the energy decay to exp(-z'^2/4κ), limiting BK to z' up to about 2 sqrt(κc) A^(1/6).","key_machinery":"The central machinery is the $t$-channel unitarity representation (Eq. 1) of the scattering amplitude as a sum over dipole densities $\\rho^P_n$ and $\\rho^T_n$, together with the specific projectile densities of Eq. (17), taken from the author's earlier paper, which are engineered so that the fan sum (Eq. 23a) reproduces the analytic BK solution. For the nucleus, each nucleon develops its own dipole cascade, giving the product structure of Eq. (26)-(27); introducing the integral representation $\\Gamma(\\omega+k)=\\int_0^\\infty dt\\,t^{\\omega+k-1}e^{-t}$ turns the combinatorial sums over $k_i$ into closed-form $\\beta$-function terms (Eqs. 39-41, 46). The load-bearing step is the saddle-point evaluation of the $\\omega$-integrals: a configuration with $j$ nucleon cascades, of which $j'$ actually meet the projectile dipoles, contributes $\\exp[-(j'/(j'+1))z'^2/(2\\kappa)]$, and the steepest descent selects $j'=1$.","core_discovery":"The discovery is that the summed large-Pomeron-loop $S$-matrix for a dipole on a nucleus factorizes as $S(z')=C S_A(b)\\exp(-C S_A(b))\\exp(-z'^2/(4\\kappa))$ in the rough cylindrical model, with the same Gaussian-in-$z'$ exponent as the dipole-dipole amplitude. The derivation splits the target into $j$ nucleon cascades; the general term (Eq. 49a) contains exponentials $\\exp[-(j'/(j'+1))z'^2/(2\\kappa)]$ from the saddle points, and the term with $j'=1$ — one nucleon cascade meeting the projectile, all others contributing single dipoles — dominates at large $z'$. This yields the exponent $1/(4\\kappa)$. The BK fan-diagram sum, by contrast, corresponds to $j'=j$ and gives $\\exp(-z'^2/(2\\kappa))$; the two amplitudes cross at $z'\\approx 2\\sqrt{\\kappa c}A^{1/6}$. Since the loop contribution stays larger, BK cannot describe the energy dependence at ultrahigh energy.","pith_inferences":["If the density ansatz is right, the dominance of the $j'=1$ configuration suggests that the single-nucleon-fluctuation channel, rather than the many-nucleon average, controls the asymptotic energy dependence; this may imply a universal high-energy amplitude for any target with the same exponent.","The crossover scale $A^{1/6}$ implies that heavier nuclei delay the breakdown of BK to higher rapidities; an electron-ion collider could look for an $A$-dependent turnover in the energy dependence of the cross section.","The paper does not derive Eq. (17) from first principles; supplying such a derivation, or computing the fate of the assumption that $C^j(z)$ tends to a constant, would be the natural next step and would also open the same method to gluon production and exclusive final states.","One testable consequence of the $j'=1$ dominance is that the contribution of configurations with more than one interacting nucleon cascade is exponentially suppressed at large $z'$; measuring final-state multiplicities or azimuthal correlations might expose which cascade configuration actually dominates."],"forward_implications":["At ultrahigh energies the dipole-nucleus $S$-matrix decays as $\\exp(-z'^2/(4\\kappa))$, staying larger than the BK prediction, so the BK equation underestimates the amplitude beyond $z'\\approx 2\\sqrt{\\kappa c}A^{1/6}$.","The BK equation remains trustworthy only in the limited range $z' \\le 2\\sqrt{\\kappa c}A^{1/6}$; for heavy nuclei this bound is roughly $18+\\ln(Q^2 R_N^2)$, beyond currently accessible energies, so BK stays usable in practice at existing colliders.","The enhanced Pomeron-loop contribution, not the fan diagrams, controls the ultrahigh-energy regime; saturation physics must include loop resummation to maintain $s$-channel unitarity.","The asymptotic dipole-nucleus amplitude acquires the same energy dependence as the dipole-dipole amplitude, so the rare-fluctuation exponent $1/(4\\kappa)$ is universal at very high rapidity."],"supporting_citations":[{"why":"Supplies the projectile dipole densities $\\rho^P_n$ of Eq. (17), the key input the loop summation is built on.","marker":"[58]"},{"why":"Provides the analytic BK solution (Eq. 5) that the densities are engineered to reproduce at the fan-diagram level.","marker":"[25]"},{"why":"Introduces the dipole densities via the generating functional and the unitarity formula Eq. (1) for the scattering amplitude.","marker":"[20]"},{"why":"Gives the earlier rare-fluctuation estimate $\\exp(-z'^2/4\\kappa)$ that the $j=1$ term of this paper recovers.","marker":"[60]"},{"why":"Argues the BK equation violates $s$-channel unitarity, motivating the enhanced-loop summation.","marker":"[22]"},{"why":"Companion argument that non-linear equations need loop modifications in the saturation region.","marker":"[23]"},{"why":"Derives the fan-diagram sum leading to BK and defines the saturation momentum scale used throughout.","marker":"[2]"},{"why":"Defines the BFKL equation and Pomeron Green's function underlying all the dipole cascades.","marker":"[1]"},{"why":"Establishes the coordinate-space BK equation and the form of $t$-channel unitarity used in the amplitude.","marker":"[5]"}],"fun_headline_variants":["Large Pomeron loops break BK at ultrahigh energies","BK only valid up to A^{1/6} scale beyond which loops win","Summing large Pomeron loops: BK fails at high energy","Pomeron loop sums change energy decay, limit BK validity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the projectile dipole densities of Eq. (17), imported from the author's earlier paper, are the true QCD dipole densities beyond the BK approximation; if they misrepresent multi-Pomeron fluctuations, the claimed change from a $1/(2\\kappa)$ to a $1/(4\\kappa)$ energy exponent does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Large Pomeron loops break BK at ultrahigh energies","BK only valid up to A^{1/6} scale beyond which loops win","Summing large Pomeron loops: BK fails at high energy","Pomeron loop sums change energy decay, limit BK validity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2784,"prompt_tokens":880,"completion_tokens":1904,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":1831}},"tokens_in":496,"tokens_out":1904,"duration_ms":13464,"temperature":1.0,"reasoning_tokens":1831,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T15:14:51.307517+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the projectile dipole densities $\\rho^P_n$ of Eq. (17) directly from a numerical solution of the full evolution hierarchy without assuming the specific $\\Gamma(\\omega+n)/\\Gamma(\\omega)$ form, and check whether the summed $S$-matrix still switches from the BK-like $\\exp(-z'^2/(2\\kappa))$ at intermediate energies to $\\exp(-z'^2/(4\\kappa))$ at large $z'$; if the factorial moments of the density deviate from the assumed form, the predicted exponent and crossover shift.","supporting_citations":[{"cited_title":"3: Summing large Pomeron loops for dipole-nucleus scattering","cited_arxiv_id":null,"evidence_quote":"Derives the fan-diagram sum leading to BK and defines the saturation momentum scale used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the BFKL equation and Pomeron Green's function underlying all the dipole cascades."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the coordinate-space BK equation and the form of $t$-channel unitarity used in the amplitude."}],"review_version":1}