**Oleg Teryaev 1,2**

<sup>1</sup> Joint Institute for Nuclear Research, 141980 Dubna, Russia; teryaev@jinr.ru

<sup>2</sup> Institute of Engineering and Physics, Dubna International University, 141980 Dubna, Russia

Received: 15 June 2020; Accepted: 18 August 2020; Published: 24 August 2020

**Abstract:** The energy-momentum relocalization in classical and quantum theory is addressed with specific impact on non-perturbative QCD and hadronic structure. The relocalization is manifested in the existence of canonical and symmetric (Belinfante and Hilbert) energy momentum tensors (EMT). The latter describes the interactions of hadrons with classical gravity and inertia. Canonical EMT, in turn, is naturally emerging due to the translation invariance symmetry and appears when spin structure of hadrons is considered. Its relation to symmetric Hilbert and Belinfante EMTs requires the possibility to neglect the contribution of boundary terms for the classical fields. For the case of quantum fields this property corresponds to the absence of zero-momentum poles of matrix element of the axial current dual to the spin density. This property is satisfied for quarks manifesting the symmetry counterpart of *UA*(1) problem and may be violated for gluons due to QCD ghost pole.

**Keywords:** gravity; relocalization; topology; boundary; poles

## **1. Introduction**

The space-time symmetry related to the energy-momentum and angular momentum conservation is manifested in field theory as the appearance of energy-momentum and spin currents (see [1] and Ref. therein). Their definition is not unique. The structure of Lagrangian immediately defines, after the application of the Noether theorem, the canonical densities. Passing to the quantum operators and their matrix elements one may analyse how the fundamental fields are manifested in the spin structure of elementary particles.

From the other side, the interaction of particles with gravity involves the symmetruc Hilbert tensor resulting from the variation with respect to metric and the symmetry property naturally emerges after the absorption of spin density to the orbital one using the Belinfante procedure [1].

The interplay of both forms of Energy-Momentum Tensor (EMT) is especially important for hadrons, as in the absence of mathematically rigorous confinement theory their spin structure is an important problem of non-perturbative Quantum Chromodynamics (QCD). The same non-perturbative effects are responsible for the most of the visible mass of the Universe, and, therefore, for its gravitational interaction.

These interactions of hadrons with gravity (and inertia, due to equivalence principle) are encoded in their gravitational formfactors [2–6] (see also [7] and Ref. therein). They define the macrospic properties of all objects, and, as it appeared more recently, the responce of hadrons to fastest ever rotation and acceleration emerging in heavy-ion collisions (see [8] and Ref. therein). Indeed, the angular velocity of quark-gluon matter in the non-central heavy-ion collisions corresponds to the change of the velocity of the order of speed of light *c* at the distance of order of Compton wavelength *lC*, *ω* ∼ *c*/*lC*, which is some 25 orders of magnitude larger than angular velocity of Earth rotation. By coincidence, the acceleration *<sup>a</sup>* of this matter which is of the order of *<sup>a</sup>* <sup>∼</sup> *<sup>c</sup>*2/*lC* is larger than Earth's gravity *<sup>g</sup>* by almost the same factor.

One may wonder, why the highly non-inertial frame formed by quark-gluon matter in heavy-ion collisions can have any impact on the observables measured by the detector located in the laboratory frame. The non-inertial matter will play a role if its interaction with hadrons may be considered as a quantum measurement which is certainly true if particle spin (essentially quantum object!) or Hawking–Unruh radiation [8] are considered. Let us also note here that the main outcome of equivalence principle (EP) for spin motion in the gravitational field, the equality of classical and quantum rotators (orbital and spin angular momenta) precession frequencies (see [7] and Ref. therein), becomes trivial for the rotating frame (like Earth) if spin is considered just as some vector remaining constant [9] in the inertial frame and rotating in the frame of the Earth like Foucault pendulum. The non-trivial meaning would emerge if the quantum measurement of the spin in the rotating frame is considered and its similarity to pendulum and orbital angular momentum (AM) is the manifestation of EP.

The interaction with gravity, as it was already mentioned, is described by the symmetric Hilbert EMT representing (like Belinfante and any symmetric EMT) the angular momentum (AM) as the orbital one. This form of AM allows one to derive the EP as low-energy theorem (see [7] and Ref. therein) by making use of momentum and angular momentum conservation. Like in QED, the global symmetry puts the restrictions for the interaction (defined by local one) for small momenta. In distinction from QED, where only terms of zero order in momenta are fixed, while the linear ones (momenta) are dynamical, in gravity, because of more complicated gauge group, the momenta are also fixed and anomalous gravitomagnetic moment is absent, which is another formulation [2] of EP.

At the same time, in the analysis of hadronic spin structure the canonical expressions naturally appear. The spin of fundamental fields plays the special role in the physical interpretation of QCD hadronic structure. The interplay of various forms of EMT and AM are discussed in detail in the problems of hadronic structure [10] and heavy-ion collisions [11].

The procedure of relocalization [1] changing the local quantities but preserving the conserved (angular) momenta, requires the possibility to discard the surface terms, which is usually assumed. At the same time, their consideration when generalized to the case of quantum operators was earlier found [12–14] to lead to some non-trivial constraints for matrix elements of axial currents. Here we develop these ideas and put them into modern context. As a result, we find the constraints for the zero-mass poles in matrix elements of singlet axial current (dual to quark spin desnity) leading to an amazing interplay between general symmetry properties of relocalization and very specific topologocal QCD dynamics.

#### **2. Boundary Terms in Coordinate and Momentum Space**

Let us start with the following expression for the quark–gluon angular momentum density

$$M^{\mu,\nu\rho} = \frac{1}{2} \epsilon^{\mu\nu\rho\sigma} J\_{\mathfrak{F}\rho\sigma} + \mathbf{x}^{\nu} T^{\mu\rho} - \mathbf{x}^{\rho} T^{\mu\nu}.\tag{1}$$

The first term in the right hand side (r.h.s.) is just the canonical quark spin tensor dual to singlet axial current. Note that the energy–momentum tensor here accumulates also the quark orbital momentum as well as the total gluon angular momentum. We may proceed further along this way and express the quark spin in the orbital form with the simultaneous change of the energy–momentum tensor to the one suggested by Belinfante long ago

$$M\_B^{\mu,\nu\rho} = \mathbf{x}^{\nu} T\_B^{\mu\rho} - \mathbf{x}^{\rho} T\_B^{\mu\nu}. \tag{2}$$

As the conservation of the angular momentum

$$
\partial\_{\mu}M\_{B}^{\mu,\nu\rho} = 0\tag{3}
$$

immediately leads to the symmetry of *Tμρ* (so that symmetric Hilbert tensor may also be considered in such a role), the latter implies that

$$
\epsilon\_{\mu\nu\rho\mathfrak{u}} M\_B^{\mu,\nu\rho} = 0.\tag{4}
$$

One might conclude that the totally antisymmetric quark spin tensor is somehow cancelled and does not contribute to the total angular momentum [15]. This is also the manifestation of the general belief that axial current and angular momentum represent the different aspects of spin structure. Still, it appears possible to extract some quantitative information about their interplay.

The Belinfante (and Hilbert) EMT lead to the same AM as canonical one so that:

$$\int d^3 \mathbf{x} M\_{\mathbf{B}}^{0,\nu\rho} = \int d^3 \mathbf{x} M^{0,\nu\rho}. \tag{5}$$

We assume (which was done only tacitly in [12]) that also the stronger condition is valid

$$\int d^3 \mathbf{x} M\_{\mathbf{B}}^{\mu,\nu\rho} = \int d^3 \mathbf{x} M^{\mu,\nu\rho} \, , \tag{6}$$

so that

$$
\epsilon\_{\mu\nu\rho\alpha} \int d^3 \mathbf{x} \mathcal{M}^{\mu\nu\rho} = 0. \tag{7}
$$

Substituting here the definition (1) one get

$$\int d^3 \mathbf{x} \left( \Im f\_5^a(\mathbf{x}) + 2 \epsilon^{\mu \nu \rho a} \mathbf{x}\_{\nu} T\_{A, \mu a}(\mathbf{x}) \right) = 0,\tag{8}$$

where *Tμα <sup>A</sup>* = (*Tμα* <sup>−</sup> *<sup>T</sup>αμ*)/2 is the antisymmetric part of energy momentum tensor responsible for the separate non-conservation of orbital and spin AM. The conservation of total AM results in the relation

$$\frac{1}{4} \epsilon^{\mu\nu\rho\sigma} \partial\_{\mu} f^{5}\_{\sigma} = T^{\rho\nu}\_{A} \,. \tag{9}$$

Note that we consider the weak gravitational fields and the derivatives, as well as quantum states in what follows, correspond to flat space. The consideration of strong fields may be achieved by applying the Dirac equation in curved space (see [7] and Ref. therein).

Set of Equations (8) and (9) allows one to exclude either spin (*J*5) or orbital (*TA*) AM. The latter is easier, as one can use the local Equation (9). Furthermore, axial current operator is related to many observables.

Excluding EMT antisymmetric part by making use of the conservation of AM density (9), one get in the case of classical fields:

$$\left(g\_{\rho\nu}g\_{a\mu} - g\_{\rho\mu}g\_{a\nu}\right)\int d^3x\partial^\rho(f\_5^a\mathbf{x}^\nu) = 0. \tag{10}$$

This is in fact the way to represent the surface terms whose neglect is necessary to apply the Belinfante procedure.

Passing to the most interesting case of quantum operators one should switch the Equations (8) and (after incorporating the AM conservation at operator level) (10) between particle (nucleon) states with the momenta *P* and *P* + *q*.

$$\int (\mathcal{g}\_{\rho\nu}\mathcal{g}\_{a\mu} - \mathcal{g}\_{\rho\mu}\mathcal{g}\_{a\nu}) \int d^3 \mathbf{x} \langle P|\partial^{\rho}(f^{\mu}\_5 \mathbf{x}^{\nu})|P+q\rangle = 0. \tag{11}$$

Expressing the local operator *J<sup>α</sup>* <sup>5</sup> (*x*) = *exp*(*iPx*<sup>ˆ</sup> )*J<sup>α</sup>* <sup>5</sup> (0)*exp*(−*iPx*<sup>ˆ</sup> ) by action of shift operator *exp*(*iPx*ˆ ) allows one to perform the integration resulting in appearance of *δ*(*q*). Now *<sup>∂</sup> <sup>∂</sup>x<sup>μ</sup>* is substituted by <sup>−</sup>*iq<sup>μ</sup>* and *<sup>x</sup><sup>μ</sup>* by *<sup>i</sup> <sup>∂</sup> <sup>q</sup><sup>μ</sup>* acting on that *<sup>δ</sup>*3(*q*). The latter is, by definition, equal, up to a sign, to the derivative acting on the matrix element. As a result, one obtains the following constraint:

$$q^{\mu} \frac{\partial}{\partial q^{a}}|\_{0} \langle P|J\_{5}^{a}|P+q\rangle = q^{a} \frac{\partial}{\partial q^{a}}|\_{0} \langle P|J\_{5}^{\mu}|P+q\rangle. \tag{12}$$

It is a quantum counterpart of (10) and it is natural that surface terms in coordinate space correspond to zero momenta. To make it more clear, let us multiply both sides by *qμ*:

$$q^2 \frac{\partial}{\partial q^\alpha} \langle P | J\_5^a | P + q \rangle = (q^\beta \frac{\partial}{\partial q^\beta} - 1) q\_\gamma \langle P | J\_5^\gamma | P + q \rangle. \tag{13}$$

This equality is obviously valid up to the second and higher powers of *q*. Note that the differential operator in the r.h.s. subtracts the terms linear in *q* from the divergence matrix element proportional to *sq* for the pure kinematical reasons.

What can be dangerous is the pole for *<sup>q</sup>*<sup>2</sup> <sup>→</sup> 0 which naturally appears for anomalous axial current already in perturbation theory for massless fermions [16]. For massless quarks, due to t'Hooft consistency principle, in the case of non-singlet currents these poles correspond to the exchange of the massless Goldstone mesons.

The exception is provided by singlet channel where *η* - remains massive manifesting the famous *UA*(1) problem and the correspondent pole is absent (see, e.g., [17] and Ref. therein):

$$
\langle P, S \vert J \mathfrak{z}\_{\varbeta}(0) \vert P + q, S \rangle = 2MS\_{\varbeta}G\_1 + q\_{\varbeta}(Sq)G\_2,\tag{14}
$$

$$q^2 G\_2|\_0 = 0.\tag{15}$$

The *G*<sup>2</sup> pole term, if present, would provide a contribution linear in *q* to the l.h.s.

Therefore, solution of *UA*(1) problem provides simulatenously the necessary dynamical mechanism for relocalization of massless quarks spin. This, in turn, leads to relation of conservation laws and canonical EMT with Belinfante and Hilbert EMT, supporting the emergence of equivalence principle as low-energy theorem [7].

#### **3. Problems with Relocalization for Gluons**

The situation is changed in the case of gluons. The relevant matrix element of topological current

$$
\langle P, S \vert K\_{\mu}^{\tilde{\sigma}}(0) \vert P + q, S \rangle = 2MS\_{\mu}G\_{1}(q^{2}) + q\_{\mu}(Sq)G\_{2}(q^{2}), \tag{16}
$$

$$q^2 G\_2(q^2)|\_{0} \neq 0,\tag{17}$$

contains the contribution *G*˜2 of the relevant Kogut–Susskind ghost (or instanton [18]) pole [19] which is fully responsible [20] for the value of the forward matrix element of anomaly-free quark gluon current *J μ* <sup>5</sup> <sup>−</sup> *<sup>K</sup>μ*.

The consideration of topological current as dual to spin is naturally supported by the studies of bosonic anomalies in gravitational field [21] which may be relevant also for consideration of rotating quark-gluon matter in heavy-ion collisions [22]. Therefore one has a contradiction between kinematics of Relocalization Invariance (RI) requiring the absence of surface terms (corresponding to zero-mass poles of matrix elements) and instanton-type dynamics requiring their presence. The possible outcomes are the following

(i) If RI is indeed violated the coupling of nucleons to gravity (described by the formfactors of Belinfante EMT [13,23]) may be unconstrained by the form of conservation laws in terms of canonical EMT. In extreme case, assuming that just the canoncical form is related to translational invariance, this might result in the violation of Equivalence Principle for nucleons at several percent level which may be tested experimentally and is probably already excluded by the data. (ii) One may assume "Hadronic censorship" leading to the absence of the ghost pole: in this case the matrrix element

$$> 0$$

in the chiral limit. Bearing in mind the smallness of gluon spin one should mostly attribute the quark spin to the (predominantly strange) quark mass. This may explain the relative smallness of quark spin ("Spin Crisis") and may be checked, say, by lattice calculations of pseudoscalar quark densities.

(iii) The simplest solution would be the impossibility to separate spin and orbital momenta of gluons in the meaningful way.

#### **4. Discussion**

We found the relation between general space-time symmetry responsible for interactions of hadrons with gravity and the specific QCD dynamics. As a result, the quark spin relocalization is supported by solution of *UA*(1) problem and the same non-perturbative dynamics may spoil the extraction of totally antisymmetric gluon spin density.

The future studies, besides the exploration of mentioned in the previous section alternatives may include following developments:


**Funding:** This research was funded by RFBR grant 18-02-01107.

**Acknowledgments:** I am intebted to Yu.N. Obukhov, G.Yu. Prokhorov, A.Ya. Silenko and V.I. Zakharov for useful discussions.

**Conflicts of Interest:** The author declares no conflict of interest.

#### **References**


c 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Symmetry* Editorial Office E-mail: symmetry@mdpi.com www.mdpi.com/journal/symmetry

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18