Search

We present our investigation of SU(2) gauge theory with 8 flavours, and SU(3) gauge theory with 12 flavours. For the SU(2) case, at strong bare coupling, $\beta \lesssim 1.45$, the distribution of the lowest eigenvalue of the Dirac operator can be described by chiral random matrix theory for the Gaussian symplectic ensemble. Our preliminary result indicates that the chiral phase transition in this theory is of bulk nature. For the SU(3) theory, we use high-precision lattice data to perform the step-scaling study of the coupling, $g_{{\rm GF}}$, in the Gradient Flow scheme. We carefully examine the reliability of the continuum extrapolation in the analysis, and conclude that the scaling behaviour of this SU(3) theory is not governed by possible infrared conformality at $g_{{\rm GF}}^{2} \lesssim 6$., Comment: 7 pages, 4 figures. Contribution to the proceedings of the 33rd International Symposium on Lattice Field Theory, July 14-18, 2015, Kobe, Japan

We perform the step-scaling investigation of the running coupling constant, using the gradient-flow scheme, in SU(3) gauge theory with twelve massless fermions in the fundamental representation. The Wilson plaquette gauge action and massless unimproved staggered fermions are used in the simulations. Our lattice data are prepared at high accuracy, such that the statistical error for the renormalised coupling, g_GF, is at the subpercentage level. To investigate the reliability of the continuum extrapolation, we employ two different lattice discretisations to obtain g_GF. For our simulation setting, the corresponding gauge-field averaging radius in the gradient flow has to be almost half of the lattice size, in order to have this extrapolation under control. We can determine the renormalisation group evolution of the coupling up to g^2_GF ~ 6, before the onset of the bulk phase structure. In this infrared regime, the running of the coupling is significantly slower than the two-loop perturbative prediction, although we cannot draw definite conclusion regarding possible infrared conformality of this theory. Furthermore, we comment on the issue regarding the continuum extrapolation near an infrared fixed point. In addition to adopting the fit ansatz a'la Symanzik for performing this task, we discuss a possible alternative procedure inspired by properties derived from low-energy scale invariance at strong coupling. Based on this procedure, we propose a finite-size scaling method for the renormalised coupling as a means to search for infrared fixed point. Using this method, it can be shown that the behaviour of the theory around g^2_GF ~ 6 is still not governed by possible infrared conformality., Comment: 24 pages, 6 figures; Published version; Appendix A added for tabulating data; One reference included; Typos corrected

We present preliminary result for the step-scaling study of the coupling constant with the Yang-Mills gradient flow, in the twelve-favour SU(3) gauge theory. In this work, the lattice simulation is performed using unimproved staggered fermions and the Wilson plaquette gauge action, from which the gradient flow is also implemented. Imposing twisted boundary condition a'la t'Hooft and Parisi, our calculation is performed at zero fermion mass. The renormalised coupling constant is extracted via the computation of the energy density. In order to examine the reliability of the continuum extrapolation, we investigate this coupling constant using two different lattice discretisations. Our result shows that in order to control the systematic effects in the continuum extrapolation, it is necessary to implement a large enough gradient-flow time. In the current calculation, the gauge-field averaging radius corresponding to the flow time has to be as large as 40% of the lattice size., Comment: 7 pages, contribution to the proceedings of Lattice 2014

We present the investigation of the strong bare-coupling regime of SU(2) lattice gauge theory with 8 fermion flavors in the fundamental representation. The simulations are performed with unimproved staggered fermions and the plaquette gauge action. One bulk phase transition is observed through the measurement of the plaquette. The results of cold-start and hot-start simulations, as well as the hysteresis study, indicate the order of this transition can be weakly first order. Using the smeared Polyakov loops, and a method inspired by the constraint effective potential, we study the vacuum structure near the confining-deconfining phase transition. The Dirac operator eigenvalue spectrum is investigated, where further analysis is needed to clarify the properties of the chiral phase structure., Comment: 7 pages, 10 figures, talk presented at the 32nd International Symposium on Lattice Field Theory - Lattice 2014, 23-28 June, 2014, Columbia University New York, NY, USA

We present new results of our ongoing project on the investigation of the phase structure of the Higgs-Yukawa model at small and large bare Yukawa couplings. The critical exponents of the second order bulk phase transitions of this model are determined from finite-size analyses and compared to the pure O(4)-model to test for triviality and the possibility of having a non-Gaussian fixed point. In addition, we will present a first study of Higgs boson masses and fermion correlation functions., Comment: 7 pages, contribution to Lattice 2013

We present new data on our ongoing project on the investigation of the phase structure of the Higgs-Yukawa model at large bare Yukawa couplings. The data presented last year are extended in terms of statistics, the number of bare Yukawa couplings at existing, and new larger volumes. In addition, this study is extended by a finite temperature project at the physical top quark mass m_t =175 GeV and a hypothetical fourth generation top quark with a mass of m_t' =700 GeV ., Comment: 7 pages, 4 figures, Poster presented at The XXXth International Symposium on Lattice Field Theory - Lattice 2012, June 24-29, 2012, Cairns, Australia

Non-perturbative numerical lattice studies of the Higgs-Yukawa sector of the standard model with exact chiral symmetry are reviewed. In particular, we discuss bounds on the Higgs boson mass at the standard model top quark mass, and in the presence of heavy fermions. We present a comprehensive study of the phase structure of the theory at weak and very strong values of the Yukawa coupling as well as at non-zero temperature., Comment: 23 pages, 13 figures. Invited review for the special issue of Advances in High Energy Physics (AHEP) on Very Heavy Quarks at the LHC. Published version

We present details of a lattice study of infrared behaviour in SU(3) gauge theory with twelve massless fermions in the fundamental representation. Using the step-scaling method, we compute the coupling constant in this theory over a large range of scale. The renormalisation scheme in this work is defined by the ratio of Polyakov loops in the directions with different boundary conditions. We closely examine systematic effects, and find that they are dominated by errors arising from the continuum extrapolation. Our investigation suggests that SU(3) gauge theory with twelve flavours contains an infrared fixed point., Comment: 29 pages, 15 figures, 4 tables. Minor revision. Published version

In this article, we present an ongoing lattice study of the Higgs-Yukawa model, in the regime of strong-Yukawa coupling, using overlap fermions. We investigated the phase structure in this regime by computing the Higgs vacuum expectation value, and by exploring the finite-size scaling behaviour of the susceptibility corresponding to the magnetisation. Our preliminary results indicate the existence of a second-order phase transition when the Yukawa coupling becomes large enough, at which the Higgs vacuum expectation value vanishes and the susceptibility diverges., Comment: 7 pages, 3 figures, 1 table. Talk presented at the XXIX International Symposium on Lattice Field Theory - Lattice 2011, Lake Tahoe, California, USA

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 propose a new method for Hybrid Monte Carlo (HMC) simulations with odd numbers of dynamical fermions on the lattice. It employs a different approach from polynomial or rational HMC. In this method, gamma-five hermiticity of the lattice Dirac operators is crucial and it can be applied to Wilson, domain-wall, and overlap fermions. We compare HMC simulations with two degenerate flavors and (1 + 1) degenerate flavors using optimal domain-wall fermions. The ratio of the efficiency, (number of accepted trajectories) / (simulation time), is about 3:2. The relation between pseudofermion action of chirally symmetric lattice fermions in four-dimensional(overlap) and five-dimensional(domain-wall) representation are also analyzed., Comment: 13 pages

We present the first GPU-based conjugate gradient (CG) solver for lattice QCD with domain-wall fermions (DWF). It is well-known that CG is the most time-consuming part in the Hybrid Monte Carlo simulation of unquenched lattice QCD, which becomes even more computational demanding for lattice QCD with exact chiral symmetry. We have designed a CG solver for the general 5-dimensional DWF operator on NVIDIA CUDA architecture with mixed-precision, using the defect correction as well as the reliable updates algorithms. We optimize our computation by even-odd preconditioning in the 4D space-time lattice, plus several innovative techniques for CUDA kernels. For NVIDIA GeForce GTX 285/480, our CG solver attains 180/233 Gflops (sustained)., Comment: 7 pages, talk presented at The XXVIII International Symposium on Lattice Field Theory, Lattice 2010, June 14-19, 2010, Villasimius, Italy

We construct positive-definite pseudofermion actions for one fermion flavor in lattice field theory, for Wilson and domain-wall fermions respectively. The positive definiteness of these actions ensures that they can be simulated with the Hybrid Monte Carlo (HMC) method. For lattice QCD with optimal domain-wall quarks, we compare the efficiency of HMC simulations of 2-flavor and (1+1)-flavor, and find that the efficiency ratio is about 3:2., Comment: 7 pages. Talk presented at the XXVII International Symposium on Lattice Field Theory, July 26-31 2009, Peking University, Beijing, China

We present an overview of our project of simulation of unquenched lattice QCD with optimal domain-wall quarks, using a GPU cluster currently constituting of 16 units of Nvidia Tesla S1070 plus 64 graphic cards with Nvidia GTX285 (total 128 GPUs with 128 Teraflops peak), attaining sustained computing power of 15.36 Teraflops. The first production run in two-flavor QCD is on-going, using the Iwasaki gauge action on a set of lattices with sizes $ 16^3 \times (32,10,8,6,4) \times (16,32) $ at the lattice spacing $ a \sim 0.1$ fm, with eight sea quark masses down to $ m_\pi \simeq 200 $ MeV. We outline our simulation algorithm, and describe the present status of the production run. Preliminary results of pseudoscalar mass and decay constant are also presented., Comment: 7 pages, 4 figures, talk presented at The XXVII International Symposium on Lattice Field Theory - LAT2009, July 26-31 2009, Peking University, Beijing, China

We perform the first study of treating b, c, and s quarks as Dirac fermions in lattice QCD with exact chiral symmetry. On a 32^3 60 lattice with 1/a ~ 7.68 GeV, we compute point-to-point quark propagators, and measure the time-correlation functions for mesons with quark contents b_bbar, c_bbar, s_bbar, and c_cbar. The lowest-lying meson mass spectra, the pseudoscalar decay constants, and the b and c quark masses are determined., Comment: 7 pages, 1 figure, Proceedings of Lattice 2007, Regensburg, Germany, July 30 - August 4, 2007

We present the first study of treating b, c, and s quarks as Dirac fermions in lattice QCD with exact chiral symmetry. For 100 gauge configurations generated with single-plaquette action at beta = 7.2 on the 32^3 x 60 lattice, we compute point-to-point quark propagators for 33 quark masses in the range 0.01 < m_q a < 0.85, and measure the time-correlation function of (pseudo-)scalar, (pseudo-)vector, and tensor mesons, for symmetric and asymmetric quark masses respectively. The lowest-lying mass spectra of mesons with quark contents b_bbar, c_bbar, s_bbar, and c_cbar are determined, together with the pseudoscalar decay constants. Our results are sumarized in Tables 1-5. Furthermore, we also determine the b and c quark masses in the MS_bar scheme, m_b = 4.65(5) GeV, and m_c = 1.16(4) GeV., Comment: 12 pages

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 on the results of our calculation of meson correlators in a finite volume. The calculation is carried out in the quenched approximation near the chiral limit (down to Mq = 2.6 MeV) using the overlap fermion. For these small quark masses, the scalar and pseudo-scalar correlators are well approximated with a few hundred eigenmodes. The results for both connected and disconnected correlators are compared with the theoretical predictions of quenched chiral perturbation theory., Comment: 6 pages

We investigate the effects of low-lying fermion eigenmodes on the QCD partition function in the $\epsilon$-regime. The fermion determinant is approximated by a truncated product of low-lying eigenvalues of the overlap-Dirac operator. With two flavors of dynamical quarks, we observe that the lattice results for the lowest eigenvalue distribution, eigenvalue sum rules and partition function reproduce the analytic predictions made by Leutwyler and Smilga, which strongly depend on the topological charge of the background gauge configuration. The value of chiral condensate extracted from these measurements are consistent with each other. For one dynamical quark flavor, on the other hand, we find an apparent disagreement among different determinations of the chiral condensate, which may suggest the failure of the $\epsilon$-expansion in the absence of massless Nambu-Goldstone boson., Comment: 23 pages, 9 figures

We present a quenched calculation of meson correlators with the overlap fermion at very small masses 2.6--13 MeV. In this region the pion Compton wavelength is larger than our lattice size L\simeq 1.23fm and the system is in the so-called epsilon-regime of chiral perturbation theory. We found that the scalar and pseudo-scalar correlators are precisely approximated by a few hundred low-lying fermion eigenmodes in this regime, whereas axial-vector correlator receives significant contributions from higher eigenmodes. We also measure the disconnected pseudo-scalar correlator, which is well saturated with the low-lying modes. Matching these lattice data with the one-loop expressions for the correlators in quenched chiral perturbation theory, we evaluate the decay constant F_pi and the chiral condensate Sigma as well as the parameters m_0 and alpha, which describe the artifacts of the quenched approximation., Comment: 28pages,31figures,v2: references added, typos fixed

We investigate the effects of low-lying fermion modes on the QCD partition function in the epsilon-regime. With the overlap Dirac operator we calculate several tens of low-lying fermion eigenvalues on the quenched lattice. By partially incorporating the fermion determinant through the truncated determinant approximation, we calculate the partition function and other related quantities for Nf = 1 and compare them with the theoretical predictions obtained by Leutwyler and Smilga., Comment: 3 pages, 5 figures, Talk presented at Lattice2004(chiral), Fermilab, June 21-26, 2004