ORIGINAL RESEARCH article
Sec. High-Energy and Astroparticle Physics
A Matrix Model from String Field Theory
- Yokote Seiryo Gakuin High School, Yokote, Japan
We demonstrate that a Hermitian matrix model can be derived from level truncated open string field theory with Chan-Paton factors. The Hermitian matrix is coupled with a scalar and U(N) vectors which are responsible for the D-brane at the tachyon vacuum. Effective potential for the scalar is evaluated both for finite and large N. Increase of potential height is observed in both cases. The large N matrix integral is identified with a system of N ZZ branes and a ghost FZZT brane.
Two different candidates for the nonperturbative formulation of string theory have been known—string field theory (SFT) and matrix model. It is commonly believed that they are just different descriptions of the underlying theory. Therefore, it is important to investigate the relationship between these formulations. However, a few examples have been known to reduce open SFT (OSFT)  to certain matrix models. First is OSFT for topological A or B model  which reduces to Chern-Simons matrix model  or ordinary Hermitian matrix model  respectively. Second is OSFT for (2,1) minimal string  which reduces to Kontsevich matrix model . In addition, less direct examples for c = 1 , c = 0 , and critical  strings have been known. In those examples, each matrix model is obtained from different OSFT associated with particular boundary conformal field theory (BCFT). A systematic way to derive different matrix models from OSFT in fixed background has not yet known. Finding such method is important to study the background independence of OSFT.
In this paper, we present an example for such method. Our idea is simple: instead of varying BCFT, we start with a general setup in critical string, and approximate the string field to its first few components. An example we will study is the level truncation in the universal sector  of critical (D = 26) OSFT in which the string field is approximated at level n as
A possible reason that such method has not yet been examined is the lack of understanding of Chan-Paton factors in OSFT. Although Chan-Paton factors can be introduced to OSFT consistently, their origin has not yet been explained. Recently, Erler and Maccaferri proposed a new construction of classical solutions of OSFT in terms of the regularized boundary conditions changing (BCC) projectors . Their construction covers a wide range of backgrounds including multiple D-branes, which is the main interest of present paper. Following their work, Kishimoto, Masuda, Takahashi, and Takemoto (KMTT) demonstrated that N Chan-Paton factors naturally arise from the decomposition of string field in terms of the regularized BCC projectors . SFT expanded around the multiple D-branes solution was interpreted as a system of N + 1 D-branes. We will employ their formulation as our foundation.
This paper is organized as follows. Section 2 introduces KMTT formulation. In Section 3, we derive a one-matrix model in terms of the level truncation. In Section 4, we derive the partition function of the model both in finite and large N. Section 5 summarizes our results and further discussions are given.
2. KMTT Formulation
where ΨT is the tachyon vacuum solution and Σa and are the regularized BCC projectors which obey
Here QT is the kinetic operator at the tachyon vacuum defined by
The index a in Σa or is a Chan-Paton factor which labels D-branes. The solution (Equation 3) carries degrees of freedom required for multiple D-branes. We further expect that non-Abelian gauge gauge symmetry is realized in OSFT expanded around (Equation 3). KMTT  introduced the decomposition of the string field which realizes this:
where χ, χa, and χab are component string fields. By expanding OSFT action around Ψ0 according to Equation (6), they derived the matrix action
where QT0 is the kinetic operator defined by QT0Ψ = QBΨ + ΨTΨ. The trace in (Equation 8) runs over indices of ϕab, and identical indices in Equation (10) are to be summed over3. KMTT  claimed that (Equation 7) describes N + 1 D-branes rather than N D-branes. In addition to the N unstable D-branes, there is a “D-brane” at the tachyon vacuum described by χ. According to their interpretation, ϕab connects two unstable D-branes a and b while χa connects D-brane a with the D-brane at the tachyon vacuum. χ represents fluctuation on the D-brane at the tachyon vacuum.
A remarkable feature of the action (Equation 7) is the presence of the “vector” sector described by χa and , which has not been found in literature. The action S3 is quadratic for the vectors therefore can be integrated out in the path integral. With assuming suitable gauge fixing, we perform path integral for χa and and obtain a determinant factor
in the partition function. In the following sections, we will evaluate this determinant by truncating the string field rather than imposing conventional gauge condition such as linear gauge .
3. Truncation to Matrices
Although the string fields presented in previous section carry Chan-Paton factors, they also depend on infinitely many labels which distinguish each state in a BCFT. In order to make our analysis tractable, we introduce an approximation in which dynamical variables reduce to matrices. Let us consider the string field in the universal sector4 truncated at level n, where the level is defined by eigenvalue of the kinetic operator of OSFT. Denoting the level k base element of the string field ψk, we truncate each component in Equation (6) as
where is a Hermite matrix, is a complex vector, is its complex conjugate, and t(k) is a real number. Note that these component fields do not depend on any other variables. As a concrete example of such truncation, we choose the expansion examined in Zeze  for dressed 0 gauge . In this case, the level k base is given by
where K, B, c are the elements of the KBc algebra . We would like to study level 1 truncation in which string fields are given by
A matrix action obtained by this truncation is rather complex since it includes both level 0 and 1 fields. However, it is possible to reduce degrees of freedom further. It is know that level 1 fields cannot have cubic terms since 
Thus, level 1 fields are quadratic in the action and can be integrated out using equations of motion. With the help of the explicit values of products of the basis ψ0 and ψ1 evaluated in Zeze , a nontrivial solution of equations of motion for , and t(1) turns out to be
By this assignment, four string fields in Equation (12) become proportional to Erler-Schnabl solution :
Thus, the four fields5 can be written as
Although only few fields are included in the approximation, we expect that it well captures the essence of the dynamics of D-branes, as is the case in the conventional level truncation analysis. An evidence that supports our expectation is that the truncation interpolates between two analytic solutions: one is the Erler-Maccaferri's multiple D-branes denoted as while the other is the perturbative vacuum, also denoted as (δab, 0, 0, −1).
Let us derive the truncated action. The requirement that ΨT reproduces the correct value of the D-brane tension is represented by
Applying these to Equation (7), we find
where we have rescaled the open string coupling as
For later convenience, we shift M to M + 1 and omit the prime in g′. Then, ignoring constant shift, we obtain an action
where . In this way, the truncation (Equation 12) reduces OSFT action to a cubic matrix action coupled with a scalar and complex vectors6. The partition function for Equation (25) can be obtained by employing the standard technique of matrix model7. As readily found in Equation (25), the action is invariant under U(N) transformation , where U denotes an U(N) matrix. This is a residual gauge symmetry of KMTT action . This U(N) symmetry can be fixed by diagonalizing Mab to its eigenvalues λa with inserting Vandermonde determinant in the partition function. After performing Gaussian integral for ξ, we obtain a partition function
This can describes a system of N + 1 particles moving in the potential W. An eigenvalue λa feels repulsive forces from other eigenvalues through the last term of Equation (27). Also, it is attracted toward −t − 1 due to the second term in Equation (27).
It is interesting to compare our result with the matrix formulation of c < 1 noncritical string theories. It has been recognized that Hermitian one-matrix model serves nonperturbative definition of (2, 2k + 1) minimal string theory. As an example, let us consider the action studied in Kutasov et al. :
where ψ and are fermionic vectors. They can be integrated out so that insert a factor
in the matrix integral. In double scaling limit, this determinant is identified with a FZZT brane [24, 25], and the fermionic vectors are identified with fermionic strings between a FZZT brane and a stuck of N ZZ branes . In contrast, our action (Equation 25) contains bosonic vectors ξ and rather than fermionic ones. Integration with respect to them yields a factor
in the partition function. Comparison between Equation (29) and Equation (30) naturally identifies the determinant Equation (30) as a ghost FZZT brane  which cancels the effect of a FZZT brane8. Unfortunately, corresponding observable in minimal string theory is not yet identified.
4. Effective Potential
In this section, we will evaluate the effective potential for t in terms of saddle point equations for eigenvalues. In large N limit of t'Hooft expansion, saddle point configurations are leading contributions to the matrix integral. Even for finite N, saddle point configurations also offer a good approximation to the matrix integral when g is small. In either case, saddle point equations are obtained from the variation of Equation (27):
4.1. N = 1
Let us study the dynamics of the model at finite N. We begin with N = 1, which corresponds to a system of an unstable D-brane and another D-brane at the tachyon vacuum. Denoting λa as λ, the potential (Equation 27) is given by
and saddle point equations are
By combining above two equations, we obtain a coupling independent equation
which has two roots λ = t and λ = −t − 1. The latter is not appropriate since it hits the singularity (t + λ + 1)−1 in the partition function. Therefore, we choose t = λ as our solution. Given this choice, Equation (34) and Equation (35) reduce to single equation
This can be rewritten to a cubic equation,
whose roots correspond to saddle points. It is easily understood that while there are three roots for small g, two of them disappear beyond critical value of g. Let us discuss further details as follows. Small g expansion of these roots is
It is useful to show positions of these roots in a plot of the potential W. Figures 1, 2 are such plots for different values of g. Figure 1 shows that saddle points for small g; saddle points are placed within the region (−1, 0) with equal intervals.
As g become larger, two larger roots become closer and annihilate beyond the critical value
Figure 2 shows a placement of a root after the annihilation.
Analysis made here is merely a classical approximation and not full quantum treatment. Fortunately, we can integrate out one variable in the potential without relying on the saddle point approximation since present example is enough simple. Let us change the variables as
Then, the partition function can be written as
where we performed Gaussian integral for v in the second line. Thus, we obtain one dimensional effective potential
We note that the saddle point approximation corresponds to setting v = 0 in Equation (42), since we our choice of the saddle point, λ = t corresponds to it. Therefore,
By comparing Equation (44) and Equation (45), we find that they only differ in the coefficient of the logarithmic term which is negligible for small g. Thus, the saddle point approximation works well for small coupling.
Finally, let us compare the shape of the potential for different values of g. Figure 3 is a plot of Equation (44) for small and large g. It is observed that the stable vacuum around u = 0 disappears for large g due to dominance of the logarithmic term. This phenomenon corresponds to the annihilation of saddle points which has already been observed in the saddle point analysis. Such dependence of the effective potential on the coupling is consistent with the fact that the D-brane describes the system at small coupling.
Figure 3. Effective potential for u, where solid line is a plot for g = 0.1 while dashed line is for g = 1.0.
4.2. N ≥ 2
Next we proceed to N ≥ 2, where the equations of motions are given by Equations (31) and (32). We again perform saddle point approximation by finding roots of these equations. The roots can be obtained numerically for each value of g. Numerical solutions for g = 0.01 and up to N = 3 are shown in Tables 1–3. We pick real roots only, and each solution is specified as (t, λ1, …, λN). Due to the symmetry which exchanges λs, we only need to specify single configuration for each value of t among N permutations of λ. We also evaluate the value of effective potential V(λ, t) normalized by the value of D-brane tension 1/(6g2).
We observe an interesting pattern in these results. Each root is located at one of a ‘site’ which is obtained by dividing the region −1 ≤ λ ≤ 0 into N + 1 intervals, i.e.,
We also observe that the minimum value of t is given by −2/(N + 1) while values of λs can reach −1. We also observed that the sum of all eigenvalues are close to its value of potential, i.e.,
Finally, it is also interesting to see that the maximum value of the potential is close to 2 for all N.
We would like to close this subsection with a summary of our result:
• The stable vacuum is lost for large coupling. This indicates breakdown of D-brane description. The critical value of the coupling can be determined by saddle point equations.
• At small g, the local maximum of the effective potential is about twice higher than that of original one.
4.3. Large N
Let us evaluate the partition function in large N limit following with the standard method of matrix model (hep-th/0410165). We begin with the partition function
and introduce the t'Hooft coupling
Then, saddle point equations read
In large N limit, eigenvalues are described by a continuous distribution ρ(λ) and summation for eigenvalues is replaced with integration
Then, saddle point equations are replaced with
At leading order, the second term in Equation (54) becomes negligible. Thus, we obtain saddle point equations
The latter equation (Equation 56) solves the planar limit of a cubic matrix model whose solution can be found elsewhere . It is convenient to introduce the resolvent
The equation (Equation 56) can be replaced with an equation for ω(z). Once ω(z) is obtained, Equation (55) can be solved by finding solutions of
The “one-cut” solution for the resolvent in large N limit is known to be 
where a and b are endpoints of the branch cut, which define a support for the eigenvalue density ρ(λ). Requiring ω(z) ~ z−1 at infinity, we obtain equations
which solve a and b as functions of μ. These equations are conveniently rewritten in terms of parameters σ = a + b and as
The latter equation restricts σ inside −2 ≤ σ ≤ 0. Further, Equation (62) tells us that there are no real roots of Equation (62) within −2 < σ < 1. Therefore, σ is constrained within
Let us choose a brunch which starts from σ = 0. From Equation (62), the maximum value of μ reads
Given these ingredients, Equation (55) can be rewritten into
Integrating this equation yields an effective potential for t,
A plot of effective potential (Equation 67) is shown in Figure 4. The local minimum at t = 0 and the local maximum around t = −1 are observed for any value of μ. This is consistent with the t'Hooft limit, in which g is small so that the D-brane description of the system holds. On the other hand, the height of the local maximum is always higher than that of the original potential W. The height is maximum for μ = 0; in this case, the potential is twice higher than W. While the height decreases as μ increases, it remains higher than W even at a maximum value of μ = μc. The increase of the height from the original one can be understood from the particle description of eigenvalues. Eigenvalues filled in the bottom of W(t) pull the tachyon t by attractive force. It makes harder for t to climb the potential wall, thus increases the height of the potential height. Such dependence of the effective potential on μ is quite different from that of the probe eigenvalue model investigated in Hanada et al.  where a plateau along the eigenvalue distribution is observed. The existence of the plateau is explained by the fact that a probe eigenvalue cannot be distinguished from others. Our model has no plateau since t can be distinguished from other eigenvalues.
Figure 4. Effective potential at large N. The value of the potential is normalized by 1/(6μ). The original potential W(t) is shown together.
5. Conclusions and Discussions
In this paper, we proposed a systematic method to derive matrix models from level truncated OSFT. Obtained matrix model contains U(N) vectors and a scalar in addition to Hermite matrix. We have evaluated the effective potential of the scalar both for finite and large N. Increase of the potential height was observed at small coupling. In Section 3, we have interpreted our model as a system of ZZ branes and a ghost FZZT brane.
We would like to discuss further issues to be explored. First, we would like to present an alternative but rather heuristic interpretation of our result. Let us go back to the inverse of the determinant :
The basic idea is that this quantity can be regarded as a propagator with ϕ and χ insertions. Recall that this factor is obtained by an integrating out χa and which connect a D-brane with the tachyon vacuum where the world-sheet boundary disappears. Therefore, it is natural to think that Equation (70) amounts to a disk amplitude with single world-sheet boundary of a D-brane. Schematically, such disk amplitude can be written
As is well known, small t limit of such amplitude corresponds to closed string propagation . Therefore, the inverse determinant (Equation 70) encodes gravitational force between N D-branes and the tachyon vacuum9.
Second issue is abut the level truncation. We have seen that the approximated OSFT action at first few levels yields a one-matrix model which can be interpreted as c < 1 noncritical string theory . It is interesting to improve the approximation by including higher level fields to obtain multi-matrix models. We speculate that that the improved matrix action continues to be dual to some closed string theory on more nontrivial background. Finally, we will recover original OSFT with infinitely many matrices in infinite level limit. We also expect that this matrix model describes a critical string theory in nontrivial background through AdS/CFT like duality . Thus, the level truncated OSFT offers a way to describe closed strings in approximated geometry.
Last issue is the matrix description of OSFT based on the left-right splitting of open strings which have been examined in past [33–36]. Our model based on the KMTT decomposition looks quite differently from these models. However, as mentioned in the previous paragraph, our model will recover full OSFT in infinite level. Thus, the left-right splitting models and our model both describe same OSFT. We expect that all ingredients of KMTT, including Chan-Paton factors and regularized BCC projectors, are embedded in the left-right type matrix model in quite nontrivial manner.
Together with recent development which deals with different backgrounds as classical solutions , our result presents further evidence for SFT as a formulation of nonperturbative string theory. We hope that further developments in this direction will shed light on the landscape of string theory.
The author confirms being the sole contributor of this work and approved it for publication.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. ^Each component field carries Chan-Paton indices.
2. ^As we will see later, extra vectors and scalars couple with the matrix. Here we omit them for simplicity.
3. ^We will ignore a constant shift in the action hereafter since it is not relevant for remaining discussions.
5. ^This truncation makes sense if we replace ΨT with another representation of the tachyon vacuum solution, since we do not require its explicit form in following analysis.
7. ^For example, see hep-th/0410165.
9. ^This interpretation is consistent with the attractive force between λ and t observed in Section 3.
29. Hanada M, Hayakawa M, Ishibashi N, Kawai H, Kuroki T, Matsuo Y, et al. Loops versus matrices: the nonperturbative aspects of noncritical string. Prog Theor Phys. (2004) 112:131–81. doi: 10.1143/PTP.112.131
Keywords: string field theory, D-branes, matrix model, string theory, tachyon Condensation
Citation: Zeze S (2016) A Matrix Model from String Field Theory. Front. Phys. 4:39. doi: 10.3389/fphy.2016.00039
Received: 20 March 2016; Accepted: 22 August 2016;
Published: 07 September 2016.
Edited by:Osvaldo Civitarese, National University of La Plata, Argentina
Reviewed by:Pietro Antonio Grassi, University of Eastern Piedmont, Italy
Nicolás Esteban Grandi, National Scientific and Technical Research Council, Argentina
Copyright © 2016 Zeze. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Syoji Zeze, email@example.com