Source author record

Keitaro Nagata

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

25 published item(s)

preprint2016arXiv

Entanglement entropy for pure gauge theories in 1+1 dimensions using the lattice regularization

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.

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.

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.

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.

preprint2016arXiv

Test for a universal behavior of Dirac eigenvalues in the complex Langevin method

We apply the complex Langevin (CL) method to a chiral random matrix theory (ChRMT) at non-zero 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 correct convergence case, while it follows the Ginibre ensemble for 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 correct convergence case.

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.

preprint2015arXiv

Entanglement in Four-Dimensional SU(3) Gauge Theory

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 Hadronic scale, the entropic $c$-function decreases smoothly, and it is consistent with zero within error bars at $l \gtrsim 0.9$ fm.

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.

preprint2015arXiv

Lee-Yang zero distribution of high temperature QCD and Roberge-Weiss phase transition

Canonical partition functions and Lee-Yang zeros of QCD at finite density and high temperature are studied. Recent lattice simulations have confirmed that the free energy of QCD is a quartic function of quark chemical potential at temperature slightly above pseudo-critical temperature $T_c$, as in the case with a gas of free massless fermions. We present analytic derivation of the canonical partition functions and Lee-Yang zeros for this type of free energy using the saddle point approximation. We also perform lattice QCD simulation in a canonical approach using the fugacity expansion of the fermion determinant, and carefully examine its reliability. By comparing the analytic and numerical results, we conclude that the canonical partition functions follow the Gaussian distribution of the baryon number, and the accumulation of Lee-Yang zeros of these canonical partition functions exhibit the first-order Roberge-Weiss phase transition. We discuss the validity and applicable range of the result and its implications both for theoretical and experimental studies.

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.

preprint2014arXiv

A filtering technique for the temporally reduced matrix of the Wilson fermion determinant

The Wilson fermion determinant can be written in the form of a series expansion in fugacity $ξ=\exp(μ/T)$, provided that the eigenmodes of the temporally reduced operator are obtained. Since the calculation of all eigenmodes rapidly becomes prohibitive for larger volumes, we develop a method to calculate only the low-energy eigenmodes of the reduced matrix using a matrix filetering technique. This provides a basis for an approximation to neglect uninteresting ultraviolet contributions.

preprint2014arXiv

Heavy quark potential at finite imaginary chemical potential

We investigate chemical-potential ($μ$) dependence of the static-quark free energies in both the real and imaginary $μ$ regions, using the clover-improved two-flavor Wilson fermion action and the renormalization-group improved Iwasaki gauge action. Static-quark potentials are evaluated from Polyakov-loop correlators in the deconfinement phase and the imaginary $μ=iμ_{\rm I}$ region and extrapolated to the real $μ$ region with analytic continuation. As the analytic continuation, the potential calculated at imaginary $μ=iμ_{\rm I}$ is expanded into a Taylor-expansion series of $iμ_{\rm I}/T$ up to 4th order and the pure imaginary variable $iμ_{\rm I}/T$ is replaced by the real one $μ_{\rm R}/T$. At real $μ$, the 4th-order term weakens $μ$ dependence of the potential sizably. Also, the color-Debye screening mass is extracted from the color-singlet potential at imaginary $μ$, and the mass is extrapolated to real $μ$ by analytic continuation. The screening mass thus obtained has stronger $μ$ dependence than the prediction of the leading-order thermal perturbation theory at both real and imaginary $μ$.

preprint2013arXiv

Color screening potential at finite density in two-flavor lattice QCD with Wilson fermions

We investigate chemical-potential (μ) dependence of static-quark free energies in both the real and imaginary μregions, performing lattice QCD simulations at imaginary μand extrapolating the results to the real μregion with analytic continuation. Lattice QCD calculations are done on a 16^{3}\times 4 lattice with the clover-improved two-flavor Wilson fermion action and the renormalization-group improved Iwasaki gauge action. Static-quark potential is evaluated from the Polyakov-loop correlation functions in the deconfinement phase. As the analytic continuation, the potential calculated at imaginary μ=iμ_{\rm I} is expanded into a Taylor-expansion series of iμ_{\rm I}/T up to 4th order and the pure imaginary variable iμ_{\rm I}/T is replaced by the real one μ_{\rm R}/T. At real μ, the 4th-order term weakens μdependence of the potential sizably. At long distance, all of the color singlet and non-singlet potentials tend to twice the single-quark free energy, indicating that the interactions between heavy quarks are fully color-screened for finite μ. For both real and imaginary μ, the color-singlet q{\bar q} and the color-antitriplet qq interaction are attractive, whereas the color-octet q{\bar q} and the color-sextet qq interaction are repulsive. The attractive interactions have stronger μ/T dependence than the repulsive interactions. The color-Debye screening mass is extracted from the color-singlet potential at imaginary μ, and the mass is extrapolated to real μby analytic continuation. The screening mass thus obtained has stronger μdependence than the prediction of the leading-order thermal perturbation theory at both real and imaginary μ.

preprint2012arXiv

Development of Lattice QCD Tool Kit on Cell Broadband Engine Processor

We report an implementation of a code for SU(3) matrix multiplication on Cell/B.E., which is a part of our project, Lattice Tool Kit on Cell/B.E.. On QS20, the speed of the matrix multiplication on SPE in single precision is 227GFLOPS and it becomes 20GFLOPS {this vaule was remeasured and corrcted.} together with data transfer from main memory by DNA transfer, which is 4.6% of the hardware peak speed (460GFLOPS), and is 7.4% of the theoretical peak speed of this calculation (268.77GFLOPS). We briefly describe our tuning procedure.

preprint2012arXiv

EoS of finite density QCD with Wilson fermions by Multi-Parameter Reweighting and Taylor expansion

The equation of state (EoS), quark number density and susceptibility at nonzero quark chemical potential $μ$ are studied in lattice QCD simulations with a clover-improved Wilson fermion of 2-flavors and RG-improved gauge action. To access nonzero $μ$, we employ two methods : a multi-parameter reweighting (MPR) in $μ$ and $β$ and Taylor expansion in $μ/T$. The use of a reduction formula for the Wilson fermion determinant enables to study the reweighting factor in MPR explicitly and heigher-order coefficients in Taylor expansion free from errors of noise method, although calculations are limited to small lattice size. As a consequence, we can study the reliability of the thermodynamical quantities through the consistency of the two methods, each of which has different origin of the application limit. The thermodynamical quantities are obtained from simulations on a $8^3\times 4$ lattice with an intermediate quark mass($m_{\rm PS}/m_{\rm V}=0.8)$. The MPR and Taylor expansion are consistent for the EoS and number density up to $μ/T\sim 0.8$ and for the number susceptibility up to $μ/T \sim 0.6$. This implies within a given statistics that the overlap problem for the MPR and truncation error for the Taylor expansion method are negligible in these regions. In order to make MPR methods work, the fluctuation of the reweighting factor should be small. We derive the equation of the reweighting line where the fluctuation is small, and show that the equation of the reweighting line is consistent with the fluctuation minimum condition.

preprint2012arXiv

On early onset of quark number density at zero temperature

We study a longstanding problem in lattice QCD at low temperature and nonzero quark chemical potential on an onset of the quark number density at $μ=m_π/2$. We introduce a physical parametrization of the eigenvalues in the reduction formula of the fermion determinant. It is shown that the parametrization reduces the quark number density operator to an expression with the Fermi distribution of the quark. For each configuration, the eigenvalues of the reduced matrix correspond to one-particle energy states of a quark. The gap of the eigenspectrum of the reduced matrix corresponds to the gap of the energy states, which causes the $μ$-independence of the fermion determinant for small $μ$ at T=0. Once $μ$ exceeds the gap, the quark number density becomes nonzero for each configuration, which causes the early onset of the quark number density.

preprint2012arXiv

Towards extremely dense matter on the lattice

QCD is expected to have a rich phase structure. It is empirically known to be difficult to access low temperature and nonzero chemical potential $μ$ regions in lattice QCD simulations. We address this issue in a lattice QCD with the use of a dimensional reduction formula of the fermion determinant. We investigate spectral properties of a reduced matrix of the reduction formula. Lattice simulations with different lattice sizes show that the eigenvalues of the reduced matrix follow a scaling law for the temporal size $N_t$. The properties of the fermion determinant are examined using the reduction formula. We find that as a consequence of the $N_t$ scaling law, the fermion determinant becomes insensitive to $μ$ as $T$ decreases, and $μ$-independent at T=0 for $μ<m_π/2$. The $N_t$ scaling law provides two types of the low temperature limit of the fermion determinant: (i) for low density and (ii) for high-density. The fermion determinant becomes real and the theory is free from the sign problem in both cases. In case of (ii), QCD approaches to a theory, where quarks interact only in spatial directions, and gluons interact via the ordinary Yang-Mills action. The partition function becomes exactly $Z_3$ invariant even in the presence of dynamical quarks because of the absence of the temporal interaction of quarks. The reduction formula is also applied to the canonical formalism and Lee-Yang zero theorem. We find characteristic temperature dependences of the canonical distribution and of Lee-Yang zero trajectory. Using an assumption on the canonical partition function, we discuss physical meaning of those temperature dependences and show that the change of the canonical distribution and Lee-Yang zero trajectory are related to the existence/absence of $μ$-induced phase transitions.

preprint2011arXiv

Imaginary Chemical Potential Approach for the Pseudo-Critical Line in the QCD Phase Diagram with Clover-Improved Wilson Fermions

The QCD phase diagram is studied in the lattice QCD simulation with the imaginary chemical potential approach. We employ a clover-improved Wilson fermion action of two-flavors and a renormalization-group improved gauge action, and perform the simulation at an intermediate quark mass on a $8^3\times 4$ lattice. The QCD phase diagram in the imaginary chemical potential $μ_I$ region is investigated by performing the simulation for more than 150 points on the $(β,μ_I)$ plane. We find that the Roberge-Weiss phase transition at $μ_I/T=π/3$ is first order and its endpoint is second order, which are identified by the phase of the Polyakov loop. We determine the pseudo-critical line from the susceptibility of the Polyakov loop modulus. We find a clear deviation from a linear dependence of the pseudo-critical line on $μ_I^2$.

preprint2011arXiv

QCD Phase Diagram with Imaginary Chemical Potential

We report our recent results on the QCD phase diagram obtained from the lattice QCD simulation. The location of the phase boundary between hadronic and QGP phases in the two-flavor QCD phase diagram is investigated. The imaginary chemical potential approach is employed, which is based on Monte Carlo simulations of the QCD with imaginary chemical potential and analytic continuation to the real chemical potential region.

preprint2010arXiv

On the two-boson exchange corrections to parity-violating elastic electron-proton scattering

The details of the calculation of the two-boson exchange effects in the parity-violating elastic $ep$ scattering within a simple hadronic model, including both the nucleon and $Δ(1232)$-resonance intermediate states, are presented. We examine the sensitivity of our results with respect to choice of form factors. We emphasize the importance to use correct relations relating N to $Δ$ and $Δ$ to N transition vertex functions. The $NΔ$ Coulomb quadrupole transition is found to play important role at higher $Q^2\geq 3.0$ GeV$^2$. We also elucidate the relation between our results and the well-known result on the $γZE$ effect given by Marciano and Sirlin (MS). The effect of the nucleon contribution $δ_N$ to parity-asymmetry $A_{PV}$, is found to be in general, larger than the corresponding $Δ$ contribution $δ_Δ$ except at extreme forward angles. The corrections to the extracted values of the strange form factors $G^{s}_{E}+βG^{s}_{M}$ from the HAPPEX, A4, and G0 data are also presented. The total TBE corrections to the extracted values of $G^{s}_{E}+βG^{s}_{M}$ in recent experiments of HAPPEX G0, and A4 are, depending on kinematics, found to be small except in a few cases where they range from -20.6% to $48.3 %$.

preprint2010arXiv

Quartet of spin-3/2 baryons in chiral multiplet $(1, 1/2) \oplus (1/2, 1) $ with mirror assignment

We study the possible existence of chiral partners in the spin-$\thalf$ sector of the baryon spectrum. We consider a quartet scheme where four spin-3/2 baryons, $P_{33}$, $D_{33}$, $D_{13}$ and $P_{13}$, group into higher-dimensional chiral multiplets $(1, \half)\oplus (\half,1)$ with a mirror assignment. With an effective $SU(2)_R\times SU(2)_L$ Lagrangian, we derive constraints imposed by chiral symmetry together with the mirror assignment on the masses and coupling constants of the quartet. Using the effective Lagrangian, we try to find a set of baryons suitable for the chiral quartet. It turns out that two cases reasonably agree with the mass pattern of the quartet: ($Δ(1600)$, $Δ(1940)$, $N(1520)$, $N(1720)$) and ($Δ(1920)$, $Δ(1940)$, $N(2080)$, $N(1900)$).

preprint2010arXiv

Wilson Fermion Determinant in Lattice QCD

We present a formula for reducing the rank of Wilson fermions from $4 N_c N_x N_y N_z N_t$ to $4 N_c N_x N_y N_z$ keeping the value of its determinant. We analyse eigenvalues of a reduced matrix and coefficients $C_n$ in the fugacity expansion of the fermion determinant $\sum_n C_n (\exp(μ/T))^n$, which play an important role in the canonical formulation, using lattice QCD configurations on a $4^4$ lattice. Numerically, $\log |C_n|$ varies as $N_x N_y N_z$, and goes easily over the standard numerical range; We give a simple cure for that. The phase of $C_n$ correlates with the distribution of the Polyakov loop in the complex plain. These results lay the groundwork for future finite density calculations in lattice QCD.

preprint2009arXiv

Contribution of Two-Boson-Exchange with $Δ(1232)$ Excitation to Parity-Violating Elastic Electron-Proton Scattering

We study the leading electroweak corrections in the precision measurement of the strange form factors. Specifically, we calculate the two-boson-exchange (TBE), two-photon-exchange (TPE) plus $γZ$-exchange ($γZ$E), corrections with $Δ(1232)$ excitation to the parity-violating asymmetry of the elastic electron-proton scattering. The interplay between nucleon and $Δ$ contributions is found to depend strongly on the kinematics, as $δ_Δ$ begins as negligible at backward angles but becomes very large and negative and dominant at forward angles, while $δ_N$ always stays positive and decreases monotonically with increasing $ε$. The total TBE corrections to the extracted values of $G^{s}_{E}+βG^{s}_{M}$ in recent experiments of HAPPEX and G0 are, depending on kinematics, found to be large and range between 13% to -75% but small in the case of A4 experiments.

preprint2009arXiv

Wilson fermions with imaginary chemical potential

We study the phase structure of imaginary chemical potential. We calculate the Polyakov loop using clover-improved Wilson action and renormalization improved gauge action. We obtain a two-state signals indicating the first order phase transition for $β= 1.9, μ_I = 0.2618, κ=0.1388$ on $8^3\times 4$ lattice volume We also present a result of the matrix reduction formula for the Wilson fermion.