{"id":"3ebd685c-430b-4504-9ff3-26b5f80adab8","arxiv_id":"2601.12489","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A branching parton model with fixed splitting probability w gives xf(x)∼x^{-δ}, δ=ln(2w)/ln2, and, if an ad hoc freeze is imposed, a saturated parton density at very small x.","lead":"This toy model treats a fast proton as a cascade of partons that split in half and sometimes merge, and derives a power-law growth of the parton density at small momentum fraction x. The paper also claims to find saturation of the density at very small x, but that part rests on an added stabilization assumption rather than on the model's equations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Saturation is not derived: §4's n_s estimate uses eN_n≈(2w)^n although the fusion term in Eq (7) is O(1/v) there, and the plateau is imposed by hand after the model disclaims knowledge beyond n_s.","rationale":"Read in good faith, the paper's power-law section is internally consistent: under the stated equal-halving and independent-splitting assumptions, Eq. (6) follows from (5), and the momentum sum (1) checks. The paper is explicit that w is not determined by the model. The weak point is the saturation announcement. The derivation of the saturation scale requires knowing \\tilde N_n, but the model never solves (7); the paper replaces \\tilde N_n by its no-fusion value (2w)^n in Eq. (20). This replacement is invalid at the relevant scale because the fusion term in (7) is then O(1/v), not negligible. And the plateau is not a consequence of any equation: the model's own evolution equation gives a turnover and decline, and stabilization is posited. Since the abstract and conclusions advertise saturation as the main non-perturbative finding, the central claim is unsupported. The reader's rationale already identified this; the reader's declared weakest assumption (equal halving) is a different, less fatal issue. Hence the rejection stands; I would keep the verdict unchanged.","tokens_in":937,"tokens_out":2379,"duration_ms":160863,"concrete_test":"Numerically solve the boundary-value problem (7)-(11) for w=0.54, v=0.01 and v=0.001, \\bar n=13, using the fusion sum (9) (and, for comparison, its leading term v C_N). From the solution \\tilde N_n, compute N_n via (12). Then check: (i) is \\tilde N_n at the n_s defined by \\tilde N_n = (1-w+v)/(2v) close to (2w)^{n_s}? (ii) does N_n reach (21) and remain flat for n>n_s, or does it peak and decline toward zero (as Eq. (18) suggests)? If \\tilde N_n differs significantly from (2w)^{n_s}, or if no plateau appears without imposing 'complete stabilization', Eqs. (20)-(22) fail.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central advertised result is saturation: a plateau dN/dy = (1-w+v)^2/(4v) for x ≤ x_s ~ v^{1/δ_w}, Eqs. (21)-(22). This is not derived from the model. Eq. (20) locates n_s by inserting the no-fusion form \\tilde N_n → (2w)^n into (19). At that n_s, \\tilde N_n ~ 1/v, and in the exact recurrence (7) the fusion source N(v,\\tilde N_{n+1}) ~ v \\tilde N_{n+1}^2/2 ~ 1/(2v), the same parametric order as the splitting term 2w \\tilde N_{n-1}. So the approximation \\tilde N_n ≈ (2w)^n is not uniformly valid up to n_s; the saturation scale is effectively assumed, not computed. Furthermore, even granting the maximum, Eq. (18) is a parabola: N_n increases up to n_s and then decreases to zero. The paper then adds 'the simplest option is complete stabilization' and admits 'our model provides no information about the behavior of the system' beyond n_s. The plateau is superimposed by hand. Hence 'saturation is detected' overstates the mathematical content.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a two-parameter parton model for the small-x PDF in the non-perturbative region. Partons are organized into generations; each splitting is assumed to divide the parent's longitudinal momentum exactly in half, so generation n has x_n = 2^{-n}. With a constant splitting probability w and no fusion, the mean occupation of generation n is N_n=(1-w)(2w)^n, leading to x f(x)=(1-w)x^{-\\delta_w} with \\delta_w=\\ln(2w)/\\ln 2. Parton fusion with probability v is then added through the nonlinear recurrence Eq. (7). The paper derives an upper bound from probability conservation, analyzes a truncated balance equation, and proposes that the density saturates at n_s with x_s~v^{1/\\delta_w} and dN/dy~(1-w+v)^2/(4v), after which complete stabilization is assumed. The final section compares these results qualitatively with BFKL, GLR, and BK expectations.","tokens_in":9692,"tokens_out":14185,"duration_ms":144104,"significance":"If the saturation derivation were valid, the paper would offer a simple, non-perturbative argument that small-x saturation is a parton-density phenomenon rather than a virtuality-scale phenomenon. The pure-splitting part is internally consistent: momentum conservation is explicitly checked in Eq. (1), and the conversion from discrete generations to a continuous PDF is transparent. The paper is also candid in places: it states that Eq. (7) is not solved and that the model provides no information beyond n_s. However, the central advertised result is not actually derived. The power-law exponent is essentially a reparametrization of the assumed binary-splitting rule, the saturation scale is obtained by inserting the no-fusion solution into a relation where fusion is not negligible, and the plateau is imposed by hand. These are load-bearing defects, so the paper cannot be recommended for publication in its present form.","major_comments":[{"comment":"The saturation scale is not derived. Equation (19) defines n_s by eN_{n_s}=(1-w+v)/(2v), but eN_n is the solution of Eq. (7), which the paper states is unknown. The text then says “Assuming that eN_{ns}→(2w)^{ns} as v→0,” i.e. it substitutes the no-fusion form into (19). At the resulting n_s one has (2w)^{n_s}∼1/v. In Eq. (7) the fusion source is at least v eN_{n+1}(eN_{n+1}-1)/2 ∼ 1/(2v), the same parametric order as the splitting term 2w eN_{n-1}. Thus the no-fusion approximation is not uniformly valid up to n_s, and x_s∼v^{1/δ_w} is an assumption, not a consequence of the recurrence.","section":"§4, Eqs. (19)–(20)"},{"comment":"The saturation plateau is imposed, not derived. Eq. (18) gives N_n as a parabola in eN_n; once eN_n exceeds eN_{ns}, N_n decreases to zero. The text calls this decreasing option “formally consistent” but unsatisfactory, then selects “the simplest option is complete stabilization,” and explicitly admits “our model provides no information about the behavior of the system” beyond n_s. Therefore the constant dN/dy≈(1-w+v)^2/(4v) in Eq. (22) is an external stabilization hypothesis. The Abstract’s statement that “the phenomenon of saturation of the parton density is detected” overstates what the model establishes.","section":"§4, after Eq. (18); Abstract"},{"comment":"The power-law exponent is built into the discrete kinematics. Since x_n=2^{-n} and N_n=(1-w)(2w)^n by construction, substituting n=ln(1/x)/ln 2 gives x f(x)=(1-w)x^{-δ_w} with δ_w=ln(2w)/ln 2. The numerical comparison in §2 then fixes w by assuming δ_w=0.1 or 0.3 from the soft/BFKL Pomeron. Thus the model does not independently predict the small-x slope; it re-expresses an assumed input. This does not invalidate the pure-splitting kinematics, but it substantially reduces the force of the first advertised result.","section":"§2, Eq. (6)"},{"comment":"The approximations leading to Eq. (21) are not controlled. Eq. (17) discards ΔN with the justification that ΔN=O(v^2). For occupation numbers near the would-be maximum, eN_n=O(1/v), and the discarded terms are v^2 C_{eN_n-2}+⋯=O(1), not O(v^2) uniformly. While this is subleading relative to N_{ns}=O(1/v), it means that Eq. (21) is obtained without a quantitative error estimate. More importantly, the preceding substitution of the no-fusion form for eN_n into (19) is the step that actually determines n_s, and that substitution is unjustified, as noted above.","section":"§3–§4, Eqs. (17)–(21)"}],"minor_comments":[{"comment":"“PDF at small xin” and similar spacing issues in the abstract and body should be corrected.","section":"Abstract"},{"comment":"“in the letter case” should be “in the latter case.”","section":"§4"},{"comment":"“parton distribution fusion” appears to be a typo for “parton distribution function.”","section":"§5"},{"comment":"The notation C_N is ambiguous: it is defined as N(N-1)/2, but the compact sum uses C_{N-2k}. Spell out that C_m is the binomial coefficient m(m-1)/2 or write inom{m}{2} for clarity.","section":"§3, Eqs. (8)–(9)"},{"comment":"The mixed notation eN_{ns} and n_s should be clarified; it is easy to misread eN_{ns} as eN_n times s.","section":"§4"}],"recommendation":"reject","confidential_remarks":"The paper is clearly written and the pure-splitting part is sound, but the central saturation claim is asserted rather than derived. The abstract and conclusions are stronger than the mathematical content of §4 supports. I do not see a local edit that would make the saturation result valid; it would require either solving or numerically controlling Eq. (7) or reframing the paper as a toy model with an externally imposed plateau, which would remove the main advertised result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe pure-splitting part of this paper is a tidy exercise: with the explicit equal-halving assumption, the average occupation number (2w)^n gives x f(x) ~ x^{-δ_w}, δ_w = ln(2w)/ln 2, and momentum conservation checks out exactly. That part is internally consistent and clearly presented. The author also deserves credit for stating plainly that w is a free parameter and that fixing it from pomeron intercepts is illustrative, not predictive.\n\nThe soft spot is the saturation claim, and it is load-bearing. Section 4 does not solve the nonlinear recurrence (7). It drops ΔN to get the parabola (18), locates the maximum at \\tilde N = (1-w+v)/(2v), and identifies the corresponding n_s by inserting the no-fusion form \\tilde N_n ≈ (2w)^n. The problem is that at that n_s, \\tilde N is O(1/v), and the fusion term in (7) is also O(1/v) — same parametric order as the splitting term. So the no-fusion approximation is not valid up to n_s, and the saturation scale is effectively assumed, not computed. After the maximum, the text admits 'our model provides no information about the behavior of the system' and then imposes the plateau as 'the simplest option.' Calling that 'saturation is detected' overstates what has been shown. The stress-test note is right on this.\n\nThere is also a circularity in the power-law exponent: δ_w is defined via w, and w is then chosen to reproduce δ_w = 0.1 or 0.3 from pomeron intercepts. So the slope is an input, not an output. The equal-halving assumption is strong; asymmetric splittings would change the x-dependence. The paper acknowledges some of this but not the circularity.\n\nWhere that leaves the paper: as a toy model of the branching cascade, it is a clean and honest demonstration, though the branching multiplicity itself goes back to Polyakov. As a derivation of saturation, it fails. The qualitative discussion of why saturation need not be tied to Q^2 is worth thinking about. But the central advertised result is unsupported.\n\nI would bring this to a reading group only as an example of overclaiming from a toy model. I would not cite it in my own work. For peer review: I would desk reject it in its current form. If the author replaces the assumed plateau with a controlled treatment of (7) — even a numerical solution — then it would deserve another look.","headline":"Tidy branching-cascade toy for the small-x power law, but the saturation result is assumed, not derived.","tokens_in":10108,"tokens_out":4034,"would_cite":false,"duration_ms":40073,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives a power-law small-x parton density from a binary splitting cascade, with the exponent fixed by the splitting probability, and shows that adding parton fusion saturates the density at a value set by the fusion probability.","keywords":["parton distribution function","small-x physics","parton model","branching cascade","parton fusion","parton saturation","non-perturbative QCD","generation ladder"],"falsifier":"Measure the gluon density x f(x) over a range of small x where fusion is negligible; the model requires x f(x/2)/(x f(x)) = 2w to be constant for every halving in x. If this ratio drifts with x instead of staying flat, the fixed-exponent power law fails; likewise, if the per-rapidity density dN/dy keeps growing with decreasing x past the predicted saturation point instead of slowing to a plateau, the fusion-saturation mechanism is falsified.","tokens_in":9148,"feed_emoji":"⚛️","tokens_out":6614,"duration_ms":69908,"temperature":0.7,"pith_summary":"The paper tries to show that the small-x parton density of a fast-moving proton can be understood without perturbation theory, just by counting partons in a branching cascade. In the cascade, each splitting divides longitudinal momentum in half and occurs with probability w, so the average count in generation n is (2w)^n; translating generation number into momentum fraction gives x f(x) = (1-w)x^{-δ_w} with δ_w = ln(2w)/ln2. Adding fusion of same-generation partons with probability v makes the recurrence nonlinear and forces the density per unit rapidity to saturate at (1-w+v)^2/(4v) once x reaches about v^{1/δ_w}. If the derivation is sound, the small-x exponent and the saturation plateau are tied to two model parameters, and saturation arises from a critical density rather than from a virtuality scale. A sympathetic reader would care because this offers a non-perturbative, mechanism-level explanation of saturation that can be compared with data.","feed_headline":"Splitting probability sets the small-x density exponent","feed_subtitle":"Adding parton fusion caps the density at a saturation plateau whose height is fixed by the fusion probability","key_machinery":"The binary generation ladder: partons are indexed by generation n, with longitudinal fraction x_n=2^{-n}, and the conversion identity Δx_n=x_n Δn turns discrete occupation numbers into the continuous relation x f(x)=N_n. The fusion term in the balance equation and the parabolic expression for the number remaining per generation identify the maximum-density point and the plateau height.","core_discovery":"In the author's model, partons are organized by generations: the parent proton parton is generation 0, and every splitting produces two daughters in the next generation, each carrying exactly half the parent's longitudinal momentum. With splitting probability w, the average number of partons that appear in generation n is (2w)^n, and after the 1-w fraction that do not split is removed, the average number remaining is (1-w)(2w)^n. Because x_n=2^{-n} and the interval Δx_n equals x_n for one generation step, the continuous PDF satisfies x f(x_n)=(1-w)(2w)^n, hence x f(x)=(1-w)x^{-δ_w}, δ_w=ln(2w)/ln2, over the moderate-small-x region. When fusion of two partons in the same generation is include","pith_inferences":["The author does not connect the splitting probability w to a measured quantity; a testable extension would identify w with the effective intercept extracted from deep-inelastic data and then predict the high-energy multiplicity exponent δ_w=ln(2w)/ln2, which could be checked in proton-proton and electron-proton data.","Because the model treats partons as longitudinally identical and ignores transverse structure, it predicts the same saturation height for protons of different transverse size; future measurements of the small-x gluon density in protons versus larger nuclei could test this density-only mechanism against a transverse-area mechanism.","The discrete identity x f(x_{n+1})/x f(x_n)=2w gives a direct experimental falsifier: measuring the gluon density ratio at x and x/2 across a decade should be a constant; any x-dependence would force the model to allow asymmetric splitting.","If fusion is stronger than assumed, the number remaining per generation rises and then falls to zero; the paper leaves open whether the system then stabilizes into a saturated state or collapses, an explicit model extension not settled here."],"forward_implications":["With w>1/2, the total parton multiplicity grows as (P/μ)^{δ_w}, a power law in momentum; for w=1/2 it reduces to logarithmic growth ~ln(P/μ), marking the transition to branching cascades.","In the moderate-x regime, x f(x) is a pure power law, so f(x)~x^{-1-δ_w}; the exponent is tied to the splitting probability and lies in the range -1 < -δ_w < 0.","When fusion is included, saturation appears at x_s∼v^{1/δ_w}, with saturated density per unit rapidity (1-w+v)^2/(4v), both controlled by the fusion probability v.","The initial proton momentum needed for slow partons to form a saturated medium is estimated as P_s/μ∼v^{-1/δ_w}.","The model reproduces, at a qualitative level, the saturation phenomenon predicted by previous perturbative QCD analyses, but explains it as a density effect rather than a virtuality effect."],"fun_headline_variants":["Parton splitting odds set small-x exponent","Fusion creates saturation plateau in parton density","Nonperturbative PDF: exponent from splitting, cap from fusion","Small-x density: exponent from splitting, saturation from fusion"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The derivation assumes every splitting divides the parent's longitudinal momentum exactly in half, so that all partons of generation n share x=2^{-n} and the conversion Δx_n=x_n Δn holds; if real splittings are asymmetric or momentum fractions spread within a generation, the power exponent δ_w and the saturation scale change.","fun_headline_variants_meta":{"raw":{"variants":["Parton splitting odds set small-x exponent","Fusion creates saturation plateau in parton density","Nonperturbative PDF: exponent from splitting, cap from fusion","Small-x density: exponent from splitting, saturation from fusion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1259,"prompt_tokens":681,"completion_tokens":578,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":425,"completion_tokens_details":{"reasoning_tokens":514}},"tokens_in":425,"tokens_out":578,"duration_ms":6253,"temperature":1.0,"reasoning_tokens":514,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:46:42.808425+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gluon density x f(x) over a range of small x where fusion is negligible; the model requires x f(x/2)/(x f(x)) = 2w to be constant for every halving in x. If this ratio drifts with x instead of staying flat, the fixed-exponent power law fails; likewise, if the per-rapidity density dN/dy keeps growing with decreasing x past the predicted saturation point instead of slowing to a plateau, the fusion-saturation mechanism is falsified.","supporting_citations":[],"review_version":1}