Monte Carlo simulations of a supersymmetric matrix model of dynamical compactification in non perturbative string theory
The IKKT or IIB matrix model has been postulated to be a non perturbative definition of superstring theory. It has the attractive feature that spacetime is dynamically generated, which makes possible the scenario of dynamical compactification of extra dimensions, which in the Euclidean model manifests by spontaneously breaking the SO(10) rotational invariance (SSB). In this work we study using Monte Carlo simulations the 6 dimensional version of the Euclidean IIB matrix model. Simulations are found to be plagued by a strong complex action problem and the factorization method is used for effective sampling and computing expectation values of the extent of spacetime in various dimensions. Our results are consistent with calculations using the Gaussian Expansion method which predict SSB to SO(3) symmetric vacua, a finite universal extent of the compactified dimensions and finite spacetime volume.
SUSY matrix model of dynamical compactification
\FullConferenceThe 30th International Symposium on Lattice Field Theory
June 24 â 29, 2012
Superstring theory is a natural candidate of a unified theory of all interactions, including gravity. It lacks a non-perturbative definition, which would allow us to address dynamically questions such as the preferred vacuum where our universe sits and determine its properties such as classical spacetime dimensionality. The large limit of the 10 dimensional IKKT or IIB matrix model has been proposed to provide such a definition . Using dualities, which suggest that the five types of superstring theory are actually limits of a unique theory, the IIB matrix model is expected to describe the unique underlying theory despite its explicit connection to perturbative type IIB superstring theory.
The IIB matrix model has a series of attractive features: The model
has a unique scale, which raises the possibility for the theory to
choose a unique vacuum. Spacetime and matter content can arise
dynamically from the distribution of the eigenvalues of the bosonic
matrices, which makes possible the realization of the scenario of
dynamical compactification of extra dimensions. The motivation for
studying such a scenario in the Euclidean model comes from noticing
that lower dimensional configurations are stationary points of the
fluctuating complex part of the action of the model
. Using the
Gaussian Expansion Method (GEM), one can show that the SO(3) symmetric
vacuum has the smallest free energy density compared to higher or
lower dimensional configurations . Dynamical
therefore, occurs by spontaneously breaking the SO(10) rotational
symmetry of the model (SSB). Moreover, a universal scale for the
small dimensions is computed
Monte Carlo simulations of the Euclidean IIB matrix model, as well as simpler related models, can confirm these results from first principles and shine light into the mechanism of dynamical compactification of extra dimensions [4, 5]. Such simulations are hindered by a severe complex action problem and are quite difficult. Some results, however, have been obtained by studying related lower dimensional toy models using the so called factorization method. This method, proposed originally in  has been tested also in Random Matrix theory and finite density QCD . It attempts to sample the most important configurations that contribute to the partition function, first by numerically solving equations that compute the maxima of the distribution functions of appropriately chosen observables that are strongly correlated to the fluctuating phase  and then by simulating a constrained system in the region of the solutions. This makes importance sampling possible by also taking into account the suppression of configurations by the fluctuations of the phase, together with the suppression caused by the measure and the real part of the action in the partition function. The overlap problem is thus solved and the complex action problem reduced due to milder fluctuations of the phase within the sampled configurations and the use of scaling properties that permit extrapolations to larger systems.
In this work, we present Monte Carlo simulation results of a 6 dimensional version of the IIB Matrix Model. Contrary to , where the oneloop approximation was used, the full model is simulated. This is necessary since, as we will show, the short distance non perturbative dynamics of the eigenvalues of the matrices play a crucial role in generating the scales of dynamical compactification. We show that in the absence of the complex part of the action, no SSB of the SO(6) rotational invariance occurs. We apply the factorization method in order to compute the eigenvalues of spacetime’s “moment of inertia” tensor and compare our results with those obtained in  using the GEM. In particular, we test the GEM predictions of SSB to SO(3) symmetric vacua, a universal scale for compactified dimensions and finite spacetime volume. Our results are consistent with these scenaria.
2 The Model
We study a 6 dimensional version of the IIB matrix model defined by the partition function
The model can be formally obtained by the dimensional reduction of the
Euclidean 6 dimensional U() Super Yang Mills theory
to zero dimensions: The matrices ,
are traceless and hermitian and transform as O(6) vectors. The
matrices , are traceless
with grassmannian entries and transform as Weyl spinors. The model
turns out to have supersymmetry, which leads to the
interpretation of the eigenvalues of defining the points of
the 6 dimensional (Euclidean) spacetime
where is the square of the extent of the SO(6) symmetric configurations, and . The SO(6) symmetric configurations dominate in the phase quenched model discussed below where SSB does not occur . Eq. (2) expresses the constant volume property of the SO() vacua and implies that the volume of spacetime in the 6 dimensional Euclidean IIB matrix model is finite .
Monte Carlo simulations are performed by integrating out the fermions in (1) first, obtaining
where and . The determinant turns out to be generically complex creating a very strong complex action problem in the simulations. We applied the factorization method , which first amounts to considering the phase quenched model . In this model, we computed the phase quenched expectation values , and defined the normalized eigenvalues . Deviation of from 1, is the result of the suppression of the dominant configurations of by the fluctuations of the phase . Due to the strong correlations of with , we consider the distribution functions [6, 8]
The value of in the large limit is determined by the minimum of the “free energy” which, given Eq. (4), is the large limit of a solution to
Using the dominant solution of (5) as an estimator of for finite , has the advantages that the overlap problem is solved by simulating and the complex action problem reduced, since is found to scale for relatively small and it is possible to extrapolate solutions to larger values of . Moreover, the numerical errors of the solution do not grow exponentially with [6, 7].
3 Simulations and Results
Monte Carlo simulations are performed on the system where and , are parameters. The rational hybrid Monte Carlo method (RHMC) is used in the simulations. We take large enough, so that is sharply peaked at and the results are independent of . We use the estimators , and .
For we obtain the phase quenched model . We simulate this system and show the results for , in the left plot of Fig. 1. In the large limit, these values converge to the same value as predicted by GEM . No SSB of SO(6) occurs in the absence of a fluctuating phase as expected [2, 5].
In order to compute the solution of (5) in the large limit,
it is important to use the scaling properties of the functions
and for small and in the large limit.
We find that
, is almost constant for large . Thus the
extended direction becomes decorrelated with the phase and there is no
need to constrain the large eigenvalues as in . We also
observe good scaling of for . For
solving (5), the small scaling behavior 
is important. We compute the coefficients and
, as shown in the left plot of Fig. 2. For
the expected small scaling is
and implies the existence of a hard core
potential suppressing the small region, as it was first found in
. This is a strictly non perturbative effect
which is absent e.g. in the one loop
approximation. Around a scaling is found to dominate
Finally we attempt to compare the free energies and by numerically computing . Their value at the solution of Eq. (5) compares the free energies of the SO(3) and SO(2) vacua respectively. The calculation is done as described in  by computing
at . The term vanishes in the large limit. The comparison of the free energies can be read off the right plot of Fig. 3 where the results are shown. Although we cannot draw a definite conclusion, we see no inconsistency with a dominating SO(3) vacuum as predicted by GEM.
In summary, we simulated from first principles a 6 dimensional version of the IIB matrix model. We studied the scenario of dynamical compactification by SSB of the eigenvalues of the tensor and found that our results are consistent with the quantitative predictions of the GEM. The importance of the fluctuating phase in inducing SSB was demonstrated first by confirming that no SSB occurs in the phase quenched model. A strong complex action problem was found to hinder ordinary Monte Carlo simulations and the factorization method was used in order to determine and sample the region favored by the competing effects of the fluctuating phase, the real part of the action and the density of states. The distribution function of the eigenvalues were considered and by using the scaling properties of the factors and we computed the scale of the compactified dimensions. The results were found to be consistent with the unique –independent numerical value computed using the GEM and with the constant volume property . A strong, non perturbative, hard core potential against the collapse of the eigenvalues is shown to be generated dynamically which plays a crucial role in obtaining non trivial solutions in the large limit, generating scales and that are comparable and giving finite spacetime volume.
Recent results in , however, have shown that
it is possible to study the Lorentzian IIB matrix model and
obtain an expanding, large, 3 dimensional space which arises after a
- Note that the actual length of each dimension is, by definition, and respectively.
- Albeit a non classical spacetime with a fuzzy geometry since the dominant configurations of cannot be simultaneously diagonalizable .
- The constant is not necessary in the calculations described below.
- Note that in the one loop approximation, the scaling holds for all values of .
- In fact, it turned out that the complex action problem can be totally avoided in the Lorentzian model in a highly nontrivial manner .
- N. Ishibashi, H. Kawai, Y. Kitazawa and A. Tsuchiya, Nucl. Phys. B 498, 467 (1997); H. Aoki, S. Iso, H. Kawai, Y. Kitazawa and T. Tada, Prog. Theor. Phys. 99 (1998) 713;for a review see e.g. H. Aoki et. al., Prog. Theor. Phys. Suppl. 134 (1999) 47; T. Azuma, arXiv:hep-th/0401120.
- J. Nishimura, G. Vernizzi, JHEP 0004, 015 (2000); Phys. Rev. Lett. 85, 4664-4667 (2000).
- J. Nishimura and F. Sugino, JHEP 0205 (2002) 001; H. Kawai et. al., Nucl. Phys. B 647, 153 (2002); Prog. Theor. Phys. 109, 115 (2003); T. Aoyama, J. Nishimura, T. Okubo, Prog. Theor. Phys. 125, 537-563 (2011); J. Nishimura, T. Okubo, F. Sugino, arXiv:1108.1293 [hep-th]; Prog. Theor. Phys. 114, 487 (2005).
- W. Krauth, H. Nicolai and M. Staudacher, Phys. Lett. B 431 (1998) 31; W. Krauth and M. Staudacher, Phys. Lett. B 435 (1998) 350; Phys. Lett. B 453 (1999) 253; T. Hotta, J. Nishimura and A. Tsuchiya, Nucl. Phys. B 545 (1999) 543.
- J. Ambjørn, K. N. Anagnostopoulos, W. Bietenholz, T. Hotta and J. Nishimura, JHEP 0007 (2000) 013; JHEP 0007 (2000) 011; J. Ambjørn, K. N. Anagnostopoulos, W. Bietenholz, F. Hofheinz and J. Nishimura, Phys. Rev. D 65 (2002) 086001;
- K.N. Anagnostopoulos, J. Nishimura, Phys. Rev. D 65, 105012 (2002).
- J. Ambjørn, K. N. Anagnostopoulos, J. Nishimura and J. J. M. Verbaarschot, JHEP 0210 (2002) 062; Phys. Rev. D 70 (2004) 035010; V. Azcoiti, G. Di Carlo, A. Galante and V. Laliena, Phys. Rev. Lett. 89 (2002) 141601; Z. Fodor, S. D. Katz and C. Schmidt, JHEP 0703 (2007) 121; S. Ejiri, Phys. Rev. D 77, 014508 (2008); K. N. Anagnostopoulos, T. Azuma, J. Nishimura, PoS LATTICE2010, 167 (2010).
- K.N. Anagnostopoulos, T. Azuma, J. Nishimura, Phys. Rev. D83, 054504 (2011); J. High Energy Phys. 10 (2011) 126.
- S. -W. Kim, J. Nishimura and A. Tsuchiya, Phys. Rev. Lett. 108 (2012) 011601; Phys. Rev. D 86 (2012) 027901; JHEP 1210 (2012) 147; PoS (Lattice2012) 230.