Volume 158B, number 5
PHYSICS LETTERS
29 August 1985
G O L D S T O N E B O S O N S AND SCALAR G L U O N I U M Harald G O M M and Joseph S C H E C H T E R
PhysicsDepartment,SyracuseUniversity,Syracuse,NY13210,~USA Received 26 December 1984; revised manuscript received 15 March 1985
We study the model of Goldstone bosons interacting with a scalar gluonium field in such a way that the trace anomaly equation is automatically satisfied. It is argued that the mixing of gluonium and quarkonium fields plays an important role. A possible decoupling of the Goldstone fields from a physical scalar-glueballis argued to be reasonable only if the glueball'smass is extremely low.
The properties of a scalar glueball (see ref. [ 1] for an up to date review) are of great interest for understanding QCD. Recently the possibility of relating the properties of such a particle to the trace (or dilation) anomaly equation has been explored in simple effective lagrangian models [ 2 - 5 ] . Both the case [2,3] of approximating pure (no "matter") QCD by a single scalar field and the case [2,5] where spin 0 matter fields belonging to the (3, 3*) representation of chiral SU(3) were present were treated. The former case was described by a unique lagrangian (up to terms with two derivatives) while the latter was noted to permit an intrinsic choice as to whether the field postulated to dominate the trace of the energy momentum tensor Ouu at low energies corresponds to a physical glueball or gets eliminated in terms of singlet scalar quarkonium (a situation which appears [6] to hold for the 7?' particle in the pseudoscalar gluonium channel). Very recently, l_Anik [5] has claimed that the models ofref. [2] do not satisfy a low energy theorem [7] on the matrix element of Ouu between pion states. He constructs a model designed to satisfy this theorem with the unusual result that the (zero mass) pions decouple from the glueball. If this result could be shown to be true and model independent it would clearly be spectacular, implying possibly very large SU(3) breaking in scalar gluonium decay. This has led us to examine the models of ref. [2] in more detail. We find that those models do in fact satisfy the low energy theorem of ref. [7], although it may not 0370-2693/85/$ 03.30 © Elsevier Science PUblishers B.V. (North-Holland Physics Publishing Division)
be obvious at first sight. The model of l_Anik is seen to correspond to a very special case of a model in ref. [2] which only seems reasonable physically if the scalar glueball is (roughly) of the order of (or smaller than) the pion mass. Otherwise scalar gluonium is expected to decay very quickly into two pseudoscalars. An interesting feature which emerges from our analysis is that if a physical glueball dominates Ouu at low energies, an upper bound may be set for the lightest particle in the scalar channel. The precise value of this bound depends on the numerical value of the QCD gluon condensate. In our treatment we use the linear rather than the (often more convenient) non-linear chiral realization to describe the spin zero quarkonium fields. The physics of the non-linear model corresponds to taking the physical scalar quarkonium mass to infinity. While this is a good approximation at energies comparable to the pion mass it is clearly very dubious when dealing with a scalar gluonium state which is generally expected to have a mass comparable to that of scalar quarkonium. For orientation first consider the case when no quark matter is present. With the convenient normalization of a scalar gluonium field H = - [[3(g)/g] X Tr(Fuv Fur), Fur being the QCD gauge field and fl(g) the renormalization group function, we want to find a model which satisfies the anomaly equation [8]
~ D = Ouu =H.
(1) 449
Volume 158B, number 5
PHYSICS LETTERS
Here Du is the dilation current. The appropriate lagrangian is [2]
£ = -½ aH-3/2(a ° H ) 2 - } H l n (H/A4),
(2)
where A is a constant of mass dimension one and a is a dimensionless constant. The dimension one gluebaU field h is defined in terms of the fluctuation around the minimum ((/-/) = A4/e) as H = (H) + Zh. The quantity Z is given in terms of the glueball mass m b by Z = 2rnh(H)l/2 while a = (H)1/2[4m2h . (H) is - 4 times the vacuum energy density and is thus given by 0.00177 GeV 4 according to the naive interpretation of the bag model [9] or 0.0135 to 0.034 GeV 4 according to the QCD sum rule [10] approach. The model may be presented in a number of equivalent ways; for example the substitution H = (H) eX converts it to the model of ref. [3]. Now consider the situation [2] when scalar and pseudoscalar mesons are added [2]. These belong to the matrix M which transforms as M-+ ULMU~R under left unitary (UL) and right unitary (UR) transformations. We will discuss two possibilities: (model A):
£A = --½ Tr(ao m ~ m t )
V2(M,H),
(3)
f'B = £A - ½aH-3/2(~ H)2"
(4)
-
(Model B):
VI(M,H )
-
The "potential" terms V1 and V2 are assumed to contain no derivatives. VI(M, H) is both chiral and scale invariant. For example, with the conventional assignment of dimension 1 to M,
V 1 = c 1 H 1/2 Tr(MM t ) + c 2 Tr [(MM*) 2 ] + c 3 H - 1/2 Tr [(MM "~)3 ] + ....
(5)
V2(M, H) is chiral invariant and is constructed in such a way [2] that eq. (1) holds. (For examlile the second term of eq. (2) by itseld is a possible choice.) We are neglecting quark masses and also possible mixing with ~q?qq fields• We assume that the chiral symmetry is spontaneously broken. This puts an important implicit restriction on the total potential V = V1 + V2. If V were scale invariant, chiral symmetry could not be spontaneously broken since there would be no way for the pion decay constant Frr to be computed out of dimensionless quantities. In the present model F~r is some multiple of A defined in eq. (2). Model A 450
29 August 1985
contains no "kinetic" term for H, which then gets eliminated in terms of scalar singlet quarkonium by its equation of motion. In this model the tightest important state for the operator T r ( F ~ Fur ) is of 7:lq type. We will work here with a simplified " t o y " theory in which the only matter fields present are the pion and a scalar singlet o. (The more general results will be indicated at the end.) Thus we expand M = (o + ixdp)/x/~. It is actually unnecessary to specify the form of V. We need just the chiral symmetry relation
(~ V/ao)¢ k - ( ~ V/ag~k ) o = 0,
(6)
and the scale anomaly equation
H=-4V
+Op.~V/adp+ oaV/ao + 4 H a V / ~ H .
(7)
First, it is interesting to demonstrate that the models here do satisfy the low energy theorem [7] for zero mass pions: (4P0 p0)l/2@+(p,)10uuln+(p))
= - 2 p ' p ' + higher order terms.
(8)
Eq. (8) holds because the general form of the matrix • t t .[. t element of0ov ]s ot[(p - p )o (p - p )v 2 p ' p 6or ] .1. . • • f3(n + p t /~O ~(n' - + p t J~V , while the normalization xt- - i f d3x04# =Po gives fl = 1[2 + O[(p'p')21 and the fact that 0or is a chiral singlet implies that the matrix element should vanish when p ' ~ 0 so a = -/L We see that the derivation of (8) just uses chiral symmetry and is independent of QCD; it is also irrelevant which version of the energy momentum tensor is taken. Now consider model A in eq. (3). We can assume that H has been eliminated so we are simply left with a chiral (but not scale) symmetric theory constructed out of ~and o. It is trivial to see that the trace of the canonical energy momentum tensor satisfies (8), but here we should employ the "new4mproved" tensor [11 ] since it is the one for which ao D o = 0uu. The trace of this tensor is easily found to be
Oou = - 4 V + #O"[]# + o [] o.
(9)
The first two terms on the right make no contribution between massless pion states. All is not lost since the third term ~ (o)D o = (F,r[v~)r-1 o contributes through a o-pole. One needs the mrrr coupling constant which is obtained by differentiating (6) with respect to both o and ~ and evaluating it at the mini-
Volume 158B, number 5
PHYSCIS LETTERS
[see e q. (2)]. Diagonalization of (13) yields
mum:
gcr~rTr= (03 V/O°O1r+O~r-) = v~m2a I F ,
(10)
where m 17 2 = (02 V/Oo 2) is the o-mass and F~r "~ 135 MeV. Using (10) then gives for the matrix element of (9) between pion states
m2(p _ p,)2/[m2 + (lO _ p,)2] = - 2 p ' p ' + higher order terms, in agreement with eq. (8). In the case of model B, we have verified eq. (8) by taking an additional term 4H OV/OH in (9) into account. However we shall not give the details. It should be remarked that in model A the o ~article has the tremendous width P(o ~ rrrr) ~ 3m~/ 16zrF2 (about 1.7 GeV for m o ~ 0.8 GeV) characteristic of the old sigma model and not [1 ] qualitatively dissimilar to the experimental situation. This was discussed many years ago by Crewther [7] and by Ellis [7]. Now let us discuss model B [eq. (4)] which turns out to have a fascinating structure. Both a scalar matter field o and a scalar glue field h are present. These are required to mix, giving physical fields Op and hp : o~
(
cos0
(h/= _sin0
sin0)
.
COS0 / ( ~ p ) P
(11)
The structure of this mixing follows from the scale anomaly equation (7); differentiating in one case with respect to H and in another case with respect to o yields the two equations:
1 = ((o)/Z) (0 2 V/OoOh) + (4(H)/Z 2) (0 2 V/0h2), 0 = (o) (0 2 V/Oo 2) + (4(H)/Z) (0 2 V/OhOo).
(12)
This determines the o - h mass squared mixing matrix to be
A
-B
}, (13)
-B
29 August 1985
C + B2/A
/
where A = (02 V]Oo2), B = - ( 0 2 V/OoOh) = ZF~r A~ 4X/r2(H) and C = Z2/4(I-1). Notice that if A, the "bare" scalar quarkonium squared mass (very roughly m 2) and (/-/) the strength of the gluon condensate are considered known, (13) depends only on the parameter Z which is proportional to a "bare" gluonium mass
tan 20 = 2B/(A - C - B2/A), m2(op,hp) = ~(A + C+B2[A + 2B[sin 20).
(14)
Differentiating (7) with respect to ¢ twice shows that the bare gluon field h must have a coupling to two pions, related to goTr~rin (I 0) by gh~r. = -ZFTrgo"/ 4X/~(H). The coupling constants of the physical particles are then
gap.Zr = ga.Tr [cos0 + (B/A )sinO ] ,
ghp.. =go.. [sin0 - (BIA )cos O] .
(15)
It is interesting to consider the formal limit of our expression when the bare o squared mass,A, gets very large (non-linear o-model limit). Then m2(Op) goes to infinity as A. (1 + B2/A 2), m2(hp) ~ C/(1 + B2/A 2) and sin 0 ~ B/(A 2 + B2)1/2. Note that gl~,r~r -* 0; the scalar field of finite mass decouples fronfthe Goldstone bosons. At first this appears to be a mechanism for hiding scalar gluonium. However it does not seem physically realistic to consider the bare quarkonium mass (around 1 GeV) to be very large compared to the gluonium mass. Nevertheless one would like to be certain that there is no scalar glueball lurking at extremely low energies (say several MeV up to around 2m~r). It seems that its non-observation in the decay process K + -* rr+hp is a possible strong argument against it. We can obtain an upper bound on the lighter of the two scalar masses in model B. This is most readily gotten by allowing Z -* ~ in the expression m2(hp) C/(1 + B2/A 2) just obtained whenA -* oo: mlighter < 2x/~(H)I/2/F.
(16)
The LHS is written as mlighter since the same bound holds even if m(hp) > m(Op). It is amusing that this bound depends only on the glue and matter condensates in the effective models and is independent of the bare masses. In the present "toy model" this bound is not restrictive if the QCD sum rule determinations of 0rt) are used. If, for definiteness we fixA = 1 GeV 2 and (H) = 0.0135 GeV 4 we find that the heavier scalar particle will always have a width greater than its mass when its mass is less than 1.6 GeV. Corresponding to heavier particle masses of 1.02 GeV, 1.12 and 1.60 GeV the lighter particle masses (widths) are predicted to be 451
Volume 158B, number 5
PHYSICS LETTERS
0.44 (0.013), 0.74 (0.48) and 0.88 (1.98) GeV. It seems instructive to compare our model with ref. [5] by making the non-linear substitution M = ½F~ U, U = exp (2i$/F~) in eqs. (3) and (4). The kinetic term becomes the standard form of the nonlinear model - (~~ F~) 2 Tr(a~ Ua~ U t ) and the potential term collapses drastically. For example (5) becomes (noting U U t = 1), V1
=1 c F 2 H 1 / 2 + t__c F 6 H - 1 / 2
41
Ir
643
7r
+ ....
In other words, the scale invariance of V1 is lost. As discussed previously, this is absolutely necessary in order that chiral symmetry break spontaneously. Fur. thermore it is trivially seen that the Goldstone fields decouple from the glue field H in this limit. Thus the model of ref. [5] is just a (probably unrealistic) special case of the earlier model in ref. [2]. After the work o f this paper was completed we received a report by Ellis and L~lnik [12] which presents a somewhat different model along the lines of ref. [5]. They also find a large width for heavy scalar gluonium to decay into two pions. Their model is however different from the present one. It uses the non4inear rather than the linear sigma model and does not include the effects of glue-matter mixing. An interesting aspect of the present model (with a linearly transforming matter multiplet) is that, as briefly noted after eq. (5), it displays a scenario for the way in which, according to the usual assumption, QCD makes the dimensional parameter F,r emerge as a function of a characteristic scale [say A in (2)]. The effective lagrangian which enforces the trace anomaly is not only a measure of confinement (since it gives a non-perturbative vacuum of lower energy than the perturbative one) but, because of glue-matter mixing, also provides a scale so that (o> = F~/x/~ 4= 0 can emerge as the minimum of a chiral invariant potential. In the present note we have concentr~ited on some conceptual aspects of the gluonium-Goldstone boson system using a " t o y " model which exhibits a number of interesting and characteristic features. Elsewhere we will study a more realistic model and compare with experiment. Here we would just like to remark that the mass mixing matrix which results when one includes three flavors and assigns scale dimension three in a consistent way to the Ftq mesons is identical to (13) apart from numerical factors. This lowers the upper bound on the fighter mass to 452
29 August 1985
m2< < y 16 ~ -(=H()6 6 3
MeV) 2,
(H) = 0.0135 GeV 4,
9r
= (1048 MeV) 2 ,
(H) = 0.0338 GeV 4.
This bound may be of relevance. We would like to thank John Donoghue for some helpful remarks. This work was supported by the US Department of Energy under Contract No. DE-AC0276ER03533. References
[1] S.R. Sharpe, M.R. Pennington and R.L. Jaffe, Phys. Rev. D30 (1984) 1013. [2] J. Scheehter, Phys. Rev. D21 (1980) 3393; A. Salomone, J. Schechter and T. Tudron, Phys. Rev. D23 (1981) 1143. [3] A.A. Migdal and M.A. Shifman, Phys. Lett. l14B (1982) 445. [4] J.M. Cornwall and A. Soni, Phys. Rev. D7 (1984) 1424. [5] J. Lfinik, Phys. Lett. 144B (1984) 439. [61 C. Rosenzweig, J. Schechter and C.G. Trahern, Phys. Rev. D21 (1980) 3388; P. Di Vecchia and G. Veneziano, Nucl. Phys. B171 (1981) 253; P. Nath and R. Arnowitt, Phys. Rev. D23 (1981) 473; E. Witten, Ann. Phys. (NY) 128 (1980) 363; A. Aurilia, Y. Takahashi and P. Townsend, Phys. Lett. 95B (1980) 265. [7] M. GeU-Mann, Proc. 3rd Hawaii Topical Conf. in Particle physics, 1969 (Western Periodicals, CA, 1969); J. Ellis, Nucl. Phys. B22 (1970) 478; R.J. Crewther, Phys. Lett. 33B (1970) 305; R. Jackiw, Phys. Rev. D3 (1971) 1347, 1356; M. Voloshin and V. Zakharov, Phys. Rev. Lett. 45 (1980) 688; V.A. Novikov and M.A. Shifman, Z. Phys. C8 (1981) 43. [8] N.K. Nielsen, NucL Phys. B210 (1977) 212; J.C. Collins, A. Duncan and S.D. Joglekar, Phys. Rev. D16 (1977) 438; see also S.L. Adler, J.C. Collins and A. Duncan, Phys. Rev. D15 (1977) 1712; M.S. Chanowitz and J. Ellis, Phys. Rev. 7 (1973) 2490; R.J. Crewther, Phys. Rev. Lett: 28 (1972) 1421. [9] T. De Grand, R.L. Jaffe, K. Johnson and J. Kiskis, Phys. Rev. D12 (1975) 2060. [10] M.A. Shifman, A.I. Vainshtain and V.I. Zakharov, Nucl. Phys. B147 (1979) 385,448. [11] C.G. CaUan, S. Coleman and R. Jackiw, Ann. Phys. 59 (1970) 42. [12] J. Ellis and J. Lfirdk, CERN report TH 4030/84 (October 1984); see also J. Ellis, Nuel. Phys. B22 (1970) 478.