Nuclear Physics B 728 (2005) 1–30
Gauge/gravity duality for interactions of spherical membranes in 11-dimensional pp-wave
Hok Kong Lee a, Tristan McLoughlin a, Xinkai Wu b
a California Institute of Technology, Pasadena, CA 91125, USA b Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA Received 3 March 2005; received in revised form 19 July 2005; accepted 5 September 2005 Available online 12 September 2005
Abstract We investigate the gauge/gravity duality in the interaction between two spherical membranes in the 11- dimensional pp-wave background. On the supergravity side, we find the solution to the field equations at locations close to a spherical source membrane, and use it to obtain the light-cone Lagrangian of a spherical probe membrane very close to the source, i.e., with their separation much smaller than their radii. On the gauge theory side, using the BMN matrix model, we compute the one-loop effective potential between two membrane fuzzy spheres. Perfect agreement is found between the two sides. Moreover, the one-loop effective potential we obtain on the gauge theory side is valid beyond the small-separation approximation, giving the full interpolation between interactions of membrane-like objects and that of graviton-like objects. 2005 Elsevier B.V. All rights reserved.
1. Introduction
In this paper we investigate the gauge/gravity duality in the maximally supersymmetric eleven-dimensional pp-wave background. The gauge theory side is represented by a one- dimensional matrix theory first proposed in [1] while the gravity side is represented by eleven- dimensional supergravity. Similar investigations have been performed in flat space for various M-theory objects [3–9]. In particular, we will be interested in the two-body interactions of point like gravitons and spherical M2-branes.
E-mail addresses: [email protected] (H.K. Lee), [email protected] (T. McLoughlin), [email protected] (X. Wu).
0550-3213/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.nuclphysb.2005.09.001 2 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
In pp-wave there is an interesting new connection we can make between a graviton and an M2-brane. Under the influence of the 3-form background whose strength is proportional to a pa- rameter µ, each stable M2-brane curls up into a sphere, with its radius r0 proportional to µ and − its total momentum in the x direction, i.e., we have r0 ∼ µP−. If one takes the limit of µ → 0 while keeping the total momentum P− fixed, then the radius of the sphere goes to zero and we get a point like graviton. If on the other hand, we increase P− simultaneously such that the radius goes to infinity, then the end product will be a flat membrane instead. Thus the gravitons and the flat membranes of flat space can both be regarded as different limits of spherical membranes in pp-wave. One could make another observation by considering two spherical membranes sepa- rated by a distance z.Ifz r0, then one expects the membranes to interact like two point-like gravitons. If z r0 then the interaction should be akin to that between flat membranes. There- fore by computing the interactions between membranes of arbitrary radii and separation, we can take different limits to understand the interactions of both gravitons and membranes. In this paper we will compare the light-cone Lagrangian of supergravity with the effective potential of matrix theory. With a slight abuse in terminology, we will sometimes refer to both as the effective potential. On the supergravity side, we will use a linear approximation when solving the field equations. This means all the metric components computed this way will be 2 proportional to no higher than the first power of κ11, higher order effects such as recoiling and other back reactions could then be neglected. On the matrix theory side, the interaction between two spherical membranes begins at one loop, and for the purpose of comparing to linearized supergravity, only the one-loop effective potential is needed. A more detailed discussion can be found in [2]. Due to the complexity of the field equations of 11D supergravity, we will only compute the effective potential of the supergravity side in the near membrane limit (z r0). The graviton limit (z r0) was already computed in a previous paper [2]. On the matrix theory side, however, it is possible to compute the expression for general z and r0, and by taking the appropriate limits, we are able to find perfect agreement with supergravity. In other words, the matrix theory result provides a smooth interpolations between the near membrane limit and the graviton limit. We also compare our results with the those of Shin and Yoshida [10,11] and where there is overlap we again have perfect agreement. This paper is organized as follows. We begin in Section 2 with a brief discussion of the theories on both sides of the duality. The relations between the parameters on the two sides as well as a discussion on the special limits of interests can be found here. In Section 3 we first define the probe membrane’s light-cone Lagrangian on the supergravity side for a general background, and then consider it in the pp-wave background perturbed by a source. Next in Section 4 the linearized field equations are first diagonalized for an arbitrary static source, and then solved in the special case of a spherical membrane source in the near membrane limit (z r0). The metric and the three-form potential are then used to compute the supergravity light-cone Lagrangian. Section 5 is devoted to computation on the matrix theory side. The one-loop effective potential is found by integrating over all the fluctuating fields. The near membrane limit of the resulting potential is compared with the supergravity light-cone Lagrangian and agreement is found. In Section 6 we compute the one-loop interpolating effective potential on the matrix theory side for arbitrary separation that takes us between the membrane limit and the graviton limit. The result is compared with our earlier work [2] as well as with that of Shin and Yoshida [10,11].Thisis followed by a discussion in Section 7. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 3
2. The two sides of the duality
2.1. The spherical membranes
The nonzero components of the maximally supersymmetric eleven-dimensional pp-wave met- ric and the four-form field strength are given by: 3 9 2 1 i 2 1 a 2 g+− = 1,g++ =−µ x + x ,g= δ , (1) 9 36 AB AB i=1 a=4 F123+ = µ. (2) The index convention throughout this paper is: µ,ν,ρ,... take the values +, −, 1,...,9; A,B,C,... take the values 1,...,9; i,j,k,... take the values 1,...,3; and a,b,c,... take the values 4,...,9. The matrix theory action in such a background is known [1]: 9 1 (M3R)2 9 S = dt Tr D XA 2 + iψT D ψ + XA,XB 2 2R t t 4R A=1 A,B=1 9 1 µ 2 3 µ 2 9 − M3R ψT γ A ψ,XA + − Xi 2 − Xa 2 = 2R 3 = 6 = A 1 i 1 a 4 µ (M3R)µ 3 − i ψT γ ψ − i XiXj Xk , (3) 4 123 3R ij k i,j,k=1 I I where Dt X = ∂t X − i[X0,X ]. The constant M in the action above is the eleven-dimensional Planck mass, and R is the compactification radius in the X− light-like direction in the DLCQ formalism. The supersymmetric configurations of this action have been well studied, see for example [12]. Although in the end we will add a perturbation in the X4 to X9 directions to make it nonsupersymmetric, we begin with the following configuration: µ Xi = Ji,Xa = 0, (4) 3MR3 where [Ji,Jj ]=iij k Jk. If Ji is an N × N irreducible matrix, then it represents a single sphere with radius r0 given by:
1 µ µN r = Tr X2 = N 2 − 1 ≈ . (5) 0 N i 6M3R 6M3R The last approximation is taken with N large, such that the matrix theory correction to the radius (denoted as δr0) is negligible compared to r0. In other words, we choose N large enough that δr0 ∼ 1 → 0. The reader is reminded that the purpose of this paper is to compare matrix theory r0 N2 to the predictions of supergravity, so we are not interested in any finite N effects which are related to matrix theory corrections to supergravity. If Ji is reducible it represents multiple concentric spherical membranes, and the radius of each irreducible component can be found in the same way as for the case above. 4 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
The configuration we are going to use is one where Ji is the direct sum of two irreducible components. This represents two membranes of different radius. In Section 5 we will put in nontrivial Xa to break the supersymmetry, which physically corresponds to placing one of the membranes away from the origin. On the supergravity side, the bosonic part of the Lagrangian of a probe membrane in a general background is 1 L Xµ,∂ Xµ =−T − det(g ) − ij k A ∂ Xµ∂ Xν∂ Xρ , (6) i ij 6 µνρ i j k µ ν with T being the membrane tension, and gij ≡ Gµν∂iX ∂j X being the pullback metric. − One can define the light-cone Lagrangian Llc via a Legendre transformation in the X direc- tion (see Section 3). In the unperturbed pp-wave background, we have 1 A A Π− 2 1 i 2 1 a 2 (L ) = Π−∂ X ∂ X − µ X + X lc pp 2 0 0 2 9 36 2 T A B B A 2 µ i j k − ∂1X ∂2X − ∂1X ∂2X + T ij k ∂1X ∂2X X , (7) 4Π− 3 − with Π− being the momentum density in the X direction. Looking at (Llc)pp, above, we see that one solution to the equations of motion in the unperturbed pp-wave is a spherical membrane i i = 2 a = + = − = = µΠ− at rest with X X r0 , X 0, X t, X 0, with its radius being r0 3T sin θ (note that we 0,1,2 Π− take the worldvolume coordinates σ to be t,θ,φ; also recall that it is sin θ that is a constant on the worldvolume). In terms of the total momentum P− = dθ dφΠ−, this radius is
µP− r = . (8) 0 12πT Later in this paper, we will take the background to be the pp-wave perturbed by a source membrane: Gµν = (Gµν)pp + hµν , and Aµνρ = (Aµνρ)pp + aµνρ . Then Llc = (Llc)pp + δLlc, and we shall compare δLlc with the one-loop effective potential on the matrix theory side (while (Llc)pp agrees with the tree level Lagrangian on the matrix theory side, of course). Identifying the radius calculated on the two sides and P− with N/R, we get the following relation between the tension of the membrane and the Planck mass: 2πT = M3. (9) 5 Another relation we will use to compare the two side is the identification κ2 = 16π [2,13], 11 M9 and for convenience we will define the parameter α = 1 that will appear often on the matrix M3R theory side.
2.2. The effective potential
As already stated, on the supergravity side the quantity of interest is the probe light-cone Lagrangian. On the matrix theory side, the one-loop effective potential is the relevant quantity, and is defined through the Euclideanized effective action with all quadratic fluctuations integrated out. In the following, the subscripts M and E will be used to denote Minkowski and Euclidean signature, respectively. Euclideanization is carried out by defining the following Euclideanized quantities:
τE = iτM , SE =−iSM . (10) H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 5
The Euclideanized action SE will then be expanded about a certain background (to be speci- fied in Section 5) and the quadratic fluctuations integrated out to produce the one-loop effective action Γ . The effective action on the matrix theory side is then defined via the relation: E
ΓE =− dτE Veff. (11)
The minus sign in front of the integral is slightly unconventional, but it is put there for the convenience of comparison with supergravity. It is chosen such that the tree level part of the effective potential is simply the light-cone Lagrangian (Llc)pp rather than −(Llc)pp.AfterVeff is computed, the result could then be analytically continued back into Minkowski signature by replacing vE → ivM . In order to facilitate the comparison of the two sides, it is useful to first examine the form of the action. The supergravity effective action is given in Eq. (95). For the simple case when the two membranes have the same radius and in the limit where their separation in the X4 to X9 direction, z, is small, i.e., z r0, the supergravity effective action is given in Eq. (95). Putting w, the difference of the two radii to zero, and not keeping track of the exact coefficients, we could rewrite the supergravity result in term of matrix theory parameters:
v4 v2 µ2 V = α + + . (12) eff µ2z5 z3 z Hence we see that a comparison with supergravity means looking at order (α)1 on the matrix theory side.
2.3. The membrane limit and the graviton limit
We will first state the two limits we are interested in:
• membrane limit: z 1, (13) αµ z N, (14) αµ • graviton limit: z 1, (15) αµ z N, (16) αµ whereweusedz to denote the separation of the two spherical membranes in the X4 to X9 directions, and recall α = 1 . M3R The membrane limit is derived from the condition: 1 z 1. (17) N r0 The first inequality ensures that the effect of nonzero z is greater than any matrix theory correc- tions to supergravity, which we are not interested in. The second inequality ensures we are at the = αµN near membrane limit. Using r0 6 (see Eq. (5)), we arrive at the limit as stated. 6 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
The graviton limit is when z is much greater than r0, so that approximately the two spheres interact like two point like gravitons. We still enforce the condition 1 z for comparison with N r0 supergravity but reverse the second inequality in the membrane limit to z 1 to arrive at the r0 graviton limit stated above. In this paper, on the supergravity side we will calculate the light-cone Lagrangian only in the membrane limit. The light-cone Lagrangian in the graviton limit was already computed in an earlier work [2]. On the matrix theory side, the one-loop effective potential can be computed for general z. The two limits of this potential are then compared with the supergravity side and we will find perfect agreement. Later in the paper we will use matrix theory to find explicitly the potential that interpolates between the two limits.
3. The supergravity light-cone Lagrangian
The light-cone Lagrangian Llc of the probe is basically its Lagrangian Legendre transformed − in the x degree of freedom. Here we briefly derive the Llc for a probe membrane in a pp-wave background perturbed by some source. (The reader is referred to Section 4.1 of [14] for the detailed discussion of the light-cone Lagrangian of point particle probes and membrane probes in terms of Hamiltonian systems with constraints.) Denote the background metric and three-form as Gµν(x), Aµνρ(x), respectively, and the membrane embedding coordinates as Xµ(σ i), with σ i,i= 0, 1, 2 being the world-volume coor- dinates. The bosonic part of the membrane Lagrangian density is given by 1 L Xµ,∂ Xµ =−T − det(g ) − ij k A ∂ Xµ∂ Xν∂ Xρ , (18) i ij 6 µνρ i j k µ ν with T being the membrane tension, and gij ≡ Gµν∂iX ∂j X being the pullback metric. The momentum density is L ≡ ∂ =− − 0k µ − ν ρ Πλ λ T det(gij )g ∂kX Gλµ Aλνρ∂1X ∂2X . (19) ∂(∂0X ) Define ˜ ν ρ 0k µ Πλ ≡ Πλ − TAλνρ∂1X ∂2X =−T − det(gij )g ∂kX Gλµ. (20)
There are the following three primary constraints
λξ ˜ ˜ 2 λ φ0 ≡ G ΠλΠξ + T det(grs) = 0,φr ≡ Πλ∂r X , (21) where r, s = 1, 2 label the spatial world-volume coordinates. Now use the light-cone coordinates {x+,x−,xA}, and assume static, i.e., x+-independent, − A − A + 0 background Gµν(x ,x ) and Aµνρ(x ,x ). Use the light-cone gauge X = σ . Then the equa- + 0 tion of motion for X fixes the Lagrange multiplier c for the constraint φ0 1 c0 = . (22) +ξ ˜ 2G Πξ
+ 0 The constraint φ0 = 0 (in which X is set to σ ) can be used to solve for Π+ as a function of − A − A (X ,X ,∂r X ,∂r X ,Π−,ΠA). H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 7
− A − A A We shall also express ΠA as a function of (X ,X ,∂r X ,∂r X ,Π−,∂0X ) by inverting the equations of motion for XA: ∂XA GAξ Π˜ ∂XA = ξ + cr , (23) 0 +ν ˜ r ∂σ G Πν ∂σ where we have used the expression for c0 given above and the cr ’s are the Lagrange multipliers for the φr ’s. − A − A We can first define the light-cone Hamiltonian Hlc(X ,X ,∂r X ,∂r X ,Π−,ΠA) ≡−Π+, and then define the light-cone Lagrangian − A − A A A A Llc X ,X ,∂r X ,∂r X ,Π−,∂0X ≡ ΠA∂0X − Hlc = ΠA∂0X + Π+. (24) Now the background in which the probe membrane moves is the pp-wave perturbed by some source: Gµν = (Gµν)pp + hµν , and Aµνρ = (Aµνρ)pp + aµνρ , where the quantities with subscript pp are those of the unperturbed pp-wave background, and hµν,aµνρ are the perturbations caused by the source. We only need the light-cone Lagrangian to linear order in the perturbation. Solving φ0 = 0gives − A − A Π+ X ,X ; ∂r X ,∂r X ; Π−,ΠA = (Π+)pp + δΠ+, (25) where − 2 1 2 T A B B A 2 (Π+)pp = −g++Π− + ΠAΠA + ∂1X ∂2X − ∂1X ∂2X 2Π− 2 µ + T ∂ Xi ∂ Xj Xk (26) 3 ij k 1 2 and − = 1 − ˜ 2 − − ˜ δΠ+ h−−(Π+)pp 2Π−(h+− g++h−−)(Π+)pp 2Π− ˜ ν ρ ν ρ ˜ − 2T(Π+) a− ∂ X ∂ X − 2TΠ−a+ ∂ X ∂ X − 2h− Π (Π+) pp νρ 1 2 νρ 1 2 A A pp 2 2 ν ρ − Π− h++ + (g++) h−− − 2g++h+− + 2TΠ−g++a−νρ ∂1X ∂2X ν ρ − h Π Π − 2Π−(h+ − g++h− )Π − 2TΠ a ∂ X ∂ X AB A B A A A A Aνρ 1 2 2 + T δ det(grs) , (27) where − 2 ˜ 1 2 T A B B A 2 (Π+)pp = −g++Π− + ΠAΠA + ∂1X ∂2X − ∂1X ∂2X . (28) 2Π− 2 The Lagrange multipliers cr ’s are r = r + r = r c c pp δc δc , (29) r where we have taken (c )pp to be zero. Inverting Eq. (23) gives − A − A A ΠA X ,X ; ∂r X ,∂r X ; Π−,∂0X = (ΠA)pp + δΠA, (30) A r where (ΠA)pp = Π−∂0X , and δΠA depends on δc . Now A A Llc = (ΠA)pp∂0X + δΠA∂0X + (Π+)pp + δΠ+. (31) 8 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
However, as can be seen easily, the second term in the above expression is canceled by the part containing δΠA in the third term. So Llc does not contain δΠA, thus being independent of the δcr ’s. Hence one finds
Llc = (Llc)pp + δLlc, (32) with 1 A A Π− 2 1 i 2 1 a 2 (L ) = Π−∂ X ∂ X − µ X + X lc pp 2 0 0 2 9 36 2 T A B B A 2 µ i j k − ∂1X ∂2X − ∂1X ∂2X + T ij k ∂1X ∂2X X (33) 4Π− 3 and δLlc = δΠ+, where δΠ+ is given by Eq. (27) with all the ΠA’s in it replaced by (ΠA)pp = A Π−∂0X (the reason being that δΠ+ contains hµν and aµνρ and is thus already first-order). We r would like to emphasize again that Llc is independent of the choice of the δc ’s. To use the general expression for Llc given in the previous paragraph in a specific problem, one just has to substitute in the hµν ’s and aµνρ ’s found by solving the supergravity field equations for the specific source, as well as the specific probe membrane embedding XA(σ i).
4. Supergravity computation
4.1. Diagonalizing the field equations for arbitrary static source
In this subsection we present the diagonalization of the linearized supergravity equations of motion for arbitrary static sources (this subsection is basically Section 4.2 of [14]). There is, of course, no highbrow knowledge involved here: we are just solving the linearized Einstein equations and Maxwell equations, which are coupled; and by “diagonalization” we basically just mean the prescription by which we get a decoupled Laplace equation for each component of the metric and three-form perturbations. The unperturbed background is the 11D pp-wave, and we only consider static, i.e., x+-independent, field configurations, thanks to the fact that the sources considered are taken to be static, i.e., with x+-independent stress tensor and three-form current. Since we leave the source arbitrary, what we will present here are the left-hand side of the lin- earized equations. These are tensors whose computation is straightforward though a bit tedious: the reason we present them here is because they are necessary when solving the field equations, and to the best of our knowledge have not been explicitly given elsewhere. A somewhat related problem is the diagonalization of the equations of motion when the source is absent. This requires field configurations with x+-dependence. One work along this line is [15]. Roughly speaking, borrowing the language of electromagnetism, what is considered in [15] are electromagnetic waves in vacuum, while what we are considering here are electrostatics and magnetostatics for arbitrary static sources. The nonzero components up to (anti)symmetry of the Christoffel symbol, Riemann tensor, and Ricci tensor of the 11D pp-wave are
A 1 − 1 Γ++ =− ∂ g++,Γ= ∂ g++, 2 A +A 2 A 1 1 R+ + =− ∂ ∂ g++,R++ =− ∂ ∂ g++. (34) A B 2 A B 2 C C H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 9
(We usually do not substitute the explicit expression of g++, unless that brings significant sim- plification to the resulting formula.) Now let us perturb the pp-wave background by adding a source. Denote the metric perturba- tion δgµν by hµν , and the gauge potential perturbation by δAµνρ = aµνρ . hµν,aµνρ are treated as rank-two and rank-three tensors, respectively, the covariant derivative ∇ acting on them is defined using the connection coefficient of the unperturbed pp-wave background, and indices are raised/lowered, traces taken using the background metric gµν . ¯ ≡ − 1 ≡ µν Let us deal with the Einstein equations first. Define hµν hµν 2 gµνh, where h g hµν . Without the source, the Einstein equation is 1 R − Rg − κ2 [T ] = 0 (35) µν 2 µν 11 µν A and we recall that the stress tensor of the gauge field is
1 1 [T ] = F F λξρ − g F ρσλξF . µν A 2 µλξρ ν µν ρσλξ (36) 12κ11 8 The source perturbs the Einstein equation to
1 δ R − Rg − κ2 δ[T ] = κ2 [T ] (37) µν 2 µν 11 µν A 11 µν S with [Tµν]S standing for the stress tensor of the source. As usual, it helps to proceed in an organized manner, grouping different terms in the above − 1 =−1 ∇σ ∇ ¯ + + perturbed Einstein equations. One finds, δ(Rµρ 2 Rgµρ ) 2 σ hµρ Kµρ Qµρ , and 2 [ ] = + ∇σ ∇ ¯ κ11δ Tµν A Nµν Lµν , where the explicit expressions of the symmetric tensors σ hµν , Kµν , Qµν , Nµν , and Lµν can be obtained after some work. Their definitions and components are given below1:
σ ¯ •∇∇σ hµν : σ ¯ µν ¯ ¯ 1 ¯ ∇ ∇ h++ = g ∂ ∂ h++ + −(∂ ∂ g++)h+− + (∂ g++∂ g++)h−− σ µ ν A A 2 A A ¯ ¯ + 2[∂Ag++∂−h+A − ∂Ag++∂Ah+−], (38)
σ ¯ µν ¯ 1 ¯ ¯ ∇ ∇ h+− = g ∂ ∂ h+− − (∂ ∂ g++)h−− + ∂ g++∂−h− σ µ ν 2 A A A A ¯ − ∂Ag++∂Ah−−, (39)
σ ¯ µν ¯ 1 ¯ ¯ ¯ ∇ ∇ h+ = g ∂ ∂ h+ − (∂ ∂ g++)h− + ∂ g++∂−h − ∂ g++∂−h+− σ C µ ν C 2 A A C A AC C ¯ − ∂Ag++∂Ah−C, (40) σ ¯ µν ¯ ∇ ∇σ h−− = g ∂µ∂νh−−, (41) σ ¯ µν ¯ ¯ ∇ ∇σ h−C = g ∂µ∂νh−C − ∂Cg++∂−h−−, (42) σ ¯ µν ¯ ¯ ¯ ∇ ∇σ hCD = g ∂µ∂νhCD − ∂Cg++∂−h−D − ∂Dg++∂−h−C. (43)
1 µν 2 Notice that ∂+ will never appear because we only consider the static case; also note g ∂µ∂ν =−g++∂− + ∂A∂A for static configurations. 10 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
• Kµν : Its definition is 1 1 1 K ≡ R ξ h¯ + R ξ h¯ + Rσ ξ h¯ + g Rξσh¯ − Rh¯ . (44) µρ 2 µ ξρ ρ ξµ µρ σξ 2 µρ ξσ 2 µρ Its components are given by
1 ¯ 1 ¯ 1 ¯ K++ = − ∂A∂Ag++ h+− + g++h−− + (∂A∂B g++)hAB , (45) 2 2 2 1 ¯ K+− = − ∂A∂Ag++ h−−, (46) 2 1 ¯ 1 ¯ K+ = − ∂ ∂ g++ h− + − ∂ ∂ g++ h− , (47) A 4 C C A 2 A B B K−− = 0, (48) = K−A 0, (49) 1 1 ¯ K = ∂ ∂ g++ − δ ∂ ∂ g++ h−−. (50) AB 2 A B 2 AB C C • ≡ 1 ∇ +∇ − 1 ∇α ≡∇β ¯ Qµν : Its definition is Qµρ 2 ( µqρ ρqµ) 2 gµρ qα, where qα hβα. As one can recognize, Qµρ contains the arbitrariness of making different gauge choices when solv- ing the Einstein equation, where one makes a gauge choice by specifying the qµ’s. The components of Qµρ are 1 1 Q−− = ∂−q−,Q− = (∂−q + ∂ q−), Q−+ = (g++∂−q− − ∂ q ), A 2 A A 2 A A 1 1 Q = (∂ q + ∂ q ) − δ (∂−q+ − g++∂−q− + ∂ q ), AB 2 A B B A 2 AB A A 1 Q+ = ∂ q+ − (∂ g++)q− , A 2 A A 1 1 Q++ = (∂ g++)q − g++(∂−q+ − g++∂−q− + ∂ q ). (51) 2 A A 2 A A • 2 [ ] Nµν : It is defined to be the part of κ11δ Tµν A that contains only the metric perturbation, but not the three-form gauge potential perturbation. Its components are given by 3 9 2 1 ¯ 1 ¯ 1 ¯ 1 ¯ N++ = µ h+− + g++h−− − h + h , 3 12 3 ii 6 aa i=1 a=4 2 2 µ ¯ µ ¯ N+− = h−−,N+ = h− ,N+ = 0, 4 i 2 i b N−− = 0,N−i = 0,N−b = 0, 2 2 µ ¯ µ ¯ N =− δ h−−,N= 0,N= δ h−−. (52) ij 4 ij ib ab 4 ab • 2 [ ] Lµν : This is defined to be the part of κ11δ Tµν A that contains only the three-form perturba- tion, but not the metric perturbation. Its components are given by
1 L++ = µ δF + − g++δF − ,L+− = 0, 123 2 123 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 11
µ µ L+ = δF+ −,L+ = δF , i 4 ij k jk b 2 123b L−− = 0,L−i = 0,L−b = 0, µ µ µ L = δ δF −,L= δF −,L=− δ δF −. (53) ij 2 ij 123 ib 4 ij k bj k bd 2 bd 123 Next let us deal with the Maxwell equation. In the absence of the source, it is 1 √ η˜ µ1...µ11 √ ∂ −gF λµ1µ2µ3 − √ F F = 0, (54) −g λ 1152 −g µ4...µ7 µ8...µ11 where η˜ is either +1or−1 depending on the choice of convention, which one can fix later by requiring the consistency of the conventions for the equations and the solutions under consider- ation. When the source is present, we add its current J µ1µ2µ3 to the left-hand side of the above equation, and get 1 √ η˜ µ1...µ11 δ √ ∂ −gF λµ1µ2µ3 − √ F F = J µ1µ2µ3 . (55) −g λ 1152 −g µ4...µ7 µ8...µ11 We can write the left-hand side of the above equation as the sum of two totally antisymmetric tensors Zµ1µ2µ3 + Sµ1µ2µ3 , where Zµ1µ2µ3 is defined to be the part that contains the metric perturbation only, and Sµ1µ2µ3 is defined to be the part that contains the three-form perturbation only. One finds +−i ¯ +−b +ij ¯ ¯ Z = µij k ∂j h−k,Z= 0,Z= µij k (∂−h−k − ∂kh−−), + + Z ib = 0,Zbc = 0, 3 −ij 1 ¯ ¯ − ¯ ¯ Z = µ ∂ h − h− − h − ∂ h , ij k k 3 ii b kb i=1 − − Z ib = µ ∂ h¯ ,Zbc = 0, ij k j kb 3 ij k 1 ¯ ¯ − ¯ ¯ Z =−µ ∂− h − h− − h − ∂ h− , ij k 3 ii b b i=1 ij b ¯ ¯ Z = µij k (∂−hkb − ∂kh−b), Zibc = 0,Zbce = 0, (56) and +−A µν µ µ S = g ∂ ∂ a−+ + ∂ g++∂−a − − ∂− ∇ a + + ∂ ∇ a +− , (57) µ ν A B BA µ A A µ +AB µν µ µ µ S = g ∂µ∂νa−AB − ∂− ∇ aµAB + ∂A ∇ aµ−B − ∂B ∇ aµ−A , (58) −AB µν +AB S = g ∂ ∂ a+ − g++S µ ν AB µ + (∂Ag++)(∂−a−+B ) + ∂A ∇ aµ+B − ∂A(aEB−∂Eg++) −[A ↔ B] ˜ η −ABµ ...µ 123+ − (∂ g++)δF − − µ 4 7 δF , (59) D D AB 24 µ4...µ7 ABE µν S = g ∂ ∂ a − (∂ g++)(∂−a− ) − (∂ g++)(∂−a− ) − (∂ g++)(∂−a− ) µ ν ABE A BE B EA E AB µ µ µ − ∂A ∇ aµBE − ∂B ∇ aµEA − ∂E ∇ aµAB η˜ + − µ ABEµ4...µ7123 δF . (60) 24 µ4...µ7 12 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
µ µ µ µ Notice that S 1 2 3 contains ∇ aµρλ and its derivatives. Those terms correspond to the gauge freedom for the three-form gauge potential. Now that we have collected the expressions for the various tensors, we are ready to diagonalize the field equations. Recall that the Einstein equation is 1 − ∇σ ∇ h¯ + K + Q − N − L = κ2 [T ] (61) 2 σ µν µν µν µν µν 11 µν S and the Maxwell equation is
Zµ1µ2µ3 + Sµ1µ2µ3 = J µ1µ2µ3 . (62) The right-hand side of these equations is given by specifying the source that we consider (recall 2 that the three-form current J is of order κ11), hence we only need to concentrate on diagonalizing the left-hand sides. As will be seen shortly, it is useful to define “level” for tensors: lower +/upper − indices con- tribute +1 to the level; lower −/upper + indices contribute −1 to level; and the upper A/lower A indices contribute zero to the level. We shall see that the field equations should be solved in ascending order of their levels. The following is the detailed prescription of the diagonalization procedure. Let us use the shorthand notation (E.E.)µν for the lower (µν) component of the Ein- µ µ µ stein equation, and (M.E.) 1 2 3 for the upper (µ1µ2µ3) component of the Maxwell equation.
• At level −2.
The only field equation at this level is (E.E.)−−, which reads, upon using the expressions of σ ¯ the various tensors ∇ ∇σ hµν,Kµν,Qµν,..., etc. that we have given above
1 µν ¯ 2 − g ∂ ∂ h−− + Q−− = κ [T−−] . (63) 2 µ ν 11 S ¯ This equation can be immediately solved for h−− after specifying the source term and the gauge choice term Q−−.
• At level −1.
We have (E.E.)−A, which reads 1 µν ¯ ¯ 2 − g ∂ ∂ h− − (∂ g++)(∂−h−−) + Q− = κ [T− ] , (64) 2 µ ν A A A 11 A S ¯ ¯ which can now be solved for h−A,usingtheh−− found previously. Also at this level is (M.E.)+AB , which reads, µν µ µ µ ¯ ¯ g ∂µ∂νa−ij − ∂− ∇ aµij + ∂i ∇ aµ−j − ∂j ∇ aµ−i + µij k (∂−h−k − ∂kh−−) + = J ij , µν µ µ µ +ib g ∂ ∂ a− − ∂− ∇ a + ∂ ∇ a − − ∂ ∇ a − = J , µ ν ib µib i µ b b µ i µν µ µ µ +bc g ∂µ∂νa−bc − ∂− ∇ aµbc + ∂b ∇ aµ−c − ∂c ∇ aµ−b = J , (65) µ from which we can find a−AB , upon specifying the gauge choice ∇ aµρλ for the three-form and ¯ ¯ using the h−A and h−− found previously. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 13
• At level 0.
+−A ABE At this level we have (E.E.)+−, (M.E.) , (E.E.)AB , and (M.E.) . (E.E.)+− is of the form
1 µν ¯ − g ∂ ∂ h+− = known terms. (66) 2 µ ν (From now on, we will not bother writing down the detailed equations; “known terms” refers µ to the gauge choice terms Qµν , ∇ aµρλ , source terms, and terms containing previously found ¯ hµν ’s and aµνρ ’s, one can write those down by looking up the expressions given earlier for the ¯ +−A various tensors.) Hence solving it we get h+− and solving (M.E.) gives a−+A. ABE (E.E.)AB and (M.E.) are coupled, so a little more work is needed. The following are the bce details. First notice that the only unknown in (M.E.) is abce, hence solving this equation we bce µν find abce.((M.E.) contains the usual term g ∂µ∂νabce and also a term of the form ∂−adf g which comes from the F ∧ F in the Maxwell equation, hence it is not quite a Laplace equation. But, that being said, one should not have any difficulty solving it.) ibc µν ij b (M.E.) is of the form g ∂µ∂νaibc = known terms, solving which gives aibc. (M.E.) and (E.E.)kb are coupled in the following manner µν ¯ g ∂µ∂νaij b + µij k ∂−hkb = known terms, 1 µν ¯ 1 − g ∂ ∂ h + µ ∂−a = known terms. (67) 2 µ ν kb 4 klm lmb ¯ Decoupling these two equations is quite easy. Let us take a12b and h3b as the representative case. One sees that these two equations can be recombined to give µν ¯ g ∂ ∂ + iµ∂− (h + ia ) = known terms, µ ν 3b 12b µν ¯ g ∂µ∂ν − iµ∂− (h3b − ia12b) = known terms. (68) ¯ + ¯ − ¯ Solving these equations gives (h3b ia12b) and (h3b ia12b), and in turn h3b and a12b. (M.E.)ij k is coupled to (E.E.) and (E.E.) through the quantity H ≡ 2 3 h¯ − ij bd 3 i=1 ii 1 9 ¯ 3 a=4 haa in the following manner µν g ∂µ∂νa123 + µ∂−H = known terms, 1 µν ¯ 1 − g ∂ ∂ h + µδ ∂−a = known terms, 2 µ ν ij 2 ij 123 1 µν ¯ 1 − g ∂ ∂ h − µδ ∂−a = known terms. (69) 2 µ ν bd 2 bd 123 Combining the last two equations gives
µν −g ∂µ∂νH + 4µ∂−a123 = known terms. (70) Recombining this with first equation, we get µν g ∂ ∂ + 2iµ∂− (H + 2ia ) = known terms, µ ν 123 µν g ∂µ∂ν − 2iµ∂− (H − 2ia123) = known terms, (71) ¯ solving which gives H and a123. Using the resulting expression for a123 one can then find hij ¯ ABE and hbd . Thus we are done with (E.E.)AB and (M.E.) . 14 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
• At level 1.
−AB µν (M.E.) is of the form g ∂µ∂νa+AB = known terms, solving which gives a+AB ,
1 µν ¯ ¯ (E.E.)+ is of the form − g ∂ ∂ h+ = known terms, solving which gives h+ . A 2 µ ν A A • At level 2.
1 µν ¯ ¯ (E.E.)++ is of the form − g ∂ ∂ h++ = known terms, solving which gives h++. 2 µ ν Thus we have diagonalized the whole set of Einstein equations and Maxwell equations.
4.2. Application to a spherical membrane source using the near-membrane expansion
Now let us apply the general formalism of the previous subsection to the case of interest, with the source being a spherical membrane sitting at the origin of the transverse directions, 1 2 + 2 2 + 3 2 = 2 4 = 9 = + = − = i.e., having (X ) (X ) (X ) r0 , X 0,...,X 0, and X t,X 0. The gauge µ choice we shall take is: qα = 0 (hence all the Qαβ ’s vanish); and ∇ aµρλ = 0. The nonzero components of the stress tensor and three-form current for this source are given by
−1 2 − 4 9 µr0 µr0 [T−−] = Tδ(x )δ(r − r )δ x ···δ x , [T+−] =− [T−−] , s 0 3 s 3 s
2 i j 4 µr0 x x ij µr0 [T ] = − δ [T−−] , [T++] = [T−−] ij s 2 s s s (72) 3 r0 3 and
k +ij = 2 − µr0 x [ ] J κ11( 2) ij k T−− s, (73) 3 r0 √ where r ≡ xixi . √ Now let us explain√ what we mean by “near-membrane expansion”. Define w ≡ r − r0, z ≡ xaxa, and ξ ≡ w2 + z2, which parameterize the distance away from the source membrane. We shall assume that w,z,ξ are of the same order of magnitude. The near-membrane expansion is an expansion in ξ/r0. When one sits very close to membrane, one just sees a flat membrane, which is the zeroth order of the expansion. As one moves away from the membrane, one begins to feel the curvature of membrane, which gives the higher order corrections in the expansion. One should also note that the zeroth order of this expansion is just the flat space limit: µ → 0, r0 →∞, with µr0 kept finite. It is instructive to see how the zeroth order works. At this order, in the Einstein equations and ¯ Maxwell equations, the effective source terms (which are of the forms (∂g++)hµν , etc.) arising µνρ from the various tensors Kµν , Nµν , Lµν , Z , etc. are less singular then the delta-functional µνρ sources [Tµν]s and J , and can thus be thrown away. Then the resulting equations are trivially i µνρ decoupled. Also, at this order, we can treat the x ’s in [Tµν]s and J as constant vectors. One finds (using the subscript 0 to denote “zeroth order”)
k ¯ µr0 x ¯ [h−A]0 = 0, [a−ij ]0 = ij k [h−−]0, [a−ib]0 = 0, [a−bc]0 = 0, 3 r0 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 15
2 ¯ µr0 ¯ [h+−] =− [h−−] , [a−+ ] = 0, 0 3 0 A 0
2 i j ¯ µr0 x x ij ¯ [h ] = − δ [h−−] , ij 0 2 0 3 r0 [h¯ ] = 0, [h¯ ] = 0, [a ] = 0, ib 0 bc 0 ABE 0 3 k µr0 x ¯ [a+ij ]0 =− ij k [h−−]0, [a+ib]0 = 0, [a+bc]0 = 0, 3 r0
4 ¯ ¯ µr0 ¯ [h+ ] = 0, [h++] = [h−−] , (74) A 0 0 3 0 ¯ where [h−−]0 satisfies 2 2 2 2 µr0 2 ∂ ∂ ∂ ¯ − k− + + +···+ [h−−] 3 (∂w)2 (∂x4)2 (∂x9)2 0 − 1 µr 1 =− 0 κ2 T δ(w)δ x4 ···δ x9 (75) πR 3 11 1 (where we have multiplied the right-hand side of the equation by 2πR due to the Fourier trans- form along the x− direction), and is given by − µr0 2 ¯ exp( 3 k−ξ) µr0 µr0 [h−−] = ∆ 3 + 3 k−ξ + k−ξ , (76) 0 ξ 5 3 3
κ2 T ≡ 11 with ∆ 4 µr0 . 16π R( 3 ) Plugging the above zeroth order solution of the field equations into the light-cone Lagrangian L δ lc given in Eq. (27), for a spherical probe membrane with radius r0, sitting at rest in the 1 2 + 2 2 + 3 2 = 2 1, 2, 3 directions, and moving about in the 4 through 9 directions: (X ) (X ) (X ) r0 , 4 9 + = − = = µΠ− X (t), . . . , X (t), and X t,X 0, one finds (using the facts that T 3r sin θ , and we can = 0 set r0 r0 because we are looking at the zeroth order) 1 ¯ ˙ a ˙ a 2 δL = Π−[h−−] X X . (77) lc 8 0 ¯ It is worth noting that, keeping only the leading order term in k− in [h−−]0,Eq.(77) becomes the v4 Lagrangian for the case of longitudinal momentum transfer between two membranes in the flat space, given in [16]. Now let us go on to consider higher orders in the near-membrane expansion. Since in this paper we do not consider longitudinal momentum transfer, we shall set k− = 0 (which makes many fields equations decouple). Denote
2 ≡ ∂A∂A (when acting on functions of (w, z))
= 20 + δ2, (78) 2 2 with 2 ≡ ∂ + ∂ + 5 ∂ = ∂2 + ∂ ∂ being the zeroth order Laplace operator, and δ2 ≡ 0 ∂w2 ∂z2 z ∂z w a a 2 ∂ being curvature correction to it. r0+w ∂w 16 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
• At level −2.
(E.E.)−−
−1 ¯ 1 2 4 9 µr0 2h−− =− κ T δ(w)δ x ···δ x . (79) πR 11 3 Let ¯ ¯ ¯ h−− =[h−−]0 + δh−−, (80) ¯ 3 with [h−−] ≡ ∆ , which satisfies 0 ξ 5
−1 ¯ 1 2 4 9 µr0 2 [h−−] =− κ T δ(w)δ x ···δ x . (81) 0 0 πR 11 3 Then (E.E.)−− becomes ¯ ¯ ¯ 20δh−− + δ2[h−−]0 + δ2δh−− = 0. (82)
Now let us look at the order of magnitude of each term in the above equation. Notice that 20 ∼ 1 2 ∼ 1 ∼ 1 ¯ ∼ 1 2 , δ . The second term is thus ∆ 6 , which tells us δh−− ∆ 4 . Solving the ξ r0ξ r0ξ r0ξ equation iteratively we find
2 3 4 ¯ ¯ w w w w δh−− =[h−−] − + − + (83) 0 r 2 3 4 0 r0 r0 r0 and thus 5 ¯ 1 r0 ξ h−− = 3∆ + O . (84) 5 r + w 5 ξ 0 r0 ξ 5 We did not compute the O( 5 ) terms, because we are only interested in the part of the solution r0 that is singular as ξ → 0. Solving the other field equations is similar, so we just present the results ξ 5 below, omitting the O( 5 ) symbol. r0
• At level −1. ¯ h−A = 0,
a−bc = 0, a− = 0, ib 2 + 2 2 2 2 k 1 1 w 1 w z 1 wz w z a− = x µ∆ 1 − + − + . (85) ij ij k 5 r 2 3 4 ξ 2 0 6 r0 2 r0 r0
• At level 0. 2 2 2 ¯ µr0 3∆ 5 w 7 z r0 h+− =− 1 + + , 5 2 2 + 3 ξ 4 r0 8 r0 r0 w a−+A = 0, H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 17
aABD = 0, (86) 1 1 5 w 17 w2 1 z2 3 w3 1 wz2 h¯ = xkxbµ2∆ − − + − − + kb 5 1 2 2 3 3 ξ 2 4 r0 12 r 12 r 2 r 4 r 0 0 0 0 3 w4 1 w2z2 + − , (87) 2 r4 2 r4 0 0 xixj µr 2 3 w z2 w3 wz2 w4 w2z2 z4 h¯ = 0 ∆ − − + + − − + ij 2 5 1 2 3 2 3 4 4 4 r 3 ξ r0 r r r r r r 0 0 0 0 0 0 µr 2 3 1 w 7 w2 1 z2 1 w3 3 wz2 1 w4 − δij 0 ∆ + + − − + + 5 1 2 2 3 3 4 3 ξ 2 r0 12 r 24 r 4 r 8 r 4 r 0 0 0 0 0 1 w2z2 1 z4 − + , (88) 24 r4 3 r4 0 0 1 1 w 7 w2 19 z2 7 w3 19 wz2 7 w4 19 w2z2 h¯ = δ (µr )2∆ − − + + − − . bd bd 0 5 r 2 2 3 3 4 4 ξ 2 0 18 r0 72 r0 18 r0 72 r0 18 r0 72 r0 (89) • At level +1.
a+bc = 0, a+ = 0, ib 2 2 2 3 2 4 k µr0 1 w w 5 z 3 w 7 wz 3 w a+ =− x µ∆ + − − + + − ij ij k 5 1 2 2 2 3 3 4 3 ξ r0 r 4 r 2 r 4 r 2 r 0 0 0 0 0 5 w2z2 1 z4 − + , 4 4 (90) 4 r0 2 r0 ¯ h+A = 0. (91)
• At level +2. 4 2 2 3 2 4 ¯ µr0 3 5 w 31 w 1 z 17 w 1 wz 1 w h++ = ∆ + + − + − + 5 1 2 2 3 3 4 3 ξ 2 r0 12 r 24 r 12 r 12 r 3 r 0 0 0 0 0 23 w2z2 17 z4 + + . 4 4 (92) 24 r0 32 r0 (A note in passing: as it turns out, the choice of convention for η˜ in the Maxwell equation (54) does not matter, since it drops out when solving the field equations for our specific source.) = + 1 2 + 2 2 + Again, let the probe membrane have a radius r0 r0 w, with the trajectory√ (X ) (X ) 3 2 = 2 4 9 + = − = ≡ ˙ a ˙ a (X ) r0 , X (t), . . . , X (t), and X t,X 0. We shall take v X X to be of order = 1 µz (recall that the supersymmetric circular orbit has v 6 µz; so here we are considering gener- ically nonsupersymmetric orbits that can be regarded as deformations of the supersymmetric µΠ− circular one). Plugging in the supergravity solution given above into Eq. (27),usingT = , 3r0 sin θ and in the end keeping only the part of δLlc that is singular as ξ → 0, we find the probe’s Llc to 18 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 be
Llc = (Llc)pp + δLlc (93) with (Llc)pp being the action in the unperturbed pp-wave background, and µ4(4w2z2 + 7z4) − 72µ2v2(2w2 + 5z2) + 3888v4 δL = Π−∆ . (94) lc 10368ξ 5 L → ∼ 1 Notice that the above δ lc’s singular behavior as ξ 0 is homogeneous: ξ (since v is of order µz). The expression in the numerator: µ4(4w2z2 + 7z4) − 72µ2v2(2w2 + 5z2) + 3888v4 nicely factorizes into (36v2 − z2µ2)[108v2 − (4w2 + 7z2)µ2], which shows that for the special case of = 1 the supersymmetric circular orbit v 6 µz considered in [10], δLlc vanishes as expected. To be more precise, so far we have been talking about Lagrangian density. Since the membrane worldvolume is taken to be a unit sphere, the Lagrangian is given by κ2 T 2 (36v2 − z2µ2)[108v2 − (4w2 + 7z2)µ2] δL = dθ dφδL = 11 lc lc 4π 3Rµ2 1152ξ 5 α (36v2 − z2µ2)[108v2 − (4w2 + 7z2)µ2] = , (95) µ2 1152ξ 5 = µΠ− wheretogetthefirstlineweusedT to eliminate Π− in terms of T and r0, and set 3r0 sin θ r ≈ r in the end to remove higher w order curvature correction. To get the second line we used 0 0 r0 2 the expressions for κ11, T , and α in terms of M and R given at the end of Section 2.1. Here we would like to make a brief comparison of the above membrane result with the gravi- ton result given in [2]. First of all, the membrane result contains the variable w (the difference in radius between the probe membrane and the source membrane), which has no counterpart in the graviton case. Secondly, in terms of the x1 through x3 directions, the two membranes are sitting at rest at the i origin; this corresponds to setting x = 0 and vi = 0 in the graviton case. If we set w = 0inEq.(95), i.e., consider two membranes of the same size, then ξ = z, and α (δL ) = 36v2 − z2µ2 108v2 − 7z2µ2 , (96) lc membrane 1152µ2z5 i while for the gravitons, upon setting Np = Ns = N, x = 0, and vi = 0 in Eq. (19) of [2],we have α3N 2 (δL ) = 36v2 − z2µ2 140v2 − 7z2µ2 . (97) lc graviton 5736z7 When comparing the above two expressions, note the difference between the numerators: (108v2 − 7z2µ2) for membrane, and (140v2 − 7z2µ2) for graviton. Also notice their differ- ent power law dependence on z: 1 for membrane and 1 for graviton, which cannot be undone z5 z7 by integrating over the membrane (r0, the radius of the membrane, has nothing to do with z,the separation of the membranes in the X4 to X9 directions).
5. Matrix theory computation—the membrane limit
Shin and Yoshida have previously calculated the one-loop effective action for membrane fuzzy spheres extended in the first three directions and having periodic motion in a sub-plane of the H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 19 remaining six transverse directions. In Ref. [10] they considered the case of supersymmetric µ circular motion for an arbitrary radius and angular frequency 6 (this orbit preserves eight super- symmetries and was first found by [1]). Here we generalize that analysis to orbits which are not required to satisfy the classical equations of motion and are in general nonsupersymmetric. The procedures are: expanding the action to quadratic order in fluctuations, writing the fields in terms of the matrix spherical harmonics introduced in [17],2 diagonalizing the mass matrices of the bosons, fermions, and ghosts, and finally summing up the masses to get the one-loop effective action. In doing so, we shall adopt the notations of [10]. We shall consider the background BI 0 I = (1) B I , (98) 0 B(2) where
i µ i a B = J × ,B= 0 · 1N ×N , (1) 3 (1)N1 N1 (1) 1 1 i µ i i a a B = J × + x (t)1N ×N ,B= x (t)1N ×N (99) (2) 3 (2)N2 N2 2 2 (2) 2 2 with the J i ’s being su(2) generators. The above background has the interpretation of one spher- ical membrane (labeled by the subscript (1)) sitting at the origin, and the other spherical mem- brane (labeled by the subscript (2)) moving along the arbitrary orbit given by {xi(t), xa(t)}. The fluctuations of the bosonic fields are given by Z0 Φ0 ZI ΦI = (1) I = (1) A 0 † 0 ,Y I † I . (100) (Φ ) Z(2) (Φ ) Z(2) As it turns out, the part of the bosonic action containing the diagonal fluctuations Z0,ZI does not contain any new terms in addition to those given in [10], so we do not write it out here. (It will be the same for the fermionic and ghost parts of the action; it is the off-diagonal fluctuations that give the one-loop interaction potential between the two membranes.) The action for the off-diagonal fluctuations is µ 2 µ S = dt Tr −Φ˙ 02 + x2Φ02 + J i ◦ Φ02 − 2 xi Re J i ◦ Φ0 Φ0 † OD 3 3
µ 2 µ 2 + Φ˙ i2 − x2Φi2 − Φi + iij k J j ◦ Φk2 + J i ◦ Φi2 3 3 µ + 2 xi Re J i ◦ Φ0 Φ0 † − 2ix˙i Φ0 †Φi − Φi †Φ0 3
µ2 µ 2 + Φ˙ a2 − x2 + Φa2 − J i ◦ Φa2 2 6 3 µ + 2 xi Re J i ◦ Φa Φa † − 2ix˙a Φ0 †Φa − Φa †Φ0 , (101) 3 where x2 ≡ xixi + xaxa, and dots mean time-derivatives.
2 For the computation below, one only needs the transformation of the matrix spherical harmonics under SU(2),how- ever, for detailed construction of the matrix spherical harmonics, see, e.g., Appendix A of [18]. 20 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
We specialize to the case where, xi = 0,x8 = b,x9 = vt, with (x8)2 + (x9)2 denoted by z2. The effective potential will be computed by summing over the masses of the fermionic and ghost fluctuations and then subtracting the masses of the bosonic fluctuations. This method, which we will refer to as the sum over mass method, is the same as the one used in [2]. Although the above trajectory has the form of a straight line with constant velocity in the (x8,x9) plane, the a 2 a 2 final expression of Veff in terms of z ≡ (x ) and v ≡ (x˙ ) should suffice for the purpose of comparing with supergravity for arbitrary orbits xa(t). One may ask whether the sum over mass formula is valid when the masses of the fluctuations are time-dependent (one origin for such a time-dependence is the acceleration of the trajectory). In Section 4.2 and Appendix A of [2], it was carefully shown that, in the case of two-graviton interaction in the pp-wave background, the sum over mass formula was sufficient in computing the terms that could occur on the super- gravity side. The time-dependence in the masses of the fluctuations will give terms of the form of matrix theory corrections to supergravity (i.e., terms that dominate at extreme short distances and cannot be observed in supergravity), which do not concern us since we are only interested in the comparison with supergravity. Here we expect a similar argument to hold in the case of two-membrane interaction. The rather nontrivial agreement with supergravity presented at the end of this section and also the agreement with the work of Shin and Yoshida [11] in Section 6.2 confirm the validity of the sum over mass method. We expand the fields in terms of matrix spherical harmonics 1 + − 2 (N1N2) 1 j × 0,I = 0,I N1 N2 Φ φjmYjm . (102) = 1 | − | m=−j j 2 N1 N2 For our choice of background the masses of modes in the i = 1, 2, 3 and a = 4, 5, 6, 7, 8 direc- tions are the same as those in [10]. For the gauge field and a = 9 direction we have 1 + − 2 (N1N2) 1 2 ˙ µ S = dt −φ0 2 + z2 + j(j + 1) φ0 2 jm 3 jm 1 j= |N1−N2| 2 µ 2 µ 2 + 9˙ 2 − 2 + + + 0 2 φjm z j(j 1) φjm 6 3 − 0 ∗ 9 − 9 ∗ 0 2iv φjm φjm φjm φjm , (103) where there is an implicit sum over −j m j. It is straightforward to diagonalize the mass matrix for these modes. Combining all contributions from bosonic fluctuations we get a bosonic effective potential3 given by, 1 (N +N )−2 2 12 µ 2 V B =− (2j + 1) z2 + (j + 1)2 eff 3 = 1 | − |− j 2 N1 N2 1 1 (N +N ) 2 1 2 µ 2 − (2j + 1) z2 + j 2 3 = 1 | − |+ j 2 N1 N2 1
3 Note that our convention differs from that of [10] by an overall minus sign. See Section 2. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 21 1 (N +N )−1 2 12 µ 2 − (2j + 1) z2 + j(j + 1) 3 = 1 | − | j 2 N1 N2 1 (N +N )−1 2 12 µ 2 1 2 − 5(2j + 1) z2 + j + 3 2 j= 1 |N −N | 2 1 2 1 (N +N )−1 2 12 µ 2 1 1 µ 4 − (2j + 1) z2 + j(j + 1) + + + 16v2 3 8 2 6 j= 1 |N −N | 2 1 2 µ 2 1 1 µ 4 + z2 + j(j + 1) + − + 16v2 . (104) 3 8 2 6 For the fermions we start with the action given in [10]
µ L = Tr iΨ †Ψ˙ − Ψ †γ I Ψ,BI − i Ψ †γ 123Ψ (105) F 4 with
Ψ χ = (1) Ψ † . (106) χ Ψ(2) Again the action for the diagonal fluctuations has no new terms, and for the off-diagonal part we decompose the SO(9) spinor χ using the subgroup SO(3) × SO(6) ∼ SU(2) × SU(4) preserved by pp-wave, 16 → (2, 4) + (2¯, 4¯)
1 χAα χ = √ . (107) 2 χˆ Aα
Substituting this into LF , 1 µ L = Tr i χ† Aαχ˙ − χ† Aαχ − χ† Aα σ i β J i ◦ χ F Aα 4 Aα 3 α Aβ + ˆ † Aα ˙ˆ + 1 ˆ † Aα ˆ + µ ˆ † Aα i β i ◦ˆ i χ χ Aα χ χAα χ σ αJ χAβ 4 3 + xi − χ† Aα σ i β χ + χˆ † Aα σ i β χˆ α Aβ α Aβ + a † Aα a ˆ B + ˆ † α a AB † x χ ρAB χα χ A ρ χBα , (108) with the σ i ’s and ρa’s being the gamma matrices for SO(3) and SO(6), respectively. We now consider our specific background and expand in modes,
1 (N +N )−3/2 2 1 2 1 3 L = iπ† π˙ + iπˆ † πˆ˙ − j + π † π −ˆπ † πˆ F jm jm jm jm 3 4 jm jm jm jm j= 1 |N −N |−1/2 2 1 2 + a † a ˆ +ˆ† a † x πjmρ πjm πjm ρ πjm 22 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
1 (N +N )−1/2 2 1 2 1 1 + iη† η˙ + iηˆ† η˙ˆ − j + η† η −ˆη† ηˆ jm jm jm jm 3 4 jm jm jm jm j= 1 |N −N |+1/2 2 1 2 + a † a ˆ +ˆ† a † x ηjmρ ηjm ηjm ρ ηjm , (109) where again there is an implicit sum over m. This can be diagonalized and contributes to the effective action 1 (N +N )−3/2 2 1 2 µ 2 3 2 V F = 2 (2j + 1) z2 + j + + v eff 3 4 j= 1 |N −N |−1/2 2 1 2
µ 2 3 2 + z2 + j + − v 3 4 1 (N +N )−1/2 2 1 2 µ 2 1 2 + 2 (2j + 1) z2 + j + + v 3 4 j= 1 |N −N |+1/2 2 1 2
µ 2 1 2 + z2 + j + − v . (110) 3 4 There is also the contribution from the ghosts, however, this has no new terms for our choice of background, 1 (N +N )−1 2 12 µ 2 V G = 2 (2j + 1) z2 + j(j + 1) (111) eff 3 = 1 | − | j 2 N1 N2 and the total effective action is the sum of the three pieces = B + F + G Veff Veff Veff Veff. (112) We introduce the variables N and u,
N1 = N + 2u, N2 = N, (113) where u will be related to the difference in radii of the two spheres and from now on we will restore α using dimensional analysis. Define
4 2 2 1 αµ 2 2 2 2 η± = z ± + 16α v ,ν± = z ± αv, (114) 2 6 and assuming that u 1 (rather than assuming u 0 because the lower limit of the first summa- tioninEq.(104) has to be nonnegative) we can write the effective action so that j always starts from 0 and finishes at N − 1. N−1 2 1 1 αµ 2 V =− (2j + 2u − 1) z2 + (j + u)2 eff α 3 j=0 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 23 2 1 αµ 2 + (2j + 2u + 3) z2 + (j + u + 1)2 3 2 1 αµ 2 + (2j + 2u + 1) z2 + (j + u)(j + u + 1) 3 2 1 αµ 2 + 5(2j + 2u + 1) z2 + (j + u + 1/2)2 3 1 2 2 2 αµ 1 + (2j + 2u + 1) η+ + (j + u)(j + u + 1) + 3 8 1 2 2 2 αµ 1 + η− + (j + u)(j + u + 1) + 3 8 1 2 2 2 αµ 2 − 2(2j + 2u) ν+ + (j + u + 1/4) 3 1 2 2 2 αµ 2 + ν− + (j + u + 1/4) 3 1 2 2 2 αµ 2 − 2(2j + 2u + 2) ν+ + (j + u + 3/4) 3 1 2 2 2 αµ 2 + ν− + (j + u + 3/4) 3 2 1 αµ 2 − 2(2j + 2u + 1) z2 + (j + u)(j + u + 1) . (115) 3 Let us write the above summation as
N−1 Veff = V(j) (116) j=0 and use Euler–Maclaurin sum formula to convert it into integrals
N 1 1 V = V(j) dj − V(N) + V(0) + V (N) − V (0) eff 2 12 0 1 − V (N) − V (0) +···. (117) 720 It is useful to first see what we are expecting from such an integral. From Eq. (115),wesee that a typical term of V(j) is roughly of the form: 1 j j 2 1 z 2 v V(j) = j z2 + αv + α2µ2j 2 = µN 2 + + . (118) α N N N 2 αµ αµ2 24 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
Plugging the above form of V into Eq. (117) and defining new variables γ = j/N and ζ = (( z )2 + v ), not keeping track of the exact coefficients, we have: αµ αµ2 1 ζ 2 ζ 2 V = µN 3 dγ γ γ 2 + + µN 2 1 + eff N N 0 2 1 2 1 2 ζ µ 3 2 ζ + µN ∂γ γ γ + + ∂ γ γ + +··· N N γ N 0 0 ζ 1 ζ 1 ζ = µN 3 F + F + F +··· , (119) 0 N N 1 N N 2 2 N where F are functions of ζ = 1 (( z )2 + v ) originating from the nth derivative of γ (we n N N αµ αµ2 ζ shall see in the next paragraph why we use N as the argument for Fn). Note that each Fn is 1 weighted by a factor Nn . First we look at the term F0. Assuming that in a power expansion of F0(x) there exists a term 4 x3, after expanding in α it would contribute to V a v4 term of the form µN 3( ζ )3 ∼ αv , which eff N µ2z5 has the correct form of membrane interactions and has the correct order of N (see Section 2). 1 If in the power expansion of F0(x) there also exists a term x , then it would contribute to Veff 4 3 ζ ∼ N2α3v4 a v term of the form µN ( N ) z7 . This is the same expression as graviton interactions. The details of the coefficients and whether such terms exist in the power expansion of course has to be seen by actually performing the integration, however as we will see later, the integrals do produce terms of the correct interactions, in the both the membrane and the graviton limit. Next we look at Fn for n>0. The whole argument in the last paragraph goes through, except 1 now every term is weighted by an extra factor of Nn . For example, the membrane like interaction 4 produced by F looks like 1 αv . This factor of 1 means that this term is in fact a matrix n Nn µ2z5 Nn theory correction to supergravity, because it vanishes as N →∞. Therefore, we could see that in converting the summation into a series of integrals, only the first one is needed for comparison with supergravity. All the other Fn with n>0 produces only matrix theory corrections which is not the subject of interest here. Now we go back to Eq. (115). As discussed above, we only need the integral part F0, which we calculate using MATHEMATICA. After calculating the integrals, u is replaced by the supergravity = αµ variable w 3 u, and the answer is expanded first in large N, keeping only the leading order (which is order N 0), and then expanded in small α, keeping only the (α)1 order (which according to Section 2.2 is the appropriate order to be compared to supergravity in the membrane limit). Finally the answer is converted back into Minkowski signature by sending v2 →−v2. We obtain the following one-loop potential in the membrane limit4:
(36v2 − z2w2)(108v2 − (4w2 + 7z2)µ2) V = α . (120) eff 1152(w2 + z2)5/2µ2
4 This is only the effective potential in the membrane limit because when we expand in large N and keeping only the lowest order we essentially send the radius of the spheres r ∼ αµN to infinity, thus in terms of z we have z 1. Note 0 r0 also that the order the limits are taken is important. If the small α limit is taken before the large N limit, implicitly we z would be assuming αµ N which is the graviton limit. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 25
We find exact agreement when this expression is compared with the supergravity light-cone Lagrangian given in Eq. (95).
6. Interpolation of the effective potential
6.1. The interpolation
The purpose of this section is to find the effective potential that interpolates between the mem- brane limit and the graviton limit. Due to the complexity of the field equations on the supergravity side, this problem could only be attacked on the matrix theory side. On the matrix theory side, however, there is a subtlety that needs to be taken into account before such a potential could be found. Here we will analyze only the v4 term, the µ2v2 term as well as µ4 term can be studied in exactly the same way. From the supergravity side, what we wish to find is more or less clear. Near the membrane, when z 1, we expect an expansion like: r0
αv4 z z 2 z 3 z 4 V = + + + + +··· eff 2 5 1 (121) µ z r0 r0 r0 r0 and far away from the membrane, when z 1, we expect an expansion: r0
αv4 r2 r r 2 r 3 r 4 V = 0 1 + 0 + 0 + 0 + 0 +··· . (122) eff µ2z5 z2 z z z z
2 3 4 4 r2 We have used N α v = αv 0 to rewrite the graviton result so that it looks similar to the mem- z7 µ2z5 z2 brane effective potential. Therefore, we are basically looking for a function C(x) which appears in the effective potential in the following way:
αv4 z V = C . eff 2 5 (123) µ z r0 As one can see, both the graviton and the membrane actions are of this form. C(x) should have the appropriate limit at x → 0 and x →∞to give the correct potential at the membrane and the graviton limit respectively. Physically it represents curvature corrections due to the finite size of the spherical membrane. So now we go to the matrix theory side to try to find C(x). The subtlety is that the effective potential on the matrix theory side not only includes curvature corrections but also the matrix theory corrections to supergravity which we are not interested in, not to mention that the sum over mass formula is incapable of deducing such matrix theory corrections exactly [2]. For our purpose of comparing with supergravity, therefore, all matrix theory corrections to supergravity should be thrown away. Such corrections appear on the matrix theory side as 1/N corrections, mixed together with the curvature corrections, and it looks something like:
αv4 z 1 z 1 z V = C + C + C +··· . eff 2 5 0 1 2 2 (124) µ z r0 N r0 N r0
In such an expansion, only the C0 term should be kept. The reader is cautioned that naively sending N to infinity will not give us the correct interpolating potential because r0 also depends on N and such limit would only result in the effective potential in the membrane limit. 26 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
There are many ways such matrix corrections could appear. For example, in a typical matrix theory computation we may get terms of the form: αv4 V ∼ z2 + α2µ2 . (125) eff µ2z7 Rewriting the above gives:
αv4 α2µ2N 2 1 αv4 r2 1 V ∼ 1 + ∼ 1 + 0 . (126) eff µ2z5 z2 N 2 µ2z5 z2 N 2 The second term can now be identified as a matrix theory correction and is irrelevant to us. To isolate the curvature corrections (which we want) from the matrix theory corrections (which we do not want), we look at Eq. (124) more carefully. We can see that since 1 = αµ ,all N 6r0 the matrix theory corrections will appear at higher order in α. Therefore to get the interpolating effective potential from the matrix theory side, we can follow the steps below:
(1) change the summation over j in Eq. (115) into an integral over j from 0 to N; 6r0 3w (2) replace N by αµ , and u by αµ; (3) expand in small α and keep only the lowest order (which shall turn out to be order α1, with all lower orders vanishing). This is the interpolating effective potential.
With MATHEMATICA, the interpolating effective potential for two spheres of the same radius (w = 0) in Minkowski signature is found to be: α 36v2 − z2µ2 Veff = µ2z5 1152(4r2 + z2)5/2 ! " 0 # × 2 − 5 + 4 2 + 2 + 2 2 2 + 2 + 4 2 + 2 108v z 16r0 4r0 z 8r0 z 4r0 z z 4r0 z " " # − 2 2 4 2 + 2 + 4 − + 2 + 2 z µ 112r0 4r0 z 7z z 4r0 z " ##$ + 2 2 − + 2 + 2 8r0 z 2z 7 4r0 z . (127)
2 2 2 We see that Veff always carries a factor of 36v − z µ , meaning the effective potential van- = µz ishes whenever v 6 . This is expected because such configurations correspond to circular orbits which preserve half of the supersymmetries. Expanding this potential in the membrane limit of large r0 we get: α(3888v4 − 360v2z2µ2 + 7z4µ4) α(36v2 − z2µ2)(108v2 − 7z2µ2) V = = . (128) eff 1152µ2z5 1152µ2z5 This result is of course identical to the matrix theory result (with w = 0) in Section 5 where the membrane limit was taken in advance. In the graviton limit of small r0, after replacing r0 by αµN/6, we have: N 2α3(720v4 − 56v2z2µ2 + z4µ4) N 2α3(20v2 − z2µ2)(36v2 − z2µ2) V = = . (129) eff 768z7 768z7 The two limits of the effective action can then be compared with that of the supergravity side. Indeed from the expressions (96) and (97) we see that we have perfect agreement. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 27
In above we have only given the expression of the interpolating potential for two membranes with the same radius (w = 0). We have also found the interpolating potential with w included, using the same steps given above. However we choose to omit the rather lengthy expression here for brevity.
6.2. Comparison with Shin and Yoshida
As mentioned in Section 5, Ref. [10], considered the case of supersymmetric circular motion µ with angular frequency 6 and found a flat potential, which agrees with what we have found (see the comment under Eq. (127)). In a subsequent paper, [11], the authors considered the case of a slightly elliptical orbit with separation z, µt µt µt z = (r + )cos2 + (r − )sin2 = r2 + 2 + 2r cos (130) 2 6 2 6 2 2 3 and velocity v, µ µt µt v = (r + )sin2 + (r − )cos2 6 2 6 2 6 µ µt = r2 + 2 − 2r cos , (131) 6 2 2 3 where is the small expansion parameter for the eccentricity of the orbit. They considered the large separation limit, r2 0, and found an effective action, Eq. (1.2) of [11], 2 4 3 4 2 35 1 385(αµN) 1 1 Γeff = dt α µ N − 4 − 27 · 3 r7 211 · 33 N 2 r9 2 2 35 1 11(αµN)2 1 1 = 4 dt α3µ4N 2 − − 7 1 2 9 (132) 384 r2 36 4N r2 4 9 after expanding to O( ) and O(1/r2 ). Note that in the equation above we have restored µ and α, and set N1 = N2 = N, r1 = 0. To compare this with our result we substitute the above expressions of z(t) and v(t) into our interpolating potential equation (127) and expand in the parameters 1/r2 and . As we are only comparing effective actions we average our potential over one period of oscillation. We find for our time-averaged potential,
105 1 r2 35 1 11(αµN)2 1 V = α4µ2 r2 − 0 = α34µ4N 2 − eff 0 7 11 9 7 9 (133) 32 r2 r2 384 r2 36 r2 = αµN (where to reach the final line we used r0 6 ) which agrees with Eq. (132) after throwing away the 1 matrix theory corrections in the latter. In calculating the matrix theory effective potential N2 in Section 5 we assumed a constant velocity. However, as we can see by this comparison, as long as we ignore matrix theory corrections, it leads to the correct result. Thus we see that we can consistently neglect acceleration terms in the effective potential as discussed earlier. 28 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
7. Discussion and future directions
One immediate generalization of the work reported in this paper is to consider more compli- cated probe configurations (recall that in this paper we have restricted our attention to a probe membrane which is spherical and has no velocity in the x1 through x3 directions). For exam- ple, we could consider deforming the probe membrane so that it is no longer a perfect sphere. It means that the coordinates of the membrane XA will now be some general functions of θ and φ, and in particular it no longer has to appear only as a point in the 4 through 9 directions. We can also give the probe membrane a nonzero velocity in the x1 through x3 direction. These general- izations give more interesting dynamics and can be fairly easily carried out. On the supergravity side, this just requires putting the relevant probe configuration into the light-cone Lagrangian; on A the matrix theory side, this requires replacing the background configurations B(2)’s used in this paper with some more general configurations. In this paper we take the source to be a static spherical membrane preserving all of the sixteen linearly realized supersymmetries. As we know, there are other interesting membrane configu- rations that are permitted by the pp-wave background, for example, the static hyperbolic brane which preserves eight supersymmetries, given by [19]. It would be interesting to investigate gauge/gravity duality in the case where these static objects are sources and interact with some graviton or membrane probe. On the supergravity side, we can apply our general formalism of diagonalizing the field equations as well as the near-membrane expansion to find the metric and three-form perturbations produced by these sources. On the matrix theory side, for the aforemen- tioned hyperbolic membrane, one could use the unitary representations of SO(2, 1) worked out in [20–22] to expand the fluctuating fields and compute the one-loop effective potential. Moreover, it is evidently desirable to push our computation on the supergravity side beyond the near-membrane expansion, finding the solution to the field equations which will then give the δLlc to be compared with the fully interpolating potential (127) found on the matrix theory side. The membrane-interaction we described in this paper has no longitudinal momentum transfer. However, as we can see, our supergravity computation can be quite easily generalized to nonzero longitudinal momentum transfer, by not setting k− to zero when solving the field equations (e.g., see the comments made after Eq. (77)). On the gauge theory side, this requires an instanton computation. We hope to say more on this in the near future. There are also M-theory pp-wave backgrounds with less supersymmetries, and matrix theories in these pp-waves backgrounds have been proposed [23,24]. It would be interesting to investigate the gauge/gravity duality in these less supersymmetric settings. Looking slightly further afield it would be interesting to consider the extension of our results to the tiny graviton matrix theory. The BMN matrix theory used in this paper can be derived from the regularized M2-brane action, specifically the M2-brane in the 11-dimensional pp-wave background, and we may interpret it as a theory of N spherical M2-branes of very small size. The tiny graviton matrix theory of [25] is quite analogous, however instead of spherical M2-branes it describes the dynamics of very small D3-branes on the ten-dimensional pp-wave background. This is then conjectured to be a dual description of DLCQ strings in this background. Evidence for this conjecture was provided in the original paper [25] and subsequently in [26].In[25] the zero energy vacuum configurations were found and it was shown that there are solutions consisting of concentric fuzzy spheres much like those analyzed in this paper (though of course they are now fuzzy three-spheres). It should be possible to modify our calculations to study the analogous dynamics of spherical branes in the tiny graviton matrix theory and compare with what is found from the type IIB supergravity. H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30 29
Finally we note that while this paper was being written, the work by Lin, Lunin, and Mal- dacena (LLM) [27] appeared, in which all 1/2 BPS excitations of M-theory plane wave were found in a quite elegant manner. We expect that the singular supergravity solution for the spher- ical membrane presented in this work can be obtained from a certain regular LLM geometry by taking the appropriate limit and reducing to the linearized level, and it would be interesting to see how this can be done in detail.5 Moreover, one should be able to use LLM’s results to com- pute the interaction between a probe (e.g. a graviton, or a spherical membrane) and an arbitrary 1/2 BPS source, and compare the answer with that of a loop computation in the pp-wave matrix theory. However, how to represent these general 1/2 BPS sources on the matrix theory side is not very clear to us at this point.
Acknowledgements
The authors would like to thank John Schwarz for helpful advice and support during the sum- mer. We would also like to thank I. Swanson for useful discussions and Oleg Lunin for explaining the work [27] to us. This work was supported by DOE grant DE-FG03-92ER40701. X. Wu is also supported in part by DOE grant DE-FG01-00ER45832 and NSF grant PHY/0244811.
References
[1] D. Berenstein, J.M. Maldacena, H. Nastase, Strings in flat space and pp waves from N = 4 super-Yang–Mills, JHEP 0204 (2002) 013, hep-th/0202021. [2] H.K. Lee, X. Wu, Two-graviton interaction in pp-wave background in matrix theory and supergravity, Nucl. Phys. B 665 (2003) 153, hep-th/0301246. [3] K. Becker, M. Becker, A two-loop test of M(atrix) theory, Nucl. Phys. B 506 (1997) 48, hep-th/9705091. [4] K. Becker, M. Becker, On graviton scattering amplitudes in M-theory, Phys. Rev. D 57 (1998) 6464, hep- th/9712238. [5] K. Becker, M. Becker, Complete solution for M(atrix) theory at two loops, JHEP 9809 (1998) 019, hep-th/9807182. [6] K. Becker, M. Becker, J. Polchinski, A.A. Tseytlin, Higher order graviton scattering in M(atrix) theory, Phys. Rev. D 56 (1997) 3174, hep-th/9706072. [7] Y. Okawa, T. Yoneya, Multi-body interactions of D-particles in supergravity and matrix theory, Nucl. Phys. B 538 (1999) 67, hep-th/9806108. [8] Y. Okawa, T. Yoneya, Equations of motion and Galilei invariance in D-particle dynamics, Nucl. Phys. B 541 (1999) 163, hep-th/9808188. [9] Y. Okawa, Higher-derivative terms in one-loop effective action for general trajectories of D-particles in matrix theory, Nucl. Phys. B 552 (1999) 447, hep-th/9903025. [10] H. Shin, K. Yoshida, One-loop flatness of membrane fuzzy sphere interaction in plane-wave matrix model, Nucl. Phys. B 679 (2004) 99, hep-th/0309258. [11] H. Shin, K. Yoshida, Membrane fuzzy sphere dynamics in plane-wave matrix model, hep-th/0409045. [12] J.H. Park, Supersymmetric objects in the M-theory on a pp-wave, JHEP 0210 (2002) 032, hep-th/0208161. [13] K. Becker, M. Becker, J. Polchinski, A.A. Tseytlin, Higher order graviton scattering in M(atrix) theory, Phys. Rev. D 56 (1997) 3174, hep-th/9706072. [14] X. Wu, PhD thesis, available at http://etd.caltech.edu/etd/available/etd-05252004-230238/. [15] T. Kimura, K. Yoshida, Spectrum of eleven-dimensional supergravity on a pp-wave background, Phys. Rev. D 68 (2003) 125007, hep-th/0307193. [16] J. Polchinski, P. Pouliot, Membrane scattering with M-momentum transfer, Phys. Rev. D 56 (1997) 6601, hep- th/9704029. [17] K. Dasgupta, M.M. Sheikh-Jabbari, M. Van Raamsdonk, Matrix perturbation theory for M-theory on a PP-wave, JHEP 0205 (2002) 056, hep-th/0205185.
5 Note added: see [28] for the computation of the membrane giant graviton supergravity solution at the full nonlinear order. 30 H.K. Lee et al. / Nuclear Physics B 728 (2005) 1–30
[18] S.R. Das, J. Michelson, A.D. Shapere, Fuzzy spheres in pp-wave matrix string theory, Phys. Rev. D 70 (2004) 026004, hep-th/0306270. [19] D. Bak, Supersymmetric branes in PP wave background, Phys. Rev. D 67 (2003) 045017, hep-th/0204033. [20] L.J. Dixon, M.E. Peskin, J. Lykken, N = 2 superconformal symmetry and SO(2, 1) current algebra, Nucl. Phys. B 325 (1989) 329. [21] J.M. Maldacena, H. Ooguri, Strings in AdS(3) and SL(2,R) WZW model. I, J. Math. Phys. 42 (2001) 2929, hep- th/0001053. [22] B. Morariu, A.P. Polychronakos, Quantum mechanics on noncommutative Riemann surfaces, Nucl. Phys. B 634 (2002) 326, hep-th/0201070. [23] M. Cvetic,ˇ H. Lü, C.N. Pope, M-theory pp-waves, Penrose limits and supernumerary supersymmetries, Nucl. Phys. B 644 (2002) 65, hep-th/0203229. [24] N. Iizuka, Supergravity, supermembrane and M(atrix) model on pp-waves, Phys. Rev. D 68 (2003) 126002, hep- th/0211138. [25] M.M. Sheikh-Jabbari, Tiny graviton matrix theory: DLCQ of IIB plane-wave string theory, a conjecture, JHEP 0409 (2004) 017, hep-th/0406214. [26] M.M. Sheikh-Jabbari, M. Torabian, Classification of all 1/2 BPS solutions of the tiny graviton matrix theory, JHEP 0504 (2005) 001, hep-th/0501001. [27] H. Lin, O. Lunin, J. Maldacena, Bubbling AdS space and 1/2 BPS geometries, JHEP 0410 (2004) 025, hep- th/0409174. [28] D. Bak, S. Siwach, H.U. Yee, 1/2 BPS geometries of M2 giant gravitons, hep-th/0504098. Nuclear Physics B 728 (2005) 31–54
Low-energy photon–photon collisions to two loops revisited
J. Gasser a,M.A.Ivanovb,M.E.Sainioc,d
a Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland b Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna (Moscow region), Russia c Helsinki Institute of Physics, PO Box 64, 00014 University of Helsinki, Finland d Department of Physical Sciences, University of Helsinki, Finland Received 27 June 2005; accepted 1 September 2005 Available online 26 September 2005
Abstract In view of ongoing experimental activities to determine the pion polarizabilities, we have started to recalculate the available two-loop expressions in the framework of chiral perturbation theory, because they have never been checked before. We make use of the chiral Lagrangian at order p6 now available, and of improved techniques to evaluate the two-loop diagrams. Here, we present the result for the neutral pions. The cross section for the reaction γγ → π0π0 agrees with the earlier calculation within a fraction of a percent. We present analytic results for the dipole and quadrupole polarizabilities, and compare the latter with a recent evaluation from data on γγ → π0π0. 2005 Elsevier B.V. All rights reserved.
PACS: 11.30.Rd; 12.38.Aw; 12.39.Fe; 13.60.Fz
Keywords: Chiral perturbation theory; Two-loop diagrams; Pion polarizabilities; Compton-scattering
1. Introduction
We consider the process γγ → π 0π 0 in the framework of chiral perturbation theory (ChPT) [1–3]. The one-loop calculation of the scattering amplitude was performed in Refs. [4,5], and the two-loop amplitude was worked out in [6]. Because the effective Lagrangian at order p6 was not available at that time, the ultraviolet divergences were evaluated in the MS scheme,
E-mail address: [email protected] (J. Gasser).
0550-3213/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.nuclphysb.2005.09.010 32 J. Gasser et al. / Nuclear Physics B 728 (2005) 31–54 then dropped and replaced with a corresponding polynomial in the external momenta. The three new counterterms which enter at this order in the low-energy expansion were estimated with resonance saturation. Whereas such a procedure is legitimate from a technical point of view, it does not make use of the full information provided by chiral symmetry. The evaluation of the two-loop amplitude involving charged pions was performed later by Burgi [7]. Over the last ten years, considerable progress has been made in this field, both in theory and experiment. As for theory, the Lagrangian at order p6 has been constructed [8,9], and its divergence structure has been determined [10]. This provides an important check on the above calculations: adding the counterterm contributions from the p6 Lagrangian to the MS amplitude evaluated in [6] and in [7] must provide a scale independent result. Also in the theory, improved techniques to evaluate the two-loop diagrams that occur in these amplitudes have been developed [11]. The improvement arises mainly in the evaluation of diagrams with four external legs, where the techniques of Ref. [11] allow one to extract the ultraviolet divergences by use of simple recursion relations. We are now able to present the final result for the two-loop amplitudes in a rather compact form (in Refs. [6,7], the result was presented partly in numerical form only, because the algebraic expressions were too long to be published). Concerning experiment, quadrupole polarizabilities [12] for the neutral pions have recently been determined from data on γγ → π 0π 0 [13]. Further, the charged pion polarizabilities (α − β)π+ have been determined at the Mainz Microtron MAMI [14], with a result that is at variance with the two-loop calculation presented in [7]. Last but not least, there is an ongoing experiment by the COMPASS Collaboration at CERN to measure the charged pion and kaon polarizabilities [15,16]. In view of these developments, and because the two-loop expressions for the polarizabilities had never been checked, we decided to recalculate these amplitudes, using the improved tech- niques of Ref. [11] to evaluate the integrals, and invoking the chiral Lagrangian at order p6 [9, 10]. As the calculation in the case of neutral pions involves considerably less diagrams, and be- cause the Fortran code for these amplitudes is still available to us for checks, we have decided to start the program with a re-evaluation of these amplitudes. This is the main purpose of the present work. The evaluation of the corresponding expressions for the charged pions and for the kaons is underway and will be presented elsewhere [17]. The article is organized as follows. Section 2 contains the necessary kinematics of the process γγ → π 0π 0. To make the article self contained, we summarize in Section 3 the necessary ingre- dients from the effective Lagrangian framework. In Section 4, we display the Feynman diagrams and discuss their evaluation. Section 5 contains a concise representation of the two Lorentz in- variant amplitudes that describe the scattering matrix element. In Section 6, we compare the present work with the previous calculation [6], while Section 7 contains explicit expressions for the dipole and quadrupole polarizabilities valid at next-to-next-to-leading order in the chiral ex- pansion, together with a numerical analysis and a comparison with an evaluation from data on γγ → π 0π 0 [13]. The summary and an outlook are given in Section 8. Finally, several technical aspects of the calculation are relegated to Appendices A–E.
2. Kinematics
The matrix element for the reaction
0 0 γ(q1)γ (q2) → π (p1)π (p2) (2.1) J. Gasser et al. / Nuclear Physics B 728 (2005) 31–54 33 is given by 0 0 4 (4) N π (p1)π (p2) out γ(q1)γ (q2) in = i(2π) δ (Pf − Pi)T , (2.2) with T N = e2µνV , 1 2 µν −i(q1x+q2y) 0 0 Vµν = i dxe π (p1)π (p2) out Tjµ(x)jν(y)|0. (2.3)
2 Here jµ is the electromagnetic current, and α = e /4π 1/137. We consider real photons, 2 = · = qi 0, with i qi 0. The decomposition of the correlator Vµν into Lorentz invariant ampli- tudes reads
Vµν = A(s, t, u)T1µν + B(s,t,u)T2µν + C(s,t,u)T3µν + D(s,t,u)T4µν, 1 T = sg − q q , 1µν 2 µν 1ν 2µ 2 T2µν = 2s∆µ∆ν − ν gµν − 2ν(q1ν∆µ − q2µ∆ν),
T3µν = q1µq2ν,
T4µν = s(q1µ∆ν − q2ν∆µ) − ν(q1µq1ν + q2µq2ν),
∆µ = (p1 − p2)µ, (2.4) where 2 2 2 s = (q1 + q2) ,t= (p1 − q1) ,u= (p2 − q1) ,ν= t − u (2.5) are the standard Mandelstam variables. The tensor Vµν satisfies the Ward identities µ = ν = q1 Vµν q2 Vµν 0. (2.6) The amplitudes A and B are analytic functions of the variables s,t and u, symmetric under crossing (t, u) → (u, t). The amplitudes C and D do not contribute to the process considered here, because i · qi = 0. It is useful to introduce in addition the helicity amplitudes 4 − 2 8(Mπ tu) H++ = A + 2 4M − s B, H+− = B. (2.7) π s The helicity components H++ and H+− correspond to photon helicity differences λ = 0, 2, re- | = 3 0 (3) − spectively. With our normalization of states p1 p2 2(2π) p1δ (p1 p2), the differential cross section for unpolarized photons in the centre-of-mass system is
γγ→π0π0 2 dσ α s 2 2 = β(s)H(s,t), H(s,t) =|H++| +|H+−| , (2.8) dΩ 64 = − 2 with β(s) 1 4Mπ /s. The relation between the helicity amplitudes M+± in Ref. [13] and the amplitudes used here is
M++(s, t) = 2παH++(s, t), M+−(s, t) = 16παB(s,t). (2.9) The physical regions for the reactions γγ → π 0π 0 and γπ0 → γπ0 are displayed in Fig. 1, where we also indicate with dashed lines the nearest singularities in the amplitudes A and B. These singularities are generated by two-pion intermediate states in the s,t and u channel. 34 J. Gasser et al. / Nuclear Physics B 728 (2005) 31–54
Fig. 1. Mandelstam plane with three related physical regions. s-channel: γγ → ππ, t-andu-channel: γπ → γπ.We → → = 2 denote the threshold for γγ ππ(γπ γπ) by A (P ). The dashed lines at s,t,u 4Mπ indicate the presence of branch-points in the amplitude, generated by two-pion intermediate states.
3. The effective Lagrangian and its low energy constants
The effective Lagrangian consists of a string of terms. Here, we consider QCD with two flavours, in the isospin symmetry limit mu = md =ˆm. At next-to-next-to-leading order (NNLO), one has
Leff = L2 + L4 + L6. (3.1)
The subscripts refer to the chiral order. The expression for L2 is
F 2 L = D UDµU † + M2 U + U † , 2 4 µ e D U = ∂ U − i(QU − UQ)A ,Q= diag(1, −1), (3.2) µ µ µ 2 where e is the electric charge, and Aµ denotes the electromagnetic field. The quantity F denotes the pion decay constant in the chiral limit, and M2 is the leading term in the quark mass expansion 2 2 = 2 + ˆ ··· of the pion (mass) , Mπ M (1 O(m)). Further, the brackets denote a trace in flavour space. In Eq. (3.2), we have retained only the terms relevant for the present application, i.e., we have dropped additional external fields. We choose the unitary 2 × 2matrixU in the form √ 2 0 + 2 π √π 2π U = σ + iπ/F, σ + = 1 × ,π= − . (3.3) F 2 2 2 2π −π 0 The Lagrangian at NLO has the structure [2]