Physics Letters B 318 (1993) 438-446 North-Holland
PHYSICS LETTERS B
On the calculation of effective actions by string methods Michael G. Schmidt and Christian Schubert 1,2 Institut fiir Theoretische Physik, Universitiit Heidelberg, Philosophenweg 16, D-69120 Heidelberg, FRG Received 20 August 1993
Editor: P.V. Landshoff Strassler's formulation of the string-derived Bern-Kosower formalism is reconsidered with particular emphasis on effective actions and form factors. Two- and three-point form factors in the nonabelian effective action are calculated and compared with those obtained in the heat kernel approach of Barvinsky, Vilkovisky et al. We discuss the FockSchwinger gauge and propose a manifestly covariant calculational scheme for one-loop effective actions in gauge theory.
One of the main differences between string theory and particle theory is the fact that in string theory the full (on-shell) S-matrix can be calculated using the Polyakov path integral, i.e. in first quantization. As string theory reduces to particle theory in the limit of infinite string tension, the same should be possible in principle also for the panicle S-matrix. And indeed, by a painstaking analysis of this limit, Bern and Kosower [ 1 ] have been able to derive from string theory a new set of Feynman rules for the calculation of one-loop amplitudes in conventional field theory which, though quite different from the conventional ones, appear to be completely equivalent [2 ]. Those rules are particularly well-suited for calculations in gauge theories, as they combine contributions of different Feynman diagrams into gauge invariant structures [3]. Strassler [4] recently succeeded in deriving essentially the same set of rules without explicit reference to string theory, using one-dimensional path integrals with Green functions adapted to the circle. Those path integrals are the particle theory analogue of the genus one Polyakov path integral. He further demonstrated the advantages of the new formalism for the calculation of both abelian and non-abelian effective actions [5]. In this paper, we will somewhat reformulate this method and apply it to the computation of form factors, i.e. to the organization of the higher derivative terms appearing in the effective action. It will turn out that the twoand three- point form factors for nonabelian gauge theory may be computed quite efficiently. In the three-point case, the calculation will be non-covariant. A manifestly covariant method of calculation will be outlined and applied to the two-point case. Let us begin with Strassler's [4] worldline path integral expression for the one-loop effective action o f a Dirac spinor of mass m minimally coupled to an (abelian or non-abelian) background gauge field A: oo
o T
x.N'Wrfl~xV¢exp[-f dr((1)dc2+½¢~+igAu~U-½igg~,UFu~,~)].
(1)
o
1 Partially supported by funds provided by the Deutsche Forschungsgemeinschaft. 2 E-mail address:
[email protected].
438
0370-2693/93/$ 06.00 © 1993 Elsevier Science Publishers B.V. All rights reserved.
Volume 318, number 3
PHYSICS LETTERSB
9 December 1993
In this construction, the x~(z)'s are the periodic functions from the circle with circumference T into Ddimensional Euclidean spacetime, and the ¥~ (z)'s their antiperiodic supersymmetric partners, obeying
{~,~}
(2)
= g"~.
JV is the path integral normalization factor, T
Ac-l = Trf Z~x~ exp[-/ dT((~)~c2 + ½~')],
(3)
0
and 8 the worldline metric; in the following we set 8 = 2. Note that our normalization includes the spinor degrees of freedom, which makes for a factor of 4 compared to [4,5]. In the nonabelian case path ordering is implied. Insertions into the free path integral will be contracted using the one-dimensional Green functions adapted to the (anti-)periodicity conditions, (XI~(.C1)XV(.~2)) =
(Vv(Zt)V"(Z2)) =
_gaVGB(zt, Z2) ½g~Gr(zl, z2)
('el ~T2)2],
= _gVV [[ Zt_.C2 I = ½ga"sign(zl- z2).
(4)
As a warm-up, let us reconsider the classical case of the l-loop effective action for spinor electrodynamics in a constant abelian background [6]. Here choosing an appropriate gauge one may write Av = ½xPFpv, the full action remains quadratic, and the calculation of F [A ] can be reduced to a calculation of determinants as usual. However, in the present formalism those determinants involve just integral operators on the circle: Denoting determinants by ] I and the derivative with respect to z by 0, all we have to calculate is oo
T
F[A]=-2f-~-[4nT]-n/2e-m2rj~/Dxlg~/exp[-/dz(lyc2+ 0
½g~b+ ½igx~Fg.k"-igggF~.g~)]
0
I :1
O0
I- 102 + ½igF @ 0 [
JTo
i oim' 1/2 [ ½ 0 -
igF[ -1/2
2igF ~ 002-11-1/211 - 2igF ~ tg-I I1/2
=
o oO
- 2 [ ~-[4nT]-n/Ee-m2r]I - igF ® [GB ]1-1/2]I - igF ® [Ge ][1/2 t#
0 OO
= -2
O0
[4nT]
~
v^r'| -IZ.., - - n
0
~
(-ig)nTnTrFn(Tr[GB]n -Tr[GF]n) "
(5)
n=l
Here we used the Log Det = Trace Log identity in the last line, and with [G] we mean the integral operator associated to the Green function G. The T-dependence of those integral operators has been extracted by a rescaling r ~ Tr, which leaves us with operators on the unit circle. Dots on a G will always signify a derivative with respect to the first variable; in particular, 0n(zl, z2) = sign(rl - z2)
2 ( T I -- r 2 )
T
0B(Zl,'~2) = 2~(Zl-
Z2)
2 T"
(6) 439
Volume 318, number 3
PHYSICSLETTERSB
9 December 1993
The operator traces may be computed using the Fourier bases {exp[2nikz][k E Z} for the periodic and {exp[2ni(k + ½)z] Ik E Z} for the antiperiodic functions. The result is T r [ 0 n ] 2n =
2(-1)n((2n)~ -2n, Tr[GF] 2n =
2(2 2n-
1)(-1)n((2n)~z -2n,
(7)
in terms of the Riemann (-function (traces of odd powers vanish). Finally oo
oo ~-(gT)Znn
F[A]=--2f-~[47rT]-O/2e-m2Texp[--n~ 0
o¢
F[A]
--
1 )((2n)zr-2"TrF2n].
(8)
=
This is the generating functional for eq. (3.23) in [ 5 ] #1 (see e.g. [8])
!
(2 2n-1
/ dSe-i~,.~r ~ .cosh(gas)cos(gbs) s
= 8n2j
[g a°sinh(gas)sin(gbs )
,2.
It may be transformed into the standard representation
I]
(9)
~2 ,
0
where b 2 = ½ [ - ( E Z - B 2)+ by diagonalizing the matrix
~/(E2-Bz) 2+ 4(E. B) 2],
(1o)
Fu~, showing
T r F 2n = 2[(a2) n + ( - b 2 ) " ]
(11)
and using the Taylor expansions of In [x cot (x) ], In [x coth (x) ]. After this amusing calculation, let us proceed to our main objective, which is the calculation of form factors in the nonabelian theory. That means we want to calculate (taking the case of a (complex) scalar loop first)
F[A] =
[4nT]-~/2e-"2rXTrfVxexp - f dz(~.~2 + igA,Y') 0
(12)
0
up to a fixed power in the background field, but to all orders in the covariant derivatives Du = 0~ - igA u of the field (for the history of this subject and for possible physical applications see ref. [9] and references therein). To begin with, we introduce the loop center of mass x0, writing
xU(z) = x~ + yU(r),
(13)
with T
f
dzyU(z) = 0,
(14)
0
~l We believe that the recursion relations derived in [5] for the above traces should be related to the well-known recursion formula for the Bernoulli numbers. #2 Analogous methods have been already used by Metsaev and Tseytlin [7] to obtain a generalized Schwinger formula for the bosonic string and, by taking the infinite string tension limit of this generalization, for pure YM theory. 440
Volume 318, number 3
9 December 1993
PHYSICSLETTERSB
and extract the integral over the center of mass from the path integral:
(15) Next we Taylor expand the interaction part ~UA. (x (z)) with respect to
xo, use ~u
JcUAu (x) = yUeV°Au (Xo)
= ))u to write (16)
and expand the path-ordered interaction exponential to get cx~
oo
"rl = T
F[A] =Trf-~-[nnT]-O/2e-m2r/dOx°Z (-ig)nTn 0
n=0
/
z2
"~n- I
dz2fdx3.., f dzn
0
0
0 T
× . A f f " D y y ' u l ( T l , O. Y .( q ). ° o. ) a ( lm )(xo)
~"Unln'"~Y(tn)O(n)A(n) -.-~, ~n,,~u, (x0)exp [ - f
dzly "2] .
(17)
0
We have labeled the background fields Aut,..., Au., and the first z integration has been eliminated by using the freedom of choosing the point 0 somewhere on the loop. Instead of restricting the colour trace by fixing the cyclic order, we have introduced an explicit factor of l/n. The single terms in this expansion may now be computed by Wick contractions in the one-dimensional worldline field theory, using the above bosonic Green function on the circle and formulas familiar from string theory,
(ey(q)ot,)ev(~2)oc2))= e-an(~,,~)oc,)o(2), For the term quadratic in
(j~u (~'1) ey(~'2)o(2)) = --GB ('t'l, T2)0~2 ) •
A u, this leads to a single
(18)
z integration,
T T f
t':, t,'r _ x2~/~l ~#21.-GB(T,'t2)O(I)O(2).4(I)A(2) d'r2[Gt)(T, z2)g u~u2_ 'O'Bl.l,t2! tJ(2)tJ(i)Jw -"a/zI -'a/~2 ,
(19)
0
which using the scaling properties of Gn, Gn, Go and partial integration with respect to z2 may be transformed into 1 T2f
~ltl ~.u2 ~--TGB(I,u2)O(I)O(2)A(1)A(2) du2GB (1, u2)2 (glqt,(20(l)O(2) _ ~(2)~(I) . . . . Ut ~-u2 ,
(20)
0
where Gn now refers to the unit circle, Gn ( 1, u2) = u2 ( 1 - u2). But ,q#1,0#2 x A(I)A(2) (glqlz20(l)O(2) -- ~(2)~(1) )za/~1 -ca/t2
is just the "abelian" part of 1Fu(l)F(2)u~,and we need not see more of this tensor to determine its form factor: Anticipating covariantization and performing further partial integrations both with respect to u2 (to get rid of the G2) and to Xo (yielding 0<1)0(2) = -0(2) 2) to get an expression standardized in the sense of [9], we can already predict the following result for the full two-point one-loop effective action: oo
~[A]=
Tr/-~-[4~zT]-O/2e-mZTT2(-g2)/dDxoFu.(xo)F~Cal(,)Ffl'(Xo),
(21)
0
441
Volume 318, number 3
PHYSICSLETTERSB
9 December 1993
with 1
F~(~)=
E
- ~
-
,
f(~) =
/
due -"~-~)~,
~=-TD~D
~.
(22)
0
This is the same expression as has been reached by [9 ] using covariant perturbation theory. Before proceeding to the three-point case, let us say something about the systematics of the calculation in general. Every 17(zi) has to be contracted, either with another 17(zk), or with an eY(~)°"L The most basic case is if there are only contractions of the second kind, as the result will be just one or a product of several closed chains of GB' s,
aB(Til,Ti2)aB(Ti2, ri3)...aB(Tin, r il)Cl x"lt'2 1 (il). ..O.#i(in),
(23)
multiplied by a universal factor of
exp[--i~
(24)
from the contraction of the exponentials among themselves. Contractions of the first kind would result in Ga's, which might be reverted into GB's by partial integrations later (this problem has already been discussed at length in refs. [1-5 ] ). From the systematic point of view, however, we prefer to eliminate the Ga's before they arise: Before contracting a 17(zi) with a 17(Zk), one performs a partial integration in Zk. This leaves a new 1?(rk) after contraction of the 1?(ri) in the main part of the partial integration, and hence the (1?y) chain construction works in this case, too. But there are further terms both from the boundary of the rk integration and from differentiating the Zk-dependent boundary of Tk+~. Altogether we obtain for one contraction step of the first kind a factor of #3
--GB ('fi, "t'k)g lai#k [--1?2(rk)O (2k) "1- J(rk -
rk-1) --J(rk)
-- J ( r k -
rk+l)].
(25)
Whereas the first term of this expression is part of a continuing chain, the other ones may be represented pictorially as the end of a (1?y) chain at Zg (# rio, if contraction has started with 17(rio)). As one can infer already from the two-point calculation, it is useful to always combine the first term with the double contraction
(P(ri)e y~'k)a~)(eyl,,)o.~1?(rk)). The result of one particular total contraction will be a product of both open and closed chains of Gn's, where open chains start freely and may but need not necessarily end on other chains (closed or open). Treating the three-point path integral in this way, we obtain (after the usual rescaling r = Tu) the u-integrals 1
u2
2 1 ~ ) )GBt2 . 2 (O(1)GB13 v" du2 f du3{ [ (gaVO(I)O(2) - 0(2)0(1 + O~2)GB23)+cyclicterms]
f 0
0
_ GB12dB13dB23 [gas (0~1)0(2)0(3) x exp [--GB12C9(1)0(2) --
_ 0~2)0(1)0(3)) ..[. C . t .
Gsl300)O(3 ) -
(1)
~q2 ~/t .qv 2 /z v + u(3)U(l)t/(2 ) -- 0(2)0(3)0(1)] } (2)
GB230(2)tgt3)]A 2 (xo)A~
( x o ) A ( 3 ) ( X o ) -1- boundaryterms,
(26)
where "boundary terms" stands for all terms involving one of the J-functions in the contraction formula (25). Those terms reduce to a single u-integration and may be written as functions of f (all)), f (0~2) ) and f (Oq(3)), 2 #3 J-functions with undefined 442
Zk±l should be skipped.
Volume 318, number 3
PHYSICS LETTERS B
9 December 1993
where f is the basic two-point function defined in eq. (22) (part of them in fact does not contribute to the three-point form factors, but to the covariantization of the two-point form factor). According to the analysis of [9 ], on the three-point level there are - up to powers of covariant derivatives, Bianchi identities and partial integrations in x0 - two independent invariant tensor structures, which may be chosen as
Tl = TrF~uF~F~u,
T2 = TrF~'PDUFu,~D~F~p.
(27)
The terms in the second square bracket of formula (26) combine to produce simply the "abelian" part of -T1, while the terms in the first square bracket contribute to both T~ and T2. Collecting all contributions to T~ we obtain for its coefficient the u-integral 1
u2
f du2f 0
du3exp[-GBI20(1)O(2)-GB130(I)O(3)-GB230(2)O(3)]
0 ×
{(~B12dB13dB23
1
.2 ][GB12(-GBI3
.2 .2 • "1- 6 B 2 3 ) -I- GBI3(GBI2 d- 6 B 2 3 ) d- G B 2 3 ( G B 1 2 - G B 1 3 ) ] } .
(28)
To compare this coefficient with the results of [ 9 ], we again standardize by partially integrating both with respect to u and Xo, introduce Tl~(i)ltl-~(i)
and performe the change of variables 31 = u 2 - u 3 ,
32 = 1 - ( u l - u 3 ) ,
33 = u l - u 2 ,
(29)
leading to OG
Fr,[A]= f -~-[4nT]-O/2e-m2rT3(-ig3) f d°xoFr~(,l,~2,~3)TrFAF£Fy(xo), 0
~al
4
f~ 1 (~1,~2,~3) =
3-~
[27
,2] (d + --) + 12-- + 8-~- F(~1,~2,~3) + boundary terms.
(30)
Here we have defined 27 = ~ 1 + ~ 2 + ~ 3 ,
"=~1~2~3,
A =~2+~22+~2_2~1~2_2~1~3_2~2~3,
(31)
and the basic three-point function F (~1, ~2, ~3) = f d3at~(1
- Oil - or2 -
OL3)exp[--alOt2~3 -- ot20t3~1 --
Ot3Ot1~2] ,
(32)
i /
u~0
originating from the universal contraction of exponentials alone. Again we have assumed that the higher order calculations will produce the missing "nonabelian" parts of the tensor TI, as should be granted by the consistency of the Bem-Kosower formalism. Expression (30) for the form factor Fr~ ~a of Ti now may be easily identified with the one found in [ 9 ] by covariant perturbation theory # 4 . #4 The fact that superparticle path integrals can be an efficient alternative to heat kernel methods is well-known from the calculation of index densities [ 10 ]. 443
Volume 318, number 3
PHYSICS LETTERS B
9 December 1993
Now let us return to the case of the spinor loop, eq. (1). At the two-point level, there is one additional contraction to compute compared to the scalar case, namely J~' (1) 1 ( X o ) ~b¢vl (~/AI ('C I ) e y(zl )0(1) ./AlV
(z,)g/~2('r2)eY(~2)°12)F~222(Xo)~"2(r2)) ,
(33)
giving a u-integral 1
I T 2 f dU2GF( 1, u2 )2 (glqV2gVl#
2 -- o~'m~2~,w2~,,~e- r ° B (uI'u2)O(1)O(2)F(I)/t 1vt F.(2)u2~, 2.
(34)
0
Using Ge ( 1, u2) = 1 and taking the global factor of - 2 in ( 1 ) into account, we obtain the following relationship between the scalar and spinor two-point form factors: F2spin(~)
=
-2[F~(~)-
¼f(~)].
(35)
In the three-point case, we have two additional contraction possibilities: The first one is a contraction of (tl~,/tl a.y(Tl )0(1) /~ (1) , v I j~.,/~2a,Y(r2)o(2)F (2) j~tv 2 .,,,U3eY(:3)ota)F (3) illV3% ,p. ~ .z //lV 1 IV W ~ /~2v2 '/I Y" kt3v3 "t" / ,
(36)
and it simply produces a contribution of -~F (~1, ~2, ~3) to Tt. The second one, ( o / / / I ~ y ( I ' I ) 0 ( I ) iL'(I) ' " v l W ~ z /ZlV l Or*
~¢IZ2eY(~2)°(2)F(22)2~l/IJ2y'lt3eY(T3)°(3)Akt3)
(37)
contributes to Tt and T2, with a contribution to Tt of -½ F (~1, ~2, ~3) and thus cancelling the first one. The final result for the T~-form factor for spinor loops becomes simply F~'!spin ( ~ 1 , ~ 2 , ~ 3 )
= - 2 F T Iseal (~b~2,~3).
(38)
The form factors for T2 will be treated in a more detailed publication [ 11 ], as well as those of the cyclic four-point tensor. Now let us rather sketch briefly how this type of calculations might be done in a manifestly covariant way. For fixed x0 the Fock-Schwinger gauge [ 12 ]
y~Aa~(xo + y ( z ) ) = 0
(39)
is very convenient for writing down manifestly gauge-invariant expressions for the effective action. In this gauge, 1
drlqFflu(x° + qy)yP.
Aa~(xo + Y) = /
(40)
0
Following the proof of [ 12 ], F~a~ can be covariantly Taylor-expanded as a F~(Xo + fly) = (enY'D)abFb~z(Xo),
(41)
leading to 1
A.(xo + y) =
drpl yPe~rDFp~(xo) = ½yPFp. + ~y~ypD.Fp~ + . . . . 0
444
(42)
Volume 318, number 3
PHYSICS LETTERS B
9 December 1993
Expanding again our starting expression (12) for the scalar effective action in powers of g and using the F o c k Schwinger gauge relation (42), we arrive at ~'I=T
r A _ , - r ' l - D / 2 ^ - m 2 T f ...ID~ ~
F[F] =
t..~l~ j
~
j
u .~oz..."
0
/ x
n=0
1
1
d~hrh..,
0
(--lg)
r/ r/ ~
n
"c2
dr2
T 0
"On_ 1
dl3.., 0
drn 0
0
,qn{l)... ODin)
0 T
f
..
exp
[-f 0
]
dr ¼¢2 F ~ t % , . . . f ; ; u , T r ( T a , . . . T a , ) ,
(43)
where the D {i) act on Fpa~i and in this sense commute in the prefactor. (1/rh)O/OD~pi~) creates the y:~ of relation (42). Note that the contraction of a certain y(z~) in e ~y{~)t)") with y ( r k ) ( k J= i) and the commutation of - r ~ ( i ) factors in the exponential are noncommuting operations in the nonabelian case. y (r~)~, The systematics of the contractions are the same as before. However, the contraction of the exponentials now leads to
[E
exp-
~" i
th~lkGB(ri, Zk)D
,,,D ,k,],
(44)
J
which is symbolical in the case of nonabelian gauge theories: Each polynomial in D~i) for fixed i has to be written in all possible orderings, and the resulting expression has to be normalized by the number of possible orderings. This is also the way further factors D {i) have to be handled. In nonabelian gauge theory the ordering of T a matrices and of the various covariant D ~i) acting on the Fu~ in the ith position requires reduction to a set of standard invariants using D-commutators (leading to further F ~ ' s ) and Bianchi identities, an admittedly painful procedure. Formfactors of n fields F ~ in the covariant formalism just described contain one further integration th per field. In the case of two fields F ~ we have checked that the additional integrations can be performed and that one arrives at the a-form of ref. [9]. In the case with three Fu~ this seems to be a more serious problem and is under consideration. For application in Q C D and in the standard electroweak theory of course besides the complex scalar fields we need gauge fields and fermionic fields in the loop and also Higgs fields in the background. In this context the above representation of the fermions by spinning superparticles in the loop [ 13 ] appears to be very economical. However if we want to calculate fluctuation determinants, e.g. for a discussion of radiative corrections to instantons, sphalerons, bounce identities, a different technique seems to have advantages: To calculate the bilinear part of the quantum Lagrangian in a gauge field - Higgs fields background in a convenient quantum gauge (for the background field we can have a separate gauge). The 't H o o f t - F e y n m a n background gauge seems to be most practical - see refs. [14,15 ]. The corresponding set of operators has the form - D 2 (A) + V (~b, F ~ ), with D ( A ) being some collection of covariant derivatives in some gauge group representations and V a general potential matrix built out of background scalar field ~, field strength F,~, and covariant derivatives thereof (in the case of the electroweak theory the main matrix operator depends on 3 × 3 + 4 = 13 fields [ 14] ). Log Det ( - D 2 + V) is then evaluated with the methods discussed above. Proceeding as in the pure gauge case, there will be also factors e y~ui)t)"~ V {i) (~, F ~ ) in the trace. The resulting expressions are lengthy but straightforward to write down and to evaluate for a finite number of V and D, if a convention for a set of independent operators is fixed. This procedure does not have the cancellations of that of ref. [ 14] with plane wave insertions in Tr Log, and it is 445
Volume 318, number 3
PHYSICS LETTERS B
9 December 1993
within the range o f the method to extend their results beyond the third order #5 (or beyond the 6th order in the case without background field [ 15 ] ). We would like to thank M. Reuter and S. Theisen for helpful discussions. C.S. also thanks J. Lauer for discussions and the Deutsche Forschungsgemeinschaft for financial support during part o f this work. #5 W o r k in progress.
References [1] Z. Bern and D.A. Kosower, Phys. Rev. Lett. 66 (1991) 1669; Nucl. Phys. B 379 (1992) 451. [2] Z. Bern and D.C. Dunbar, Nucl. Phys. B 379 (1992) 562. [3] Z. Bern, L. Dixon and D.A. Kosower, SLAC-PUB-5947; Z. Bern, D.C. Dunbar and T. Shimada, Phys. Lett. B 312 (1993) 277; Z. Bern, UCLA-93-TEP-5, and references therein. [4] M.J. Strassler, Nucl. Phys. B 385 (1992) 145. [ 5 ] M.J. Strassler, SLAC-PUB-5978. [6] J. Schwinger, Phys. Rev. 82 (1951) 664. [7] R.R. Metsaev and A.A. Tseytlin, Nucl. Phys. B 298 (1988) 109. [8] C. Itzykson and J. Zuber, Quantum field theory (McGraw-Hill, 1980). [9] A.O. Barvinsky and G.A. Vilkovisky, Nucl. Phys. B 282 (1987) 163; B 333 (1990) 471, 512; G.A. Vilkovisky, CERN-TH 6392/92; A.O. Barvinsky, Yu. V. Gusev, V.V. Zhytnikov and G.A. Vilkovisky, MANITOBA-93-0274. [ 10] L. Alvarez-Gaumd, Commun. Math. Phys. 90 (1983) 161. [ 11 ] M.G. Schmidt and C. Schubert, in preparation. [12] M.A. Shifman, Nucl. Phys. B 173 (1980) 12, and references therein. [13] L. Brink, P. Di Vecchia and P. Howe, Nucl. Phys. B 118 (1977) 76; A.P. Balachandran, P. Salomonson, B. Skagerstam and J. Winnberg, Phys. Rev. D 15 (1977) 2308; F.A. Berezin and M.S. Marinov, Ann. Phys. (NY) 104 (1977) 336; R. Casalbuoni, Phys. Lett. B 62 (1976) 49; A. Barducci, R. Casalbuoni and L. Lusanna, Nuovo Cim. 35A (1976) 377; J.W. van Holten, in: Proc. Sem. on Math. structures in field theories (1986-87), CWI syllabus, Vol. 26 (1990) 109. [14] L. Carson and L. McLerran, Phys. Rev. D 41 (1990) 647. [15] L. Carson, Phys. Rev. D 42 (1990) 2853.
446