The energy-momentum tensor of spin-1 hadrons: formalism
Wim Cosyn1,
Sabrina Cotogno ORCID: orcid.org/0000-0002-7440-80712,
Adam Freese3 &
Cédric Lorcé2
The European Physical Journal C volume 79, Article number: 476 (2019) Cite this article
A preprint version of the article is available at arXiv.
We provide the complete decomposition of the local gauge-invariant energy-momentum tensor for spin-1 hadrons, including non-conserved terms for the individual parton flavors and antisymmetric contributions originating from intrinsic spin. We state sum rules for the gravitational form factors appearing in this decomposition and provide relations for the mass decomposition, work balance, total and orbital angular momentum, mass radius, and inertia tensor. Generalizing earlier work, we derive relations between the total and orbital angular momentum and the Mellin moments of twist-2 and 3 generalized parton distributions, accessible in hard exclusive processes with spin-1 targets. Throughout the work, we comment on the unique features in these relations originating from the spin-1 nature of the hadron, being absent in the lower spin cases.
Avoid the common mistakes
In recent years, a lot of attention has been given to the energy-momentum tensor (EMT) as a fundamental object of study in hadronic physics and QCD, see e.g. [1, 2] and references therein. Hadronic matrix elements of the (local) EMT operator for quarks and gluons are parametrized in terms of gravitational form factors (GFFs),Footnote 1 just as hadronic matrix elements of the charge current operator are parametrized in terms of electromagnetic form factors. The GFFs encode properties that are of great interest, such as the hadronic mass and angular momentum sum rules and their spatial distributions [1, 4,5,6], and can shed light on novel properties of hadrons, such as the way that stress and shear forces are distributed within them [2, 4, 7, 8]. The topic places itself in the more general quest for a thorough understanding of the hadron structure, in which the spatial and momentum distributions of quarks and gluons play a crucial role.
Much of the recent literature on the hadron structure and QCD EMT has focused on either spin-1 / 2 [2, 5, 6, 8,9,10] or spin-0 [11,12,13] systems, see also [14,15,16,17,18] for recent lattice studies. In the former case, the proton has been the predominant object of study. This is natural, as nucleons are often understood as the primary building blocks of nuclear matter, and understanding the mass and spin decomposition of the proton is a necessary step in understanding the origin of most visible mass. In the latter case, the pion is of great interest not only for the unique role it plays in dynamical chiral symmetry breaking, which is believed to be the origin of the majority of hadron mass, but also because of its simplicity as a system and the ability to study its properties (including components of its EMT) on the lattice [19].
On the other hand, the hadronic physics community has in general dedicated little attention to the internal structure of hadrons of spin higher than 1/2. From a theoretical point of view, a full picture of higher spin hadrons and nuclei is desirable because it would serve in elucidating QCD dynamics: spin-1 (and higher spin) systems carry information on non-nucleonic degrees of freedom, i.e. the dynamics beyond quarks and gluons confined within the individual nucleons [20].
Information on spin-1 hadrons would allow us to thoroughly study such different parton contributions and dynamics in the spirit, for instance, of the theoretical calculations of the gravitational form factors for vector mesons in holographic QCD [21] and on the lattice [22]. Being nearly the only experimentally available hadronic spin-1 target, the deuteron has attracted a fair amount of attention over the past decades. It is the simplest bound state of more than one nucleon and, therefore, it has been of prime importance to unravel the nature of the nuclear binding. On the other hand, its internal structure and dynamics are the ultimate effect of the interactions between the elementary constituents, and this makes the deuteron a promising avenue towards understanding how QCD produces the force that binds nucleons together in nuclei [23]. After the first measurement by the HERMES collaboration of a tensor polarized collinear structure function of the deuteron [24], the so-called \(b_1\) function defined in [25], it became clearer that going beyond the single-nucleon formulation is needed to describe the experimental data, especially in specific regions of the parton momentum range [26,27,28,29,30].
The same arguments hold for the study of the gluonic content of higher spin hadrons, which requires once again to account for additional gluon functions in momentum and coordinate space that are exclusive to tensor polarized structures and therefore related to spin-1 or higher. This fact has stimulated a recent interest in the theoretical [31, 32] and lattice community [33, 34]. The deuteron is thus expected to play a major role in the 12 GeV program at Jefferson Lab (JLab) dedicated to spin-1 targets [35].
As a fundamental entity encoding the spatial and mechanical properties of hadrons, the EMT of a spin-1 system such as the deuteron contains much of this dynamical information that is of interest to the nuclear physics community. This information is encoded in GFFs familiar from the EMT of spin-0 and spin-1/2 systems, but also within a host of additional form factors novel to spin-1 systems. This is analogous to the spin-1 electromagnetic current containing one more form factor than the spin-1/2 current, and has a similar origin. A spin-1 system has an additional degree of freedom, which can manifest itself in higher multipole moments (in this case, a quadrupole moment) or a tensor polarization mode. A full understanding of spin-1 systems requires a complete categorization of all the independent Lorentz structures that can appear in its EMT, and an elucidation of the physical significance of the GFFs that appear with these structures.
Expressions for the decomposition of the EMT for spin-1 hadrons have appeared earlier in Refs. [21, 36, 37]. In the present work, we provide the complete EMT decomposition that also includes all non-conserved terms (appearing incompletely in [37]) and we study the properties of and the relations between the GFFs that parametrize the local operator for the EMT. More specifically, we derive new sum rules and present expressions for the mass and angular momentum decomposition of a spin-1 hadron in terms of the new structures for quarks and gluons.
The results in this paper may be relevant for the experiments at JLab and a future EIC [38] and for the proposed fixed-target projects @LHC [39], where different polarized hadrons and nuclei can be employed. Current data for spin-1 tomography is rather scarce, with HERMES having measured deeply virtual Compton scattering (DVCS) on the deuteron with both unpolarized [40] and polarized targets [41]. In these measurements, hadrons were not detected in the final state, but simulations were used to select a sample of enhanced coherent deuteron contribution. More recently, Jefferson Lab has measured deeply virtual \(\pi ^0\) production on the deuteron [42] and a recent letter of intent allows for coherent deuteron DVCS measurements [43]. The latter should also be possible in Hall B using the ALERT detector [44]. Finally, generalized distribution amplitudes (GDAs) for the rho-rho meson pair, accessible in the crossed reaction \(\gamma ^*\gamma \rightarrow \rho \rho \) [45,46,47], can be related to the rho-meson GFFs similarly to the pion case [12] and could also potentially be studied at Belle II.
This work is organized as follows. In Sect. 2, we give a full decomposition of the most general form that the EMT of a spin-1 hadron can take. This section also contains sum rules that follow immediately from energy-momentum conservation. In Sect. 3, we calculate the multipole moments of matrix elements of the spin-1 EMT. Mass and angular momentum decompositions are derived in this section, along with additional sum rules and a work-energy balance relation. Section 4 explores the connections between the EMT and Mellin moments of twist-2 and twist-3 generalized parton distributions. Finally, in Sect. 5, we summarize our results. In addition, in Appendix A the form factors counting technique is reviewed, Appendix B and C contain additional information on Lorentz projectors and the polarization bilinears useful to obtain the parametrizations of the EMT, and Appendix D displays the covariant parametrization of the GPD correlator.
Decomposition of the energy-momentum tensor
The goal of this section is to construct the most general possible parametrization for the EMT of an on-shell, spin-1 hadron in terms of GFFs. A variety of definitions exists for the EMT in QCD (for a review, see [1]), but here we work with the gauge-invariant kinetic form of the QCD EMT \(T^{\mu \nu }=T^{\mu \nu }_q+T^{\mu \nu }_g\), where
with \( \overset{\leftrightarrow }{D^{\mu }} =(\overset{\rightarrow }{\partial ^{\mu }} -\overset{\leftarrow }{\partial ^{\mu }}) -2igA^{\mu } \). Due to the presence of spin, the QCD EMT is in general not symmetric under exchange of its free Lorentz indices, with the entirety of the asymmetry in the quark contribution. The EMT \(T^{\mu \nu }\) is a conserved current, with the symmetric and antisymmetric components being separately conserved. Accordingly, we consider the general form of the EMT in two layers: the symmetric component of the EMT and the full asymmetric EMT.
Symmetric EMT
For a spin-1 system, there are only six possible independent rank-2 Lorentz structures that are symmetric, \(\mathsf {P}\)-even, \(\mathsf {T}\)-even, consistent with the hermiticity property, Lorentz-covariant, linear in each of the initial and final state polarization vectors, and conserved [21, 36]. This is fewer than the seven Lorentz structures that arise from the (1, 1) representation of the Lorentz group (see Appendix A), meaning one of those Lorentz structures is non-conserved. The conserved symmetric EMT takes the following form
$$\begin{aligned}&\langle p^{\prime }, \lambda ^{\prime } \mid T_{\mu \nu }(0) \mid p, \lambda \rangle \nonumber \\&\quad = -\,2P_\mu P_\nu \left[ (\epsilon ^{\prime *} \epsilon ) \mathcal {G}_1(t) - \frac{(\Delta \epsilon ^{\prime *} )(\Delta \epsilon )}{2M^2} \mathcal {G}_2(t)\right] \nonumber \\&\qquad -\, \frac{1}{2}(\Delta _\mu \Delta _\nu {-} \Delta ^2 g_{\mu \nu }) \left[ (\epsilon ^{\prime *} \epsilon ) \mathcal {G}_3(t) {-} \frac{(\Delta \epsilon ^{\prime *} )(\Delta \epsilon )}{2M^2} \mathcal {G}_4(t)\right] \nonumber \\&\qquad +\, P_{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon ) - \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \mathcal {G}_5(t) \nonumber \\&\qquad +\, \frac{1}{2} \left[ \Delta _{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon ) + \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \right. \nonumber \\&\qquad \left. -\, {\epsilon }^{\prime *}_{\{\mu }\epsilon _{\nu \}} \Delta ^2 - g_{\mu \nu }(\Delta \epsilon ^{\prime *})(\Delta \epsilon )\right] \mathcal {G}_6(t), \end{aligned}$$
where M is the hadron mass, \(P=(p'+p)/2\) is the average four-momentum, \(t=\Delta ^2\) with \(\Delta =p'-p\) is the four-momentum transfer, and for each four-vector a, b one has \(a_{\{\mu } b_{\nu \}}=(a_\mu b_\nu +a_\nu b_\mu )/2\). Energy and momentum must be conserved in a closed system, so in this decomposition of the symmetric EMT a sum over all partons is understood. The partial EMT for quarks and gluons does not have to be conserved however, so there are three additional independent Lorentz structures that can appear
$$\begin{aligned}&\langle p^{\prime }, \lambda ^{\prime } \mid T^a_{\mu \nu }(0) \mid p, \lambda \rangle \nonumber \\&\quad =-\,2P_\mu P_\nu \left[ (\epsilon ^{\prime *} \epsilon ) \mathcal {G}^a_1(t)-\frac{(\Delta \epsilon ^{\prime *} ) (\Delta \epsilon )}{2M^2} \mathcal {G}^a_2(t)\right] \nonumber \\&\qquad {-}\, \frac{1}{2}(\Delta _\mu \Delta _\nu {-} \Delta ^2 g_{\mu \nu }) \left[ \!(\epsilon ^{\prime *} \epsilon )\mathcal {G}^a_3(t) -\frac{(\Delta \epsilon ^{\prime *} )(\Delta \epsilon )}{2M^2} \mathcal {G}^a_4(t)\!\right] \nonumber \\&\qquad +\,P_{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon ) - \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \mathcal {G}^a_5(t)\nonumber \\&\qquad +\, \frac{1}{2} \left[ \Delta _{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon )+ \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \right. \nonumber \\&\qquad \left. -\, \epsilon _{\{\mu }^{\prime *}\epsilon _{\nu \}} \Delta ^2 - g_{\mu \nu }(\Delta \epsilon ^{\prime *})(\Delta \epsilon )\right] \mathcal {G}^a_6(t)\nonumber \\&\qquad +\,\epsilon _{\{\mu }^{\prime *}\epsilon _{\nu \}} M^2\mathcal {G}^a_7(t) + g_{\mu \nu } M^2 (\epsilon '^*\epsilon ) \mathcal {G}^a_8(t)\nonumber \\&\qquad +\, \frac{1}{2}g_{\mu \nu }(\Delta \epsilon '^*)(\Delta \epsilon ) \mathcal {G}^a_9(t), \end{aligned}$$
where \(a=q,g\). Summing over all partons, we should recover (2) which implies the following sum rules
$$\begin{aligned} \sum _{a=q,g} \mathcal {G}_i^a(t) = 0\quad \text {for}\quad i=7,8,9. \end{aligned}$$
Note that we have named the GFFs to agree with the conventions in [37], although we find an additional non-conserved pure trace GFF, in agreement with the counting in Appendix A. It should also be noted that unlike the total GFFs \(\mathcal {G}_i(t)\), the partial GFFs \(\mathcal {G}_i^a(t)\) are usually scale and scheme dependent.
Asymmetric EMT
When the constituents of a system possess intrinsic angular momentum, the EMT is in general expected to be asymmetric, see e.g. [1, 48, 49] for recent discussions. This is a simple consequence of the conservation of the generalized angular momentum \(\partial _\mu M^{\mu \alpha \beta }=0\), with \(M^{\mu \alpha \beta }=x^{\alpha } T^{\mu \beta }-x^{\beta } T^{\mu \alpha }+S^{\mu \alpha \beta }\) and where \(S^{\mu \alpha \beta }\) is the intrinsic generalized angular momentum tensor, which when combined with the conservation of the EMT implies that \(T^{\alpha \beta }-T^{\beta \alpha }=-\partial _\mu S^{\mu \alpha \beta }\). In agreement with the counting in Appendix A, we find only two antisymmetric Lorentz structures satisfying all the constraints. The most general form of the EMT is therefore
$$\begin{aligned}&\langle p^{\prime }, \lambda ^{\prime } \mid T^a_{\mu \nu }(0) \mid p, \lambda \rangle \nonumber \\&\quad =-\,2P_\mu P_\nu \left[ (\epsilon ^{\prime *} \epsilon ) \mathcal {G}^a_1(t)-\frac{(\Delta \epsilon ^{\prime *} ) (\Delta \epsilon )}{2M^2} \mathcal {G}^a_2(t)\right] \nonumber \\&\qquad {-}\, \frac{1}{2}(\Delta _\mu \Delta _\nu {-} \Delta ^2 g_{\mu \nu }) \left[ \!(\epsilon ^{\prime *} \epsilon )\mathcal {G}^a_3(t) -\frac{(\Delta \epsilon ^{\prime *} )(\Delta \epsilon )}{2M^2} \mathcal {G}^a_4(t)\!\right] \nonumber \\&\qquad +\,P_{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon ) - \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \mathcal {G}^a_5(t)\nonumber \\&\qquad +\,\frac{1}{2} \left[ \Delta _{\{\mu }\left( \epsilon ^{\prime *}_{\nu \}} (\Delta \epsilon )+ \epsilon _{\nu \}} (\Delta \epsilon ^{\prime *}) \right) \right. \nonumber \\&\qquad \left. -\, \epsilon _{\{\mu }^{\prime *}\epsilon _{\nu \}} \Delta ^2 - g_{\mu \nu }(\Delta \epsilon ^{\prime *})(\Delta \epsilon )\right] \mathcal {G}^a_6(t)\nonumber \\&\qquad +\,\epsilon _{\{\mu }^{\prime *}\epsilon _{\nu \}} M^2 \mathcal {G}^a_7(t)+ g_{\mu \nu } M^2 (\epsilon '^*\epsilon ) \mathcal {G}^a_8(t)\nonumber \\&\qquad +\, \frac{1}{2}g_{\mu \nu }(\Delta \epsilon '^*)( \Delta \epsilon ) \mathcal {G}^a_9(t)\nonumber \\&\qquad +\,P_{[\mu }\left( \epsilon ^{\prime *}_{\nu ]} (\Delta \epsilon ) - \epsilon _{\nu ]} (\Delta \epsilon ^{\prime *}) \right) \mathcal {G}^a_{10}(t)\nonumber \\&\qquad +\,\Delta _{[\mu }\left( \epsilon ^{\prime *}_{\nu ]} (\Delta \epsilon )+ \epsilon _{\nu ]} (\Delta \epsilon ^{\prime *}) \right) \mathcal {G}^a_{11}(t), \end{aligned}$$
where \(a_{[\mu } b_{\nu ]}=(a_\mu b_\nu -a_\nu b_\mu )/2\). Since one of the two new tensors is non-conserved, energy-momentum conservation imposes the additional sum rule
$$\begin{aligned} \sum _{a=q,g} \mathcal {G}_{11}^a(t) = 0. \end{aligned}$$
This is an interesting new feature of the spin-1 target, since a spin-0 target has no antisymmetric part and a spin-1/2 target has only a conserved contribution. This has to do with the fact that the intrinsic generalized angular momentum tensor for a scalar field vanishes, \(S^{\mu \alpha \beta }_0=0\), and is completely antisymmetric for a Dirac field, as \( S^{\mu \alpha \beta }_{1/2} =\frac{1}{2}\,\epsilon ^{\mu \alpha \beta \lambda } \overline{\psi }\gamma _\lambda \gamma _5\psi \) with \(\epsilon _{\,0123}=+1\). In the case of a massive vector field \(V^{\mu }\), the intrinsic generalized angular momentum tensor reads \(S^{\mu \alpha \beta }_1=-2F^{\mu [\alpha }V^{\beta ]}\), so that \(\partial _\alpha \partial _\mu S^{\mu \alpha \beta }_1\ne 0\), opening the possibility of having a non-vanishing intrinsic energy dipole moment beside intrinsic angular momentum [48].
Since we will be working with the kinetic form [1] of the QCD EMT, the QCD equations of motion imply that the antisymmetric part of the EMT can be expressed in terms of the axial-vector current as follows [1, 50, 51]
$$\begin{aligned} \overline{\psi }\gamma ^{[\mu }i\overset{\leftrightarrow }{D^{\nu }}]\psi =-\frac{1}{2}\epsilon ^{\mu \nu \rho \sigma }\partial _\rho (\overline{\psi } \gamma _\sigma \gamma _5\psi ). \end{aligned}$$
The matrix elements of the axial-vector current being parametrized as [52, 53]
$$\begin{aligned}&\langle p^{\prime }, \lambda ^{\prime } \mid \overline{\psi }(0)\gamma _\mu \gamma _5\psi (0) \mid p, \lambda \rangle \nonumber \\&\quad =-2i\epsilon _{\mu \alpha \beta P}\left( \epsilon ^{\prime *\alpha }\epsilon ^{\beta }\tilde{G}_1(t)\right. \nonumber \\&\qquad \left. +\frac{\Delta ^{\alpha }[\epsilon ^{\prime *\beta }(\Delta \epsilon ) -\epsilon ^{\beta }(\Delta \epsilon ^{\prime *})]}{M^2}\tilde{G}_2(t)\right) \end{aligned}$$
with the notation \(\epsilon _{\mu \alpha \beta P}=\epsilon _{\mu \alpha \beta \lambda }P^{\lambda }\), we find from considering the matrix elements of (7)
$$\begin{aligned} \mathcal {G}_{10}^q(t)&= -\tilde{G}_1(t)+\frac{t}{M^2}\tilde{G}_2(t), \end{aligned}$$
(9a)
$$\begin{aligned} \mathcal {G}_{11}^q(t)&= 0, \end{aligned}$$
and \(\mathcal {G}_{10}^g(t)=\mathcal {G}_{11}^g(t)=0\). Note that the vanishing of the antisymmetric part of \(T^g_{\mu \nu }\) has to do with the impossibility of writing down the gluon spin contribution in a form that is both local and gauge invariant [1]. We thus find that the antisymmetric part of the EMT for a spin-1 hadron is conserved.Footnote 2
Multipole moments of the energy-momentum tensor
Much of the interesting information about a hadron's mechanical properties that is contained in the EMT is encoded by the multipole moments of the EMT matrix elements. These include static observables such as mass, angular momentum, the inertia tensor, and so on, but additionally include information about how each of these decomposes into quark and gluon contributions. The decomposition of hadron mass and angular momentum into quark and gluon contributions – and the latter also into spin and orbital angular momentum components – has been a major focus of recent literature on the EMT. As with much of the other literature on the QCD EMT, this focus has been primarily directed towards spin-1/2 systems (predominantly the proton) and spin-0. In this section, we elaborate on the mechanical properties of spin-1 hadrons encoded by the multipole moments of their EMT, including both properties that are analogous to the lesser-spin cases and those that are new to spin-1.
Mass decomposition and balance equation
The mass decomposition and balance equation associated with a spin-1 target are obtained in terms of the following properly normalized matrix element of the EMT [5, 8]
$$\begin{aligned} \big \langle \big \langle \int \text {d}^3r\, T^a_{\mu \nu }(0,\varvec{r})\big \rangle \big \rangle&\equiv \frac{\langle p, \lambda ^{\prime } \mid \int \text {d}^3r\, T^a_{\mu \nu }(0,\vec {r}) \mid p, \lambda \rangle }{\langle p, \lambda \mid p, \lambda \rangle }\nonumber \\&=\frac{1}{2p^0}\,\langle p, \lambda ^{\prime } \mid T^a_{\mu \nu }(0) \mid p, \lambda \rangle . \end{aligned}$$
Using the covariant expression for the density matrix of a spin-1 system [54, 55] (see also App. C for more details)
$$\begin{aligned} \epsilon _\beta \epsilon ^*_\alpha =-\frac{1}{3}P_{\beta \alpha } +\frac{i}{2M}\epsilon _{\beta \alpha \mathcal {S} p}-\mathcal {T}_{\beta \alpha }, \end{aligned}$$
where the projector onto the subspace orthogonal to \(p^{\mu }\) is given by
$$\begin{aligned} P_{\mu \nu }=g_{\mu \nu }-\frac{p_\mu p_\nu }{M^2}, \end{aligned}$$
and the covariant vector and tensor polarizations by
$$\begin{aligned} \mathcal {S}^{\mu }(p)&=-\epsilon ^{\mu \rho \sigma \lambda } \text {Im}(\epsilon _\rho \epsilon ^*_\sigma )\frac{p_\lambda }{M}, \end{aligned}$$
$$\begin{aligned} \mathcal {T}_{\mu \nu }(p)&=-\frac{1}{3}P_{\mu \nu } -\text {Re}(\epsilon _\mu \epsilon ^*_\nu ), \end{aligned}$$
we find
$$\begin{aligned} \langle p, \lambda ^{\prime } \mid T^a_{\mu \nu }(0) \mid p, \lambda \rangle= & {} 2p_\mu p_\nu \left[ \mathcal {G}^a_1(0)+\frac{1}{6}\mathcal {G}^a_7(0)\right] \nonumber \\&-2g_{\mu \nu }M^2\left[ \frac{1}{2}\mathcal {G}^a_8(0)+\frac{1}{6}\mathcal {G}^a_7(0)\right] \nonumber \\&-\mathcal {T}_{\mu \nu } M^2\mathcal {G}^a_7(0). \end{aligned}$$
The first two Lorentz structures do not depend on the spin and are indeed common to all targets. The last Lorentz structure is new. It is related to the target tensor polarization and therefore does not appear in the case of spin-0 or spin-1/2 targets. Because of Poincaré invariance, the forward matrix element of the total EMT has to assume the formFootnote 3
$$\begin{aligned} \langle p, \lambda ^{\prime } \mid T_{\mu \nu }(0) \mid p, \lambda \rangle = 2p_\mu p_\nu , \end{aligned}$$
from which we conclude that
$$\begin{aligned} \sum _{a=q,g} \mathcal {G}_1^a(0) = 1 \end{aligned}$$
using the constraints in Eq. (4).
The Lorentz-invariant coefficients in Eq. (15) can be interpreted in terms of proper internal energy and pressure-volume work [5, 8]. In the target rest frame, the partial internal energy is given byFootnote 4
$$\begin{aligned} U_a=\frac{1}{2M}\langle p, \lambda ^{\prime } \mid T_a^{00}(0) \mid p, \lambda \rangle =\left[ \mathcal {G}^a_1(0)-\frac{1}{2}\mathcal {G}^a_8(0)\right] M, \end{aligned}$$
and the partial isotropic pressure–volume work by
$$\begin{aligned} W_a=\frac{\delta ^{ij}}{6M}\langle p, \lambda ^{\prime } \mid T_a^{ij}(0) \mid p, \lambda \rangle =\left[ \frac{1}{2}\mathcal {G}^a_8(0)+\frac{1}{6}\mathcal {G}^a_7(0)\right] M. \end{aligned}$$
The new feature of a spin-1 target is the presence of a partial pressure–volume work anisotropy
$$\begin{aligned} W^{ij}_a= & {} \frac{1}{2M}\left\langle p, \lambda ^{\prime } \mid T_a^{ij}(0) \mid p, \lambda \right\rangle -\delta ^{ij}W_a\nonumber \\= & {} \mathcal {T}^{ij}\left[ -\frac{1}{2}\mathcal {G}^a_7(0)\right] M \end{aligned}$$
associated with the tensor polarization. The mass decomposition then takes the form
$$\begin{aligned} M=\sum _{a=q,g}U_a \end{aligned}$$
and the balance equations read
$$\begin{aligned} \sum _{a=q,g} W_a=0,\quad \sum _{a=q,g} W^{ij}_a=0. \end{aligned}$$
Angular momentum decomposition
As explained in detail in [1, 56], higher spatial moments of the energy-momentum distribution are ambiguous if defined naively as \(\langle \langle \int \text {d}^3r\,r^j T^a_{\mu \nu }(0,\vec {r})\rangle \rangle \). The reason for this is because information about the spatial distribution is lost in the forward limit \(\Delta \rightarrow 0\). Spatial distributions can only be defined in frames where no energy is transfered to the system \(\Delta ^0=\vec {P}\cdot \vec {\Delta }/P^0=0\). In this work, we will only consider the Breit frame \(\vec {P}=\vec {0}\) where three-dimensional spatial distributions of the EMT are defined as [4, 6, 8]
$$\begin{aligned} \langle T_a^{\mu \nu }\rangle (\vec {r}){\equiv }\int \frac{\text {d}^3\Delta }{(2\pi )^3}\,e^{-i\vec {\Delta }\cdot \vec {r}}\,\frac{1}{2P^0}\left\langle \frac{\vec {\Delta }}{2}, \lambda ^{\prime } \Bigg | T_a^{\mu \nu }(0) \Bigg | {-}\frac{\vec {\Delta }}{2}, \lambda \right\rangle \end{aligned}$$
with \(P^0=\sqrt{M^2+\frac{\vec {\Delta }^2}{4}}\). The dipole moment of the spatial distribution is then given by
$$\begin{aligned}&\int \text {d}^3r\,r^j \langle T_a^{\mu \nu }\rangle (\vec {r})\nonumber \\&\quad =\left\{ -i\nabla ^j_\Delta \left[ \frac{1}{2P^0}\left\langle \frac{\vec {\Delta }}{2}, \lambda ^{\prime } \Bigg | T_a^{\mu \nu }(0) \Bigg | -\frac{\vec {\Delta }}{2}, \lambda \right\rangle \right] \right\} _{\Delta =0}.\nonumber \\ \end{aligned}$$
Using the Breit-frame expansion of the polarization four-vector bilinear derived in Appendix C, we find
$$\begin{aligned} \int \text {d}^3r\,r^j \langle T_a^{\mu \nu }\rangle (\vec {r})= & {} \frac{1}{2}g^{0\{\mu }\epsilon ^{\nu \} j\mathcal {S} 0}\left[ \mathcal {G}^a_5(0)+\frac{1}{2}\mathcal {G}^a_7(0)\right] \nonumber \\&+\frac{1}{2}g^{0[\mu }\epsilon ^{\nu ] j\mathcal {S} 0} \,\mathcal {G}^a_{10}(0). \end{aligned}$$
Clearly, the only non-vanishing dipole moment is associated with the momentum distribution \(\langle T_a^{0k}\rangle (\vec {r})\) and is orthogonal to the vector polarization of the target \(\mathcal {S}^{\mu }=(0,\vec {s})\). It simply originates from the parton orbital angular momentum (OAM)Footnote 5
$$\begin{aligned} L^i_a= & {} \epsilon ^{ijk}\int \text {d}^3r\,r^j \langle T_a^{0k}\rangle (\vec {r})\nonumber \\= & {} \frac{s^i}{2}\left[ \mathcal {G}^a_5(0)+\frac{1}{2}\mathcal {G}^a_7(0)+\mathcal {G}^a_{10}(0)\right] . \end{aligned}$$
In QCD, we then find that the parton total angular momentum (AM) is given in the target rest frame by
$$\begin{aligned} J^i_a=\epsilon ^{ijk}\int \text {d}^3r\,r^j \langle T_a^{\{0k\}}\rangle (\vec {r})=\frac{s^i}{2}\left[ \mathcal {G}^a_5(0)+\frac{1}{2}\mathcal {G}^a_7(0)\right] . \end{aligned}$$
We naturally recover \(J^i_q=L^i_q+S^i_q\) with the quark spin contribution being given by
$$\begin{aligned} S^i_q=\frac{1}{2}\int \text {d}^3r\,\langle \overline{\psi }\gamma ^i\gamma _5\psi \rangle (\vec {r})=\frac{s^i}{2}\tilde{G}_1(0)=-\frac{s^i}{2}\mathcal {G}^q_{10}(0). \end{aligned}$$
For gluons we simply have \(J^i_g=L^i_g\) because no local gauge-invariant definition of the gluon spin does exist, see e.g. [1] for a recent detailed discussion.
Beside the term proportional to \(\mathcal {G}^a_5\) already obtained in [37], we also find a contribution from the \(\mathcal {G}^a_7\) GFF describing the non-conserved part of the EMT. Interestingly, such a contribution cannot appear for spin-0 and 1/2 targets, since in these cases the non-conserved terms are necessarily of the form of a pure trace and hence decoupled from AM. Summing over all the partons, the \(\mathcal {G}_7\) contribution drops out according to Eq. (4) and we get the AM constraint
first derived in [21].
Mass radius and inertia tensor
Beside mass, two other important quantities characterizing the energy distribution can be defined, namely the mass radius and the inertia tensor. Both are expressed in terms of the following second-order moments
$$\begin{aligned} C^{ij}_a= & {} \int \text {d}^3r\,r^i r^j \langle T_a^{00}\rangle (\vec {r})\nonumber \\= & {} \left\{ -\nabla ^i_\Delta \nabla ^j_\Delta \left[ \frac{1}{2P^0} \left\langle \frac{\vec {\Delta }}{2}, \lambda ^{\prime } \Bigg | T_a^{00}(0) \Bigg | -\frac{\vec {\Delta }}{2}, \lambda \right\rangle \right] \right\} _{\Delta =0}.\nonumber \\ \end{aligned}$$
The mass radius defined as \(R_M=\sqrt{\langle r^2\rangle }\) with
$$\begin{aligned} \langle r^2\rangle =\frac{1}{M}\sum _{a=q,g} C^{ij}_a\delta ^{ij} \end{aligned}$$
gives an idea of the spatial extension of the energy distribution. The inertia tensor [58, 59] defined as
$$\begin{aligned} I^{ij}_a= C^{kl}_a\left( \delta ^{kl}\delta ^{ij} -\delta ^{ik}\delta ^{jl}\right) \end{aligned}$$
allows one to determine the moment of inertia \(I_a^{\vec {n}} =I^{ij}_a n^in^j\) of the system about an arbitrary axis \(\vec {n}\) passing through the center of mass, which coincides in the Breit frame with the origin [48]. It is related to the mass quadrupole moment
$$\begin{aligned} Q^{ij}_a= C^{kl}_a\left( \delta ^{ik}\delta ^{jl}-\frac{1}{3}\delta ^{kl} \delta ^{ij}\right) =- I^{kl}_a\left( \delta ^{ik}\delta ^{jl}-\frac{1}{3} \delta ^{kl}\delta ^{ij}\right) \end{aligned}$$
which measures the deviation from a spherical distribution of the energy.
Using once more the Breit-frame expansion of the polarization four-vector bilinear derived in Appendix C, we find
$$\begin{aligned} C^{ij}_a=\frac{1}{M}\left[ \delta ^{ij}\left( \mathcal {A}_a(0) +\frac{1}{3}\mathcal {B}_a(0)\right) -\mathcal {T}^{ij}\mathcal {B}_a(0)\right] \end{aligned}$$
$$\begin{aligned} \mathcal {A}_a(t)= & {} -\frac{1}{4}\left( \mathcal {G}^a_1(t)+2\mathcal {G}^a_3(t) +\frac{1}{2}\mathcal {G}^a_8(t)\right) \nonumber \\&+2M^2\frac{\mathrm d}{\mathrm dt} \left( \mathcal {G}^a_1(t)-\frac{1}{2}\mathcal {G}^a_8(t)\right) , \end{aligned}$$
(35a)
$$\begin{aligned} \mathcal {B}_a(t)= & {} -\mathcal {G}^a_1(t)-\mathcal {G}^a_2(t)+\mathcal {G}^a_5(t) +\frac{1}{2}\mathcal {G}^a_6(t)\nonumber \\&+\frac{1}{4}\mathcal {G}^a_7(t) +\frac{1}{2}\mathcal {G}^a_8(t)-\frac{1}{2}\mathcal {G}^a_9(t). \end{aligned}$$
(35b)
The squared mass radius, inertia tensor and mass quadrupole moment are then given by
$$\begin{aligned}&\langle r^2\rangle =\frac{1}{M^2}\sum _{a=q,g}\left[ 3\mathcal {A}_a(0) +\mathcal {B}_a(0)\right] , \end{aligned}$$
$$\begin{aligned}&I^{ij}_a=\frac{1}{M}\left[ 2\delta ^{ij}\left( \mathcal {A}_a(0)+\frac{1}{3} \mathcal {B}_a(0)\right) +\mathcal {T}^{ij}\mathcal {B}_a(0)\right] , \end{aligned}$$
$$\begin{aligned}&Q^{ij}_a=-\frac{1}{M}\mathcal {T}^{ij}\mathcal {B}_a(0). \end{aligned}$$
As expected, the quadrupole moment in spin-1 hadrons is different from zero due to the presence of the tensor polarization.
GFFs as moments of GPDs
The connection between the EMT and partonic distributions has long been a topic of consideration (see e.g. [60,61,62]). Generalized parton distributions (GPDs) in particular allow for GFFs to be extracted from their Mellin moments. Since GPDs parametrize the non-perturbative structure contributing to hard reactions such as deeply virtual compton scattering (DVCS) and virtual meson production, they are the most promising avenue for experimentally exploring the form factors appearing in the EMT decomposition.
Leading-twist GPDs have been extensively studied for their polynomiality relations [62]. A specific case allows the second Mellin moments of helicity-independent twist-2 GPDs to be related to the non-trace GFFs appearing in the symmetric component of the EMT. Such relations have been studied not only for spin-0 [12, 63, 64] and spin-1/2 [8, 61, 62], but also spin-1 systems [21, 22, 37, 65].
Beyond leading twist, the Penttinen–Polyakov–Shuvaev–Strikman (PPSS) sum rule [66] relates the second moment of a twist-3 GPD to the orbital angular momentum carried by quarks in spin-1/2 hadrons, and it has also been shown [1] that twist-4 GPDs contain information about the non-conserved GFF in a nucleon. In this section we will derive a spin-1 analogue of the PPSS sum rule.
We proceed to derive sum rules for second Mellin moments of GPDs up to twist-3.Footnote 6 Consider the quark and gluon vector correlators.Footnote 7
$$\begin{aligned} V_{\lambda \lambda ^{\prime }}^{\mu ,q}= & {} \langle p^{\prime }, \lambda ^{\prime } | \mathcal {O}^{\mu }_{qV}| p, \lambda \rangle \nonumber \\= & {} \frac{1}{2}\int _{-\infty }^{\infty } \frac{\mathrm {d}\kappa }{2\pi } e^{i(Pn)\kappa x}\nonumber \\&\times \left\langle p^{\prime }, \lambda ^{\prime } | \overline{\psi }\left( -\frac{n\kappa }{2}\right) \gamma ^{\mu } \left[ {-}\frac{n\kappa }{2},\frac{n\kappa }{2}\right] \psi \left( \frac{n\kappa }{2}\right) | p, \lambda \right\rangle , \end{aligned}$$
$$\begin{aligned} V_{\lambda \lambda ^{\prime }}^{\mu ,g}= & {} \langle p^{\prime }, \lambda ^{\prime } | \mathcal {O}^{\mu }_{gV}| p, \lambda \rangle \nonumber \\= & {} \frac{\delta ^{\phantom {\{}\mu }_{\{\alpha } n^{\phantom {\mu }}_{\beta \}}}{x(Pn)}\int _{-\infty }^{\infty } \frac{\mathrm {d}\kappa }{2\pi } e^{i(Pn)\kappa x}\nonumber \\&\times \left\langle p^{\prime }, \lambda ^{\prime } | \mathrm {Tr}\left\{ \left[ \frac{n\kappa }{2}, -\frac{n\kappa }{2}\right] \right. \right. \nonumber \\&\times \left. \left. F^{\alpha \lambda }\left( -\frac{n\kappa }{2}\right) \left[ {-}\frac{n\kappa }{2},\frac{n\kappa }{2}\right] F_\lambda ^{\phantom {\lambda }\beta } \left( \frac{n\kappa }{2}\right) \right\} | p, \lambda \right\rangle , \end{aligned}$$
where n is a light-like four-vector and [y, z] denotes a straight Wilson line joining the spacetime points y and z. These correlators enter the description of deeply virtual Compton scattering and can be parametrized up to twist 3 as followsFootnote 8 [53, 68]
$$\begin{aligned} V_{\lambda \lambda '}^{\mu ,a}= & {} -\frac{P^{\mu }}{(Pn)} (\epsilon ^{\prime *} \epsilon ) \; H_1^a + \frac{\epsilon ^{\prime *\mu } (\Delta \epsilon ) -\epsilon ^{\mu } (\Delta \epsilon ^{\prime *} )}{2(Pn)}\, H_2^a\nonumber \\&+ \frac{P^{\mu }}{(Pn)}\frac{(\Delta \epsilon ^{\prime *} )(\Delta \epsilon )}{2M^2} \, H_3^a - \frac{\epsilon ^{\prime *\mu } (\Delta \epsilon ) +\epsilon ^{\mu } (\Delta \epsilon ^{\prime *} )}{2(Pn)}\, H_4^a \nonumber \\&+\left[ M^2 \frac{\epsilon ^{\prime *\mu }(n \epsilon ) +\epsilon ^{\mu } (n \epsilon ^{\prime *}) }{2(Pn)^2} + \frac{1}{3} \frac{P^{\mu }}{(Pn)}(\epsilon ^{\prime *} \epsilon ) \right] H_5^a\nonumber \\&-\frac{\Delta ^{\mu }_T}{(Pn)}(\epsilon ^{\prime *}\epsilon ) \; G_1^a -\frac{\Delta ^{\mu }_T}{(Pn)}\frac{(n\epsilon ^{\prime *}) (\Delta \epsilon ) -(n\epsilon ) (\Delta \epsilon ^{\prime *} )}{2(Pn)}\, G_2^a \nonumber \\&+\frac{\Delta ^{\mu }_T}{(Pn)}\frac{(\Delta \epsilon ^{\prime *} ) (\Delta \epsilon )}{2M^2} \; G_3^a\nonumber \\&-\frac{\Delta ^{\mu }_T}{(Pn)} \frac{(n\epsilon ^{\prime *}) (\Delta \epsilon ){+}(n\epsilon ) (\Delta \epsilon ^{\prime *} )}{2(Pn)}\, G_4^a\nonumber \\&+\frac{\Delta ^{\mu }_T}{(Pn)}\left[ M^2 \frac{(n \epsilon ^{\prime *}) (n \epsilon )}{(Pn)^2} + \frac{1}{3}(\epsilon ^{\prime *} \epsilon ) \right] G_5^a \nonumber \\&+\frac{\epsilon ^{\prime *\mu }_T (\Delta \epsilon ) -\epsilon ^{\mu }_T (\Delta \epsilon ^{\prime *} )}{2(Pn)} \, G_6^a - \frac{\epsilon ^{\prime *\mu }_T (\Delta \epsilon ) {+}\epsilon ^{\mu }_T (\Delta \epsilon ^{\prime *} )}{2(Pn)}\, G_7^a\nonumber \\&+ M^2 \frac{\epsilon ^{\prime *\mu }_T(n \epsilon ) +\epsilon ^{\mu }_T (n \epsilon ^{\prime *}) }{2(Pn)^2}\, G_8^a\nonumber \\&+ M^2 \frac{\epsilon ^{\prime *\mu }_T(n \epsilon ) -\epsilon ^{\mu }_T (n \epsilon ^{\prime *}) }{2(Pn)^2}\, G_9^a . \end{aligned}$$
The first five terms (\(H^a_i\)) correspond to the twist-2 GPDs, and the remaining nine (\(G^a_i\)) are purely twist-3. In the quark sector, the twist-3 GPDs satisfy the relation \(\int \mathrm {d}x\,G_i^q=0\) as a consequence of the charge current conservation. We suppressed the dependence of the GPDs on the parton longitudinal momentum x, longitudinal momentum transfer \(\xi =-(\Delta n)/2(Pn)\), and squared momentum transfer \(t=\Delta ^2\) for conciseness of notation, and made use of the Sudakov decomposition of four-vectors \((n^2=\bar{n}^2=0,\, n\bar{n}=1)\)
$$\begin{aligned} z^{\mu } = (z n) \bar{n}^{\mu } + (z \bar{n}) n^{\mu } + z_T^{\mu }. \end{aligned}$$
The second Mellin moment of the light-front string operators are related to the EMT up to twist 3 as follows
$$\begin{aligned} \int _{-1}^1\mathrm {d}x\,x \mathcal {O}^{\mu }_{qV}&=\frac{1}{4(P n)^{2}}\overline{\psi }(0) \gamma ^{\mu }(i\overset{\leftrightarrow }{D}n) \psi (0)=\frac{T_q^{\mu n}}{2(Pn)^2}, \end{aligned}$$
$$\begin{aligned} \int _{-1}^1\mathrm {d}x\,x\mathcal {O}^{\mu }_{gV}&=\frac{1}{(Pn)^2}\mathrm {Tr}\left[ F^{\mu \lambda }(0) F_{\lambda n}(0)\right] =\frac{T_g^{\mu n}}{2(Pn)^2}. \end{aligned}$$
Taking the off-forward matrix element on both sides allows us to relate seven of the GFFs to moments of leading-twist vector GPDs [37, 65]. Comparing the Mellin moment of Eq. (40) with the decomposition in Eq. (2) in the symmetric frame \(P^{\mu }_T=0\), we find the following relations for quarks and gluons at twist 2
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x \left[ H^{a}_1(x,\xi ,t) - \frac{1}{3}H^{a}_5(x,\xi ,t)\right] = \mathcal {G}^{a}_1(t) + \xi ^2 \mathcal {G}^{a}_3(t), \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x H^{a}_2(x,\xi ,t) = \mathcal {G}^{a}_5(t), \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x H^{a}_3(x,\xi ,t) = \mathcal {G}^{a}_2(t) + \xi ^2 \mathcal {G}^{a}_4(t), \end{aligned}$$
(43c)
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x H^{a}_4(x,\xi ,t) = \xi \mathcal {G}^{a}_6(t), \end{aligned}$$
(43d)
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x H^{a}_5(x,\xi ,t) = -\frac{t}{4M^2}\mathcal {G}^{a}_6(t) + \frac{1}{2}\mathcal {G}^{a}_7(t). \end{aligned}$$
(43e)
Since the GFFs appear as second Mellin moments of GPDs, they are special cases of generalized form factors, which correspond to arbitrary moments of GPDs. We note the following correspondence between the GFFs as defined in this work, and the \(s=2\) generalized form factors appearing in Ref. [65]:
$$\begin{aligned}&A^{a}_{2,0}(t) = \mathcal {G}^{a}_1(t), \end{aligned}$$
$$\begin{aligned}&B^{a}_{2,0}(t) = \mathcal {G}^{a}_5(t), \end{aligned}$$
$$\begin{aligned}&C^{a}_{2,0}(t) = \mathcal {G}^{a}_2(t), \end{aligned}$$
$$\begin{aligned}&2 D^{a}_{2,1}(t) = -\mathcal {G}^{a}_6(t), \end{aligned}$$
$$\begin{aligned}&E^{a}_{2,1}(t) = -\frac{t}{4M^2}\mathcal {G}^{a}_6(t) + \frac{1}{2}\mathcal {G}^{a}_7(t), \end{aligned}$$
$$\begin{aligned}&4 F^{a}_{2}(t) = \mathcal {G}^{a}_3(t), \end{aligned}$$
(44f)
$$\begin{aligned}&4 H^{a}_{2}(t) = \mathcal {G}^{a}_4(t). \end{aligned}$$
We also find the following relations for quarks at twist 3
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_1(x,\xi ,t) = -\frac{\xi }{2} \mathcal {G}^{q}_3(t)\nonumber \\&\quad =-\frac{1}{4}\frac{\partial }{\partial \xi } \int _{-1}^1\mathrm {d}x\, x \left[ H^{q}_1(x,\xi ,t) - \frac{1}{3}H^{q}_5(x,\xi ,t)\right] , \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_2(x,\xi ,t) = 0, \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_3(x,\xi ,t) = -\frac{\xi }{2} \mathcal {G}^{q}_4(t)\nonumber \\&\quad =-\frac{1}{4}\frac{\partial }{\partial \xi }\int _{-1}^1 \mathrm {d}x\, x H^{q}_3(x,\xi ,t), \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_4(x,\xi ,t) = -\frac{1}{4}\mathcal {G}^{q}_6(t)\nonumber \\&\quad =-\frac{1}{4}\frac{\partial }{\partial \xi } \int _{-1}^1\mathrm {d}x\, x H^{q}_4(x,\xi ,t), \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_6(x,\xi ,t) =-\frac{1}{2}\left[ \mathcal {G}^{q}_5(t) + \mathcal {G}^{q}_{10}(t)\right] \nonumber \\&\quad =-\frac{1}{2}\int _{-1}^1\mathrm {d}x\, x H^{q}_2(x,\xi ,t) +\frac{1}{2}\left[ \tilde{G}_1(t) -\frac{t}{M^2}\tilde{G}_2(t)\right] , \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, x G^{q}_7(x,\xi ,t) =-\frac{\xi }{2}\mathcal {G}^{q}_6(t)\nonumber \\&\quad =-\frac{1}{2}\int _{-1}^1\mathrm {d}x\, x H^{q}_4(x,\xi ,t), \end{aligned}$$
$$\begin{aligned}&\int _{-1}^1\mathrm {d}x\, xG^{q}_8(x,\xi ,t) =0, \end{aligned}$$
(45i)
where we have used Eq. (9b).
Based on Eq. (27), we find that the total quark or gluon AM in a state with maximal vector polarization along the z-direction can be expressed in terms of twist-2 GPDs as follows
$$\begin{aligned} J^z_a=\int _{-1}^1\mathrm {d}x\, \frac{x}{2}\left[ H^{a}_2(x,0,0)+H^{a}_5(x,0,0)\right] , \end{aligned}$$
which is nothing but Ji's relation [61] for spin-1 targets. Summing over quark and gluon contributions, we recover the spin sum rule derived in [21]
$$\begin{aligned} J^z=\sum _{a=q,g}\int _{-1}^1\mathrm {d}x\, \frac{x}{2}\,H^{a}_2(x,0,0)=1. \end{aligned}$$
Unlike the case of spin-1/2 targets [66, 69, 70], the quark OAM (26) requires not only a pure twist-3 GPD but also a twist-2 GPD
$$\begin{aligned} L^z_q=\int _{-1}^1\mathrm {d}x\, x\left[ \frac{1}{2}H^{q}_5(x,0,0)-G^{q}_6(x,0,0)\right] . \end{aligned}$$
This twist-2 GPD contribution is associated with the tensor polarization and is therefore absent in the case of spin-1/2 targets.
Quark and gluon contributions to mass and pressure-volume work involve trace terms and hence twist-4 GPDs. Only the partial pressure-volume work anisotropy can be related to a twist-2 GPD
$$\begin{aligned} W^{ij}_a=-\mathcal {T}^{ij}M\int _{-1}^1\mathrm {d}x\, xH^{a}_5(x,0,0). \end{aligned}$$
Summing over quarks and gluons, the mass sum rule and the balance equations imply the following constraints
$$\begin{aligned} \sum _{a=q,g}\int _{-1}^1\mathrm {d}x\, xH^{a}_1(x,0,0)&=1, \end{aligned}$$
$$\begin{aligned} \sum _{a=q,g}\int _{-1}^1\mathrm {d}x\, xH^{a}_5(x,0,0)&=0. \end{aligned}$$
These relations can be interpreted as statements of energy-momentum conservation for collinear parton distribution functions (PDFs). Using the notation of [31, 71],Footnote 9 the unpolarized and tensor-polarized PDFs are given respectively by \(f_1^a(x) = H_1^a(x,0,0)\) and \(f^a_{1LL}(x) = H^a_5(x,0,0)\). In the case of quarks, they enter the deep inelastic structure functions \(F_1\) and \(b_1\) at leading order and leading twist as
$$\begin{aligned} F_1(x,Q^2)= & {} \frac{1}{2}\sum _{q} e_q^2 \left[ H_1^q(x,0,0;\mu ^2=Q^2) \right. \nonumber \\&\left. - H_1^q(-x,0,0;\mu ^2=Q^2)\right] , \end{aligned}$$
$$\begin{aligned} b_1(x,Q^2)= & {} \frac{1}{2} \sum _{q} e_q^2 \left[ H_5^q(x,0,0;\mu ^2=Q^2) \right. \nonumber \\&\left. - H_5^q(-x,0,0;\mu ^2=Q^2)\right] . \end{aligned}$$
The corresponding gluon PDFs mix with the quark ones and contribute to the DIS structure functions at higher order in the \(\alpha _s\) perturbative expansion. Thus Eq. (50) is a statement of the momentum sum rule for PDFs. The collinear structure functions for a scattering off a tensor polarized targets were first introduced in [20, 25, 72] and the separate contributions of quarks and gluons to Eq. (51) were previously discussed in [73, 74].
In this work, we found the most general form that the asymmetric, gauge-invariant kinetic energy-momentum tensor (EMT) of a spin-1 hadron can take. Expressions were given for both the full EMT and the partial EMT due to a single parton type. We explored the physical meaning of the gravitational form factors appearing in this EMT, including sum rules imposed by conservation of momentum and angular momentum, the decomposition of spin-1 hadron mass, and multipole moments of the EMT. We also explored connections between the gravitational form factors and other functions describing partonic structure, such as axial form factors and generalized parton distributions up to twist three.
The spin-1 EMT was found to contain many more gravitational form factors than the corresponding spin-0 or spin-1/2 EMTs. A total of 11 form factors are present in the EMT decomposition, with 9 of these in the symmetric part and 2 in the antisymmetric part. Among them, 6 structures have no analogues in the lower-spin cases and are related to the presence of tensor polarization modes. They contribute to features new to spin-1 hadrons such as a quadrupole moment and possibly a non-zero intrinsic energy dipole moment aside from the intrinsic angular momentum.
The structure of the spin-1 EMT is rich, and there remains much to be explored. The pressure and shear force distributions encoded within it, and how these differ from the simpler spin-0 case, are worthy of detailed study. It is also worth investigating how the EMT of a composite spin-1 hadron compares to that of an elementary spin-1 particle, such as a photon or one of the heavy electroweak gauge bosons. These topics will be the subject of future work, along with illustrative model calculations.
Experimentally, measurements of coherent hard exclusive processes with deuteron targets are possible at JLab and the future EIC with forward detectors. Extraction of the chiral-even vector GPDs from these measurements would then constrain the deuteron gravitational form factors through the Mellin moments of these GPDs. Similarly, extraction of GDAs for the rho-rho meson pair from the crossed reaction \(\gamma ^*\gamma \rightarrow \rho \rho \) at Belle would constrain the rho meson gravitational form factors.
Shortly after the present work was completed and made available, the independent work [75] appeared, dealing with the EMT of spin-1 hadrons. The results are consistent with ours and the document contains an especially useful comparison between the different nomenclature for the GFFs present in the literature.
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no data sets were generated or analyzed during the current study.]
We call the form factors "gravitational" since the EMT is usually understood as the source of gravitational interactions. However, we do not measure in practice these form factors through gravity and we do not know the exact form that a theory of quantum gravity will take. Thus, it is unclear whether gravitation sees the symmetric, Belinfante-improved EMT as in general relativity or an asymmetric EMT as in Einstein-Cartan theory [3]. Despite this, we refer to all form factors appearing in either the symmetric or asymmetric EMT as "gravitational" form factors.
We note that this conclusion holds for the gauge-invariant kinetic EMT, but may not hold for the (non-gauge-invariant) canonical EMT. This is because the gluon contribution to the canonical EMT is asymmetric, and due to the fact that \({\partial _\alpha \partial _\mu S^{\mu \alpha \beta }_1\ne 0}\), the divergences of the symmetric and antisymmetric parts of the canonical gluon EMT are not expected to separately vanish.
This is consistent with the total four-mometum of the system being given by \(p^{\mu }=\langle \langle \int \text {d}^3r\, T^{0\mu }(0,\vec {r})\rangle \rangle \).
One can of course avoid having recourse to the rest frame and use instead the projectors \(p^{\alpha } p^{\beta }/M^2\) and \(P_{\alpha \beta }\) [5].
Since we consider the local gauge-invariant EMT, we are dealing with the kinetic form of OAM, see e.g. [1, 57] for more details.
It is possible to consider up to twist-4 for quarks and twist-6 for gluons, but contributions beyond twist-3 are less promising for future phenomenological studies.
The Ji convention [62] is used in this work for the gluon GPDs to simplify all the following formulas. Note that \(2x H_{g,\mathrm{Ji}} = H_{g,\mathrm{Diehl}}\) [67].
For the complete covariant parametrization, see Appendix D.
Note that Ref. [31] uses \((Pn)=1\).
Interestingly, since \(\sigma ^{\mu \nu }=\frac{i}{2}[\gamma ^{\mu },\gamma ^{\nu }]\) we can interpret the Clifford algebra \(\{\gamma ^{\mu },\gamma ^{\nu }\}=2g^{\mu \nu }\mathbb I\) as the condition of vanishing quadrupole.
E. Leader, C. Lorcé, Phys. Rept. 541, 163 (2014). arXiv:1309.4235 [hep-ph]
ADS Google Scholar
M.V. Polyakov, P. Schweitzer, Int. J. Mod. Phys. A 33, 1830025 (2018). arXiv:1805.06596 [hep-ph]
E. Cartan, Ann. Sci. Ecole Norm. Sup. 40, 325 (1923)
M.V. Polyakov, Phys. Lett. B 555, 57 (2003). arXiv:hep-ph/0210165 [hep-ph]
ADS MathSciNet Google Scholar
C. Lorcé, Eur. Phys. J. C 78, 120 (2018a). arXiv:1706.05853 [hep-ph]
C. Lorcé, L. Mantovani, B. Pasquini, Phys. Lett. B 776, 38 (2018). arXiv:1704.08557 [hep-ph]
M.V. Polyakov, A.G. Shuvaev, (2002). arXiv:hep-ph/0207153 [hep-ph]
C. Lorcé, H. Moutarde, A.P. Trawiński, Eur. Phys. J. C 79, 89 (2019). arXiv:1810.09837 [hep-ph]
M.V. Polyakov, H.-D. Son, JHEP 09, 156 (2018). [JHEP 18, 156 (2020)]. arXiv:1808.00155 [hep-ph]
I.V. Anikin, (2019). arXiv:1902.00094 [hep-ph]
J. Hudson, P. Schweitzer, Phys. Rev. D 96, 114013 (2017). arXiv:1712.05316 [hep-ph]
S. Kumano, Q.-T. Song, O.V. Teryaev, Phys. Rev. D 97, 014020 (2018). arXiv:1711.08088 [hep-ph]
K. Tanaka, Phys. Rev. D 98, 034009 (2018). arXiv:1806.10591 [hep-ph]
C. Alexandrou, M. Constantinou, K. Hadjiyiannakou, K. Jansen, C. Kallidonis, G. Koutsou, A. Vaquero Avilés-Casco, C. Wiese, Phys. Rev. Lett. 119, 142002 (2017). arXiv:1706.02973 [hep-lat]
Y.-B. Yang, M. Gong, J. Liang, H.-W. Lin, K.-F. Liu, D. Pefkou, P. Shanahan, Phys. Rev. D 98, 074506 (2018a). arXiv:1805.00531 [hep-lat]
Y.-B. Yang, J. Liang, Y.-J. Bi, Y. Chen, T. Draper, K.-F. Liu, Z. Liu, Phys. Rev. Lett. 121, 212001 (2018b). arXiv:1808.08677 [hep-lat]
P.E. Shanahan, W. Detmold, Phys. Rev. D 99, 014511 (2019). arXiv:1810.04626 [hep-lat]
P. E. Shanahan, W. Detmold (2018). arXiv:1810.07589 [nucl-th]
T. Horn, C.D. Roberts, J. Phys. G 43, 073001 (2016). arXiv:1602.04016 [nucl-th]
R.L. Jaffe, A. Manohar, Phys. Lett. B 223, 218 (1989)
Z. Abidin, C.E. Carlson, Phys. Rev. D 77, 095007 (2008). arXiv:0801.3839 [hep-ph]
W. Detmold, D. Pefkou, P.E. Shanahan, Phys. Rev. D 95, 114515 (2017). [Erratum: PhysRevD.99.019904]. arXiv:1703.08220 [hep-lat]
W. Boeglin, M. Sargsian, Int. J. Mod. Phys. E 24, 1530003 (2015). arXiv:1501.05377 [nucl-ex]
A. Airapetian et al., (HERMES), Phys. Rev. Lett. 95, 242001 (2005). arXiv:hep-ex/0506018 [hep-ex]
P. Hoodbhoy, R.L. Jaffe, A. Manohar, Nucl. Phys. B 312, 571 (1989)
N.N. Nikolaev, W. Schafer, Phys. Lett. B 398, 245 (1997). [Erratum: Phys. Lett.B407,453(1997)]. arXiv:hep-ph/9611460 [hep-ph]
K. Bora, R.L. Jaffe, Phys. Rev. D 57, 6906 (1998). arXiv:hep-ph/9711323 [hep-ph]
J. Edelmann, G. Piller, W. Weise, Z. Phys. A 357, 129 (1997). arXiv:nucl-th/9701026 [nucl-th]
G.A. Miller, Phys. Rev. C 89, 045203 (2014). arXiv:1311.4561 [nucl-th]
W. Cosyn, Y.-B. Dong, S. Kumano, M. Sargsian, Phys. Rev. D 95, 074036 (2017). arXiv:1702.05337 [hep-ph]
D. Boer, S. Cotogno, T. van Daal, P.J. Mulders, A. Signori, Y.-J. Zhou, JHEP 10, 013 (2016). arXiv:1607.01654 [hep-ph]
S. Cotogno, T. van Daal, P.J. Mulders, JHEP 11, 185 (2017). arXiv:1709.07827 [hep-ph]
W. Detmold, P.E. Shanahan, Phys. Rev. D 94, 014507 (2016). [Erratum: Phys. Rev.D95,no.7,079902(2017)]. arXiv:1606.04505 [hep-lat]
F. Winter, W. Detmold, A.S. Gambhir, K. Orginos, M.J. Savage, P.E. Shanahan, M.L. Wagman, Phys. Rev. D 96, 094512 (2017). arXiv:1709.00395 [hep-lat]
K. Slifer, Proceedings, Tensor Polarized Solid Target Workshop: Newport News, USA, March 10–12, 2014. J. Phys. Conf. Ser. 543, 012003 (2014)
B.R. Holstein, Phys. Rev. D 74, 084030 (2006). arXiv:gr-qc/0607051 [gr-qc]
S.K. Taneja, K. Kathuria, S. Liuti, G.R. Goldstein, Phys. Rev. D 86, 036008 (2012). arXiv:1101.0581 [hep-ph]
A. Accardi et al., Eur. Phys. J. A 52, 268 (2016). arXiv:1212.1701 [nucl-ex]
D. Kikoła, M.G. Echevarria, C. Hadjidakis, J.-P. Lansberg, C. Lorcé, L. Massacrier, C.M. Quintans, A. Signori, B. Trzeciak, Few Body Syst. 58, 139 (2017). arXiv:1702.01546 [hep-ex]
A. Airapetian et al., (HERMES), Nucl. Phys. B 829, 1 (2010). arXiv:0911.0095 [hep-ex]
A. Airapetian et al., (HERMES), Nucl. Phys. B 842, 265 (2011). arXiv:1008.3996 [hep-ex]
M. Mazouz et al., (Jefferson Lab Hall A), Phys. Rev. Lett. 118, 222002 (2017). arXiv:1702.00835 [hep-ex]
A. Camsonne, M. Carmignotto, E. Fuchey, R. Montgomery, A. Puckett, Z. Ye, Z. Ye, et al., (Jefferson Lab Hall A SBS), Deeply virtual compton scattering using the tagged deeply inelastic scattering experimental setup (2018). LOI12-18-002
W. Armstrong et al. (2017). arXiv:1708.00888 [nucl-ex]
P. Achard et al., (L3),Phys. Lett. B 568, 11 (2003). arXiv:hep-ex/0305082 [hep-ex]
I.V. Anikin, B. Pire, O.V. Teryaev, Phys. Rev. D 69, 014018 (2004). arXiv:hep-ph/0307059 [hep-ph]
I.V. Anikin, B. Pire, O.V. Teryaev, Phys. Lett. B 626, 86 (2005). arXiv:hep-ph/0506277 [hep-ph]
C. Lorcé, Eur. Phys. J. C 78, 785 (2018b). arXiv:1805.05284 [hep-ph]
W. Florkowski, R. Ryblewski, (2018). arXiv:1811.04409 [nucl-th]
I.I. Balitsky, V.M. Braun, Nucl. Phys. B 311, 541 (1989)
C. Lorcé, JHEP 08, 045 (2015). arXiv:1502.06656 [hep-ph]
T. Frederico, E.M. Henley, S.J. Pollock, S. Ying, Phys. Rev. C 46, 347 (1992)
E.R. Berger, F. Cano, M. Diehl, B. Pire, Phys. Rev. Lett. 87, 142302 (2001). arXiv:hep-ph/0106192 [hep-ph]
H. Joos, G. Kramer, Z. Phys. 178, 542 (1964)
D. Zwanziger, Phys. Rev. B 137, 1535 (1965)
B.L.G. Bakker, E. Leader, T.L. Trueman, Phys. Rev. D 70, 114001 (2004). arXiv:hep-ph/0406139 [hep-ph]
K.-F. Liu, C. Lorcé, Eur. Phys. J. A 52, 160 (2016). arXiv:1508.00911 [hep-ph]
K.R. Symon, Mechanics, 3rd edn. (Addison-Wesley Pub. Co, Reading, MA, 1971)
H. Goldstein, Mechanics, 2nd edn. (Addison-Wesley, Reading, MA, 1980)
R.L. Jaffe, A. Manohar, Nucl. Phys. B 337, 509 (1990)
X.-D. Ji, Phys. Rev. Lett. 78, 610 (1997). arXiv:hep-ph/9603249 [hep-ph]
X.-D. Ji, J. Phys. G G24, 1181 (1998). arXiv:hep-ph/9807358 [hep-ph]
M.V. Polyakov, C. Weiss, Phys. Rev. D 60, 114017 (1999). arXiv:hep-ph/9902451 [hep-ph]
H.-D. Son, H.-C. Kim, Phys. Rev. D 90, 111901 (2014). arXiv:1410.1420 [hep-ph]
W. Cosyn, A. Freese, B. Pire, Phys. Rev. D 99, 094035 (2019). arXiv:1812.01511 [hep-ph]
M. Penttinen, M.V. Polyakov, A.G. Shuvaev, M. Strikman, Phys. Lett. B 491, 96 (2000). arXiv:hep-ph/0006321 [hep-ph]
M. Diehl, Phys. Rept. 388, 41 (2003). arXiv:hep-ph/0307382 [hep-ph]
I.V. Anikin, R.S. Pasechnik, B. Pire, O.V. Teryaev, Eur. Phys. J. C 72, 2055 (2012). arXiv:1112.1849 [hep-ph]
D.V. Kiptily, M.V. Polyakov, Eur. Phys. J. C 37, 105 (2004). arXiv:hep-ph/0212372 [hep-ph]
Y. Hatta, S. Yoshida, JHEP 10, 080 (2012). arXiv:1207.5332 [hep-ph]
A. Bacchetta, P.J. Mulders, Phys. Rev. D 62, 114004 (2000). arXiv:hep-ph/0007120 [hep-ph]
X. Artru, M. Mekhfi, Z. Phys. C 45, 669 (1990)
A.V. Efremov, O.V. Teryaev, Sov. J. Nucl. Phys. 36, 557 (1982). [Yad. Fiz.36,950(1982)]
A.V. Efremov, O.V. Teryaev, in Dubna Deuteron - 93 (1994). arXiv:hep-ph/9910555 [hep-ph]
M.V. Polyakov, B.-D. Sun, (2019). arXiv:1903.02738 [hep-ph]
X.-D. Ji, R.F. Lebed, Phys. Rev. D 63, 076005 (2001). arXiv:hep-ph/0012160 [hep-ph]
S. Liuti, in Talk at INT Workshop (February 6–17, 2012). http://www.int.washington.edu/talks/WorkShops/int_12_49W/People/Liuti_S/Liuti.pdf
S. Liuti, in 23rd International Symposium on Spin Physics (SPIN 2018) Ferrara, Italy, September 10–14, 2018 (2019) arXiv:1903.08760 [hep-ph]
J. Krause, J. Math. Phys. 18, 889 (1977)
D.E. Fahnline, Am. J. Phys. 50, 818 (1982). https://doi.org/10.1119/1.12748
C. Lorcé, Phys. Rev. D 97, 016005 (2018). arXiv:1705.08370 [hep-ph]
W.N. Polyzou, W. Glöckle, H. Witala, Few Body Syst. 54, 1667 (2013). arXiv:1208.5840 [nucl-th]
We are grateful to Bernard Pire and Maxim Polyakov for useful discussions. The work of SC and CL was supported by the Agence Nationale de la Recherche under the projects No. ANR-18-ERC1-0002 and ANR-16-CE31-0019. AF was supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, Contract No. DE-AC02-06CH11357 and an LDRD initiative at Argonne National Laboratory under Project No. 2017-058-N0.
Department of Physics and Astronomy, Ghent University, Proeftuinstraat 86, 9000, Ghent, Belgium
Wim Cosyn
Centre de Physique Théorique, École polytechnique, CNRS, Université Paris-Saclay, 91128, Palaiseau, France
Sabrina Cotogno & Cédric Lorcé
Argonne National Laboratory, Lemont, IL, 60439, USA
Adam Freese
Sabrina Cotogno
Cédric Lorcé
Correspondence to Sabrina Cotogno.
Appendix A: Form factor counting
Let us consider the total number of GFFs we may expect to appear in the decomposition of the spin-1 EMT. This will depend on the EMT definition we use. First, let us consider the most general possible rank-2 tensor that is even under charge conjugation. It decomposes into several representations of the Lorentz group, namely (0, 0), \((1,0)\oplus (0,1)\), and (1, 1), which are respectively 1-, 6-, and 9-dimensional representations. These representations have \(J^{PC}\) quantum number decompositions given in Table 1.
Table 1 \(J^{PC}\) quantum numbers of several representations of the Lorentz group, constrained to be C-even
With these \(J^{PC}\) numbers, we can use the method of Ji & Lebed to count form factors [76].
We consider a crossed channel matrix element \(\langle h\bar{h}|\hat{\mathcal {O}}|0\rangle \), and write out all the possible \(J^{PC}(L)\) quantum numbers that \(h\bar{h}\) can have up to \(J=2\). We then count the matches between these \(J^{PC}\) and those in our hypothetical rank-2 tensor. Such lists can be found in Table 2.
Table 2 Limited list of \(J^{PC}_L(S)\) quantum numbers possible for a \(h\bar{h}\) state. Only \(J^{PC}\) numbers present in Table 1 are included here. The lists for spin-0 and spin-1 are from [65, 77, 78], while the spin-1/2 list is from [76]
With these lists in hand, we can count the number of form factors that are contributing to the general rank-2 tensor by each of the representations it decomposes into. The numbers are given in Table 3.
Table 3 The number of form factors in a decomposition of a rank-2 tensor belonging to any of the representations listed in the table
Each of the rows in Table 3 tells us something meaningful about form factor counts.
The first row gives us trace terms, which can be thought of as "non-conserving" terms. There is one such term for both spin-0 and spin-1/2, and they are associated with the \(\bar{C}(t)\) GFFs. For spin-1, there are two non-conserving terms of this kind.
The second row tells us the number of GFFs appearing in the antisymmetric part of the EMT.
The third row tells us the number of GFFs in the decomposition of a traceless symmetric rank-2 tensor, and they correspond to the number of GFFs that are known to appear in the second Mellin moments of twist-2 vector GPDs.
Appendix B: Parametrization with Lorentz projectors
The matrix elements of the most general local, gauge-invariant EMT for a massive spin-1 target can be written as
$$\begin{aligned} \langle p^{\prime }, \lambda ^{\prime } \mid T_{\mu \nu }(0) \mid p, \lambda \rangle =\epsilon ^{*\alpha }(p^{\prime },\lambda ^{\prime }) T_{\mu \nu ,\alpha \beta }(P,\Delta )\epsilon ^{\beta }(p,\lambda ), \end{aligned}$$
(B1)
where the effective vertex \(T_{\mu \nu ,\alpha \beta }(P,\Delta )\) is a Lorentz tensor constructed out of the invariant tensors \(g_{\mu \nu }\) and \(\epsilon _{\mu \nu \alpha \beta }\), and the available four-vectors P and \(\Delta \). Because of the onshell relations \((p^{\prime }\epsilon ^{\prime *})=(p\epsilon )=0\), we choose to discard terms involving \(P_\alpha \) or \(P_\beta \) as they are not independent from those involving \(\Delta _\alpha \) or \(\Delta _\beta \). Finally, discrete symmetries impose further constraints on the effective vertex. Invariance under parity and time reversal implies that
$$\begin{aligned} T_{\mu \nu ,\alpha \beta }(P,\Delta )=T_{\bar{\mu }\bar{\nu },\bar{\alpha }\bar{\beta }} (\bar{P},\bar{\Delta })=T^*_{\bar{\mu }\bar{\nu },\bar{\alpha }\bar{\beta }}(\bar{P}, \bar{\Delta }), \end{aligned}$$
where a bar over a four-vector or a Lorentz index stands for the parity-transformed object \(\bar{a}^{\mu }=a^{\bar{\mu }}=(a^0,-\varvec{a})\). The hermiticity property leads to the last constraint
$$\begin{aligned} T_{\mu \nu ,\alpha \beta }(P,\Delta )=T^*_{\mu \nu ,\beta \alpha }(P,-\Delta ). \end{aligned}$$
It follows from these constraints that the effective vertex is real, and that its symmetric part under the exchange \(\alpha \leftrightarrow \beta \) involves only even powers of \(\Delta \), whereas the antisymmetric part involves only odd powers of \(\Delta \). Note also that the Levi-Civita tensor need not be considered since it is the only object with negative intrinsic parity and since the product \(\epsilon _{\mu \nu \alpha \beta }\epsilon _{\rho \sigma \tau \lambda }\) can be rewritten in terms of the metric only.
Interestingly, we can find a complete parametrization of the effective vertex in terms of the Lorentz projectors onto (0, 0), \((1,0)\oplus (0,1)\), and (1, 1) representations
$$\begin{aligned}&(I_{\mu \nu })_{\alpha \beta }=\frac{1}{4}g_{\mu \nu }g_{\alpha \beta }, \end{aligned}$$
(B4a)
$$\begin{aligned}&(D_{\mu \nu })_{\alpha \beta }=\frac{1}{2}(g_{\mu \alpha }g_{\nu \beta } -g_{\mu \beta }g_{\nu \alpha }), \end{aligned}$$
(B4b)
$$\begin{aligned}&(Q_{\mu \nu })_{\alpha \beta }=\frac{1}{2}(g_{\mu \alpha }g_{\nu \beta } +g_{\mu \beta }g_{\nu \alpha })-\frac{1}{4}g_{\mu \nu }g_{\alpha \beta }. \end{aligned}$$
(B4c)
We recognize in particular the generators of Lorentz transformations \((M^{\mu \nu })^{\alpha \beta }=2i(D^{\mu \nu })^{\alpha \beta }\) associated with the vector representation \((\frac{1}{2},\frac{1}{2})\). Taking into account all the constraints discussed above, we find that the most general effective vertex can be written as a linear combination of the following eleven independent Lorentz structures (omitting the \(\alpha ,\beta \) indices for convenience)
$$\begin{aligned}&I_{\mu \nu },P_\mu I_{\nu P},\Delta _\mu I_{\nu \Delta },\nonumber \\&P_\mu D_{\nu \Delta },P_\nu D_{\mu \Delta },\nonumber \\&Q_{\mu \nu },\Delta _\mu Q_{\nu \Delta },\Delta _\nu Q_{\mu \Delta }, g_{\mu \nu }Q_{\Delta \Delta },P_\mu P_\nu Q_{\Delta \Delta }, \Delta _\mu \Delta _\nu Q_{\Delta \Delta }.\nonumber \\ \end{aligned}$$
At first sight, it seems that these are at odds with the counting in Appendix A, since we found three structures involving the projector onto the (0, 0) representation and six involving the projector onto the (1, 1) representation. In fact, the counting in Appendix A corresponds to the projections applied on the pair of indices \(\mu \nu \). The structures corresponding to the (0, 0) part of the effective vertex are then the only two pure trace terms \(I_{\mu \nu }\) and \(g_{\mu \nu }Q_{\Delta \Delta }\). The structures corresponding to the (1, 1) part of the effective vertex are those obtained from removing the trace of the seven other symmetric terms \(P_\mu I_{\nu P}\), \(\Delta _\mu I_{\nu \Delta }\), \(Q_{\mu \nu }\), \(\Delta _\mu Q_{\nu \Delta }\), \(\Delta _\nu Q_{\mu \Delta }\), \(P_\mu P_\nu Q_{\Delta \Delta }\), and \(\Delta _\mu \Delta _\nu Q_{\Delta \Delta }\).
This technique has the advantage of showing the connection with lower spin targets. For spin-0 targets we have \(I_{\mu \nu }=\frac{1}{4} g_{\mu \nu }\), \(D_{\mu \nu }=0\) and \(Q_{\mu \nu }=0\), and for spin-1 / 2 targets we have \(I_{\mu \nu }=\frac{1}{4} g_{\mu \nu }\mathbb I\), \(D_{\mu \nu }=i\sigma _{\mu \nu }\) and \(Q_{\mu \nu }=0\).Footnote 10 According to (B5) this leads respectively to 3 and 5 GFFs, in agreement with Appendix A.
Appendix C: Polarization four-vector bilinears
Matrix elements of spin-1 targets can be written in terms of polarization four-vector bilinears contracted with some Lorentz tensors
$$\begin{aligned}&\langle p^{\prime }, \lambda ^{\prime } \mid O^{\mu _1\cdots \mu _n}(0) \mid p, \lambda \rangle \nonumber \\&\quad =\epsilon ^*_\alpha (p^{\prime },\lambda ^{\prime }) O^{\mu _1\cdots \mu _n,\alpha \beta }(P,\Delta )\epsilon _\beta (p,\lambda ). \end{aligned}$$
In the forward limit, the polarization four-vector bilinear reduces to the covariant density matrix which can be written as [54, 55]
$$\begin{aligned}&\epsilon _\beta (p,\lambda )\epsilon ^*_\alpha (p,\lambda ^{\prime }) =\rho _{\beta \alpha }(p)\nonumber \\&\quad =\left[ \frac{1}{3}\mathbb I +\frac{1}{2}\mathcal {S}^{\mu }(p)\Sigma _\mu +\mathcal {T}^{\mu \nu }(p)\Sigma _{\mu \nu }\right] _{\beta \alpha }, \end{aligned}$$
where the vector and tensor polarization operators read
$$\begin{aligned}&\Sigma _\mu =\frac{1}{2M}\epsilon _{\mu \rho \sigma p}M^{\rho \sigma }, \end{aligned}$$
$$\begin{aligned}&\Sigma _{\mu \nu }=-\Sigma _{\{\mu }\Sigma _{\nu \}}+\frac{1}{3} \left( g_{\mu \nu }-\frac{p_\mu p_\nu }{M^2}\right) \Sigma \cdot \Sigma , \end{aligned}$$
and are orthogonal to \(p^{\mu }\) as required by the onshell relation \((p\epsilon )=0\). Clearly, the covariant vector and tensor polarizations can be taken so as to satisfy the same properties as the associated operators
$$\begin{aligned}&p_\mu \mathcal {S}^{\mu }(p)=0, \end{aligned}$$
$$\begin{aligned}&\mathcal {T}^{\nu \mu }(p)= \mathcal {T}^{\mu \nu }(p), \end{aligned}$$
$$\begin{aligned}&p_\mu \mathcal {T}^{\mu \nu }(p)=0, \end{aligned}$$
$$\begin{aligned}&g_{\mu \nu }\mathcal {T}^{\mu \nu }(p)=0. \end{aligned}$$
They carry the dependence on the polarizations \(\lambda \) and \(\lambda '\), which has been omitted for convenience.
The generators of Lorentz transformations in the \((\frac{1}{2},\frac{1}{2})\) representation being given by
$$\begin{aligned} (M^{\rho \sigma })_{\alpha \beta }=i(\delta ^{\rho }_\alpha \delta ^{\sigma }_\beta -\delta ^{\rho }_\beta \delta ^{\sigma }_\alpha ), \end{aligned}$$
the unit, vector and tensor polarization operators take the explicit form
$$\begin{aligned}&(\mathbb I)_{\beta \alpha }=-P_{\beta \alpha }, \end{aligned}$$
$$\begin{aligned}&(\Sigma _\mu )_{\beta \alpha }=\frac{i}{M}\epsilon _{\mu \beta \alpha p}, \end{aligned}$$
$$\begin{aligned}&(\Sigma _{\mu \nu })_{\beta \alpha }=-P_{\mu \{\beta }P_{\alpha \}\nu } +\frac{1}{3}P_{\mu \nu }P_{\beta \alpha }, \end{aligned}$$
where \(P_{\mu \nu }\) is the projector onto the subspace orthogonal to \(p^{\mu }\) given by
$$\begin{aligned} P_{\mu \nu }=g_{\mu \nu }-\frac{p_\mu p_\nu }{M^2}. \end{aligned}$$
All these operators are orthogonal to each other, so that the normalization, covariant vector and tensor polarizations can simply be obtained as follows
$$\begin{aligned}&\text {Tr}[\rho (p)\mathbb I]=1, \end{aligned}$$
$$\begin{aligned}&\text {Tr}[\rho (p)\Sigma ^{\mu }]=-\mathcal {S}^{\mu }(p), \end{aligned}$$
$$\begin{aligned}&\text {Tr}[\rho (p)\Sigma ^{\mu \nu }]=\mathcal {T}^{\mu \nu }(p). \end{aligned}$$
Denoting the rest-frame four-momentum as \(k^{\mu }=Mg^{\mu 0}\), the general polarization four-vector bilinear can be written as
$$\begin{aligned} \epsilon ^{\beta }(p,\lambda )\epsilon ^{*\alpha }(p^{\prime },\lambda ^{\prime }) =(\Lambda _{\mathrm{can}})^{\beta }_{\phantom {\beta }\rho } \,(\Lambda ^{\prime }_{\mathrm{can}})^{\alpha }_{\phantom {\alpha }\sigma } \,\rho ^{\rho \sigma }(k), \end{aligned}$$
where the standard canonical rotationless Lorentz boost tensor reads
$$\begin{aligned}&(\Lambda _{\mathrm{can}})^{\mu }_{\phantom {\mu }\nu }\nonumber \\&\quad =\begin{pmatrix} \frac{p^0}{M}&{}\quad \frac{p^1}{M}&{}\quad \frac{p^2}{M}&{}\quad \frac{p^3}{M}\\ \frac{p^1}{M}&{}\quad 1{+}\frac{(p^1)^2}{M(p^0+M)}&{}\quad \frac{p^1p^2}{M(p^0+M)}&{}\quad \frac{p^1p^3}{M(p^0+M)}\\ \frac{p^2}{M}&{}\quad \frac{p^1p^2}{M(p^0+M)}&{}\quad 1{+}\frac{(p^2)^2}{M(p^0+M)}&{}\quad \frac{p^2p^3}{M(p^0+M)}\\ \frac{p^3}{M}&{}\quad \frac{p^1p^3}{M(p^0+M)}&{}\quad \frac{p^2p^3}{M(p^0+M)}&{}\quad 1+\frac{(p^3)^2}{M(p^0+M)}, \end{pmatrix}\nonumber \\ \end{aligned}$$
or in a more covariant form [79, 80]
$$\begin{aligned} (\Lambda _{\mathrm{can}})^{\mu }_{\phantom {\mu }\nu }=\delta ^{\mu }_\nu -\frac{(p+k)^{\mu }(p+k)_\nu }{k\cdot (p+k)}+\frac{2p^{\mu } k_\nu }{M^2}. \end{aligned}$$
Denoting for convenience the rest-frame covariant vector and tensor polarizations as
$$\begin{aligned} \mathcal {S}^{\mu }\equiv \mathcal {S}^{\mu }(k),\qquad \mathcal {T}^{\mu \nu } \equiv \mathcal {T}^{\mu \nu }(k), \end{aligned}$$
we find that the general polarization four-vector bilinear can be expressed as
$$\begin{aligned}&\epsilon ^{\beta }(p,\lambda )\epsilon ^{*\alpha }(p^{\prime },\lambda ^{\prime })\nonumber \\&=\frac{1}{3}\left[ -g^{\beta \alpha }+\frac{p^{\beta } p^{\prime \alpha }}{M^2} -\frac{(p+k)^{\beta }\Delta ^{\alpha }}{k\cdot (p+k)}\right. \nonumber \\&\qquad \left. +\frac{\Delta ^{\beta }(p^{\prime }+k)^{\alpha }}{k\cdot (p^{\prime }+k)} +\frac{\Delta ^2}{2}\frac{(p+k)^{\beta }(p^{\prime }+k)^{\alpha }}{[k\cdot (p+k)][k\cdot (p^{\prime }+k)]}\right] \nonumber \\&\qquad +\frac{i}{2M}\left[ \epsilon ^{\beta \alpha \mathcal {S} k} -\epsilon ^{\beta p^{\prime }\mathcal {S} k}\frac{(p^{\prime }+k)^{\alpha }}{k\cdot (p^{\prime }+k)}\right. \nonumber \\&\qquad \left. -\frac{(p+k)^{\beta }}{k\cdot (p+k)}\epsilon ^{p\alpha \mathcal {S} k} +\epsilon ^{pp^{\prime }\mathcal {S}k} \frac{(p+k)^{\beta } (p^{\prime }+k)^{\alpha }}{[k\cdot (p+k)][k\cdot (p^{\prime }+k)]}\right] \nonumber \\&\qquad -\mathcal {T}^{\beta \alpha }+\mathcal {T}^{\beta p^{\prime }} \frac{(p^{\prime }+k)^{\alpha }}{k\cdot (p^{\prime }+k)} +\frac{(p+k)^{\beta }}{k\cdot (p+k)}\mathcal {T}^{p\alpha }\nonumber \\&\qquad -\mathcal {T}^{pp^{\prime }}\frac{(p+k)^{\beta }(p^{\prime }+k)^{\alpha }}{[k\cdot (p+k)][k\cdot (p^{\prime }+k)]}. \end{aligned}$$
In the forward limit, we recover the standard expression for the covariant density matrix
$$\begin{aligned} \rho _{\beta \alpha }(p)=-\frac{1}{3}P_{\beta \alpha } +\frac{i}{2M}\epsilon _{\beta \alpha \rho \sigma } \mathcal {S}^{\rho }(p)p^{\sigma }-\mathcal {T}_{\beta \alpha }(p), \end{aligned}$$
$$\begin{aligned} \mathcal {S}^{\mu }(p)=(\Lambda _{\mathrm{can}})^{\mu }_{\phantom {\mu }\rho } \mathcal {S}^{\rho },\qquad \mathcal {T}^{\mu \nu }(p) =(\Lambda _{\mathrm{can}})^{\mu }_{\phantom {\mu }\rho } (\Lambda _{\mathrm{can}})^{\nu }_{\phantom {\nu }\sigma }\mathcal {T}^{\rho \sigma }. \end{aligned}$$
As noticed in the case of Dirac bilinears [81], in the forward limit \(\Delta =0\) the dependence on the rest-frame four-momentum \(k^{\mu }\) can be absorbed into the p-dependent covariant polarization. However, in the general off-forward case \(\Delta \ne 0\) this is not possible anymore and the k-dependence remains explicit. This dependence is unavoidable in a relativistic theory and comes from the fact that, because of the non-commutativity of boosts, canonical polarization of a massive particle has to be defined from a boost of the rest-frame polarization to the frame of interest [82].
In the Breit frame \(\vec {P}=\vec {0}\), we find that the general polarization four-vector bilinear reduces to
$$\begin{aligned}&\epsilon ^{\beta }(p,\lambda )\epsilon ^{*\alpha } (p^{\prime },\lambda ^{\prime })\nonumber \\&\quad =\frac{1}{3}\left[ -g^{\beta \alpha } +\left( 1+\frac{\Delta ^2}{4M^2}\right) \frac{k^{\beta } k^{\alpha }}{M^2}-\frac{k^{[\beta }\Delta ^{\alpha ]}}{M^2}+\frac{\Delta ^{\beta }\Delta ^{\alpha }}{4M^2}\right] \nonumber \\&\qquad +\frac{i}{2M}\left[ \epsilon ^{\beta \alpha \mathcal {S}k} -\frac{k^{\{\beta }\epsilon ^{\alpha \}\Delta \mathcal {S} k}}{M^2} +\frac{\Delta ^{[\beta }\epsilon ^{\alpha ] \Delta \mathcal {S} k}}{4M^2}\right] \nonumber \\&\qquad {-}\mathcal {T}^{\beta \alpha }{-}\frac{k^{[\beta } \mathcal {T}^{\alpha ]\Delta }}{M^2}{+}\frac{\Delta ^{\{\beta } \mathcal {T}^{\alpha \}\Delta }}{4M^2}+\frac{\mathcal {T}^{\Delta \Delta }}{4M^2} \frac{k^{\beta } k^{\alpha }}{M^2}+\mathcal {O}(\Delta ^3).\nonumber \\ \end{aligned}$$
Appendix D: Covariant parametrization of Generalized Parton Distributions
Just like for local operators, hadronic matrix elements of non-local operators can be parametrized in a covariant form if one includes in the list of available four-vectors the direction of non-locality, namely the lightlike four-vector n in the case of parton distributions like GPDs. The most general parametrization of the quark GPD correlator (39a39b), that respects the constraints imposed by hermiticity, parity and time-reversal, reads
$$\begin{aligned} V_{\lambda \lambda ^{\prime }}^{\mu } = \frac{1}{(Pn)}\epsilon ^*_\alpha (p^{\prime },\lambda ^{\prime }) \mathcal {V}^{\mu ,\alpha \beta }(P,\Delta ,n)\,\epsilon _\beta (p,\lambda ) \end{aligned}$$
(D1)
$$\begin{aligned}&\mathcal {V}^{\mu ,\alpha \beta }(P,\Delta ,n)\nonumber \\&\quad = g^{\alpha \beta }\left[ P^{\mu } F_1+\Delta ^{\mu } F_2+\frac{M^2n^{\mu }}{(Pn)}\,F_3 \right] \nonumber \\&\qquad +\frac{\Delta ^{\alpha }\Delta ^{\beta }}{M^2} \left[ P^{\mu } F_4+\Delta ^{\mu } F_5 +\frac{M^2n^{\mu }}{(Pn)}\,F_6 \right] \nonumber \\&\qquad +\frac{M^2n^{\alpha }n^{\beta }}{(Pn)^2} \left[ P^{\mu } F_7+\Delta ^{\mu } F_8+\frac{M^2n^{\mu }}{(Pn)}\,F_9 \right] \nonumber \\&\qquad +\frac{n^{\{\alpha } \Delta ^{\beta \}}}{(Pn)} \left[ P^{\mu } F_{10}+\Delta ^{\mu } F_{11}+\frac{M^2n^{\mu }}{(Pn)}\,F_{12} \right] \nonumber \\&\qquad +\frac{n^{[\alpha } \Delta ^{\beta ]}}{(Pn)} \left[ P^{\mu } F_{13}+\Delta ^{\mu } F_{14}+\frac{M^2n^{\mu }}{(Pn)}\,F_{15} \right] \nonumber \\&\qquad +g^{\mu \{\alpha } \left[ \Delta ^{\beta \}} F_{16}+\frac{M^2n^{\beta \}}}{(Pn)}\,F_{17} \right] \nonumber \\&\qquad +g^{\mu [\alpha } \left[ \Delta ^{\beta ]} F_{18}+\frac{M^2n^{\beta ]}}{(Pn)}\,F_{19} \right] , \end{aligned}$$
and \(F_i = F_i(x,\xi ,t)\). Since \(\mathcal {V}^{\mu , \alpha \beta }(P,\Delta ,n)\) must be invariant under a rescaling of the lightlike direction \(n\mapsto \alpha n\), factors of M / (Pn) appear whenever necessary. Additional factors of the hadron mass M have also been included to keep GPDs dimensionless. The relation between \(F_i\) and the standard GPD basis \(H_i,G_i\) of (40) can easily be obtained by projection onto the twist-2 and 3 parts.
Cosyn, W., Cotogno, S., Freese, A. et al. The energy-momentum tensor of spin-1 hadrons: formalism. Eur. Phys. J. C 79, 476 (2019). https://doi.org/10.1140/epjc/s10052-019-6981-3
DOI: https://doi.org/10.1140/epjc/s10052-019-6981-3