On the hadron mass decomposition

We argue that the standard decompositions of the hadron mass overlook pressure effects, and hence should be interpreted with great care. Based on the semiclassical picture, we propose a new decomposition that properly accounts for these pressure effects. Because of Lorentz covariance, we stress that the hadron mass decomposition automatically comes along with a stability constraint, which we discuss for the first time. We show also that if a hadron is seen as made of quarks and gluons, one cannot decompose its mass into more than two contributions without running into trouble with the consistency of the physical interpretation. In particular, the so-called quark mass and trace anomaly contributions appear to be purely conventional. Based on the current phenomenological values, we find that in average quarks exert a repulsive force inside nucleons, balanced exactly by the gluon attractive force.


Introduction
According to the Standard Model of Particle Physics, the masses of almost all known elementary particles are generated through the Brout-Englert-Higgs (BEH) mechanism. The current light quark masses obtained from such a mechanism correspond however to only about 1% of the nucleon mass. Therefore, the mass of the ordinary matter around us essentially finds its origin in the strong interactions which confine quarks and gluons inside hadrons, and not the BEH mechanism [1,2].
Understanding the origin of the hadron mass in Quantum Chromodynamics (QCD) represents a formidable challenge owing to the relativistic, quantum and non-perturbative nature of the problem. The QCD lagrangian in the chiral limit appears to be scale invariant at the classical level. This implies in particular that all the hadron masses should vanish in that limit. Scale invariance is however an anomalous symmetry, in the sense that it appears to be broken at the quantum level by radiative corrections. This gives rise, via dimensional transmutation, to a dimensionful parameter Λ QCD ≈ 0.2 GeV in the theory [3,4], and therefore to a non-trivial spectrum.
Lattice QCD calculation of hadron masses obtained from the analysis of correlations functions in Euclidean time is in remarkable agreement with the experimental spectrum [5,6]. Unfortunately, this method gives little insight on how these masses arise from quark and gluon contributions.
In order to address the question of the origin of the hadron mass, one should rather start from the energy-momentum tensor (EMT) of QCD. At the classical level, it reads where ψ is the quark field, G aµν is the gluon field strength, 2igA aµ t a is the symmetric non-abelian covariant derivative, and η µν is the Minkowski metric. Note that, as argued e.g. in [7,8], there is no need in Particle Physics to require the EMT to be symmetric in its Lorentz indices. This aspect will however not affect the purpose of the present paper and will therefore not be discussed further. An important property of the total EMT is that it is conserved ∂ µ T µν = 0, which is a consequence of the invariance of the theory under space-time translations. At the quantum level, it does not therefore require any global renormalization [9]. The trace of the renormalized EMT is given by [10,11,12,13,14,15] T µ µ = β(g) 2g G 2 + (1 + γ m ) ψmψ, where m is the quark mass matrix and γ m is its anomalous dimension. The first term is known as the trace anomaly, which is a pure quantum effect as indicated by the β function factor. In the chiral limit m → 0, this term prevents the trace of the EMT to vanish, and hence prevents QCD to become a scale-invariant theory with trivial spectrum. Mass decompositions can be obtained from the expectation value of the EMT. In the next Section, we briefly present the two standard decompositions used in the literature, and point out a couple of issues regarding their physical interpretation. To the best of our knowledge, some of these issues have not been addressed before, in particular regarding pressure effects. Based on the semiclassical picture, we propose in Section 3 a new decomposition free of these issues. Then we discuss in Section 4 the new picture in more detail and we estimate the various contributions using our current phenomenological knowledge. Finally, we summarize our results in Section 5.
where the state is normalized according to The hadron mass M can then be defined via the trace of the EMT [17,18,16,20,19,21,22,23] This expression is quite appealing owing to its manifest Poincaré invariance. In particular, it relies on the identity P 2 = M 2 , which defines mass as the square root of one of the Casimir operators associated with the Poincaré group. Using the explicit expression (2) for the trace of the EMT operator it is tempting to conclude that most of the light hadron masses comes from gluons, based on the fact that the second term on the right hand side is known to be a rather small contribution (except for the pion) [17,24,25,26,27]. Note that this partonic picture applies whenever the renormalization scale is much larger than 2 GeV. For a renormalization scale below 0.5 GeV, the picture may look quite different and closer to a non-relativistic picture where the (effective) constituents are quite massive and the (residual) binding energy is small and negative [25,26].
The physical interpretation of the origin of hadron mass based on Eq. (5) suffers however from a couple of caveats. Indeed, it seems a bit odd from a physical point of view that the definition of mass involves the trace of the EMT only at a single point. Moreover, the precise form of Eq. (3), and hence of Eq. (5), depends on how hadron states are normalized. A proper physical interpretation should in principle not depend on our conventions for the normalization of states. It is therefore more natural to define mass via the normalized expectation value of some spatially extended operator It is understood that such an expectation value contains implicitly the subtraction of the corresponding vacuum matrix element. Note in particular that the precise form of Eq. (3) was fixed by the requirement that T 0µ = P µ . It follows that the proper expectation value of the trace of the EMT reads independently of how states are normalized. The presence of the P 0 factor indicates that this expectation value is frame dependent. In our opinion, frame dependence is not really a problem since physical quantities often find a simple interpretation only in some particular set of reference frames. For instance, using the light-front form of dynamics [28], Eq. (8) becomes T µ µ = M 2 /P + which can be interpreted in the symmetric frame (defined by P ⊥ = 0 ⊥ ) as twice the light-front energy T µ µ | P ⊥ =0 ⊥ = 2P − . Frame dependence 1 seems to be a generic feature associated with many decompositions of a system. If one insists on manifest Lorentz invariance 2 , it turns out that Lorentz non-invariant quantities defined in a given frame can be formally expressed as Lorentz scalars [7,29]. By construction, the latter will give the same result in any frame, and will coincide with the original noninvariant quantities in the given frame. Typical examples are the proper time and length which are Lorentz-invariant quantities giving respectively the time and length measured in the rest frame of the system. In other words, it is in principle always possible to express quantities in a Lorentz-invariant form, but the price to pay is that the physical interpretation necessarily singles out a privileged frame. In the present case, the frame dependence of Eq. (8) originates from the fact that we are considering the operator d 3 r T µ µ (r), which is not Lorentz invariant unlike T µ µ (0). It is however in line with the standard expectation that mass effects in hadronic matrix elements should become negligible in the ultra-relativistic regime P 0 ≫ M . If one insists on preserving manifest Lorentz invariance, one can alternatively consider the integral over the proper volume P 0 M d 3 r T µ µ (r). Frame dependence will disappear from the formal expressions but it will remain in their physical interpretation.
Beside the question of frame dependence, a more important point we would like to stress is that the physical interpretation should be directly based on the quantum operator. Providing an interpretation at the level of the matrix element may be misleading. For example, in Eq. (5) the hadron mass is connected to the trace of the EMT at the level of matrix elements, but the connection remains obscure at the operator level, especially when the operator is decomposed into several contributions. Indeed, the trace T µ µ involves beside energy T 0 0 the spatial components T i i usually associated with normal stresses. It seems also counterintuitive that the contribution from the quark mass to the hadron mass is not proportional to m, but rather to √ m. In fact, there exist (infinitely) many operators whose forward matrix elements are proportional to some power of the hadron mass. Following the same logic as for the trace operator, one could then in principle argue that these operators can be used to provide alternative definitions (and hence alternative decompositions) of the hadron mass, raising the question of deciding which one is the "correct" one. This happens simply because the hadron mass gives the natural scale of the matrix elements. The only way out is to use a quantum operator with a clear connection to the concept of mass.
In view of all these caveats, we feel that claims about the physical origin of hadron mass based on Eq. (5) are somewhat misleading. Note also that the above remarks can similarly be applied to the Gell-Mann-Oakes-Renner formula [30].

Ji's decomposition
About 20 years ago, Ji proposed a decomposition of the hadron mass analogous to the virial theorem for a harmonic oscillator and the hydrogen atom [31,32], see also [33]. As a first step, he decomposed the QCD EMT into traceless and trace parts T µν =T µν +T µν (9) withT µν = 1 4 η µν T α α . This choice is of course not unique 3 , but it is well suited for treating the trace anomaly contribution. Using the relativistic normalization (4) for the states, the corresponding forward matrix elements read 3 The general trace part is defined asT µν c =T µν + cT µν with c an arbitrary constant.
and are renormalization scale independent, as required by Lorentz symmetry.
In the Hamiltonian formalism, the mass of a system is defined as the total energy in the rest frame. It then follows that In other words, the trace of the EMT appears to contribute only to 25% of the hadron mass. Ji proceeded with a decomposition of the traceless part into quark and gluon contributions, and of the trace part into quark mass and trace anomaly contributions The corresponding individual forward matrix elements are not protected by Lorentz symmetry and so the coefficients a(µ 2 ) and b(µ 2 ) depend generally on the renormalization scale µ. It follows from Eq. (13) that the Hamiltonian operator can naturally be decomposed as where Owing to the QCD equations of motion, the quark contribution can be put in the form 4 where D = ∇ + igA. Rearranging a bit the quark mass term between H ′ q and with Ji arrived at the following decomposition of the hadron mass in the rest frame where In particular, using the parametrization in Eqs. (14)- (17), the individual contributions can be expressed as In the chiral limit, one has b = 0 and the gluon energy arising from the trace anomaly corresponds precisely to the vacuum energy introduced phenomenologically in the MIT bag model [34]. This decomposition is quite popular in hadronic physics [2] and motivated many lattice QCD calculations, see e.g. [35,36,37] for recent related works.
Ji's decomposition is sometimes criticized because it is performed in a specific frame and applies to massive states only [25,26]. As we already argued in the previous section, the frame dependence is not really a problem but a general feature. Note also that Ji's decomposition can formally be put in a covariant form by considering the Lorentz-invariant quantity T 0µ u µ , where u µ ≡ P µ /M is interpreted as the hadron four-velocity. Although this form is frame independent, its physical interpretation becomes simple only in the rest frame, where u µ = (1, 0). For a massless state, since there is no rest frame, one can consider instead the energy decomposition in any frame. In this case, there is no contribution from the trace part since P 2 = 0 and so only the coefficient a(µ 2 ) is needed 5 .
The actual problem with Ji's decomposition is that the separation of the EMT into traceless and trace parts is not inconsequential for the physical interpretation of the individual contributions 6 . Namely, although T 00 ,T 00 andT 00 all have the dimension of energy densities, they actually correspond to different thermodynamic potentials as we will show in the following. By focusing on the µ = ν = 0 component in the rest frame, Ji's decomposition does not make any distinction between the Lorentz tensors P µ P ν and M 2 η µν , and therefore disregards pressure effects.
To sum up, although all terms in Ji's decomposition (25) can formally be defined and evaluated, they cannot however be interpreted as pure mass contributions. In the next section, we will explain how to achieve a proper mass decomposition based on a more covariant treatment using the semiclassical picture.

A new decomposition
Consider some decomposition of the EMT at the operator level T µν = i T µν i . From the general parametrization of the matrix element of a rank-two tensor in a state with mass M and spin j ≤ 1 2 [39,40,7], we obtain in the forward limit where A i (0) andC i (0) are two energy-momentum form factors evaluated at zero momentum transfer and depending on the renormalization scale µ. They are related to Ji's coefficients as follows Four-momentum conservation implies that as one can see from Eq. (3). In order to identify the physical meaning of T µν i , we consider its expectation value averaged over the hadron volume V = M V /P 0 , with V denoting the hadron proper volume, In the hadron rest frame, we find where are naturally interpreted as partial proper (internal) energy density and isotropic pressure, respectively. Note that ε i and p i are average quantities over the hadron proper volume, and hence are proportional to the average proper energy density of the hadron M/V . In a generic frame, we can write which has the same structure as the EMT of an element of perfect fluid in relativistic hydrodynamics [41]. Using now the relation (32), we find that Ji's coefficients correspond to the notions of enthalpy density h i and interaction measure I i Remembering that ε i and p i are average quantities over the proper volume, the thermodynamic potentials are then simply defined as where U i is the partial internal energy, H i = U i +W i is the partial enthalpy, and W i = p i V is the partial pressure-volume work. The latter quantity represents the energy associated with the occupation of the volume V by the subsystem i pushing (p i > 0) or pulling (p i < 0) the rest of the system. Since the internal energy represents the sum of all the kinetic and potential energies, the proper mass decomposition we were looking for is Moreover, the total pressure-volume work should vanish for a stable system These two conditions are naturally equivalent to Eq. (33).

Discussion
In the previous section, we observed that the generic form of the EMT for a hadron with spin j ≤ 1 2 is characterized by two Lorentz scalar quantities similarly to that of an element of perfect fluid. We are of course not claiming that hadrons consist of a set of perfect fluids 7 , let alone that a hydrodynamical description is quite difficult to justify in this context. Forward matrix elements of the EMT, like the ones we considered, allow us to determine only static mechanical properties of the state. In the semiclassical limit, these matrix elements can effectively be thought of as a continuum description averaged over time [34,44,45]. We simply adopted the terminology of Continuum Mechanics to identify the physical meaning of the various energy-momentum form factors.
Decomposing the EMT into various contributions T µν = i T µν i amounts in practice to an effective (time-averaged) coupled multifluid picture of the hadron [46], where each species of the mixture is regarded as a separate continuum coexisting with the continuums made up of other species. Each continuum is then described by its own (partial) energy density ε i and pressure p i . Since these quantities are averaged over the hadron proper volume, they trivially obey a barotropic equation of state for species i Depending on the value of the parameter w i , this equation of state describes an element of effective perfect fluid from phantom to ekpyrotic matter [47,48] w i < −1 phantom matter

Back to the old decompositions
Considering the trace of Eq. (35) shows that the trace decomposition in Eq. (6) does not represent a decomposition of the total energy of the system, but rather a decomposition of the interaction measure Each individual contribution involves, beside internal energy U i , the partial pressure-volume work W i . For a stable system, the total pressurevolume work i W i = 0 has to vanish, explaining why one obtains at the end just the total mass M of the system 8 .
Except possibly for the pion, the first term on the right hand side of Eq. (6) corresponding to the gluon contribution is expected to dominate in light hadrons. Contrary to what is often claimed in the literature [17,24,25,26,27], this does not mean however that gluons are at the origin of most of the light hadron masses. It actually indicates that In non-relativistic systems, one can safely neglect pressures in front of energy densities, leading then to ε g ≫ |ε q |, but in relativistic systems like light hadrons, one may expect that pressures become comparable to energy densities. Since stability of the system imposes p q = −p g , assuming an approximate equal sharing of the energy between quarks and gluons ε q ≈ ε g > 0 leads to the conclusion that p q > 0, p g < 0.
We will see later that such a scenario is indeed supported by the current phenomenology. So, in average, quarks are responsible for the repulsive force and gluons for the attractive force inside light hadrons. The two antagonist forces balance each other so as to lead to a stable bound state. This is reminiscent of the structure of stars, where the repulsive force caused by the gas pressure is balanced by the attractive gravitational force 9 .
Hadrons are bound states of two types of constituents, namely quarks and gluons. We can therefore decompose the QCD EMT into quark and gluon contributions T µν = T µν q + T µν g , where although there is some arbitrariness in doing so owing to quark-gluon interactions [7]. Here, the quark mass term is part of T µν q and the trace anomaly is part of T µν g . Accordingly, it seems therefore natural to adopt an effective coupled two-fluid picture where T µν q,g are treated as separate EMT and hence described in terms of their own energy densities ε q,g and pressures p q,g .
Let us now come back to Ji's decomposition. Looking at the expressions (27)-(30), we observe that the individual contributions to the hadron mass correspond to different linear combinations of the coefficients a and b. In the effective coupled two-fluid picture, they correspond therefore to different linear combinations of energy density and pressure owing to Eq. (38). This is of course not physically acceptable, in the sense that one is adding apples and oranges. It is however mathematically correct since the sum over species gives at the end just the hadron mass M , once again thanks to the stability constraint i W i = 0. The only way to make physical sense out of Ji's decomposition within the effective coupled two-fluid picture is to define the quark mass and trace anomaly contributions ε m,a from the onset, and to split the quark and gluon energies as follows with ε q−m ≡ ε q − ε m and ε g−a ≡ ε g − ε a , while keeping the pressure part unchanged. This is consistent, but the price to pay is a loss of covariance since different components of the EMT are not treated in the same way. What Ji did actually in Eq. (13) is a covariant decomposition of the QCD EMT into four parts. He considered with Eqs. (19)-(22) that each quantitȳ T 00 q ,T 00 g ,T 00 m , andT 00 a should represent an energy density, and hence implicitly adopted an effective coupled four-fluid picture of the hadron. The problem with such a description is that e.g. gluons carying kinetic and potential energies are considered as different from those involved in the trace anomaly, in the sense that they are effectively treated as two distinct continuums with their own equations of state. Our opinion is that this is not acceptable from a physical point of view. Since hadrons are composed of quarks and gluons, an effective description of hadrons in terms of only two continuums is more natural.

Virial decomposition
In the effective coupled two-fluid picture, the hadron mass is decomposed into quark and gluon internal energies M = U q + U g , where U q,g = T 00 q,g | P =0 with T µν q,g given by Eqs. (47) and (48). The question now is whether one can further decompose into parton kinetic and potential energies, quark mass and trace anomaly contributions, similarly to Ji's decomposition.
Let us start with the gluon sector. Since the only contribution to the trace of the gluon EMT comes from the trace anomaly, it seems natural to write For the quark sector, we reshuffle the quark mass terms after the separation into traceless and trace parts The constant c is then chosen such thatT 00 q takes the form ψ † (−iD · α)ψ upon using the QCD equations of motion, leading us to Now, contrary to Ji, we will not treatT µν i andT µν i as separate EMT, but rather as mere parts of the EMT T µν i . We then have T 00 In other words, we see that traceless parts contribute to 3 4 and trace parts to 1 4 of the internal energy content. Summing over all the species i naturally leads to Eq. (12). Keeping only the internal energy contributions leads to the following finer decomposition of the hadron mass in the effective coupled two-fluid picture M =Ũ q +Ǔ q +Ū g +Û g (55) with For comparison, Ji's decomposition reads explicitly and differs from our finer decomposition by a reshuffling of the energy contributions within the quark and gluon sectors, separately. We refrain from interpretingǓ q as quark mass contribution andÛ g as trace anomaly contribution. Indeed, in the classical limit γ m → 0 we havě U q = U q which should also include quark kinetic and potential energies. Moreover, the contributionÛ g does not vanish when the trace anomaly is set to zero T µν g = 0. By keeping only the internal energy contributions, we lost the direct connection with the matrix elements and hence a simple physical interpretation. At best, our finer hadron mass decomposition (55) can be seen as some sort of virial decomposition [31,32,33].
A similar treatment can be applied to the pressure-volume work. In the effective coupled two-fluid picture, the partial pressure-volume works of quarks and gluons satisfy 0 = W q +W g , where W q,g = T 33 q,g | P =0 . Treating once again T µν i andT µν i as separate EMT, we get In this case, traceless parts contribute to 1 4 and trace parts to 3 4 of the pressurevolume work. Summing over all the species i, we obtain the analogue of Eq. (12) for pressure-volume work Keeping only the pressure-volume work contributions leads to the following finer decomposition of the hadron stability constraint in the effective coupled two-fluid picture 0 =W q +W q +W g +Ŵ g with For comparison, a similar decomposition in the effective coupled four-fluid picture reads with and differs from our finer decomposition by a reshuffling of the work contributions within the quark and gluon sectors, separately. Once again, we refrain from interpretingW q as quark mass contribution andŴ g as trace anomaly contribution. At best, Eq. (63) can be seen as some sort of virial decomposition.
The above finer decompositions are not based on the nature of the constituents but on Lorentz symmetry and this has a dramatic impact on the physical picture. Indeed, we were not able find a simple physical interpretation for the individual contributions to the decompositions (55) and (63), because the effective coupled two-fluid picture forced us to discard some of the terms in the matrix elements for consistency. In order to avoid this, one has to adopt a picture where the number of effective fluids is at least equal to the number of contributions. In the present case, we need at least an effective four-fluid picture, leading us back directly to Ji's decomposition. We have however already argued that the four-fluid picture is not natural from the physical point of view, since gluons carrying kinetic and potential energies are the same as those involved in the trace anomaly. Similarly, quarks carrying kinetic and potential energies are the same as those characterized by mass m.
Another indication that a decomposition based on Lorentz symmetry is not natural from the physical point of view is the following. Any EMT with the structure (35) and generic equation of state p i = w i ε i can be represented as a particular combination of two EMT with fixed equations of state. This can be seen as some sort of decomposition onto a basis. For example, through a decomposition of the EMT into traceless and trace parts T µν are characterized byp i = −ε i . In Ji's decomposition, the gluon contribution is indeed divided into kinetic and potential energies treated as a pure radiation w g = 1 3 , and trace anomaly treated as a cosmological constant w a = −1. Since the choice of a basis is not unique, the decomposition based on Lorentz symmetry is purely conventional, and hence artificial. In fact, Lorentz symmetry has already been used to provide a physical interpretation of the various components of the EMT, namely by distinguishing in our case energy density ε i from pressure p i in Eq. (35). One can then easily understand why using Lorentz symmetry again to perform a decomposition amounts to choosing an arbitrary basis of EMT with fixed equations of state not determined by the physics of the problem.
To summarize, although a decomposition of the hadron mass or pressurevolume work into four contributions can formally be achieved in the effective coupled two-fluid picture, the individual terms cannot be interpreted in a simple way. If the hadron is seen as made of quarks and gluons, the only unambiguous decompositions are Any further decomposition based on Lorentz symmetry, like e.g. separation of traceless and trace parts, is artificial and purely arbitrary. The only physically acceptable further decompositions are those based on constituent properties, like e.g. flavor or polarization [52,53].

Phenomenology
All that is needed to characterize the quark and gluon EMT (31) are the form factors A q,g (0) andC q,g (0) in the forward limit. Using the energy-momentum sum rules (33), we can reduce this set to e.g. the quark form factors A q (0) andC q (0), or equivalently Ji's coefficients a and b owing to Eq. (32). In the following, we will consider only the proton case and fix the renormalization scale to µ = 2 GeV.
The parameter a, which is also interpreted as the average fraction of hadron momentum carried by quarks, can be extracted from deep-inelastic leptonproton scatterings. A recent global analysis with leading-order parametrization obtained a = 0.546 ± 0.005 [54]. The parameter b is related to the scalar charge of the proton and has been estimated to b = 0.113 ± 0.010 by Gao et al. [2], based on a recent determination of the pion-nucleon σ-term [55], a recent lattice calculation of the strangeness content [56] and neglecting the heavy quark contributions. It has also been suggested to extract the parameter b from the trace anomaly using quarkonium-hadron scattering close to threshold [21,57]. For completeness, the anomalous quark mass dimension is approximatively given by γ m ≈ −0.15 for n f = 3 active flavors [58]. The various decompositions obtained with these values are depicted in Figs. 1-4.
In figure 1 we represent the trace decomposition given by Eq. (6) and determined only by the parameter b. As already discussed, it is largely dominated by the gluon contribution and is at the origin of the claim that most of the nucleon mass comes from gluons [17,24,25,26,27]. The trace of the EMT being given by T µ µ = T 00 − j T jj , what the quark and gluon contributions in Eq. (6) do actually represent are the combinations U q,g −3W q,g . Since the total pressure-volume work vanishes, one artificially supresses one of the contributions in favor of the other with pressure effects. As we will see below, it turns out that W q = −W g > 0 which explains why the quark contribution appears to be much smaller than the gluon contribution. Because of these pressure effects, the sole trace of the QCD EMT does not provide enough information to determine the actual quark and gluon contributions to the hadron mass.
Ji's decomposition (25) is represented in Fig. 2. The M q and M m contributions slightly differ from ones in the pie chart drawn in [2] because we did not neglect the anomalous quark mass dimension γ m . All the information about the hadron EMT in the forward limit, parametrized by Ji's coefficients a and b, is included in this decomposition. As argued in the previous sections, a decomposition of the hadron mass into four contributions is however artificial. In the effective coupled two-fluid picture, the four contributions appear to be combinations of internal energies and pressure-volume works In order to avoid adding apples and oranges, Ji required that each term must represent some internal energy, which amounts to adopting an effective coupled four-fluid picture. In such a picture, no new information arises from the corresponding decomposition of the pressure-volume work (65) represented in Fig. 3, since the individual contributions are directly determined by the energy ones as a consequence of the conventional splitting into traceless and trace parts 10 . In the effective coupled two-fluid picture, they are explicitly given by The information about the energy content and the stability of the system are a priori independent. By treating the traceless and trace parts as separate EMT, Ji effectively combined both information into a single decomposition. While mathematically correct, such a decomposition is however not satisfactory from the physical point of view because gluons involved in the trace anomaly are treated as if they were distinguishable from those carrying kinetic and potential energies (they are indeed associated with different equations of state in Ji's approach).
Since hadrons are made of quarks and gluons, it seems more natural to decompose into two contributions only. Ji's coefficients a and b appear to characterize two different hadron properties of the hadron, namely mass and pressure-volume work. Accordingly, instead of a single decomposition of the hadron mass into four contributions like Ji's, we propose to consider separately the decompositions of hadron mass and pressure-volume work into two contributions, given by Eq. (70) and represented in Fig. 4. We do not require +89% +11% gluon quark 0% 100% Fig. 1 Trace decomposition of the nucleon mass (6)  any splitting into traceless and trace parts, and hence we do not fix a priori the relation between energy density and pressure. Note that all the standard energy conditions are satisfied [59,60] null energy condition ǫ i + p i ≥ 0, weak energy condition ǫ i ≥ 0 and ǫ i + p i ≥ 0, dominant energy condition ǫ i ≥ |p i |, strong energy condition ǫ i + p i ≥ 0 and ǫ i + 3p i ≥ 0.
One is of course free to define quark mass and trace anomaly contributions, but these turn out to be purely conventional and to lead to physical conumdrums when treated covariantly.

Conclusions
We used forward matrix elements of the energy-momentum tensor to characterize static mechanical properties of hadrons. The components of such an energy-momentum tensor can be interpreted semi-classically in terms of parton energy density and pressure averaged over time and the hadron proper volume. We showed that, because of pressure effects, the physical interpretation of the standard decompositions of the hadron mass have to be considered with care. In the trace decomposition, a pressure-volume work contribution artificially emphasizes the role played by gluons. In Ji's decomposition, the splitting of the quark and gluon energy-momentum tensors into traceless and trace parts mixes the information about the hadron mass budget with the pressure-volume work budget. This can be understood by the fact that Lorentz symmetry is already used to provide a physical interpretation of the various components of the energy-momentum tensor. It cannot be used a second time to define separate mass contributions without modifying the physical picture. In particular, from the point of view of physical interpretation, using Ji's decomposition amounts to treating gluons involved in the trace anomaly and those carrying kinetic and potential energies as separate entities with different equations of state, which is not physically acceptable.
Since hadrons are made of quarks and gluons, it is natural to decompose their mass into two contributions only. Any further decomposition that is not based on the properties of the constituents will be somewhat arbitrary and will mix internal energy with pressure-volume work. This mixing is mathematically harmless because the total pressure-volume work vanishes for a stable system, but it is a problem for the physical interpretation of the individual contributions.
We proposed a new picture, where the hadron mass and pressure-volume work budgets are kept separate and expressed in terms of the sole quark and gluon contributions which are physically unambiguous. In particular, unlike Ji's decomposition we do not fix a priori the equations of state for quarks and gluons.
Finally, we quantitatively compared the different decompositions based on recent phenomelogical estimates. It turned out that, as expected for highly relativistic systems, pressure-volume work contributions are of the same order of magnitude as internal energy contributions. In particular, quarks are responsible in average for the repulsive force and gluons for the attractive force inside nucleons.
Having clearly identified the pressure contributions opens many interesting applications related to compact stars. For example, determining the quark and gluon equations of state inside a nucleon may give important clues about the internal structure of compact stars.