掲載種別：研究論文（国際会議プロシーディングス）  

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.

Scopus

researchmap

掲載種別：研究論文（国際会議プロシーディングス）  

We apply the complex Langevin method (CLM) to overcome the sign problem in 4D SU(2) gauge theory with a theta term extending our previous work on the 2D U(1) case. The topology freezing problem can be solved by using open boundary conditions in all spatial directions, and the criterion for justifying the CLM is satisfied even for large θ as far as the lattice spacing is sufficiently small. However, we find that the CP symmetry at θ = π remains to be broken explicitly even in the continuum and infinite-volume limits due to the chosen boundary conditions. In particular, this prevents us from investigating the interesting phase structures suggested by the't Hooft anomaly matching condition. We also try the so-called subvolume method, which turns out to have a similar problem. We therefore discuss a new technique within the CLM, which enables us to circumvent the topology freezing problem without changing the boundary conditions.

Scopus

researchmap

掲載種別：研究論文（学術雑誌）   出版者・発行元：Springer Science and Business Media {LLC}  

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.

DOI： 10.1007/JHEP04(2022)179

Scopus

arXiv

researchmap

掲載種別：研究論文（学術雑誌）   出版者・発行元：Springer Science and Business Media {LLC}  

The tensor renormalization group method is a promising approach to lattice field theories, which is free from the sign problem unlike standard Monte Carlo methods. One of the remaining issues is the application to gauge theories, which is so far limited to U(1) and SU(2) gauge groups. In the case of higher rank, it becomes highly nontrivial to restrict the number of representations in the character expansion to be used in constructing the fundamental tensor. We propose a practical strategy to accomplish this and demonstrate it in 2D U(N) and SU(N) gauge theories, which are exactly solvable. Using this strategy, we obtain the singular-value spectrum of the fundamental tensor, which turns out to have a definite profile in the large-N limit. For the U(N) case, in particular, we show that the large-N behavior of the singular-value spectrum changes qualitatively at the critical coupling of the Gross-Witten-Wadia phase transition. As an interesting consequence, we find a new type of volume independence in the large-N limit of the 2D U(N) gauge theory with the θ term in the strong coupling phase, which goes beyond the Eguchi-Kawai reduction.

DOI： 10.1007/JHEP12(2021)011

Scopus

arXiv

researchmap

掲載種別：研究論文（学術雑誌）   出版者・発行元：Springer Science and Business Media {LLC}  

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.

DOI： 10.1007/JHEP09(2020)023

Scopus

arXiv

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Oxford University Press ({OUP})  

<jats:title>Abstract</jats:title>
<jats:p>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.</jats:p>

DOI： 10.1093/ptep/ptaa042

Scopus

arXiv

researchmap