---
title: Transverse-Momentum Dihadron Fragmentation Functions
url: https://www.emergentmind.com/topics/transverse-momentum-sensitive-dihadron-fragmentation-functions
type: topic
---

# Transverse-Momentum Dihadron Fragmentation Functions

Searching arXiv for the cited and closely related papers on transverse-momentum-sensitive dihadron fragmentation functions.
arXiv search query: "transverse momentum sensitive dihadron fragmentation functions interference fragmentation helicity-dependent G1 perp"
Transverse-momentum-sensitive dihadron fragmentation functions are fragmentation correlators for a parton that hadronizes into two observed hadrons plus unobserved remnants, with explicit sensitivity to the pair’s invariant mass, relative momentum, and transverse-momentum geometry. In the standard notation of dihadron fragmentation, the observed pair carries total momentum $P_h=P_1+P_2$, relative momentum $R=(P_1-P_2)/2$, invariant mass $M_h$, and azimuth $\phi_R$; these variables make the two-hadron system more informative than a single hadron because the orientation of the hadron-pair plane can correlate with the spin of the fragmenting quark. In the collinear limit, after integrating over the parent parton transverse momentum, the formalism reduces to dihadron fragmentation functions (DiFFs) such as the unpolarized $D_1$ and the chiral-odd interference fragmentation function $H_1$, while in the transverse-momentum-dependent setting the fragmentation correlator retains additional angular structure and Fourier moments that encode helicity- and transversity-sensitive observables [1012.0054].

## 1. Kinematics and correlator structure

The basic collinear dihadron correlator discussed for $\pi^+\pi^-$ production is
$$
\Delta^q(z,\cos\theta,M_h^2,\phi_R)=\frac{z|\vec R|}{16M_h}\int d^2\vec k_T\, dk^+\, \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^{ik\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 .
$$
Here $z=P_h^-/k^-$ is the light-cone momentum fraction carried by the hadron pair. Once the transverse momentum of the parent parton is integrated out, the framework is collinear factorization rather than transverse-momentum-dependent factorization, and the cross section becomes a product of PDFs and DiFFs rather than a convolution in transverse momentum [1012.0054].

At leading twist, the collinear projection leaves two functions in the discussion,
$$
D_1^q(z,\cos\theta,M_h^2)=4\pi\,\mathrm{Tr}\!\left[\Delta^q(z,\cos\theta,M_h^2,\phi_R)\gamma^-\right],
$$
and
$$
\epsilon_T^{ij} R_{Tj} M_h\, H_1^q(z,\cos\theta,M_h^2) =4\pi\,\mathrm{Tr}\!\left[\Delta^q(z,\cos\theta,M_h^2,\phi_R)\, i\sigma^{i-}\gamma_5\right].
$$
In the paper’s language, $D_1$ is the unpolarized DiFF and $H_1$ is the chiral-odd, transverse-spin-dependent DiFF, also called the interference fragmentation function because it arises from interference of different partial waves, especially $s$- and $p$-waves [1012.0054].

The explicitly transverse-momentum-sensitive formulation keeps the pair transverse momenta before integration. A quark dihadron fragmentation function can then be defined as
$$
D_{q\to h_j h_k}(y_j,y_k,\vec P_{j\perp},\vec P_{k\perp}) = \int \frac{dr^+\, d^2\vec r_\perp}{8\,y_j y_k (2\pi)^5} \,e^{ik^- r^+/2} \sum_X \operatorname{Tr}\!\left[ \gamma^-\, \langle 0|\psi_q(r^+,0,\vec r_\perp)|h_j,h_k;X\rangle \langle h_j,h_k;X|\bar\psi_q(0)|0\rangle \right],
$$
with analogous gluon definitions. In this form the nonperturbative object depends on the hadron longitudinal momentum fractions and on the hadron transverse momenta in the parton-collinear frame [2507.11495].

## 2. Spin-dependent content: transversity, helicity, and azimuthal structure

The central spin-dependent collinear DiFF is the interference function $H_1^{\sphericalangle}$, which converts the transverse polarization of the fragmenting quark into an azimuthal asymmetry of the hadron-pair plane. In semi-inclusive deep-inelastic scattering on a transversely polarized target,
$$
e(l)+N^\uparrow(P)\to e(l')+\pi^+(P_1)+\pi^-(P_2)+X,
$$
the relevant single-spin asymmetry is
$$
A_{UT}^{\sin(\phi_R+\phi_S)\sin\theta}(x,y,z,M_h^2) \propto -\frac{|\bm R|}{M_h}\, \frac{\sum_q e_q^2\, h_1^q(x)\, H_{1,q}^{sp}(z,M_h^2)} {\sum_q e_q^2\, f_1^q(x)\, D_{1,q}^{ss+pp}(z,M_h^2)}.
$$
This is the key transversity relation: the numerator contains the transversity distribution $h_1^q(x)$ multiplied by the chiral-odd interference DiFF $H_1^{sp}$, and the denominator contains the unpolarized PDF $f_1^q(x)$ and the unpolarized DiFF $D_1^{ss+pp}$ [1012.0054].

In the TMD dihadron formalism, the fragmentation sector is richer because the two-hadron final state can carry its own polarization and partial-wave content. The TMD dihadron cross section is significantly more complex than that for single-hadron SIDIS, and the final state can be analyzed in partial waves $(\ell,m)$. The paper explicitly states that the collinear interference fragmentation function $H_1^{\varangle}$ is associated with the $(\ell=1,m=1)$ partial wave, that the Collins-like TMD moments with $(2,\pm 2)$ correspond to specific partial-wave angular modulations, and that $H_1^{\varangle}$ receives both $sp$ and $pp$ interference contributions [1308.2989].

A distinct TMD DiFF is the helicity-dependent function $G_1^\perp$, which encodes a correlation between the longitudinal polarization of a fragmenting quark and the azimuthal/transverse-momentum structure of a detected hadron pair. The relevant leading-twist part of the $e^+e^- \to (h_1h_2)(\bar h_1\bar h_2)X$ cross section contains
$$
\mathcal{F}\!\Big[ D_1^a \bar{D}_1^{\bar{a}}\Big] - \mathcal{F}\!\Big[\frac{(R\times k)_3}{M_h^2} \frac{(\bar R\times \bar k)_3}{\bar M_h^2} G_1^{\perp a}\bar G_1^{\perp \bar a}\Big],
$$
where the convolution
$$
\mathcal{F}[w D^a \bar D^{\bar a}] = \int d^2k\, d^2\bar k\; \delta^2(k_T+\bar k_T-q_T)\, w(k,\bar k,R,\bar R)\, D^a\,\bar D^{\bar a}
$$
shows that the azimuthal structure is inseparable from the transverse-momentum convolution unless an appropriate weighting is introduced [1712.06384].

## 3. Fourier moments, weighted asymmetries, and the resolution of null searches

To expose the angular structure of TMD DiFFs, the unintegrated functions are expanded in Fourier cosine series in the relative angle between $k_T$ and $R_T$,
$$
D^a(z,\xi,k^2,R^2,\cos\theta_{kR}) = \frac{1}{\pi}\sum_{n=0}^\infty \frac{\cos(n\theta_{kR})}{1+\delta^{0,n} D^{a,[n]}(z,\xi,k^2,R^2),
$$
with a corresponding decomposition for $G_1^\perp$. The integrated unpolarized DiFF is then
$$
D_1^{a}(z,M_h^2) = z^2 \int d^2k \int d\xi\; D_1^{a,[0]}(z,\xi,k^2,R^2).
$$
These Fourier moments are the objects selected by weighted observables [1712.06384].

A central result is that the earlier BELLE-style unweighted asymmetry used to search for $G_1^\perp$ should vanish when the correct cross section is used:
$$
\langle \cos(2(\phi_R-\phi_{\bar R})) \rangle = 0.
$$
More generally, any unweighted moment depending only on the hadron-pair relative azimuths fails to isolate the required structure. The explanation is that the old angular weight did not properly break the transverse-momentum convolution between the quark and antiquark fragmentation functions [1712.06384].

The new proposal is to weight the $e^+e^-$ cross section by the virtual photon transverse momentum squared, $q_T^2$, multiplied by a trigonometric combination of azimuthal angles. For the simplest case,
$$
\left\langle \frac{q_T^2\left(3\sin(\phi_{qR})\sin(\phi_{q\bar R})+\cos(\phi_{qR})\cos(\phi_{q\bar R})\right)} {M_h\bar M_h} \right\rangle
$$
becomes
$$
\frac{12\alpha^2 A(y)}{\pi Q^2} \sum_{a,\bar a} e_a^2 \Big(G_1^{\perp a,[0]}-G_1^{\perp a,[2]}\Big) \Big(\bar G_1^{\perp\bar a,[0]}-\bar G_1^{\perp\bar a,[2]}\Big).
$$
The paper defines the integrated helicity-dependent DiFF as
$$
G_1^{\perp a}(z,M_h^2) \equiv G_1^{\perp a,[0]}(z,M_h^2)-G_1^{\perp a,[2]}(z,M_h^2),
$$
with
$$
G_1^{\perp a,[n]}(z,M_h^2) \equiv z^2\int d^2k\int d\xi\; \left(\frac{k^2}{2M_h^2}\right)\frac{|R|}{M_h} \,G_1^{\perp a,[n]}(z,\xi,k^2,R^2).
$$
The corresponding asymmetry,
$$
A_{e^+e^-}^{\Rightarrow}(z,\bar z,M_h^2,\bar M_h^2) = 4\, \frac{ \sum_{a,\bar a} e_a^2\, G_1^{\perp a}(z,M_h^2)\, G_1^{\perp\bar a}(\bar z,\bar M_h^2) }{ \sum_{a,\bar a} e_a^2\, D_1^{a}(z,M_h^2)\, D_1^{\bar a}(\bar z,\bar M_h^2) },
$$
is a collinear expression obtained after the $q_T^2$ weighting removes the transverse convolution. The same paper proposes an analogous weighted SIDIS asymmetry,
$$
A_{SIDIS}^{\Rightarrow}(x,z,M_h^2) = S_L\, \frac{ \sum_a g_1^a(x)\, z\, G_1^{\perp a}(z,M_h^2) }{ \sum_a f_1^a(x)\, D_1^a(z,M_h^2) },
$$
which provides a universality test because the same $G_1^\perp$ enters both $e^+e^-$ and SIDIS [1712.06384].

## 4. Partial waves, resonances, and measured angular modulations in SIDIS

The transverse-momentum-sensitive dihadron final state is not a simple scalar final state. It can contain non-resonant and resonant contributions, scalar, vector, and higher partial waves, interference between partial waves, and polarization information of intermediate states such as the $\rho$ resonance. This is why the TMD dihadron cross section is more complicated than the single-pseudoscalar case even before one turns to phenomenology [1308.2989].

Within the partial-wave expansion, the paper identifies two important classes of transverse-target modulations. The collinear $(1,1)$ moment related to transversity and the Collins function coincides with
$$
\sin(\phi_R+\phi_S)\sin\vartheta .
$$
The TMD $(2,\pm2)$ partial waves are associated with
$$
\sin((1\mp2)\phi_h \pm 2 \phi_R + \phi_S)\sin^2\vartheta ,
$$
and for $\pi\pi$ dihadrons are written as
$$
\sin(3\phi_h-2\phi_R+\phi_S)\sin^2\vartheta
$$
and
$$
\sin(-\phi_h+2\phi_R+\phi_S)\sin^2\vartheta .
$$
These modulations isolate different fragmentation mechanisms tied to the internal state of the hadron pair [1308.2989].

The HERMES analysis used data taken at HERA with a $27.6~\mathrm{GeV}$ lepton beam, a transversely polarized hydrogen target, and data from 2002–2005. Event selection included
$$
Q^2 > 1~\mathrm{GeV}^2,\quad W^2 > 10~\mathrm{GeV}^2,\quad 0.023 < x < 0.4,\quad 0.2 < y < 0.95,\quad 0.2 < z < 0.8,\quad 0.05~\mathrm{GeV} < P_{h\perp} < 1.60~\mathrm{GeV}.
$$
Hadron identification used event-level algorithms, the RICH detector, and lepton-hadron separation above 98%. The analysis corrected the acceptance in parameter space of the fit amplitudes, not in histogram space, used four angular variables, and included 42 angular moments: 24 unpolarized twist-2 and twist-3 moments and 18 transverse-target Collins/Sivers-type moments. The paper also notes the Monte Carlo generator TMDGen, a TMD spectator model for dihadron fragmentation, and systematic uncertainties from smearing, radiative effects, beam-charge differences, and RICH particle ID [1308.2989].

The reported results are that the $(1,1)$ amplitudes can be used in global collinear transversity fits and help flavor separation, while the $(2,\pm2)$ amplitudes test the Lund/Artru and Gliske Gluon Radiation expectations. In particular, both $(2,\pm2)$ amplitudes are consistent with zero outside the $\rho$ peak, the $(2,-2)$ amplitude is also consistent with zero in the $\rho$ region, and the $(2,2)$ amplitude in the $\rho$-mass region for $\pi^\pm\pi^0$ is nonzero and negative. The paper interprets this as consistent with the predicted opposite-sign Collins effect for vector mesons relative to pseudoscalar mesons [1308.2989].

## 5. Phenomenological extraction of unpolarized and polarized DiFFs

A major phenomenological point is that the unpolarized DiFF $D_1$ must be determined before one can reliably extract the polarized interference function $H_1^{\sphericalangle}$. In the BELLE-based parameterization discussed for $\pi^+\pi^-$ pairs, the dependence on $(z,M_h)$ cannot be factorized: the data show that the shape in $M_h$ changes with $z$, and vice versa. The hadron-pair spectrum is also not purely continuum-like. In the invariant-mass region $2m_\pi < M_h \lesssim 1.5$ GeV, the dominant channels are $\rho\to\pi^+\pi^-$, $\omega$ production through both $\omega\to\pi^+\pi^-$ and $\omega\to\pi^+\pi^-\pi^0$ with the $\pi^0$ unobserved, and a continuum or “incoherent” $\pi^+\pi^-$ background [1012.0054].

Because BELLE provides flavor-separated Monte Carlo histograms, the parameterization distinguishes between light flavors $uds$ and charm $c$. The charm contribution is non-negligible at BELLE energies, especially in the continuum. Separate parameterizations are therefore required for $\rho$, $\omega$, and continuum contributions, and separately for $uds$ and $c$ where appropriate. The fitted unpolarized function is written schematically as
$$
f_{D_1}^a(z,M_h)=2M_h\, z^2\, D_1^a(z,M_h^2),
$$
and the fits achieved $\chi^2/\text{d.o.f.}$ values around $1.25$ for $\rho$, $1.3$ for $\omega$, $1.4$ for the $uds$ background, and $1.55$ for the charm background. The key conclusion is that a genuinely two-dimensional dependence on $(z,M_h)$ is necessary and that a factorized form does not work [1012.0054].

On the polarized side, the same paper states that $H_1^{\sphericalangle}$ had been computed only in a spectator model at that stage, while the available BELLE information was preliminary and the data were binned in $(z,\bar z)$ and $(M_h,\bar M_h)$ with large uncertainties. For that reason the extraction of a detailed functional form was not unique, and a simplified factorized ansatz was adopted,
$$
f_{H_1^{\sphericalangle}(z,M_h,\bar z,\bar M_h)\propto f(z)f(\bar z)\,g(M_h)g(\bar M_h),
$$
even though the expected physical origin is $sp$-wave interference and therefore a more structured dependence [1012.0054].

This phenomenology suggests two recurring constraints on the subject. First, transverse-momentum-sensitive dihadron fragmentation is strongly shaped by resonance structure and flavor separation. Second, clean extraction of spin-dependent DiFFs requires independent control of unpolarized fragmentation and of the functional dependence on pair kinematics.

## 6. Model realizations, QCD evolution, and later extensions

The NJL-jet model provides a fully hadronization-based realization of transverse-momentum-sensitive dihadron fragmentation. In that framework, hadronization is a sequential quark-jet process in which the initial quark emits a hadron, the remnant quark continues, and further hadrons are emitted iteratively. The key extension is that the Monte Carlo simulation tracks the transverse momenta of all produced hadrons, assigning each hadron a light-cone fraction $z_i=P_i^-/k^-$, a transverse momentum $\mathbf{P}_i^\perp$, and a mass $M_{h_i}$. The pair invariant mass is then constrained by the hadron masses, the momentum fraction ratio, and the relative transverse momentum of the pair [1310.1917].

The unpolarized DFF is interpreted as a number density,
$$
D_q^{h_1h_2}(z,M_h^2)\,\Delta z\,\Delta M_h^2 = \left\langle N_q^{h_1h_2}(z,z+\Delta z;M_h^2,M_h^2+\Delta M_h^2)\right\rangle,
$$
and extracted by simulating quark hadronization, recording all primary hadrons, forming all possible hadron pairs, binning them in $z$ and $M_h^2$, and averaging over many events. The paper uses $N_{\text{Sims}} = 10^{10}$, $\Delta z = 0.002$, and $\Delta M_h^2 = 2\,\text{MeV}^2$, and reports that about 8 hadron-emission steps are enough for convergence in the region $z \gtrsim 0.2$ [1310.1917].

A major result of the NJL-jet study is that vector-meson decays have a dramatic effect on DFFs, much larger than on ordinary single-hadron FFs, because of combinatorial enhancement in pair counting. The vector mesons included are $\rho$, $\omega$, $K^*$, and $\phi$, with only strong decays considered. The model finds visible $\rho^0$, $\omega$, $K^*$, and $\phi$ structures in the relevant channels and large differences between primary-only and full-final-state DFFs. The evolution of the DFFs from the model scale $Q_0 \simeq 0.2\,\text{GeV}$ to $Q^2 = 4\,\text{GeV}^2$ and $Q^2 = 100\,\text{GeV}^2$ is performed with the LO fixed-$M_h^2$ equation
$$
\frac{d}{d\log Q^2} D_q^{h_1h_2}(z,M_h^2,Q^2) = \frac{\alpha_s(Q^2)}{2\pi} \sum_i \int_z^1 \frac{du}{u}\, P_{qi}(u)\, D_i^{h_1h_2}\!\left(\frac{z}{u},M_h^2,Q^2\right),
$$
using QCDNUM. The qualitative effect is that the overall magnitude decreases, strength shifts to smaller $z$, and the shape in $M_h^2$ remains largely unchanged while resonance peaks remain visible [1310.1917].

Later developments extend the role of transverse-momentum-sensitive dihadron fragmentation beyond standard SIDIS and $e^+e^-$ transversity observables. In small-$x$ DIS at NLO, back-to-back dihadron production in the Color Glass Condensate generates a Sudakov double logarithm with coefficient
$$
-\frac{\alpha_s}{2\pi}\left(C_F+\frac{N_c}{2}\right),
$$
rather than the dijet value $-\frac{\alpha_s N_c}{4\pi}$. To preserve the universality of the Sudakov soft factor associated with the Weizsäcker-Williams TMD gluon distribution, the collinear fragmentation functions are promoted to TMD fragmentation functions and resummed through Collins-Soper-Sterman evolution [2405.19404].

A further extension relates $e^+e^- \to h_1 h_2 X$ factorization to the small-angle energy-energy correlator. In the collinear region, the EEC jet function can be written as a weighted integral over transverse-momentum-sensitive dihadron fragmentation functions plus a single-hadron term, and the same formalism is used to describe the confinement transition region where $Q^2 z \sim \Lambda_{\rm QCD}^2$ [2507.11495].

These developments indicate that transverse-momentum-sensitive dihadron fragmentation functions are not only tools for accessing transversity and helicity in hadronization, but also universal nonperturbative inputs in problems involving back-to-back hadron production, Sudakov resummation, and small-angle hadronic energy correlators.

Source: https://www.emergentmind.com/topics/transverse-momentum-sensitive-dihadron-fragmentation-functions