BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jun Nishimura

Jun Nishimura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-thhep-latgr-qccond-mat.stat-mechhep-phphysics.comp-phastro-ph.COcond-mat.str-elastro-phastro-ph.HEnucl-thquant-ph

Catalog footprint

What is connected

69works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

69 published item(s)

preprint2023arXiv

Numerical studies on the finite-temperature CP restoration in 4D SU(N) gauge theory at $θ=π$

Recent studies on the 't Hooft anomaly matching condition have suggested a nontrivial phase structure in 4D SU($N$) gauge theory at $θ=π$. In the large-$N$ limit, it has been found that CP symmetry at $θ=π$ is broken in the confined phase, while it restores in the deconfined phase, which is indeed one of the possible scenarios. However, at small $N$, one may find other situations that are consistent with the consequence of the anomaly matching condition. Here we investigate this issue for $N=2$ by direct lattice calculations. The crucial point to note is that the CP restoration can be probed by the sudden change of the tail of the topological charge distribution at $θ=0$, which can be seen by simulating the theory at imaginary $θ$ without the sign problem. Our results suggest that the CP restoration at $θ=π$ occurs at temperature higher than the deconfining temperature unlike the situation in the large-$N$ limit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Backpropagating Hybrid Monte Carlo algorithm for fast Lefschetz thimble calculations

The Picard-Lefschetz theory has been attracting much attention as a tool to evaluate a multi-variable integral with a complex weight, which appears in various important problems in theoretical physics. The idea is to deform the integration contour based on Cauchy's theorem using the so-called gradient flow equation. In this paper, we propose a fast Hybrid Monte Carlo algorithm for evaluating the integral, where we "backpropagate" the force of the fictitious Hamilton dynamics on the deformed contour to that on the original contour, thereby reducing the required computational cost by a factor of the system size. Our algorithm can be readily extended to the case in which one integrates over the flow time in order to solve not only the sign problem but also the ergodicity problem that occurs when there are more than one thimbles contributing to the integral. This enables, in particular, efficient identification of all the dominant saddle points and the associated thimbles. We test our algorithm by calculating the real-time evolution of the wave function using the path integral formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Complex Langevin studies of the emergent space-time in the type IIB matrix model

The type IIB matrix model has been proposed as a non-perturbative definition of superstring theory since 1996. We study a simplified model that describes the late time behavior of the type IIB matrix model non-perturbatively using Monte Carlo methods, and we use the complex Langevin method to overcome the sign problem. We investigate a scenario where the space-time signature changes dynamically from Euclidean at early times to Lorentzian at late times. We discuss the possibility of the emergence of the (3+1)D expanding universe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Signature change of the emergent space-time in the IKKT matrix model

The IKKT matrix model (or the type IIB matrix model) is known as a promising candidate for a nonperturbative formulation of superstring theory in ten dimensions. As a most attractive feature, the model admits the emergence of (3+1)-dimensional space-time associated with the spontaneous breaking of the (9+1)-dimensional Lorentz symmetry. Numerical confirmation of such a phenomenon has been attempted for more than two decades. Recently it has been found that the sign problem, the main obstacle in simulating this model, can be overcome by the complex Langevin method. It has been shown that the Lorenzian version of the model is smoothly connected with the Euclidean version, in which the SO(10) symmetry is found to be spontaneously broken to SO(3). Here we propose to add a Lorentz invariant "mass" term to the original model and discuss a scenario that (3+1)-dimensional expanding space-time with Lorentzian signature appears at late times. Some numerical results supporting this scenario are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A new phase in the Lorentzian type IIB matrix model and the emergence of continuous space-time

The Lorentzian type IIB matrix model is a promising candidate for a non-perturbative formulation of superstring theory. In previous studies, Monte Carlo calculations provided interesting results indicating the spontaneous breaking of SO(9) to SO(3) and the emergence of (3+1)-dimensional space-time. However, an approximation was used to avoid the sign problem, which seemed to make the space-time structure singular. In this talk, we report our results obtained by using the complex Langevin method to overcome the sign problem instead of using this approximation. In particular, we discuss the emergence of continuous space-time in a new phase, which we discovered recently.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Relationship between the Euclidean and Lorentzian versions of the type IIB matrix model

The type IIB matrix model was proposed as a non-perturbative formulation of superstring theory in 1996. We simulate a model that describes the late time behavior of the IIB matrix model by applying the complex Langevin method to overcome the sign problem. We clarify the relationship between the Euclidean and the Lorentzian versions of the type IIB matrix model in a recently discovered phase. By introducing a constraint, we obtain a model where the spacetime metric is Euclidean at early times, whereas it {\it dynamically} becomes Lorentzian at late times.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Complex Langevin analysis of 2D U(1) gauge theory on a torus with a $θ$ term

Monte Carlo simulation of gauge theories with a $θ$ term is known to be extremely difficult due to the sign problem. Recently there has been major progress in solving this problem based on the idea of complexifying dynamical variables. Here we consider the complex Langevin method (CLM), which is a promising approach for its low computational cost. The drawback of this method, however, is the existence of a condition that has to be met in order for the results to be correct. As a first step, we apply the method to 2D U(1) gauge theory on a torus with a $θ$ term, which can be solved analytically. We find that a naive implementation of the method fails because of the topological nature of the $θ$ term. In order to circumvent this problem, we simulate the same theory on a punctured torus, which is equivalent to the original model in the infinite volume limit for $ |θ| < π$. Rather surprisingly, we find that the CLM works and reproduces the exact results for a punctured torus even at large $θ$, where the link variables near the puncture become very far from being unitary.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Complex Langevin analysis of the spontaneous breaking of 10D rotational symmetry in the Euclidean IKKT matrix model

The IKKT matrix model is a promising candidate for a nonperturbative formulation of superstring theory, in which spacetime is conjectured to emerge dynamically from the microscopic matrix degrees of freedom in the large-$N$ limit. Indeed in the Lorentzian version, Monte Carlo studies suggested the emergence of (3+1)-dimensional expanding space-time. Here we study the Euclidean version instead, and investigate an alternative scenario for dynamical compactification of extra dimensions via the spontaneous symmetry breaking (SSB) of 10D rotational symmetry. We perform numerical simulations based on the complex Langevin method (CLM) in order to avoid a severe sign problem. Furthermore, in order to avoid the singular-drift problem in the CLM, we deform the model and determine the SSB pattern as we vary the deformation parameter. From these results, we conclude that the original model has an SO(3) symmetric vacuum, which is consistent with previous results obtained by the Gaussian expansion method (GEM). We also apply the GEM to the deformed matrix model and find consistency with the results obtained by the CLM.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Complex Langevin studies of the dynamical compactification of extra dimensions in the Euclidean IKKT matrix model

The type IIB matrix model, also known as the IKKT matrix model, is a promising candidate for a nonperturbative formulation of superstring theory. In this talk we study the Euclidean version of the IKKT matrix model, which has a "sign problem" due to the Pfaffian coming from integrating out the fermionic degrees of freedom. To study the spontaneous breaking of the SO(10) rotational symmetry, we apply the Complex Langevin Method (CLM) to the Euclidean IKKT matrix model. We conclude that the SO(10) symmetry is broken to SO(3), in agreement with the previous studies by the Gaussian Expansion Method (GEM). We also apply the GEM to the deformed model and find consistency with the CLM result. These are proceedings of Takehiro Azuma's talk at Asia-Pacific Symposium for Lattice Field Theory (APLAT 2020) on August 4-7, 2020, based on the paper arXiv:2002.07410.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical Compactification of Extra Dimensions in the Euclidean IKKT Matrix Model via Spontaneous Symmetry Breaking

The IKKT matrix model has been conjectured to provide a promising nonperturbative formulation of superstring theory. In this model, spacetime emerges dynamically from the microscopic matrix degrees of freedom in the large-N limit, and Monte Carlo simulations of the Lorentzian version provide evidence of an emergent (3+1)-dimensional expanding space-time. In this talk, we discuss the Euclidean version of the IKKT matrix model and provide evidence of dynamical compactification of the extra dimensions via the spontaneous symmetry breaking (SSB) of the 10D rotational symmetry. We perform numerical simulations of a system with a severe complex action problem by using the complex Langevin method (CLM). The CLM suffers from the singular-drift problem and we deform the model in order to avoid it. We study the SSB pattern as we vary the deformation parameter and we conclude that the original model has an SO(3) symmetric vacuum, in agreement with previous calculations using the Gaussian expansion method (GEM). We employ the GEM to the deformed model and we obtain results consistent with the ones obtained by CLM.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

New perspectives on the emergence of (3+1)D expanding space-time in the Lorentzian type IIB matrix model

The type IIB matrix model is a promising candidate for a nonperturbative formulation of superstring theory. In the Lorentzian version, in particular, the emergence of (3+1)D expanding space-time was observed by Monte Carlo studies of this model. Here we provide new perspectives on the (3+1)D expanding space-time that have arised from recent studies. First it was found that the matrix configurations generated by the simulation are singular in that the submatrices representing the expanding 3D space have only two large eigenvalues associated with the Pauli matrices. This problem was conjectured to occur due to the approximation used to avoid the sign problem in simulating the model. In order to confirm this conjecture, the complex Langevin method was applied to overcome the sign problem instead of using the approximation. The results indeed showed a clear departure from the Pauli-matrix structure, while the (3+1)D expanding behavior remained unaltered. It was also found that classical solutions obtained within a certain ansatz show quite generically a (3+1)D expanding behavior with smooth space-time structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The emergence of expanding space-time and intersecting D-branes from classical solutions in the Lorentzian type IIB matrix model

The type IIB matrix model is a promising candidate for a nonperturbative formulation of superstring theory. As such, it is expected to explain the origin of space--time and matter at the same time. This has been partially demonstrated by the previous Monte Carlo studies on the Lorentzian version of the model, which suggested the emergence of (3+1)-dimensional expanding space--time. Here we investigate the same model by solving numerically the classical equation of motion, which is expected to be valid at late times since the action becomes large due to the expansion of space. Many solutions are obtained by the gradient descent method starting from random matrix configurations, assuming a quasi-direct-product structure for the (3+1)-dimensions and the extra 6 dimensions. We find that these solutions generally admit the emergence of expanding space--time and a block-diagonal structure in the extra dimensions, the latter being important for the emergence of intersecting D-branes. For solutions corresponding to D-branes with appropriate dimensionality, the Dirac operator is shown to acquire a zero mode in the limit of infinite matrix size.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gauge cooling for the singular-drift problem in the complex Langevin method --- an application to finite density QCD

We study full QCD at finite density and low temperature with light quark mass using the complex Langevin method. Since the singular drift problem turns out to be mild on a $4^3 \times 8$ lattice we use, the gauge cooling is performed only to control the unitarity norm in this exploratory study. We report on our preliminary data obtained from the complex Langevin simulation up to certain Langevin time. While the data are still noisy due to lack of statistics, the onset of the baryon number density seems to occur at larger $μ$ than half the pion mass, which is the value for the phase quenched QCD. The validity of our simulation is tested by the recently proposed criterion based on the probability distribution of the drift term.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gauge cooling for the singular-drift problem in the complex Langevin method --a test in Random Matrix Theory for finite density QCD

Recently, the complex Langevin method has been applied successfully to finite density QCD either in the deconfinement phase or in the heavy dense limit with the aid of a new technique called the gauge cooling. In the confinement phase with light quarks, however, convergence to wrong limits occurs due to the singularity in the drift term caused by small eigenvalues of the Dirac operator including the mass term. We propose that this singular-drift problem should also be overcome by the gauge cooling with different criteria for choosing the complexified gauge transformation. The idea is tested in chiral Random Matrix Theory for finite density QCD, where exact results are reproduced at zero temperature with light quarks. It is shown that the gauge cooling indeed changes drastically the eigenvalue distribution of the Dirac operator measured during the Langevin process. Despite its non-holomorphic nature, this eigenvalue distribution has a universal diverging behavior at the origin in the chiral limit due to a generalized Banks-Casher relation as we confirm explicitly.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large-scale computation of the exponentially expanding universe in a simplified Lorentzian type IIB matrix model

The type IIB matrix model is a conjectured nonperturbative formulation of superstring theory. Recent studies on the Lorentzian version of the model have shown that only three out of nine spatial directions start to expand after some critical time. On the other hand, due to the unbounded action of the Lorentzian model, one has to introduce infrared (IR) cutoffs in order to make the partition function finite. In this work we investigate whether the effects of the IR cutoffs disappear in the infinite volume limit. For that purpose, we study a simplified model with large matrix size up to $N=256$ by Monte Carlo simulation. First we confirm the exponentially expanding behavior of the "universe". Then we generalize the form of the IR cutoffs by one parameter, and find that the results become universal in some region of the parameter. It is suggested that the effects of IR cutoffs disappear in this region, which is confirmed also from the studies of Schwinger-Dyson equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical tests of the gauge/gravity duality conjecture for D0-branes at finite temperature and finite N

According to the gauge/gravity duality conjecture, the thermodynamics of gauge theory describing D-branes corresponds to that of black branes in superstring theory. We test this conjecture directly in the case of D0-branes by applying Monte Carlo methods to the corresponding gauge theory, which takes the form of the BFSS matrix quantum mechanics. In particular, we take the continuum limit by extrapolating the UV cutoff to infinity. First we perform simulations at large N so that string loop corrections can be neglected on the gravity side. Our results for the internal energy exhibit the temperature dependence consistent with the prediction including the α' corrections. Next we perform simulations at small N but at lower temperature so that the α' corrections can be neglected on the gravity side. Our results are consistent with the prediction including the leading string loop correction, which suggests that the conjecture holds even at finite N.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the condition for correct convergence in the complex Langevin method

The complex Langevin method (CLM) provides a promising way to perform the path integral with a complex action using a stochastic equation for complexified dynamical variables. It is known, however, that the method gives wrong results in some cases, while it works, for instance, in finite density QCD in the deconfinement phase or in the heavy dense limit. Here we revisit the argument for justification of the CLM and point out a subtlety in using the time-evolved observables, which play a crucial role in the argument. This subtlety requires that the probability distribution of the drift term should fall off exponentially or faster at large magnitude. We demonstrate our claim in some examples such as chiral Random Matrix Theory and show that our criterion is indeed useful in judging whether the results obtained by the CLM are trustable or not.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spontaneous symmetry breaking induced by complex fermion determinant --- yet another success of the complex Langevin method

In many interesting systems, the fermion determinant becomes complex and its phase plays a crucial role in the determination of the vacuum. For instance, in finite density QCD at low temperature and high density, exotic fermion condensates are conjectured to form due to such effects. When one applies the complex Langevin method to such a complex action system naively, one cannot obtain the correct results because of the singular-drift problem associated with the appearance of small eigenvalues of the Dirac operator. Here we propose to avoid this problem by adding a fermion bilinear term to the action and extrapolating its coefficient to zero. We test this idea in an SO(4)-invariant matrix model with a Gaussian action and a complex fermion determinant, whose phase is expected to induce the spontaneous breaking of the SO(4) symmetry. Our results agree well with the previous results obtained by the Gaussian expansion method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The argument for justification of the complex Langevin method and the condition for correct convergence

The complex Langevin method is a promising approach to the complex-action problem based on a fictitious time evolution of complexified dynamical variables under the influence of a Gaussian noise. Although it is known to have a restricted range of applicability, the use of gauge cooling made it applicable to various interesting cases including finite density QCD in certain parameter regions. In this paper, we revisit the argument for justification of the method. In particular, we point out a subtlety in the use of time-evolved observables, which play a crucial role in the previous argument. This requires that the probability of the drift term should fall off exponentially or faster at large magnitude. We argue that this is actually a necessary and sufficient condition for the method to be justified. Using two simple examples, we show that our condition tells us clearly whether the results obtained by the method are trustable or not. We also discuss a new possibility for the gauge cooling, which can reduce the magnitude of the drift term directly.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The complex Langevin analysis of spontaneous symmetry breaking induced by complex fermion determinant

In many interesting physical systems, the determinant which appears from integrating out fermions becomes complex, and its phase plays a crucial role in the determination of the vacuum. An example of this is QCD at low temperature and high density, where various exotic fermion condensates are conjectured to form. Another example is the Euclidean version of the type IIB matrix model for 10d superstring theory, where spontaneous breaking of the SO(10) rotational symmetry down to SO(4) is expected to occur. When one applies the complex Langevin method to these systems, one encounters the singular-drift problem associated with the appearance of nearly zero eigenvalues of the Dirac operator. Here we propose to avoid this problem by deforming the action with a fermion bilinear term. The results for the original system are obtained by extrapolations with respect to the deformation parameter. We demonstrate the power of this approach by applying it to a simple matrix model, in which spontaneous symmetry breaking from SO(4) to SO(2) is expected to occur due to the phase of the complex fermion determinant. Unlike previous work based on a reweighting-type method, we are able to determine the true vacuum by calculating the order parameters, which agree with the prediction by the Gaussian expansion method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A renormalization group method for studying the early universe in the Lorentzian IIB matrix model

We propose a new method for studying the early universe in the Lorentzian version of the IIB matrix model, which is considered to be a nonperturbative formulation of superstring theory. This method is based on the idea of renormalization group, and it enables us to study the time-evolution of the universe for much longer time than in the previous work, which showed that the SO(9) rotational symmetry is spontaneously broken down to SO(3) after a "critical time". We demonstrate how this method works in a simplified model, which is expected to capture the behaviors of the original model when the space is not so large. In particular, we present clear evidence that the three-dimensional space expands exponentially after the critical time in this simplified model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Justification of the complex Langevin method with the gauge cooling procedure

Recently there has been remarkable progress in the complex Langevin method, which aims at solving the complex action problem by complexifying the dynamical variables in the original path integral. In particular, a new technique called the gauge cooling was introduced and the full QCD simulation at finite density has been made possible in the high temperature (deconfined) phase or with heavy quarks. Here we provide a rigorous justification of the complex Langevin method including the gauge cooling procedure. We first show that the gauge cooling can be formulated as an extra term in the complex Langevin equation involving a gauge transformation parameter, which is chosen appropriately as a function of the configuration before cooling. The probability distribution of the complexified dynamical variables is modified by this extra term. However, this modification is shown not to affect the Fokker-Planck equation for the corresponding complex weight as far as observables are restricted to gauge invariant ones. Thus we demonstrate explicitly that the gauge cooling can be used as a viable technique to satisfy the convergence conditions for the complex Langevin method. We also discuss the "gauge cooling" in 0-dimensional systems such as vector models or matrix models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo studies of dynamical compactification of extra dimensions in a model of nonperturbative string theory

The IIB matrix model has been proposed as a non-perturbative definition of superstring theory. In this work, we study the Euclidean version of this model in which extra dimensions can be dynamically compactified if a scenario of spontaneously breaking the SO(10) rotational symmetry is realized. Monte Carlo calculations of the Euclidean IIB matrix model suffer from a very strong complex action problem due to the large fluctuations of the complex phase of the Pfaffian which appears after integrating out the fermions. We employ the factorization method in order to achieve effective sampling. We report on preliminary results that can be compared with previous studies of the rotational symmetry breakdown using the Gaussian expansion method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New insights into the problem with a singular drift term in the complex Langevin method

The complex Langevin method aims at performing path integral with a complex action numerically based on complexification of the original real dynamical variables. One of the poorly understood issues concerns occasional failure in the presence of logarithmic singularities in the action, which appear, for instance, from the fermion determinant in finite density QCD. We point out that the failure should be attributed to the breakdown of the relation between the complex weight that satisfies the Fokker-Planck equation and the probability distribution associated with the stochastic process. In fact, this problem can occur in general when the stochastic process involves a singular drift term. We show, however, in a simple example that there exists a parameter region in which the method works although the standard reweighting method is hardly applicable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Power-law expansion of the Universe from the bosonic Lorentzian type IIB matrix model

Recent studies on the Lorentzian version of the type IIB matrix model show that (3+1)D expanding universe emerges dynamically from (9+1)D space-time predicted by superstring theory. Here we study a bosonic matrix model obtained by omitting the fermionic matrices. With the adopted simplification and the usage of a large-scale parallel computer, we are able to perform Monte Carlo calculations with matrix size up to $N=512$, which is twenty times larger than that used previously for the studies of the original model. When the matrix size is larger than some critical value $N_{\rm c}\simeq 110$, we find that (3+1)D expanding universe emerges dynamically with a clear large-$N$ scaling property. Furthermore, at sufficiently late times, we observe a power-law behavior $t^{1/2}$ of the spatial extent with respect to time $t$, which is reminiscent of the expanding behavior of the Friedmann-Robertson-Walker universe in the radiation dominated era. We discuss possible implications of this result on the original model including fermionic matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing a generalized cooling procedure in the complex Langevin simulation of chiral Random Matrix Theory

The complex Langevin method has been attracting much attention as a solution to the sign problem since the method was shown to work in finite density QCD in the deconfined phase by using the so-called gauge cooling procedure. Whether it works also in the confined phase with light quarks is still an open question, though. In order to shed light on this question, we apply the method to the chiral Random Matrix Theory, which describes the epsilon regime of finite density QCD. Earlier works reported that a naive implementation of the method fails to reproduce the known exact results and that the problem can be solved by choosing a suitable coordinate. In this work we stick to the naive implementation, and show that a generalized gauge cooling procedure can be used to avoid the problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Understanding the problem with logarithmic singularities in the complex Langevin method

In recent years, there has been remarkable progress in theoretical justification of the complex Langevin method, which is a promising method for evading the sign problem in the path integral with a complex weight. There still remains, however, an issue concerning occasional failure of this method in the case where the action involves logarithmic singularities such as the one appearing from the fermion determinant in finite density QCD. In this talk, we point out that this failure is due to the breakdown of the relation between the complex weight which satisfies the Fokker-Planck equation and the probability distribution generated by the stochastic process. In fact, this kind of failure can occur in general when the stochastic process involves a singular drift term. We show, however, in simple examples, that there exists a parameter region in which the method works although the standard reweighting method is hardly applicable.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Realizing three generations of the Standard Model fermions in the type IIB matrix model

We discuss how the Standard Model particles appear from the type IIB matrix model, which is considered to be a nonperturbative formulation of superstring theory. In particular, we are concerned with a constructive definition of the theory, in which we start with finite-N matrices and take the large-N limit afterwards. In that case, it was pointed out recently that realizing chiral fermions in the model is more difficult than it had been thought from formal arguments at N=infinity and that introduction of a matrix version of the warp factor is necessary. Based on this new insight, we show that two generations of the Standard Model fermions can be realized by considering a rather generic configuration of fuzzy S^2 and fuzzy S^2 * S^2 in the extra dimensions. We also show that three generations can be obtained by squashing one of the S^2's that appear in the configuration. Chiral fermions appear at the intersections of the fuzzy manifolds with nontrivial Yukawa couplings to the Higgs field, which can be calculated from the overlap of their wave functions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recent developments in the type IIB matrix model

We review recent developments in the type IIB matrix model, which was conjectured to be a nonperturbative formulation of superstring theory. In the first part we review the recent results for the Euclidean model, which suggest that SO(10) symmetry is spontaneously broken. In the second part we review the recent results for the Lorentzian model. In particular, we discuss Monte Carlo results, which suggest that (3+1)-dimensional expanding universe emerges dynamically. We also discuss some results suggesting the emergence of exponential expansion and the power-law expansion at later times. The behaviors at much later times are studied by the classical equation of motion. We discuss a solution representing 3d expanding space, which suggests a possible solution to the cosmological constant problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct test of the AdS/CFT correspondence by Monte Carlo studies of N=4 super Yang-Mills theory

We perform nonperturbative studies of N=4 super Yang-Mills theory by Monte Carlo simulation. In particular, we calculate the correlation functions of chiral primary operators to test the AdS/CFT correspondence. Our results agree with the predictions obtained from the AdS side that the SUSY non-renormalization property is obeyed by the three-point functions but \emph{not} by the four-point functions investigated in this paper. Instead of the lattice regularization, we use a novel regularization of the theory based on an equivalence in the large-N limit between the N=4 SU(N) theory on RxS^3 and a one-dimensional SU(N) gauge theory known as the plane-wave (BMN) matrix model. The equivalence extends the idea of large-N reduction to a curved space and, at the same time, overcomes the obstacle related to the center symmetry breaking. The adopted regularization preserves 16 SUSY, which is crucial in testing the AdS/CFT correspondence with the available computer resources. The only SUSY breaking effects, which come from the momentum cutoff $Λ$ in R direction, are made negligible by using sufficiently large $Λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Holographic description of quantum black hole on a computer

The discovery of the fact that black holes radiate particles and eventually evaporate led Hawking to pose the well-known information loss paradox. This paradox caused a long and serious debate since it claims that the fundamental laws of quantum mechanics may be violated. A possible cure appeared recently from superstring theory, a consistent theory of quantum gravity: if the holographic description of a quantum black hole based on the gauge/gravity duality is correct, the information is not lost and quantum mechanics remains valid. Here we test this gauge/gravity duality on a computer at the level of quantum gravity for the first time. The black hole mass obtained by Monte Carlo simulation of the dual gauge theory reproduces precisely the quantum gravity effects in an evaporating black hole. This result opens up totally new perspectives towards quantum gravity since one can simulate quantum black holes through dual gauge theories.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo studies of the spontaneous rotational symmetry breaking in dimensionally reduced super Yang-Mills models

It has long been speculated that the spontaneous symmetry breaking (SSB) of SO(D) occurs in matrix models obtained by dimensionally reducing super Yang-Mills theory in D=6,10 dimensions. In particular, the D=10 case corresponds to the IIB matrix model, which was proposed as a nonperturbative formulation of superstring theory, and the SSB may correspond to the dynamical generation of four-dimensional space-time. Recently, it has been shown by using the Gaussian expansion method that the SSB indeed occurs for D=6 and D=10, and interesting nature of the SSB common to both cases has been suggested. Here we study the same issue from first principles by a Monte Carlo method in the D=6 case. In spite of a severe complex-action problem, the factorization method enables us to obtain various quantities associated with the SSB, which turn out to be consistent with the previous results obtained by the Gaussian expansion method. This also demonstrates the usefulness of the factorization method as a general approach to systems with the complex-action problem or the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo studies on the expanding behavior of the early universe in the Lorentzian type IIB matrix model

The type IIB matrix model is a conjectured nonperturbative formulation of superstring theory. Recently the Lorentzian version of the model has been studied by Monte Carlo simulation, and it has been shown that only three out of nine spatial directions start to expand after a critical time. We extend this work by investigating the expanding behavior for much longer time. We find that the 3d space expands exponentially for some period of time, which may be interpreted as inflation. We also simulate a simplified model, which is expected to capture some qualitative features of the original model at much later times. We observe that the exponential expansion eventually changes into a power-law t^(1/2) behavior, which agrees with the expanding behavior of the Friedmann-Robertson-Walker (FRW) universe in the radiation dominated era.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Realizing chiral fermions in the type IIB matrix model at finite N

We discuss how chiral fermions can appear in the type IIB matrix model, which is considered to be a nonperturbative formulation of superstring theory. In particular, we are concerned with a constructive definition of the theory, in which we start with a finite-N configuration and take the large-N limit later on. We point out that there exists a certain necessary condition which the structure of the extra dimensions should satisfy. As an example, we consider a previous proposal using intersecting branes and show that chiral fermions can indeed be realized in four dimensions by introducing a matrix counterpart of warped space-time. This is remarkable in view of the well-known difficulty in realizing chiral fermions in lattice gauge theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

(3+1)-dimensional expanding universe from a Lorentzian matrix model for superstring theory in (9+1)-dimensions

We study the Lorentzian version of the type IIB matrix model as a nonperturbative formulation of superstring theory in (9+1)-dimensions. Monte Carlo results show that not only space but also time emerges dynamically in this model. Furthermore, the real-time dynamics extracted from the matrices turns out to be remarkable: 3 out of 9 spatial directions start to expand at some critical time. This can be interpreted as the birth of our Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Expanding (3+1)-dimensional universe from a Lorentzian matrix model for superstring theory in (9+1)-dimensions

We reconsider the matrix model formulation of type IIB superstring theory in (9+1)-dimensional space-time. Unlike the previous proposal in which the Wick rotation was used to make the model well-defined, we regularize the Lorentzian model by introducing infrared cutoffs in both the spatial and temporal directions. Monte Carlo studies reveal that the two cutoffs can be removed in the large-N limit and that the theory thus obtained has no parameters other than one scale parameter. Moreover, we find that three out of nine spatial directions start to expand at some "critical time", after which the space has SO(3) symmetry instead of SO(9).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Late time behaviors of the expanding universe in the IIB matrix model

Recently we have studied the Lorentzian version of the IIB matrix model as a nonperturbative formulation of superstring theory. By Monte Carlo simulation, we have shown that the notion of time ---as well as space---emerges dynamically from this model, and that we can uniquely extract the real-time dynamics, which turned out to be rather surprising: after some "critical time", the SO(9) rotational symmetry of the nine-dimensional space is spontaneously broken down to SO(3) and the three-dimensional space starts to expand rapidly. In this paper, we study the same model based on the classical equations of motion, which are expected to be valid at later times. After providing a general prescription to solve the equations, we examine a class of solutions, which correspond to manifestly commutative space. In particular, we find a solution with an expanding behavior that naturally solves the cosmological constant problem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local field theory from the expanding universe at late times in the IIB matrix model

Recently we have shown that (3+1)-dimensional expanding universe appears dynamically and uniquely from the Lorentzian version of the IIB matrix model, which is considered as a nonperturbative formulation of superstring theory. Similarly, it is possible that the Standard Model appears uniquely at the electroweak scale from the same model at late times. In order to pursue such a possibility, we discuss how to derive the effective local field theory for the massless modes that appear at late times in the same formulation. As a concrete example, we consider the massless modes associated with the spontaneous breaking of (9+1)-dimensional Poincare symmetry and supersymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Monte Carlo studies of 3d N=6 SCFT via localization method

We perform Monte Carlo study of the 3d N=6 superconformal U(N)*U(N) Chern-Simons gauge theory (ABJM theory), which is conjectured to be dual to M-theory or type IIA superstring theory on certain AdS backgrounds. Our approach is based on a localization method, which reduces the problem to the simulation of a simple matrix model. This enables us to circumvent the difficulties in the original theory such as the sign problem and the SUSY breaking on a lattice. The new approach opens up the possibility of probing the quantum aspects of M-theory and testing the AdS_4/CFT_3 duality at the quantum level. Here we calculate the free energy, and confirm the N^{3/2} scaling in the M-theory limit predicted from the gravity side. We also find that our results nicely interpolate the analytical formulae proposed previously in the M-theory and type IIA regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical studies of the ABJM theory for arbitrary N at arbitrary coupling constant

We show that the ABJM theory, which is an N=6 superconformal U(N)*U(N) Chern-Simons gauge theory, can be studied for arbitrary N at arbitrary coupling constant by applying a simple Monte Carlo method to the matrix model that can be derived from the theory by using the localization technique. This opens up the possibility of probing the quantum aspects of M-theory and testing the AdS_4/CFT_3 duality at the quantum level. Here we calculate the free energy, and confirm the N^{3/2} scaling in the M-theory limit predicted from the gravity side. We also find that our results nicely interpolate the analytical formulae proposed previously in the M-theory and type IIA regimes. Furthermore, we show that some results obtained by the Fermi gas approach can be clearly understood from the constant map contribution obtained by the genus expansion. The method can be easily generalized to the calculations of BPS operators and to other theories that reduce to matrix models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observations of High Energy Cosmic-Ray Electrons from 30 GeV to 3 TeV with Emulsion Chambers

We have performed a series of cosmic-ray electron observations using the balloon-borne emulsion chambers since 1968. While we previously reported the results from subsets of the exposures, the final results of the total exposures up to 2001 are presented here. Our successive experiments have yielded the total exposure of 8.19 m^2 sr day at the altitudes of 4.0 - 9.4 g/cm^2. The performance of the emulsion chambers was examined by accelerator beam tests and Monte-Carlo simulations, and the on-board calibrations were carried out by using the flight data. In this work we present the cosmic-ray electron spectrum in the energy range from 30 GeV to 3 TeV at the top of the atmosphere, which is well represented by a power-law function with an index of -3.28+-0.10. The observed data can be also interpreted in terms of diffusive propagation models. The evidence of cosmic-ray electrons up to 3 TeV suggests the existence of cosmic-ray electron sources at distances within ~1 kpc and times within ~1x10^5 yr ago.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The origin of space-time as seen from matrix model simulations

The AdS/CFT correspondence, or more generally the gauge/gravity duality, is a remarkable conjecture obtained from superstring theory with various D-brane backgrounds. According to this conjecture, a higher-dimensional curved space-time emerges from supersymmetric gauge theory in lower-dimensional flat space-time. In the first part of this article, we review Monte Carlo studies of U(N) supersymmetric gauge theories, which confirmed the gauge/gravity duality for various observables. In particular, Monte Carlo results for thermodynamic quantities enable us to understand the microscopic origin of the black hole entropy associated with the dual geometry. We also discuss results for Wilson loops and correlation functions, which agree nicely with the predictions from the gravity side. In the second part, we review recent developments in a nonperturbative formulation of superstring theory, which may be regarded as a counterpart of the lattice gauge theory in QCD. In particular, we discuss Monte Carlo results for the Lorentzian matrix model, which suggest that (3+1)-dimensional expanding universe emerges dynamically from type IIB superstring theory in (9+1) dimensions if one treats the theory nonperturbatively.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A general approach to the sign problem - the factorization method with multiple observables

The sign problem is a notorious problem, which occurs in Monte Carlo simulations of a system with the partition function whose integrand is not real positive. The basic idea of the factorization method applied on such a system is to control some observables in order to determine and sample efficiently the region of configuration space which gives important contribution to the partition function. We argue that it is crucial to choose appropriately the set of the observables to be controlled in order for the method to work successfully in a general system. This is demonstrated by an explicit example, in which it turns out to be necessary to control more than one observables. Extrapolation to large system size is possible due to the nice scaling properties of the factorized functions, and known results obtained by an analytic method are shown to be consistently reproduced.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A practical solution to the sign problem in a matrix model for dynamical compactification

The matrix model formulation of superstring theory offers the possibility to understand the appearance of 4d space-time from 10d as a consequence of spontaneous breaking of the SO(10) symmetry. Monte Carlo studies of this issue is technically difficult due to the so-called sign problem. We present a practical solution to this problem generalizing the factorization method proposed originally by two of the authors (K.N.A. and J.N.). Explicit Monte Carlo calculations and large-N extrapolations are performed in a simpler matrix model with similar properties, and reproduce quantitative results obtained previously by the Gaussian expansion method. Our results also confirm that the spontaneous symmetry breaking indeed occurs due to the phase of the fermion determinant, which vanishes for collapsed configurations. We clarify various generic features of this approach, which would be useful in applying it to other statistical systems with the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Direct test of the gauge-gravity correspondence for Matrix theory correlation functions

We study correlation functions in (0+1)-dimensional maximally supersymmetric U(N) Yang-Mills theory, which was proposed by Banks et al. as a non-perturbative definition of 11-dimensional M-theory in the infinite-momentum frame. We perform first-principle calculations using Monte Carlo simulations, and compare the results against the predictions obtained previously based on the gauge-gravity correspondence from 10 dimensions. After providing a self-contained review on these predictions, we present clear evidence that the predictions in the large-N limit actually hold even at small N such as N=2 and 3. The predicted behavior seems to continue to the far infrared regime, which goes beyond the naive range of validity of the 10D supergravity analysis. This suggests that the correlation functions also contain important information on the M-theory limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Expanding universe as a classical solution in the Lorentzian matrix model for nonperturbative superstring theory

Recently we have shown by Monte Carlo simulation that expanding (3+1)-dimensional universe appears dynamically from a Lorentzian matrix model for type IIB superstring theory in (9+1)-dimensions. The mechanism for the spontaneous breaking of rotational symmetry relies crucially on the noncommutative nature of the space. Here we study the classical equations of motion as a complementary approach. In particular, we find a unique class of SO(3) symmetric solutions, which exhibits the time-dependence compatible with the expanding universe. The space-space noncommutativity is exactly zero, whereas the space-time noncommutativity becomes significant only towards the end of the expansion. We interpret the Monte Carlo results and the classical solution as describing the behavior of the model at earlier time and at later time, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonperturbative studies of supersymmetric matrix quantum mechanics with 4 and 8 supercharges at finite temperature

We investigate thermodynamic properties of one-dimensional U(N) supersymmetric gauge theories with 4 and 8 supercharges in the planar large-N limit by Monte Carlo calculations. Unlike the 16 supercharge case, the threshold bound state with zero energy is widely believed not to exist in these models. This led A.V. Smilga to conjecture that the internal energy decreases exponentially at low temperature instead of decreasing with a power law. In the 16 supercharge case, the latter behavior was predicted from the dual black 0-brane geometry and confirmed recently by Monte Carlo calculations. Our results for the models with 4 and 8 supercharges indeed support the exponential behavior, revealing a qualitative difference from the 16 supercharge case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Systematic study of the SO(10) symmetry breaking vacua in the matrix model for type IIB superstrings

We study the properties of the space-time that emerges dynamically from the matrix model for type IIB superstrings in ten dimensions. We calculate the free energy and the extent of space-time using the Gaussian expansion method up to the third order. Unlike previous works, we study the SO(d) symmetric vacua with all possible values of d within the range $2 \le d \le 7$, and observe clear indication of plateaus in the parameter space of the Gaussian action, which is crucial for the results to be reliable. The obtained results indeed exhibit systematic dependence on d, which turns out to be surprisingly similar to what was observed recently in an analogous work on the six-dimensional version of the model. In particular, we find the following properties: i) the extent in the shrunken directions is given by a constant, which does not depend on d; ii) the ten-dimensional volume of the Euclidean space-time is given by a constant, which does not depend on d except for d = 2; iii) The free energy takes the minimum value at d = 3. Intuitive understanding of these results is given by using the low-energy effective theory and some Monte Carlo results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing the AdS/CFT correspondence by Monte Carlo calculation of BPS and non-BPS Wilson loops in 4d N=4 super-Yang-Mills theory

We test the AdS/CFT correspondence by calculating Wilson loops in N = 4 super Yang-Mills theory on R*S^3 in the planar limit. Our method is based on a novel large-N reduction, which reduces the problem to Monte Carlo calculations in the plane-wave matrix model or the BMN matrix model, which is a 1d gauge theory with 16 supercharges. By using the gauge-fixed momentumspace simulation, we obtain results respecting 16 supersymmetries. We report on the Monte Carlo results for the BPS circular Wilson loop, which reproduce the exact result up to strong coupling. As a future prospect, we calculate a track-shapedWilson loop from the gravity side, which shows that a clear test of the AdS/CFT for the non-BPS case is also feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards an Effective Importance Sampling in Monte Carlo Simulations of a System with a Complex Action

The sign problem is a notorious problem, which occurs in Monte Carlo simulations of a system with a partition function whose integrand is not positive. One way to simulate such a system is to use the factorization method where one enforces sampling in the part of the configuration space which gives important contribution to the partition function. This is accomplished by using constraints on some observables chosen appropriately and minimizing the free energy associated with their joint distribution functions. These observables are maximally correlated with the complex phase. Observables not in this set essentially decouple from the phase and can be calculated without the sign problem in the corresponding "microcanonical" ensemble. These ideas are applied on a simple matrix model with very strong sign problem and the results are found to be consistent with analytic calculations using the Gaussian Expansion Method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Study of the Complex Action Problem in a Simple Model for Dynamical Compactification in Superstring Theory Using the Factorization Method

The IIB matrix model proposes a mechanism for dynamically generating four dimensional space--time in string theory by spontaneous breaking of the ten dimensional rotational symmetry $\textrm{SO}(10)$. Calculations using the Gaussian expansion method (GEM) lend support to this conjecture. We study a simple $\textrm{SO}(4)$ invariant matrix model using Monte Carlo simulations and we confirm that its rotational symmetry breaks down, showing that lower dimensional configurations dominate the path integral. The model has a strong complex action problem and the calculations were made possible by the use of the factorization method on the density of states $ρ_n(x)$ of properly normalized eigenvalues $\tildeλ_n$ of the space--time moment of inertia tensor. We study scaling properties of the factorized terms of $ρ_n(x)$ and we find them in agreement with simple scaling arguments. These can be used in the finite size scaling extrapolation and in the study of the region of configuration space obscured by the large fluctuations of the phase. The computed values of $\tildeλ_n$ are in reasonable agreement with GEM calculations and a numerical method for comparing the free energy of the corresponding ansatze is proposed and tested.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-lattice simulation of supersymmetric gauge theories as a probe to quantum black holes and strings

In the past decade we have witnessed remarkable developments in the gauge-gravity duality, which suggested a new approach to superstring theory and quantum space-time. In this context it is important to study supersymmetric large-N gauge theories in the strongly coupled regime. I will summarize the results and insights obtained so far by non-lattice simulations. A simple example of the gauge-gravity duality is the one between 1d U(N) gauge theory with 16 supercharges and the so-called black 0-brane solution in type IIA supergravity. In order for this duality to be valid, one has to take the 't Hooft large-N limit and to take the strong coupling limit on the gauge theory side. The gauge theory can be regularized by fixing the gauge completely thanks to one dimension, and by introducing a Fourier mode cutoff. One can then use the standard RHMC algorithm to simulate the system. The energy calculated as a function of the temperature was compared with the results obtained from the gravity side based on the black hole thermodynamics. This confirmed the gauge-gravity duality with high accuracy and provided the microscopic origin of the black hole thermodynamics. From the calculation of the Wilson loop, one obtains the Schwarzschild radius of the dual geometry. One can actually use the present 1d model with supersymmetric mass deformation to study \mathcal{N}=4 super Yang-Mills theory on R \times S^3 based on a novel large-N reduction, which generalizes the original idea of Eguchi and Kawai. It is remarkable that we can now simulate the 4d superconformal field theory, which appears in the most typical case of the gauge-gravity duality known as the AdS/CFT correspondence. In particular, no fine-tuning is required unlike previous proposals based on the lattice regularization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spontaneous breaking of the rotational symmetry in dimensionally reduced super Yang-Mills models

We investigate the spontaneous breaking of the SO(D) symmetry in matrix models, which can be obtained by the zero-volume limit of pure SU(N) super Yang-Mills theory in D = 6, 10 dimensions. The D = 10 case corresponds to the IIB matrix model, which was proposed as a non-perturbative formulation of type IIB superstring theory, and the spontaneous breaking corresponds to the dynamical compactification of space-time suggested in that model. First we study the D = 6 case by the Gaussian expansion method, which turns out to yield clearer results than the previous results for the D = 10 case for certain technical reasons. By comparing the free energy of the SO(d) symmetric vacua for d = 2, 3, 4, 5, we conclude that the breaking SO(6) \to SO(3) actually occurs. We find that the extent of space-time in the shrunken directions is almost independent of d. In units of this universal scale, the extended directions seem to have large but still finite extents depending on d. We show that these results for the extent of space-time can be explained quantitatively by an argument based on the low-energy effective theory. With these new insights, we reconsider the previous results for the IIB matrix model, and find that they are also consistent with our argument based on the low-energy effective theory. Thus we arrive at comprehensive understanding and some quantitative predictions concerning the nature of the spontaneous symmetry breaking taking place in these models. The space-time picture that emerges from the IIB matrix model and its implication on possible interpretations of the model are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Supersymmetry non-renormalization theorem from a computer and the AdS/CFT correspondence

We perform Monte Carlo calculation of correlation functions in 4d N=4 super Yang-Mills theory on R*S^3 in the planar limit. In order to circumvent the well-known problem of lattice SUSY, we adopt the idea of a novel large-N reduction, which reduces the calculation to that of corresponding correlation functions in the plane-wave matrix model or the BMN matrix model. This model is a 1d gauge theory with 16 supersymmetries, which can be simulated in a manner similar to the recent studies of the D0-brane system. We study two-point and three-point functions of chiral primary operators at various coupling constant, and find that they agree with the free theory results up to overall constant factors. The ratio of the overall factors for two-point and three-point functions agrees with the prediction of the AdS/CFT correspondence.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dominance of a single topological sector in gauge theory on non-commutative geometry

We demonstrate a striking effect of non-commutative (NC) geometry on topological properties of gauge theory by Monte Carlo simulations. We study 2d U(1) NC gauge theory for various boundary conditions using a new finite-matrix formulation proposed recently. We find that a single topological sector dictated by the boundary condition dominates in the continuum limit. This is in sharp contrast to the results in commutative space-time based on lattice gauge theory, where all topological sectors appear with certain weights in the continuum limit. We discuss possible implications of this effect in the context of string theory compactifications and in field theory contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Monte Carlo studies of Matrix theory correlation functions

We study correlation functions in (0+1)-dimensional maximally supersymmetric U(N) gauge theory, which represents the low-energy effective theory of D0-branes. In the large-N limit, the gauge-gravity duality predicts power-law behaviors in the infrared region for the two-point correlation functions of operators corresponding to supergravity modes. We evaluate such correlation functions on the gauge theory side by the Monte Carlo method. Clear power-law behaviors are observed at N=3, and the predicted exponents are confirmed consistently. Our results suggest that the agreement extends to the M-theory regime, where the supergravity analysis in 10 dimensions may not be justified a priori.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Putting M theory on a computer

We propose a non-lattice simulation for studying supersymmetric matrix quantum mechanics in a non-perturbative manner. In particular, our method enables us to put M theory on a computer based on its matrix formulation proposed by Banks, Fischler, Shenker and Susskind. Here we present Monte Carlo results of the same matrix model but in a different parameter region, which corresponds to the 't Hooft large-N limit at finite temperature. In the strong coupling limit the model has a dual description in terms of the N D0-brane solution in 10d type IIA supergravity. Our results provide highly nontrivial evidences for the conjectured duality. In particular, the energy (and hence the entropy) of the non-extremal black hole has been reproduced by solving directly the strongly coupled dynamics of the D0-brane effective theory.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The instability of intersecting fuzzy spheres

We discuss the classical and quantum stability of general configurations representing many fuzzy spheres in dimensionally reduced Yang-Mills-Chern-Simons models with and without supersymmetry. By performing one-loop perturbative calculations around such configurations, we find that intersecting fuzzy spheres are classically unstable in the class of models studied in this paper. We also discuss the large-N limit of the one-loop effective action as a function of the distance of fuzzy spheres. This shows, in particular, that concentric fuzzy spheres with different radii, which are identified with the 't Hooft-Polyakov monopoles, are perturbatively stable in the bosonic model and in the D=10 supersymmetric model.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Monte Carlo approach to nonperturbative strings -- demonstration in noncritical string theory

We show how Monte Carlo approach can be used to study the double scaling limit in matrix models. As an example, we study a solvable hermitian one-matrix model with the double-well potential, which has been identified recently as a dual description of noncritical string theory with worldsheet supersymmetry. This identification utilizes the nonperturbatively stable vacuum unlike its bosonic counterparts, and therefore it provides a complete constructive formulation of string theory. Our data with the matrix size ranging from 8 to 512 show a clear scaling behavior, which enables us to extract the double scaling limit accurately. The ``specific heat'' obtained in this way agrees nicely with the known result obtained by solving the Painleve-II equation with appropriate boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Probability distribution of the index in gauge theory on 2d non-commutative geometry

We investigate the effects of non-commutative geometry on the topological aspects of gauge theory using a non-perturbative formulation based on the twisted reduced model. The configuration space is decomposed into topological sectors labeled by the index nu of the overlap Dirac operator satisfying the Ginsparg-Wilson relation. We study the probability distribution of nu by Monte Carlo simulation of the U(1) gauge theory on 2d non-commutative space with periodic boundary conditions. In general the distribution is asymmetric under nu -> -nu, reflecting the parity violation due to non-commutative geometry. In the continuum and infinite-volume limits, however, the distribution turns out to be dominated by the topologically trivial sector. This conclusion is consistent with the instanton calculus in the continuum theory. However, it is in striking contrast to the known results in the commutative case obtained from lattice simulation, where the distribution is Gaussian in a finite volume, but the width diverges in the infinite-volume limit. We also calculate the average action in each topological sector, and provide deeper understanding of the observed phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The index of the overlap Dirac operator on a discretized 2d non-commutative torus

The index, which is given in terms of the number of zero modes of the Dirac operator with definite chirality, plays a central role in various topological aspects of gauge theories. We investigate its properties in non-commutative geometry. As a simple example, we consider the U(1) gauge theory on a discretized 2d non-commutative torus, in which general classical solutions are known. For such backgrounds we calculate the index of the overlap Dirac operator satisfying the Ginsparg-Wilson relation. When the action is small, the topological charge defined by a naive discretization takes approximately integer values, and it agrees with the index as suggested by the index theorem. Under the same condition, the value of the index turns out to be a multiple of N, the size of the 2d lattice. By interpolating the classical solutions, we construct explicit configurations, for which the index is of order 1, but the action becomes of order N. Our results suggest that the probability of obtaining a non-zero index vanishes in the continuum limit, unlike the corresponding results in the commutative space.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A non-perturbative study of 4d U(1) non-commutative gauge theory -- the fate of one-loop instability

Recent perturbative studies show that in 4d non-commutative spaces, the trivial (classically stable) vacuum of gauge theories becomes unstable at the quantum level, unless one introduces sufficiently many fermionic degrees of freedom. This is due to a negative IR-singular term in the one-loop effective potential, which appears as a result of the UV/IR mixing. We study such a system non-perturbatively in the case of pure U(1) gauge theory in four dimensions, where two directions are non-commutative. Monte Carlo simulations are performed after mapping the regularized theory onto a U(N) lattice gauge theory in d=2. At intermediate coupling strength, we find a phase in which open Wilson lines acquire non-zero vacuum expectation values, which implies the spontaneous breakdown of translational invariance. In this phase, various physical quantities obey clear scaling behaviors in the continuum limit with a fixed non-commutativity parameter $θ$, which provides evidence for a possible continuum theory. The extent of the dynamically generated space in the non-commutative directions becomes finite in the above limit, and its dependence on $θ$ is evaluated explicitly. We also study the dispersion relation. In the weak coupling symmetric phase, it involves a negative IR-singular term, which is responsible for the observed phase transition. In the broken phase, it reveals the existence of the Nambu-Goldstone mode associated with the spontaneous symmetry breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Dynamical aspects of the fuzzy CP$^{2}$ in the large $N$ reduced model with a cubic term

``Fuzzy CP^2'', which is a four-dimensional fuzzy manifold extension of the well-known fuzzy analogous to the fuzzy 2-sphere (S^2), appears as a classical solution in the dimensionally reduced 8d Yang-Mills model with a cubic term involving the structure constant of the SU(3) Lie algebra. Although the fuzzy S^2, which is also a classical solution of the same model, has actually smaller free energy than the fuzzy CP^2, Monte Carlo simulation shows that the fuzzy CP^2 is stable even nonperturbatively due to the suppression of tunneling effects at large N as far as the coefficient of the cubic term ($α$) is sufficiently large. As αis decreased, both the fuzzy CP$^2$ and the fuzzy S^2 collapse to a solid ball and the system is essentially described by the pure Yang-Mills model (α= 0). The corresponding transitions are of first order and the critical points can be understood analytically. The gauge group generated dynamically above the critical point turns out to be of rank one for both CP^2 and S^2 cases. Above the critical point, we also perform perturbative calculations for various quantities to all orders, taking advantage of the one-loop saturation of the effective action in the large-N limit. By extrapolating our Monte Carlo results to N=\infty, we find excellent agreement with the all order results.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Absence of a fuzzy $S^4$ phase in the dimensionally reduced 5d Yang-Mills-Chern-Simons model

We perform nonperturbative studies of the dimensionally reduced 5d Yang-Mills-Chern-Simons model, in which a four-dimensional fuzzy manifold, ``fuzzy S$^{4}$'', is known to exist as a classical solution. Although the action is unbounded from below, Monte Carlo simulations provide an evidence for a well-defined vacuum, which stabilizes at large $N$, when the coefficient of the Chern-Simons term is sufficiently small. The fuzzy S$^{4}$ prepared as an initial configuration decays rapidly into this vacuum in the process of thermalization. Thus we find that the model does not possess a ``fuzzy S$^{4}$ phase'' in contrast to our previous results on the fuzzy S$^{2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Gaussian expansion analysis of a matrix model with the spontaneous breakdown of rotational symmetry

Recently the gaussian expansion method has been applied to investigate the dynamical generation of 4d space-time in the IIB matrix model, which is a conjectured nonperturbative definition of type IIB superstring theory in 10 dimensions. Evidence for such a phenomenon, which is associated with the spontaneous breaking of the SO(10) symmetry down to SO(4), has been obtained up to the 7-th order calculations. Here we apply the same method to a simplified model, which is expected to exhibit an analogous spontaneous symmetry breaking via the same mechanism as conjectured for the IIB matrix model. The results up to the 9-th order demonstrate a clear convergence, which allows us to unambiguously identify the actual symmetry breaking pattern by comparing the free energy of possible vacua and to calculate the extent of ``space-time'' in each direction.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

On Existence of Nontrivial Fixed Points in Large $N$ Gauge Theory in More than Four Dimensions

Inspired by a possible relation between large $N$ gauge theory and string theory, we search for nontrivial fixed points in large $N$ gauge theory in more than four dimensions. We study large $N$ gauge theory through Monte Carlo simulation of the twisted Eguchi-Kawai model in six dimensions as well as in four dimensions. The phase diagram of the system with the two coupling constants which correspond to the standard plaquette action and the adjoint term has been explored.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Two-loop Renormalization in Quantum Gravity near Two Dimensions

We study two--loop renormalization in $(2+ε)$--dimensional quantum gravity. As a first step towards the full calculation, we concentrate on the divergences which are proportional to the number of matter fields. We calculate the $β$ functions and show how the nonlocal divergences as well as the infrared divergences cancel among the diagrams. Although the formalism includes a subtlety concerning the general covariance due to the dynamics of the conformal mode, we find that the renormalization group allows the existence of a fixed point which possesses the general covariance. Our results strongly suggest that we can construct a consistent theory of quantum gravity by the $ε$ expansion around two dimensions.

0 citations0 reviewsNo code linkedOpen access