{"id":"787a6cb4-19dc-4a26-8f18-703dd4c6aa18","arxiv_id":"hep-ph/0306279","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Solving the Balitsky-Kovchegov equation with impact parameter numerically, the authors find power-law tails in impact parameter that make the cross section violate the Froissart bound despite nonlinear saturation.","lead":"A numerical solution of the Balitsky-Kovchegov equation with full impact-parameter dependence shows that the interaction probability develops power-law tails in impact parameter, so the total cross section grows faster than the Froissart bound allows. The result matters because it shows that saturation alone does not restore unitarity: confinement, which the equation omits, must be added by hand.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The Froissart-violation claim rests on an exponential fit to R_BD(Y) over Y ≤ 11 with no grid-convergence or error analysis, in a regime the authors themselves call pre-asymptotic; the qualitative conclusion is carried mainly by the analytic argument of [19].","rationale":"The paper is a careful, re-implementable numerical study, and its qualitative picture — a kernel-generated power tail in b, a finite saturation window in r, and a black-disc radius growing exponentially in rapidity — is coherent and consistent with the analytic work of [19, 20]. The reader's verdict of CONDITIONAL is defensible and I would not move it. I agree with the reader that the quantitative exponents are uncalibrated, but I disagree with the specific mechanism proposed as the weakest assumption. The reader argues that the Gaussian profile S(b) = 10 exp(−b²/2) (Eq. 15) has support at every b, so the numerics cannot separate kernel-driven from initial-condition-driven tails. The paper largely closes that gap itself: Fig. 6 shows the input falling exponentially while the Y = 0.1 output already falls as 1/b^3.6, and Fig. 7 explicitly decomposes the right-hand side into short-range (|z − b| < r₀) and long-range (|z − b| > r₀) contributions, showing the latter dominates at large b. The analytic estimate in Eq. 22, r²πR₀²(Y)/b⁴, likewise makes the power form a property of the kernel, with only the coefficient set by the target profile. So the deferred companion paper [31] is a reasonable request for a systematic scan, but its absence does not by itself undercut the power-law claim. Where I would sharpen the reader's concern is elsewhere: the Froissart conclusion depends on the asymptotic behaviour of λ_BD, and the extraction window (Y ≤ 11) is one the authors describe as still pre-asymptotic. Combined with the absence of any grid- or step-size convergence study for a computation whose large-b amplitudes span nine orders of magnitude, this is the place where the numerical support is thinnest. That is a calibration gap, not an error, and the analytic argument of [19] independently supports the conclusion, so the verdict should remain CONDITIONAL rather than move toward rejection. Credit is due for the internal cross-check λ_s = 2λ_BD from conformal invariance (Eq. 30), which the quoted numbers do satisfy, and for the explicit lower-bound logic σ ≥ σ_BD that makes the Froissart argument valid even though only the black-disc piece is computed.","tokens_in":16107,"tokens_out":4438,"duration_ms":72538,"concrete_test":"Rerun the evolution to Y = 11 with the ln b grid doubled to 200 points and ΔY halved, and extend the evolution to Y = 14–15 on the finer grid; extract γ and λ_BD from both runs. If λ_BD continues to drift upward beyond the quoted 0.6–1.0 range, or if γ shifts by more than about 0.3, the exponential fit is pre-asymptotic and the numerical basis for abstract:s5 weakens, leaving the claim resting on the analytic argument of [19] alone.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract:s5) requires the total cross section to grow faster than Y². Note first that with the reported tail exponent γ ≈ 2–3 > 2, the integral ∫ b db N converges at large b, so the total cross section is dominated by the black-disc region and σ ≈ 2πR²_BD (Eq. 26). Everything therefore hinges on Eq. 27, R_BD = R₀ exp(λ_BD α_s Y). That exponent is extracted from Y ≤ 11 on a (100 r × 100 b × 20 cosθ) grid with linear interpolation, with no stated error bars, no grid-refinement study, and no ΔY sensitivity analysis. The paper itself states in Sec. 5.3 that λ_BD 'increases rapidly for dipole sizes smaller than the given range above, mostly due to the pre-asymptotic effects' — a direct admission that the fitted exponent has not stabilized. The internal consistency relation λ_s = 2λ_BD (Eq. 30) does hold within the quoted ranges (2 × 0.9–1.0 versus 1.7–2.0), which is genuine support, but it constrains only the ratio, not whether either exponent is asymptotic. Separately, the tail amplitudes used to fix γ reach 10⁻⁹ (Fig. 6), precisely where a 100-point linear-interpolation grid in ln b is coarsest relative to the dynamic range. The numerical evidence is therefore consistent with, but does not by itself establish, asymptotic exponential growth. On the reader's specific framing: the Gaussian profile's support at all b is largely answered by the paper's own Figs. 6 and 7 — the input falls as exp(−b²/2), yet at Y = 0.1 the output already falls as 1/b^3.6, and Fig. 7 separates short- from long-range contributions. The power form is kernel-generated; what remains initial-condition-dependent is the coefficient, not the power.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript presents a numerical solution of the Balitsky-Kovchegov equation retaining the full impact-parameter dependence, using a Glauber-Mueller initial condition with a Gaussian profile S(b) = 10 exp(-b^2/2) and fixed alpha_s = 0.2. The authors report four main results: (i) because the target has finite extent, the amplitude decreases at large dipole sizes r, so the saturated region is bounded in r between 1/Q_s(b,Y) and R_H(b,Y); (ii) the saturation scale acquires a b-dependence whose tail falls as a power Q_s ~ 1/b^gamma with gamma about 1.6-2.0; (iii) a power-law tail N ~ 1/b^gamma with gamma about 2-3 emerges immediately in the evolution, driven by long-range contributions in the kernel (Fig. 7); and (iv) the black-disc radius grows as R_BD = R_0 exp(lambda_BD alpha_s Y) with lambda_BD about 0.6-1.0, so the dipole-nucleus cross section grows faster than the Froissart bound despite the nonlinear saturation term. The authors attribute the Froissart violation to the power tail and note that confinement, absent from the BK kernel, would be needed to restore unitarity.","tokens_in":16495,"tokens_out":9759,"duration_ms":178429,"significance":"If correct, the result is significant: it shows that saturation in the b-dependent BK equation is local in impact parameter but not global, so that unitarity of the total cross section requires physics outside the equation. The paper is a careful and honest numerical study that goes beyond the local approximation of earlier work, and it is consistent with the analytic arguments of Ref. [19]. Strengths include the full three-variable discretization, the explicit iterative scheme in Eqs. (12)-(13), the consistency check lambda_s = 2 lambda_BD from conformal invariance (Eq. 30), and the clean demonstration in Fig. 7 that the long-range kernel contribution dominates at large b. The authors acknowledge the main limitations (no running coupling, no confinement, alternative initial conditions deferred to Ref. [31]).","major_comments":[{"comment":"The headline claim of abstract:s5, that the cross section violates the Froissart bound, rests on the exponential growth of the black-disc radius in Eq. (27), R_BD = R_0 exp(lambda_BD alpha_s Y). The exponent lambda_BD is extracted from Y <= 11 with no stated error bars, no grid-refinement study, and no Delta-Y sensitivity analysis. The authors themselves state in Sec. 5.3 that lambda_BD 'increases rapidly for dipole sizes smaller than the given range above, mostly due to the pre-asymptotic effects,' which is a direct admission that the fitted exponent has not stabilized. A convergence study, varying the grid resolution, the rapidity step, and the fit window, is needed to establish that lambda_BD is asymptotic rather than pre-asymptotic. Without it, the numerical evidence is consistent with, but does not by itself establish, the Froissart-violation claim.","section":"Sec. 5.3, Eq. (27)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The causal attribution of the Froissart violation is imprecise. With the reported tail exponent gamma about 2-3 > 2, the integral of b db N converges at large b, so the total cross section is dominated by the black-disc region and sigma is approximately 2 pi R_BD^2 (Eq. 26); the violation therefore comes from the exponential growth of R_BD, not directly from the power tail. The role of the power tail is to feed the evolution at all b and thereby drive the exponential growth of R_BD, and to ensure that the tail region contributes at the same rate as the black disc rather than being exponentially suppressed. The sentence in Sec. 5.3 that the violation 'is a consequence of the power tails in impact parameter N ~ b^{-gamma}, in the same way as the steep exponential profile ~ exp(-b^2) leads to the linear dependence ~ Y' conflates these two steps. The argument should be made precise: state explicitly that for gamma > 2 the tail integral converges, that the Y-scaling is set by R_BD(Y), and that the power tail is what makes R_BD grow exponentially rather than as a power of Y.","section":"Sec. 5.3 and Conclusions","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The robustness of the extracted exponents gamma and lambda_BD to the choice of initial profile is deferred to Ref. [31], listed as 'work in progress.' Because the Gaussian profile S(b) = 10 exp(-b^2/2) has support at every b, the reader cannot exclude that the asymptotic tail is influenced by the initial condition. Fig. 7 is a good partial answer, since at Y = 0.1 the output already falls as 1/b^3.6 while the input falls as exp(-b^2/2), but a compact-support or steeper profile would settle the question. Either include such a test in this manuscript or state explicitly that the asymptotic exponents may depend on the initial profile and that Ref. [31] will address it.","section":"Sec. 3 and Sec. 5.1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The numerical accuracy criterion is not specified. The text states 'if the relative difference |(N^(1)_xy - N^0_xy)/N^0_xy| < epsilon (some accuracy)' but epsilon is never given, and the convergence of the iterative scheme in Eqs. (12)-(13) is not demonstrated; the statement that 'usually only couple of iterations are needed' is not quantified. Since the headline exponents are extracted from this solution, the accuracy of the solver should be documented, for example by showing that the results are stable under a factor-of-two change in grid resolution and in Delta-Y.","section":"Sec. 3","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"'Kovch egov' and 'Kovner and Wiedeman' should read 'Kovchegov' and 'Wiedemann'; 'considereta tions' in Sec. 1 should read 'considerations.' These are typesetting artifacts but should be cleaned.","section":"Sec. 2 and Sec. 5.1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"For large b there are two solutions of Eq. (19), corresponding to small and large dipole sizes; the paper should state which root is used to define Q_s(b,Y).","section":"Eq. (19)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The tail amplitudes reach 10^-9, at the edge of the 100-point linear-interpolation grid in ln b; a comment on numerical reliability in this region would help.","section":"Fig. 6","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The derivation of lambda_s = 2 lambda_BD via Eq. (29) is heuristic ('roughly'); it should be labeled as a plausibility argument rather than a derivation.","section":"Eq. (30)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that lambda_BD = 0.9-1.0 'is probably closer to the asymptotic one' is speculative; either support it or soften it.","section":"Sec. 5.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Clarify that 'the resulting cross section' is the dipole-nucleus (or gamma*-nucleus) cross section and that the Froissart bound is being applied to it.","section":"Abstract and Sec. 5.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that the amplitude has a maximum at r = 2b should be explained physically (it follows from the geometry of the finite target), not merely noted.","section":"Fig. 3 right panel","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that the crossover between short- and long-range contributions 'is determined by the form of the initial condition' slightly undercuts the 'instant emergence' language; consider rephrasing.","section":"Sec. 5.1","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Worth your time if you work on saturation with b-dependence. This is the first full numerical solution of the BK equation keeping all three variables (ln r, ln b, cos θ), and it delivers three things the local-approximation literature does not have: a finite saturation window 1/Qs(b,Y) < r < R_H(b,Y), a b-dependent saturation scale whose large-b tail is power-like, Qs ~ 1/b^γ with γ ≈ 1.6–2, and a black-disc radius growing exponentially in rapidity. The headline claim at abstract:s5 — that the b-dependent BK cross section violates Froissart despite the nonlinear term — follows from that last number and lines up with the Kovner–Wiedemann analytic argument in [19]. I think it is correct in substance. What is done well: the method is described concretely enough to re-implement (grid, iteration scheme, αs = 0.2, Glauber–Mueller input), and the conformal-invariance consistency check λ_s = 2λ_BD actually holds within the quoted ranges. Fig. 7 is the right diagnostic — it splits short- from long-range contributions and shows the long-range piece owns the power tail. Fig. 6 makes the strongest single point: an exponential input exp(−b²/2) already comes out as 1/b^3.6 after one small step. So the power is kernel-generated; only the coefficient is initial-condition-dependent. I would push back on the worry that the Gaussian profile undermines the claim — the paper answers that itself, and the deferral to [31] matters for the coefficient, not the power. Where it is soft: the exponents are ranges, not numbers with errors. No grid-refinement study, no ΔY sensitivity, and the authors themselves say λ_BD has not stabilized for small r (pre-asymptotic). The tail amplitudes in Fig. 6 reach 1e-9, which is where a 100-point linear grid in ln b is thinnest. So exponential growth is consistent with the data but not established asymptotically. That is a calibration gap, not a flaw in the qualitative conclusion. Send it out; a referee should press on error bars and grid convergence, not on the physics.","headline":"First full three-variable numerical BK solution with impact parameter; the Froissart-violation conclusion is right, but the exponents are quoted as uncalibrated ranges.","tokens_in":17159,"tokens_out":1470,"duration_ms":23090},"pith_extraction":{"msc":[],"pacs":["12.38.-t","13.60.Hb","24.85.+p"],"model":"selfhosted/deepseek-v4.1-flash","headline":"Solving the Balitsky-Kovchegov equation with the full impact-parameter dependence shows that the amplitude develops a power-law tail in impact parameter, so the total dipole-nucleus cross section grows faster than the Froissart bound allows","keywords":["Balitsky-Kovchegov equation","colour dipole model","impact parameter dependence","saturation scale","geometric scaling","Froissart bound","small-x QCD","black disc radius"],"falsifier":"Repeat the evolution with an initial profile of compact support in impact parameter, or with a kernel cutoff that removes the long-range region $|z-b| > r_0$, and check whether the $1/b^{gamma}$ tail and the exponential growth $R_{BD} sim exp(lambda_{BD} alpha_s Y)$ survive; if the tail tracks the support of the input instead, the Froissart-violation claim does not hold.","tokens_in":17450,"feed_emoji":"⚛️","tokens_out":2562,"duration_ms":45998,"temperature":0.7,"pith_summary":"The Balitsky-Kovchegov (BK) equation is the standard mean-field description of how a small colour dipole scatters off a dense nuclear target as energy rises: linear BFKL growth is tamed by a nonlinear term that saturates the amplitude at one. Almost all earlier work solved it while pretending the target is an infinite uniform slab, so the answer depended only on the dipole size $r$. This paper keeps the impact parameter $b$ and solves the equation numerically on a three-dimensional grid in $(r, b, cos theta)$ with a Glauber-Mueller initial profile. The authors find that the exponential fall-off in $b$ they start with is not preserved: the evolution replaces it with a power tail $N sim 1/b^{gamma}$, which makes the black disc radius grow exponentially with rapidity and pushes the integrated cross section past the Froissart bound. If that is right, the nonlinear term enforces unitarity only locally in $b$, and global unitarity has to come from physics the equation does not contain.","feed_headline":"Saturation alone cannot save the Froissart bound","feed_subtitle":"With impact parameter included, the BK amplitude grows a power-law tail in b and the cross section outruns unitarity.","key_machinery":"The load-bearing object is the BK equation itself, the large-$N_c$ closed equation for the two-point dipole amplitude $N(x,y;Y)$, whose kernel describes a parent dipole splitting into two with a measure invariant under Möbius transformations of the transverse plane. The paper exploits two consequences of that symmetry structure. Translation invariance is deliberately broken by the initial condition, which lets the $b$-dependence propagate; and conformal invariance of the kernel implies that the amplitude depends on the combination $r r_0/b^2$, which yields the relation $lambda_s = 2 lambda_{BD}$ between the growth rates of the saturation scale and of the black disc radius, a relation the numerics confirm. The saturation scale is extracted by solving $langle N(r=1/Q_s, b, theta, Y) rangle_theta = 1/2$, and the black disc radius by the same condition with $b$ as the unknown.","core_discovery":"The paper's central claim is that the b-dependent BK equation violates the Froissart bound even though it contains the nonlinear term responsible for saturation. Two numerical facts carry it. First, the amplitude as a function of dipole size is not monotone: because the initial profile is localised in $b$, a very large dipole with both ends outside the target has $N to 0$, so saturation occupies only a finite window $1/Q_s(b,Y) < r < R_H(b,Y)$ whose width grows with rapidity. Second, the steep Gaussian profile in $b$ is washed out immediately, and already after one small rapidity step the amplitude behaves as $N sim 1/b^{3.6}$, softening to $1/b^{gamma}$ with $gamma$ between 2 and 3 at higher rapidity, with the growth at large $b$ matching the hard Pomeron intercept. The authors trace this to long-range contributions in the kernel, where one end of the parent dipole sits in the strong-field region and the other far outside it. Fitting the black disc radius to $R_{BD} = R_0 exp(lambda_{BD} alpha_s Y)$ gives $lambda_{BD} approx 0.6$ to $1.0$, far faster than the linear-in-$Y$ growth an exponential profile would give, hence the bound is violated.","pith_inferences":["The result reframes saturation as a local statement in impact parameter: the nonlinear term caps $N$ at each $b$, but nothing in the equation prevents the black disc from spreading sideways, so global unitarity is an emergent constraint on the kernel rather than a property of the equation.","This suggests a sharp numerical test that the paper defers: rerun with a compact-support profile in $b$ and check whether the $1/b^{gamma}$ tail still appears at the same rate. If the tail survives, it is kernel-driven; if it tracks the input support, the Froissart-violation claim is an artefact of the Gaussian profile.","The same mechanism should show up in the onium-onium Monte Carlo amplitudes, where multiple scattering is included but the target is not a uniform slab, so comparing the extracted $lambda_{BD}$ across those two approaches would test whether the effect is universal or specific to the mean-field truncation.","If confinement-modified kernels restore a faster-than-power fall-off in $b$, the practical consequence is a ceiling on how much of the observed small-$x$ rise can be attributed to BFKL-type evolution at large impact parameter, which is directly relevant to interpreting forward-physics and diffractive data."],"forward_implications":["The saturation region in dipole size is bounded on both sides, $1/Q_s(b,Y) < r < R_H(b,Y)$, and that window widens with rapidity, so geometric scaling in $r Q_s(b,Y)$ holds only for $r$ well below the outer scale $R_H$.","The impact-parameter dependence of the saturation scale factorises at high rapidity, $Q_s(b,Y) = exp(lambda_s alpha_s Y) g(b)$ with $lambda_s approx 1.7$ to $2.0$ nearly independent of $b$, and the large-$b$ tail becomes a power $Q_s sim 1/b^{gamma}$ with $gamma approx 1.6$ to $2.0$ rather than the input Gaussian.","The growth rates satisfy $lambda_s = 2 lambda_{BD}$, so the exponential expansion of the black disc is not an independent parameter but a consequence of conformal invariance of the kernel.","Because the power tail comes from long-range Coulomb-like contributions, restoring the Froissart bound requires breaking the conformal symmetry of the kernel, for instance by a naive cutoff that drops the long-range piece of the integration region."],"supporting_citations":[{"why":"Supplies the BK equation itself, the object being solved numerically.","marker":"[10]"},{"why":"Supplies the Balitsky hierarchy from which the large-$N_c$ closed equation is obtained.","marker":"[7]"},{"why":"Provides the analytical argument that long-range kernel contributions generate power tails in $b$ and violate the Froissart bound, and the comparison value of $lambda_{BD}$ the numerics are checked against.","marker":"[19]"},{"why":"Earlier numerical BK studies in the translationally invariant (no impact parameter) approximation, the baseline the present results are contrasted with.","marker":"[14]"},{"why":"Prior numerical work establishing saturation and geometric scaling in the local approximation, including the running of the saturation scale.","marker":"[16]"},{"why":"Source of the impact-parameter-dependent saturation momentum $Q_s(b,Y)$ used to define and interpret the extracted scale.","marker":"[30]"},{"why":"Glauber multiple-scattering formula that fixes the form of the initial condition $N_0(r,b) = 1 - exp(-r^2 S(b))$.","marker":"[29]"},{"why":"Onium-onium Monte Carlo study with multiple scattering whose exponential growth of the cross section is cited as independent support for the Froissart-violation conclusion.","marker":"[23]"},{"why":"Saddle-point solution of the dipole BFKL equation in transverse space, used to derive the conformal relation $lambda_s = 2 lambda_{BD}$.","marker":"[37]"}],"fun_headline_variants":["Adding impact parameter breaks the Froissart bound in BK","BK numerics: saturation holds only in a finite dipole window","Long-range kernel terms give BK a power-law tail in b","Numerical BK with b-dependence outruns unitarity","Impact-parameter BK: saturation scale grows, bound does not hold"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The demonstration starts from a Gaussian profile in impact parameter that is nonzero everywhere, so long-range contributions are present from the first evolution step, and the paper's own check of whether the kernel rather than the initial condition generates the power tail is deferred to a companion paper listed as work in progress.","fun_headline_variants_meta":{"raw":{"variants":["Adding impact parameter breaks the Froissart bound in BK","BK numerics: saturation holds only in a finite dipole window","Long-range kernel terms give BK a power-law tail in b","Numerical BK with b-dependence outruns unitarity","Impact-parameter BK: saturation scale grows, bound does not hold"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1302,"prompt_tokens":913,"completion_tokens":389,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":316}},"tokens_in":913,"tokens_out":389,"duration_ms":7898,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":316,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T09:08:50.081784+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"Repeat the evolution with an initial profile of compact support in impact parameter, or with a kernel cutoff that removes the long-range region $|z-b| > r_0$, and check whether the $1/b^{gamma}$ tail and the exponential growth $R_{BD} sim exp(lambda_{BD} alpha_s Y)$ survive; if the tail tracks the support of the input instead, the Froissart-violation claim does not hold.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the BK equation itself, the object being solved numerically."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Balitsky hierarchy from which the large-$N_c$ closed equation is obtained."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the analytical argument that long-range kernel contributions generate power tails in $b$ and violate the Froissart bound, and the comparison value of $lambda_{BD}$ the numerics are checked against."},{"cited_title":"Comments on parton saturation at small x in large nuclei","cited_arxiv_id":"hep-ph/0010041","evidence_quote":"Earlier numerical BK studies in the translationally invariant (no impact parameter) approximation, the baseline the present results are contrasted with."},{"cited_title":"Golec-Biernat, L","cited_arxiv_id":null,"evidence_quote":"Prior numerical work establishing saturation and geometric scaling in the local approximation, including the running of the saturation scale."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the impact-parameter-dependent saturation momentum $Q_s(b,Y)$ used to define and interpret the extracted scale."},{"cited_title":"Glauber, Phys","cited_arxiv_id":null,"evidence_quote":"Glauber multiple-scattering formula that fixes the form of the initial condition $N_0(r,b) = 1 - exp(-r^2 S(b))$."},{"cited_title":"Salam, Nucl","cited_arxiv_id":null,"evidence_quote":"Onium-onium Monte Carlo study with multiple scattering whose exponential growth of the cross section is cited as independent support for the Froissart-violation conclusion."},{"cited_title":"Mueller, Nucl","cited_arxiv_id":null,"evidence_quote":"Saddle-point solution of the dipole BFKL equation in transverse space, used to derive the conformal relation $lambda_s = 2 lambda_{BD}$."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}