Source author record

Shinji Ejiri

Shinji Ejiri 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

30works
5topics
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

30 published item(s)

preprint2023arXiv

Critical point in heavy-quark region of QCD on fine lattices

We perform a finite-size scaling analysis of the critical point in the heavy-quark region of QCD at nonzero temperature. Our previous analysis on the Binder cumulant at $N_t=4$ is extended to finer lattices with $N_t=6$ and $8$. The aspect ratio is also extended up to $15$ to suppress the non-singular contribution. High-precision analysis of the Binder cumulant is realized by an efficient Monte-Carlo simulation with the hopping-parameter expansion (HPE). Effects of higher-order terms in the HPE are incorporated by the reweighting method.

preprint2022arXiv

Canonical partition function and center symmetry breaking in finite density lattice gauge theories

We study the nature of the phase transition of lattice gauge theories at high temperature and high density by focusing on the probability distribution function, which represents the probability that a certain density will be realized in a heat bath. The probability distribution function is obtained by creating a canonical partition function fixing the number of particles from the grand partition function. However, if the Z_3 center symmetry, which is important for understanding the finite temperature phase transition of SU(3) lattice gauge theory, is maintained on a finite lattice, the probability distribution function is always zero, except when the number of particles is a multiple of 3. For U(1) lattice gauge theory, this problem is more serious. The probability distribution becomes zero when the particle number is nonzero. This problem is essentially the same as the problem that the expectation value of the Polyakov loop is always zero when calculating with finite volume. In this study, we propose a solution to this problem. We also propose a method to avoid the sign problem, which is an important problem at finite density, using the center symmetry. In the case of U(1) lattice gauge theory with heavy fermions, numerical simulations are actually performed, and we demonstrate that the probability distribution function at a finite density can be calculated by the method proposed in this study. Furthermore, the application of this method to QCD is discussed.

preprint2022arXiv

Scope and convergence of the hopping parameter expansion in finite temperature QCD with heavy quarks around the critical point

Hopping parameter expansion is a useful tool to investigate heavy dynamical quarks in lattice QCD, while the range of its applicability has been sometimes questioned. We study the convergence and the valid range of the hopping parameter expansion in the determination of the critical point (critical quark mass) of QCD with heavy quarks at finite temperature and density. On lattices with sufficiently large spatial extent, the terms in the hopping parameter expansion are classified into Wilson loop terms and Polyakov-type loop terms. We first study the case of the worst convergence in which all the gauge link variables are unit matrices and thus the Wilson loops and the Polyakov-type loops get their maximum values. We perform explicit calculation up to more than 100th order of the hopping parameter expansion. We show that the hopping parameter expansion is convergent up to the chiral limit of free Wilson quarks. We then perform a Monte-Carlo simulation to measure correlation among Polyakov-type loop terms up to the 20th order of the hopping parameter expansion. In previous studies, strong correlation between the leading order Polyakov loop term and the next-to-leading order bent Polyakov loop terms was reported and used to construct an effective theory to incorporate the next-to-leading order effect by a shift of the leading order coupling parameter. We establish that the strong correlation among Polyakov-type loop terms holds also at higher orders of the hopping parameter expansion, and extend the effective theory to incorporate higher-order effects up to high orders. Using the effective theory, we study the truncation error of the hopping parameter expansion. We find that the previous next-to-leading order result of the critical point for $N_t=4$ are well reliable. For $N_t \ge 6$, we need to incorporate higher-order effects in the effective theory.

preprint2021arXiv

Finite-size scaling around the critical point in the heavy quark region of QCD

Finite-size scaling is investigated in detail around the critical point in the heavy-quark region of nonzero temperature QCD. Numerical simulations are performed with large spatial volumes up to the aspect ratio $N_s/N_t=12$ at a fixed lattice spacing with $N_t=4$. We show that the Binder cumulant and the distribution function of the Polyakov loop follow the finite-size scaling in the $Z(2)$ universality class for large spatial volumes with $N_s/N_t \ge 9$, while, for $N_s/N_t \le 8$, the Binder cumulant becomes inconsistent with the $Z(2)$ scaling. To realize the large-volume simulations in the heavy-quark region, we adopt the hopping parameter expansion for the quark determinant: We generate gauge configurations using the leading order action including the Polyakov loop term for $N_t=4$, and incorporate the next-to-leading order effects in the measurements by the multipoint reweighting method. We find that the use of the leading-order configurations is crucially effective in suppressing the overlapping problem in the reweighting and thus reducing the statistical errors.

preprint2021arXiv

Latent heat and pressure gap at the first-order deconfining phase transition of SU(3) Yang-Mills theory using the small flow-time expansion method

We study latent heat and the pressure gap between the hot and cold phases at the first-order deconfining phase transition temperature of the SU(3) Yang-Mills theory. Performing simulations on lattices with various spatial volumes and lattice spacings, we calculate the gaps of the energy density and pressure using the small flow-time expansion (SFtX) method. We find that the latent heat $Δε$ in the continuum limit is $Δε/T^4 = 1.117 \pm 0.040$ for the aspect ratio $N_s/N_t=8$ and $1.349 \pm 0.038$ for $N_s/N_t=6$ at the transition temperature $T=T_c$. We also confirm that the pressure gap is consistent with zero, as expected from the dynamical balance of two phases at $T_c$. From hysteresis curves of the energy density near $T_c$, we show that the energy density in the (metastable) deconfined phase is sensitive to the spatial volume, while that in the confined phase is insensitive. Furthermore, we examine the effect of alternative procedures in the SFtX method - the order of the continuum and the vanishing flow-time extrapolations, and also the renormalization scale and higher-order corrections in the matching coefficients. We confirm that the final results are all very consistent with each other for these alternatives.

preprint2020arXiv

Calculation of PCAC mass with Wilson fermion using gradient flow

We calculate the PCAC mass for $(2+1)$ flavor full QCD with Wilson-type quarks. We adopt the Small Flow-time eXpansion (SFtX) method based on the gradient flow which provides us a general way to compute correctly renormalized observables even if the relevant symmetries for the observable are broken explicitly due to the lattice regularization, such as the Poincáre and chiral symmetries. Our calculation is performed on heavy $u, d$ quarks mass ($m_π/m_ρ\simeq0.63$) and approximately physical $s$ quark mass with fine lattice $a \simeq 0.07$~fm. The results are compared with those computed with the Schrödinger functional method.

preprint2020arXiv

End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD

We study the end point of the first-order deconfinement phase transition in two and 2+1 flavor QCD in the heavy quark region of the quark mass parameter space. We determine the location of critical point at which the first-order deconfinement phase transition changes to crossover, and calculate the pseudo-scalar meson mass at the critical point. Performing quenched QCD simulations on lattices with the temporal extents Nt=6 and 8, the effects of heavy quarks are determined using the reweighting method. We adopt the hopping parameter expansion to evaluate the quark determinants in the reweighting factor. We estimate the truncation error of the hopping parameter expansion by comparing the results of leading and next-to-leading order calculations, and study the lattice spacing dependence as well as the spatial volume dependence of the result for the critical point. The overlap problem of the reweighting method is also examined. Our results for Nt=4 and 6 suggest that the critical quark mass decreases as the lattice spacing decreases and increases as the spatial volume increases.

preprint2016arXiv

Critical point search from an extended parameter space of lattice QCD at finite temperature and density

Aiming to understand the phase structure of lattice QCD at nonzero temperature and density, we study the phase transitions of QCD in an extended parameter space, where the number of flavor and quark masses are considered as parameters. Performing simulations of 2 flavor QCD and using the reweighting method, we investigate (2+Nf) flavor QCD at finite density, where two light flavors and Nf massive flavors exist. Calculating probability distribution functions, we determine the critical surface terminating first order phase transitions in the parameter space of the light quark mass, the heavy quark mass and the chemical potential. Through the study of the many flavor system, we discuss the phase structure of QCD at finite density.

preprint2016arXiv

Equation of state in (2+1)-flavor QCD with gradient flow

The energy-momentum tensor and equation of state are studied in finite-temperature (2+1)-flavor QCD with improved Wilson quarks using the method proposed by Makino and Suzuki based on the gradient flow. We find that the results of the gradient flow are consistent with the previous results using the $T$-integration method at $T \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 280$ MeV ($N_t\mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}}10$), while a disagreement is found at $T \mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}} 350$ MeV ($N_t \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 8$) presumably due to the small-$N_t$ lattice artifact. We also report on the results on the renormalized chiral condensate and its disconnected susceptibility using the method of Hieda and Suzuki. The results show a clear signal of the expected chiral restoration crossover even with Wilson-type quarks which violate the chiral symmetry explicitly.

preprint2016arXiv

Exploring the nature of chiral phase transition in two-flavor QCD using extra heavy quarks

Chiral phase transition of two flavor QCD at finite quark masses is known to be crossover except near the chiral limit, but it can turn to a first order transition when adding many extra flavors. This property is used to explore the nature of the phase transition of massless two flavor QCD using lattice numerical simulations. The extra heavy flavors being incorporated in the form of the hopping parameter expansion through the reweighting, the number of the extra flavors and their masses appear only in a single parameter, defined by $h$. We determine the critical value of the parameter, at which the first order and the crossover regions are separated, and examine its dependence on the two flavor mass. The lattice calculations are carried out at $N_t$=4, and show that the critical value does not depend on the two flavor mass in the range we have studied ($0.46 \le m_π/m_ρ\le 0.66$) and appears to remain finite and positive in the chiral limit, suggesting that the phase transition of massless two flavor QCD is of second order.

preprint2016arXiv

Latent heat at the first order phase transition point of SU(3) gauge theory

We calculate the energy gap (latent heat) and pressure gap between the hot and cold phases of the SU(3) gauge theory at the first order deconfining phase transition point. We perform simulations around the phase transition point with the lattice size in the temporal direction Nt=6, 8 and 12 and extrapolate the results to the continuum limit. We also investigate the spatial volume dependence. The energy density and pressure are evaluated by the derivative method with non-perturabative anisotropy coefficients. We adopt a multi-point reweighting method to determine the anisotropy coefficients. We confirm that the anisotropy coefficients approach the perturbative values as Nt increases. We find that the pressure gap vanishes at all values of Nt when the non-perturbative anisotropy coefficients are used. The spatial volume dependence in the latent heat is found to be small on large lattices. Performing extrapolation to the continuum limit, we obtain $ Δε/T^4 = 0.75 \pm 0.17 $ and $ Δ(ε-3 p)/T^4 = 0.623 \pm 0.056.$

preprint2016arXiv

Many flavor approach to study the nature of chiral phase transition of two-flavor QCD

We perform lattice numerical simulations to study the phase transition of QCD at finite temperature to clarify the nature of the transition of massless two flavor QCD. We investigate QCD with two light and Nf heavy quarks instead of two-flavor QCD, and focus on the light quark mass dependence of the critical heavy mass, below which the transition is of first order. The heavy quarks are incorporated into two flavor configurations in the form of the hopping parameter expansion through the reweighting technique. The nature of the transition is identified by the shape of the constraint effective potential at the critical temperature. Our result indicates that the critical heavy mass remains finite in the chiral limit of the two flavors, suggesting the phase transition of massless two-flavor QCD is of second order.

preprint2016arXiv

Temperature dependence of topological susceptibility using gradient flow

We study temperature dependence of the topological susceptibility with the $N_{f}=2+1$ flavors Wilson fermion. We have two major interests in this paper. One is a comparison of gluonic and fermionic definitions of the topological susceptibility. Two definitions are related by the chiral Ward-Takahashi identity but their coincidence is highly non-trivial for the Wilson fermion. By applying the gradient flow both for the gauge and quark fields we find a good agreement of these two measurements. The other is a verification of a prediction of the dilute instanton gas approximation at low temperature region $T_{pc}< T<1.5T_{pc}$, for which we confirm the prediction that the topological susceptibility decays with power $χ_{t}\propto(T/T_{pc})^{-8}$ for three flavors QCD.

preprint2015arXiv

Analytic continuation of finite density QCD with heavy quarks in the strong coupling region

Complex nature of finite density QCD with heavy quarks in the strong coupling region is studied. For this purpose, we consider the effective potential as a function of Polyakov line, and study thermodynamic singularities and associated Stokes boundaries in the complex chemical potential plane. We also perform an explicit analytic continuation of the first order transition and crossover lines in the complex chemical potential plane.

preprint2015arXiv

First order transition regions in the quark masses and chemical potential parameter space of QCD

We investigate the phase transitions of (2+Nf)-flavor QCD, where two light flavors and Nf massive flavors exist, aiming to understand the phase structure of (2+1)-flavor QCD. Performing simulations of 2-flavor QCD with improved staggered and Wilson fermions and using the reweighting method, we calculate probability distribution functions in the many-flavor QCD. Through the shape of distribution functions, we determine the critical surface terminating first order phase transitions in the parameter space of the light quark mass, heavy quark mass and the chemical potential, and find that the first order region becomes larger with Nf. We then study the critical surface at finite density for large Nf and the first order region is found to become wider with the increasing chemical potential. On the other hand, the light quark mass dependence of the critical mass of heavy quarks seems weak in the region we investigated. The result of this weak dependence suggests that the critical mass of heavy quark remains finite in the chiral limit of 2-flavors and there exists a second order transition region on the line of the 2-flavor massless limit above the tri-critical point. Moreover, we extend the study of 2-flavor QCD at finite density to the case of a complex chemical potential and investigate the singularities where the partition function vanishes, so-called Lee-Yang zeros. The plaquette effective potential is computed in the complex plane. We find that the shape of the effective potential changes from single-well on the real axis to double-well at large imaginary chemical potential and the double-well potential causes the singularities.

preprint2015arXiv

Many flavor approach to study the critical point in finite density QCD

We discuss the QCD critical point at finite density through the study of many flavor QCD, in which two light flavors and Nf massive flavors exist. Performing simulations of QCD with two flavors of improved Wilson fermions, we calculate probability distribution functions of many flavor QCD at finite temperature and density. The dynamical effects of massive flavors and the chemical potential are added using the reweighting technique. From the shape of the distribution functions, we determine the critical surface separating the first order transition and crossover regions in the parameter space of the light and massive quark masses and the chemical potentials. It is found that the critical massive quark mass becomes larger as the chemical potential increases in (2+Nf) flavor QCD. The indication to the (2+1) flavor QCD is then discussed.

preprint2014arXiv

Complex singularities around the QCD critical point at finite densities

Partition function zeros provide alternative approach to study phase structure of finite density QCD. The structure of the Lee-Yang edge singularities associated with the zeros in the complex chemical potential plane has a strong influence on the real axis of the chemical potential. In order to investigate what the singularities are like in a concrete form, we resort to an effective theory based on a mean field approach in the vicinity of the critical point. The crossover is identified as a real part of the singular point. We consider the complex effective potential and explicitly study the behavior of its extrema in the complex order parameter plane in order to see how the Stokes lines are associated with the singularity. Susceptibilities in the complex plane are also discussed.

preprint2013arXiv

End Point of a First-Order Phase Transition in Many-Flavor Lattice QCD at Finite Temperature and Density

Towards the feasibility study of the electroweak baryogenesis in realistic technicolor scenario, we investigate the phase structure of (2+Nf)-flavor QCD, where the mass of two flavors is fixed to a small value and the others are heavy. For the baryogenesis, an appearance of a first order phase transition at finite temperature is a necessary condition. Using a set of configurations of two-flavor lattice QCD and applying the reweighting method, the effective potential defined by the probability distribution function of the plaquette is calculated in the presence of additional many heavy flavors. Through the shape of the effective potential, we determine the critical mass of heavy flavors separating the first order and crossover regions and find it to become larger with Nf. We moreover study the critical line at finite density and the first order region is found to become wider as increasing the chemical potential. Possible applications to real (2+1)-flavor QCD are discussed.

preprint2013arXiv

Nature of finite temperature and density phase transitions in many-flavor QCD

We investigate the phase structure of (2+Nf)-flavor QCD, where two light flavors and Nf massive flavors exist, to discuss the feasibility of the electroweak baryogenesis in realistic technicolor scenario and to understand properties of finite density QCD. Because an appearance of a first order phase transition at finite temperature is a necessary condition for the baryogenesis, it is important to study the nature of finite temperature phase transition. Applying the reweighting method, the probability distribution function of the plaquette is calculated in the many-flavor QCD. Through the shape of the distribution function, we determine the critical mass of heavy flavors terminating the first order region, and find it to become larger with Nf. We moreover study the critical line at finite density and the first order region is found to become wider as increasing the chemical potential. We discuss how the properties of real (2+1)-flavor QCD at finite temperature and density can be extracted from simulations of many-flavor QCD.

preprint2013arXiv

Phase structure of hot dense QCD by a histogram method

We study the phase structure of QCD at high temperature and density by lattice QCD simulations adopting a histogram method. The quark mass dependence and the chemical potential dependence of the nature of phase transition are investigated focusing on the probability distribution function (histogram). The shape of the distribution function changes with the quark mass and chemical potential. Through the shape of the distribution, the critical surface which separates the first order transition and crossover regions in the heavy quark region is determined for the (2+1)-flavor case. Moreover, we determined the critical point at finite density for two-flavor QCD with an intermediate quark mass, using a Gaussian approximation of the complex phase distribution of the quark determinant. The chemical potential dependence of the critical quark mass is also evaluated in the situation where two light quarks and many massive quarks exist. We find that the first order transition region becomes wider with the chemical potential in the many-flavor QCD.

preprint2013arXiv

Singularities around the QCD critical point in the complex chemical potential plane

We consider thermodynamic singularities appearing in the complex chemical potential plane in the vicinity of QCD critical point. In order to investigate what the singularities are like in a concrete form, we resort to an effective theory based on a mean field approach. We study the behavior of extrema of the real part of the complex effective potential in the complex order parameter plane.

preprint2012arXiv

Ab initio study of QCD thermodynamics on the lattice at zero and finite densities

The WHOT-QCD Collaboration is pushing forward a series of lattice studies of QCD at finite temperatures and densities using improved Wilson quarks. Because Wilson-type quarks require more computational resources than the more widely adopted staggered-type quarks, various theoretical and computational techniques have to be developed and applied. In this paper, we introduce the fixed-scale approach armed with the T-integration method, the Gaussian method based on the cumulant expansion, and the histogram method combined with the reweighting technique. Adopting these methods, we have carried out the first study of finite-density QCD with Wilson-type quarks and the first calculation of the equation of state with 2+1 flavors of Wilson-type quarks. We present results of these studies and discuss perspectives towards a clarification of the properties of 2+1 flavor QCD at the physical point, at finite temperatures and densities.

preprint2012arXiv

Phase structure of finite density QCD with a histogram method

We study the phase structure of QCD in the $T-μ$ plane using a histogram method and the reweighting technique by performing phase quenched simulations of two-flavor QCD with RG-improved gauge action and O($a$) improved Wilson quark action. Taking the effects of the complex phase of the quark determinant using the cumulant expansion method, we calculate the probability distribution function of plaquette and phase-quenched determinant as a function of $T$ and $μ$. We discuss the order of the QCD phase transition consulting the shape of the probability distribution function.

preprint2010arXiv

Lattice QCD at finite temperature and density

We study the phase structure of QCD at finite temperature and density by numerical simulations on a lattice. The most important point for the numerical study at finite density is treatment of the sign problem. We propose a method to avoid the sign problem, which is based on a cumulant expansion of the complex phase in the density of state method combined with the reweighting method. Using the method, we study the critical point terminating a first order phase transition line in lattice QCD at high temperature and density.

preprint2009arXiv

Singularities of QCD in the complex chemical potential plane

We study the thermodynamic singularities of QCD in the complex chemical potential plane by a numerical simulation of lattice QCD, and discuss a method to understand the nature of the QCD phase transition at finite density from the information of the singularities. The existence of singular points at which the partition function (Z) vanishes is expected in the complex plane. These are called Lee-Yang zeros or Fisher zeros. We investigate the distribution of these singular points using the data obtained by a simulation of two-flavor QCD with p4-improved staggered quarks. The convergence radius of a Taylor expansion of ln Z in terms of the chemical potential is also discussed.

preprint2009arXiv

Third moments of conserved charges as probes of QCD phase structure

The third moments of conserved charges, the baryon and electric charge numbers, and energy, as well as their mixed moments, carry more information on the state around the QCD phase boundary than previously proposed fluctuation observables and higher order moments. In particular, their signs give plenty of information on the location of the state created in relativistic heavy ion collisions in the temperature and baryon chemical potential plane. We demonstrate this with an effective model.

preprint2009arXiv

Third moments of conserved charges in QCD phase diagram

We point out that the third moments of conserved charges, the baryon and electric charge numbers, and energy, as well as their mixed moments, change their signs around the QCD phase boundary in the temperature and baryon chemical potential plane. These signs can be measured in relativistic heavy ion collisions, and will give clear information on the phase structure of QCD and the state of the system in the early stage of relativistic heavy ion collisions. The behaviors of these moments on the temperature axis and at small quark chemical potential can be analyzed in lattice QCD simulations. We emphasize that the third moments obtained on the lattice, together with the experimental results, will provide a deep understanding about the QCD phase diagram and the location of the state created in heavy ion collisions.

preprint2006arXiv

Lee-Yang zero analysis for the study of QCD phase structure

We comment on the Lee-Yang zero analysis for the study of the phase structure of QCD at high temperature and baryon number density by Monte-Carlo simulations. We find that the sign problem for non-zero density QCD induces a serious problem in the finite volume scaling analysis of the Lee-Yang zeros for the investigation of the order of the phase transition. If the sign problem occurs at large volume, the Lee-Yang zeros will always approach the real axis of the complex parameter plane in the thermodynamic limit. This implies that a scaling behavior which would suggest a crossover transition will not be obtained. To clarify this problem, we discuss the Lee-Yang zero analysis for SU(3) pure gauge theory as a simple example without the sign problem, and then consider the case of non-zero density QCD. It is suggested that the distribution of the Lee-Yang zeros in the complex parameter space obtained by each simulation could be more important information for the investigation of the critical endpoint in the $(T, μ_q)$ plane than the finite volume scaling behavior.

preprint2000arXiv

A Study on the Non-perturbative Existence of Yang-Mills Theories with Large Extra Dimensions

Pure lattice SU(2) Yang-Mills theory in five dimensions is considered, where an extra dimension is compactified on a circle. Monte-Carlo simulations indicate that the theory possesses a continuum limit with a non-vanishing string tension if the compactification radius is smaller than a certain value $R_M$ which is $O(1/10)$ of the inverse of the square root of the string tension. We verify non-perturbatively the power-law running of gauge coupling constant. Our method can be applied to the investigation of continuum limits in other higher-dimensional gauge theories.

preprint1996arXiv

Monopole action and monopole condensation in SU(3) lattice QCD

Effective monopole actions for various extended monopoles are derived from vacuum configurations after abelian projection in the maximally abelian gauge in $T=0$ and $T\ne 0$ $SU(3)$ lattice QCD. The actions obtained appear to be independent of the lattice volume adopted. At zero temperature, monopole condensation is seen to occur from energy-entropy balance in the strong coupling region. Larger $β$ is included in the monopole condensed phase as more extended monopoles are considered. The scaling seen in the $SU(2)$ case is not yet observed. The renormalization flow diagram suggests the existence of an infrared fixed point. A hysteresis behavior is seen around the critical temperature in the case of the $T\ne 0$ action.