General approach to the sign problem: 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 observable. 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.

title: Towards an Effective Importance Sampling in Monte Carlo Simulations of a System with a Complex Action
abstract: 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.

title: Matrix model with manifest general coordinate invariance
abstract: We present a formulation of a matrix model which manifestly possesses the general coordinate invariance when we identify the large N matrices with differential operators. In order to build a matrix model which has the local Lorentz invariance, we investigate how the so (9,1) Lorentz symmetry and the u ( N ) gauge symmetry are mixed together. We first analyze the bosonic part of the model, and we find that the Einstein gravity is reproduced in the classical low-energy limit. We then present a proposal to build a matrix model which has N =2 SUSY and reduces to the type IIB supergravity in the classical low-energy limit.

title: Monte Carlo simulations of a supersymmetric matrix model of dynamical compactification in non perturbative string theory
abstract: 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.

title: New states of gauge theories on a circle
abstract: A bstractWe study a one-dimensional large-N U(N) gauge theory on a circle as a toy model of higher dimensional Yang-Mills theories at finite temperature. To investigate the profile of the thermodynamical potential in this model, we evaluate a stochastic time evolution of several states, and find that an unstable confinement phase at high temperature does not decay to a stable deconfinement phase directly. Before it reaches the deconfinement phase, it develops to several intermediate states. These states are characterised by the expectation values of the Polyakov loop operators, which wind the temporal circle different times. We reveal that these intermediate states are the saddle point solutions of the theory, and similar solutions exist in a wide class of SU(N) and U(N) gauge theories on S1 including QCD and pure Yang-Mills theories in various dimensions. We also consider a Kaluza-Klein gravity, which is the gravity dual of the one-dimensional gauge theory on a spatial S1, and show that these solutions may be related to multi black holes localised on the S1. Then we present a connection between the stochastic time evolution of the gauge theory and the dynamical decay process of a black string though the Gregory-Laflamme instability.

title: Monte Carlo studies of the spontaneous rotational symmetry breaking in dimensionally reduced super Yang-Mills models
abstract: A bstractIt 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.

title: Propagators for a Quantized Scalar Field in Some Isotropic Universe
abstract: To make clearer the relation examined in a previous paper between Chitre and Hartle's path-integral quantization procedure and our canonical one for a massive scalar field in a special isotropic universe, various propagators such as the 4-dimensional commutation func­ tion, the elementary solution and the associated Feynman propagator are explicitly deter­ mined. Owing to the determination, there arises a new changing procedure superior to the previous one from the isotropic expanding universe to the Minkowski space-time. Thus it is shown without any ambiguity that both the path-integral quantization procedure and our canonical one leading to the same pair-creation o£ particles as Chitre and Hartle's cannot recover the usual theory for a free field in the limit when the cosmic expansion disappears. (Notice, however, that our quantization procedure permits another case which can recover the Minkowski theory for a free field in the above limit, while it gives rise to no pair-creation.) Propagators for a massless scalar field are also dealt with in connection with the problem of conformal anomaly arising in the vacuum expectation value of the stress-energy tensor for our quantized scalar field, which is common to any quantized field in a curved space-time.

title: Exact Disk Solution of the Einstein Equation
abstract: The 4-soliton solution with complex pole trajectories to the vacuum Einstein equation is studied. The solution describes a static and axisymmetric gravitational field with asymptotic flatness, which is produced by a thin disk. We obtain the energy-momentum tensor of the disk to reveal its properties. The soliton technique 1 ) to solve the gravitational field equation gave us a new viewpoint on the exact solutions of the vacuum Einstein equation. There is an infinite series of solutions which are characterized by the soliton number n. In the stationary case with axial symmetry the series includes the Kerr-NUT solution as a 2-soliton solution and the Tomimatsu-Sato solution as a certain limit of 4-soliton solution. In the static case 2 ),3) there is an infinite series of the Weyl solutions, which is interpreted as a nonlinear superposition of multi-Weyl solutions. Two of the present authors 3 ),4) showed that the soliton technique gives not only these series of solutions but a series of ring-type solutions by adopting complex pole trajectories in the static case. The first non-trivial example of this series is the Weyl solution expressed in the oblate spheroidal coordinates,5) which has a ring-type curvature singularity. In our classification this solution is interpreted as the 2-soliton solution with non-complex-conjugate pole trajectories. The next example is the 4-soliton solution in which there are two types. One is just a nonlinear superposition of the above 2-soliton solutions, which has two ring-type singularities with different radii. The other, which is more important, is the 4-soliton solution constructed from two pairs of complex conjugate pole trajectories. This was referred to as the two­ ring solution, because the curvature invariants have finite limits at two rings although the metric coefficients diverge there. This fact contrasts with the fact that the curvature invariants diverge at the ring in the 2-soliton case. The purpose of this paper is to study the detailed structure of the two-ring solution. It will be shown that the solution should be called a disk solution rather than a ring solution. The properties of the disk will be also clarified by calculating the energy-momentum tensor. In the next section we give the solution in an explicit form and investigate its mathematical structure. From this we find that the plane containing the rings is separated into a regular region inside the inner ring and a region with the discontinuity of curvature outside the inner ring. We shall call this outer region a "disk" because of its non-zero energy density as will be.discussed in § 3. The disk is further divided into two parts by the outer ring. In § 3 we obtain the

title: On Anisotropic and Homogeneous Cosmological Models in the Renormalized Theory of Gravitation

title: Equilibrium Condition of Dilaton Black Holes with Electric Charge

title: De Sitter metric as a self-consistent solution of the back reaction problem
abstract: The de Sitter metric is a solution of the Einstein equation even when quantum effects of matter fields are included. However between the curvature and the cosmological constant Λ changes drastically. For some values of Λ, there are two or three solutions, while there is none in some other cases. There can be even a positive curvature solution for negative Λ. These features depend on mass and spin of matter fields and the matter-curvature coupling.