Source author record

Yoshinobu Kuramashi

Yoshinobu Kuramashi 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

37works
8topics
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

37 published item(s)

preprint2022arXiv

Metal-insulator transition in (2+1)-dimensional Hubbard model with tensor renormalization group

We investigate the doping-driven metal-insulator transition of the (2+1)-dimensional Hubbard model in the path-integral formalism with the tensor renormalization group method. We calculate the electron density $\langle n\rangle$ as a function of the chemical potential $μ$ choosing three values of the Coulomb potential with $U=80$, $8$, and $2$ as representative cases of the strong, intermediate, and weak couplings. We have determined the critical chemical potential at each $U$, where the Hubbard model undergoes the metal-insulator transition from the half-filling plateau with $\langle n\rangle=1$ to the metallic state with $\langle n\rangle > 1$. Our results indicate that the model exhibits the metal-insulator transition over the vast region of the finite coupling $U$.

preprint2022arXiv

Nucleon isovector tensor charge from lattice QCD with physical light quarks

We present preliminary results for the axial, scalar and tensor charges of the nucleon measured in 2+1 flavor QCD with the physical light quarks ($m_π=135$ MeV). Our simulations are carried out with gauge configurations generated by the PACS Collaboration with the stout-smeared $O(a)$ improved Wilson fermions and Iwasaki gauge action at a single lattice spacing of $0.085\ (\mathrm{fm})$. There are two lattice ensembles of the PACS gauge configurations, which have physical lattice sizes over $(10\ \mathrm{fm})^4$ and $(5\ \mathrm{fm})^4$, respectively. We compute the nucleon three-point correlation functions in the axial, scalar, and tensor channels. For the renormalization, we use the Rome-Southampton method as the intermediate scheme in order to evaluate the renormalization constants for the scalar and tensor currents in fully nonperturbative manner. We then evaluate the renormalized values of the scalar and tensor charges ($g_S$ and $g_T$) in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV with a help of the continuum perturbation theory for the matching between two schemes. We compare our preliminary results of $g_S$ and $g_T$ with those of other collaboration results.

preprint2022arXiv

Tensor renormalization group study of (3+1)-dimensional $Z_2$ gauge-Higgs model at finite density

We investigate the critical endpoints of the (3+1)-dimensional $Z_2$ gauge-Higgs model at finite density together with the (2+1)-dimensional one at zero density as a benchmark using the tensor renormalization group method. We focus on the phase transition between the Higgs phase and the confinement phase at finite chemical potential along the critical end line. In the (2+1)-dimensional model, the resulting endpoint is consistent with a recent numerical estimate by the Monte Carlo simulation. In the (3+1)-dimensional case, however, the location of the critical endpoint shows disagreement with the known estimates by the mean-field approximation and the Monte Carlo studies. This is the first application of the tensor renormalization group method to a four-dimensional lattice gauge theory and a key stepping stone toward the future investigation of the phase structure of the finite density QCD.

preprint2021arXiv

The lower moments of nucleon structure functions in lattice QCD with physical quark masses

We present results for the nucleon structure functions and form factors obtained from 2+1 flavor lattice QCD with physical light quark masses ($m_π=135$ MeV) in a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared ${\mathscr{O}}(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β=1.82$ and $2.00$ corresponding to lattice spacings of $0.085$ fm and $0.064$ fm respectively. The lower moments of structure functions, $\langle x \rangle_{u-d}$ and $\langle x \rangle_{Δu - Δd}$ given by the twist-2 operators being properly renormalized, are evaluated in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV only at $β=1.82$, since the renormalization factors at $β=2.00$ have not yet determined nonperturbatively in the RI/MOM scheme. Instead, at two lattice spacings, we evaluate appropriate ratios of $g_{A}/g_{V}$ and $\langle x \rangle_{u-d}/\langle x \rangle_{Δu -Δd}$, which are not renormalized in the continuum limit. These quantities thus can be directly compared with the experimental data without the renormalization.

preprint2020arXiv

$K_{l3}$ form factors at the physical point on (10.9 fm)$^3$ volume

We present the calculation of the $K_{l3}$ form factors with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on a large volume of (10.9 fm)$^3$ at one lattice spacing of $a = 0.085$ fm. We extract the form factors from 3-point functions with three different time separations between the source and sink operators to confirm suppression of excited state contributions. The form factors are calculated in very close to the zero momentum transfer, $q^2 = 0$, thanks to the large volume, so that stable interpolations to $q^2 = 0$ are carried out. Using our form factors, we obtain the form factor at $q^2 = 0$, $f_+(0) = 0.9603(16)(^{+14}_{\ -4})(44)(19)(1)$, where the first, second, and fifth errors are statistical, systematic errors from fit functions and the isospin breaking effect, respectively. The third and fourth errors denote the finite lattice spacing effects estimated from the renormalization factor and contribution beyond the leading order SU(3) chiral perturbation theory (ChPT). The result of $f_+(0)$ yields the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, $|V_{us}| = 0.2255(13)(4)$, where the first error comes from our calculation and the second from the experiment. This value is consistent with the ones determined from the unitarity of the CKM matrix and the $K_{l2}$ decay within one standard deviation, while it is slightly larger than recent lattice calculations by at most 1.5 $σ$. Furthermore, we evaluate the shape of the form factors and the phase space integral from our results. We confirm that those results are consistent with the experiment, and also $|V_{us}|$ determined with our phase space integral agrees with the one in the above.

preprint2020arXiv

$K_{l3}$ form factors in $N_f = 2+1$ QCD at physical point on large volume

We present our results of the $K_{l3}$ form factors on the volume whose spatial extent is more than $L=$10 fm, with the physical pion and kaon masses using the stout-smearing clover $N_f = 2+1$ quark action and Iwasaki gauge action at $a^{-1}\approx2.3$ GeV. The $K_{l3}$ form factor at zero momentum transfer is obtained from fit based on the next-to-leading (NLO) formula in SU(3) chiral perturbation theory. We estimate systematic errors of the form factor, mainly coming from the finite lattice spacing effect. We also determine the value of $|V_{us}|$ by combining our result with the experiment and check the consistency with the standard model prediction. The result is compared with the previous lattice calculations.

preprint2020arXiv

Nature of the phase transition for finite temperature $N_{\rm f}=3$ QCD with nonperturbatively O($a$) improved Wilson fermions at $N_{\rm t}=12$

We study the nature of the finite temperature phase transition for three-flavor QCD. In particular we investigate the location of the critical endpoint along the three flavor symmetric line in the light quark mass region of the Columbia plot. In the study, the Iwasaki gauge action and the nonperturvatively O($a$) improved Wilson-Clover fermion action are employed. We newly generate data at $N_{\rm t}=12$ and set an upper bound of the critical pseudoscalar meson mass in the continuum limit $m_{\rm PS,E}\lesssim 110$MeV.

preprint2020arXiv

Nucleon form factors and root-mean-square radii on a (10.8 fm$)^4$ lattice at the physical point

We present the nucleon form factors and root-mean-square (RMS) radii measured on a (10.8 fm$)^4$ lattice at the physical point. We compute the form factors at small momentum transfer region in $q^2\le 0.102$ GeV$^2$ with the standard plateau method choosing four source-sink separation times $t_{\rm sep}$ from 0.84 to 1.35 fm to examine the possible excited state contamination. We obtain the electric and magnetic form factors and their RMS radii for not only the isovector channel but also the proton and neutron ones without the disconnected diagram. We also obtain the axial-vector coupling and the axial radius from the axial-vector form factor. We find that these three form factors do not show large $t_{\rm sep}$ dependence in our lattice setup. On the other hand, the induced pseudoscalar and pseudoscalar form factors show the clear effects of the excited state contamination, which affect the generalized Goldberger-Treiman relation.

preprint2020arXiv

Restoration of chiral symmetry in cold and dense Nambu--Jona-Lasinio model with tensor renormalization group

We analyze the chiral phase transition of the Nambu--Jona-Lasinio model in the cold and dense region on the lattice developing the Grassmann version of the anisotropic tensor renormalization group algorithm. The model is formulated with the Kogut--Susskind fermion action. We use the chiral condensate as an order parameter to investigate the restoration of the chiral symmetry. The first-order chiral phase transition is clearly observed in the dense region at vanishing temperature with $μ/T\sim O(10^3)$ on a large volume of $V=1024^4$. We also present the results for the equation of state.

preprint2020arXiv

Tensor renormalization group study of two-dimensional U(1) lattice gauge theory with a $θ$ term

We make an analysis of the two-dimensional U(1) lattice gauge theory with a $θ$ term by using the tensor renormalization group. Our numerical result for the free energy shows good consistency with the exact one at finite coupling constant. The topological charge density generates a finite gap at $θ=π$ toward the thermodynamic limit. In addition finite size scaling analysis of the topological susceptibility up to $V=L\times L=1024\times 1024$ allows us to determine the phase transition at $θ=π$ is the first order.

preprint2019arXiv

Investigation of complex $ϕ^{4}$ theory at finite density in two dimensions using TRG

We study the two-dimensional complex $ϕ^{4}$ theory at finite chemical potential using the tensor renormalization group. This model exhibits the Silver Blaze phenomenon in which bulk observables are independent of the chemical potential below the critical point. Since it is expected to be a direct outcome of an imaginary part of the action, an approach free from the sign problem is needed. We study this model systematically changing the chemical potential in order to check the applicability of the tensor renormalization group to the model in which scalar fields are discretized by the Gaussian quadrature. The Silver Blaze phenomenon is successfully confirmed on the extremely large volume $V=1024^2$ and the results are also ensured by another tensor network representation with a character expansion.

preprint2016arXiv

A performance evaluation of CCS QCD Benchmark on the COMA (Intel(R) Xeon Phi$^{TM}$, KNC) system

The most computationally demanding part of Lattice QCD simulations is solving quark propagators. Quark propagators are typically obtained with a linear equation solver utilizing HPC machines. The CCS QCD Benchmark is a benchmark program solving the Wilson-Clover quark propagator, and is developed at the Center for Computational Sciences (CCS), University of Tsukuba. We optimized the benchmark program for a \Intel \XeonPhi (Knights Corner, KNC) system named "COMA (PACS-IX)" at CCS Tsukuba under the Intel Parallel Computing Center program. A single precision BiCGStab solver with the overlapped Restricted Additive Schwarz (RAS) preconditioner was implemented using SIMD intrinsics, OpenMP and MPI in the offload mode. With the reverse-offloading technique, we could reduce the communication and offloading overheads. We observed a performance of $\sim 200$ GFlops sustained for the Wilson-Clover hopping matrix multiplication on the lattice sizes larger than $24^3\times 32$ on a sinlge card of the COMA system. A good weak scaling perofmace was observed on the local lattice sizes larger than $24^3\times 32$.

preprint2016arXiv

Critical endline of the finite temperature phase transition for 2+1 flavor QCD around the SU(3)-flavor symmetric point

We investigate the critical endline of the finite temperature phase transition of QCD around the SU(3)-flavor symmetric point at zero chemical potential. We employ the renormalization-group improved Iwasaki gauge action and non-perturbatively $O(a)$-improved Wilson-clover fermion action. The critical endline is determined by using the intersection point of kurtosis, employing the multi-parameter, multi-ensemble reweighting method to calculate observables off the SU(3)-symmetric point, at the temporal size $N_{\rm T}$=6 and lattice spacing as low as $a \approx 0.19$ fm. We confirm that the slope of the critical endline takes the value of $-2$, and find that the second derivative is positive, at the SU(3)-flavor symmetric point on the Columbia plot parametrized with the strange quark mass $m_s$ and degenerated up-down quark mass $m_{\rm l}$.

preprint2016arXiv

Update on Nf=3 finite temperature QCD phase structure with Wilson-Clover fermion action

We present an update of the finite temperature phase structure analysis for three flavor QCD. In the study the Iwasaki gauge action and non-perturvatively O($a$) improved Wilson-Clover fermion action are employed. We discuss finite size scaling analysis including mixings of magnetization-like and energy-like observables. Preliminary results are shown of the continuum limit of the critical point using newly generated data at Nt=8,10, including estimates of the critical pseudo-scalar meson mass and critical temperature.

preprint2015arXiv

Curvature of the critical line on the plane of quark chemical potential and pseudo scalar meson mass for three-flavor QCD

We investigate the phase structure of three-flavor QCD in the presence of finite quark chemical potential $μ/T\lesssim1.2$ by using the non-perturbatively $O(a)$ improved Wilson fermion action on lattices with a fixed temporal extent $N_{\rm t}=6$ and varied spatial linear extents $N_{\rm s}=8,10,12$. Especially, we focus on locating the critical end point that characterizes the phase structure, and extracting the curvature of the critical line on the $μ$-$m_π$ plane. For Wilson-type fermions, the correspondence between bare parameters and physical parameters is indirect. Hence we present a strategy to transfer the bare parameter phase structure to the physical one, in order to obtain the curvature. Our conclusion is that the curvature is positive. This implies that, if one starts from a quark mass in the region of crossover at zero chemical potential, one would encounter a first-order phase transition when one raises the chemical potential.

preprint2015arXiv

Phase structure of $N_{\rm f}=3$ QCD at finite temperature and density by Wilson-Clover fermions

We investigate the phase structure of 3-flavor QCD in the presence of finite quark chemical potential by using Wilson-Clover fermions. To deal with the complex action with finite density, we adopt the phase reweighting method. In order to survey a wide parameter region, we employ the multi-parameter reweighting method as well as the multi-ensemble reweighting method. Especially, we focus on locating the critical end point that characterizes the phase structure. It is estimated by the kurtosis intersection method for the quark condensate. For Wilson-type fermions, the correspondence between bare parameters and physical parameters is indirect, thus we present a strategy to transfer the bare parameter phase structure to the physical one. We conclude that the curvature with respect to the chemical potential is positive. This implies that, if one starts from a quark mass in the region of crossover at zero chemical potential, one would encounter a first-order phase transition when one raises the chemical potential.

preprint2015arXiv

Scalar correlators near the 3-flavor thermal critical point

We investigate screening masses at both sides of the first order finite temperature transition with 3 quark flavors using the nonperturbatively improved clover fermion action and the Iwasaki gauge action. We have developed the method of hierarchical truncations with stochastic probing to accelerate the noise estimator for evaluating quark loops at every spatial lattice slices. At parameter values we study, the flavor singlet scalar meson has a screening mass about half of the pion screening mass. It becomes lighter as the system approaches the critical endpoint.

preprint2015arXiv

Study of quark mass dependence of binding energy for light nuclei in 2+1 flavor lattice QCD

We investigate the formation of light nuclei with the nuclear mass number less than or equal to four in 2+1 flavor QCD using a non-perturbative improved Wilson quark and Iwasaki gauge actions. The quark mass is decreased from our previous work to the one corresponding to the pion mass of 0.30 GeV. In each multi-nucleon channel, the energy shift of the ground state relative to the assembly of free nucleons is calculated on two volumes, whose spatial extents are 4.3 fm and 5.8 fm. From the volume dependence of the energy shift, we distinguish a bound state of multi nucleons from an attractive scattering state. We find that all the ground states measured in this calculation are bound states. As in the previous studies at larger $m_π$, our result indicates that at $m_π= 0.30$ GeV the effective interaction between nucleons in the light nuclei is relatively stronger than the one in nature, since the results for the binding energies are larger than the experimental values and a bound state appears in the dineutron channel, which is not observed in experiment. Possible sources of systematic error in this calculation are discussed.

preprint2014arXiv

Critical behavior of lattice Schwinger model with topological term at $θ=π$ using Grassmann tensor renormalization group

Lattice regularized Schwinger model with a so-called $θ$ term is studied by using the Grassmann tensor renormalization group. We perform the Lee-Yang and Fisher zero analyses in order to investigate the phase structure at $θ=π$. We find a first order phase transition at larger fermion mass. Both of the Lee-Yang zero and Fisher zero analyses indicate that the critical endpoint at which the first order phase transition terminates belongs to the Ising universality class.

preprint2014arXiv

Critical end point of Nf=3 QCD at finite temperature and density

We investigate the phase structure of 3-flavor QCD in the presence of finite quark chemical potential $aμ=0.1$ by using the Wilson-Clover fermion action. Especially, we focus on locating the critical end point that characterizes the phase structure. We do this by the kurtosis intersection method for the quark condensate. For Wilson-type fermions, the correspondence between bare parameters and physical parameters is indirect. Hence we present a strategy to transfer the bare parameter phase structure to the physical one.

preprint2014arXiv

Critical endpoint of finite temperature phase transition for three flavor QCD

We investigate the critical endpoint of finite temperature phase transition of $N_f=3$ QCD at zero chemical potential. We employ the renormalization-group improved Iwasaki gauge action and non-perturbatively $O(a)$-improved Wilson-clover fermion action. The critical endpoint is determined by using the intersection point of kurtosis for the temporal size $N_t$=4, 6, 8. Spatial sizes of $N_l$=6-16 ($N_t$=4), 10-24 ($N_t$=6), and 12-24 ($N_t$=8) are employed. We find that $N_t$=4 is out of the scaling region. Using results for $N_t$=6 and 8, and making linear extrapolations in $1/N_t^2$, we obtain $\sqrt{t_0}T_{\rm E}=0.0975(14)(8)$, $\sqrt{t_0}m_{\rm PS,E}=0.2254(52)(105)$ and $m_{\rm PS,E}/T_{\rm E}=2.311(63)(13)$, where the first error is statistical error, the second error is systematic error, and $m_{\rm PS}$ is the pseudo scalar meson mass. If one uses $1/\sqrt{t_0}=1.347(30)$ GeV reported by Borsanyi et al., one finds $T_{\rm E}=131(2)(1)(3)$ MeV, $m_{\rm PS,E}=304(7)(14)(7)$ MeV and $m_{\rm PS,E}/m_{\rm PS,E}^{\rm phys, sym}=0.739(17)(34)(17)$, where the third error comes from error of $\sqrt{t_0}$ and $m_{\rm PS}^{\rm phys, sym}=\sqrt{(m_π^2+2m_K^2)/3}$. Our current estimation of $\sqrt{t_0}m_{\rm PS,E}$ in the continuum limit is about 25% smaller than the SU(3) symmetric point.

preprint2014arXiv

Grassmann Tensor Renormalization Group Approach to One-Flavor Lattice Schwinger Model

We apply the Grassmann tensor renormalization group to the lattice regularized Schwinger model with one-flavor of the Wilson fermion. We study the phase diagram in the $(β,κ)$ plane performing a detailed analysis of the scaling behavior of the Lee-Yang zeros and the peak height of the chiral susceptibility. Our results strongly indicate that the whole range of the phase transition line starting from $(β,κ)=(0.0,0.380665(59))$ and ending at $(\infty,0.25)$ belongs to the two-dimensional Ising universality class similarly to the free fermion case.

preprint2013arXiv

Finite size scaling study of $N_{\text{f}}=4$ finite density QCD on the lattice

We explore the phase space spanned by the temperature and the chemical potential for 4-flavor lattice QCD using the Wilson-clover quark action. In order to determine the order of the phase transition, we apply finite size scaling analyses to gluonic and quark observables including plaquette, Polyakov loop and quark number density, and examine their susceptibility, skewness, kurtosis and Challa-Landau-Binder cumulant. Simulations were carried out on lattices of a temporal size fixed at $N_{\text{t}}=4$ and spatial sizes chosen from $6^3$ up to $10^3$. Configurations were generated using the phase reweighting approach, while the value of the phase of the quark determinant were carefully monitored. The $μ$-parameter reweighting technique is employed to precisely locate the point of the phase transition. Among various approximation schemes for calculating the ratio of quark determinants needed for $μ$-reweighting, we found the Taylor expansion of the logarithm of the quark determinant to be the most reliable. Our finite-size analyses show that the transition is first order at $(β, κ, μ/T)=(1.58, 0.1385, 0.584\pm 0.008)$ where $(m_π/m_ρ, T/m_ρ)=(0.822, 0.154)$. It weakens considerably at $(β, κ, μ/T)=(1.60, 0.1371, 0.821\pm 0.008)$ where $(m_π/m_ρ, T/m_ρ)=(0.839, 0.150)$, and a crossover rather than a first order phase transition cannot be ruled out.

preprint2013arXiv

Multi-nucleon bound states in $N_f=2+1$ lattice QCD

We report on our on-going effort to calculate the properties of light nuclei directly from quarks and gluons based on lattice QCD. After briefly introducing our motivations and aims, we describe our strategy of fixing the strange quark mass at its physical value and approaching the physical point for the up and down quark masses step by step from the region of heavy quark masses. A successful calculation for the pion mass of 0.51GeV is reviewed, and the status for a lighter pion mass of 0.30GeV is reported.

preprint2013arXiv

Zeros of QCD partition function from finite density lattices

Partition function zeros steer the critical behavior of a system. Studying four-flavor lattice QCD at finite temperature and density with the Wilson-clover fermion action and the Iwasaki gauge action using a phase-quenched fermion determinant, we combine statistics from multiple chemical potentials to improve sampling of the configuration space, and aim at unraveling the movement of zeros in finite systems. Preparing for further investigations, we discuss methods and criteria used to sieve through complex parameter space spanned by $(\text{Re}μ, \text{Im}μ)$ and $(\text{Re}μ, \text{Im}β)$, and present statistically robust zeros of the partition function.

preprint2012arXiv

Bound states of multi-nucleon channels in N_f=2+1 lattice QCD

We calculate the energies for multi-nucleon ground states with the nuclear mass number less than or equal to 4 in 2+1 flavor QCD at the lattice spacing of a = 0.09 fm employing a relatively heavy quark mass corresponding to m_pi = 0.51 GeV. We investigate the volume dependence of the energy shift of the ground state and the state of free nucleons to distinguish a bound state from attractive scattering states. From the investigation we conclude that ^4He, ^3He, deuteron and dineutron are bound at m_pi = 0.51 GeV. We compare their binding energies with those in our quenched studies and also with some recent investigations.

preprint2012arXiv

Finite size scaling for 4-flavor QCD with finite chemical potential

We explore the phase diagram spanned by the temperature and the chemical potential for 4-flavor QCD by the phase-reweighting approach. In order to determine the order of phase transition, we perform finite size scaling studies for various quantities, for example, susceptibility, kurtosis and Challa-Landau-Binder cumulant. At the parameter (beta=1.60, kappa=0.1371, c_sw=1.9655 and N_T=4), where the Kentucky group reported a first-order phase transition in their canonical simulation, we observe that the transition is consistent with being of first order.

preprint2012arXiv

Helium nuclei, deuteron and dineutron in 2+1 flavor lattice QCD

We calculate the binding energies for multi-nucleon bound states with the nuclear mass number less than or equal to 4 in 2+1 flavor QCD at the lattice spacing of a = 0.09 fm employing a relatively heavy quark mass corresponding to m_pi = 0.51 GeV. To distinguish a bound state from attractive scattering states, we investigate the volume dependence of the energy shift between the ground state and the state of free nucleons by changing the spatial extent of the lattice from 2.9 fm to 5.8 fm. We conclude that ^4He, ^3He, deuteron and dineutron are bound at m_pi = 0.51 GeV. We compare their binding energies with those in our quenched studies and also with several previous investigations.

preprint2011arXiv

On the phase of quark determinant in lattice QCD with finite chemical potential

We investigate the phase of the quark determinant with finite chemical potential in lattice QCD using both analytic and numerical methods. Applying the winding number expansion and the hopping parameter expansion to the logarithm of the determinant, we show that the absolute value of the phase has an upper bound that grows with the spatial volume but decreases exponentially with an increase in the temporal extent of the lattice. This analytic but approximate result is confirmed with a numerical study in four-flavor QCD in which the phase is calculated exactly. Since the phase is well controlled on lattices with larger time extents, we try the phase reweighting method in a region beyond $μ/T=1$ where the Taylor expansion method cannot be applied. Working in four-flavor QCD, we find a first-order like behavior on a $6^3\times 4$ lattice at $μ/T\approx 0.8$ which was previously observed by Kentucky group with the canonical method. We also show that the winding number expansion has a nice convergence property beyond $μ/T=1$. We expect that this expansion is useful to study the high density region of the QCD phase diagram at low temperatures.

preprint2008arXiv

PACS-CS results for 2+1 flavor lattice QCD simulation on and off the physical point

We report on the PACS-CS project focusing on a direct simulation of 2+1 flavor QCD on the physical point and chiral analysis of meson and baryon masses off the physical point with both the SU(2) and SU(3) chiral perturbation theories. Configurations are generated with the O(a)-improved Wilson quark action and the Iwasaki gauge action. The up-down quark is simulated by employing the DDHMC algorithm with several improvements and the UV-filtered PHMC algorithm is implemented for the strange quark. We investigate the convergence behaviors of the SU(2) and SU(3) chiral expansions up to NLO for the pseudoscalar meson sector, where the up-down quark mass ranges from 3 MeV to 24 MeV and the strange quark mass is chosen around the physical value. The fit results for the low energy constants are compared with those recently obtained by other groups. We also discuss the importance of the direct simulation at the physical point by comparing the physical quantities measured on the physical point with those estimated by the extrapolation method.

preprint2004arXiv

Perturbative Determination of Mass Dependent $O(a)$ Improvement Coefficients for the Vector and Axial Vector Currents with a Relativistic Heavy Quark Action

We carry out a perturbative determination of mass dependent renormalization factors and $O(a)$ improvement coefficients for the vector and axial vector currents with a relativistic heavy quark action, which we have designed to control $m_Qa$ errors by extending the on-shell $O(a)$ improvement program to the case of $m_Q \gg Λ_{\rm QCD}$ with $m_Q$ the heavy quark mass. We discuss what kind of improvement operators are required for the heavy-heavy and the heavy-light cases under the condition that the Euclidean rotational symmetry is not retained anymore because of the $m_Qa$ corrections. Our calculation is performed employing the ordinary perturbation theory with the fictitious gluon mass as an infrared regulator. We show that all the improvement coefficients are determined free from infrared divergences. Results of the renormalization factors and the improvement coefficients are presented as a function of $m_Q a$ for various improved gauge actions as well as the plaquette action.

preprint2003arXiv

One-loop calculation of mass dependent ${\cal O}(a)$ improvement coefficients for the relativistic heavy quarks on the lattice

We carry out the one-loop calculation of mass dependent ${\cal O}(a)$ improvement coefficients in the relativistic heavy quark action recently proposed, employing the ordinary perturbation theory with the fictitious gluon mass as an infrared regulator. We also determine renormalization factors and improvement coefficients for the axial-vector current at the one-loop level. It is shown that the improvement coefficients are infrared finite at the one-loop level if and only if the improvement coefficients in the action are properly tuned at the tree level.

preprint2001arXiv

Nucleon Decay Matrix Elements from Lattice QCD

We present a GUT-model-independent calculation of hadron matrix elements for all dimension-six operators associated with baryon number violating processes using lattice QCD. Our results cover all the matrix elements required to estimate the partial lifetimes of (proton,neutron)$\to$($π,K,η$) +(${\bar ν},e^+,μ^+$) decay modes. We point out the necessity of disentangling two form factors that contribute to the matrix elements; previous calculations did not make the separation, which led to an underestimate of the physical matrix elements. With a correct separation, we find that the matrix elements have values 3-5 times larger than the smallest estimates employed in phenomenological analyses of the nucleon decays, which gives stronger constraints on GUT models. We also find that the values of the matrix elements are comparable with the tree-level predictions of chiral Lagrangian.

preprint2001arXiv

Relativistic Heavy Quarks on the Lattice

Lattice QCD should allow quantitative predictions for the heavy quark physics from first principles. Up to now, however, most approaches have based on the nonrelativistic effective theory, with which the continuum limit can not be taken in principle. In this paper we investigate feasibility of relativistic approaches to the heavy quark physics in lattice QCD. We first examine validity of the idea that the use of the anisotropic lattice could be advantageous to control the $m_Q a$ corrections. Our perturbative calculation, however, reveals that this is not true. We instead propose a new relativistic approach to handle heavy quarks on the isotropic lattice. We explain how power corrections of $m_Q a$ can be avoided and remaining uncertainties are reduced to be of order $(aΛ_{QCD})^2$.

preprint1999arXiv

Perturbative renormalization factors of quark operators for domain-wall QCD

We calculate one-loop renormalization factors of several quark operators including bilinear, three- and four-quark operator for domain-wall fermion action. Since Green functions are constructed for external physical quark fields, our renormalization method is simple and can be easily applied to calculation of any quark operators. Our results show that these renormalized quark operators preserve several chiral properties of continuum massless QCD, which can be understood by the property of external quark line propagator.