Spherical transverse M5-branes from the plane wave matrix model

We consider matrix theoretical description of transverse M5-branes in M-theory on the 11-dimensional maximally supersymmetric pp-wave background. We apply the localization to the plane wave matrix model (PWMM) and show that the transverse spherical fivebranes with zero light cone energy in M-theory are realized as the distribution of low energy moduli of the $SO(6)$ scalar fields in PWMM.


Introduction
Matrix models are conjectured to give nonperturbative formulations of M-theory [1]. This formulation is expected to realize a second quantization of M-theory, which contains all the fundamental objects in the theory. However, the description of states with M5-branes in the matrix models has not been established yet. Understanding this problem will shed light on the matrix-model formulation of M-theory.
In this paper, we focus on M-theory defined on the maximally supersymmetric ppwave solution of the 11-dimensional supergravity and consider the description of certain M5-branes living in this geometry in terms of the matrix model. On this background, there exist stable spherical M2-and M5-branes with zero light cone energy. According to the matrix-model conjecture, objects with zero light cone energy should be realized as vacuum states in the corresponding matrix model. Hence, these spherical branes should also be realized as certain vacuum states in the matrix model. In this paper, we investigate this relation in detail by using the localization method.
The matrix model for M-theory on the pp-wave background is called the plane wave matrix model (PWMM) [2]. This model is given by a mass deformation of the BFSS matrix model [1], where the mass parameter is proportional to the three form flux on the pp-wave geometry. Because of the mass deformation, PWMM possesses many discretely degenerate vacua, unlike the BFSS matrix model. The relation between these vacua and objects with vanishing light cone energy in M-theory was proposed in [2,3]. Here, in particular, the vacua corresponding to the above mentioned spherical M5-brane and its multiple generalization were also specified. For the case of a single M5-brane, this correspondence was tested by comparing the BPS protected mass spectra of PWMM with that of the M5-brane [3].
Let us review this proposal in more detail. The vacua of PWMM, which preserve all the supersymmetry, are given by the fuzzy sphere [4] and are labeled by N -dimensional representations of the SU (2) Lie algebra, where N is the matrix size of PWMM. Generally, the classical vacuum configuration in PWMM takes the form of Any N -dimensional representation gives a supersymmetric vacuum and, in general, the representation is reducible. Then, one can make an irreducible decomposition: (1.2) Here, L are fixed while n s go to infinity, smooth two-spheres are realized from the fuzzy sphere. One can naturally expect that this smooth sphere is the spherical M2-brane with zero light cone energy.
On the other hand, in [2], the spherical M5-brane was conjectured to be realized as the trivial vacuum of PWMM, where all the fields are vanishing. This is the case where the representation in (1.1) is a direct sum of N trivial representations. Furthermore, the conjecture was generalized to the case of multiple spherical M5-branes [3]. In these conjectures, the M5-branes are considered to be realized in the limit such that n s are fixed and N (s) 2 go to infinity in (1.2). In order to describe this limit more precisely, let us introduce Young diagrams associated with the partition of (1.2). In the decomposition (1.2), we assume that n 1 > n 2 > · · · > n Λ without loss of generality. Then we consider a Young diagram which consists of where R is the radius of the light like circle 2 . Note that the total light cone momentum is given by p + = Λ s=1 (n s − n s+1 )p + s and this is equal to N/R. Note also that N 5 := max{n s |s = 1, 2, . . . } = n 1 corresponds to the total number of M5-branes. This conjecture is highly nontrivial. For example, let us consider the simplest partition with Λ = 1, n 1 = 1, N (1) 2 = N , which corresponds to the trivial vacuum of PWMM.
At the classical level, the vacuum configuration is just vanishing, so that we can not see any structure of the M5-brane. For example, it looks seemingly impossible to reproduce geometric information of the spherical M5-brane (the radius etc.) from the trivial configuration. Nevertheless, the conjecture claims that a single spherical M5-brane is realized in the trivial vacuum. 1 For s = Λ, we define n Λ+1 := 0. 2 As we will see in the next section, the radius of a single (i.e. not coincident) M5-brane is proportional to (p + ) 1/4 . Thus, larger p + gives a larger radius. Though this relation had never been derived for coincident M5-branes, our results discussed below shows that this is also true for coincident M5-branes. Fig. 1 is based on this picture, so that the sth stack has a larger radius than (s − 1)th stack.
To bridge this gap, one needs to recall that M-theory is conjectured to be realized in an appropriate large-N limit of PWMM, where the coupling constant also becomes very large as the matrix size N goes to infinity. Thus, one has to deal with the strongly coupled regime of PWMM, in order to understand the description of M5-branes. In the strong coupling region, there must be a large quantum fluctuation around the classical vacuum configuration. Thus, typical configurations of matrices will be very different from the classical configuration. There is a possibility that the spherical M5-branes are formed as a typical configuration of matrices in the strong coupling region of PWMM 3 .
In this paper, we investigate this possibility by directly studying the strong coupling regime of PWMM. The limit we consider is for any s, t = 1, 2, · · · , Λ. This limit corresponds to Λ stacks of M5-branes with different radii as shown in Fig. 1. In addition, we also scale the coupling constant of PWMM in such a way that the M5-branes decouple with the bulk gravity and only the degrees of freedom on the M5-branes become relevant [3]. This decoupling limit turns out to be the strong coupling limit in the 't Hooft limit of PWMM, as we will describe in the next section. In this decoupling limit, we apply the localization to PWMM and reduce some BPS correlation functions to certain eigenvalue integrals. By evaluating the eigenvalue integral, we argue that the eigenvalue distribution of the low energy modes of the SO (6) scalar fields forms Λ stacks of spherical shells and coincides with the expected configuration of the spherical M5-branes in M-theory 4 . In particular, we show that, for a single M5brane, the radius of the shell completely agrees with the value computed by using the classical Dirac-Nambu-Goto action of a single M5-brane. This result strongly supports the proposal of [3] and shows that PWMM indeed contains the multiple M5-brane states.
We also apply the same argument to M2-branes and show that the spherical M2-brane can be described in a similar way using the eigenvalue integral.
This paper is organized as follows. In section 2, we review M-theory on the pp-wave background. We show that there exist spherical M2-and M5-branes with zero light cone 3 See also [5] for the description of M5-branes in a different matrix model. 4 A part of this result was briefly reported in the letter [6] for the case of concentric M5-branes. In this paper, we not only describe the technical details of [6] but also generalize the result of [6] to the most general configurations of the spherical M5-branes. energy on this background. We also compute the radii of these objects. In section 3, we review PWMM. In section 4, we apply the localization to PWMM and evaluate the moduli distribution of scalar fields. We show that the distribution agrees with the configuration of the spherical M5-branes. In section 5, we consider the case of M2-branes. In section 6, we summarize our results and discuss the low energy theory of PWMM.

M-theory on the pp-wave background
In this section, we review M-theory on the maximally supersymmetric plane wave background in the 11-dimensional supergravity. The background geometry is given by x a x a dx + dx + , where µ is the flux parameter of the three form field 5 . We will see that spherical M2-brane and M5-brane exist as the lowest energy states with respect to the light cone Hamiltonian.
We refer the method in [7] for the calculation in this section.

Spherical M2-brane
We first consider a single M2-brane in the background (2.1). The bosonic part of the M2-brane action is given by the Dirac-Nambu-Goto action plus a Chern-Simons term as Here, h αβ is the induced metric, for the embedding function X µ (σ). The overall constant T M2 in (2.2) is the tension of M2-brane given by where l p stands for the Planck length. By introducing a symmetric auxiliary field γ αβ , we rewrite the action into the Polyakov type: This action has a diffeomorphism symmetry for the worldvolume coordinates σ α = (σ 0 , σ 1 , σ 2 ) of the membrane. If we consider an M2-brane with topology R × Σ, where R is the time direction and Σ is a Riemann surface, we can fix this symmetry by putting where a, b = 1, 2 and the determinant is taken in this 2 × 2 subspace. ν is a constant which will be related to the light cone momentum of the M2-brane below. Then, the action becomes Here, in the second line, we have introduced a canonical Poisson bracket on the membrane defined by {f, g} = ab (∂ a f )(∂ b g) for each fixed σ 0 . In terms of the Poisson bracket, the Chern-Simons term can be written as The gauge fixing condition (2.6) as well as the equation of motion of the auxiliary field produce the following constraints: From the second constraint, it also follows that {g µνẊ µ , X ν } = 0.
(2.10) Thus, the system is reduced to the theory (2.7) with these constraints imposed.
The constraints (2.9) can be explicitly solved in the light cone gauge, Here, we have defined X ± by We then consider the Hamilton formalism. We denote by P µ the canonical conjugate momentum of X µ . The total light cone momentum is then given by where we have chosen the spacial coordinates such that they have a volume d 2 σ = 4π.
This relates the constant ν to the light cone momentum. The Hamiltonian is given by where, V 2 is the volume of the unit sphere, V 2 = 4π. The remaining constraint (2.10) is written in terms of the transverse components X A as Now, let us consider a vacuum configuration, which minimizes the Hamiltonian (2.14).
Note that the potential for X i forms a perfect square, (2.16) From this, we find that the vacuum configuration is given by where x i are the embedding function of the unit sphere in R 3 satisfying The radius is also determined as The configuration (2.17) obviously has the spherical shape. Thus we see that, in M-theory on the pp-wave background, there exists a spherical zero energy M2-brane with the radius given by (2.19).

Spherical M5-brane
Then, let us consider a single M5-brane. We start from the bosonic part of the action, where dC 6 = * F 4 and the tension is written as (2.21) We can apply the computation in the previous subsection to (2.20). Then, we can obtain the light-cone Hamiltonian for the M5-brane, Here, V 5 is the volume of the unit 5-dimensional sphere, V 5 = π 3 . The curly bracket with five entries in (2.22) is the 5-dimensional analogue of the Poisson bracket defined by We notice that the potential terms of X a forms a perfect square, (2.24) Thus, we find that the vacuum configuration is given by a spherical fivebrane of the form, where x a are the embedding function of the unit 5-sphere into R 6 satisfying x a x a = 1, {x a 1 , · · · , x a 5 } = a 1 a 2 ···a 6 x a 6 . (2.26) The radius of the fivebrane is determined as (2.27)

Decoupling limits
In this paper, we focus on the limits in which the radii of the spherical M2-and M5branes become very large and only the degrees of freedom on these branes survive for low energy physics [3].
Let us introduce the radius r = √ x i x i of the two sphere on which the M2-brane is wrapping. The metric (2.1) is written as where · · · represents the other terms which are irrelevant in this discussion and we have where R is the radius of the original compactified circle along the light-like direction. Since the shift ofx + is much smaller than that ofx − in the large-r limit, this can be effectively regarded as a spatial compactification near the large M2-brane 6 . From the structure of the metric (2.28), we find that the physical radius of the M-circle is given byR ∼ R/(µr).
In the perspective of the type IIA superstring theory, the spherical M2-brane wrapping on the two-sphere in (2.28) corresponds to a D2-brane. The gauge coupling constant on D2-branes is given by g 2 Y M ∼ g s l −1 s , where g s and l s are the string coupling and the string length. By translating this into the M-theory parameters using the standard dictionary, g s ∼ (R/l p ) 3/2 and l s ∼ (l 3 p /R) 1/2 , one can express the coupling constant as g 2 Y M ∼R 2 /l 3 p . In the limit where the radius of the D-branes becomes large, it is convenient to rescale the metric, so that the parameter which controls the theory on D2-brane is given by the dimensionless coupling constant g 2 Y M r M2 , where r M2 is the radius of the M2-brane. By 6 One can also take another coordinate (x + ,x − ) such that the metric becomes canonical Minkowski metric. In this coordinate, the both shifts ofx ± are given by R/(µr) and this looks like a light cone compactification. However, note that from d dx + ∼ 1 µr d dx + , we see that the energy alongx + direction is given by 1 µr H. Similarly, the momentum alongx − direction is µrp + + 1 µr H. In the limit discussed below, both r and p + becomes large, so that the energy is much smaller than the spatial momentum. Thus, after all, this can be indeed regarded as a spatial compactification.
using (2.19), The coupling constant can be expressed as (2.30) We are interested in the case where l p and µ are fixed. Moreover, in order to have an interacting theory on the D2-branes in the r M2 → ∞ limit, we would like to fix the coupling constant (2.30). Then, the decoupling limit of the D2-branes is given by The fixed quantity in (2.31) measures the size of the M-circle for each fixed p + , so that the M2-brane in 11-dimension is realized in the limit where R 2 p + becomes large. The limit (2.31) can be written in terms of the parameters of the matrix model. The D0-brane charge (the matrix size) N is related to the light cone momentum by p + = N/R and the gauge coupling of D0-branes is given as g 2 ∼ R 3 l −6 p . Thus, the limit (2.31) is translated to (2.32) Then, the decoupling limit of M2-brane is given by sending g 2 N to infinity. Next, we consider the decoupling limit of the spherical M5-brane. The theory on NS5-branes is known as the little string theory. This theory is characterized by the string tension proportional to 1/l 2 s . We can apply the above argument for D2-branes to the little string theory. Here, the fixed quantity is replaced by the tension of the little string which is made dimensionless by using the radius of the M5-brane (2.27): (2.33) Thus, the decoupling limit of NS5-brane is given by The M5-branes in 11-dimension are realized by further taking R 4 p + to be large.
In terms of the parameters of the matrix model, the decoupling limit of NS5-brane This is just the 't Hooft limit of the matrix model 7 . The M5-brane limit corresponds to the strong coupling limit with respect to the 't Hooft coupling g 2 N .
For multiple M5-branes, the radius of each M5-brane should become large to decouple from the gravity. Furthermore, if there are some stacks of M5-branes with different radii as shown in Fig. 1, the distances between the nearest stacks should also become large. The limit realizing this situation is such that the all radii become large with the same order.
Since the radius of each M5-brane is proportional to a positive power of the light cone momentum, the decoupling limit for the multiple fivebranes should be given by (2.35) The plane wave matrix model In this section, we review the plane wave matrix model (PWMM) [2].
The Hamiltonian of PWMM is obtained by the matrix regularization of the Hamiltonian (2.14) of a single M2-brane [9]. In the matrix regularization, real functions on the world volume f (σ a ) are linearly mapped to N × N Hermitian matrices, in such a way that integrals and the Poisson algebra of functions are consistently mapped to traces and the commutator algebra of the corresponding matrices, respectively. Namely, under this mapping, we have For example, let us consider the case where the spatial world volume is a unit sphere embedded in R 3 . The image of the embedding function x i which satisfies (2.18) is given by the N -dimensional irreducible representation of the SU (2) generators, The normalization is chosen so that ix 2 i = 1 N holds in the large-N limit. One can check that (3.1) is satisfied by (3.2) for sufficiently large N .
By applying the matrix regularization to the Hamiltonian (2.14), we obtain the bosonic part of the Hamiltonian of PWMM as 9) are now N × N matrices, which correspond to the images of P A (σ a ) and X i (σ a ) in (2.14). The constraint (2.15) is replaced by In obtaining (3.3), we have also rescaled the momenta as P A → N 4π P A . The rescaled momenta correspond to the canonical momenta of X A in PWMM. When one quantizes the theory of M2-brane (2.14), one has the canonical commutation relation 8 , Without the rescaling, according to (3.1), this would be mapped to The rescaling just removes the factor N 4π on the right-hand side and makesP Aij the canonically normalized momenta ofX A ij . We consider vacua of PWMM. Noticing that the potential for X i forms a perfect square, we find that the Hamiltonian is minimized when For later convenience, we introduce the action of PWMM. We first rescale the matrices as Then, the bosonic action of PWMM can be written in a simple form as 9 Here, the coupling constant is related to the original parameters by and the covariant derivative is defined by The gauge field A is introduced to take the constraint (3.4) into account. In the A = 0 gauge, the Gauss law constraint reproduces (3.4).

Localization in PWMM
We consider a complex scalar field in PWMM defined by The real and imaginary parts of φ are given by an SO(3) scalar and an SO(6) scalar, respectively, up to the time dependent rotation. When one makes a double Wick-rotation for the time and Y 9 directions, one can construct (four) supercharges which leave φ invariant. This allows us to exactly compute the expectation values of operators made of only φ by using the localization method [10].
In order to perform the localization, one first needs to define the boundary conditions in the Euclidean time direction. Since we are interested in PWMM expanded around a 9 We have also rescaled the time coordinate appropriately.
fixed vacuum, the appropriate boundary condition is such that all the fields approach to the vacuum configuration as the Euclidean time goes to ±∞. With this boundary condition, the path integral of PWMM defines the theory around the fixed background.
For the theory around the generic vacuum (1.2), the result of the localization obtained in [11][12][13] is summarized below. See appendix A for the detail of localization. We have the following equality: Hermitian matrix. The expectation value · · · on the left-hand side of (4.2) is taken with respect to the original action of PWMM expanded around the background (1.2). On the other hand, the expectation value · · · M M on the right-hand side of (4.2) is taken with respect to the following matrix integral: 2 ) are eigenvalues of M s and Z 1−loop is the one-loop determinant, which arises in the 1-loop calculation of the localization. Z 1−loop is given by The prime on the last product means that the second factor in the numerator with s = t, J = 0 and i = j is not included in the product.
Note that the right-hand side of (4.2) does not depend on the time coordinates t a . So this relation implies that the correlator on the left-hand side does not depend on time.
This property can be understood from the SUSY Ward identity, as shown in [12].
We remark that, in the calculation of the localization, some possible instanton corrections are neglected [11][12][13]. This corresponds to kink-like configurations in PWMM which connect two distinct vacua [14][15][16][17]. However, the instanton amplitudes are bounded from below by N 2 /λ times the difference of the quadratic Casimirs of the two vacua. Thus, in the decoupling limit of the M5-brane, this effect is suppressed.

Coincident M5-branes from the simplest partition
To illustrate our computation, let us first consider the simplest partition with Λ = 1, namely the vacuum with According to the proposal in [3], this corresponds to N 5 coincident M5-branes. In this case, the eigenvalue integral (4.4) reduces to a one matrix model: In the decoupling limit of the M5-brane, N 2 becomes infinity, so that the saddle point approximation is valid in evaluating the eigenvalue integral (4.7). As usual, we introduce the eigenvalue distribution Here, q m represents the range of the support of ρ(x) 10 . Note that, we are interested in the decoupling limit of M5-brane where the 't Hooft coupling λ := g 2 N 2 goes to infinity. In this regime, the Gaussian attractive force of the eigenvalue integral (4.7) becomes weaker, so that q m is expected to go to infinity. If one considers the region where q m is very large compared to N 5 , one can reduce the saddle point equation of ρ(x) to where β is the Lagrange multiplier, which imposes the normalization (4.9). See appendix B for the derivation of (4.10) In the M5-brane limit, the solution to the saddle point equation is given by See appendix C.1 for the derivation of this solution 11 . Note that indeed q m becomes infinity as the 't Hooft coupling goes to infinity.
By using this solution, we can compute correlation functions of φ. For example, Note that the second term is much larger than the first term in the strong coupling regime with N 5 fixed. This originally comes from the fact that the eigenvalue distribution of M spreads over the much wider region than the distribution of L 3 in the M5-brane limit.
This property is common for any correlation function of φ, including the resolvent. In this regime, therefore, the imaginary part of φ is dominant and the real part is negligible.
In other words, the spectrum of φ lies along the imaginary axis in this limit.
Assuming that the matrices Y A become mutually commuting in the decoupling limit, one may expect that this spectrum on the imaginary axis given by ρ in (4.11) could be identified with the eigenvalue distribution of one of the SO(6) scalars. However, such identification would contradict with the discussion in [18] by Polchinski. In [18], the BFSS matrix model is considered and the trace of the square of the scalar fields Y A is shown to be bounded from below by λ 2/3 (in the notation used in this paper). And this conclusion is considered to hold also for PWMM if we assume the gauge/gravity correspondence: The dual geometry of PWMM reduces to the dual geometry of BFSS matrix model at a sufficiently large radius r O(λ 1/4 ) in the decoupling limit of M5-brane [19]. On the 11 See also [12] for another derivation using the Fermi gas method.
other hand, if one assumes that ρ in (4.11) gives the eigenvalue distribution of one of the SO(6) scalars in PWMM, this would give 1 N Tr(Y A ) 2 = O(λ 1/2 ), which is smaller than the bound in the Polchinski's argument. Thus, this leads to a contradiction and the first assumption that Y A become commuting in the decoupling limit seems to be wrong 12 .
Apart from Polchinski's argument, we can find another reasoning for the above statement, based on the gauge/gravity correspondence. The gravity dual [19] of PWMM has a typical scale λ 1/3 , which is the string scale beyond which the supergravity approximation is not valid. It is natural to expect that the matrix elements of Y A contain information of such typical scale on the gravity side, so that the scalar fields in PWMM have the Then, it is again suggested that the matrices are noncommuting even in the strongly coupled region.
The classical geometry of the supergravity and the M2/M5-branes are considered to be realized as the low energy moduli of these matrices. Roughly speaking, they will correspond to the low energy modes of the matrices and one needs to consider the low energy theory of the matrix model to find the classical geometric objects in M-theory.
The noncommuting modes, which produce the large value for 1 N Tr(Y A ) 2 = O(λ 2/3 ), have a large excitation energy, so that these modes should be frozen and irrelevant in studying the low energy theory.
Note that the complex field φ has the eigenvalue distribution of order of λ 1/4 . This is much smaller than the typical value of the noncommuting modes. From this fact, we find that φ is a good low energy field and the operators Trφ n can be considered as operators in the low energy theory. This can also be understood from our formula (4.2) of the localization. The correlation functions of φ are independent of the time coordinates and hence are invariant under taking the time averages, which projects the operators to the low energy modes (More specifically, one can eliminate the high energy modes by integrating over very short time intervals with length given by 1/C, where C is a constant much smaller than the typical energy scale for noncommuting modes but much larger than the energy scale for (4.11).). This means that the result of the localization (4.2) contains only the low energy modes.
As is discussed in [18], operators in the matrix model should be additively renormal- 12 We thank J. Maldacena for suggesting this problem and the resolution using the time average which we will discuss below. ized in the low energy theory, where the additive renormalization constants correspond to contributions from the high energy noncommuting modes. However, such additive renormalization is not needed for Trφ n . In order for the eigenvalues of φ to be of O(λ 1/4 ), the renormalization constants for Y A must cancel out in the correlators of φ. For example, this can be seen in our computation in (4.12). Since the Trφ 2 is given by the difference between Tr(Y 3 ) 2 and Tr(Y 9 ) 2 , the renormalization constants should cancel out.
The statement that φ picks up the low energy moduli of the matrices is also supported by the earlier work on the gauge/gravity correspondence for PWMM. It was shown in [12,13] that the field φ describes a system of moduli parameters on the gravity side, which is equivalent to a certain axially symmetric electrostatic system: The charge densities of the electrostatic system, which determines the geometry on the gravity side, were shown to be equivalent to the eigenvalue density of φ.
From these observations, we claim that the spectrum of φ is identified with the low energy moduli of PWMM. Furthermore, we claim that the low energy moduli in PWMM are given by commuting matrices in the decoupling limit. This can be understood as follows. Suppose that the moduli are given by noncommuting matrices and the theory on the M5-branes has some noncommutativity of the low energy moduli parameters as well as some length scale associated with the noncommutativity. The noncommutative length scale must be much smaller than the radius of the M5-brane, since otherwise the M5-brane would not be localized along the radial direction due to the nonlocality caused by the noncommutativity and hence would not be regarded as 1+5 dimensional object. Then, let us consider the length scale λ 1/4 of the low energy moduli computed from the localization. This scale corresponds to the scale of the M5-brane radius if one takes the rescaling (3.8) into account. Thus, the length scale of the low energy moduli must be much larger than the noncommutative scale. Therefore, even if the moduli have noncommutativity, this effect must be much smaller than the value of the moduli themselves in the decoupling limit. Thus, we can ignore the noncommutativity and can regard the moduli as just commuting matrices. Note that this conclusion is consistent with our result of the localization (4.2). Here, the moduli distribution is given by the distribution of 2L 3 +iM in (4.7), and L 3 and M are indeed mutually commuting variables.
The commutativity of the low energy moduli matrices might be general phenomena which occur in the strong coupling limit. As observed in [20], in some matrix models with commutator interactions, commuting matrices indeed arise in the strong coupling limit. A possible mechanism is as follows. For Yang-Mills type matrix models, one can rescale the matrices in such a way that the coupling constant appears in front of each commutators.
In the strong coupling limit, in order to have a finite value of the action, the values of commutators themselves must become small unless there is some cancellation with the kinetic terms. If this occurs, the matrices become commuting with each other. Though observing this phenomena directly in the current model is very difficult, this is very likely to occur in the low energy region, since in the low energy limit, the kinetic terms of the matrices are very small and there will be no chance to have a cancellation between the kinetic terms and the commutator terms.
Thus, we identify the real and imaginary parts of φ in the formula (4.2) with the low energy moduli for Y 3 and Y 9 , respectively. In particular, ρ in (4.11) is identified with the moduli of Y 9 . Recall that, in the decoupling limit of M5-brane, we have seen that the spectrum of φ becomes pure imaginary. Hence, with the suitable normalization of matrices (namely, going back to the original normalization in (3.8)), one finds that the moduli of the SO(6) scalar have a wide distribution while the moduli of the SO(3) scalars collapse to the origin in the decoupling limit of the M5-brane. Now, let us consider the description of the spherical M5-brane. We consider the SO (6) symmetric uplift of the distribution [21,22] of the moduli of a single SO(6) scalar. The uplifted distributionρ is defined as the solution of for any n, where r = x 2 a is the distance from the origin. The normalization factor on the right-hand side is chosen so thatρ represents a density function before the rescaling (3.8). For the density ρ in (4.11), the unique solution to (4.13) is (4.14) The radius r 0 is given by (4.15) For N 5 = 1, the shape of the density function of the SO(6) moduli agrees with the shape of the spherical M5-brane. In particular, the radius shows a perfect agreement with the M5-brane: r 0 = r M5 . Therefore, we conclude that the spherical M5-brane is indeed realized as the low energy moduli distribution of the SO(6) scalar fields in PWMM.
The N 5 -dependence of the radius agrees with the expected form in [3] based on the perturbative expansion in PWMM.

Multiple M5-branes from generic partitions
Let us generalize the above calculation to the case of the general partition (1.2). According to [3], this corresponds to Λ stacks of M5-branes with different radii as shown in Fig. 1.
As we discussed in section 2.3, to make the M5-branes decouple from the bulk gravity, we consider the limit (2.35) such that the large-N limit is taken as in (1.4).

(4.22)
Note that the coefficient (n s − n s+1 ) ofρ s is just the number of M5-branes in the sth stack. From these properties,ρ s can be naturally identified with the density function for an M5-brane in the sth stack.
Obviously, the SO(6) symmetric uplift of {ρ s } is given by Λ stacks of the spherical shells. By taking the rescaling (3.8) into account, the sth stack has the radius where p + s is defined in (1.3). Thus, we have shown that, as shown in Fig. 1, the generic partition indeed describes concentric stacks of M5-branes with radii given by (4.23).
branes in the type IIA superstring theory can be understood even at the level of the classical action. However, it it still nontrivial whether we can observe the emergence in the strong coupling region of PWMM. Here, we study the emergence of M2-branes in the decoupling limit of the M2-branes. In this section, we only consider the simplest partition (4.6) but the generalization is straightforward.
In the M2-brane limit, where N 5 goes to infinity, the one-loop determinant of the eigenvalue integral (4.7) converges to the hyperbolic tangent function. Thus, we obtain Note that the model depends only on N 2 and g 2 /N 5 .
The typical value of the eigenvalues of this model should depend on N 2 and g 2 /N 5 .
Then, in the decoupling limit of M2-branes, the typical value is much smaller than N 5 .
This implies that, in the result of the localization (4.2), the eigenvalue distribution of M is much narrower than that of L 3 . Hence the spectrum of φ lies on the real axis in this limit.
This implies that the moduli of Y 3 are given by the classical vacuum configuration 2L 3 while the moduli of Y 9 collapse to the origin. It is easy to see that the SO(3) symmetric uplift of this configuration gives the two-sphere and the radius agrees with that of the spherical M2-brane on the supergravity side for N 2 = 1. Thus, we see that the spherical M2-brane is also realized as the moduli of SO (3) scalars.
However, we should notice that, unlike the decoupling limit of M5-branes, the instanton corrections could contribute to the partition function in the M2-brane limit 13 . If this is the case, since the result of the localization does not include the instanton corrections [11][12][13], our computation is not correct. Then, in order for our computation to make sense, we need to consider the limit where the number of M2-branes goes to infinity.
In this limit, the instanton effects will be suppressed. Thus, at least in the large-N 2 case, the result of the localization shows the emergence of the spherical M2-branes in PWMM.
When N 2 is large, we can find an exact solution for the eigenvalue distribution of (5.1) and can check that the typical value of the eigenvalues is proportional to (λ/N 5 ) 1/3 for large λ/N 5 . See appendix D.
Of course, there is still a possibility that the instantons do not affect our computation.
For instance, this happens if there exists a fermionic zero mode at the saddle points of instantons in the localization computation. This needs a further analysis of the localization saddles in PWMM.

Summary and discussion
In this paper, we tested a conjecture on the description of spherical M5-branes in the matrix model formulation of M-theory. We considered the plane wave matrix model (PWMM), which is expected to describe the M-theory on the maximally supersymmetric 11-dimensional plane wave geometry.
We first reviewed that, in the M-theory, there exist spherical M2-and M5-branes with zero light cone energy. These spherical branes are considered to be described as certain vacuum states in PWMM. This relation between the spherical branes and the vacua of PWMM is stated in [3]. In particular, it is conjectured that a single spherical M5-brane corresponds to the trivial vacuum of PWMM.
Through a direct computation in PWMM using the localization, we showed that the spherical M2-and M5-branes are formed by the distribution of the moduli of SO(3) and SO(6) scalar fields, respectively. This result strongly supports the proposal in [3].
As we discussed in section 4.2, we can assume that the moduli in PWMM are given by commuting matrices in the decoupling limit of the M5-branes. Here, let us consider a possible effective theory of these commuting matrices in the decoupling limit. We require the theory to have the SO(6) symmetry and to be able to reproduce our result of the localization. For the case of coincident M5-branes, a possible solution to these requirements is given by a commuting matrix model with 6 matrices defined by 14 whereρ is the distribution of moduli y a i (a = 4, 5, · · · , 9, i = 1, 2, · · · , N ) for the SO (6) 14 The same model was also considered in different contexts [22][23][24].
scalars Y a ,ρ The Lagrange multiplier β is introduced to impose the normalization condition onρ. The second term in (6.1) is understood as the Vandermonde determinant i<j | y i − y j | 2 for the commuting matrices. We fix the parameter m in (6.1) as so that the model reproduces the result of the localization below. In the 't Hooft limit, the WKB approximation becomes exact. The saddle point equation is given by The solution to this equation is obtained in [22][23][24] aŝ ). It would be interesting to find more general commuting matrix model, which reproduces our result for the general partition. In addition, we also need to investigate whether some low energy excitations can also be reproduced from the commuting matrix model or not.
Finding a good description of the low energy theory should be one of the most important problem in understanding the description of the classical geometry in the matrix theory. We hope that our result gives a clue to this problem.

A Localization in PWMM
In this appendix, we perform the localization and derive the formula (4.2). In this appendix, following the method in [10], we use a Lorentzian signature obtained by a double Wick rotation for the time-direction and the direction of one of the SO(6) scalar fields.
To use some 10 dimensional notation, we relabel the SO(3) scalar fields in PWMM as , the SO(6) scalar fields as (Y 4 , · · · , Y 9 ) → (Y 5 , · · · , Y 10 ) and the gauge field as A → Y 1 . The double Wick rotation is performed for the Y 1 and Y 10 directions and hence, the Y 1 's direction is Euclidean and Y 10 's direction is Lorentzian. We also use Y 0 to express the scalar field in the Lorentzian signature, which is related to Y 10 by Y 0 = iY 10 .
A.1 Off-shell supersymmetry of PWMM In the above notation, the full action of PWMM can be written in the 10-dimensional notation as Here, Ψ is the 10-dimensional Majorana Weyl spinor with 16 components and we use the gamma matrices defined in [11]. We have also used the following notation: where i, j, k = 2, 3, 4 and a, b = 5, 6, · · · , 10. In order to realize the off-shell supersymmetries, we further add seven auxiliary fields The theory has the off-shell supersymmetry, See [11] for the definition ofΓ a . The parameter has to satisfy the Killing spinor equation of PWMM and the closure of the supersymmetry requires ν I to satisfy The following spinors give a solution to these conditions: where η 1 is any 4-component constant vector. We use η 1 = (1, 0, 0, 0) in the following computation.

A.2 Saddle point of the localization
To perform the localization, we add an exact term tδ s V to the action, where After some calculation, one can find that the bosonic part of δ s V is calculated to be where c := cosh τ and S is defined by The derivatives D M are defined in (A.2). F ± ab stands for the selfdual and anti-selfdual part of F ± ab in the subspace a, b = 1, 2, 3, 4 or a, b = 5, 6, 7, 8. After the Wick rotation, Y 0 = iY 10 and K i = iK (E) i (i = 1, 2, · · · , 7), the bosonic part δ s V | bos becomes a sum of positive-definite terms.
We consider the theory around a fixed vacuum (1.2). Then, we impose the boundary condition such that all fields approach to the vacuum configuration at τ → ±∞. Then, taking the temporal gauge Y 1 = 0, we find that the saddle point configuration is given bŷ
We also introduce the ghost fields, (C, C 0 ,C,C 0 , b, b 0 , a 0 ,ã 0 ), where (b, b 0 , a 0 ,ã 0 ) are bosonic and (C,C, C 0 ,C 0 ) are fermionic fields. The fields with subscript 0 shall contain only zero modes for both τ direction and the fuzzy sphere directions. They are defined through the following BRS transformations, The commutator in the above equation shall express the anti-commutator for fermionic variables. The square of δ B is a gauge transformation with parameter a 0 , Then Q = δ s + δ B has the following action: One can easily show that Q 2 is given as The gauge-fixing and ghost actions are defined by where F corresponds to the gauge fixing condition. We use

A.4 1-loop determinants
Let us perform the 1-loop calculation around the saddle point (A.10). We first redefine the fields asX and divide the fields to four groups as They form doublets under the action of Q as where R := Q 2 is given by the sum of the U (1) and gauge transformations as shown in (A.20).
Then we expand the full action S P W + tQ(V + V gh ) around the saddle point configuration (A.10) as Z i →Ẑ i + Z i and Z i →Ẑ i + Z i . Then the quadratic part of the fluctuations in V + V gh is schematically written as where D ij (i, j = 0, 1) are some linear differential operators. Thus, the quadratic part of the action takes the form Hence, the one-loop integral produces the determinants, Here, the determinants should be taken in the appropriate functional spaces of the fluctuations. Recall that we adopted the boundary condition such that all fields go to the vacuum configuration as τ → ∞. This implies that the fluctuations should vanish at infinities.
Note that D 10 is a linear map from V Z 0 to V Z 1 and commutes with R. Then the determinants in (A.29) cancel between ImD 10 ⊂ V Z 1 and ImD * 10 ⊂ V Z 0 , where D * 10 is the adjoint of D 10 . Hence, the 1-loop determinant reduces to Furthermore, since R and D 10 commute, the kernel and the cokernel are given by direct sums of the eigenspaces of R. Thus, we can express the 1-loop determinant as where V r i and V r i are the restrictions of the kernel and the cokernel to the eigenspace of R with eigenvalue r i , respectively. Therefore, the remaining task is to evaluate r i and the index dimV r i − dimV r i in each eigenspace.
By integrating the ghost fieldã 0 , we obtain the constraint a 0 = −2φ. At the saddle point, this is equal to −2iM + 4L 4 . Thus, r i is given by the sum of eigenvalues of [−2iM + 4L 4 , ] and the diagonal U (1) charge.
By studying the structure of D 10 for each supersymmetry multiplet, we can easily compute the index. The result is as follows [11]. The contribution from the hypermultiplet, which contains Y 5 , · · · , Y 8 , is given by Λ s,t=1 (A.32) The contribution from the vector multiplet, which contains Y 1 , · · · , Y 4 , Y 9 , is given by Combining these contributions with the Vandermonde determinant for diagonalizing M , we obtain the 1-loop determinant (4.5). See below for the derivation of these 1-loop determinants.

A.5 Derivation of 1-loop determinants
The relevant part of the action is given by Z 1 D 10 Z 0 . In terms of the component fields, this can be written explicitly as where P denotes the path ordering. The general solution to the differential equation Df = 0 is then given by where

Hypermultiplet
Let us consider the hypermultiplet. We use complex combinations, We can read off the action of D 10 from (A.34). If (W 1 , W 2 ) is an element of kerD 10 , we have where s = sinh τ and c = cosh τ . To analyze the structure of these equation, we use the fuzzy spherical harmonics, which behave nicely under the adjoint action of L i . See [25][26][27] for the definition. For the vacuum of the form (  2J−J and their complex conjugates contribute to the index. The eigenvalues of R for these modes are r = 2(±(2J + 1) + i(q si − q tj )), and thus we obtain (A.32).

Vector multiplet
Next, we consider the vector multiplet. We first calculate dim(kerD 10 ). For {Y A ,ã 0 , b 0 |A = 1, 2, 3, 4, 9} ∈ kerD 10 , we have To simplify the equations, let us consider the limit τ → ±∞ in (A.48). Since F → 0 in this limit, we obtain b 0 = 0. Noticing that b 0 has only the constant mode, by using (A.48) again, we find that F should be vanishing for arbitrary point on R, namely, for J = 0. Therefore, the relation Y 4 = sY 9 holds for any J.
We make some redefinitions asC and also introduce complex fields, Υ ± = Υ 6 ± iΥ 7 . With this notation, we can write (A.58) as

B The saddle point equation
In this appendix, we derive (4.10). We start with the effective action of (4.7) , (B.1) Here, λ = g 2 N 2 is the 't Hooft coupling and β is the Lagrange multiplier for the normalization of the eigenvalue density (4.9). The saddle point equation is obtained by differentiating S ef f with respect to ρ(q) and is given by Here the integral of q should be understood as the principal value.
We first consider the following identity, log tanh 2 πx 2 = log πx 2 which follows from the infinite product expression of the hyperbolic sine and cosine functions. By using this identity, we find that the second term in (B.2) can be written as up to a constant term. We ignore the constant term since it can always be absorbed by a redefinition of β. In the regime where N 5 is finite but λ is very large, q m also becomes very large. To see the q m -dependence clearly, let us rescale the variables as q = q m ξ. From the fact that q m log tanh 2 πq m ξ 2 → −πδ(ξ) (q m → ∞), (B.5) we find that the first term in (B.4) is equal to πρ(q) in this limit. In the second term in (B.4), we approximate the discrete sum with a continuous integral by replacing J/q m → η and ∞ J=N 5 → q m ∞ N 5 /qm dη. Then, the second term can be evaluated as Thus, in the strongly coupled regime, the saddle point equation (B.2) is reduced to (4.10).

C Solving the saddle point equation
In this appendix, we construct solutions of the saddle point equations of the eigenvalue integrals obtained by the localization.
The equation (C.6) is identical with the equation of motion of the quartic one matrix model. Hence, the solution takes the same form as the quartic matrix model, where the resolvent is written as We substitute (C.6) into this expression. Then, the asymptotic behavior of the resolvent, ω(z) → 1 z (z → ∞), gives three conditions, which enable us to express a, b and β in terms of q m : Thus, the resolvent is finally determined as The eigenvalue density is given by the discontinuity of (C.9) as ρ(q) = 1 π In the region (C.11), the minimum of S ef f is realized at q 4 m = 8λ. (C.14) By substitute this into (C.10), we obtain (4.11).

C.2 For the generic partition
Here, we construct a solution to (4.18) in the decoupling limit of M5-branes.
(C. 25) Again, this is the same form as (C.5), so that the solution forρ r is given by the same form as (C.10).
Finally we determine q r . By using the equation of motion for ρ s , the on-shell can be computed as where λ s = g 2 s r=1 N (r) 2 . Thus, the minimum is given by q s = (8λ s ) 1/4 . Thus, we obtained (4.19).

D Eigenvalue distribution in the M2-brane limit
In this appendix, we solve the eigenvalue integral (5.1) for large N 2 . Putting a = g 2 N 2 /N 5 , we consider the scaling limit such that N 5 , N 2 , a → ∞, N 5 /N 2 → ∞ and N 2 /a → ∞.
We again assume that the typical value of eigenvalues is very large in this limit, since the Gaussian attractive force becomes weak.
We introduce the eigenvalue density ρ(q) as (4.8). The effective action of (5.1) is written in terms of ρ(q) as . Note that this result is consistent with our assumption that the typical value of the eigenvalues is very large in the strong coupling region.