Search

We study the implementation of the even-odd Wilson fermion matrix for lattice QCD simulations on the A64FX architecture. Efficient coding of the stencil operation is investigated for two-dimensional packing to SIMD vectors. We measure the sustained performance on the supercomputer Fugaku at RIKEN R-CCS and show the profiler result of our code, which may signal an unexpected source of slow-down in addition to the detailed efficiency of each part of the code., Comment: 10 pages, contribtuion to the International Workshop on Arm-based HPC: Practice and Experience (IWAHPCE-2023), held in conjunction with The International Conference on High Performance Computing in Asia-Pacific Region (HPC Asia 2023), Singapore, Feb 27 - March 2, 2023

Bridge++ is a general-purpose code set for lattice QCD simulations aiming at a readable, extensible, and portable code while keeping practical high performance. The new version 2.0 employs machine-dependent optimization, enabling flexible data layout in float/double precision, while it was fixed layout and only with the double precision in previous versions. We report the performance on supercomputer Fugaku with Arm A64FX-SVE architecture by Fujitsu., Comment: 10 pages, contribution to the 39th International Symposium on Lattice Field theory (LATTICE2022), 8th-13th August, 2022, Bonn, Germany. The new version of Bridge++ mentioned in this article is now availabe at https://bridge.kek.jp/Lattice-code/

We investigate color superconductivity on the lattice using the gap equation for the Cooper pair condensate. The weak coupling analysis is justified by choosing the physical size of the lattice to be smaller than the QCD scale, while keeping the aspect ratio of the lattice small enough to suppress thermal excitations. In the vicinity of the critical coupling constant that separates the superconducting phase and the normal phase, the gap equation can be linearized, and by solving the corresponding eigenvalue problem, we obtain the critical point and the Cooper pair condensate without assuming its explicit form. The momentum components of the condensate suggest spatially isotropic s-wave superconductivity with Cooper pairs formed by quarks near the Fermi surface. The chiral symmetry in the massless limit is spontaneously broken by the Cooper pair condensate, which turns out to be dominated by the scalar and the pseudo-scalar components. Our results provide useful predictions, in particular, for future lattice simulations based on methods to overcome the sign problem such as the complex Langevin method., Comment: 30 pages, 15 figures, 1 table, v2: A.3 modified, Ref. [57] added

We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial condition is linearly unstable to convection. We find in our 2D simulation that convection is instigated there indeed and expands inward after being full-fledged. The convection is then settled to a quasi-steady state in $\sim100\,\mathrm{ms}$, being sustained by the negative $Y_e$ gradient, which is in turn maintained by neutrino emissions. It enhances the luminosities and mean energies of all species of neutrinos compared to 1D. Taking advantage of the Boltzmann solver, we analyze the possible occurrence of the neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in the regions, where $Y_e$ is lower, and that the growth rate can be as high as $\sim 10^{-1}\,{\mathrm{cm}^{-1}}$., Comment: accepted to ApJ

We report an implementation of a multigrid solver for the Clover fermion on supercomputer Fugaku, which uses A64FX CPU with Arm architecture. On Fugaku, a highly optimized implementation of BiCGStab solver with domain decomposed preconditioner for the Clover fermion, called QCD Wide SIMD library (QWS), is available. We use the preconditioner in QWS as a smoother of the multigrid solver. Our implementation shows reasonably good scaling and performance. The code is developed by using Bridge++ code framework and its extension., Comment: 8 pages, 3 figures, talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We discuss the flavor number dependence of QCD at low temperature and high density by the complex Langevin method. In our previous work, the complex Langevin method is confirmed to satisfy the criterion for correct convergence in certain regions, such as $\mu_{\rm q} / T = 5.2-7.2$ on $8^3 \times 16$ and $\mu_{\rm q} / T = 1.6-9.6$ on $16^3 \times 32$ using $N_{\rm f} = 4$ staggered fermion at $\beta = 5.7$. We extend this study to more realistic flavor cases, $N_{\rm f} = 2, 2 + 1, 3$, using Wilson fermions. We present the flavor number dependence of the validity regions of the complex Langevin method and the quark number., Comment: 9 pages, 5 figures, presented at the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We develop a new method to investigate color superconductivity (CSC) on the lattice based on the Thouless criterion, which amounts to solving the linearized gap equation without imposing any ansatz on the structure of the Cooper pairs. We perform explicit calculations at the one-loop level with the staggered fermions on a $8^3 \times 128$ lattice and the Wilson fermions on a $4^3 \times 128$ lattice, which enables us to obtain the critical $\beta(=6/g^2)$ as a function of the quark chemical potential $\mu$, below which the CSC phase is expected to appear. The obtained critical $\beta$ has sharp peaks at the values of $\mu$ corresponding to the discretized energy levels of quarks similarly to what was observed in previous studies on simplified effective models. From the solution to the linearized gap equation, one can read off the flavor and spatial structures of the Cooper pairs at the critical $\beta$. In the case of massless staggered fermion, in particular, we find that the chiral $\mathrm{U}(1)$ symmetry of the staggered fermions is spontaneously broken by the condensation of the Cooper pairs., Comment: 9 pages, 4 figures, submitted to the proceedings for the 38th International Symposium on Lattice Field Theory - Lattice2021

It is expected that the color superconductivity (CSC) phase appears in QCD at low temperature and high density. On the basis of the lattice perturbation theory, a possible parameter region in which the CSC occurs has been predicted. In this work, we perform complex Langevin simulation on an $8^3\times 128$ lattice using four-flavor staggered fermions. We find, in particular, that the quark number has plateaux with respect to the chemical potential similar to our previous study, indicating the formation of the Fermi sphere. A diquark-antidiquark operator, which is an order parameter of color superconductivity, is formulated on the lattice using the U(1) noise. Our result for this operator is found to fluctuate violently when the Fermi surface coincides with the energy levels of quarks. We also discuss partial restoration of the chiral symmetry at high density., Comment: 9 pages, 4 figures, talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

A portable implementation of elaborated algorithm is important to use variety of architectures in HPC applications. In this work we implement and benchmark an algebraic multi-grid solver for Lattice QCD on three different architectures, Intel Xeon Phi, Fujitsu A64FX, and NVIDIA Tesla V100, in keeping high performance and portability of the code based on the object-oriented paradigm. Some parts of code are specific to an architecture employing appropriate data layout and tuned matrix-vector multiplication kernels, while the implementation of abstract solver algorithm is common to all architectures. Although the performance of the solver depends on tuning of the architecture-dependent part, we observe reasonable scaling behavior and better performance than the mixed precision BiCGSstab solvers., Comment: 17 pages, 7 figures. The 21st International Conference on Computational Science and Its Applications, ICCSA 2021

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., Comment: 6 pages, 6 figures. Talk by I.Kanamori at the XXXII IUPAP Conference on Computational Physics (CCP2021), Coventry, UK, 1-5 August 2021

We present results on the world's first over 100 PFLOPS single precision lattice QCD quark solver on the japanese new supercomputer Fugaku. We achieve a factor 38 time speedup from the supercomputer K on the same problem size, $192^4$, with 102 PFLOPS, 10% floating-point operation efficiency against single precision floating-point operation peak. The evaluation region is the single precision BiCGStab for a Clover-Wilson Dirac matrix with Schwarz Alternating Procedure domain decomposition preconditioning using Jacobi iteration for the local domain matrix inversion., Comment: 27 pages, 8 figures

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame., Comment: 28pages, 22figures, accepted for publication in the ApJ

The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics., Comment: 89pages, Prog. Theor. Exp. Phys. accepted

We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our code's capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code., Comment: 23 pages, 19 figures, accepted version

We demonstrate that the complex Langevin method (CLM) enables calculations in QCD at finite density in a parameter regime in which conventional methods, such as the density of states method and the Taylor expansion method, are not applicable due to the severe sign problem. Here we use the plaquette gauge action with $\beta = 5.7$ and four-flavor staggered fermions with degenerate quark mass $m a = 0.01$ and nonzero quark chemical potential $\mu$. We confirm that a sufficient condition for correct convergence is satisfied for $\mu /T = 5.2 - 7.2$ on a $8^3 \times 16$ lattice and $\mu /T = 1.6 - 9.6$ on a $16^3 \times 32$ lattice. In particular, the expectation value of the quark number is found to have a plateau with respect to $\mu$ with the height of 24 for both lattices. This plateau can be understood from the Fermi distribution of quarks, and its height coincides with the degrees of freedom of a single quark with zero momentum, which is 3 (color) $\times$ 4 (flavor) $\times$ 2 (spin) $=24$. Our results may be viewed as the first step towards the formation of the Fermi sphere, which plays a crucial role in color superconductivity conjectured from effective theories., Comment: 17 pages, 4 figures; published version

We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there., Comment: 29 pages, 22 figures, submitted to ApJ

Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The protoneutron star seen in the FS model is more compact than that in the LS model because the existence of multinuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the motion of matter determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semitransparent region., Comment: 30 pages, 24 figures, 1 explosion, published in ApJ

We explore the QCD phase diagram at finite density with four-flavor staggered fermions using the complex Langevin method, which is a promising approach to overcome the sign problem. In our previous work on an $8^3 \times 16$ lattice at $\beta = 5.7$ with the quark mass $m = 0.01$, we have found that the baryon number density has a clear plateau as a function of the chemical potential. In this study, we use a $16^3 \times 32$ lattice to reduce finite volume effects and find that the plateau structure survives. Moreover, the number of quarks in the plateau region turns out to be 24, which is exactly the same as the one obtained previously on the $8^3 \times 16$ lattice. We provide a simple interpretation of this number, which suggests that the Fermi sphere is starting to form., Comment: 7 pages, 2 figures, submitted to the proceedings for 37th International Symposium on Lattice Field Theory - Lattice2019

Monte Carlo studies of QCD at finite density suffer from the sign problem, which becomes easily uncontrollable as the chemical potential $\mu$ is increased even for a moderate lattice size. In this work we make an attempt to approach the high density low temperature region by the complex Langevin method (CLM) using four-flavor staggered fermions with reasonably small quark mass on a $8^3 \times 16$ lattice. Unlike the previous work on a $4^3 \times 8$ lattice, the criterion for correct convergence is satisfied within a wide range of $\mu$ without using the deformation technique. In particular, the baryon number density exhibits a plateau behavior consistent with the formation of eight baryons, and it starts to grow gradually at some $\mu$., Comment: 7 pages, 8 figures, talk presented at the 36th Annual International Symposium on Lattice Field Theory, July 22-28, 2018, East Lansing, MI, USA

Exploring the phase diagram of QCD at finite density is a challenging problem since first-principle calculations based on standard Monte Carlo methods suffer from the sign problem. As a promising approach to this issue, the complex Langevin method (CLM) has been pursued intensively.In this work, we investigate the applicability of the CLM in the vicinity of the deconfinement phase transition using the four-flavor staggered fermions. In particular, we look for a signal of the expected first order phase transition within the validity region of the CLM., Comment: 7 pages, 8 figures, proceedings for the 36th Annual International Symposium on Lattice Field Theory (LATTICE2018)

We investigate implementation of lattice Quantum Chromodynamics (QCD) code on the Intel AVX-512 architecture. The most time consuming part of the numerical simulations of lattice QCD is a solver of linear equation for a large sparse matrix that represents the strong interaction among quarks. To establish widely applicable prescriptions, we examine rather general methods for the SIMD architecture of AVX-512, such as using intrinsics and manual prefetching, for the matrix multiplication. Based on experience on the Oakforest-PACS system, a large scale cluster composed of Intel Xeon Phi Knights Landing, we discuss the performance tuning exploiting AVX-512 and code design on the SIMD architecture and massively parallel machines. We observe that the same code runs efficiently on an Intel Xeon Skylake-SP machine., Comment: 17 pages, 9 figures, talk given by I.K. at the Workshop Large Scale Computational Physics (LSCP 2018) in the 18th International Conference on Computational Science and its Applications (ICCSA 2018), 2-5 July 2018, Melbourne. arXiv admin note: text overlap with arXiv:1712.01505

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed., Comment: 24 pages, 23 figures, 0 explosion, published in ApJ

We investigate implementation of lattice Quantum Chromodynamics (QCD) code on the Intel Xeon Phi Knights Landing (KNL). The most time consuming part of the numerical simulations of lattice QCD is a solver of linear equation for a large sparse matrix that represents the strong interaction among quarks. To establish widely applicable prescriptions, we examine rather general methods for the SIMD architecture of KNL, such as using intrinsics and manual prefetching, to the matrix multiplication and iterative solver algorithms. Based on the performance measured on the Oakforest-PACS system, we discuss the performance tuning on KNL as well as the code design for facilitating such tuning on SIMD architecture and massively parallel machines., Comment: 8 pages, 12 figures. Talk given at LHAM'17 "5th International Workshop on Legacy HPC Application Migration" in CANDAR'17 "The Fifth International Symposium on Computing and Networking" and to appear in the proceedings

We report the performance of Wilson and Domainwall Kernels on a new Intel Xeon Phi Knights Landing based machine named Oakforest-PACS, which is co-hosted by University of Tokyo and Tsukuba University and is currently fastest in Japan. This machine uses Intel Omni-Path for the internode network. We compare performance with several types of implementation including that makes use of the Grid library. The code is incorporated with the code set Bridge++., Comment: 8 pages, 9 figures, Proceedings for the 35th International Symposium on Lattice Field Theory (Lattice 2017)

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods., Comment: Accepted for publication in ApJ

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 $\mu$ 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., Comment: 7 pages, 4 figures. Proceedings of the 34th International Symposium on Lattice Field Theory (Lattice 2016)

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe., Comment: Accepted for publication in ApJS

We numerically study the SU(2) gauge theory with two dynamical flavors of the domain-wall fermions in fundamental representation. The meson spectra and the residual mass are measured on three lattice volumes and at two values of gauge coupling so as to investigate the finite volume effect. On generated configurations, eigenvalues of the overlap fermion operator are determined and compared to the random matrix theory. To quantify the effect of violation of the exact chiral symmetry, we measure the correlation between the eigenvectors of the domain-wall and the overlap operators., Comment: 7 pages, 9 figures. Proceedings of 31st International Symposium on Lattice Field Theory (LATTICE 2013), July 29 - August 3, 2013, Mainz, Germany

We examine the axial U(1) symmetry near and above the finite temperature phase transition in two-flavor QCD using lattice QCD simulations. Although the axial U(1) symmetry is always violated by quantization, (i.e.) the chiral anomaly, the correlation functions may manifest effective restoration of the symmetry in the high temperature phase. We explicitly study this possibility by calculating the meson correlators as well as the Dirac operator spectral density near the critical point. Our numerical simulations are performed on a $16^3\times 8$ lattice with two flavors of dynamical quarks represented by the overlap fermion formalism. Chiral symmetry and its violation due to the axial anomaly is manifestly realized with this formulation, which is a prerequisite for the study of the effective restoration of the axial U(1) symmetry. In order to avoid discontinuity in the gauge configuration space, which occurs for the exactly chiral lattice fermions, the simulation is confined in a fixed topological sector. It induces finite volume effect, which is well described by a formula based on the Fourier transform from the $\theta$-vacua. We confirm this formula at finite temperature by calculating the topological susceptibility in the quenched theory. Our two flavor simulations show degeneracy of the meson correlators and a gap in the Dirac operator spectral density, which implies that the axial U(1) symmetry is effectively restored in the chirally symmetric phase., Comment: 28 pages, 8 figures

We consider the simulation of finite temperature QCD with two flavors of dynamical overlap fermions in order to study the suppression of the axial U(1) symmetry breaking at the chiral phase transition point. As a preliminary study, pure gauge simulations are performed to investigate how fixing the topology affects physical quantities like the topological susceptibility, $\chi_t$, at finite temperature, showing that it is possible to reconstruct known results from the fixed topology sector. First results on the degeneracy of meson correlators in the high temperature QGP sector are shown., Comment: 7 pages, 4 figures. XXIX International Symposium on Lattice Field Theory, July 10 - 16 2011, Squaw Valley, Lake Tahoe, California

Incorporated with twisted boundary condition, Polyakov loop correlators can give a definition of the renormalized coupling. We employ this scheme for the step scaling method (with step size s = 2) in the search of conformal fixed point of SU(3) gauge theory with 12 massless flavors. Staggered fermion and plaquette gauge action are used in the lattice simulation with six different lattice sizes, L/a = 20, 16, 12, 10, 8 and 6. For the largest lattice size, L/a = 20, we used a large number of Graphics Processing Units (GPUs) and accumulated 3,000,000 trajectories in total. We found that the step scaling function sigma (u) is consistent with u in the low-energy region. This means the existence of conformal fixed point. Some details of our analysis and simulations will also be presented., Comment: 7 pages, 7 figures

We study the infrared properties of SU(3) gauge theory coupled to 12 massless Dirac fermions in the fundamental representation. The renormalized running coupling constant is calculated in the Twisted Polyakov loop scheme on the lattice. From the step-scaling analysis, we find that the infrared behavior of the theory is governed by a non-trivial fixed point., Comment: This paper has been withdrawn due to the major revision from a letter into a full paper, which will be submitted as a new paper. In the new paper, the author list will also change

We measure the renormalized coupling in the Twisted Polyakov loop scheme for SU(3) gauge theory coupled with $N_f=12$ fundamental fermions. To find the infrared fixed point of this theory, we focus on the step scaling function for the renormalized coupling. We take the continuum limit using the linear function of $(a/L)^2$ and a constant fit function. We find that there is a sizeable systematic error due to the choice of the continuum extrapolation function in the low energy region. We will give some directions to reduce the systematic errors., Comment: 7pages, To appear in the proceedings for the 28th International Symposium on Lattice Field Field: Lattice 2010

We present our preliminary study of the SU(2) gauge theory with 8 flavors of fermions in fundamental representation. This theory could be a candidate of the gauge theory with conformal fixed point. By using Wilson/Polyakov loop in a finite volume with twisted boundary conditions, we study the renormalization group flow of the gauge coupling constant. Our calculation gives consistent result with the perturbative prediction of the running coupling in the weak coupling region. We investigate a possible signal for conformal behavior in the strong coupling region., Comment: 7 pages, 5figures, Presented at The XXVIII International Symposium on Lattice Field Theory, Lattice2010

We present some preliminary results of the project on finite temperature QCD with overlap fermions at KEK. We performed a series of simulations to assess the effects of fixing the topological sector at finite temperature and we will show the first calculations of topological susceptibility and meson masses for quenched and full QCD., Comment: 7 pages. 5 figures. Contribution to the proceedings of The XXVIII International Symposium on Lattice Field Theory, Lattice2010, June 14-19, 2010, Villasimius, Italy

We present a non-perturbative study of the running coupling constant in the Twisted Polyakov Loop (TPL) scheme. We investigate how the systematic and statistical errors can be controlled {\it via} a feasibility study in SU(3) pure Yang-Mills theory. We show that our method reproduces the perturbative determination of the running coupling in the UV. In addition, our numerical result agrees with the theoretical prediction of this coupling constant in the IR. We also present our preliminary results for $N_f=12$ QCD, where an IR fixed point may be present., Comment: 8pages, Proceedings at the XXVII International Symposium on Lattice Field Theory, July 26-31 2009, Peking University, Beijing, China

We propose a new renormalization scheme of the running coupling constant in general gauge theories using the Wilson loops. The renormalized coupling constant is obtained from the Creutz ratio in lattice simulations and the corresponding perturbative coefficient at the leading order. The latter can be calculated by adopting the zeta-function resummation techniques. We perform a benchmark test of our scheme in quenched QCD with the plaquette gauge action. The running of the coupling constant is determined by applying the step-scaling procedure. Using several methods to improve the statistical accuracy, we show that the running coupling constant can be determined in a wide range of energy scales with relatively small number of gauge configurations., Comment: 30pages, figs and comments added,reference added(v3)

We propose a new method to compute the running coupling constant of gauge theories on the lattice. We first give the definition of the running coupling in the new scheme using the Wilson loops in a finite volume, and explain how the running of the coupling constant is extracted from the measurement of the volume dependence. The perturbative calculation of the renormalization constant to define the scheme is also given at the leading order. As a benchmark test of the new scheme we apply the method the case of the quenched QCD. We show the preliminary result from our numerical simulations which are carried out with plaquette gauge action for various lattice sizes and bare lattice couplings. With techniques to improve the statistical accuracy, we show that we can determine the non-perturbative running of the coupling constant in a wide range of the energy scale with relatively small number of gauge configurations in our scheme. We compare our lattice data of the running coupling constant with perturbative renormalization group evolution at one- and two-loop order, and confirm the consistency between them at high energy., Comment: 13 pages, 6 figures, based on talks given by E. Itou and M. Kurachi at The XXVI International Symposium on Lattice Field Theory, July 14 - 19, 2008, Williamsburg, Virginia, USA; barely visible changes in Figs. 5 and 6

The $B^* B\pi$ coupling is a fundamental parameter of chiral effective Lagrangian with heavy-light mesons and can constrain the chiral behavior of $f_B$, $B_B$ and the $B\to \pi l \nu$ form factor in the soft pion limit. We compute the $B^* B \pi $ coupling with the static heavy quark and the $O(a)$-improved Wilson light quark. Simulations are carried out with $n_f=2$ unquenched $12^3\times 24$ lattices at $\beta=1.80$ and $16^3\times 32$ lattices at $\beta=1.95$ generated by CP-PACS collaboration. To improve the statistical accuracy, we employ the all-to-all propagator technique and the static quark action with smeared temporal link variables following the quenched study by Negishi {\it et al.}. These methods successfully work also on unquenched lattices, and determine the $B^*B\pi$ coupling with 1--2% statistical accuracy on each lattice spacing., Comment: 19pages,26figures

The B* B pi coupling is a fundamental parameter of chiral effective Lagrangian with heavy-light mesons and can constrain the B->pi l nu form factor in the soft pion limit which will be useful for precise determination of |Vub|. We compute the B* B pi coupling with the static heavy quark and the O(a)-improved Wilson light quark. Simulations are carried out with n_f=2 unquenched 12^3x24 lattices at beta=1.80 generated by CP-PACS collaboration. Following the quenched study by Negishi et al., we employ the all-to-all propagator with 200 low eigenmodes as well as HYP smeared link to improve the statistical accuracy., Comment: 7 pages, 5 figures, talk presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4, 2007, Regensburg, Germany Journal-ref: PoS(LATTICE 2007)365

I report the status of the dynamical overlap fermion project by the JLQCD and TWQCD collaborations. So far, the simulations have been completed with two flavors of overlap sea quarks in a wide range of sea quark mass corresponding the $p$-regime and the $\epsilon$-regime on a $16^3\times 32$ lattice at $a=$0.12 fm. More recently, runs with 2+1 flavors of sea quarks have also started. This talk mainly discusses the physics results on the $N_f=2$ lattice after describing the lattice formulation and algorithms., Comment: 15 pages, 14 figures, contribution to 25th International Symposium on Lattice Field Theory, Regensburg, Germany, 30 Jul - 4 Aug 2007

We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the numerical costs. Our main target is a $N_f=2$ simulation with overlap fermion on a $16^3\times 32$ lattice at lattice spacing around 0.12 fm., Comment: 7 pages, 8 figures, talk given at Lattice 2006 (Algorithms, Machines, and Networks), to appear in Proceedings Of Science, POS(LAT2006)031

We report our exploratory study on the matching condition of HQET with QCD including 1/m corrections. We introduce a new observable from the dependence of the heavy-light effective energy on the twisted boundary condition parameter $\theta$, which could be used to match the kinetic term $\vec{D}^2/(2m)$. Carrying out quenched QCD simulations for fixed lattice spacing in small volumes with O(a)-improved Wilson fermions, we study the 1/m dependence of this observable, from which the static limit and 1/m coefficient can be extracted. We also compare our preliminary result with HQET., Comment: 6 pages, 2 figures, poster presented at Lattice 2005 (Heavy quark physics), PoS style

We make remarks on the Maximum Entropy Method (MEM) for studies of the spectral function of hadronic correlators in finite temperature lattice QCD. We discuss the virtues and subtlety of MEM in the cases that one does not have enough number of data points such as at finite temperature. Taking these points into account, we suggest several tests which one should examine to keep the reliability for the results, and also apply them using mock and lattice QCD data., Comment: 6 pages, 3 figures, poster presented at Lattice 2005 (Nonzero temperature and density)

We study charmonium properties below and above $T_c$ up to 1.8$T_c$, on quenched anisotropic lattices. Information of the spectral functions is extracted using the maximum entropy method and the constrained curve fitting. We also calculate the color singlet and averaged free energies and evaluate the charmonium spectrum with the potential model analysis. The relation between the lattice result of the spectral function analysis and the potential model is discussed., Comment: Lattice2004(non-zero), 3 pages, 4 figures