Source author record

Kazuyuki Kanaya

Kazuyuki Kanaya 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

22works
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

22 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.

preprint2023arXiv

The International Linear Collider: Report to Snowmass 2021

The International Linear Collider (ILC) is on the table now as a new global energy-frontier accelerator laboratory taking data in the 2030s. The ILC addresses key questions for our current understanding of particle physics. It is based on a proven accelerator technology. Its experiments will challenge the Standard Model of particle physics and will provide a new window to look beyond it. This document brings the story of the ILC up to date, emphasizing its strong physics motivation, its readiness for construction, and the opportunity it presents to the US and the global particle physics community.

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

General purpose lattice QCD code set Bridge++ 2.0 for high performance computing

Bridge++ is a general-purpose code set for a numerical simulation of lattice QCD aiming at a readable, extensible, and portable code while keeping practically high performance. The previous version of Bridge++ is implemented in double precision with a fixed data layout. To exploit the high arithmetic capability of new processor architecture, we extend the Bridge++ code so that optimized code is available as a new branch, i.e., an alternative to the original code. This paper explains our strategy of implementation and displays application examples to the following architectures and systems: Intel AVX-512 on Xeon Phi Knights Landing, Arm A64FX-SVE on Fujitsu A64FX (Fugaku), NEC SX-Aurora TSUBASA, and GPU cluster with NVIDIA V100.

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.

preprint2020arXiv

Four quark operators for kaon bag parameter with gradient flow

To study the CP-violation using the $K_0-\bar{K}_0$ oscillation, we need the kaon bag parameter which represents QCD corrections in the leading Feynman diagrams. The lattice QCD provides us with the only way to evaluate the kaon bag parameter directly from the first principles of QCD. However, a calculation of relevant four quark operators with theoretically sound Wilson-type lattice quarks had to carry a numerically big burden of extra renormalizations and resolution of extra mixings due to the explicit chiral violation. Recently, the Small Flow-time eXpansion (SFtX) method was proposed as a general method based on the gradient flow to correctly calculate any renormalized observables on the lattice, irrespective of the explicit violations of related symmetries on the lattice. To apply the SFtX method, we need matching coefficients, which relate finite operators at small flow-times in the gradient flow scheme to renormalized observables in conventional renormalization schemes. In this paper, we calculate the matching coefficients for four quark operators and quark bi-linear operators, relevant to the kaon bag parameter.

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

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

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.

preprint2013arXiv

The Legacy of Hideki Yukawa, Sin-itiro Tomonaga, and Shoichi Sakata: Some Aspects from their Archives

Hideki Yukawa, Sin-itiro Tomonaga and Shoichi Sakata pioneered nuclear and particle physics and left enduring legacies. Their friendly collaboration and severe competition laid the foundation to bring up the active postwar generation of nuclear and particle physicists in Japan. In this presentation we illustrate milestones of nuclear and particle physics in Japan from 1930's to mid-1940's which have been clarified in Yukawa Hall Archival Library, Tomonaga Memorial Room and Sakata Memorial Archival Library.

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.

preprint2011arXiv

Finite Temperature QCD on the Lattice -- Status 2010

In the last couple of years, there has been big progress in finite temperature QCD on the lattice. Large-scale dynamical simulations of 2+1 flavor QCD with various improved staggered quark actions have been started to produce results for various thermodynamic quantities which are extrapolated to the continuum limit at around physical quark masses, and thus are capable for a direct comparison with experiment. At the same time, the theoretical uneasiness with staggered-type lattice quarks motivated several groups to accelerate studies with Wilson-type quarks and lattice chiral quarks. In this review, I discuss these important developments in finite temperature QCD made in the past year.

preprint2008arXiv

Search for the Charmonium Dissociation Temperature with Variational Analysis in Lattice QCD

Charmonium dissociation temperatures are studied in a quenched anisotropic lattice QCD with standard plaquette gauge action and O(a) improved Wilson fermion action. Simulations are carried out at temperatures in the range 0.88T_c to 2.3T_c. From the meson correlators, we first subtract the contribution of constant mode, which was reported to mislead the analysis of, in particular, P wave signals. We then calculate effective masses and Bethe-Salpeter wave functions for ground (1S, 1P) and excited states (2S, 2P) using the multi-state variational analysis. To distinguish between bound states and scattering states, we apply two methods: First, we compare effective masses for charmonium correlation functions with finite spatial momenta under different spatial boundary conditions. Since the scattering state energies are sensitive to the boundary conditions, we expect finite volume effects when the charmonium dissociates. Second, we study the Bethe-Salpeter wave function, which should become broad when the charmonium bound state turns into a scattering state. With both methods, we have fond no clear evidences of dissociation for the ground and exited charmonium states up to 2.3T_c so far.

preprint2001arXiv

Elementary Particles on a Dedicated Parallel Computer

Numerical simulation of lattice QCD is now an indispensable tool for investigating non-perturbative properties of quarks at low energies. In this paper, I introduce a recent lattice QCD calculation of light hadron spectrum and quark masses performed on a dedicated parallel computer, CP-PACS. From a simulation of QCD with two flavors of dynamical quarks, we find $m_{ud}^{\bar{\rm MS}}(2{\rm GeV}) = 3.44^{+0.14}_{-0.22}$ MeV for the averaged mass of $u$ and $d$ quarks to reproduce the experimental masses of $π$ and $ρ$ mesons, and $m_{s}^{\bar{\rm MS}}(2{\rm GeV}) = 88^{+4}_{-6}$ MeV or $90^{+5}_{-11}$ MeV to reproduce the $K$ or $ϕ$ meson mass. These values are about 20--30% smaller than the previous estimates using the quenched approximation.

preprint1997arXiv

Developments in Finite Temperature QCD on the Lattice with Dynamical Quarks

I report on some recent developments in nonperturbative studies of finite temperature QCD with dynamical quarks on the lattice. I discuss new studies of improved lattice actions and their application to finite temperature QCD. I also summarize the status of lattice investigations about the order of the finite temperature QCD transition for the case of two flavors of degenerate light quarks, using both staggered and Wilson lattice fermions.

preprint1993arXiv

Interface tension in SU(3) lattice gauge theory at finite temperatures on an $N_t=2$ lattice

The surface tension $σ$ of the confined-deconfined interface is calculated in pure $SU(3)$ lattice gauge theory at finite temperatures employing the operator and integral methods on a lattice of a size $8^2\times N_z\times 2$ with $N_z=16$ and 40. Analyses of non-perturbative corrections in asymmetry response functions strongly indicate that the use of one-loop values for the response functions lead to an overestimate of $σ$ in the operator method. The operator method also suffers more from finite-size effects due to a finite thickness of the interface, leading us to conclude that the integral method yields more reliable values for $σ$. Our result with the integral method $σ/T_c^3=0.134(16)$ is consistent with earlier results and also with that obtained with a transfer matrix method. Result is also reported on $σ$ obtained on a lattice $18^2\times 48\times 4$ with the integral method.