arXiv:hepph/9410274v1 12 Oct 1994
A QCD ANALYSIS OF THE MASS STRUCTURE OF THE NUCLEON?
Xiangdong Ji
Center for Theoretical Physics Laboratory for Nuclear Science
and Department of Physics Massachusetts Institute of Technology
Cambridge, Massachusetts 02139
MITCTP: #2368
HEPPH: #9410274
(Submitted to: Physical Review Letters
October 1994)
Abstract
From the deepinelastic momentum sum rule and the trace anomaly of the energymomentum tensor, I derive a separation of the nucleon mass into the contributions of the quark and gluon kinetic and potential energies, the quark masses, and the trace anomaly.
Typeset using REVTEX
?This work is supported in part by funds provided by the U.S. Department of Energy (D.O.E.) under cooperative agreement #DFFC0294ER40818.
1
The nucleon derives its mass (939 MeV) from the quarkgluon dynamics of its underlying structure. However, due to the complexity of the lowenergy Quantum Chromodynamics (QCD), a more detailed understanding of the nucleon mass seems di?cult. The lattice QCD is successful in reproducing the measured mass from the fundamental lagrangian [1], but the approach provides little insight on how the number is partitioned between the nucleon’s quark and gluon content. Years after the advent of QCD, our knowledge of the nucleon’s mass structure mostly comes from models: nonrelativistic quark models, Bag models, the Skyrme model, string models, the NambuJonaLasinio model, to just name a few. Though all the models are made to ?t the mass of the nucleon, they di?er considerably on the account of its origin. Depending on di?erent facets of QCD the models are created to emphasize, the interpretations of the nucleon mass often go opposite extremes.
In this Letter I show that an insight on the mass structure of the nucleon can be produced within QCD with the help of the deepinelastic momentum sum rule and the trace anomaly. The result is a separation of the nucleon mass into the contributions from the quark, antiquark, gluon kinetic and potential energies, the quark masses, the gluon trace anomaly. Numerically, the only large uncertainty is the size of P mss?sP , the strange scalar charge of the nucleon. Some implications of this breakup of the masses are discussed following the result.
Let me begin with the energymomentum tensor of QCD,
T ?ν
=
1 ψ?i?D(?γν)ψ 2
+
1 g?ν F 2 4
?
F
?αF
ν α
,
(1)
where ψ is the quark ?eld with color, ?avor, and Dirac indices; F ?ν is the gluon ?eld strength
with color derivative
indices and F 2 = ?D? = →D? ? ←D?,
wFiαthβ F→Dαβ?;
and all = →?? +
implicit indices are summed igA? and ←D? = ←?? ? igA?,
over. The where A?
covariant = A?a ta is
the gluon potential. The symmetrization of the indices ? and ν in the ?rst term is indicated
by (?ν). Eq. (1) is quite formal, for it contains neither the gauge ?xing and ghost terms, nor
the trace anomaly. The ?rst type of terms are BRSTexact [2] and have vanishing physical
matrix elements according to the JoglekarLee theorems [3]. I will add the trace anomaly
explicitly when the renormalization issue is dealt with.
A few results about the energymomentum tensor are wellknown. First of all, it is a
symmetric and conserved tensor,
T ?ν = T ν?, ??T ?ν = 0.
(2)
Because of the second property, the tensor is a ?nite operator and does not need an overall renormalization. All the ?elds and couplings in Eq. (1) are bare and their divergence are cancelled by the standard set of renormalization constants. The only complication is the gluon part of the tensor cannot be renormalized with a vanishing trace (see below). Second, the tensor de?nes the hamiltonian operator of QCD,
HQCD = d3x T 00(0, x),
(3)
which is also ?nite and scaleindependent. Third, the matrix element of the tensor operator in the nucleon state is [4],
P T ?νP = P ?P ν/M,
(4)
2
where P is the nucleon state with momentum P ? and is normalized according to P P = (E/M)(2π)3δ3(0), and E and M is the energy and mass of the nucleon, respectively. Lastly,
the trace of the tensor is [5]
T??ν
=
1 g?ν 4
(1 + γm)ψ?mψ +
β(g) 2g
F
2
,
(5)
where m is a quark mass matrix, γm is the anomalous dimension of the mass operator, and β(g) is the βfunction of QCD. At the leading order β(g) = ?β0g3/(4π)2 and β0 = (11 ?2nf /3), where nf is the number of ?avors. The second term is called the trace anomaly and is generated in the process of renormalization.
According to the above, the mass of the nucleon is
M=
P
d3x T 00(0, x)P P P
≡ T 00 ,
(6)
in the nucleon’s rest frame. Although I formally work with the matrix elements of the nucleon, it actually is the di?erence of the nucleon matrix elements and the vacuum matrix elements that enters all the formula (the vacuum has zero measurable energy density). According to (6), a mass separation can be found through a decomposition of T ?ν into various parts, which are then evaluated with the deepinelastic momentum sum rule and the scalar charge of the nucleon. [Note that the parts of the energymomentum tensor are not separately conserved, so the breaking of the nucleon energy cannot be Lorentz covariant.]
First of all, let me decompose the T ?ν into traceless and trace parts,
T ?ν = T??ν + T??ν ,
(7)
where T??ν is traceless. According to Eq. (4), I have,
P T??νP = (P ?P ? ? 1 M 2g?ν)/M ,
(8)
4
P T??νP = 1 g?νM.
(9)
4
Combining Eq. (6) with the above three equations, I get,
T?00 = 3 M,
(10)
4
T?00 = 1 M.
(11)
4
Thus 3/4 of the nucleon mass comes from the traceless part of the energymomentum tensor and 1/4 from the trace part. The magic number 4 is just the spacetime dimension. This decomposition, a bit like the virial theorem, is valid for any bound states in ?eld theory!
The traceless part of the energymomentum tensor can be decomposed into the contribution from the quark and gluon parts,
T??ν = T?q?ν + T?g?ν ,
(12)
where
3
T?q?ν
=
1 ψ?i?D(?γν)ψ 2
?
1 g?νψ?mψ, 4
(13)
T?g?ν
=
1 g?νF 2 4
?
F ?αF να.
(14)
Although the sum of T?q?ν and T?g?ν with bare ?elds and bare couplings is ?nite (now neglecting the trace anomaly), the individual operators are divergent and must be renormalized. Under
renormalization, they mix with each other and with other BRSTexact and the equations
of motion operators, which have vanishing physical matrix elements [2,3]. For my purpose, I regard both operators renormalized and dependent on a renormalization scale ?2. De?ne
their matrix elements in the nucleon state,
P T?q?νP
= a(?2)(P ?P ν ? 1 g?νM 2)/M, 4
(15)
P T?g?νP
= (1 ? a(?2))(P ?P ν ? 1 g?νM 2)/M, 4
(16)
where I have used Eq. (8) to get the second equation. The constant a(?2) is related to the deep inelastic sum rule [6],
1
a(?2) =
x[qf (x, ?2) + q?f (x, ?2)]dx,
(17)
f0
where the sum is over all quark ?avors and qf (x, ?2) and q?f (x, ?2) are quark momentum distributions inside the nucleon in the in?nite momentum frame or lightfront coordinate system. Again, according the Eq. (6), I ?nd the contribution to the nucleon mass,
T?q00
= 3a(?2)M, 4
(18)
T?g00
= 3(1 ? a(?2))M. 4
(19)
Finally, I turn to the trace part of the energymomentum tensor T??ν. According to Eq. (5), I decompose it into T?m?ν and T?a?ν, the mass term and trace anomaly term, respectively. Both operators are ?nite and scale independent. If I de?ne,
b = 4 T?m00 /M,
(20)
then according to Eq. (11), the anomaly part contributes,
T?a00 = 41(1 ? b)M.
(21)
Thus, the energymomentum tensor T ?ν can be separated into four gaugeinvariant parts, T?q?ν, T?g?ν, T?m?ν, and T?a?ν. They contribute, respectively, 3a/4, 3(1 ? a)/4, b/4, and (1 ? b)/4 fractions of the nucleon mass. The corresponding breakdown for the hamiltonian is, HQCD =
Hq′ + Hg + Hm′ + Ha, with
4
Hq′ =
d3x ψ?(?iD · α)ψ + 3 ψ?mψ , 4
(22)
Hg =
d3x 1 (E2 + B2), 2
(23)
Hm′ =
d3x 1 ψ?mψ, 4
(24)
Ha =
d3x
9αs 16π
(E2
? B2).
(25)
where I have consistently neglected γm and the higherorder terms in β(g). One can put them back if a higher precision analysis becomes necessary. I also have taken nf = 3. [Note that the heavy quarks do contribute to the mass term, the kinetic and potential energy term, and the trace anomaly term. However, the contributions cancel each other in the limit of mf → ∞, and thus for simplicity I neglect them.] If I rearrange the mass terms by de?ning,
Hq = d3x ψ?(?iD · α)ψ,
(26)
Hm = d3x ψ?mψ,
(27)
then the QCD hamiltonian becomes,
HQCD = Hq + Hm + Hg + Ha.
(28)
Here Hq (Eq. (26)) represents the quark and antiquark kinetic and potential energies and contributes 3(a ? b)/4 fraction of the nucleon mass. Hm (Eq. (27)) is the quark mass term and contributes b fraction of the mass. Hg (Eq. (23)) is the normal part of the gluon energy and contributes 3(1 ? a)/4 fraction of the mass. Finally, Ha (Eq. (25)) is the gluon energy from the trace anomaly. It contributes (1 ? b)/4 fraction of the mass.
To determine the decomposition numerically, I need the matrix elements a and b. The deepinelastic scattering experiments have determined a(?2) with an accuracy of a few per
cent. Using a recent ?t to the quark distributions [7],
aMS(1GeV2) = 0.55
(29)
where MS refers to the modi?ed minimal subtraction scheme. Without the heavy quarks, the matrix element b is,
bM = P muu?u + mdd?dP + P mss?sP
(30)
The ?rst term is the πN σterm apart from a small isospinviolating contribution of order 2 MeV. A most recent analysis gave a magnitude of 45±5 MeV for this term [8]. So the only unknown in our analysis is the strange scalar charge P mss?sP in the nucleon. There are model calculations for this quantity in the literature [9]. Here I choose to estimate it using two standard approaches, though both of them are not completely satisfactory.
In the ?rst approach, the strange quark mass is considered small in the QCD scale and so the chiral perturbation theory can be used to calculate the SU(3) symmetry breaking e?ects. A recent secondorder analysis on the spectra of the baryon octet combined with the measured σterm yields [8],
5
P s?sP ? 0.11 × P u?u + d?dP
(31)
? 0.77,
(32)
where in the second line, I have used (mu + md)/2 ? 7 MeV at the scale of 1 GeV2 [10]. Taking the strange quark mass to be 150 MeV at the same scale, I get,
bM ? 160 MeV,
(33)
Using Eq. (11), I have,
P

αs π
F
2P
= ?693 MeV.
(34)
In the second approach, the strange quark is considered heavy in the QCD scale. Using heavyquark expansion, it was found [11],
P mQQ?QP
=
?1 12
P  αs F 2P π
.
(35)
Thus the strange quark contribution in
P
ψ?mψ
+
β(g) 2g
F
2P
= M,
(36)
which is an explicit form of Eq. (11), cancels. From the above equation and the σterm, I ?nd,
P

αs π
F
2P
= ?740 MeV.
(37)
This yields a strange matrix element P mss?sP = 62 MeV. Together with the σterm, I determine,
bM = 107 MeV.
(38)
The complete result of the mass decomposition at the scale of ?2 = 1GeV2, together with the two numerical estimates, is shown Table 1. I have not shown the errors due to omission of higherorder perturbative e?ects and errors on the σterm and current quark masses. The total e?ect on individual numbers is about 5 to 10 MeV. Thus I have rounded up the numbers to nearest 10 MeV. The largest uncertainty is from the matrix element P mss?sP , which could be larger than the di?erence of the two estimates shown. Nevertheless, I will argue below that the total strange contribution to the nucleon mass is quite small and with a smaller uncertainty.
The following comments can be made with regard to the numerical result.
? The quark kinetic and potential energies contribute about 1/3 of the nucleon mass. Because the quark kinetic energy must be very large when con?ned within a radius of 1 fm, there must exists a large cancellation between the kinetic and potential energies. This may not be entirely surprising in the presence of strong interactions between quarks and gluons. Such strong interactions are clearly at the origin of the chiral symmetry breaking, modeled, for instance, by the NambuJonaLasinio [9].
6
? The decomposition of the quark energy into di?erent ?avors is possible. Taking the number 270 MeV (the ms → 0 limit) as an example, I ?nd the upquark energy in the proton is 250 MeV using the momentum fraction carried by up quark 0.375 [7], the down quark energy, 105 MeV, and strange quark energy, ?85 MeV. Further decomposition into valence and sea contribution cannot be made without knowledge of the separate valence and sea contributions to the scalar charge.
? The quark mass term accounts for about 1/8 of the nucleon mass. About half of which or more is carried by the strange quark. The contributions from the up and down quarks are well determined by the σterm.
? The normal gluon energy is about 1/3 of the nucleon mass and the trace anomaly part contributes about 1/4. From these two, I deduce the colorelectric and colormagnetic ?elds in the nucleon separately (take αs(1GeV)? 0.4),
P E2P = 1700 MeV,
(39)
P B2P = ?1050 MeV.
(40)
So the magnetic?eld energy is negative in the nucleon! This of course is due to a cancellation between the quark’s magnetic ?eld and that of the vacuum. The electric ?eld in the vacuum is presumably small, however, it is large and positive in the nucleon. This behavior of the color ?elds in presence of quarks is very interesting, it may help to unravel the structure of the QCD vacuum.
? In the chiral limit, the gluon energy from the trace anomaly (M/4) corresponds exactly to the vacuum energy in the MIT bag model [12]. The role of such energy in the model is to con?ne quarks. Thus we see here a clear way through which the scale symmetry breaking leads to quark con?nement. To keep the con?nement mechanism, a model must include Ha, the anomaly part of the hamiltonian.
? The strange quark contributes about ?60 MeV through trace anomaly. When adding together with the kinetic and potential energy contribution ?85 MeV and the mass term 115 MeV (the ms → 0 limit) the total strange contribution to the nucleon mass is a mere ?30 MeV. (The other limit gives a total of ?45 MeV.) The smallness of the contribution is, to a large extent, insensitive to the matrix element P mss?sP .
To summarize, we have ?nd a separation of the nucleon mass into contributions from the quark kinetic and potential energy, gluon energy, and the trace anomaly. The largest uncertainty is from the strange matrix element P mss?sP . The result has interesting implications on the quarkgluon structure of the nucleon and on the response of the QCD vacuum to color charges.
I wish to thank D. Freedman, F. Low, and K. Johnson for useful discussions and suggestions.
7
REFERENCES
[1] C. Rebbi, Lattice Gauge Theories and Monte Carlo Simulations, Addison Wesley, Reading, MA, 1988.
[2] J. C. Collins, Renormalization, Cambridge Univ. Press, Cambridge, 1984. [3] S. D. Joklekar and B. W. Lee, Ann. Phys. 97 (1976) 160. [4] R. L. Ja?e and A. Manohar, Nucl. Phys. B337 (1990) 509. [5] R. T. Crewther, Phys. Rev. Lett. 28 (1972) 1421; M. S. Chanowitz and J. Ellis, Phys.
Lett. B40 (1972) 397; J. C. Collins, A. Duncan, and S. D. Joglekar, Phys. Rev. S16 (1977) 438; N. K. Nielsen, Nucl. Phys. B120 (1977) 212. [6] F. J. Yndurain, Quantum Chromodynamics, SpringVerlag, 1983. [7] The CTEQ2 distribution, by J. Botts, J. Huston, H. L. Lai, J. G. Mor?n, J. F. Owens, J. Qiu, W. K. Tung, and H. Weerts, 1994. [8] J. Gasser, H. Leutwyler, and M. E. Sainio, Phys. Lett. B253 (1991) 252. [9] T. Hatsuda and T. Kunihiro, UTHEP270, RYUTHP 941, 1994. To be published in Phys. Rep. [10] S. Narison, Riv. Nuovo Cimento, 10 (1987) No. 2. [11] A. I. Vainshtein, V. I. Zakharov, and M. A. Shifman, JETP Lett. 22 (1975) 55; E. Witten, Nucl. Phys. B104 (1976) 445. [12] A. Chodos, R. L. Ja?e, K. Johnson, C. B. Thorn, V. F. Weisskopf, Phys. Rev. D9 (1974) 3471.
8
TABLES
TABLE I. A decomposition of the nucleon mass into di?erent contributions. The matrix elements a and b are de?ned in Eqs. (15) and (20).
mass type quark energy quark mass gluon energy trace anomaly
Hi ψ?(?iD · α)ψ
ψ?mψ
1 2
(E2
+
B2)
9αs 16π
(E2
?
B2)
Mi 3(a ? b)/4
b 3(1 ? a)/4 (1 ? b)/4
ms → 0(MeV) 270 160 320 190
ms → ∞(MeV) 300 110 320 210
9