Search

We study the $U(1)_A$ anomaly at high temperatures of $N_f=2+1$ lattice QCD with chiral fermions. Gauge ensembles are generated with M\"obius domain-wall (MDW) fermions, and the measurements are reweighted to those with overlap fermions. We report on the results for the Dirac spectra, the $U(1)_A$ susceptibility, and the topological susceptibility in the temperature range of $T=136$, $153$, $175$, and $204$ MeV, where the up and down quark masses are set to be near the physical point as well as at lighter or heavier masses., Comment: 10 pages, 4 figures, 1 table; Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

We present an update on the study of the QCD phase transition with 3 flavors of M\"obius domain wall fermions at zero baryon density. We performed simulations on lattices of size $36^3\times12\times16$ and $24^3\times12\times32$ with a variety of quark masses at a fixed lattice spacing $a=0.1361(20)$ fm, which correspond to a temperature 121(2) MeV. By analyzing the chiral condensate, chiral susceptibilitities and Binder cumulant on $36^3\times12\times16$ lattices together with the result obtained from our previous study on $24^3\times12\times16$ lattices, we identified a crossover occurring at quark mass around $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 3-4$ MeV for this temperature. Besides, we show the effects of residual chiral symmetry breaking on chiral condensate and chiral susceptibilities between $L_s=16$ and 32., Comment: 10 pages, 6 figures, Contribution to the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31 to August 4, 2023, Fermilab, USA

We investigate the finite temperature QCD phase transition with three degenerate quark flavors using M\"{o}bius domain wall fermions. To explore the order of phase transition on the lower left corner of Columbia plot and if possible, to locate the critical endpoint we performed simulations at temperatures around 181 and 121 MeV with lattice spacing $a=0.1361(20)$ fm corresponding to temporal lattice extent $N_{\tau}=8,12$ with varying quark mass for two different volumes with aspect ratios $N_{\sigma}/N_{\tau}$ ranging from 2 to 3. By analyzing the volume and mass dependence of the chiral condensate, disconnected chiral susceptibility and Binder cumulant we find that there is a crossover at $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 44\, \mathrm{MeV}$ for $\mathrm{T_{pc}}\sim$ 181 MeV, At temperature 121 MeV, the binder cumulant suggests a crossover at $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 3.7\, \mathrm{MeV}$, although a study of volume dependence would be important to confirm this., Comment: 9 pages, 5 figures, Contribution to the 39th International Symposium on Lattice Field Theory (Lattice 2022), 08-13 August 2022, Bonn, Germany

We study the $U(1)_A$ anomaly in the high-temperature phase of $N_f=2+1$ lattice QCD with chiral fermions. Gauge ensembles are generated with M\"obius domain-wall (MDW) fermions, and in the measurements the determinant is reweighted to that of overlap fermions. We report the results for the overlap Dirac spectrum, $U(1)_A$ susceptibility, and topological susceptibility at $T=204$ and $175$ MeV., Comment: 7 pages, 4 figures, 1 table; Proceedings of 38th International Symposium on Lattice Field Theory (LATTICE2021), 26-30 July 2021, Zoom/Gather@MIT, USA

Simulations for the thermodynamics of the 2+1 flavor QCD are performed employing chiral fermions. The use of M\"obius domain-wall fermions with stout-link smearing is more effective on the finer lattices where all the relevant chiral symmetries are realized more accurately. We report on the initial simulations near the (pseudo) critical point using the line of constant physics with an average $ud$ quark mass slightly heavier than physical at $a\lesssim 0.1$ fm., Comment: 9 pages, 8 figures, 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@MIT

We present preliminary results for the kaon semileptonic form factors using the PACS10 configurations, whose physical volume is more than (10 fm)$^3$ at the physical point with the lattice spacings of 0.085 and 0.064 fm. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared Clover quark action. For the continuum extrapolation, we calculate the form factors with the local and conserved vector currents. The form factors in the two lattice spacings are extrapolated to the continuum limit using a fit function based on the NLO SU(3) ChPT formula with terms corresponding to finite lattice spacing effects. The value of $|V_{us}|$ is determined using our preliminary result of the form factor at the zero momentum transfer in the continuum limit. The result is compared with recent lattice results, and also predictions of the standard model from the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix., Comment: 9 pages, 4 figures, proceedings of the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

We determine the properties of generalised parton distributions (GPDs) from a lattice QCD calculation of the off-forward Compton amplitude (OFCA). By extending the Feynman-Hellmann relation to second-order matrix elements at off-forward kinematics, this amplitude can be calculated from lattice propagators computed in the presence of a background field. Using an operator product expansion, we show that the deeply-virtual part of the OFCA can be parameterised in terms of the low-order Mellin moments of the GPDs. We apply this formalism to a numerical investigation for zero-skewness kinematics at two values of the soft momentum transfer, $t = -1.1, -2.2 \;\text{GeV}^2$, and a pion mass of $m_{\pi}\approx 470\;\text{MeV}$. The form factors of the lowest two moments of the nucleon GPDs are determined, including the first lattice QCD determination of the $n=4$ moments. Hence we demonstrate the viability of this method to calculate the OFCA from first principles, and thereby provide novel constraint on the $x$- and $t$-dependence of GPDs., Comment: 20 pages, 7 figures

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

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., Comment: 7 pages, 5 figures, Talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

We have initiated a program to compute the Compton amplitude from lattice QCD with the Feynman-Hellman method. This amplitude is related to the structure function via a Fredholm integral equation of the first kind. It is known that these types of equations are inherently ill--posed - they are, e.g., extremely sensitive to perturbations of the system. We discuss two methods which are candidates to handle these problems: the model free inversion based on singular value decomposition and one Bayesian type approach. We apply the Bayesian method to currently available lattice data for the Compton amplitude., Comment: 7 pages, 5 figures, poster presented at the 37th International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

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., Comment: 14 pages, 11 figures

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 $\sigma$. 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., Comment: 34 pages and 11 figures

We study the two-dimensional complex $\phi^{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., Comment: 17 pages, 9 figures

We investigate the critical endpoints of the finite temperature phase transition of 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 endpoints are determined by using the intersection point of kurtosis, employing the multi-parameter, multi-ensemble reweighting method. We present results for the critical endline at $N_{\rm T}$ = 6 and the continuum extrapolation for the critical endpoint of the SU(3)-flavor symmetric point., Comment: Typos corrected, references added. 8 pages, 5 figures, talk presented at 37th International Symposium on Lattice Field Theory - Lattice2019, 16 - 22 June 2019, Wuhan, China

We report our study on the critical endpoint of the finite temperature phase transition in 4-flavor QCD with Wilson-Clover fermions. Using the kurtosis intersection method, we determined the critical endpoint on lattices with $N_t$ = 4, 6 and 8. Our continuum extrapolated results show that the pseudo-scalar meson mass at the critical endpoint, $m_\mathrm{PS,E}$, for 4-flavor is clearly larger than that for 3-flavor. We also compared our results to those with staggered fermions and found that $m_\mathrm{PS,E}/T_E$ for 4-flavor with Wilson fermions might remain finite even in the continuum limit in contrast to that with staggered fermions, where $m_\mathrm{PS,E}/T_E$ is very close to zero, which suggests that the difference between Wilson and staggered fermions is at least not due to the rooting., Comment: 7 pages, 6 figures, talk presented in the 36th International Symposium on Lattice Field Theory (Lattice 2018), 22-28 July, 2018, Michigan State University, East Lansing, Michigan, USA

The tensor renormalization group attracts great attention as a new numerical method that is free of the sign problem. In addition to this striking feature, it also has an attractive aspect as a coarse-graining of space-time; the computational cost scales logarithmically with the space-time volume. This fact allows us to aggressively approach the thermodynamic limit. While taking this advantage, we study the critical coupling of the two dimensional $\phi^{4}$ theory on large and fine lattices. We present the numerical results along with the extrapolation procedure to the continuum limit and compare them with the previous ones by Monte Carlo simulations., Comment: 7 pages, 5 figures, talk presented at the 36th International Symposium on Lattice Field Theory (Lattice 2018), 22-28 July, 2018, Michigan, USA

We make a detailed analysis of the spontaneous $Z_{2}$-symmetry breaking in the two dimensional real $\phi^{4}$ theory with the tensor renormalization group approach, which allows us to take the thermodynamic limit easily and determine the physical observables without statistical uncertainties. We determine the critical coupling in the continuum limit employing the tensor network formulation for scalar field theories proposed in our previous paper. We obtain $\left[ \lambda / \mu_{\mathrm{c}}^{2} \right]_{\mathrm{cont.}} = 10.913(56)$ with the quartic coupling $\lambda$ and the renormalized critical mass $\mu_{\mathrm{c}}$. The result is compared with previous results obtained by different approaches., Comment: 21 pages, 11 figures

Early B-physics experiments have left us with a number of puzzles in heavy flavour physics. New lattice calculations (with a greater understanding of QCD effects in the Standard Model) will be needed to support the increase in experimental precision to be achieved by upcoming experiments such as Belle II. We extend the CSSM/UKQCD/QCDSF studies of SU(3) flavour breaking effects by presenting new results for the decay constants $f_B$ and $f_{B_s}$., Comment: 7 pages, 3 figures, conference proceedings for talk presented at the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

We apply the projective truncation technique to the tensor renormalization group (TRG) algorithm in order to reduce the computational cost from $O(\chi^6)$ to $O(\chi^5)$, where $\chi$ is the bond dimension, and propose three kinds of algorithms for demonstration. On the other hand, the technique causes a systematic error due to the incompleteness of a projector composed of isometries, and in addition requires iteration steps to determine the isometries. Nevertheless, we find that the accuracy of the free energy for the Ising model on a square lattice is recovered to the level of TRG with a few iteration steps even at the critical temperature for $\chi$ = 32, 48, and 64., Comment: 11 pages, 17 figures

It has become increasingly important to include one or more individual flavours of dynamical fermion in lattice QCD simulations. This is due in part to the advent of QCD+QED calculations, where isospin symmetry breaking means that the up, down, and strange quarks must be treated separately. These single-flavour pseudofermions are typically implemented as rational approximations to the inverse of the fermion matrix, using the technique known as Rational Hybrid Monte Carlo (RHMC). Over the years, a wide range of methods have been developed for accelerating simulations of two degenerate flavours of pseudofermion, while there are comparatively fewer such techniques for single-flavour pseudofermions. Here, we investigate two different filtering methods that can be applied to RHMC for simulating single-flavour pseudofermions, namely polynomial filtering (PF-RHMC), and filtering via truncations of the ordered product (tRHMC). A novel integration step-size tuning technique based on the characteristic scale is also introduced. Studies are performed on two different lattice volumes, demonstrating that one can achieve significant reductions in the computational cost of single-flavour simulations with these filtering techniques., Comment: 14 pages, 17 figures, 20 tables; to be submitted to Computer Physics Communications

Supersymmetric models with spontaneous supersymmetry breaking suffer from the notorious sign problem in stochastic approaches. By contrast, the tensor network approaches do not have such a problem since they are based on deterministic procedures. In this work, we present a tensor network formulation of the two-dimensional lattice $\mathcal{N}=1$ Wess-Zumino model while showing that numerical results agree with the exact solutions for the free case., Comment: 31 pages, 10 figures, revised version

We develop calculational method for fermionic Green functions in the framework of Grassmann higher-order tensor renormalization group. The validity of the method is tested by applying it to three-dimensional free Wilson fermion system. We compare the numerical results for chiral condensate and two-point correlation functions with the exact ones obtained by analytical methods.

Filtering algorithms for two degenerate quark flavours have advanced to the point that, in 2+1 flavour simulations, the cost of the strange quark is significant compared with the light quarks. This makes efficient filtering algorithms for single flavour actions highly desirable, in particular when considering 1+1+1 flavour simulations for QED+QCD. Here we discuss methods for filtering the RHMC algorithm that are implemented within BQCD, an open-source Fortran program for Hybrid Monte Carlo simulations., Comment: 8 pages, 3 figures, proceedings of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We present an update of BQCD, our Hybrid Monte Carlo program for simulating lattice QCD. BQCD is one of the main production codes of the QCDSF collaboration and is used by CSSM and in some Japanese finite temperature and finite density projects. Since the first publication of the code at Lattice 2010 the program has been extended in various ways. New features of the code include: dynamical QED, action modification in order to compute matrix elements by using Feynman-Hellman theory, more trace measurements, a more flexible integration scheme, polynomial filtering, term-splitting for RHMC, and a portable implementation of performance critical parts employing SIMD., Comment: Poster presented at the 35th International Symposium on Lattice Field Theory, Granada, Spain, 18-24 June 2017

We study a tensor network formulation of the two dimensional lattice $\mathcal{N}=1$ Wess-Zumino model with Wilson derivatives for both fermions and bosons. The tensor renormalization group allows us to compute the partition function without the sign problem, and basic ideas to obtain a tensor network for both fermion and scalar boson systems were already given in previous works. In addition to improving the methods, we have constructed a tensor network representation of the model including the Yukawa-type interaction of Majorana fermions and real scalar bosons. We present some numerical results., Comment: 8 pages, 4 figures, talk presented at the 35th International Symposium on Lattice Field Theory (Lattice 2017), 18-24 June 2017, Granada, Spain

We study the critical point for finite temperature Nf=3 QCD using several temporal lattice sizes up to 10. In the study, the Iwasaki gauge action and non-perturbatively O(a) improved Wilson fermions are employed. We estimate the critical temperature and the upper bound of the critical pseudo-scalar meson mass., Comment: 8 pages, 7 figures, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We study the finite temperature phase structure for three-flavor QCD with a focus on locating the critical point which separates crossover and first order phase transition region in the chiral regime of the Columbia plot. In this study, we employ the Iwasaki gauge action and the non-perturvatively O($a$) improved Wilson-Clover fermion action. We discuss the finite size scaling analysis including the mixing of magnetization-like and energy-like observables. We carry out the continuum extrapolation of the critical point using newly generated data at $N_{\rm t}=8$, $10$ and estimate the upper bound of the critical pseudo-scalar meson mass $m_{\rm PS,E} \lesssim 170 {\rm MeV}$ and the critical temperature $T_{\rm E}=134(3){\rm MeV}$. Our estimate of the upper bound is derived from the existence of the critical point as an edge of the 1st order phase transition while that of the staggered-type fermions is based on its absence., Comment: 13 pages, 6 figures

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., Comment: 7 pages, 6 figures, presented at the 34th annual International Symposium on Lattice Field Theory (Lattice2016), 24-30 July 2016, University of Southampton, UK

The use of mass preconditioning or Hasenbusch filtering in modern Hybrid Monte Carlo simulations is common. At light quark masses, multiple filters (three or more) are typically used to reduce the cost of generating dynamical gauge fields; however, the task of tuning a large number of Hasenbusch mass terms is non-trivial. The use of short polynomial approximations to the inverse has been shown to provide an effective UV filter for HMC simulations. In this work we investigate the application of polynomial filtering to the mass preconditioned Hybrid Monte Carlo algorithm as a means of introducing many time scales into the molecular dynamics integration with a simplified parameter tuning process. A generalized multi-scale integration scheme that permits arbitrary step- sizes and can be applied to Omelyan-style integrators is also introduced. We find that polynomial-filtered mass-preconditioning (PF-MP) performs as well as or better than standard mass preconditioning, with significantly less fine tuning required., Comment: 33 pages, 16 figures, 8 tables; updated to reflect accepted manuscript

The extraction of hadron matrix elements in lattice QCD using the standard two- and three-point correlator functions demands careful attention to systematic uncertainties. One of the most commonly studied sources of systematic error is contamination from excited states. We apply the variational method to calculate the axial vector current $g_{A}$, the scalar current $g_{S}$ and the quark momentum fraction $\left$ of the nucleon and we compare the results to the more commonly used summation and two-exponential fit methods. The results demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements., Comment: 14 pages, 21 figures

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}$., Comment: 9 pages, 5 figures

The extraction of hadron form factors in lattice QCD using the standard two- and three-point correlator functions has its limitations. One of the most commonly studied sources of systematic error is excited state contamination, which occurs when correlators are contaminated with results from higher energy excitations. We apply the variational method to calculate the axial vector current gA and compare the results to the more commonly used summation and two-exponential fit methods. The results demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements., Comment: 7 pages, 6 figures, talk presented at Lattice 2015, PoS (LATTICE2015)

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., Comment: 7 pages, 5 figures, The 33rd International Symposium on Lattice Field Theory 14-18 July 2015, Kobe ,Japan

We investigate the phase structure of three-flavor QCD in the presence of finite quark chemical potential $\mu/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 $\mu$-$m_{\pi}$ 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., Comment: 18 pages, 8 figures

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., Comment: 7 pages, 4 figures, and 3 tables, presented at the 32nd International Symposium on Lattice Field Theory (LATTICE2014), 23-28 June, 2014, Columbia University New York, NY

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_\pi^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., Comment: 24 pages, 14 figures

We investigate the phase structure of 3-flavor QCD in the presence of finite quark chemical potential $a\mu=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., Comment: 7 pages, 5 figures, 32nd International Symposium on Lattice Field Theory LATTICE 2014, June 23 - 28, 2014, New York, USA

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}\mu, \text{Im}\mu)$ and $(\text{Re}\mu, \text{Im}\beta)$, and present statistically robust zeros of the partition function., Comment: 7 pages, 1 table, and 8 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany. Changes include: fix a typo and clarify the caption for Figure 1; clearly state estimating error for the number of independent configurations

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 $\mu$-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 $\mu$-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 $(\beta, \kappa, \mu/T)=(1.58, 0.1385, 0.584\pm 0.008)$ where $(m_\pi/m_\rho, T/m_\rho)=(0.822, 0.154)$. It weakens considerably at $(\beta, \kappa, \mu/T)=(1.60, 0.1371, 0.821\pm 0.008)$ where $(m_\pi/m_\rho, T/m_\rho)=(0.839, 0.150)$, and a crossover rather than a first order phase transition cannot be ruled out., Comment: 35 pages

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., Comment: 7 pages, 5 figures, Proceedings of the 30th International Symposium on Lattice Field Theory, June 24 - 29, 2012, Cairns, Australia

We compute the strangeness and light-quark contributions Delta s, Delta u and Delta d to the proton spin in nf=2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a approx 0.073 fm, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the renormalization of these matrix elements which involves mixing between contributions from different quark flavours. Our main result is the small negative value MSbar scheme value Delta s = -0.020(10)(4), at a scale mu = 7.4 GeV, of the strangeness contribution to the nucleon spin. The second error is an estimate of the uncertainty, due to the missing extrapolation to the physical point., Comment: 5 pages, 2 figures, V2: minor changes, to enhance the readability of the article. References updated. "The" removed from the title

We determine the strangeness and light quark fractions of the nucleon mass by computing the quark line connected and disconnected contributions to the matrix elements m_q in lattice QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a approx 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized values \sigma_{\pi N} = 38(12) MeV at the physical point and f_{T_s}=\sigma_s/m_N= 0.012(14)^{+10}_{-3} for the strangeness contribution at our larger than physical sea quark mass., Comment: 11 pages, 7 figures, v2: 1 reference added, some corrections

We present an update of an on-going project to determine the disconnected contributions to hadronic structure, specifically, the scalar matrix element and the quark contribution to the spin of the nucleon., Comment: 7 pages, 7 figures, talk given at 28th International Symposium on Lattice Field Theory (Lattice 2010), Villasimius, Italy, 14-19 June 2010

We publish BQCD as free software under the GNU General Public License. BQCD is a Hybrid Monte-Carlo program that simulates lattice QCD with dynamical Wilson fermions. It is one of the main production programs of the QCDSF collaboration. The program can simulate 2 and 2 + 1 fermion flavours with pure, clover improved, and stout smeared fat link Wilson fermions as well as standard plaquette, and an improved (rectangle) gauge action. The single flavour is simulated with the Rational Hybrid Monte-Carlo algorithm., Comment: 7 pages, 1 figure, poster presented at the XXVIII International Symposium on Lattice Field Theory (Lattice 2010), June 14-19 2010, Villasimius, Italy

We compute moments of distribution amplitudes using gauge configurations with two flavors of clover fermions from QCDSF/DIK and operators which are optimized with respect to their behavior under the lattice symmetries. The knowledge of these quantities helps in understanding the internal structure of hadrons and in the analysis of (semi-)exclusive processes. We present results for the nucleon distribution amplitude which suggest that the asymmetries (the deviations from the asymptotic form) are smaller than indicated by sum rule calculations. Using the same approach we were also able to calculate the same quantities for the $N^\star(1535)$, the parity partner of the nucleon. These results show a stronger deviation from the asymptotic form., Comment: Presented at 26th International Symposium on Lattice Field Theory, Williamsburg, Virginia, USA. 8 pages, 2 tables, 2 figures

In recent years the investigation of hadron structure using lattice techniques has attracted growing attention. In this talk we give an overview on recent work with a focus on results for nucleon spectrum and structure from the QCDSF collaboration., Comment: Talk at the "Advanced Studies Institute on SYMMETRIES AND SPIN (SPIN-Praha-2007)", Prague, July 8 - July 14, 2007

We calculate low moments of the leading-twist and next-to-leading twist nucleon distribution amplitudes on the lattice using two flavors of clover fermions. The results are presented in the MSbar scheme at a scale of 2 GeV and can be immediately applied in phenomenological studies. We find that the deviation of the leading-twist nucleon distribution amplitude from its asymptotic form is less pronounced than sometimes claimed in the literature., Comment: 5 pages, 3 figures, 2 tables. RevTeX style. Normalization for \lambda_i corrected. Discussion of the results extended. To be published in PRL

The theta term that breaks the Strong CP symmetry is introduced in the two flavors of dynamical QCD simulation. theta is analytically continued to a pure imaginary number to make the probability of Monte Carlo positive. The Neutron's Electric Dipole Moment (NEDM) is measured on the ensemble under a uniform and week electric field. Other applications of theta terms are also discussed., Comment: 8 pages, 3 figures, talk given by Taku Izubuchi at the XXV International Symposium on Lattice Field Theory, Regensburg, Germany, 30 July - 4 August 2007

We report about on-going simulations of $N_f=2+1$ lattice QCD. We use a tadpole improved Symanzik gauge action and stout link smeared Wilson fermions with a clover term. We employ the Hasenbusch trick for the degenerate u- and d-quarks, and the RHMC algorithm for the simulation of the strange quark., Comment: 6 pages, Contribution to Lattice 2007, Regensburg, Germany, 30 July - 4 August 2007

The electromagnetic form factors provide important insight into the internal structure of the nucleon and continue to be of major interest for experiment and phenomenology. For an intermediate range of momenta the form factors can be calculated on the lattice. However, the reliability of the results is limited by systematic errors mostly due to the required extrapolation to physical quark masses. Chiral effective field theories predict a rather strong quark mass dependence in a range which was yet inaccessible for lattice simulations. We give an update on recent results from the QCDSF collaboration using gauge configurations with dynamical Nf=2, non-perturbatively O(a)-improved Wilson fermions at pion masses as low as 350 MeV., Comment: 11th International Baryons Conference, Seoul, Korea; 7 pages and 4 figures; added references