1 Introduction↩︎

From Einstein-Podolsky-Rosen’s critique questioning the completeness of quantum mechanics [1] followed by Schrödinger’s quick rejoinder and identification of entanglement [2], [3] as a defining trait of quantum mechanics in 1935, to Bell’s/CHSH inequality in the sixties [4], [5], and the first loophole-free Bell tests a few years back [6]–[8], entanglement has provided one of the clearest ways to distinguish genuine quantum correlations from correlations that can be reproduced by classical statistical mixtures in a wide range of physical systems and energy frontiers. In recent years, entanglement has begun to play an increasingly important role in high-energy physics as well [9]. A collider event is not merely a dynamical action in phase space contributing to the cross section; it also gives birth to organic quantum correlations between multiple final states [10]. It can also provide information about the correlations between different degrees of freedom for a single parton [11], [12]. The spin, momentum, and color correlations of the particles produced in a hard scattering retain detailed information about the underlying quantum nature of the interaction. When these correlations are reconstructed experimentally, collider processes become laboratories for studying quantum information in relativistic quantum field theory.
This perspective has gained significant momentum following recent studies of spin correlations in heavy-particle production [13], [14]. In particular, the observation of spin entanglement in \(t{\bar t}\) production at the LHC has demonstrated that genuine quantum correlations can survive in a strongly interacting collider environment and can be extracted from suitable angular observables [10]. These developments have opened a new direction at the interface of quantum information theory, particle physics, and quantum chromodynamics. In this emerging program, QCD radiation is not merely a complication that dilutes a clean partonic signal. Rather, it is part of the quantum dynamics that determines how entanglement is produced, transported, modified, and, in some cases, degraded.
Heavy-quark pair production provides a particularly natural setting for this program. The heavy-quark mass provides a perturbative scale, while the spin correlations of the pair can, in favorable channels, be related to measurable decay or fragmentation products. At hadron colliders, top quarks are especially clean because they decay before hadronization. At an Electron-Ion Collider, one can study a complementary class of processes in which a virtual photon probes the gluonic structure of a hadronic target. Photon-gluon fusion then produces heavy quark-(anti-quark) pairs whose spin density matrix is calculable within perturbative QCD. The EIC, therefore, offers a uniquely clean environment in which one may connect entanglement observables with the dynamics of the quarks and gluons [15]–[19].
The interface between collider physics and quantum information science provides a new perspective to explore the high energy scattering processes. Significant progress has been made in understanding quantum entanglement in electron-positron collisions [20], [21], the effects of decoherence in high-energy scattering processes [22], and the applications of quantum-information at the Electron-Ion Collider (EIC) [23]–[25].
A significant recent result was obtained by Qi, Guo, and Xiao, who studied spin entanglement and Bell nonlocality in heavy-quark pairs produced through photon-gluon fusion at an EIC [15]. They showed that, for a longitudinally polarized virtual photon, the produced quark-(anti-quark) pair is maximally entangled at leading order throughout the kinematic region. The corresponding spin-correlation matrix describes a pure Bell-type state, and the Bell-CHSH inequality is maximally violated. For transversely polarized photons, the state is generally mixed, but sizable regions of phase space still exhibit entanglement and Bell nonlocality. Their analysis identified the longitudinal photon channel as an exceptionally clean leading-order source of maximal spin entanglement in QCD at the EIC.
A follow-up development was provided by Fucilla and Hatta, in the context of diffractive heavy-quark production, by considering exclusive production through color-singlet Pomeron exchange [16]. It was found that the longitudinal photon again produces a maximally entangled and maximally Bell-violating quark-(anti-quark) pair. This agreement with the inclusive result is nontrivial, since the Pomeron represents a color-singlet gluonic exchange rather than a single gluon. It suggests that quantum-information observables may provide a new diagnostic of the QCD Pomeron and of gluonic dynamics at high energy.
This naturally raises a deeper question. Is this maximal entanglement in the longitudinal channel a robust feature of the underlying QCD dynamics, or is it an idealized leading-order result that naturally degrades once realistic QCD radiations are included?
This question is the central motivation of the present work. At leading order, the observed final state consists only of the heavy quark and anti-quark. The spin state of the pair is then pure at the amplitude level, and the longitudinal channel produces maximal entanglement. Beyond leading order, however, real gluon emission contributes. If the emitted gluon is resolved, it becomes part of the measured final state. If it is unresolved, it must be summed over. In the latter case, the experimentally relevant object is then not the pure density matrix of the full quark-antiquark-gluon system, but the reduced spin density matrix of the observed heavy-quark pair, obtained after tracing over the unobserved gluon momentum, polarization, and color. This tracing operation turns the heavy-quark pair into an open quantum subsystem.
At this point, the distinction between ordinary radiative corrections and genuine spin decoherence is essential. In the case of soft radiation, the leading soft-gluon factor is eikonal and spin independent. It changes the inclusive rate and contributes to the familiar infrared structure of QCD but, after normalization of the spin-density matrix, does not rotate or distort the leading-order spin-correlation pattern. Decoherence requires a more subtle mechanism: unresolved radiation must carry information about the spin geometry of the hard scattering. Such information first appears through subleading soft effects, where recoil and angular momentum enter the amplitude. These terms are less singular than the leading eikonal factor, but they are precisely the terms capable of changing the normalized spin density matrix.
In this paper, we show that, for the longitudinal photon, unresolved soft gluon radiation induces a dephasing mechanism. At the leading order, the spin operator of the heavy-quark pair lies in the production plane. At next to leading order, the spin-dependent soft recoil produces a small, event-by-event rotation of this plane about the direction normal to it. For a resolved emitted gluon, this rotation is coherent and unitary. However, when the gluon is unresolved, the rotation angle is averaged over the soft phase space. Even though symmetry causes the mean rotation to vanish, its variance survives. This opens a dephasing channel where the correlation normal to the production plane is preserved, while the in-plane spin correlations are suppressed. Consequently, the initially maximally entangled state slips away from maximal entanglement. In this framework, the unresolved gluon acts as an environment. It encodes information about recoil and about the orientation of the heavy-quark spin-correlation plane. Since that information is not measured, it is inaccessible to the observed subsystem. The full final state remains quantum mechanical and evolves unitarily, but the reduced state of the heavy-quark pair becomes mixed.
The paper is organized as follows. We first review the leading-order spin structure of the longitudinal photon channel and the origin of maximal entanglement as done in [15], [16]. We then discuss the real-emission contribution at next-to-leading order and explain how an unresolved gluon enters the reduced spin density matrix. Next, we isolate the spin-dependent soft-recoil component that rotates the leading-order spin-correlation plane and derive the associated dephasing channel. We obtain an analytic expression for the dephasing strength in the soft-recoil approximation and discuss its limiting behavior.

2 Leading Order spin structure↩︎

In this study, we consider the process \(\gamma^*(q)+g(p)\rightarrow q(p_1)+{\bar{q}}(p_2)+g(k)\) where the final state soft gluon is unresolved. We work in the center-of-mass frame of the observed heavy quark (anti-quark) pair in the helicity basis \(\{\hat{n},\hat{r}, \hat{k}\}\), where \(\hat{k}\) is the direction of the outgoing quark, \(\hat{n}\) is orthogonal to the production plane spanned by the incoming virtual photon direction and the outgoing quark momentum, \(\hat{r}\) is orthogonal to \(\hat{k}\) and \(\hat{n}\). The ordered set \(\{\hat{n},\hat{r}, \hat{k}\}\) is an orthonormal basis where \(\hat{k}=\left(\sin\theta\cos\phi, \sin\theta\sin\phi, \cos\theta\right)\). We define the invariant mass of \(q{\bar q}\) pair by \(M^2=\left(p_1+p_2\right)^2\), and also define the heavy-quark velocity, \(\beta=\sqrt{1-{4m^2}/{M^2}}\). The angle \(\theta\) represents the angle between the virtual-photon direction and the quark momentum. We work in the heavy-quark-pair center-of-mass frame. The observed heavy quark and antiquark momenta are \[p_1^\mu=E(1,\beta\hat{k}),\qquad p_2^\mu=E(1,-\beta\hat{k}), \qquad E=\frac{M}{2},\qquad \beta=\sqrt{1-\frac{4m^2}{M^2}} .\] Let the virtual photon define the longitudinal direction. The orthonormal triad is \[\hat{k}=(\sin\theta,0,\cos\theta),\qquad \hat{r}=(\cos\theta,0,-\sin\theta),\qquad \hat{n}=(0,1,0),\] where \(\hat{r}\) lies in the production plane and \(\hat{n}\) is normal to it. A generic soft-gluon direction relative to the outgoing heavy quark is parametrized as \[\hat{\ell} = \cos\chi\,\hat{k} +\sin\chi\cos\varphi\,\hat{r} +\sin\chi\sin\varphi\,\hat{n}, \qquad k^\mu=\omega(1,\hat{\ell}).\] Thus, \[p_1\cdot k=E\omega(1-\beta\cos\chi), \label{eq:p1dotk}\tag{1}\] which is the universal soft-recoil denominator.
In the leading order there is no final state gluon. It is therefore \(2\rightarrow2\) process as \(\gamma^*g\rightarrow q{\bar q}\). The Dirac spinors in the Pauli-spinor decomposition for the \(q{\bar q}\) pair produced in CM frame are defined as, \[\begin{align} u(p_1,\zeta)= \begin{pmatrix} a\zeta\\ b\sigma_k \zeta \end{pmatrix}, ~~~~~~~~~~~ v(p_2,\eta)=\begin{pmatrix} -b\sigma_k \eta\\ a\eta \end{pmatrix}, \end{align}\] where, \(a,b=\sqrt{E\pm m}\), \(E=M/2\) and, \(\zeta\) (and \(\eta\)) are the Pauli quark (anti-quark) spinors. The spin amplitude can be written as, \[\begin{align} {\cal M}_{\lambda_q \lambda_{\bar q}} \sim{\bar v}(p_2)\gamma^+u(p_1) =\eta^{\dagger}_{\lambda_q} X_0 ~\zeta_{\lambda_{\bar q}}, \end{align}\] where \(X_0\) is the spin-operator, acting on the Pauli-spin space, encodes the spin structure of the interaction process. For the longitudinal photon channel, in the helicity basis the leading order spin-amplitude can be represented by the spin operator, \[X_{0}={ A}(\beta,\theta)\,\sigma_r+{ B}(\beta,\theta)\,\sigma_k. \label{eq:lo95operator}\tag{2}\] Here \(\sigma_r,\sigma_k\) are the Pauli matrix projected along the basis directions \(\{\hat{r},\hat{k}\}\) and the functions \({ A}\) and \({ B}\) are real kinematic coefficients depending upon the quark velocity and the scattering angle \(\theta\). Although their overall normalization is irrelevant for the normalized spin-correlation matrix, their relative size determines the orientation of the spin operator \(X_0\) in the \(\{\hat{r},\hat{k}\}\) plane. One may notice the absence of the identity \(I\) component as well as the \(\sigma_n\) component in Eq.@eq:eq:lo95operator . Therefore, the longitudinally polarized virtual photon produces a spin operator which lies entirely in the \(\{\hat{r},\hat{k}\}\) plane or the production plane. This operator constrains the spin state of the \(q{\bar q}\) pair produced in the interaction.
It is convenient to describe the two-particle spin state in the density matrix formalism to explicitly show the spin correlation. To obtain the density matrix for a process, we take the square of the scattering amplitude, sum over the spin of the initial states and keep the spinor indices open for the final state. For the longitudinal case, the amplitude squared is proportional to \[\begin{align} {\bar u}_{\beta'}(p_1)\gamma^+v_{\beta}(p_2){\bar v_{\alpha}(p_2)}\gamma^+u_{\alpha'}(p_1). \end{align}\] The spin density matrix can further be decomposed into \(4 \times 4\) Hermitian matrix as, \[\rho=\frac{1}{4}\left[I\otimes I+ P_i \sigma_i\otimes I+\bar{ P}_jI\otimes\sigma_j+ \sum_{i,j}C_{ij}\,\sigma_i\otimes\sigma_j\right], \label{eq:density95general}\tag{3}\] where \(P\) and \(\bar{ P}\) are the particle polarization vectors of the quark and anti-quark, and \(C_{ij}\) is the spin-correlation matrix. The normalized longitudinal spin-correlation matrix at leading order written in the basis \(\{\hat{n},\hat{r},\hat{k}\}\) is found to be the form [15], \[C^L_{\rm LO}= \begin{pmatrix} 1&0&0\\ 0&C_{rr}&C_{rk}\\ 0&C_{kr}&C_{kk} \end{pmatrix}, \label{eq:clo}\tag{4}\] where, the non-zero components are given by, \[\begin{align} C_{rr}=-C_{kk} &=\frac{1-(2-\beta^2)\cos^2\theta}{1-\beta^2\cos^2\theta}, \tag{5}\\ C_{rk}=C_{kr} &=-\frac{\sqrt{1-\beta^2}\,\sin 2\theta}{1-\beta^2\cos^2\theta}. \tag{6} \end{align}\] The structure of the \(C\)-matrix reflects the planar nature of the longitudinally polarized photon interaction with the gluon. While the component \(C_{nn}=1\) indicates perfect correlation along the normal direction, the remaining \(2\times2\) block shows the nontrivial correlation in the production plane. The denominator \(1-\beta^2\cos^2\theta\) coming from the normalized squared amplitude, governs the angular and velocity dependence of the recoil.

3 Soft Emission and Universal Soft Theorem↩︎

The \(C\)-matrix in Eq.@eq:eq:clo represents a pure maximally entangled spin state. At the leading order, the spin directions of the produced \(q{\bar q}\) pair are correlated in a way that cannot be reproduced by a classical mixture of independently polarized quarks. At this order, as there is no unobserved final-state radiation, there is no information lost from the \(q{\bar q}\) system. Hence, the density matrix is pure. However, when a real gluon is emitted, it produces a larger Hilbert space, defined as, \({\cal H}_{q}\otimes {\cal H}_{\bar q}\otimes {\cal H}_{g}\). If the gluon is not observed, the relevant density matrix is not the full density matrix on this larger space and is given by the reduced density matrix given as, \(\rho_{q\bar q}^{\rm red}=\operatorname{Tr}_g\,\rho_{q\bar q g}\), where \(\operatorname{Tr}_g\) denotes a sum over the gluon polarization and color, together with an integral over the unresolved momentum region, which is the origin of decoherence.
For emission of a soft gluon of momentum \(k\), polarization \(\epsilon\), and color \(a\), the tree-level soft theorem gives \[\mathcal{M}_{q{\bar q} g}^a(k) = g_s\sum_i T_i^a \left[ \frac{p_i\cdot\epsilon^*}{p_i\cdot k} + \frac{\epsilon^*_\mu k_\nu J_i^{\mu\nu}}{p_i\cdot k} \right] \mathcal{M}_{q{\bar q}} +O(k) . \label{eq:soft-theorem}\tag{7}\] where the sum runs over the quark, anti-quark and the incoming gluon line. The angular-momentum generator decomposes into orbital and intrinsic-spin parts, \[J_i^{\mu\nu}=L_i^{\mu\nu}+S_i^{\mu\nu}, \qquad L_i^{\mu\nu} = p_i^\mu\frac{\partial}{\partial p_{i\nu}} -p_i^\nu\frac{\partial}{\partial p_{i\mu}} . \label{eq:Lmunu}\tag{8}\] The eikonal term is a scalar in the final-state spin space. This factor is spin blind: it multiplies the leading-order Pauli operator \(X_0\) by a scalar in spin space. Consequently, the leading eikonal contribution changes the real-emission rate but does not change the normalized spin-density matrix. It is therefore not the origin of decoherence in the normalized two-spin state. The part relevant for recoil-induced dephasing is therefore \[\delta X_{\rm orb}(k) = g_s\sum_i T_i^a \frac{\epsilon^*_\mu k_\nu L_i^{\mu\nu}}{p_i\cdot k} X_0(\beta,\theta). \label{eq:deltaXorb}\tag{9}\] Since, \(X_0\) depends on the hard momentum only through the scalar variables \(\beta,\theta\) and the basis matrices \(\sigma_r,\sigma_k\), the orbital generator acts by the chain rule: \[L_i^{\mu\nu}X_0 = (L_i^{\mu\nu}\beta)\partial_\beta X_0 +(L_i^{\mu\nu}\theta)\partial_\theta X_0 +A\,L_i^{\mu\nu}\sigma_r+B\,L_i^{\mu\nu}\sigma_k . \label{eq:chain-rule}\tag{10}\] Equation 10 is the key first-principles step. The subleading soft theorem does not merely multiply the Born spin structure, through \(L_i^{\mu\nu}\) it differentiates the hard kinematics and the spin basis. The term proportional to \(\partial_\beta X_0\) changes the radial position on the spin texture and can modify rates or coherent NLO spin structures. It is not a common rotation in the \((\sigma_r,\sigma_k)\) plane. The dephasing coefficient defined by 12 is controlled by the tangent projection of the remaining pieces. The stochastic rotation of the leading-order spin-correlation plane is generated by the orbital recoil term. The intrinsic spin term can contribute to the full NLO amplitude, but it does not generate the geometric tangent projection used below to define the dephasing coefficient addressed in this paper.

4 Infinitesimal rotation of the leading spin operator↩︎

The leading operator \(X_0=A\sigma_r+B\sigma_k\) may be viewed as a vector in the two-dimensional plane spanned by \(\sigma_r\) and \(\sigma_k\). A small rotation of the production plane about \(\hat{n}\) by an angle \(\delta\phi_n\) acts on the basis vectors as, \[\delta\hat{k}=\delta\phi_n\,\hat{r}, \qquad \delta\hat{r}=-\delta\phi_n\,\hat{k}. \label{eq:basis-rotation}\tag{11}\] Equivalently, the Pauli matrices transform as, \[\delta\sigma_k=\delta\phi_n\,\sigma_r, \qquad \delta\sigma_r=-\delta\phi_n\,\sigma_k.\] We note here that, depending on how one defines the basis, the overall sign of \(\delta\phi_n\) may be reversed. The physical dephasing coefficient, however, depends on \(\langle\delta\phi_n^2\rangle\) and is therefore independent of this sign convention.
The rotation of \(X^{(0)}=A\sigma_r+B\sigma_k\) about the \(\hat{n}\)-direction by an infinitesimal angle \(\delta\phi_n\) is given by the Unitary transformation, \[X_0 \to e^{-i\delta\phi_n\sigma_n/2}~X_0~e^{+i\delta\phi_n\sigma_n/2}.\] To the first order in this angle, the change in the spin operator is given by, \[\delta X_{\mathrm{rot}}=\frac{i}{2}\delta\phi_n[\sigma_n,X_0].\] Using the identity, \([\sigma_i,\sigma_j]=2i\varepsilon_{ijk}\sigma_k\) and the cyclic ordering of the basis \(\{\hat{n},\hat{r},\hat{k}\}\), the commutator in the above equation becomes, \[\begin{align} \delta X_{\rm rot} = \delta\phi_n\,(A\sigma_k-B\sigma_r) \equiv \delta\phi_n\,T_n, \end{align}\] where, \(T_n = A\sigma_k-B\sigma_r\) is in the “tangent" direction and is orthogonal to \(X_0\) in the Pauli spin operator space as \(\mathrm{Tr}(T_n^\dagger X_0) =\mathrm{Tr}\left[(A\sigma_k-B\sigma_r)(A\sigma_r+B\sigma_k)\right]=0\). The component of the real gluon emission operator \(X_g\) along the rotation tangent is isolated by projecting it onto \(T_n\) as, \[\delta\phi_n(k)= \frac{\mathrm{Tr}[T_n^\dagger X_g]}{\mathrm{Tr}[T_n^\dagger T_n]} =\frac{A\delta c_k(k)-B\delta c_r(k)}{A^2+B^2}, \label{eq:projection}\tag{12}\] where \[X_g=\delta c_0(k)I +\delta c_n(k)\sigma_n +\delta c_r(k)\sigma_r + \delta c_k(k) \sigma_k .\] The \(\sigma_n\) and \(I\) terms do not contribute to the tangent projection because their traces with \(T_n\) vanish.

4.1 Out-of-plane recoil tilt↩︎

The orbital part of the subleading soft theorem acts as a generator of changes in the hard momenta. For the outgoing heavy quark, \[S_{\mathrm{orb},1}^{(1)}X_0 =g_sT_1^a\frac{\epsilon_\mu^*k_\nu L_1^{\mu\nu}}{p_1\cdot k}X_0.\] The factor \(1/(p_1\cdot k)\) gives the characteristic soft-recoil denominator \[\frac{1}{p_1\cdot k}=\frac{1}{E\omega(1-\beta\cos\chi)}.\] The numerator \(\epsilon_\mu k_\nu L^{\mu\nu}\) generates an infinitesimal change of the angular variables defining the hard amplitude. The part relevant for dephasing is the part that changes the orientation of the production plane, not merely the magnitude of \(A\) or \(B\) but their ratios. The soft momentum component odd under reflection through the production plane is \(k_r=\omega\sin\chi\sin\phi\). A nonzero \(k_r\) tilts the recoil out of the original production plane. Therefore, the stochastic rotation angle must be proportional to \[\frac{k_r}{p_1\cdot k} \propto \frac{\sin\chi\sin\phi}{1-\beta\cos\chi}. \label{eq:angular-factor}\tag{13}\] This factor has two important properties. First, it is odd under \(\phi\to2\pi-\phi\). Hence the mean rotation vanishes. Second, its square is even and positive, so the variance survives the unresolved average. The remaining prefactor is fixed by the leading-order geometry. The rotation must vanish at threshold, where \(\beta\to0\) and the outgoing heavy quarks are produced nearly at rest. Thus, \(R_n\propto\beta\). It must also vanish in the forward and backward limits \(\theta\to0,\pi\), where the production plane is degenerate. Thus, \[R_n\propto\sin\theta.\] Finally, normalization by the leading longitudinal amplitude gives the factor \[\frac{1}{\sqrt{1-\beta^2\cos^2\theta}}.\] This is the same kinematic combination that appears in the normalized leading-order spin-correlation matrix. Combining these ingredients, the soft-recoil seed is \[\delta\phi_n^{(\perp)}= \frac{\beta\sin\theta}{\sqrt{1-\beta^2\cos^2\theta}} \frac{\sin\chi\sin\phi}{1-\beta\cos\chi}. \label{eq:Rn-derived}\tag{14}\]

4.2 In-plane recoil from the angular shift↩︎

The preceding discussion isolates the part of the out-of-plane recoil tilt of the production plane out of its Born orientation. This contribution is encoded in the recoil seed \(R_n\), and it corresponds to a direct rotation of the basis vectors \(\{\hat{r},\hat{k}\}\) about the normal direction \({\hat{n}}\). There is, however, a second recoil effect at the same order which must also be included. The leading spin operator is not only a vector in the \((\sigma_r,\sigma_k)\) operator plane; it is also a function of the Born scattering angle \(\theta\), \[\begin{align} X_0(\beta,\theta) = A(\beta,\theta)\sigma_r + B(\beta,\theta)\sigma_k . \end{align}\] Therefore, a soft recoil that changes the hard kinematics according to \(\theta \longrightarrow \theta+\delta\theta\), induces an additional variation \[\begin{align} \delta_\theta X_0=\delta\theta \partial_\theta X_0 = \delta\theta\, \left( A_\theta \sigma_r + B_\theta \sigma_k \right) \end{align}\] where, \(A_\theta\equiv \partial_\theta A\) and \(B_\theta\equiv \partial_\theta B\). This is an in-plane recoil effect: the production plane need not be tilted, but the point on the Born spin texture is shifted. To determine whether this variation contributes to dephasing, it should be projected onto the same tangent direction used above, \(T_n=A\sigma_k-B\sigma_r\). Using the projection formula, \[\begin{align} \delta\phi_n = \frac{\mathrm{Tr}\left(T_n^\dagger \delta X\right)}{\mathrm{Tr}\left(T_n^\dagger T_n\right)}, \end{align}\] the angular-shift contribution is \[\begin{align} \delta\phi_n^{(\theta)} = \frac{A B_\theta-B A_\theta}{A^2+B^2}\, \delta\theta . \label{eq:dphitheta95general} \end{align}\tag{15}\] This expression is purely geometric. It is the connection associated with motion along the Born spin-correlation curve in the \((\sigma_r,\sigma_k)\) plane. For the longitudinal-photon channel, the leading spin-correlation matrix is reproduced, up to an irrelevant common normalization, by the choice \[\begin{align} A=\sin\theta, \qquad B=\sqrt{1-\beta^2}\cos\theta . \end{align}\] Thus, Eq 15 then gives, \[\begin{align} \delta\phi_n^{(\theta)}=-\frac{\sqrt{1-\beta^2}}{1-\beta^2\cos^2\theta}\, \delta\theta. \label{eq:dphitheta95result} \end{align}\tag{16}\] The angular recoil \(\delta\theta\) is generated by the in-plane component of the unresolved gluon momentum. In the same soft-recoil approximation used to obtain \(R_n\), one may write \[\begin{align} \delta\theta = \frac{\beta\sin\theta}{\sqrt{1-\beta^2\cos^2\theta}}\, \frac{\sin\chi\cos\phi}{1-\beta\cos\chi} . \label{eq:dtheta95soft} \end{align}\tag{17}\] Thus the in-plane phase-rotation contribution becomes \[\begin{align} \delta\phi_n^{(\theta)} = - \frac{\sqrt{1-\beta^2}}{1-\beta^2\cos^2\theta} \frac{\beta\sin\theta}{\sqrt{1-\beta^2\cos^2\theta}}\, \frac{\sin\chi\cos\phi}{1-\beta\cos\chi} . \label{eq:dphitheta95soft} \end{align}\tag{18}\] The total recoil-induced rotation angle should therefore be written as \[\begin{align} \delta\phi_n(k) = \delta\phi_n^{(\perp)}(k) + \delta\phi_n^{(\theta)}(k). \label{eq:dphin95total} \end{align}\tag{19}\] The out-of-plane contribution is proportional to \(\sin\phi\), whereas the in-plane angular-shift contribution is proportional to \(\cos\phi\). Hence, for an azimuthally symmetric unresolved soft region, \[\begin{align} \left\langle \delta\phi_n^{(\perp)} \delta\phi_n^{(\theta)} \right\rangle_{\phi} \propto \int_0^{2\pi} d\phi\,\sin\phi\cos\phi=0 . \end{align}\] The two contributions therefore add incoherently in the dephasing variance, \[\begin{align} \langle \delta\phi_n^2 \rangle = \langle \left(\delta\phi_n^{(\perp)}\right)^2 \rangle + \langle \left(\delta\phi_n^{(\theta)}\right)^2 \rangle . \label{eq:dphin95variance95sum} \end{align}\tag{20}\] Since the two terms have the same soft denominator and the same unresolved energy logarithm, their relative size is determined by the kinematic prefactor as, \[\begin{align} \langle\delta\phi^2_n\rangle= \langle\left(\delta\phi_n^{(\perp)}\right)^2\rangle\left[1+ \frac{1-\beta^2}{\left(1-\beta^2\cos^2\theta\right)^2} \right], \label{eq:gamma95total95corrected} \end{align}\tag{21}\] where, \[\delta\phi_n^{(\perp)}= \frac{\beta\sin\theta}{\sqrt{1-\beta^2\cos^2\theta}} \frac{\sin\chi\sin\phi}{1-\beta\cos\chi}. \label{jmqfgdnb}\tag{22}\] Thus the in-plane angular-shift contribution is comparable to the out-of-plane tilt near threshold, while it becomes suppressed for highly relativistic heavy quarks. Away from the central region, the denominator \((1-\beta^2\cos^2\theta)^2\) can enhance the relative importance of the angular-shift contribution.

4.3 Why the spin part does not generate the dephasing coefficient?↩︎

The intrinsic spin part of the subleading soft operator is \[S_{\mathrm{spin}}^{(1)}=g_s\sum_iT_i^a \frac{\epsilon_\mu^*k_\nu S_i^{\mu\nu}}{p_i\cdot k}.\] For a fermion, using \(S^{\mu\nu}=i[\gamma^\mu,\gamma^\nu]/4\), the spatial part has the schematic Pauli reduction \[\epsilon_i k_j S^{ij}\sim (\boldsymbol{\epsilon}\times\mathbf{k})\cdot\boldsymbol{\Sigma}.\] This produces intrinsic spin insertions. Such insertions may certainly modify the amplitude. They may generate local Pauli structures proportional to \(I\), \(\sigma_n\), \(\sigma_r\), and \(\sigma_k\), and they may contribute to coherent spin-dependent NLO corrections. However, the dephasing coefficient \(\delta\phi_n\) defined above is not an arbitrary spin correction. It is the coefficient of a geometric tangent rotation of the leading operator, \(\delta X_{\mathrm{rot}}=\delta\phi_nT_n\) with \(T_n=A\sigma_k-B\sigma_r\). The defining projection is \[\delta\phi_n=\frac{\mathrm{Tr}[T_n^\dagger X_g]}{\mathrm{Tr}[T_n^\dagger T_n]}.\] The orbital part contributes to this projection because it differentiates the basis vectors \(\hat{r}\) and \(\hat{k}\) and therefore changes the orientation of the production plane. The intrinsic spin part does not act on the external angular basis. It acts on the spinor indices at fixed hard geometry. Thus, in the soft-recoil approximation used here, it does not generate the basis rotation encoded in \(T_n\). This statement may be written as \[\mathrm{Tr}\left[T_n^\dagger S_{\mathrm{spin}}^{(1)}X_0\right]_{\rm tangent}=0, \label{eq:spin-projection-zero}\tag{23}\] where the subscript emphasizes that only the component interpretable as a common geometric rotation of the leading spin-correlation plane is being projected. By contrast, \[\mathrm{Tr}\left[T_n^\dagger S_{\mathrm{orb}}^{(1)}X_0\right]_{\rm tangent}\neq0.\] The Eq 23 says only that the spin soft operator does not contribute to the particular stochastic rotation angle \(\delta\phi_n\) that defines the dephasing coefficient \(R_n\). The full NLO density matrix may include additional coherent and spin-dependent contributions from the intrinsic spin operator, which we will explore elsewhere. The dephasing channel is an open-system effect caused by averaging over unobserved recoil orientations. It is therefore controlled by the orbital recoil information carried away by the unresolved gluon. Intrinsic spin insertions are local amplitude corrections; they do not by themselves encode the stochastic orientation of the production plane.

5 Stochastic unitary evolution and the dephasing dissipators↩︎

For each fixed unresolved gluon momentum, the recoil-induced rotation is coherent. In the common-axis approximation, it is represented by \[U(k)=\exp\left[-\frac{i}{2}\delta\phi_n(k)J_n\right], \qquad J_n=\sigma_n\otimes I+I\otimes\sigma_n. \label{eq:U}\tag{24}\] This form assumes that the observed quark and antiquark spin axes are rotated by the same geometric recoil angle. A more general treatment could introduce two angles, \[J_n\delta\phi_n\to \delta\phi_q\,\sigma_n\otimes I+\delta\phi_{\bar q}\,I\otimes\sigma_n.\] We use the common-rotation approximation because the dephasing mechanism being isolated is the rotation of the observed pair’s spin-correlation plane as a whole.
For a fixed \(k\), one may now get the spin-density matrix as, \[\rho(k)=U(k)\rho^{(0)}U^\dagger(k).\] Expanding Eq. 24 to second order gives \[\begin{align} \rho(k)&=&\rho^{(0)}-\frac{i}{2}\delta\phi_n(k)[J_n,\rho^{(0)}] -\frac{1}{8}\delta\phi_n^2(k)[J_n,[J_n,\rho^{(0)}]]+\mathcal{O}(\delta\phi_n^3). \label{eq:rho-expand} \end{align}\tag{25}\] This equation is simply the Baker-Campbell-Hausdorff expansion of a unitary conjugation. Now, the unresolved average needs to be normalized. For any quantity \(O(k)\) define \[\langle O\rangle_{\mathrm{unres}} =\frac{1}{\mathcal{N}_{\mathrm{unres}}} \int_{\mathrm{unres}}\mathrm d\Phi_g\,W(k)\,O(k), \label{eq:average-def}\tag{26}\] where \(W(k)\) contains the soft-emission weight, the color factor, and the phase-space measure, and \[\mathcal{N}_{\mathrm{unres}}= \int_{\mathrm{unres}}\mathrm d\Phi_g\,W(k).\] This normalization removes the ordinary rate correction from the reduced density matrix. Without it, \(\mathrm{Tr}\rho_{\mathrm{av}}\) would not remain unity. While the virtual corrections are needed for an infra-red safe normalized observables, a present model isolates only the spin dephasing part due to soft recoil. Since \(\delta\phi_n\propto\sin\phi\), the linear term satisfies \(\langle\delta\phi_n\rangle_{\mathrm{unres}}=0. \label{eq:mean-zero}\) Therefore, the first nonvanishing correction to the normalized reduced density matrix is quadratic: \[\rho_{\mathrm{av}}=\rho^{(0)}-\frac{1}{8}\langle\delta\phi_n^2\rangle_{\mathrm{unres}} [J_n,[J_n,\rho^{(0)}]]. \label{eq:rhoav}\tag{27}\] This is the dephasing map. Using Eq. 14 , we get the dephasing factor as, see Appendix A, \[\Gamma_n(\beta,\theta)\equiv 2\langle\delta\phi_n^2\rangle= \frac{\alpha_s C_{\mathrm{eff}}}{2\pi} \ln\frac{Q}{E_{\rm cut}} \frac{\sin^2\theta}{1-\beta^2\cos^2\theta} \left[\frac{1}{2\beta}\ln\frac{1+\beta}{1-\beta}-1\right]\left[1+ \frac{1-\beta^2}{\left(1-\beta^2\cos^2\theta\right)^2} \right]. \label{eq:Gamma}\tag{28}\] Here, \(C_{\rm eff}\) denotes the effective color factor associated with the unresolved soft-emission weight and with the color projection of the observed heavy-quark pair. The cutoff \(E_{\rm cut}\) should be regarded as a schematic soft-resolution scale; in a complete phenomenological treatment it must be replaced by an infrared-safe measurement or jet-resolution function. The Eq. 27 defines the dissipative correction \[\mathcal{D}_n[\rho]=-\frac{1}{8}\langle\delta\phi_n^2\rangle[J_n,[J_n,\rho]]. \label{eq:D-box}\tag{29}\] This expression is trace preserving because the trace of any commutator vanishes: \(\mathrm{Tr}\mathcal{D}_n[\rho]=0\). It is also Hermiticity preserving because \(J_n\) is Hermitian. Using \([J,[J,\rho]]=J^2\rho+\rho J^2-2J\rho J\), Eq. 29 may be rewritten as \[\mathcal{D}_n[\rho] =\frac{\langle\delta\phi_n^2\rangle}{4} \left[J_n\rho J_n-\frac{1}{2}\{J_n^2,\rho\}\right]. \label{eq:Lindblad}\tag{30}\] Thus, the two-spin operator \(J_n\), or equivalently \(L_n=\sqrt{\gamma_n}\,J_n\) with \(\gamma_n=\frac{1}{4}\langle\delta\phi_n^2\rangle\), can be identified as the Lindblad jump operator. The corresponding soft-recoil corrected spin-correlation matrix takes the simple form \[C^{L}_{\rm soft} = \begin{pmatrix} 1 & 0 & 0\\ 0 & (1-\Gamma_n)C_{rr} & (1-\Gamma_n)C_{rk}\\ 0 & (1-\Gamma_n)C_{kr} & (1-\Gamma_n)C_{kk} \end{pmatrix}. \label{eq:Csoft}\tag{31}\] The coefficient \(\Gamma_n\) measures the dephasing strength induced by the unresolved soft recoil. Thus, the normal-axis correlation is protected, while the entire in-plane spin-correlation block is suppressed by the same soft-recoil dephasing factor. This is the characteristic signature of the unresolved-gluon-induced dephasing mechanism.

5.1 Multiple soft emissions and exponentiation of the dephasing factor↩︎

The result in Eq. (58) was obtained by keeping the first non-vanishing contribution from a single unresolved soft gluon. It is therefore a fixed-order expression in the soft-recoil variance. Now consider several unresolved soft emissions, with momenta \(k_1,k_2,\ldots,k_N\). In the soft approximation, the total recoil angle is additive, \[\Delta\phi_n=\sum_{a=1}^{N}\delta\phi_n(k_a).\] For an azimuthally symmetric unresolved region, \[\left\langle \Delta\phi_n\right\rangle=0,\] while the variance is additive, \[\left\langle \Delta\phi_n^2\right\rangle = \sum_{a=1}^{N} \left\langle \delta\phi_n^2(k_a)\right\rangle ,\] up to correlations between emissions. In the leading-logarithmic independent soft-emission approximation, these correlations are absorbed into the usual soft anomalous dimension or, equivalently here, into the effective color factor and the unresolved soft-emission weight.
For a fixed accumulated angle \(\Delta\phi_n\), the spin density matrix evolves as \[\rho(\Delta\phi_n) = \exp\left[-\frac{i}{2}\Delta\phi_n J_n\right] \rho^{(0)} \exp\left[+\frac{i}{2}\Delta\phi_n J_n\right].\] Averaging over many unresolved emissions therefore gives a dephasing channel. Equivalently, the in-plane spin-correlation components acquire the characteristic factor \[\left\langle e^{i\Delta\phi_n}\right\rangle .\] Since the accumulated recoil angle is a sum of many small unresolved contributions, its distribution is approximately Gaussian in the leading-log regime. Therefore, \[\left\langle e^{i\Delta\phi_n}\right\rangle = \exp\left[ -\frac{1}{2} \left\langle \Delta\phi_n^2\right\rangle \right].\] Using the definition \[\Gamma_n(\beta,\theta) \equiv 2\left\langle \delta\phi_n^2\right\rangle_{\rm unres},\] the exponentiated coherence-survival factor may be written as \[S_n(\beta,\theta) = \exp[-\Gamma_n(\beta,\theta)] .\] Thus the fixed-order result \[1-\Gamma_n(\beta,\theta)\] is recovered as the first term in the expansion \[e^{-\Gamma_n} = 1-\Gamma_n+\frac{\Gamma_n^2}{2}+\cdots .\]

Accordingly, the soft-recoil corrected longitudinal spin-correlation matrix can be written in the exponentiated form \[C_{\rm soft}^{L,\,{\rm exp}} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & e^{-\Gamma_n} C_{rr} & e^{-\Gamma_n} C_{rk} \\ 0 & e^{-\Gamma_n} C_{kr} & e^{-\Gamma_n} C_{kk} \end{pmatrix}.\] This form makes explicit the interpretation of \(\Gamma_n\) as a dephasing exponent. The normal-axis correlation remains protected, while the in-plane spin coherences are exponentially damped by multiple unresolved soft emissions.
We emphasize that the exponentiated expression should be understood as the leading-logarithmic soft-recoil resummation of the dephasing channel. The fixed-order expression in Eq. (58) is recovered by expanding the exponential to first order in \(\Gamma_n\). A complete phenomenological treatment would require an infrared-safe measurement or jet-resolution function and the corresponding soft anomalous dimension. In the present work these details are encoded in the effective factor \(C_{\rm eff}\) and the schematic unresolved scale \(E_{\rm cut}\) appearing in Eq. (55).

6 Estimates of Entanglement Measures and Bell Violation↩︎

We now estimate the degradation of entanglement and Bell nonlocality induced by the soft-recoil dephasing mechanism. In the preceding section, the fixed-order single-emission result was promoted to an exponentiated survival factor once multiple unresolved soft emissions are included. The in-plane spin-correlation block is therefore damped by \[\lambda(\beta,\theta) = \exp[-\Gamma_n(\beta,\theta)] , \label{eq:lambda95exp}\tag{32}\] where \(\Gamma_n\) is the soft-recoil dephasing exponent defined in Eq. (55). The exponentiated form has the fixed-order expansion \[\lambda = e^{-\Gamma_n} = 1-\Gamma_n+\mathcal{O}(\Gamma_n^2),\] and hence reproduces the single-emission result at leading order in the soft-recoil variance.
After diagonalizing the in-plane correlation block by a local spin-basis rotation, the longitudinal-channel correlation matrix is characterized by the singular values \[s_1=1, \qquad s_2=s_3=\lambda . \label{eq:singular95values95exp}\tag{33}\] The value \(s_1=1\) corresponds to the protected normal-axis correlation, while the two in-plane singular values are reduced by the same coherence-survival factor \(\lambda=e^{-\Gamma_n}\). Thus the soft-recoil mechanism produces an effectively one-parameter dephasing channel. In this approximation, the spin density matrix is locally equivalent to a Bell-diagonal state with two non-vanishing eigenvalues, \[p_1=\frac{1+\lambda}{2}, \qquad p_2=\frac{1-\lambda}{2}, \qquad p_3=p_4=0. \label{eq:bell95eigenvalues95exp}\tag{34}\] At leading order, \(\Gamma_n=0\), so that \(\lambda=1\), \(p_1=1\), and \(p_2=p_3=p_4=0\), as expected for a pure maximally entangled state.

6.1 Purity↩︎

The purity of a two-spin state with vanishing single-particle polarizations and correlation singular values \((s_1,s_2,s_3)\) is \[{\cal P} = {\rm Tr}\,\rho^2 = \frac{1}{4} \left( 1+s_1^2+s_2^2+s_3^2 \right).\] Using Eq. 33 , we obtain \[{\cal P}_{\rm exp} = \frac{1+\lambda^2}{2} = \frac{1+e^{-2\Gamma_n}}{2}. \label{eq:purity95exp}\tag{35}\] For weak dephasing, \[{\cal P}_{\rm exp} = 1-\Gamma_n+\mathcal{O}(\Gamma_n^2).\] Thus the purity decreases linearly with the dephasing exponent in the perturbative regime.

6.2 Linear Entropy↩︎

The linear entropy is defined by \[S_L=1-{\rm Tr}\,\rho^2 .\] Using Eq. 35 , one finds \[S_{L,{\rm exp}} = \frac{1-\lambda^2}{2} = \frac{1-e^{-2\Gamma_n}}{2}. \label{eq:linear95entropy95exp}\tag{36}\] For \(\Gamma_n\ll 1\), \[S_{L,{\rm exp}} = \Gamma_n+\mathcal{O}(\Gamma_n^2).\] The linear entropy therefore provides a direct leading-order measure of the mixedness generated by unresolved soft recoil.

6.3 Von Neumann Entropy↩︎

The von Neumann entropy is \[S_{\rm vN} = -\sum_i p_i\ln p_i .\] Using the eigenvalues in Eq. 34 , we obtain \[S_{{\rm vN},{\rm exp}} = -\frac{1+e^{-\Gamma_n}}{2} \ln\left( \frac{1+e^{-\Gamma_n}}{2} \right) - \frac{1-e^{-\Gamma_n}}{2} \ln\left( \frac{1-e^{-\Gamma_n}}{2} \right). \label{eq:vn95entropy95exp}\tag{37}\] This entropy vanishes when \(\Gamma_n=0\) and increases as unresolved soft radiation turns the initially pure Bell-like state into a mixed state.

6.4 Concurrence↩︎

For a Bell-diagonal two-spin state, the concurrence is \[{\cal C} = \max\{0,2p_{\rm max}-1\}.\] Since \(p_{\rm max}=p_1=(1+\lambda)/2\), we find \[{\cal C}_{\rm exp} = \max\{0,\lambda\}.\] In the present dephasing problem \(\lambda=e^{-\Gamma_n}\) is positive, and therefore \[{\cal C}_{\rm exp} = e^{-\Gamma_n}. \label{eq:concurrence95exp}\tag{38}\] For weak dephasing, \[{\cal C}_{\rm exp} = 1-\Gamma_n+\mathcal{O}(\Gamma_n^2).\] Thus the concurrence is exponentially reduced by the multiple-soft-emission dephasing exponent.

6.5 Negativity↩︎

The negativity of the same dephased Bell-type state is \[{\cal N} = \frac{1}{2}\max\{0,\lambda\}.\] Since \(\lambda=e^{-\Gamma_n}>0\), this gives \[{\cal N}_{\rm exp} = \frac{1}{2}e^{-\Gamma_n}. \label{eq:negativity95exp}\tag{39}\] At leading order, \(\Gamma_n=0\), and hence \({\cal N}_{\rm LO}=1/2\). The reduction in negativity is therefore \[\Delta{\cal N} = {\cal N}_{\rm LO}-{\cal N}_{\rm exp} = \frac{1}{2} \left( 1-e^{-\Gamma_n} \right). \label{eq:delta95negativity95exp}\tag{40}\] For \(\Gamma_n\ll1\), \[\Delta{\cal N} = \frac{\Gamma_n}{2} +\mathcal{O}(\Gamma_n^2).\]

6.6 Bell-CHSH Violation↩︎

The maximal Bell-CHSH parameter for a two-qubit state is determined by the two largest eigenvalues of \(C^T C\). Equivalently, if \(s_1,s_2,s_3\) are the singular values of the correlation matrix, then \[B_{\rm max} = 2\sqrt{s_{(1)}^2+s_{(2)}^2},\] where \(s_{(1)}\) and \(s_{(2)}\) denote the two largest singular values. For the soft-recoil corrected longitudinal state, \[s_1=1, \qquad s_2=s_3=e^{-\Gamma_n}.\] Therefore, \[B_{{\rm max},{\rm exp}} = 2\sqrt{1+e^{-2\Gamma_n}}. \label{eq:bell95exp}\tag{41}\] For weak dephasing, \[B_{{\rm max},{\rm exp}} = 2\sqrt{2} \left( 1-\frac{\Gamma_n}{2} \right) +\mathcal{O}(\Gamma_n^2).\] Thus unresolved soft radiation reduces the Bell-CHSH violation from its leading-order value \(2\sqrt{2}\). However, the violation remains above the classical bound as long as \[B_{{\rm max},{\rm exp}}>2,\] which is satisfied for any finite positive \(\lambda=e^{-\Gamma_n}\). Collecting the exponentiated multiple-soft-emission estimates, we have \[\begin{align} {\cal C}_{\rm exp} &= e^{-\Gamma_n}, \\ {\cal N}_{\rm exp} &= \frac{1}{2}e^{-\Gamma_n}, \\ {\cal P}_{\rm exp} &= \frac{1+e^{-2\Gamma_n}}{2}, \\ S_{L,{\rm exp}} &= \frac{1-e^{-2\Gamma_n}}{2}, \\ S_{{\rm vN},{\rm exp}} &= -\frac{1+e^{-\Gamma_n}}{2} \ln\left( \frac{1+e^{-\Gamma_n}}{2} \right) - \frac{1-e^{-\Gamma_n}}{2} \ln\left( \frac{1-e^{-\Gamma_n}}{2} \right), \\ B_{{\rm max},{\rm exp}} &= 2\sqrt{1+e^{-2\Gamma_n}} . \end{align}\]

These relations show that the entire degradation of entanglement, purity, and Bell violation is controlled by the single dephasing exponent \(\Gamma_n\). The normal-axis spin correlation remains protected, while the in-plane spin coherences are exponentially damped by unresolved soft radiation. In this sense, the reduced heavy-quark pair behaves as a one-parameter dephased Bell state: multiple unresolved soft emissions do not introduce several independent decoherence channels at this order, but instead exponentiate the same soft-recoil variance that appears in the single-emission calculation.

7 EIC observable: radiation-binned spin-dephasing ratio↩︎

The soft-recoil mechanism discussed above leads to a distinctive experimental signature: unresolved radiation suppresses the spin correlations lying in the production plane while leaving the correlation normal to the plane comparatively stable. This motivates an EIC-oriented observable based on comparing spin-correlation ratios in different radiation bins.

We consider semi-inclusive heavy-flavour production in deep-inelastic scattering, \[e(\ell) + p(P) \to e(\ell') + Q(p_1) + \bar Q(p_2) + X ,\] where \(q=\ell-\ell'\) is the virtual-photon momentum, \(Q^2=-q^2\), and the heavy quark and antiquark are reconstructed through heavy-flavour jets or through identified heavy-flavor hadrons. For example, one may consider \[e p \to e' + J_Q + J_{\bar Q} + X ,\] or, in more spin-sensitive channels, \[e p \to e' + \Lambda_Q + \bar\Lambda_Q + X , \qquad Q=c,b .\] The latter class is useful because the spin information of the heavy quark may be partially transferred to the heavy baryon and subsequently analyzed through its decay products. For each event, we define the heavy-quark-pair production plane using the reconstructed heavy-flavour momenta. Experimentally, the heavy-quark spin is not measured directly. Instead, one reconstructs spin-analyzing decay or fragmentation products. If \(a\) and \(b\) denote the spin analyzers associated with the heavy quark and antiquark, respectively, their angular distribution may be written schematically as \[\frac{1}{\sigma} \frac{d\sigma}{d\Omega_a d\Omega_b} = \frac{1}{(4\pi)^2} \left[ 1 + \alpha_a \alpha_b \sum_{i,j=r,n,k} C_{ij}\, \hat{a}_i \hat{b}_j \right], \label{eq:spin95analyzer95distribution}\tag{42}\] where \(\hat{a}_i\) and \(\hat{b}_j\) are the components of the analyzer directions in the event basis, and \(\alpha_a,\alpha_b\) are the corresponding analyzing powers. In ratios of spin correlations, the dependence on the analyzing powers can partially cancel, provided the same analyzer channels are used in the two radiation bins. To connect directly with unresolved soft radiation, we define an extra-radiation variable \[E_{\rm rad} = \sum_{h\notin J_Q,J_{\bar Q}} E_h^\ast . \label{eq:erad95def}\tag{43}\] The sum runs over reconstructed final-state particles or calorimeter deposits not assigned to the two heavy-flavour jets, and \(E_h^\ast\) is evaluated in a fixed analysis frame, for example the \(\gamma^\ast p\) center-of-mass frame or the reconstructed \(Q\bar Q\) rest frame. A radiation veto corresponds to requiring \[E_{\rm rad} < E_{\rm veto}.\] A smaller value of \(E_{\rm veto}\) defines a more Born-like sample with less additional radiation, while a larger value allows more unresolved or semi-soft recoil. For the present purpose it is useful to define two non-overlapping radiation bins, \[\begin{align} \text{low-radiation bin:} \qquad & 0 < E_{\rm rad} < E_{\rm cut}^{(1)}, \\ \text{higher-radiation bin:} \qquad & E_{\rm cut}^{(1)} < E_{\rm rad} < E_{\rm cut}^{(2)}, \end{align}\] with the hierarchy \[E_{\rm cut}^{(1)} < E_{\rm cut}^{(2)} \ll Q .\] The condition \(E_{\rm cut}^{(2)}\ll Q\) keeps the radiation in the soft or semi-soft regime relevant for the recoil-induced dephasing mechanism. In each bin we define the spin-dephasing ratio \[R_{\rm dep} = \frac{ \left(C_{rr}+C_{kk}\right)/2 }{ C_{nn} } . \label{eq:Rdep95def}\tag{44}\] The numerator measures the average in-plane spin correlation, while the denominator measures the correlation along the normal direction. The specific prediction of the soft-recoil dephasing mechanism is that \(C_{rr}\) and \(C_{kk}\) are suppressed by unresolved recoil, whereas \(C_{nn}\) is approximately preserved. Therefore \(R_{\rm dep}\) should decrease as one moves to a radiation bin with larger unresolved activity. The main EIC observable proposed here is the radiation-binned double ratio \[{\mathcal{D}}_{\rm EIC}^{\rm rad} = \frac{ R_{\rm dep}^{\rm higher} }{ R_{\rm dep}^{\rm low} } \, , \label{eq:D95EIC95rad}\tag{45}\] where \(R_{\rm dep}^{\rm low}\) is evaluated in the low-radiation bin and \(R_{\rm dep}^{\rm higher}\) in the higher-radiation bin. If the subleading soft recoil generates an effective dephasing exponent \(\Gamma_n\), the in-plane correlations behave schematically as \[C_{rr}, C_{kk} \longrightarrow e^{-\Gamma_n} C_{rr}, e^{-\Gamma_n} C_{kk}, \qquad C_{nn}\longrightarrow C_{nn}.\] Consequently, \[R_{\rm dep} \simeq R_{\rm dep}^{(0)} e^{-\Gamma_n},\] and the double ratio becomes \[{\mathcal{D}}_{\rm EIC}^{\rm rad} \simeq \exp\left[ -\left( \Gamma_n^{\rm higher} - \Gamma_n^{\rm low} \right) \right] . \label{eq:D95EIC95prediction}\tag{46}\] Since the higher-radiation bin contains larger unresolved recoil, one expects \[\Gamma_n^{\rm higher} > \Gamma_n^{\rm low}, \qquad \Rightarrow \qquad \mathcal{D}_{\rm EIC}^{\rm rad}<1 .\] Thus the experimentally relevant signature is not merely a reduction of the total heavy-flavour rate, but an anisotropic degradation of the spin-correlation tensor: \[C_{rr},C_{kk} \;\text{suppressed relative to}\; C_{nn} \quad \text{as unresolved radiation increases}.\]

8 Conclusion and outlook↩︎

In this work we have identified a mechanism by which QCD radiation can reduce the spin entanglement of a heavy \(Q\bar Q\) pair. At leading order, the longitudinal photon channel produces a pure two-spin state with a fixed correlation plane. Since no final-state radiation is unobserved at this order, the spin density matrix remains pure and the longitudinal channel displays maximal entanglement in the spin basis used here.
At next-to-leading order the situation changes qualitatively. Real soft-gluon emission enlarges the final-state Hilbert space. When the gluon is unresolved, the observed system is no longer the full quantum state, but the reduced density matrix obtained after summing over the emitted and unobserved gluon’s momentum, color and polarization. The leading order eikonal soft factor of this radiation is not itself responsible for spin decoherence: it is scalar in spin space and therefore cancels from normalized spin observables. The first genuine dephasing effect is instead tied to the subleading soft recoil.
The central result is that the orbital part of the subleading soft operator generates a small event-by-event rotation of the leading spin-correlation plane. The intrinsic spin part can contribute to the full NLO amplitude, and should not be discarded in a complete calculation, but it does not generate the geometric tangent projection that defines the common stochastic rotation angle in the soft-recoil channel we are discussing. After the unresolved gluon is traced over, the mean recoil angle vanishes by reflection symmetry, while the variance survives. This converts a coherent event-by-event rotation into a two-spin dephasing map generated by \[J_n=\sigma_n\otimes I+I\otimes\sigma_n .\] Equivalently, the dissipative part has the Lindblad form \[\mathcal{D}_n[\rho] =\gamma_n\left[J_n\rho J_n-\frac{1}{2}\{J_n^2,\rho\}\right],\] with \(\gamma_n\) proportional to the variance of the unresolved recoil angle. This establishes a direct connection between the subleading soft theorem and an open-quantum-system description of spin decoherence. The final normalized soft-recoil density matrix shows that the normal-axis correlation is protected, while the in-plane spin-correlation block is suppressed by the factor \(\exp(-\Gamma_n)\). The coefficient \(\Gamma_n(\beta,\theta)\) contains the heavy-quark velocity, the scattering-angle dependence and the unresolved soft energy logarithm.

9 Appendix A: Angular integrations for the soft-recoil variance↩︎

In Eq. 21 , the soft recoil angle is taken to have the angular dependence \[\delta\phi^{(\perp)}_n(k) = \frac{\beta\sin\theta}{\sqrt{1-\beta^2\cos^2\theta}} \frac{\sin\chi\sin\phi}{1-\beta\cos\chi}, \label{eq:delta-phi-app-clean}\tag{47}\] where \(\chi\) and \(\phi\) describe the unresolved gluon direction relative to the heavy-quark momentum and the \((\hat{n},\hat{r},\hat{k})\) basis. Squaring and averaging gives \[\langle(\delta\phi^{(\perp)}_n)^2\rangle = \frac{\beta^2\sin^2\theta}{1-\beta^2\cos^2\theta} \int\frac{d\omega}{\omega} \int\frac{d\Omega}{4\pi} \frac{\sin^2\chi\sin^2\phi}{(1-\beta\cos\chi)^2}. \label{eq:variance-app-clean}\tag{48}\] The energy integral over the unresolved region produces the soft logarithm \[\int_{E_{\rm cut}}^Q\frac{d\omega}{\omega}=\ln\frac{Q}{E_{\rm cut}}.\] The factorized angular integral is \[I(\beta)= \int\frac{d\Omega}{4\pi} \frac{\sin^2\chi\sin^2\phi}{(1-\beta\cos\chi)^2}.\] Using \(d\Omega=d\phi\,d\cos\chi\) and \[\int_0^{2\pi}d\phi\,\sin^2\phi=\pi,\] one obtains \[I(\beta)=\frac{1}{4}\int_{-1}^{1}du\, \frac{1-u^2}{(1-\beta u)^2}, \qquad u=\cos\chi .\] Equivalently, with \(v=1-\beta u\), \[\begin{align} I(\beta) &=\frac{1}{4\beta^3} \int_{1-\beta}^{1+\beta}dv\, \left[\frac{\beta^2-1}{v^2}+\frac{2}{v}-1\right] \\ &=\frac{1}{4\beta^3} \left[-4\beta+2\ln\frac{1+\beta}{1-\beta}\right]. \end{align}\] Therefore \[I(\beta)= \frac{1}{\beta^2} \left[\frac{1}{2\beta}\ln\frac{1+\beta}{1-\beta}-1\right]. \label{eq:Ibeta-app-clean}\tag{49}\] Combining Eqs. 48 and 49 , and absorbing the coupling, color and normalization conventions into \(\alpha_s C_{\rm eff}/(4\pi)\), gives the dephasing coefficient quoted in Eq. 28 .

10 Appendix B: Amplitudes leading to subleading Soft theorem↩︎

\[\begin{align} \nonumber i{\mathcal{M}}_1&=&\bar u(p_1)\left(i g_s \gamma^\rho T^b\right)\frac{i(-\cancel{p}_1-\cancel{k}+m)}{(p_1+k)^2-m^2}\left(i e e_q \gamma^\mu\right)\frac{i(\cancel{q}-\cancel{p}_1-\cancel{k}+m)}{(q-p_1-k)^2-m^2}\left(i g_s \gamma^\nu T^a\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\, \varepsilon_\rho^{*}(k)\,\\ \nonumber i{\mathcal{M}}_2&=&\bar u(p_1)\left(i e e_q\gamma^\mu\right)\frac{i(\cancel{p}_2+\cancel{k}-\cancel{p}+m)}{(p_2+k-p)^2-m^2}\left(i g_s \gamma^\rho T^b\right)\frac{i(\cancel{q}-\cancel{p}_1-\cancel{k}+m)}{(q-p_1-k)^2-m^2}\left(i g_s \gamma^\nu T^a\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\, \varepsilon_\rho^{*}(k)\,\\ \nonumber i{\mathcal{M}}_3&=&\bar u(p_1)\left(i e e_q \gamma^\mu\right)\frac{i(\cancel{p}_2+\cancel{k}-\cancel{p}+m)}{(p_2+k-p)^2-m^2}\left(i g_s \gamma^\nu T^a\right)\frac{i(\cancel{p}_2+\cancel{k} +m)}{(p_2+k)^2-m^2}\left(i g_s \gamma^\rho T^b\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\, \varepsilon_\rho^{*}(k)\,\\ \nonumber i{\mathcal{M}}_4&=&\bar u(p_1)\left(i g_s \gamma^\rho T^b\right)\frac{i(-\cancel{p}_1-\cancel{k}+m)}{(p_1+k)^2-m^2}\left(i g_s \gamma^\nu T^a\right)\frac{i(\cancel{p}-\cancel{p}_1-\cancel{k} +m)}{(p-p_1-k)^2-m^2}\left(i e e_q \gamma^\mu\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\, \varepsilon_\rho^{*}(k)\,\\ \nonumber i{\mathcal{M}}_5&=&\bar u(p_1)\left(i g_s \gamma^\nu T^a\right)\frac{i(\cancel{p}_2+\cancel{k}-\cancel{q}+m)}{(p_2+k-q)^2-m^2}\left(i g_s \gamma^\rho T^b\right)\frac{i(\cancel{p}-\cancel{p}_1-\cancel{k} +m)}{(p-p_1-k)^2-m^2}\left(i e e_q \gamma^\mu\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\,\varepsilon_\rho^{*}(k)\,\\ \nonumber i{\mathcal{M}}_6&=&\bar u(p_1)\left(i g_s \gamma^\nu T^a\right)\frac{i(\cancel{p}_2+\cancel{k}-\cancel{q}+m)}{(p_2+k-q)^2-m^2}\left(i e e_q \gamma^\mu\right)\frac{i(\cancel{p}_2+\cancel{k} +m)}{(p_2+k)^2-m^2}\left(i g_s \gamma^\rho T^b\right)v(p_2)\varepsilon_\mu(q)\, \varepsilon_\nu(p)\, \varepsilon_\rho^{*}(k)\,\\ \nonumber i\mathcal{M}_7 &=&\bar{u}(p_1) \left(i e e_q \gamma^\mu\right) \frac{i(\cancel{p}_2+\cancel{k}-\cancel{p}+m)}{(p_2+k-p)^2-m^2} \left(i g_s T^c\gamma^\sigma\right) v(p_2) \left(\frac{-i}{(p-k)^2} i g_s f^{abc}\right) V_{\nu\rho\sigma}^{abc}\big(p,-k,-(p-k)\big) \epsilon_\mu(q)\, \epsilon_\nu(p)\, \epsilon^{*}_{\rho}(k),\\ \nonumber i\mathcal{M}_8 &=& \bar{u}(p_1) \left(ig_s T^c\gamma^\sigma \right) \frac{i(\cancel{p}-\cancel{p}_1-\cancel{k}+m)}{(p-p_1-k)^2-m^2} \left(iee_q \gamma^\mu \right) v(p_2) \frac{-i}{(p-k)^2} \left(i g_s f^{abc}\right) V_{\nu\rho\sigma}^{abc}\big(p,-k,-(p-k)\big) \epsilon_\mu(q)\, \epsilon_\nu(p)\, \epsilon^{*}_{\rho}(k), \end{align}\]

where, \[\begin{align} V_{\nu\rho\sigma}^{abc}\big(p,-k,-(p-k)\big) &= - \Big[ (p+k)_{\sigma} g_{\nu\rho} + (-2k+p)_{\nu} g_{\rho\sigma} + (-2p+p)_{\rho} g_{\sigma\nu} \Big] . \end{align}\]

10.1 Leading Order Amplitude↩︎

These amplitudes have a soft expansion \[\begin{align} i{\cal M}_j=\frac{1}{\omega}i{\cal M}^{(-1)}_j+i{\cal M}^{(0)}_j+{\cal O}(\omega) \end{align}\]

\(\bullet~ \text{For emission from internal line}\) The amplitudes corresponding to the emission from the internal lines, do not contribute to the eikonal term.

\(\bullet~ \text{For emission from quark line}\)

\(\bullet\) For \({\cal M}_4\),

\[\begin{align} \nonumber i{\cal M}_{4}^{(-1)}&=&\bar u(p_1)\left(i g_s\gamma^\rho T^b\right)\, \frac{i\left(-\cancel{p}_1+m\right)}{\left(2p_1.k\right)^2}\,\left(i g_s \gamma^\nu T^a\right)\frac{i\left(\cancel{p}-\cancel{p}_1+m\right)}{\left(p-p_1\right)^2-m^2}\left(i e e_q\gamma^\mu\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k)\\ \end{align}\]

\(\bullet\) For \({\cal M}_1\),

\[\begin{align} \nonumber i{\cal M}_{1}^{(-1)}&=&\bar u(p_1)\left(i g_s\gamma^\rho T^b\right)\, \frac{i\left(-\cancel{p}_1+m\right)}{\left(2p_1.k\right)^2}\left(i e e_q\gamma^\mu\right)\frac{i\left(\cancel{q}-\cancel{p}_1+m\right)}{\left(q-p_1\right)^2-m^2}\left(i g_s \gamma^\nu T^a\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k)\\ \end{align}\]

Thus, we can write in the limit \(\omega\rightarrow0\), \[\begin{align} {\cal M}^{(-1)}_{4+1}=g_sT^b \left(\frac{p_1.\varepsilon_k}{p_1.k}\right){\cal M}_{LO}. \end{align}\]

\[\begin{align} i {\cal M}_{LO}={\bar u}(p_1)\left[(iee_q \gamma^\mu)S_F(q-p_1)(ig_sT^a\gamma^\nu)+(ig_sT^a\gamma^\nu)S_F(p_2-q)(iee_q \gamma^\mu\right]v(p_2)\varepsilon_\mu(q)\varepsilon_\nu(p) \end{align}\]

\(\bullet~ \text{For emission from anti- quark line}\)

\(\bullet\) For \({\cal M}_6\), \[\begin{align} \nonumber i{\cal M}_{6}^{(-1)}&=&\bar u(p_1)\left(i g_s \gamma^\nu T^a\right)\, \frac{i\left(\cancel{p}_2-\cancel{q}+m\right)}{\left(p_2-q\right)^2-m^2}\,\left(i e e_q \gamma^\mu\right)\frac{i\left(\cancel{p}_2+m\right)}{2p_2\cdot \hat{k}}\left(i g_s \gamma^\rho T^b\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k) \end{align}\]

\(\bullet\) For \({\cal M}_3\),

\[\begin{align} \nonumber i{\cal M}_{3}^{(-1)}&=&\bar u(p_1)\left(i e e_q \gamma^\mu\right)\, \frac{i\left(\cancel{p}_2-\cancel{p}+m\right)}{\left(p_2-p\right)^2-m^2}\,\left(i g_s \gamma^\nu T^a\right)\frac{i\left(\cancel{p}_2+m\right)}{2p_2\cdot \hat{k}}\left(i g_s \gamma^\rho T^b\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k)\\ \end{align}\]

Thus, we can write in the limit \(\omega\rightarrow0\), \[\begin{align} {\cal M}^{(-1)}_{6+3}=-g_s \left(\frac{p_2.\varepsilon_k}{p_2.k}\right){\cal M}_{LO}T^b. \end{align}\] \[\begin{align} i {\cal M}_{LO}={\bar u}(p_1)\left[(iee_q \gamma^\mu)S_F(p_2-p)(ig_sT^a\gamma^\nu)+(ig_sT^a\gamma^\nu)S_F(p_2-q)(iee_q \gamma^\mu\right]v(p_2)\varepsilon_\mu(q)\varepsilon_\nu(p) \end{align}\]

Thus, the total contribution to the eikonal term coming from diagrams 2,3,4 and 5 can be written as \[\begin{align} {\cal M}^{eikonal}=g_s\left(T^b\frac{p_1.\epsilon_k}{p_1.k}-\frac{p_2.\epsilon_k}{p_2.k}T^b\right){\cal M}^{LO}. \end{align}\]

\(\bullet~ \text{For emission from incoming gluon line}\)

\(\bullet\) For \({\cal M}_7\),

\[\begin{align} \nonumber i{\cal M}_{7}^{(-1)}&=&\bar u(p_1) \left(i e e_q \gamma^\mu\right)\frac{i\left(\cancel{p}_2 -\cancel{p}+m\right)}{\left(p_2-p\right)^2-m^2}\left(ig_sT^c \gamma^\sigma\right)v(p_2)\frac{i}{2p.k}\left(ig_sf^{abc}\right)V_{\nu\rho\sigma}(p,-k,-(p-k))\varepsilon_\nu(p)\varepsilon^*_\rho(k),\\ \nonumber \end{align}\]

\[\begin{align} V^{(0)}_{\nu\rho\sigma}=-\left[p_\sigma g_{\nu\rho}+ p_\nu g_{\rho \sigma} -2p_\rho g_{\sigma\nu}\right] \end{align}\]

\(\bullet\) For \({\cal M}_8\),

\[\begin{align} \nonumber i{\cal M}_{8}^{(-1)}&=&\bar u(p_1)\left(ig_sT^c \gamma^\sigma\right) \frac{i\left(\cancel{p} -\cancel{p}_1+m\right)}{\left(p-p_1\right)^2-m^2}\left(ie e_q \gamma^\mu \right)v(p_2)\frac{i}{2p.k}\left(ig_sf^{abc}\right)V_{\nu\rho\sigma}(p,-k,-(p-k))\varepsilon_\nu(p)\varepsilon^*_\rho(k),\\ \nonumber i{\cal M}_{8}^{(0)}&=&\bar u(p_1)\left(ig_sT^c \gamma^\sigma\right) \left[-\frac{-\cancel{k}}{(p-p_1)^2-m^2}+\frac{i\left(\cancel{p}-\cancel{p}_1+m\right)}{\left((p-p_1)^2-m^2\right)^2}\right]\left(i e e_q \gamma^\mu\right)v(p_2)\\ &&~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\times\frac{i}{2p.k}\left(ig_sf^{abc}\right)V_{\nu\rho\sigma}(p,-k,-(p-k))\varepsilon_\nu(p)\varepsilon^*_\rho(k) \end{align}\]

Thus, we can write in the limit \(\omega\rightarrow0\), \[\begin{align} {\cal M}^{(-1)}_{7+8}=-g_s f^{abc} \left(\frac{p.\varepsilon_k}{p.k}\right){\cal M}_{LO}. \end{align}\]

10.2 Subleading order amplitudes↩︎

\(\bullet\) For \({\cal M}_1\),

\[\begin{align} \nonumber i{\cal M}_{1}^{(0)}&=&\bar u(p_1) \left(i g_s\gamma^\rho T^b\right) \left[-\frac{i{\cancel{k}}}{2p_1\cdot k}\left(i ee_q \gamma^\mu\right) \frac{i\left(\cancel{q}-\cancel{p}_1+m\right)}{\left(q-p_1\right)^2-m^2} +i\frac{\left(-\cancel{p}_1+m\right)}{2p_1\cdot k}\left(i ee_q \gamma^\mu\right)\right.\\ \nonumber &&~~~~~~~~~~~~~~~~~~~~~~~~~~~\left. \times i\left\{-\frac{ \cancel{k}}{\left(q-p_1\right)^2-m^2}+\frac{2 k\cdot (q-p_1)}{\left(\left(q-p_1\right)^2-m^2\right)^2} \left(\cancel{q}-\cancel{p}_1+m\right) \right\}\right]\left(i g_s \gamma^\nu T^a\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k),\\ \end{align}\] where \({\cal M}_A\) is the Born amplitude and \(S^{\mu\nu}=1/4[\gamma^\mu,\gamma^\nu]\). This expression clearly shows that the first term corresponds to the orbital angular momentum part, whereas the second term contains the spin operator.

\[\begin{align} {\cal M}^{(0)}_1=ee_q g^2_s T^bT^a{\bar u}(p_1)\gamma^\mu\left[-\frac{\epsilon^*_k.p_1}{p_1.k}\frac{\partial}{\partial p^\alpha_1}+\frac{i\epsilon^*_\mu k_\nu S^{\mu\nu}}{p_1.k}\right]\left(\frac{\cancel{q}-\cancel{p}_1+m}{\left(q-p_1\right)^2-m^2}\right)\gamma^\nu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_2\), \[\begin{align} i{\cal M}_{2}^{(0)}&=&\bar u(p_1)\left(iee_q \gamma^\mu\right)\left[\frac{\cancel{p}-\cancel{p}_2+m}{(p-p_2)^2-m^2}\left(ig_s\gamma^\rho T^b\right)\frac{\cancel{p}_1-\cancel{q}+m}{(p_1-q)^2-m^2}\right]\left(ig_s\gamma^\nu T^a\right)v(p_2) \varepsilon_\mu(q)\varepsilon_\nu(p)\varepsilon^*_\rho(k).\\ \end{align}\]

\[\begin{align} {\cal M}_2=ee_q g^2_s T^bT^a{\bar u}(p_1)\gamma^\mu\left[\epsilon^*_{k\alpha}\frac{\partial}{\partial p^\alpha_1}\right]\left(\frac{\cancel{q}-\cancel{p}_1+m}{\left(q-p_1\right)^2-m^2}\right)\gamma^\nu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_3\),

\[\begin{align} \nonumber i{\cal M}_{3}^{(0)}&=&\bar u(p_1) \left(ie e_q \gamma^\mu \right) \left[\frac{i\left(\cancel{p}_2-\cancel{p}+m\right)}{\left(p_2-p\right)^2-m^2} \left(i g_s \gamma^\nu T^a\right)\frac{i{\cancel{k}}}{2p_2\cdot k} +i\left\{\frac{{\cancel{k}}}{\left(p_2-p\right)^2-m^2}-\frac{2 k\cdot (p_2-p)}{\left(\left(p_2-p\right)^2-m^2\right)^2} \left(\cancel{p}_2-\cancel{p}+m\right) \right\}\right.\\ &&~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\left.\left(i g_s \gamma^\nu T^a\right) i\frac{\cancel{p}_2+m}{2p_2\cdot k}\right]\left(i g_s \gamma^\rho T^b\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k) \\ \end{align}\]

\[\begin{align} {\cal M}^{(0)}_3=ee_q g^2_s T^aT^b{\bar u}(p_1)\gamma^\mu\left[\frac{\epsilon^*_k.p_2}{p_2.k}k_\alpha\frac{\partial}{\partial p^\alpha_2}+\frac{i\epsilon^*_\mu k_\nu S^{\mu\nu}}{p_2.k}\right]\left(\frac{\cancel{p}_2-\cancel{p}+m}{\left(p_2-p\right)^2-m^2}\right)\gamma^\nu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_4\),

\[\begin{align} \nonumber i{\cal M}_{4}^{(0)}&=&\bar u(p_1) \left(i g_s\gamma^\rho T^b\right) \left[-\frac{i{\cancel{k}}}{2p_1\cdot k} \left(i g_s \gamma^\nu T^a\right)\frac{i\left(\cancel{p}-\cancel{p}_1+m\right)}{\left(p-p_1\right)^2-m^2} +i\frac{\left(-\cancel{p}_1+m\right)}{2p_1\cdot k}\left(i g_s \gamma^\nu T^a\right)\right.\\ \nonumber &&~~~~~~~~~~~~~~~~~~~~~~~~~~~\left. \times i\left\{-\frac{ \cancel{k}}{\left(p-p_1\right)^2-m^2}+\frac{2 k\cdot (p-p_1)}{\left(\left(p-p_1\right)^2-m^2\right)^2} \left(\cancel{p}-\cancel{p}_1+m\right) \right\}\right]\left(i ee_q \gamma^\mu\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k) \\ \end{align}\]

\[\begin{align} {\cal M}^{(0)}_4=ee_q g^2_s T^bT^a{\bar u}(p_1)\gamma^\nu\left[-\frac{\epsilon^*_k.p_1}{p_1.k}\frac{\partial}{\partial p^\alpha_1}+\frac{i\epsilon^*_\mu k_\nu S^{\mu\nu}}{p_1.k}\right]\left(\frac{\cancel{p}-\cancel{p}_1+m}{\left(p-p_1\right)^2-m^2}\right)\gamma^\mu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_5\), \[\begin{align} \nonumber i{\cal M}_{5}^{(0)}&=&\bar u(p_1)\left(ig_s\gamma^\nu T^a\right)\left[i\frac{\cancel{p}_2-\cancel{q}+m}{(p_2-q)^2-m^2}\left(ig_s\gamma^\rho T^b\right)i\frac{\cancel{p}-\cancel{p}_1+m}{(p-p_1)^2-m^2}\right]\left(iee_q \gamma^\mu\right)v(p_2) \varepsilon_\mu(q)\varepsilon_\nu(p)\varepsilon^*_\rho(k).\\ \end{align}\]
\[\begin{align} {\cal M}_5=-ee_q g^2_s T^aT^b{\bar u}(p_1)\gamma^\nu\left[\epsilon^*_{k\alpha}\frac{\partial}{\partial p^\alpha_2}\right]\left(\frac{\cancel{p}_2-\cancel{q}+m}{\left(p_2-q\right)^2-m^2}\right)\gamma^\mu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_6\), \[\begin{align} \nonumber i{\cal M}_{6}^{(0)}&=&\bar u(p_1) \left(i g_s \gamma^\nu T^a\right) \left[\frac{i\left(\cancel{p}_2-\cancel{q}+m\right)}{\left(p_2-q\right)^2-m^2} \left(i e e_q \gamma^\mu\right)\frac{i{\cancel{k}}}{2p_2\cdot k} +i\left\{\frac{ {\cancel{k}}}{\left(p_2-q\right)^2-m^2}-\frac{2 k\cdot (p_2-q)}{\left(\left(p_2-q\right)^2-m^2\right)^2} \left(\cancel{p}_2-\cancel{q}+m\right) \right\}\right.\\ &&~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\left.\left(i e e_q \gamma^\mu\right) i\frac{\cancel{p}_2+m}{2p_2\cdot k}\right]\left(i g_s \gamma^\rho T^b\right)v(p_2)\epsilon_\mu(q)\epsilon_\nu(p)\epsilon_\rho^{*}(k) \\ \end{align}\]

\[\begin{align} {\cal M}^{(0)}_6=ee_q g^2_s T^aT^b{\bar u}(p_1)\gamma^\nu\left[\frac{\epsilon^*_k.p_2}{p_2.k}k_\alpha\frac{\partial}{\partial p^\alpha_2}+\frac{i\epsilon^*_\mu k_\nu S^{\mu\nu}}{p_2.k}\right]\left(\frac{\cancel{p}_2-\cancel{q}+m}{\left(p_2-q\right)^2-m^2}\right)\gamma^\mu v(p_2) \epsilon_\mu(q)\epsilon_\nu(p) \end{align}\]

\(\bullet\) For \({\cal M}_7\), \[\begin{align} i{\cal M}_{7}^{(0)}&=&\bar u(p_1) \left(ie e_q \gamma^\mu \right)\left[i\frac{\cancel{k}}{(p_2-p)^2-m^2}-i\frac{\cancel{p}_2-\cancel{p}+m}{\left((p_2-p)^2-m^2\right)^2}\right]\left(ig_sT^c \gamma^\sigma\right)v(p_2)\\ &&~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\times\frac{i}{2p.k}\left(ig_sf^{abc}\right)V_{\nu\rho\sigma}(p,-k,-(p-k))\varepsilon_\nu(p)\varepsilon^*_\rho(k) \end{align}\]

\(\bullet\) For \({\cal M}_8\), \[\begin{align} i{\cal M}_{8}^{(0)}&=&\bar u(p_1)\left(ig_sT^c \gamma^\sigma\right) \left[-\frac{i\cancel{k}}{(p-p_1)^2-m^2}+\frac{i\left(\cancel{p}-\cancel{p}_1+m\right)}{\left((p-p_1)^2-m^2\right)^2}\right]\left(i e e_q \gamma^\mu\right)v(p_2)\\ &&~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\times\frac{i}{2p.k}\left(ig_sf^{abc}\right)V_{\nu\rho\sigma}(p,-k,-(p-k))\varepsilon_\nu(p)\varepsilon^*_\rho(k) \end{align}\]

Adding the external diagrams, 1, 3, 4, and 6, we get,

\[\begin{align} {\cal M}_{ext}=g_s\left[-T^b\frac{\epsilon^*.p_1}{p_1.k}k.\partial_{p_1}{\cal M}_{Born}+\frac{\epsilon^*.p_2}{p_2.k}k.\partial_{p_2}{\cal M}_{Born}T^b+T^b\frac{i\epsilon^*_\mu k_\nu}{p_1.k}S^{\mu\nu}_1{\cal M}_{Born}+\frac{i\epsilon^*_\mu k_\nu}{p_2.k}S^{\mu\nu}_2{\cal M}_{Born}T^b\right] \end{align}\]

\[\begin{align} {\cal M}_{int}=g_s \left[T^b\epsilon^*_k.\frac{\partial {\cal M}_{A}}{\partial p_1}-\epsilon^*_k.\frac{\partial {\cal M}_{B}}{\partial p_2}T^b\right] \end{align}\]

\[\begin{align} {\cal M}^\mu={\cal M}^\mu_{ext}+{\cal M}^\mu_{int}. \end{align}\] Applying the Ward Identity, we can write, \[\begin{align} k_\mu{\cal M}^\mu=0 \end{align}\] Thus, we have \[\begin{align} k_\mu{\cal M}^\mu_{int}=- k_\mu{\cal M}^\mu_{ext}. \end{align}\] Thus, we have \[\begin{align} k_\mu{\cal M}^\mu_{int}=k_\mu\frac{\partial{\cal M}_{Born}}{\partial p_{1\alpha}}-k_\mu\frac{\partial{\cal M}_{Born}}{\partial p_{2\alpha}} \end{align}\] leading to, \[{\cal M}^{(0)}_{\rm soft} = g_s \sum_{i=1,2} {\boldsymbol{T}}_i^b \frac{\epsilon^*_{k\mu} k_\nu J_i^{\mu\nu}}{{p_i} \cdot k} {\cal M}_{Born},\]

References↩︎