{"id":"0a91c3bd-e122-4348-9081-a74e13a3df21","arxiv_id":"1908.01621","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors extend the parton-branching TMD approach with CCFM-style non-Sudakov form factors, showing that including small-x dynamics raises the small-x gluon density above the DGLAP result.","lead":"This paper is a first step toward merging two ways of calculating the internal structure of protons at very small momentum fractions. It fits the model to precision scattering data and tests how a small-x correction changes the gluon content of the proton.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's load-bearing step is the unsupported assertion in Sec. 4 that for k⊥ < q_i the 1/z non-Sudakov contribution is already covered by 1/z Δs; if this equality fails, the modified kernels double-count or miss small-x virtual corrections and Fig. 3 loses its meaning.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing step: the unproved factorization of the non-Sudakov form factor in Sec. 4. I agree that this is the central risk. The strongest claim, that PB-based TMDs can be extended to small x with angular ordering and non-Sudakov factors, depends on Eq. (4.4) being a valid rewriting of the CCFM kernels, and the paper does not establish this.\n\nI do not see an additional concern strong enough to replace this one. The LO fit section is exploratory and internally consistent; the χ2 values are reported and are comparable to NLO fits, and the q0 dependence in Figs. 1 and 2 is presented as a numerical study rather than a claim of precision. The QCDNUM agreement in the DGLAP limit gives independent support that the PB baseline is correctly implemented, although the modified kernels are not cross-checked against an independent code.\n\nThe lack of a code or benchmark distribution is a reproducibility concern, but it is secondary to the correctness question. The factorization issue is not demonstrated false by the text, so the appropriate outcome is not rejection but a conditional acceptance pending a derivation and/or numerical verification. Since the reader already reached CONDITIONAL, the verdict remains unchanged.","tokens_in":4599,"tokens_out":6720,"duration_ms":68788,"concrete_test":"Independently derive the product 1/z ~Δns ~Δs for a single branching with k⊥ < q_i, using the definitions in Eq. (4.2), and compare the resulting exponent with 1/z Δs from Eq. (4.3). If the two expressions differ by any term that is not subleading in the soft limit, the 'already covered' statement is false and Eq. (4.4) needs revision. As a numerical companion, run the PB evolution twice at μ^2 = 100 GeV^2: once with Eq. (4.4), i.e. Δns = 1 for k⊥ < q_i, and once keeping ~Δns from Eq. (4.2) in that region; a difference larger than the line width of Fig. 3 in xg(x) for x ∈ [10^-5, 10^-3] would show that the step changes the predictions and must be justified rather than asserted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The extension to small x rests on the transition from the pre-replacement CCFM quantities in Eq. (4.2) to the DGLAP-Sudakov form Δs in Eq. (4.3) and the non-Sudakov factor Δns in Eq. (4.4). The paper states, without derivation, that '1/z ΔnsΔs is already covered by 1/z Δs for k⊥ < q_i' and therefore sets Δns = 1 there. This is the load-bearing step: both the green and blue curves in Fig. 3 are produced with this choice, and the stated conclusion that non-Sudakov corrections keep the small-x gluon density above the DGLAP baseline is only meaningful if this replacement is an identity.\n\nThe definitions in Eq. (4.2) show why the step is not obviously an identity. ~Δns exponentiates the 1/z contribution integrated over q'^2 from z_{i-1} q_{i-1} to k⊥, while Δs in Eq. (4.3) exponentiates a 1/z contribution integrated from z_{i-1} q_{i-1} to q_i. For k⊥ < q_i the former range is a subset of the latter, so the product 1/z ~Δns ~Δs is not automatically equal to 1/z Δs; whether the equality holds depends on the treatment of the z-integration endpoint, the plus-prescription subtraction, and the lower bound z_{i-1} q_{i-1}. If the correct CCFM expression for k⊥ < q_i contains a residual non-Sudakov piece with a different z dependence, or if the 1/z term in Δs is not the same object as the one in ~Δns, the kernels in Eq. (4.1) either miss or double-count small-x virtual corrections. The proceedings text provides no derivation, no numerical cross-check of this specific factorization, and no code release that would allow the equality to be verified independently.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper proposes an extension of the parton-branching (PB) approach to transverse-momentum-dependent (TMD) parton densities so that it covers the small-x regime. It first presents new leading-order PB fits to inclusive HERA DIS data, reporting chi2/dof values of about 1.24 and 1.26, and studies the dependence of the fits on the soft-gluon resolution scale q0. It then modifies the PB splitting kernels by including angular ordering and non-Sudakov form factors following CCFM methods. The main numerical result is that including the non-Sudakov factor suppresses the small-x gluon density relative to the angular-ordering-only result, but the density still remains above the DGLAP-based baseline.","tokens_in":5014,"tokens_out":4362,"duration_ms":46092,"significance":"If the central factorization step in Sec. 4 is correct, the paper offers a practical route to incorporate CCFM small-x dynamics into the PB TMD framework, which would be useful for semi-hard QCD phenomenology at the LHC and future colliders. The paper's strengths include concrete LO fits with reasonable chi2/dof, a systematic step-by-step numerical study of the kernel modifications, and comparison with QCDNUM in the DGLAP limit. However, the load-bearing technical step, Eq. (4.4), is stated without derivation or numerical validation, and the quantitative small-x results are presented without uncertainty estimates. The contribution is therefore best viewed as an exploratory proposal whose central claim requires further support before the method can be considered established.","major_comments":[{"comment":"The statement that '1/z ΔnsΔs is already covered by 1/z Δs for k⊥ < qi', which leads to setting Δns = 1 in that region, is load-bearing and is not derived. The integration ranges in Eq. (4.2) differ: Δns integrates over q'^2 from z_{i-1}q_{i-1} to k⊥, while Δs in Eq. (4.3) integrates over q'^2 from z_{i-1}q_{i-1} to qi. For k⊥ < qi the product 1/z Δns Δs contains an integral over the interval [k⊥, qi] with the 1/(1-z) kernel but not with the 1/z kernel, so the claimed equality with 1/z Δs is not an identity unless additional conditions on the z-integration endpoint and plus-prescription are imposed. Please provide a derivation of this replacement or a numerical cross-check comparing results obtained with the full Eq. (4.2) product against those obtained with Eq. (4.4), since Fig. 3 and the central qualitative conclusion depend on this step.","section":"Sec. 4, Eqs. (4.2)-(4.4)"},{"comment":"The modified splitting functions in Eq. (4.1) involve singular 1/z and 1/(1-z) terms multiplied by form factors, but the paper does not specify the plus-distribution prescription or the precise z-integration limits used in the numerical implementation. Without this information the kernels are not uniquely defined and the results are not reproducible. In particular, the lower limit zM of the z integrals in Eq. (4.2) needs to be specified explicitly for both Δs and Δns, including how zM depends on q0 and on the relevant transverse momentum scale in each case.","section":"Sec. 4, Eq. (4.1)"},{"comment":"The numerical evidence for the small-x behavior is presented without uncertainty bands, and the curves are obtained from a benchmark starting distribution rather than from the LO fitted densities. As a result, the comparison among the red, green, and blue curves in Fig. 3 does not establish whether the differences are significant relative to input-distribution and parameter uncertainties. Please quantify the theoretical uncertainty, for example by varying q0, αs, and the starting distribution, or explicitly label these curves as exploratory illustrations without quantitative claims.","section":"Sec. 4, Fig. 3"}],"minor_comments":[{"comment":"The LO fits are characterized only by chi2/dof values; showing a comparison of the fitted distributions to the HERA data or to the NLO PB results would make the quality of the fits more transparent.","section":"Sec. 2"},{"comment":"Figure 1 shows the q0 dependence of the parton densities without uncertainty bands, so it is difficult to judge whether the visible small-x enhancement is statistically significant; adding fit uncertainties would strengthen the discussion.","section":"Sec. 3, Fig. 1"},{"comment":"The notation 'q,i−1' in the sentence 'DGLAP ordering ( qi > q,i−1)' appears to be a typo for q_{i-1}; the subscript notation should be made consistent throughout.","section":"Sec. 4"},{"comment":"The axes and labels in Fig. 3 are garbled (for example, the x-axis tick labels and the ratio panel), and the figure should be reformatted so that the curves are readable and the ratio is clearly defined.","section":"Sec. 4, Fig. 3"},{"comment":"The paper says 'We observe that 1/z ΔnsΔs is already covered by 1/z Δs for k⊥ < qi' without an explicit reference to the original CCFM literature for this replacement; citing the specific equations in Refs. [3,4,5] would help the reader check the validity of the step.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"This is a DIS2019 proceedings paper, so the exploratory character is understandable. The main concern is that the technical step leading to Eq. (4.4) is both central and unsupported; I would recommend requiring a derivation or a numerical validation before the method is cited as established. The paper's citation pattern is appropriate for the field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — this is a DIS2019 proceedings note from the PB TMD group. The actual new content is modest and clearly scoped: new LO PB fits to HERA1+2 data with chi2/dof around 1.25, a scan over the soft-gluon resolution scale q0, and a first implementation of CCFM-style non-Sudakov form factors in the PB evolution kernels (Eqs. 4.1-4.4). The step-by-step numerical comparison in Fig. 3 (DGLAP, full angular ordering, and angular ordering plus non-Sudakov) is genuinely informative and consistent with CCFM expectations: angular ordering raises the small-x gluon, non-Sudakov suppresses it, and the net stays above DGLAP. The reader's verdict is fair: conditional, moderate confidence.\n\nThe soft spot is exactly where the stress-test points. The transition from the CCFM form factors in (4.2) to the DGLAP Sudakov in (4.3)-(4.4) is load-bearing, and the sentence '1/z Delta_ns Delta_s is already covered by 1/z Delta_s for k_perp < q_i' is an assertion, not a derivation. It may well be right—if the 1/z term in the DGLAP Sudakov already sums the virtual corrections for k_perp < q_i, setting Delta_ns = 1 there avoids double counting. But the paper does not show the z-integration endpoints or plus-prescription details, and the definitions in (4.2) do not make the equality self-evident. Since both the green and blue curves in Fig. 3 depend on this choice, and the paper's central message is that the non-Sudakov correction keeps the gluon above the DGLAP baseline, that statement is only as solid as this factorization. No uncertainty bands on Fig. 3, no code release, no numerical cross-check of the factorization itself. Those omissions are more forgivable in a proceedings than in a full paper, but they are real.\n\nThe LO fit part is more solid. The fits are done with xFitter, the chi2 values are reasonable, and the q0 dependence is shown. That part is reproducible in principle and clearly documented. I do not think the paper overclaims—the conclusion says 'we have presented a method to incorporate CCFM effects into the PB formulation,' and 'method' is appropriate given the missing derivation.\n\nWho does this actually help? People working on PB TMDs, TMD fits, and small-x evolution. It is a useful status report and a pointer to what a longer publication should contain. I would not cite it as a definitive result, but I would bring it to a reading group interested in the small-x TMD interface.\n\nRecommendation: if this comes to a journal, send it to peer review—not to reject it, but to force the derivation or a numerical verification into the paper. As a proceedings, it is borderline but acceptable, provided the reader understands it is exploratory. My verdict matches the reader's: conditional.","headline":"A genuinely new but unproven step in PB TMD small-x extension; the key non-Sudakov factorization in Sec. 4 is asserted rather than derived, and Fig. 3 rests on it.","tokens_in":5612,"tokens_out":2806,"would_cite":false,"duration_ms":26877,"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":"The parton-branching TMD framework can be extended to small $x$ by adding angular ordering and non-Sudakov form factors to its evolution kernels.","keywords":["parton branching","TMD parton distribution functions","small x","CCFM","Sudakov form factor","non-Sudakov form factor","soft-gluon resolution scale","DGLAP evolution"],"falsifier":"A decisive check would be to run an independent numerical solution of the CCFM evolution that does not set the non-Sudakov form factor to one for $k_\\perpq_i$, and $\\Delta_{ns}=1$ for $k_\\perp