Numerical Study of MRW-Type Unintegrated Double Parton Distribution Functions from Non-Factorized DPDFs

K. Azizi; R. Kord Valeshabadi; S. Rezaie

Numerical Study of MRW-Type Unintegrated Double Parton Distribution Functions from Non-Factorized DPDFs

T0 review · 2 major / 2 minor · reviewed 2026-05-20 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read MRW-type prescriptions allow direct construction of unintegrated double parton distributions from non-factorized collinear DPDFs without piecewise splitting.

desk verdict This paper numerically extends MRW prescriptions to unintegrated double parton distributions from non-factorized GS09 inputs with a double DGLAP check, but offers limited external validation. read the letter →

arxiv 2605.15853 v2 pith:ULEMJ5KE submitted 2026-05-15 hep-ph hep-exhep-th

R. Kord Valeshabadi , S. Rezaie , K. Azizi This is my paper · ORCID

classification hep-phhep-exhep-th

keywords doublepartondistributionfunctionsunintegrateddistributionsMRWprescriptionnon-factorizedDPDFsDGLAPevolutionscatteringtransversemomentumdependence

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper numerically constructs unintegrated double parton distribution functions by adding transverse momentum dependence to non-factorized collinear DPDFs using MRW-inspired models. It employs the GS09 parametrization evolved via a double DGLAP framework validated by the momentum sum rule. The study examines three variants: double modified KMRW, double virtuality-ordered MRW, and a normalization-matched version of the latter. These can be applied to the entire DPDF including non-homogeneous terms, unlike conventional methods that require splitting into pieces. The results show preservation of normalization in one model, adjustments in others, and strong dependence on parton channels for the effects of non-factorized correlations.

What carries the argument

MRW-type prescriptions (DMKMRW, DVO-MRW, and matched DVO-MRW) for generating transverse-momentum dependence from non-factorized collinear DPDFs

What would settle it

A measurement of double parton scattering observables at the LHC showing significant normalization mismatch with the DVO-MRW predictions even after matching, or a failure of the evolved distributions to satisfy the momentum sum rule.

Watch

Extended reading notes

Core claim

The authors demonstrate that MRW-inspired prescriptions for unintegrated double parton distributions can be applied directly to the full collinear double parton distribution functions from the GS09 parametrization, including their non-homogeneous components. This avoids the piecewise treatment needed in the standard leading-order MRW construction. Through numerical evolution to unequal scales and analysis of normalization, transverse momentum dependence, and flavor sensitivity, they find that the double modified KMRW model preserves normalization while the virtuality-ordered model requires matching to do so, with non-factorized effects manifesting differently across models and being most 1

Load-bearing premise

The GS09 parametrization provides an accurate enough description of non-factorized longitudinal correlations in double parton distributions to be directly extended to transverse momentum dependence using the MRW prescriptions.

Editorial extensions

If this is right

The DMKMRW prescription preserves the normalization of the underlying DPDF by construction.
The DVO-MRW model exhibits nontrivial normalization deviations that are eliminated in its normalization-matched version.
Non-factorized longitudinal correlations affect the transverse momentum shape in virtuality-ordered models but enter primarily through the longitudinal part in DMKMRW.
Largest deviations from factorized expectations appear near the double parton kinematic boundary and in channels with valence-number or quark-antiquark correlations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Such direct constructions could simplify calculations of double parton scattering cross sections at colliders by incorporating correlations more seamlessly.
The channel-dependent sensitivity highlights the need to model specific parton flavor combinations carefully in multi-parton interaction studies.
Extending this numerical approach to other parametrizations or evolution orders might reveal additional impacts on transverse momentum distributions.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Request a human review

A listed scientist reviews the paper for a fee and the review publishes here regardless of verdict. See the reviewers or get listed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

read the letter

This paper numerically constructs unintegrated double parton distribution functions from the non-factorized GS09 collinear DPDFs by extending MRW-type prescriptions to the double-parton case. They use a double DGLAP evolution code validated by the momentum sum rule, then generate transverse momentum dependence with DMKMRW, DVO-MRW, and a matched version. What is new is the direct application to the complete DPDF without the piecewise treatment of conventional LO-MRW, plus the comparative analysis of how non-factorized longitudinal correlations show up in the unintegrated versions. The work does a reasonable job laying out the normalization behavior, flavor dependence, and where the largest effects appear, especially near the double-parton kinematic boundary. The soft spots are the lack of any comparison to experimental data or other unintegrated models, and no quantitative uncertainties or stability tests for the numerics. The assumption that the same transverse-momentum generating factors apply uniformly to the non-homogeneous component may need more scrutiny, as the stress-test suggests it could affect the kT spectrum without additional scale adjustments; the paper reports deviations but does not isolate this effect separately. That said, the sum-rule check is a positive point and the results are presented transparently. This is aimed at QCD phenomenologists modeling double parton scattering for LHC analyses. A reader looking for concrete implementations or sensitivity studies in this niche would find it useful. It deserves a serious referee because the numerical framework is solid enough to build on, even if more validation would help. I recommend sending it to peer review.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a numerical study of unintegrated double parton distribution functions (uDPDFs) constructed from the non-factorized GS09 collinear DPDFs. The collinear DPDFs are evolved to unequal scales using a numerical double DGLAP framework validated by the momentum sum rule. Three MRW-inspired prescriptions (DMKMRW, DVO-MRW, and a normalization-matched DVO-MRW variant) are applied to generate transverse-momentum dependence. The central claim is that these prescriptions can be applied directly to the full collinear DPDF, including its non-homogeneous component, thereby avoiding the piecewise treatment required in conventional LO-MRW constructions. The study examines normalization, kT spectra, flavor dependence, and sensitivity to longitudinal correlations, reporting channel-dependent non-factorized effects that are largest near the double-parton kinematic boundary.

Significance. If the direct application of the MRW factors to the non-homogeneous term is justified, the work supplies useful numerical benchmarks for modeling double parton scattering, particularly the interplay between longitudinal correlations and transverse-momentum distributions. The momentum-sum-rule validation of the double DGLAP evolution is a positive technical check, and the explicit comparison of normalization-preserving versus virtuality-ordered prescriptions quantifies model differences in a reproducible way.

major comments (2)

[momentum sum rule validation and numerical results on non-homogeneous component] The momentum-sum-rule validation is performed only at the collinear level and does not isolate whether the unintegrated non-homogeneous component recovers the input GS09 distribution after integration over transverse momentum. This check is required to confirm that the MRW-style Sudakov or virtuality-ordering factor acts consistently on the inhomogeneous term, whose evolution kernel differs from the homogeneous one.
[analysis of transverse momentum dependence and non-factorized effects] The claim that the prescriptions can be applied directly to the full DPDF (including non-homogeneous part) without additional scale-dependent adjustment rests on the assumption that the transverse-momentum generating factor is uniform regardless of generation scale. The largest reported deviations occur precisely near the double-parton kinematic boundary, where this assumption is most likely to distort the kT spectrum; no explicit test isolating the non-homogeneous contribution is presented.

minor comments (2)

[abstract and numerical results] The abstract states that the DMKMRW model is normalization preserving by construction while DVO-MRW shows deviations removed by matching, but no quantitative error estimates or cutoff sensitivity studies are provided for the reported normalization deviations.
[prescription definitions] Notation for the three prescriptions (DMKMRW, DVO-MRW, matched version) should be defined once in a dedicated subsection with explicit formulas for the Sudakov factors and matching procedure to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and indicate the revisions we plan to make.

read point-by-point responses

Referee: The momentum-sum-rule validation is performed only at the collinear level and does not isolate whether the unintegrated non-homogeneous component recovers the input GS09 distribution after integration over transverse momentum. This check is required to confirm that the MRW-style Sudakov or virtuality-ordering factor acts consistently on the inhomogeneous term, whose evolution kernel differs from the homogeneous one.

Authors: We agree with the referee that validating the integration of the unintegrated distributions back to the collinear level, specifically separating the non-homogeneous component, would strengthen the presentation. In our work, the double DGLAP evolution framework is validated using the momentum sum rule at the collinear level for the full DPDF. The MRW prescriptions are then applied to this evolved full DPDF. Since the MRW construction is such that the transverse momentum integration is designed to recover the collinear input by construction (via the Sudakov or ordering factors), we expect consistency. However, to explicitly address the concern for the inhomogeneous term, we will include in the revised manuscript a numerical demonstration that the kT-integrated uDPDFs reproduce the input GS09 DPDFs for both the homogeneous and non-homogeneous contributions. This will be added as an additional validation step. revision: yes
Referee: The claim that the prescriptions can be applied directly to the full DPDF (including non-homogeneous part) without additional scale-dependent adjustment rests on the assumption that the transverse-momentum generating factor is uniform regardless of generation scale. The largest reported deviations occur precisely near the double-parton kinematic boundary, where this assumption is most likely to distort the kT spectrum; no explicit test isolating the non-homogeneous contribution is presented.

Authors: The application of the MRW-type factors is performed on the collinear DPDF after it has been evolved to the final scales using the double DGLAP evolution. This evolution already accounts for the different kernels acting on the homogeneous and non-homogeneous terms up to the final scales. The transverse momentum dependence is then introduced via the MRW-inspired models at these final scales, consistent with how such models are applied in the single-parton case. We do not introduce additional scale-dependent adjustments because the prescriptions are scale-independent in their application once the collinear input is fixed. Regarding the deviations near the kinematic boundary, these arise from the longitudinal correlations in the GS09 DPDFs themselves. While we have not presented an explicit isolation of the non-homogeneous contribution to the kT spectra in the current version, the results for the full DPDF already incorporate all effects. In the revision, we will add a short discussion clarifying this point and, if space permits, a supplementary figure showing the kT spectra with and without the non-homogeneous term to quantify its impact. revision: partial

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: external GS09 input plus published MRW prescriptions produce numerical outputs

full rationale

The derivation chain begins with the external GS09 collinear DPDF parametrization as input, followed by numerical double DGLAP evolution validated against the momentum sum rule at the collinear level. Transverse-momentum dependence is then generated by applying published MRW-inspired prescriptions (DMKMRW, DVO-MRW and matched variant) directly to the full evolved DPDF, including its non-homogeneous term. No quantity is defined in terms of a fitted parameter that is subsequently relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled through prior author work. The central claim that the prescriptions avoid piecewise treatment is a direct numerical statement rather than an equivalence that reduces to the input by construction. The paper therefore remains self-contained against external benchmarks.

Assumptions & free parameters 0 free parameters · 2 assumptions · 0 invented entities

The study rests on the validity of the GS09 parametrization and the applicability of single-parton MRW constructions to the double-parton case; no new free parameters are introduced beyond those already present in the input DPDFs.

assumptions (2)

domain assumption The GS09 DPDFs provide a reliable description of non-factorized longitudinal correlations inside the proton.
Used directly as input for the evolution and MRW constructions.
standard math The double DGLAP evolution equations can be solved numerically to unequal scales while preserving the momentum sum rule.
Invoked to justify the scale evolution step.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Numerical Study of MRW-Type Unintegrated Double Parton Distribution Functions from Non-Factorized DPDFs." pith.science (2026). https://pith.science/paper/ULEMJ5KE

@misc{pith2026260515853,
  author       = {Pith},
  title        = {Pith review of: Numerical Study of MRW-Type Unintegrated Double Parton Distribution Functions from Non-Factorized DPDFs},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ULEMJ5KE}},
  note         = {Machine review of arXiv:2605.15853}
}

read the original abstract

Double parton scattering provides a sensitive probe of multiparton correlations inside hadrons. In this work we present a numerical study of unintegrated double parton distribution functions constructed from non-factorized collinear DPDFs. As input we use the \texttt{GS09} DPDFs and evolve them to unequal scales with a numerical double DGLAP evolution framework, which is validated through the corresponding momentum and valence-number sum rules. We investigate MRW-inspired prescriptions for generating transverse momentum dependence in the double parton case. In particular, we study the double LO-MRW approach (DLO-MRW) in the perturbative region, the double modified KMRW approach (DMKMRW), the double virtuality ordered MRW (DVO-MRW) model, and a normalization-matched version of the latter. The DMKMRW, DVO-MRW, and matched DVO-MRW prescriptions can be applied directly to the full collinear DPDF, including its nonhomogeneous component, and avoid the piecewise treatment required in the conventional LO-MRW construction. We analyze the normalization, transverse momentum dependence, flavor dependence, and sensitivity to longitudinal DPDF correlations of the resulting distributions. The DMKMRW model is normalization preserving by construction, while the DVO-MRW model shows nontrivial normalization deviations that are removed in the matched version. Non-factorized effects are strongly channel-dependent. In DMKMRW, they enter mainly through longitudinal DPDF correlations, whereas in the LO-MRW and virtuality-ordered models, they also modify the transverse momentum shape. The largest deviations occur near the double parton kinematic boundary and in channels affected by valence-number and quark-antiquark correlations.