{"id":"5b14aa08-f5ef-474f-8013-216955241427","arxiv_id":"2505.23487","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The NLO factorized QED contribution to the leading e+q->e+X subprocess of DIS is computed and shown to be infrared safe with no free parameters beyond the factorization scale.","lead":"This paper calculates the first full next-to-leading order QED correction to the hard scattering part of deep inelastic scattering, treating photon radiation from electrons and quarks on the same footing. It offers a way to include these radiative effects in future precision measurements at the Electron-Ion Collider without ad hoc cutoffs.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (20)'s QED gauge link omits the quark fractional charge e_q, so the claimed operator definition of the quark PDF/LFF is not QED gauge invariant; the factorization theorem underpinning the IR-safe subtraction is not established as written.","rationale":"After reading in good faith, the paper's strongest claim is that the NLO factorized QED hard part (Eq. 55) is IR safe and parameter-free. That claim must be read together with Eq. (36), which defines bH^(3,0) as a difference between the full partonic cross section and a set of convolution subtractions. These subtractions are exactly the collinear terms generated by the joint factorization theorem (Eq. 6). If the theorem is valid, the calculation is a self-contained demonstration of IR safety (with the usual caveat that the lengthy Appendix A coefficients are not independently checked). If the theorem is not valid, the subtraction terms are arbitrary and the hard part is not the physical coefficient. The paper asserts the theorem's proof is 'effectively the same' as Ref. [18], but does not supply the proof. Stress-testing the asserted theorem uncovers a specific problem: Eq. (20) defines the QED part of the gauge link with the elementary charge e, not the quark's charge e e_q. For the electron channel this is fine, but for quark PDFs/LFFs it fails the QED gauge-invariance test. This is not merely an aesthetic issue: the cancellation of collinear singularities in the factorization proof relies on the gauge-link phase exactly matching the charge of the field whose distribution is being defined. The 'same proof' claim therefore cannot be accepted as written. Other concerns—absence of an independent check of the coefficients, mismatch between CT18 PDFs and the joint-scheme PDFs, and the impossibility of using perturbative LDF/LFF with NLO hard parts (Sec. V)—are real but secondary; they affect the numerical illustrations, not the analytic IR-safety statement. The reader's CONDITIONAL verdict is appropriate: the core derivation is plausible and the paper is explicit about its limitations, but the factorization foundation must be either proved or corrected, and Eq. (20) must be fixed to include e_q or justified. No change to the verdict.","tokens_in":28681,"tokens_out":15057,"duration_ms":142268,"concrete_test":"Perform a U(1) gauge variation of Eq. (19): set ψ_q(z) → e^{i e e_q θ(z)} ψ_q(z), A_γ^+ → A_γ^+ + ∂^+ θ, and require the Wilson line in Eq. (20) to cancel the phase between ψ_q(0) and ψ_q(z^-). Verify that the factor e^{-i e e_q ∫ dη A_γ^+} is required; then compute the one-loop f_{q/q}^{(1,0)} from this gauge-invariant operator in a covariant gauge and compare the 1/ϵ coefficient with Eq. (38). If the coefficient changes, Eq. (20) must be corrected; if it does not, the simplified link might be a gauge artifact, but the onus is on the paper to demonstrate that.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that bH^(3,0) in Eq. (55) is completely IR safe—is conditional on the joint QED-QCD factorization theorem in Eq. (6). The paper does not prove this theorem; it asserts the proof is the same as in Ref. [18]. A concrete crack appears in the operator definition of the quark distribution (Eq. (19)) and quark-to-electron fragmentation function (Eq. (44)): the gauge link in Eq. (20) is written as P exp[-i∫dη (g_s A_g^+ t^a - e A_γ^+)], with the electron charge e but no quark fractional charge e_q. Under a QED gauge transformation the quark field ψ_q picks up a phase e^{i e e_q θ}; the Wilson line must cancel this phase to make the bilocal operator gauge invariant. As written, it cancels only a charge-1 phase, so the defined PDF/LFF are not gauge invariant for quarks. Because the factorization theorem and the universality of the subtracted collinear singularities rely on gauge-invariant operator definitions, the 'same as Ref. [18]' claim is not justified. If the gauge link should contain e e_q, then Eq. (20) is wrong as stated; if the authors intend a different convention, they need to state it. Without a correct gauge-invariant definition, the subtraction terms in Eq. (36) may not remove the correct collinear divergences, and bH^(3,0) would not be the true short-distance coefficient. This is distinct from (but more fundamental than) the absence of an independent numerical check of the lengthy coefficients in Appendix A.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first next-to-leading-order (NLO) factorized QED hard coefficient for the leading partonic subprocess e+q -> e+X of inclusive lepton-hadron DIS in the joint QED and QCD factorization framework. Section III derives the NLO QCD hard part bH^(2,1) by converting standard structure-function results into the joint-factorization language. Section IV computes the NLO QED hard part bH^(3,0), given in Eq. (55) with coefficients in Appendix A, by subtracting, order by order in alpha_em, the collinear contributions absorbed into lepton distribution functions (LDFs), lepton fragmentation functions (LFFs), and quark/photon PDFs. The paper argues that all collinear and pinch singularities cancel, so that bH^(3,0) is infrared safe and depends only on the factorization scale. Section V provides numerical estimates of the impact of collision-induced QED radiation at JLab and EIC kinematics, using CT18 PDFs and model LDF/LFF input shapes.","tokens_in":29047,"tokens_out":10964,"duration_ms":131766,"significance":"If the calculation is correct, this is a useful first step toward treating collision-induced QED radiation in inclusive DIS with the same factorization logic as QCD, avoiding ad hoc radiative-correction cuts. The organization of Eq. (55) by powers of the quark charge e_q is physically illuminating, and the explicit demonstration that the photon-PDF subtraction removes the q^2 -> 0 pinch is valuable. The paper also makes a legitimate point that the hard scale for the factorized formula is the observed-lepton transverse momentum l'_T, not Q^2, and it is candid about the nonperturbative nature of LDFs and LFFs. The manuscript is not yet fully supported, however: the operator definitions of the PDFs and LFFs have a gauge-invariance problem, the central NLO QED coefficient is presented without derivation or cross-check, and the factorization theorem on which the IR-safety claim rests is only cited, not established.","major_comments":[{"comment":"The quark PDF in Eq. (19) and the quark-to-electron LFF in Eq. (44) are not QED gauge invariant as written. Under a U(1) transformation the quark field transforms with a phase involving the fractional quark charge e_q, whereas the QED part of the gauge link in Eq. (20), P exp[-i integral (g_s A_g^+ t^a - e A_gamma^+)], is constructed to cancel only a unit-charge phase. The same defect carries over to Eq. (44), whose gauge link is also taken from Eq. (20). Because the subtraction terms in Eq. (36) are supposed to represent the collinear divergences of gauge-invariant, universal distributions, the factorization underpinning the IR-safety claim of Eq. (55) is not established as stated. Please insert the quark fractional charge e_q in the QED part of the link, or state an explicit charge convention and demonstrate that the bilocal operators are invariant under both QCD and QED gauge transformations.","section":"Sec. II, Eqs. (19) and (20); Sec. IV, Eq. (44)"},{"comment":"The central result of the paper, Eq. (55), is presented with the coefficients a_i, b_i, and c_i relegated to Appendix A and with no derivation shown. The text does not describe how the d-dimensional phase-space integrals for the diagrams in Figs. 11 and 12 are evaluated, how the 1/epsilon poles of sigma^(3,0) are distributed among the subtraction terms in Eq. (36), or how the cancellation between the two-photon-exchange diagrams and the real-interference diagrams is verified. No independent cross-check is provided against known QED radiative-correction results, even in a limiting case. As a calculation paper, this leaves the main claim unsupported. Please provide a derivation summary and at least one nontrivial consistency check, such as a demonstration of the 1/epsilon pole cancellation order by order or a comparison with a standard radiative-correction expression in a well-defined limit.","section":"Sec. IV, Eq. (55) and Appendix A"},{"comment":"The IR-safety claim for bH^(3,0) is conditional on the joint QED-QCD factorization formula in Eq. (6), but this theorem is neither proved nor precisely stated. The text says that the proof is 'effectively the same' as that in Ref. [18] and is 'straightforward to verify', and the treatment of the q^2 -> 0 pinch region is only heuristic. Since Eq. (36) is exactly the order-by-order implementation of this theorem, the identification of Eq. (55) as the correct short-distance coefficient requires either a proof sketch for the inclusive e h -> e X case, including the q^2 -> 0 region and the combined QED-QCD gauge links, or an explicit statement of the theorem from Ref. [16] and the conditions under which it applies.","section":"Sec. II, Eq. (6); Sec. IV, Eqs. (35)-(36)"}],"minor_comments":[{"comment":"There are several typographical errors: 'contributiosn' in the first paragraph of Sec. IV, 'infra-safe' in the paragraph after Eq. (55), and 'purbatively' in Sec. V. A careful proofreading pass is needed.","section":"Sec. IV and Sec. V, general"},{"comment":"The logarithm in the perturbative LDF f_{e/e}(xi) is written with (1 - zeta) in the denominator; it should presumably be (1 - xi). Please confirm and correct.","section":"Eq. (61)"},{"comment":"The notation e_l^2 is confusing, because it denotes a sum over all fermion flavors in the photon vacuum polarization, not the electron charge squared. Consider renaming this quantity, for example C_gamma or Sigma_f n_c e_f^2, to avoid confusion with the lepton charge.","section":"Eq. (56)"},{"comment":"The captions and text should state explicitly that the shaded bands in Figs. 15 and 16 are model estimates based on the ad hoc input LDFs and LFFs of Eqs. (63) and (64), not extracted universal functions, and that the horizontal axis region near the dotted vertical line corresponds to small l'_T where the factorized formula is not expected to be reliable.","section":"Figs. 15 and 16"}],"recommendation":"major_revision","confidential_remarks":"The paper is a natural continuation of the authors' program, and the central calculation may well be correct. The two issues that must be resolved before publication are the QED gauge invariance of the operator definitions in Eqs. (19) and (44) and the complete lack of derivation or cross-check for the coefficients in Appendix A. I would suggest asking the authors to provide a supplementary file containing the full calculation or a detailed derivation appendix. The novelty relative to Refs. [15,16] is the explicit NLO coefficient; without a verifiable derivation, the contribution is difficult to assess."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read the Cammarota et al. paper on factorized QED/QCD DIS. Bottom line: the calculation is likely right and worth taking seriously, but there's a technical error in the operator definition that needs to be fixed before I'd trust the formal claim.\n\nWhat's new: Eq. (55) with Appendix A is the first explicit NLO QED hard part for e+q->e+X in the joint factorization scheme. The coefficients are all written out, so the result is reproducible. The QCD part is repackaged known structure-function results, but the conversion procedure (Eqs. 15, 18) is clean and useful. The numerical section is honest: the LDF/LFF inputs are explicitly models, not fits, and they show the resulting uncertainty. I also like that they caught the divergence from using perturbative LDF/LFFs together with NLO hard parts—that's a real subtlety that would trip people up.\n\nSoft spots. The biggest one: the gauge link in Eq. (20) is missing the quark fractional charge e_q. As written, it cancels a charge-1 phase, so the quark PDF defined in Eq. (19) is not QED gauge invariant. Same issue appears in the quark LFF of Eq. (44). This matters because the factorization theorem they rely on requires gauge-invariant operator definitions. The perturbative calculations themselves use e_q correctly, so the hard coefficients are probably fine, but the formal setup is wrong as stated. It's likely a typo—the link should be P exp[-i∫(g_s A_g t^a - e e_q A_γ)]—but it has to be corrected. The paper also does not prove the factorization theorem; it cites Ref. [16] and says the proof is the same as for single-hadron production. That's acceptable as a division of labor, but it means the IR-safety claim is conditional on that theorem.\n\nTwo smaller issues. There's no independent cross-check of the Appendix coefficients—a comparison with the traditional RC result for the pure-electron-radiation term in an appropriate limit would be reassuring. And the numerics use CT18 PDFs that are not the joint QED-QCD PDFs defined in Eq. (19); the authors acknowledge the difference is small, but it's an approximation, not the exact framework.\n\nOverall: the main result is new and the coefficients are explicit enough to be tested. The gauge-link issue is a genuine defect that a referee should catch and the authors should fix. I'd send it to peer review, but I'd want the operator definition corrected and ideally a cross-check added before publication.","headline":"First NLO QED hard coefficient in the joint QED-QCD factorization scheme, likely right but with a technical gauge-link error that needs fixing.","tokens_in":29608,"tokens_out":3475,"would_cite":true,"duration_ms":34926,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V05","81V10","81T18"],"pacs":["12.38.-t","13.60.Hb","12.20.-m"],"model":"deepseek-v4-flash","headline":"The paper computes the next-to-leading-order QED and QCD short-distance hard coefficients for inclusive lepton-hadron deep inelastic scattering in a joint QED-QCD factorization, and claims the NLO QED contribution is completely infrared…","keywords":["Deep inelastic scattering","QED factorization","QCD factorization","Lepton distribution function","Lepton fragmentation function","Photon PDF","Infrared safety","NLO hard coefficients"],"falsifier":"A direct check is to repeat the NLO QED calculation with a different infrared regulator, for example a small photon mass or a photon-energy cutoff instead of dimensional regularization, and verify that the combination in Eq. (36) is finite and independent of the regulator in the limit; any leftover $1/\\epsilon$ pole or cutoff dependence in $\\widehat H^{(3,0)}_{eq\\to eX}$ would show the infrared-safety claim fails. A second, observable check is to measure the inclusive $e+p$ cross section in a kinematic region where the photon-PDF subtraction term is large and compare with the prediction using an independently extracted photon PDF.","tokens_in":28473,"feed_emoji":"⚛️","tokens_out":9184,"duration_ms":93023,"temperature":0.7,"pith_summary":"This paper tries to establish that collision-induced QED radiation in inclusive lepton-hadron deep inelastic scattering can be handled by the same factorization machinery as QCD radiation, rather than by radiative-correction factors with adjustable parameters. It presents the first next-to-leading-order calculation of the short-distance hard coefficients in a joint QED-QCD factorization for the leading electron-quark subprocess, and the central result is that the NLO QED hard part is infrared safe and depends on no free parameter except the standard factorization scale. If the result is right, analyses of DIS and related lepton-hadron observables can include QED radiation without arbitrary photon-phase-space cuts, and the same universal lepton distribution and fragmentation functions can be used across observables. The paper also gives a procedure for converting existing higher-order QCD structure-function results into this framework, and it stresses that the relevant hard scale is the observed lepton transverse momentum rather than $Q^2$.","feed_headline":"QED radiation in deep inelastic scattering tamed at NLO","feed_subtitle":"New calculation handles collision-induced QED and QCD radiation on equal footing, with no ad hoc parameters beyond the factorization scale.","key_machinery":"The load-bearing object is the joint QED-QCD factorization formula in Eq. (6), which writes the DIS cross section as convolutions of universal nonperturbative functions, namely lepton distribution functions (LDFs), lepton fragmentation functions (LFFs), and parton distribution functions (PDFs), with short-distance partonic hard coefficients $\\widehat H_{ia\\to jX}$. The derivation of the NLO QED coefficient uses the subtraction identity in Eq. (36), where the dimensionally regularized partonic cross section $\\sigma^{(3,0)}_{eq\\to eX}$ is freed of all collinear and pinch divergences by subtracting first-order LDF, LFF, and PDF terms, including a photon-PDF term whose convolution absorbs the $q^2\\to 0$ region. The hard coefficient is then organized in powers of quark charge in Eq. (55), and the MS scheme plus DGLAP-type joint evolution equations supply the only scale dependence.","core_discovery":"On the paper's own terms, the discovery is that the NLO QED contribution to the leading partonic channel of inclusive DIS, the hard coefficient $\\widehat H^{(3,0)}_{eq\\to eX}$ given in Eq. (55) with coefficients in Appendix A, is completely infrared safe after collinear-sensitive pieces are subtracted through Eq. (36). The subtraction combines the first-order lepton fragmentation function $D^{(1,0)}_{e/e}$, electron and quark distribution functions $f^{(1,0)}_{e/e}$ and $f^{(1,0)}_{q/q}$, and the photon distribution $f^{(1,0)}_{\\gamma/q}$ convoluted with the lowest-order $e\\gamma\\to eX$ hard part; the last term is what removes the perturbative pinch when the exchanged photon goes on shell and collinear to the quark. Grouping the result by quark charge shows the $e_q^2$ (electron-radiation), $e_q^3$ (two-photon-exchange/interference), and $e_q^4$ (quark-radiation) pieces are each separately finite, and the NLO QCD hard part Eq. (34) is likewise infrared safe with only the factorization-scale dependence. The claim is that no extra parameter is needed beyond the standard factorization scale in the MS scheme.","pith_inferences":["Inference: if the joint factorization theorem holds, the same subtraction pattern should extend to semi-inclusive DIS and TMD observables, where QED radiation changes the orientation of the leptonic plane relative to the hadronic plane; this could remove one model dependence from TMD extractions.","Inference: the claim that LDFs and LFFs are nonperturbative and must vanish at momentum fraction 1 is testable: global fits combining inclusive DIS, SIDIS, and electron-positron annihilation data could extract these functions and check that they agree across observables.","Inference: the photon-PDF subtraction makes a definite prediction at high energies: the photon-initiated $e\\gamma\\to eX$ channel contributes at a level fixed by the evolving photon content of the proton, so tagging low-$Q^2$ or forward-photon events would directly probe the mechanism that removes the $q^2\\to 0$ pinch.","Inference: the emphasis on $\\ell'_T$ rather than $Q^2$ as the hard scale implies that existing global PDF analyses that cut on $Q^2$ alone may include kinematic regions with sizable power corrections; repeating fits with an $\\ell'_T$-based cut could shift extracted distributions at large $y$."],"forward_implications":["NLO QED radiative effects in inclusive DIS can be included without introducing a photon-energy or photon-angle cut, so the predictive uncertainty of radiative corrections is reduced to the standard factorization scale.","The pinch singularity that appears when the exchanged photon is nearly on shell and collinear to the quark is absorbed into the hadron's photon PDF, so the same framework produces the photon-initiated contribution as part of the same calculation.","Existing and future higher-order QCD results for DIS structure functions can be converted into joint-factorization hard coefficients through Eqs. (15) and (18), so the QCD part of the framework is not limited to NLO.","The same LDFs and LFFs, with their nonperturbative endpoint behavior, must be used consistently with the hard coefficients: using perturbative delta-function endpoint distributions together with beyond-LO hard parts produces a spurious divergence, as Eq. (68) shows.","Because the true hard scale is $\\ell'^2_T=(1-y)Q^2$, predictions at large $y$ require a correspondingly larger $Q^2$, otherwise inverse-power corrections to the factorization formula are not small."],"supporting_citations":[{"why":"Classic radiative-correction scheme with free parameter that the paper's parameter-free QED factorization is contrasted with.","marker":"[8]"},{"why":"Earlier factorized treatment of radiative corrections for inelastic lepton-hadron collisions whose azimuthal-modulation discussion motivates treating QED radiation at the amplitude level.","marker":"[15]"},{"why":"Introduces the joint QED-QCD factorization formula (Eq. (6)) and the LDF/LFF framework on which this calculation is built.","marker":"[16]"},{"why":"The QCD factorization proof for inclusive single-hadron production cited as the template for the joint factorization theorem.","marker":"[18]"},{"why":"Operator definitions of parton distribution and decay functions adapted to define LDFs, LFFs, and the photon distribution.","marker":"[19]"},{"why":"Cut-vertex method used to compute the first-order quark/electron/photon distributions and the subtraction terms.","marker":"[20]"},{"why":"PDF set used in the numerical estimates of LDF/LFF evolution and cross-section ratios.","marker":"[24]"}],"fun_headline_variants":["First joint QED+QCD NLO calculation for DIS","QED and QCD radiation factorized equally at NLO","Infrared-safe DIS at NLO with no extra parameters","No ad hoc scales: NLO QED+QCD DIS calculation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the joint QED-QCD factorization formula in Eq. (6) is genuinely valid for inclusive lepton-hadron DIS; the paper assumes this theorem rather than proving it, citing [16] and an analogy to the QCD single-hadron proof [18], so if the combined gauge link or the on-shell-photon pinch region were mishandled, the subtracted hard coefficient would not be the true short-distance quantity.","fun_headline_variants_meta":{"raw":{"variants":["First joint QED+QCD NLO calculation for DIS","QED and QCD radiation factorized equally at NLO","Infrared-safe DIS at NLO with no extra parameters","No ad hoc scales: NLO QED+QCD DIS calculation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1470,"prompt_tokens":963,"completion_tokens":507,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":435}},"tokens_in":579,"tokens_out":507,"duration_ms":5485,"temperature":1.0,"reasoning_tokens":435,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:44:56.916219+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check is to repeat the NLO QED calculation with a different infrared regulator, for example a small photon mass or a photon-energy cutoff instead of dimensional regularization, and verify that the combination in Eq. (36) is finite and independent of the regulator in the limit; any leftover $1/\\epsilon$ pole or cutoff dependence in $\\widehat H^{(3,0)}_{eq\\to eX}$ would show the infrared-safety claim fails. A second, observable check is to measure the inclusive $e+p$ cross section in a kinematic region where the photon-PDF subtraction term is large and compare with the prediction using an independently extracted photon PDF.","supporting_citations":[{"cited_title":"Kripfganz, H","cited_arxiv_id":null,"evidence_quote":"Introduces the joint QED-QCD factorization formula (Eq. (6)) and the LDF/LFF framework on which this calculation is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"PDF set used in the numerical estimates of LDF/LFF evolution and cross-section ratios."}],"review_version":1}