---
title: Dihadron Fragmentation Functions – Overview
url: https://www.emergentmind.com/topics/dihadron-fragmentation-function-dff
type: topic
---

# Dihadron Fragmentation Functions – Overview

Dihadron fragmentation functions (DFFs, often DiFFs) are the two-hadron analog of ordinary fragmentation functions: they describe the correlated production of two identified hadrons from a single fragmenting parton, together with unobserved remnants. In the standard pion-pair channel, they encode the dependence on the pair total momentum, invariant mass, and internal geometry, and they become especially valuable when the hadron-pair invariant mass is small compared with the hard scale, so that the pair is naturally described as originating from the same parton rather than from two different hard partons. Their importance is twofold. First, they probe nonperturbative hadronization with a sensitivity to correlations, resonance production, decay chains, and cascade ordering that is absent in single-hadron fragmentation. Second, in collinear factorization they provide experimentally accessible spin analyzers, most notably the chiral-odd interference fragmentation function \(H_1^{\sphericalangle}\), which couples to quark transversity, as well as helicity-sensitive structures such as \(G_1^\perp\) and \(D_{1LL}\) [1012.0054, 1511.03220, 1712.06384, 2412.00394].

## 1. Kinematics, correlators, and leading functions

For a quark fragmenting into a hadron pair \(h_1 h_2\), the standard pair variables are
\[
P_h = P_1 + P_2, \qquad R = \frac{P_1-P_2}{2}, \qquad M_h^2 = P_h^2.
\]
In the \(\pi^+\pi^-\) channel one also uses
\[
z = \frac{P_h^-}{k^-} = z_1+z_2, \qquad \zeta = \frac{2R^-}{P_h^-} = \frac{z_1-z_2}{z},
\]
together with the polar angle \(\theta\) of one hadron in the pair rest frame and the azimuth \(\phi_R\) of the hadron-pair plane. A useful kinematic relation is
\[
\frac{2|{\bf R}|}{M_h} = \sqrt{1-\frac{4m_\pi^2}{M_h^2}} .
\]
This variable set makes explicit that a DFF retains both the longitudinal momentum sharing of the pair and its internal relative motion [1012.0054, 1206.1836, 1202.0323].

At operator level, the collinear quark–hadron correlator for a pair can be written as
\[
\Delta^q(z,\cos\theta,M_h^2,\phi_R) = \frac{z |\vec R|}{16\,M_h}\int d^2 \vec k_T \; d k^+\,\Delta^q(k;P_h,R) \Big|_{k^- = P_h^-/z},
\]
with
\[
\Delta^q(k,P_h,R)_{ij} = \sum_X \int \frac{d^4\xi}{(2\pi)^4}\, e^{i k\cdot \xi} \langle 0| {\cal U}^{n_+}_{(-\infty,\xi)} \,\psi_i^q(\xi)|P_h, R; X\rangle \langle P_h, R; X| \bar{\psi}_j^q(0)\, {\cal U}^{n_+}_{(0,-\infty)} |0\rangle.
\]
Dirac traces of this correlator define the leading collinear DiFFs, with \(D_1\) the unpolarized DiFF and \(H_1^{\sphericalangle}\) the polarized, chiral-odd interference DiFF [1012.0054].

At leading twist, the quark-quark correlator for two hadrons contains four DiFFs:
\[
D_1, \qquad G_1^\perp, \qquad H_1^{\sphericalangle}, \qquad H_1^\perp .
\]
After integrating over the partonic transverse momentum, the standard collinear limit retains \(D_1\) and \(H_1^{\sphericalangle}\), while \(G_1^\perp\) and \(H_1^\perp\) belong to the transverse-momentum-sensitive formulation [1511.03220, 1707.04999].

| Function | Spin structure | Typical role |
|---|---|---|
| \(D_1\) | Unpolarized quark \(\to\) unpolarized pair | Unpolarized cross sections and normalization |
| \(H_1^{\sphericalangle}\) | Transverse quark spin correlated with pair-plane orientation | Transversity in SIDIS; Artru–Collins asymmetry |
| \(G_1^\perp\) | Quark helicity correlated with transverse-momentum-sensitive handedness of the pair | Weighted \(e^+e^-\) and SIDIS helicity observables |
| \(D_{1LL}\) | Longitudinal helicity correlation between the two hadrons for an unpolarized parent parton | Neighboring dihadron helicity correlation |

For low pair invariant mass, the partial-wave content is central. The standard expansion gives
\[
D_1 \to D_{1,ss+pp}+D_{1,sp}\cos\theta + D_{1,pp}\frac{1}{4}(3\cos^2\theta-1),
\]
\[
\frac{|\bm{R}_T|}{M_h}H_1^{\sphericalangle} \to H_{1,sp}^{\sphericalangle}\sin\theta + H_{1,pp}^{\sphericalangle}\sin\theta\cos\theta.
\]
After averaging over \(\cos\theta\), the dominant surviving pieces are \(D_{1,ss+pp}\) and \(H_{1,sp}^{\sphericalangle}\). Physically, \(H_1^{\sphericalangle}\) is generated by interference between different relative orbital angular momentum states of the hadron pair, especially \(s\)- and \(p\)-wave components [1012.0054, 1511.03220].

## 2. Collinear factorization and the central observables

A defining feature of dihadron production is that, after integration over the transverse momentum of the fragmenting quark, the leading-twist description falls within collinear factorization. The resulting cross sections are products of PDFs and DFFs rather than transverse-momentum convolutions. This is the structural reason DFFs became a preferred route to transversity: the dihadron mechanism survives \(k_T\) integration, unlike the usual single-hadron Collins effect [1012.0054, 1206.1836, 1106.5897].

For SIDIS,
\[
\ell(k) + N(P) \to \ell(k') + H_1(P_1) + H_2(P_2) + X,
\]
the leading-twist cross section can be written as
\[
\frac{d\sigma}{dx \, dy\, dz\, d\phi_S\, d\phi_R\, d M_{h}^2\,d \cos{\theta}}
= \frac{\alpha^2}{x y\, Q^2}
\left\{ A(y)\,F_{UU} + |\mathbf S_T|\,B(y)\,\sin(\phi_R+\phi_S)\,F_{UT} \right\},
\]
with
\[
F_{UU} = x \sum_q e_q^2\, f_1^q(x;Q^2)\, D_1^q(z,\cos\theta,M_h;Q^2),
\]
\[
F_{UT} = \frac{|\mathbf R|\sin\theta}{M_h}\, x \sum_q e_q^2\, h_1^q(x;Q^2)\, H_1^{\sphericalangle\,q}(z,\cos\theta,M_h;Q^2).
\]
Equivalently, after the standard angular integrations, the measured asymmetry takes the form
\[
A_{\rm SIDIS}(x,z,M_h;Q) =
-\frac{B(y)}{A(y)}\,\frac{|\mathbf R|}{M_h}
\frac{\sum_q e_q^2\, h_1^q(x;Q^2)\,H_1^{\sphericalangle\,q}(z,M_h;Q^2)}
{\sum_q e_q^2\, f_1^q(x;Q^2)\,D_1^q(z,M_h;Q^2)} .
\]
The numerator contains the chiral-odd product \(h_1\,H_1^{\sphericalangle}\), while the denominator contains the unpolarized product \(f_1\,D_1\) [1012.0054, 1409.6607].

In \(e^+e^-\) annihilation, the key process is
\[
e^+e^- \to (\pi^+\pi^-)_{\rm jet1}(\pi^+\pi^-)_{\rm jet2}X.
\]
The corresponding Artru–Collins asymmetry is sensitive to the product of two interference DiFFs:
\[
A_{e^+e^-} \propto
\frac{|\mathbf R|}{M_h}\,\frac{|\bar{\mathbf R}|}{\bar M_h}
\frac{\sum_q e_q^2\,H_1^{\sphericalangle\,q}(z,M_h;Q^2)\,H_1^{\sphericalangle\,\bar q}(\bar z,\bar M_h;Q^2)}
{\sum_q e_q^2\,D_1^q(z,M_h;Q^2)\,D_1^{\bar q}(\bar z,\bar M_h;Q^2)} .
\]
This observable calibrates the polarized DiFF in \(e^+e^-\), which can then be used as input in SIDIS to access transversity [1012.0054, 1202.0323, 1409.6607].

The same mechanism extends to hadronic collisions. In \(p p^\uparrow \to (h_1h_2)X\), a \(\sin(\phi_{S_B}-\phi_R)\) modulation appears in collinear kinematics, again proportional to transversity multiplied by \(H_1^{\sphericalangle}\). This supports the standard universality picture for DiFFs across \(e^+e^-\), SIDIS, and hadron-hadron reactions [1511.03220].

## 3. Extraction programs and the transversity route

The modern phenomenology of DiFFs developed as a two-step program. First, \(D_1\) and \(H_1^{\sphericalangle}\) are constrained from \(e^+e^-\to(\pi^+\pi^-)(\pi^+\pi^-)X\). Second, the same functions are inserted into SIDIS on transversely polarized targets to extract the transversity PDF \(h_1^q(x)\). This strategy produced the first extraction of the \(u\)- and \(d\)-flavor transversity distributions in the framework of collinear factorization [1206.1836, 1106.5897].

The practical complication of the early analyses was that the unpolarized dihadron cross section had not been measured directly. Consequently, the unpolarized DiFF \(D_1\) was first parametrized to reproduce pion-pair yields from PYTHIA tuned to Belle kinematics, while \(H_1^{\sphericalangle}\) was fitted to Belle’s measured Artru–Collins asymmetry [1202.0323, 1206.1836, 1409.6607]. Within this framework, the integrated moments
\[
n_q(Q^2)=\int dz\, dM_h\, D_1^q(z,M_h;Q^2),
\qquad
n_q^\uparrow(Q^2)=\int dz\, dM_h\, \frac{|\mathbf R|}{M_h}\, H_1^{\sphericalangle\,q}(z,M_h;Q^2)
\]
enter the SIDIS asymmetry in a particularly transparent way:
\[
A_{\text{SIDIS}(x,Q^2)}
=
\frac{\sum_q e_q^2\, h_1^q(x,Q^2)\, n_q^\uparrow(Q^2)}
{\sum_q e_q^2\, f_1^q(x,Q^2)\, n_q(Q^2)} .
\]
Under the usual isospin and charge-conjugation assumptions for \(\pi^+\pi^-\), proton and deuteron data provide two independent flavor combinations and permit the separation of \(h_1^{u_v}\) and \(h_1^{d_v}\) [1206.1836, 1409.6607].

A decisive methodological update was the Belle reanalysis with the replica method. In that study, the Belle asymmetry points were replicated 100 times, each replica was independently fitted, and the uncertainty bands were defined by the central 68% of the replica ensemble. The point of this procedure was that the fit tends to push the function toward its bounds, making Gaussian or Hessian assumptions unreliable. The resulting transversity extraction, combined with more precise COMPASS proton data, produced a more realistic estimate of the uncertainties and suggested that the valence up component was smaller and had a narrower error band than in the earlier extraction, while the down component was largely unchanged because the deuteron input was unchanged [1409.6607].

The phenomenological impact of this program extends beyond the extraction of \(h_1\). In the flexible scenario of the 2015 analysis, the isovector tensor charge was quoted as
\[
g_T = 0.81 \pm 0.44 \qquad \text{at } Q^2=4~\text{GeV}^2,
\]
and the same DiFF-based transversity determination was used to infer the bounds
\[
|\epsilon_T| < 0.00162,
\qquad
|\epsilon_T| < 0.00139 \quad \text{at } 90\%\ \text{C.L.}
\]
for two different uncertainty treatments. This established DiFF phenomenology as part of the interface between hadron structure and low-energy searches for physics beyond the Standard Model [1608.04429].

## 4. Hadronization dynamics, cascade models, and resonance effects

DFFs are more sensitive than single-hadron FFs to the internal structure of the hadronization cascade. In the NJL-jet model, the quark emits hadrons sequentially and emitted hadrons do not reinteract. The corresponding dihadron fragmentation function \(D_q^{h_1h_2}(z_1,z_2)\) satisfies a recursive equation whose first two terms are “driving terms” and whose third term accounts for the case in which both observed hadrons are produced after the first step of the cascade:
\[
D^{h_1,h_2}_q(z_1,z_2)
=
\hat d_q^{h_1}(z_1)\,
\frac{D_{q_1}^{h_2}\!\left(\frac{z_2}{1-z_1}\right)}{1-z_1}
+
\hat d_q^{h_2}(z_2)\,
\frac{D_{q_2}^{h_1}\!\left(\frac{z_1}{1-z_2}\right)}{1-z_2}
+
\sum_Q \int_{z_1+z_2}^{1}\frac{d\eta}{\eta^2}\,
\hat d_q^Q(\eta)\,
D_Q^{h_1,h_2}\!\left(\frac{z_1}{\eta},\frac{z_2}{\eta}\right).
\]
The physical interpretation is direct: at large \(z_1\) or \(z_2\), the driving terms dominate because one of the hadrons is likely to have been produced in the first emission step; at smaller \(z\), the higher-order cascade term becomes more significant. In some channels, particularly for strange-quark initiated fragmentation, the driving term can vanish and the entire DFF is generated recursively [1201.1351, 1202.4036].

Monte Carlo implementations exploit a number-density interpretation:
\[
D_q^{h_1 h_2}(z, M_h^2)\,\Delta z\,\Delta M_h^2
=
\left\langle
N_q^{h_1 h_2}(z, z+\Delta z;\, M_h^2, M_h^2+\Delta M_h^2)
\right\rangle .
\]
Operationally, one bins all hadron pairs produced in a large sample of simulated fragmentation chains and extracts the DFF from the average pair multiplicity in each \((z,M_h^2)\) bin. This procedure can be applied to primary hadrons only, or to the full final state after strong resonance decays [1307.8125, 1310.1917].

A major result of the hadronization studies is that vector meson decays drastically affect pion DFFs. In the benchmark channel \(u\to\pi^+\pi^-\), including pions from vector meson decays makes the DFF much larger, typically up to an order of magnitude larger, than the primary-pion result. The enhancement is not merely a resonance correction but a combinatorial effect: once a vector meson decays, each decay product can pair with other hadrons in the event. In the restricted \(\pi\)- and \(\rho\)-only example discussed in the literature, the listed secondary channels generate a total of 17 possible \(\pi^+\pi^-\) pairs, compared with only the single direct primary channel. The invariant-mass dependence correspondingly acquires a pronounced \(\rho^0\) peak near \(M_h^2\approx(0.77~\mathrm{GeV})^2\) and a low-\(M_h^2\) enhancement from \(\omega\to\pi^+\pi^-\pi^0\) when the \(\pi^0\) is unobserved. The fact that the same qualitative behavior appears in both the NJL-jet model and PYTHIA 8.1 with Lund string fragmentation shows that this is not a model artifact [1307.8125, 1310.1917].

These studies also show that naive factorization assumptions in the \((z,M_h)\) dependence are generally unsafe. In the Belle Monte Carlo analysis of the unpolarized cross section, the invariant-mass shape changes across \(z\)-bins, and the continuum decreases with \(z\) differently for light and charm flavors. That finding motivated genuinely two-dimensional parameterizations of \(D_1\) rather than products of separate \(z\)- and \(M_h\)-dependent functions [1012.0054].

## 5. Evolution equations, perturbative regimes, and helicity-sensitive extensions

At low scales, model DFFs must be evolved to experimental scales. In the early NJL-jet program, the model input was defined at \(Q_0^2 = 0.2~\mathrm{GeV}^2\), and leading-order evolution to \(Q^2=4~\mathrm{GeV}^2\) or \(100~\mathrm{GeV}^2\) was found to shift strength toward lower \(z\), suppressing the high-\(z\) region while broadly preserving the basic \(M_h^2\) structure. For DFFs integrated over invariant mass, the LO evolution equations contain both a homogeneous term and an inhomogeneous term built from single-hadron FFs, reflecting the possibility that the two observed hadrons are produced after a perturbative splitting into two daughter partons [1207.4267, 1202.4036, 1310.1917].

This inhomogeneous structure became conceptually explicit in the 2024 analysis of neighboring dihadron helicity correlations. For the unpolarized DiFF,
\[
\frac{d D_{1,i}^{h_1 h_2}}{d\ln \mu_f^2}
=
\text{homogeneous DGLAP term}
+
\text{source term from } D_1^{h_1} D_1^{h_2},
\]
whereas the correlated helicity DiFF \(D_{1LL}^{h_1h_2}\) obeys an analogous equation in which the source term involves the longitudinal spin-transfer FF \(G_{1L}\):
\[
\frac{d D_{1LL,i}^{h_1 h_2}}{d\ln \mu_f^2}
=
\text{homogeneous DGLAP term}
+
\text{source term from } G_{1L}^{h_1} G_{1L}^{h_2}.
\]
This means that even if a DiFF is set to zero at a low scale, perturbative splittings regenerate it at higher scales. In the specific \(\Lambda\bar\Lambda\) study, the helicity-correlation observable
\[
{\cal C}_{LL,i}^{\Lambda\bar\Lambda}
=
\frac{D_{1LL,i}^{\Lambda\bar\Lambda}}{D_{1,i}^{\Lambda\bar\Lambda}}
\]
was shown to be sensitive to the flavor structure of \(G_{1L}\), especially for strange quarks and gluons, thereby opening a new way to probe circularly polarized gluon hadronization in unpolarized collisions [2412.00394].

A complementary helicity-sensitive DiFF is \(G_1^\perp\), which correlates the longitudinal polarization of a fragmenting quark with the transverse-momentum structure of the hadron pair. A crucial clarification of the literature is that the original BELLE search found no signal because the previously proposed unweighted asymmetry actually vanishes. The corrected observables are weighted asymmetries: in \(e^+e^-\), a \(q_T^2\)-weighted sine/cosine combination of azimuths isolates products of Fourier moments of \(G_1^\perp\); in SIDIS, a \(P_{h\perp}\sin(\phi_{qR})/M_h\) weight yields a collinear asymmetry proportional to \(g_1(x)\,G_1^\perp(z,M_h^2)\). Quark-jet Monte Carlo studies generate a nonvanishing \(G_1^\perp\), but also predict that the effect is small and becomes relatively suppressed as the number of hadron emissions increases [1712.06384, 1707.04999].

In the large-invariant-mass regime, DFFs admit a perturbative treatment. For \(M_{hh}\gg\Lambda_{\rm QCD}\), the hadron pair is described as a perturbative splitting followed by ordinary single-hadron fragmentation. In this regime,
\[
D_1 \sim \frac{1}{M_{hh}^2},
\qquad
H_1^{\sphericalangle} \sim \frac{1}{M_{hh}^3},
\]
so the corresponding SIDIS transverse-spin asymmetry behaves like \(1/M_{hh}\). The same analysis showed that \(H_1^{\sphericalangle}\) and the Collins fragmentation function are driven by the same underlying twist-3 collinear fragmentation correlators, and that in the intermediate region \(\Lambda_{\rm QCD}\ll M_{hh}\ll Q\), collinear factorization in terms of DiFFs and collinear factorization in terms of single-hadron FFs give the same result [1101.3273].

## 6. Interpretation, conceptual debates, and recent directions

A recurrent conceptual issue is how literally one should interpret DFFs as number densities. The current position of the field is qualified rather than absolute. DFFs do support a legitimate number-density interpretation in the same sense as ordinary fragmentation functions: as renormalized, operator-based densities for identified hadrons in restricted regions of phase space, tied to appropriate factorization theorems. What is not settled is the stronger claim that a literal global multiplicity sum rule,
\[
\sum_{h_1,h_2}\int dz_2\int dz_1\int d^2\mathbf P_{1\perp}\int d^2\mathbf P_{2\perp}\, D^{h_1h_2/i}
=
\langle N(N-1)\rangle,
\]
should be taken as a defining first-principles constraint. The objection is that the naive derivation fails because the fragmenting quark state and the purely hadronic asymptotic states are orthogonal in the usual completeness relation. The practical conclusion is that factorization, not a disputed multiplicity sum rule, should be the guiding criterion for acceptable DFF definitions [2404.02281].

This debate matters because slightly different operator definitions and phase-space parametrizations have been used across the DiFF literature. The present view is that such differences do not automatically invalidate a number-density interpretation, provided the corresponding hard parts and evolution kernels are adjusted consistently. A plausible implication is that the DiFF sector should be treated much like the rest of QCD factorization: renormalized distributions are “probability-like” or quasiprobability objects, but not literal probabilities once ultraviolet subtraction and scheme dependence are taken seriously [2404.02281].

The scope of DiFF phenomenology has also broadened substantially. At subleading twist, the same \(H_1^{\sphericalangle}\) that enables transversity extraction provides the cleanest access to the poorly known twist-3 PDF \(e(x)\) through beam-spin asymmetries in dihadron SIDIS, while related observables also involve \(h_L(x)\). In the 2015 review framework, this placed DiFFs at the center of a program linking hadron-pair production to dynamical chiral symmetry breaking, the nucleon scalar structure, and higher-twist quark-gluon dynamics [1511.03220, 1608.04429].

Two recent developments extend DiFFs into new territories. First, a 2025 extraction of unpolarized \(\pi^+\pi^-\) DFFs used Belle measurements of the differential unpolarized cross section rather than relying only on Monte Carlo yields. That analysis fitted 344 data points after cuts, employed NNLO coefficient functions implemented in APFEL++, compared a 71-parameter physics-informed fit with a 205-parameter neural-network parametrization, and found that the neural-network fit achieved a lower \(\chi^2/N_{\text{data}}\) but also exhibited larger uncertainty bands, especially for the essentially unconstrained gluon DFF [2509.11855]. Second, the 2025 EEC study showed that in the collinear limit of the energy-energy correlator, the relevant jet function can be written directly in terms of single-hadron FFs and transverse-momentum-sensitive DFFs. In that formulation, DFFs become the universal nonperturbative functions governing the confinement transition region \(Q^2 z \sim \Lambda_{\rm QCD}^2\), thereby linking hadronization in \(e^+e^-\to h_1h_2X\) to the nonperturbative structure of energy correlators [2507.11495].

Taken together, these developments define the present status of DFFs. They are simultaneously hadronization observables, spin analyzers, inputs to collinear extractions of transversity and higher-twist PDFs, test objects for universality across \(e^+e^-\), SIDIS, and hadronic collisions, and increasingly precise nonperturbative functions in their own right. The central lesson of the field is not merely that two hadrons carry more information than one, but that the pair invariant mass, relative momentum, partial-wave content, and spin correlations provide a uniquely efficient organization of that information in QCD [1012.0054, 1511.03220, 2509.11855].

Source: https://www.emergentmind.com/topics/dihadron-fragmentation-function-dff