Keitaro Nagata

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.

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.

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.

We study the entanglement entropy (EE) for pure gauge theories in 1+1 dimensions with the lattice regularization. Using the definition of the EE for lattice gauge theories proposed in a previous paper [1] (S. Aoki, T. Iritani, M. Nozaki, T. Numasawa, N. Shiba and H. Tasaki, JHEP 1506 (2015) 187), we calculate the EE for arbitrary pure as well as mixed states in terms of eigenstates of the transfer matrix in 1+1 dimensional lattice gauge theory. We find that the EE of an arbitrary pure state does not depend on the lattice spacing, thus giving the EE in the continuum limit, and show that the EE for an arbitrary pure state is independent of the real (Minkowski) time evolution. We also explicitly demonstrate the dependence of EE on the gauge fixing at the boundaries between two subspaces, which was pointed out for general cases in the paper [1]. In addition, we calculate the EE at zero as well as finite temperature by the replica method, and show that our result in the continuum limit corresponds to the result obtained before in the continuum theory, with a specific value of the counter term, which is otherwise arbitrary in the continuum calculation. We confirm the gauge dependence of the EE also for the replica method.

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.

We investigate the quantum entanglement entropy for the four-dimensional Euclidean SU(3) gauge theory. We present the first non-perturbative calculation of the entropic |$c$|-function (⁠|$C(l)$|⁠) of SU(3) gauge theory in lattice Monte Carlo simulation using the replica method. For |$0 \leqslant l \leqslant 0.7$| fm, where |$l$| is the length of the subspace, the entropic |$c$|-function is almost constant, indicating conformally invariant dynamics. The value of the constant agrees with that perturbatively obtained from free gluons, with 20% discrepancy. When |$l$| is close to the |$\Lambda _{{\rm QCD}}^{-1} \big (\sim T_c^{-1}\big )$| scale, the entropic |$c$|-function decreases smoothly, and it is consistent with zero within error bars at |$l \gtrsim 0.9$| fm.

We apply the complex Langevin (CL) method to a chiral random matrix theory (ChRMT) at nonzero chemical potential and study the nearest neighbor spacing (NNS) distribution of the Dirac eigenvalues. The NNS distribution is extracted using an unfolding procedure for the Dirac eigenvalues obtained in the CL method. For large quark mass, we find that the NNS distribution obeys the Ginibre ensemble as expected. For small quark mass, the NNS distribution follows the Wigner surmise for the correct convergence case, while it follows the Ginibre ensemble for the wrong convergence case. The Wigner surmise is physically reasonable from the chemical potential independence of the ChRMT. The Ginibre ensemble is known to be favored in a phase-quenched QCD at finite chemical potential. Our result suggests a possibility that the originally universal behavior of the NNS distribution is preserved even in the CL method for the correct convergence case.