Proton Mass Decomposition and Hadron Cosmological Constant

Lattice results on sigma terms and global analysis of parton momentum fractions are used to give the quark and glue fractions of the proton mass and rest energy. The mass decomposition in terms of the trace of the energy-momentum tensor is renormalization group invariant. The decomposition of the rest energy from the Hamiltonian and the gravitational form factors are scheme and scale dependent. The separation of the energy-momentum tensor into the traceless part which is composed of the quark and glue parton momentum fractions and the trace part has the minimum scheme dependence. We identify the glue part of the trace anomaly $\langle H_{\beta}\rangle $ as the vacuum energy from the glue condensate in the vacuum. From the metric term of the gravitational form factors, which is the stress part of the stress-energy-momentum tensor, we find that the trace part of the rest energy, dominated by $\langle H_{\beta}\rangle$, gives a {\it constant} restoring pressure which balances that from the traceless part of the Hamiltonian to confine the hadron, much like the cosmological constant Einstein introduced for a static universe. From a lattice calculation of $\langle H_{\beta}\rangle$ in the charmonium, we deduce the associated string tension which turns out to be in good agreement with that from a Cornell potential which fits the charmonium spectrum.


Introduction
In this work, we consider the decompositions of the proton mass and rest energy into their respective quark and glue components and will use lattice results and momentum factions from a global analysis to provide numerical results for each component. Even though the mass and rest energy are equal in Einstein's equation i.e. E 0 = mc 2 , many properties associated with the mass and energy are not the same. For example, the mass is a Lorentz scalar while the energy is the time component of the 4-momentum vector. In the example of e + e − annihilation to two photons e + e − −→ γγ, the mass of the two photons from the rest energy is 2m e , not the sum of the two photon mass [1,2]. This shows that, while momentum and energy have additivity properties, mass does not. When there is mass there is energy, but not vice versa. For nonrelativistic particles, mass appears in Newtonian dynamics and gravitational interaction. In general relatively, however, the gravitational field is coupled to the energy-momentum tensor that has 10 components. 1 In the present work, we shall make a distinction between the mass and rest energy in the context of separating the quark and glue components in the proton or other hadrons. As we shall see in the following text, the mass can be obtained from the trace of the energymomentum tensor (EMT). However, the fraction of each component, although renormalization group invariant, needs to be defined in the rest frame. On the other hand, the decomposition of the proton rest energy, be it through the Hamiltonian or the forward matrix elements of the gravitational form factors, is scheme and scale dependent. We shall evaluate the matrix elements associated with the operators to examine the physical meaning of them and; moreover, check if they are accessible experimentally and/or through lattice calculations.
The paper is organized as follows. In Sec. 2, we will discuss the proton mass decomposition via the trace of the energy-momentum tensor and use lattice QCD calculation of the sigma terms for quarks with different flavors to find the fractional contribution for each flavor as well as the trace anomaly. In Secs. 3 and 4, we will describe the decomposition of the rest energy in terms of the Hamiltonian and the forward matrix elements of the gravitational form factors. They are related to the quark and glue momentum fractions of the parton distribution functions which are obtained from the global analysis of experiments. As such, they are scale dependent. There is quite a bit of interest lately in the issue on how to understand the physical meaning of each component in various ways of apportioning the proton mass. The first work was done by X. Ji who has decomposed the hadron rest energy from the Hamiltonian in terms of the quark sigma terms, the quark and glue energies and the anomaly [3,4]. Lattice calculations based on this decomposition have been carried out for the mesons [5] and the nucleon [6]. Also, combined lattice and experimental results on the momentum fractions have been used to evaluate the quark and glue components of the nucleon [7]. Examining the rest energy from the gravitational form factors, C. Lorcé has interpreted the quark and glue parts in terms of the internal energies and pressurevolume work [8,9]. Y. Hatta, A. Rajan and K. Tanaka have explored the renormalization group property of the form factorC and performed perturbative calculation of the quark and glue parts of the trace of the energy-momentum tensor [10]. The scheme dependence of the quark and glue parts of the mass and rest energy has been further examined by A. Metz, B. Pasquini and S. Rodini [11]. Each of the above-mentioned work has some different expressions and interpretations in the decomposition which we will attempt to address and clarify. We shall discuss the significance of each term in the gravitational form factors in terms of thermodynamics. In Sec. 5, we discuss the role of the scalar trace and tensor traceless parts of the rest energy in the metric term (η µν ) of the gravitational form factor which is the stress in the energy-momentum tensor. Comparing them in the energy and pressure equations, we find that the scalar part of the energy, dominated by the glue part of the trace anomaly H β , has an energy density which is a constant so that it gives rise to a constant restoring pressure to balance those from the quark and gluons and, thus, confines the proton. This shows that it is the hadron cosmological constant. As a support of this idea, we show that the string tension deduced from a lattice calculation of H β in the charmonium agrees well with that from a Cornell potential for the charmonium spectrum calculation. A summary is presented in Sec. 6.

Mass Decomposition
The decomposition of proton mass and rest energy are all based on the energy momentum tensor (EMT). The Belinfante form of the EMT is a symmetric rank two tensor The bare operators are and where T a is the SU (3) color matrix.
It is natural to consider the mass decomposition in terms of the trace of the EMT, since the matrix element of the trace gives the nucleon mass and is frame independent On the other hand, the forward matrix element of the EMT component is One can consider the rest energy from the T 00 component.
As far as the EMT trace is concerned, the classical trace is zero when the quark mass is neglected. However, this conformal symmetry is broken by a trace anomaly in QCD due to quantum corrections [12,13,14,15], where β(g) is the β function and γ m is the mass anomalous dimension. In dimensional regularization, the quark condensate comes from the quark part T µ qµ in Eq. (2) and the anomaly terms with γ m and β are from T µ g µ , i.e.
It is pointed out that the above separation into quark and glue parts of the trace is scheme dependent and they can mix under renromalization giving rise to scheme dependence [10,16,11]. However, since the EMT is conserved, i.e.
Therefore, the renormalized matrix elements of the trace are represented as where ... = P |...|P . As long as one does not take the renormalized T µ q µ and T µ a µ in Eqs. (8) and (9) as separate quark and glue contributions, but considers the renormalized (T µ µ ) R as a whole, the expression in Eq. (12) is scheme independent and renormalization group invariant. However, the matrix element is proportional to M 2 , not M itself. This raised the question and ambiguity about its proper normalization [9,10,11,17,18]. The usual definition of an expectation value P 0 depends on the frame, since P |P = (2π) 3 2P 0 δ 3 (0). One exception is when all the components are in the rest frame. In this case, In the moving frame, one should take the integral d 3 x γ to be over the proper volume so that This shows that mass is frame independent.
In this work, we shall define the fractional contribution of each term in Eq. (12) where the quark mass times the quark condensate is know as the sigma term and the fraction f N f is thus the mass fraction of the sigma term The nucleon sigma terms have been calculated on the lattice and they are tabulated in FLAG [19] . We shall use the results from the overlap fermion [20] that has exact chiral symmetry on the latticeà la Gingparg-Wilson relation [21]. In this calculation [22], several lattice ensembles with different lattice spacings and quark masses with one lattice at physical pion mass are used to address the systematic errors in the continuum and infinite volume extrapolations at the physical pion point. For QCD at hadronic scale ∼ 1 GeV, the heavy quarks are integrated out so that they are reflected in the β function. In Table 1, we give f N πN and f N s for the (2 + 1)-flavor case. The trace anomaly contribution f N a is obtained from subtracting the f f N f from unity. f N a can in principle be calculated on the lattice which involves more elaborate renormalization [23] and will have larger errors than that obtained here from the mass sum rule. As we see from Table 1 and Fig. 1(a) that the sigma terms contribute only ∼ 9% to the nucleon mass. The rest comes from the trace anomaly.
When one considers QCD as part of the standard model at the weak scale, the Higgs and heavy quarks (c, b, t) are all involved as relevant operators and can appear in external states. This is relevant to high energy processes such as DIS and high energy hadron collisions. If the dark matter, such as neutrilino, couples to the Higgs, then the heavy quark sigma terms will contribute to the scattering cross section of the dark matter on nucleus [24]. It was pointed out [25] that, to the leading order in the heavy quark expansion, the matrix element m Q P |QQ|P is related to the glue condensate in the nucleon, This is precisely the n f term in the leading α s /4π expansion of β(g) [26,27], where β 0 = (11 − 2 3 n f ), except with a negative sign. This shows that for n f heavy enough quarks, the introduction of their sigma terms can trade with the n f term in β 0 , the leading term in β(g)/2g. To study the quark mass dependence, a lattice calculation with the overlap fermion has been carried out [28]. It is found that the sigma terms for quark masses heavier that ∼ 1/2 of the charm mass are the same within errors. We take this finding to mean that the sigma terms for the charm, beauty and top quarks are the same. For the charm, it is found [28] that f N c = 0.094 (31), which is taken to be the same for f N b and f N t . The sigma terms for the heavy quarks are listed in Table 1 and the pie chart for the quark sigma terms for all six flavors and the corresponding trace anomaly term is given in Fig. 1(b). Even though the heavy quark sigma terms are listed in Table 1 and plotted in Fig. 1(b), it should not be misconstrued to imply that they are the total heavy quark contributions to the proton mass. As explained above, the heavy quarks also contribute negatively through the β function with n f flavor [cf. Eq. (17)] so that the net contribution of a heavy quark with mass M H is O(1/m H ), in accordance with the decoupling theorem [29,30]. Table 1: Decomposition of proton mass in terms of the quark sigma terms of different flavors and the trace anomaly. They are given as the percentage fractions of the proton mass. The sigma terms are obtained from lattice calculations [22,28].

Rest Energy from Hamiltonian
We shall consider the rest energy from the T 00 component of the EMT matrix elements at rest. There are several ways to examine the decomposition into quark and glue contributions. The decomposition of the proton rest energy with the Hamiltonian was first considered by X. Ji [3,4]. Since the EMT in the Belinfante form is a symmetric rank two tensor, they can be separated into the traceless and trace parts. They are in different representations of the Lorentz group. The traceless and trace parts are in the irreducible (1,1) and (0,0) representations respectively. Hence, they do not mix. The renormalized EMT is separated as where (T µν ) R = 1 4 η µν T ρ ρ . So far, this separation is scale and scheme independent.T µν can be further split into the quark and glue parts. In this case, the Hamiltonian, being the spatial integral of T 00 , i.e. H = d 3 x T 00 (x), can be written as where Here, the subscript R indicates it is renormalized and mixed if needed. For example, the separation of H q and H g from the traceless T µν entails operator mixing besides renormalization and the scale dependence is introduced. These nucleon matrix elements can be extracted from experiments and lattice calculations at the scale µ. H q (µ) is the quark mass and energy operator and H g (µ) corresponds to the glue field energy operator. Their matrix elements are related to the quark and glue momentum fractions where x f (µ) and x g (µ) are the momentum fractions for the quark with flavor f and the glue from DIS or Drell-Yan experiments at the scale µ and they satisfy momentum conservation, i.e. f x f (µ) + x g (µ) = 1. This decomposition gives a sum rule for the nucleon rest energy Here, we have defined H ... = P |H ... |P / P |P at P = 0. This is the simplest and welldefined separation of the rest energy from the Hamiltonian-each term is related to experimental and/or lattice observables, with the necessary scale dependence in the MS scheme and there is no further scheme dependence for the trace term. We shall use x f (µ) and x g (µ) from CT18 global analysis of the PDFs [31] for the rest energy decomposition. They are listed in Table 2 for µ = 2 and 250 GeV. As can be seen in Table 2, the momentum fractions of the u and d partons, which are mostly due to the valence contributions, are shifted to those of the sea-quark and glue partons as the scale µ increases from 2 to 250 GeV, the electroweak scale. We shall define the fractions and plot them in Fig. 2(a) for µ = 2 GeV and Fig. 2(b) for µ = 250 GeV. We note that f H tr can be further decomposed as in Eqs. (13) and (14). It is shown that this separation can have a Lorentz-invariant interpretation [8,32,17].
There are suggestions to separate H tr into the sigma terms and the trace anomaly [3,4,11] so that the Hamiltonian is where The quark and glue momentum fractions x f (µ) and x g (µ) from CT18 global analysis of the PDFs [31] are tabulated in the M S scheme for µ = 2 GeV in the second row and µ = 250 GeV in the third row [33]. They are given as percentages of the total proton momentum.
x f,g (250 GeV) 24.5(4) 15.2(4) 5.6 (9) 4.3(9) 2.9(4) 0 52.5(1.3) 47.5 (5) However, it was pointed out that the separation of H tr introduces a scheme dependence [10,16,11] and the separation suggested here amounts to taking a specific D2 scheme [11]. Nevertheless, this is a physically motivated separation in the sense that H m and H a are each scale invariant and they can be determined from lattice calculations as we alluded to in Sec. 2. It has also been suggested [3,4] The classical form of H E (µ) is ψ † i α · Dψ and can be interpreted as the quark kinetic and potential energy. The renormalized matrix element of H E has mixings with both H g and H m . The mixing of the traceless matrix elements is where Z qq (µ) = 1 − Z gq (µ) and Z gg (µ) = 1 − Z qg (µ) to preserve momentum conservation and µ r is a reference scale for the unmixed matrix element, such as the lattice spacing in lattice calculations. The renormalization group equation for the evolution of the moments, to first order, is where C F = 4/3. The renormalized and mixed (H E ) R (µ) is then The first term on the right-hand side is the self-renormalized operator and the terms in the square bracket are from mixing with H m and H g . The matrix element of H E (µ) in terms of x (µ) and H m turns out to be which is the same as in Refs. [3,4] so that the rest energy sum rule is satisfied.
In an effort to have a simplified decomposition with physical interpretation for each term, Metz, Pasquini and Rodini [11] have derived a decomposition based on dimensional regularization in the MS scheme. It is composed of three terms which represent the quark mass, quark energy and glue energy respectively. They have the same formal expressions for the corresponding operators as in Eqs. (31), (21), and (36) for H m , H g and H E (without the perturbative Z terms), but without H a as compared to Eq. (33). The difference appears to be due to the definition of the operators. When the matrix elements are involved, which is the sum of H g and 1 4 H a . Thus, in the end, the same sum rule is satisfied as in Eq. (38).

Rest Energy from Gravitational Form Factors
The rest energy has also been discussed in the context of the gravitational form factors. These are form factors for the EMT and contain the following terms [34,35,36,37] for the quarks and gluons where T 1q,g (0) = x q,g (µ) is the momentum fraction and T 1q,g (0) + T 2q,g (0) = 2J q,g (µ) the angular momentum fraction [36,37]. From the conservation of momentum and angular momentum, one has the following sum rules [36,37] T 1 and T 2 are named after F 1 and F 2 of the electromagnetic form factors due to their analogy to the form factor structure of the vector current, where the forward F 1 (0) is the charge and forward F 1 (0) + F 2 (0) = µ p,n the magnetic moments of the nucleon. By making a connection to the stress tensor of the continuous medium, it is shown by Polyakov [38] and Polyakov and Schweitzer [39] that D q,g (q 2 ) is related to the internal force of the hadron and encodes the shear forces and pressure distributions of the quarks and glue in the nucleon. The D term, which is D(0), has been emphasized [39] to be the least known fundamental constant of the nucleon as compared to charge, magnetic moments, g A , g P , etc. The pressure distribution of the quarks have been deduced from the experimentally measured D q (q 2 ) [40] and the pressure distributions for both the quarks and glue from D q,g (q 2 ) are calculated on the lattice [41].
In view of the similarity between the structure of the gravitational form factor in Eq. (41) and that of the stress tensor of the perfect fluid frequently used in general relativity where u µ = P µ /M , C. Lorcé identified the quark and glue parts of the EMT in Eq. (41) as two fluids [8,9] with q,g ≡ T 1q,g (0) +C q,g (0) where V is the proper volume. Therefore, the rest energy has been interpreted in terms of the following thermodynamic functions Here U q and U q are the internal energies for the quarks and glue and they are equal to the matrix elements of H q + 1 4 H m and H g + 1 4 H a in Eq. (30), respectively. The work W is W q,g = p q,g V = −C q,g (0)M.
Since the enthalpy is H = U + W , one finds the enthalpies for quarks and glue are The rest energy and total work are the sum of the quark and glue internal energies and their works In view of Eqs. (55) and (52), one finds the total pressure to be zero, i.e. P = p q + p g = 0. It is thus concluded by C. Lorce that the quark gives a positive pressure (p q > 0) and is balanced by the negative pressure from the glue part (p g < 0) to have a stability condition for the nucleon, and hadrons in general.
It has been pointed out by Y. Hatta, A. Rajan and K. Tanakain [10] that the renormalization and trace taking of individual quark and glue EMT do not commute in dimensional regularization, i.e., η µν (T µν q,g ) R = (T µ q,g µ ) R , and the perturbative mixing of η µν T µν q,g has been worked out to two loops [10] and three loops [43]. Taking this into account,C q,g in one loop becomes Using γ m 0 = 6 C F and α s (2 GeV) = 0.301 and α s (250 GeV) = 0.102, we obtainC q (2 GeV) = −0.161(10),C q (250 GeV) = −0.122 (47), which are close to the asymptotic value of −0.146 for n f = 3 in the chiral limit as calculated in Ref. [10]. Separating out the quark from the glue contributions in the trace, it introduces a scheme dependence inC(0, µ) as is in the case of decomposition with the Hamiltonian in Eq. (30). On the other hand, if one does not separate out H m from H a in the trace, the rest energy from the internal energies has the same decomposition as the rest energy in Eq. (26) from the Hamiltonian.

Trace Anomaly, vacuum energy, and Cosmological Constant
Throughout the discussion of the nucleon mass and rest energy decompositions, the role of the trace anomaly remains mysterious. Besides giving rise to scale breaking, it is not clear what physical role it plays as far as the hadron mass and structure are concerned. In the mass decomposition, the hadron mass is entirely due to the trace anomaly at the chiral limit in Eq. (12). This does not reveal any dynamical information of the trace anomaly. On the other hand, as seen from the decomposition of the rest energy from the Hamiltonian in Eqs. (26) and (30) and the gravitational form factors in Eq. (59), there are other terms besides the trace anomaly that may make them more susceptible to divulging the dynamical origin and function of the trace anomaly.
As we shall see, it would be helpful to describe the nucleon in thermodynamic principles. For that, we need to first establish that the nucleon can be treated as a statistical system. To this end, we shall look at the nucleon from the path-integral formulation of QCD in the Euclidean space. QCD in Euclidean path-integral formulation has the same form as that of the classical statistical mechanics. The vacuum-to-vacuum transition in the quantum field theory is termed the grand canonical partition function where S G and S F are the gauge and fermion actions. T is the temperature and µ is the chemical potential. There are infinite degrees of freedom in the partition function. Extensive lattice calculations have been carried out to study the QCD phase diagram in finite temperature and chemical potential (for a review of the status of the lattice calculation, see for example Ref. [44]). The canonical ensemble approach has been formulated [45,46] to address the problem with a definite baryon number. This entails the projection of the fermion determinant from Z GC (V, T, µ) with imaginary µ. Since the center Z 3 symmetry is preserved in the canonical approach, the projected quark numbers are multiples of 3, thus the canonical partition Z C (V, T, n B ) is a function of the baryon number n B . The chemical potential at fixed baryon number n B is where F n B is the free energy. For one nucleon, µ (n B = 1) is the nucleon mass, i.e. µ (n B = 1) = M when T → 0. This is the same as calculating the nucleon mass from the twopoint nucleon correlator with t → ∞ in the grand canonical ensemble in Eq. (60) at µ = 0. The quark matrix element of the nucleon can be similarly obtained with projected quark propagators [46]. A first-order phase transition at finite density and temperature is clearly seen with the Maxwell construction and the critical point is determined [47]. However, this was done on a small 4 4 lattice. Once the volume is increased, the sign problem sets in abruptly that impedes the lattice calculation.
Given the above formalism, we can discuss the nucleon thermodynamic properties from the canonical ensemble. We first notice that theC terms are associated with the metric η µν in the gravitational form factor in Eq. (41) which are like the cosmological constant, except that the sumC q +C g = 0. However, we need to examine Eq. (46) further and point out that the origin of f N f and f N a is different from that of x q and x g inC q andC g . Instead of decomposing the EMT in terms of quark and glue contributions, it turns out to be more fruitful to examine the separation of the traceless and trace contributions in Eq. (18). We first note that there is a glue condensate in the vacuum with a negative energy density (i.e. 0|G 2 |0 < 0) so is there a negative quark condensate in the vacuum (i.e. 0|ψψ|0 < 0). The former is from the conformal symmetry breaking and the latter from the chiral symmetry breaking. One can picture the hadrons as bubbles in the sea of the condensates. The hadron gas can go through a first-order phase transition to become the quark-gluon plasma at finite temperature and chemical potential, much like the bubbles turning into steam during the first-order phase transition between water and steam. As such, the trace anomaly can be considered the vacuum energy of the hadron due to the fact that the bubble with a finite proper volume is created from the surrounding condensate with negative energy. In other words, where vac = − 0|H a |0 is the vacuum energy density and V is the proper three-volume of the hadron. Similarly, the disconnected insertion part of the sigma terms arise from the quark condensates in the vacuum. Since H m is small compared to H a -especially the connected insertion part, in the nucleon at the hadronic scale as shown in Table 1 and Fig. 1(a), we shall ignore it. As the volume changes, we have, from the first law of thermodynamics at zero temperature where dQ = T dS = 0. Here the pressure P vac is which is negative. On the other hand, the other energies H E (µ) + H g (µ) in Eq. (38) or x q,g in Eqs. (27) and (28) are ascribed to the quark kinetic and potential energy and the glue field energy. Their energy densities should go down if the hadron is allowed to expand so that their energy densities are like those in the cosmological models where the matter density falls off like 1/a 3 and radiation density falls off like 1/a 4 [48], where a is a scale factor for the radius of the Universe.. This will give a positive pressure to balance that from the H a so that the nucleon and other ground state hadrons in QCD are stable. We can parametrize the rest energy (or internal energy) of the nucleon with two terms to illustrate the situation.
where the first term is from the anomaly that is proportional to the volume and the second term represents the quark kinetic and potential energies and the glue field energy as in Eqs. (38) and (26). It is proportional to V p with p < 0. Taking the derivative with respect to V gives the pressure which should be canceled between the two contributions Given that E 0 = 4E S from the Virial theorem consideration [17], where E S = 1 4 ( H a + H m ) is the scalar part of the energy, one obtains p = − 1 3 . This is the consequence of the Virial theorem for the scalar and tensor energies. One notices that Eq. (65) is exactly like the MIT bag model [49] where vac is the bag constant B which provides the confinement and mat V −1/3 corresponds to the quark and gluon normal modes in the bag cavity which are proportional to 1/R where R is the bag radius [50].
More importantly, we notice that the pressure-volume equation in Eqs. (43), (46) and (55) can be written as where E T and E S are the tensor (traceless) and scalar (trace) parts of the energy T 00 in Eq. (18), i.e. where Notice that the 1 3 factor in Eq. (67) comes from the fact that E T = 3E S and it coincides with the p = − 1 3 from the volume dependence in Eq. (66). Thus the two pressure equations, Eqs. (67) and (66), are consistent with each other so that one can conclude that the scalar energy density E S V = vac is a constant. On the other hand, the tensor energy density is not a constant and E T is proportional to V −1/3 . Since E S appears in the metric term in the EMT, it yields a constant restoring pressure to balance that from the tensor energy E T , leading to the confinement of hadrons. This is analogous to the cosmological constant Λ that Einstein introduced to his equation in general relativity [51] where R µν is the Ricci curvature tensor and R is the scalar curvature. G is the Newton's constant and the source T µν is the energy-momentum tensor. The positive constant Λ is introduced to the g µν term so that it balances the gravitational pull from a static uniform matter density ρ. Einstein found the solution of Λ to be [51,52] Later on, when the Universe was found to expand by Hubble, Einstein considered the introduction of the cosmological constant the biggest blunder of his life. The confinement mechanism with the hadrons is somewhat different from that of the the cosmological constant Λ in Eq. (71) in that the pressure from the energy density of E S is negative so that it balances the positive pressure exerted by the quark and glue kinetic energy from E T . There have been suggestions that the glue part of the trace anomaly is the vacuum energy which is responsible for the confinement of hadrons [53,54,55,4,56]. We have proved from the stress-pressure equation in Eq. (67) that the energy density of E S is indeed a constant that leads to a constant restoring pressure that is responsible for the confinement of hadrons.
Heavy quarkonium is an analogous but more familiar case. Since the lattice simulation reveals an area law of the Wilson loop that displays a linearly confining potential between the infinitely heavy quark-antiquark pair [57] and the flux tube is formed along the separated color sources [58], they are consistent with a picture of a constant vacuum energy density for H a . When the flux tube has a fixed cross section A, the static potential between the heavy quark-antiqurak pair is which is the origin of the linearly rising potential, where σ is the string tension. We can estimate the string tension σ from the G 2 term of the trace anomaly in the charmonium.
The potential between infinitely heavy quark-antiquark pair can be calculated from the rectangular Wilson loop W L (r, T ) which is a plaquette with time T and spatial extent r.
It is shown that under the renormalization group consideration, the potential and its derivative is related to the G 2 term of the trace anomaly in the quenched approximation [59,60] For the charmonium, we can check this relation to see how well it holds. In this case, the corresponding matrix element for the right-hand side of Eq. (75) is For the linear potential in Eq. (73), this implies for the potential energy. H β cc can be obtained from the lattice calculation of H m cc , Taking M cc to be the spin average of the J/Ψ and η c masses, M cc = 3069 MeV. The sigma term in the charmonium in a recent lattice calculation of the charmonium structure [61] has been obtained to be H m = 2166(1) MeV. At the scale of the lattice spacing a −1 = 1.785 GeV used in the the lattice calculation [62,61], γ m (µ = 0.1785 GeV) = 0.325. This is consistent with a fit of the charmonium spectrum with the Cornell potential that determined the effective α s = 0.491(80) [63], and this implies that the leading order γ m = 6 C F αs 4π = 0.31. With γ m = 0.325, one can deduce from Eq. (78) that H β = 199 MeV.
To estimate r in Eq. (77), we shall use the ratio r r 2 1/2 . Since the asymptotic behavior of the 1S wave function of a linear potential in the Schrödinger equation is between that of the harmonic oscillator and the Coulomb potential, we shall use the average of the above ratios from the r 2 and −1/r potentials, i.e. r r 2 1/2 = 1 2 [ r r 2 1/2 | H.O. + r r 2 1/2 | Coul ] = 0.61. In a potential model to fit the charmonioum spectrum which includes both the linear and Coulomb potentials and the spin-spin and spin-orbit interactions, r 2 1/2 = 0.21fm [64]. This gives r ∼ 0.13 fm. From Eq. (77) and H β from the lattice calculation in Eq. (78), we obtain the string tension σ = 0.153 GeV 2 . This is in very good agreement with that from the recent Cornell potential analysis of the charmonium spectrum which gives the fitted string tension σ = 0.164(11) GeV 2 . Despite this close agreement, we should caution that Eq. (75) is for infinitely heavy quarks in a pure gauge theory, while the lattice calculation [62,61] is carried out with realistic finite charm quark mass on (2+1)-flavor dynamical fermion configurations.
An effective low-energy theory in curved spacetime with broken scale invariance has been formulated which can lead to confinement [65,66]. The trace anomaly (quantum anomalous energy) contribution to the nucleon mass has been considered in a Higgs mechanism [56]. Here, we identify 1 4 H β as the vacuum energy emerged from the formation of a hadron bubble in the sea of the glue condensate. We should emphasize the fact that the vacuum energy density is a constant, so is its negative restoring pressure which inevitably results in confinement. This is similar to the case where the constant force from a linear potential confines the heavy quarkonium.
The cosmological constant in general relativity has a renewed interpretation after the Universe is found to have an accelerating expansion. The pressure from the energy-momentum tensor also gravitates and contributes to the acceleration of the expanding Universe. The Freidman equation from the Friedmann-Robertson-Walker metric is where a is a scale factor, ρ/P is the density/pressure from all matter, radiation, and dark energy. The cosmological constant, which would appear as Λ/3 on the right-hand side for historical reasons, has been subsumed into the ρ and p as vacuum energy and pressure with ρ vac = Λ/(8πG) and p vac = −Λ/(8πG). When the negative pressure from the dark energy (cosmological constant) overcomes ρ, the Universe expansion accelerates. In the case of QCD, the negative pressure P vac simply balances out the pressures from the quark kinetic energy and glue field energy and confines the quarks and gluons. In this sense, the hadron is analogous to the static universe that Einstein had imaged for the cosmological constant to achieve and E S /V or, more specifically, H β plays the role of the cosmological constant for the hadrons.
It would be interesting to observe the glue part of the trace anomaly H β experimentally [67]. It is shown by D. Kharzeev [68] that the photoproduction of J/Ψ at threshold would be a place to probe this. It has also been explored by Y. Hatta and D.L. Yang, using gauge/string duality [69]. One could also explore the trace anomaly and conformal symmetry on the lattice. Through the study of the Dirac eigenvalue density, there is evidence that there is a phase above the crossover temperature that displays infrared scale invariance [70]. It would be useful to find out what impact it may have on H β . There are efforts to look for conformal window with multiflavor simulations [71]. One could ask the same question about H β . Also, using it as an indicator, one could calculate it in nuclei to see if the conformal symmetry is partially restored.

Summary
We have considered the decomposition of the proton mass and rest energy in terms of their quark and glue components. The proton mass from the trace of the EMT has a unique division in terms of the quark sigma terms and the quantum trace anomaly that are scheme and scale independent. Since the mass is not additive in general frames, the decomposition of the proton mass in terms of expectation values can be carried out in the rest frame or a comoving frame. We use the lattice results to enumerate each term. The role of the heavy quarks is clarified.
There are different ways to decompose the rest energy, be it from the Hamiltonian or the gravitational form factors. They are scheme and scale dependent. The simplest and least scheme dependent way is to divide the Hamiltonian in terms of the trace and the traceless parts, which is the same as obtained from the forward gravitational form factors. The traceless part can be separated into quark and glue momentum fractions, measurable from DIS and Drell-Yan experiments. We use the CT18 global fitting of x q and x g to display this decomposition at µ = 2 and 250 GeV. We note that the trace anomaly introduces a scale in QCD and its mass decomposition can be expressed in terms of the sigma terms and the trace anomaly. The proton mass from the trace of the EMT is simply four times the scalar energy, i.e. M = 4E S .
One interesting feature is revealed from the decomposition of the gravitational form factors. The forward η µν (metric) termC q,g is the normal stress of the EMT which is the negative of the pressure for the system. The total pressure from the trace part of the rest energy (E S ) and the traceless part (E T ) is zero due to the conservation of the EMT. It shows that the pressure from E S is canceled by − 1 3 of that of E T . This reflects the fact that E T = 3E S . More importantly, expressed in terms of their volume dependence, it discloses that E S is linearly dependent on volume which gives a constant restoring pressure to balance the positive pressure from E T which has a volume dependence of V −1/3 . E S with constant energy density is naturally interpreted as the vacuum energy from the condensate; in particular, the glue part of the trace anomaly, which dominates E S , is due to the forming of a nucleon bubble in the sea of the glue condensate. Einstein introduced the cosmological constant in general relativity with an intention to describe a static universe. It is actually more applicable to delineate the confinement of hadrons.
The linear potential between the heavy quark-antiquark pair can be understood as due to a constant vacuum energy density and a flux tube formation between the heavy quarks. We deduce the string tension from the glue part of the the trace anomaly H β cc in a lattice calculation of the charmonium and found it to be in very good agreement with that from a Cornell potential that fits the charmonium spectrum. This further supports the notion that the glue part of the trace anomaly is the hadron cosmological constant which is responsible for hadron confinement. Studies of the trace anomaly in conformal symmetry restored phase in QCD and in conformal field theories may shed light on this issue.