BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Yoshifumi Nakamura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
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öbius 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.
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_π\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.
We present the calculation of the $K_{l3}$ form factors with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on a large volume of (10.9 fm)$^3$ at one lattice spacing of $a = 0.085$ fm. We extract the form factors from 3-point functions with three different time separations between the source and sink operators to confirm suppression of excited state contributions. The form factors are calculated in very close to the zero momentum transfer, $q^2 = 0$, thanks to the large volume, so that stable interpolations to $q^2 = 0$ are carried out. Using our form factors, we obtain the form factor at $q^2 = 0$, $f_+(0) = 0.9603(16)(^{+14}_{\ -4})(44)(19)(1)$, where the first, second, and fifth errors are statistical, systematic errors from fit functions and the isospin breaking effect, respectively. The third and fourth errors denote the finite lattice spacing effects estimated from the renormalization factor and contribution beyond the leading order SU(3) chiral perturbation theory (ChPT). The result of $f_+(0)$ yields the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, $|V_{us}| = 0.2255(13)(4)$, where the first error comes from our calculation and the second from the experiment. This value is consistent with the ones determined from the unitarity of the CKM matrix and the $K_{l2}$ decay within one standard deviation, while it is slightly larger than recent lattice calculations by at most 1.5 $σ$. Furthermore, we evaluate the shape of the form factors and the phase space integral from our results. We confirm that those results are consistent with the experiment, and also $|V_{us}|$ determined with our phase space integral agrees with the one in the above.
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.
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.
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.
We study the two-dimensional complex $ϕ^{4}$ theory at finite chemical potential using the tensor renormalization group. This model exhibits the Silver Blaze phenomenon in which bulk observables are independent of the chemical potential below the critical point. Since it is expected to be a direct outcome of an imaginary part of the action, an approach free from the sign problem is needed. We study this model systematically changing the chemical potential in order to check the applicability of the tensor renormalization group to the model in which scalar fields are discretized by the Gaussian quadrature. The Silver Blaze phenomenon is successfully confirmed on the extremely large volume $V=1024^2$ and the results are also ensured by another tensor network representation with a character expansion.
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}$.
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<x\right>$ 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.
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.
We investigate the phase structure of three-flavor QCD in the presence of finite quark chemical potential $μ/T\lesssim1.2$ by using the non-perturbatively $O(a)$ improved Wilson fermion action on lattices with a fixed temporal extent $N_{\rm t}=6$ and varied spatial linear extents $N_{\rm s}=8,10,12$. Especially, we focus on locating the critical end point that characterizes the phase structure, and extracting the curvature of the critical line on the $μ$-$m_π$ plane. For Wilson-type fermions, the correspondence between bare parameters and physical parameters is indirect. Hence we present a strategy to transfer the bare parameter phase structure to the physical one, in order to obtain the curvature. Our conclusion is that the curvature is positive. This implies that, if one starts from a quark mass in the region of crossover at zero chemical potential, one would encounter a first-order phase transition when one raises the chemical potential.
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.
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.
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.
We investigate the phase structure of 3-flavor QCD in the presence of finite quark chemical potential $aμ=0.1$ by using the Wilson-Clover fermion action. Especially, we focus on locating the critical end point that characterizes the phase structure. We do this by the kurtosis intersection method for the quark condensate. For Wilson-type fermions, the correspondence between bare parameters and physical parameters is indirect. Hence we present a strategy to transfer the bare parameter phase structure to the physical one.
We investigate the critical endpoint of finite temperature phase transition of $N_f=3$ QCD at zero chemical potential. We employ the renormalization-group improved Iwasaki gauge action and non-perturbatively $O(a)$-improved Wilson-clover fermion action. The critical endpoint is determined by using the intersection point of kurtosis for the temporal size $N_t$=4, 6, 8. Spatial sizes of $N_l$=6-16 ($N_t$=4), 10-24 ($N_t$=6), and 12-24 ($N_t$=8) are employed. We find that $N_t$=4 is out of the scaling region. Using results for $N_t$=6 and 8, and making linear extrapolations in $1/N_t^2$, we obtain $\sqrt{t_0}T_{\rm E}=0.0975(14)(8)$, $\sqrt{t_0}m_{\rm PS,E}=0.2254(52)(105)$ and $m_{\rm PS,E}/T_{\rm E}=2.311(63)(13)$, where the first error is statistical error, the second error is systematic error, and $m_{\rm PS}$ is the pseudo scalar meson mass. If one uses $1/\sqrt{t_0}=1.347(30)$ GeV reported by Borsanyi et al., one finds $T_{\rm E}=131(2)(1)(3)$ MeV, $m_{\rm PS,E}=304(7)(14)(7)$ MeV and $m_{\rm PS,E}/m_{\rm PS,E}^{\rm phys, sym}=0.739(17)(34)(17)$, where the third error comes from error of $\sqrt{t_0}$ and $m_{\rm PS}^{\rm phys, sym}=\sqrt{(m_π^2+2m_K^2)/3}$. Our current estimation of $\sqrt{t_0}m_{\rm PS,E}$ in the continuum limit is about 25% smaller than the SU(3) symmetric point.
We explore the phase space spanned by the temperature and the chemical potential for 4-flavor lattice QCD using the Wilson-clover quark action. In order to determine the order of the phase transition, we apply finite size scaling analyses to gluonic and quark observables including plaquette, Polyakov loop and quark number density, and examine their susceptibility, skewness, kurtosis and Challa-Landau-Binder cumulant. Simulations were carried out on lattices of a temporal size fixed at $N_{\text{t}}=4$ and spatial sizes chosen from $6^3$ up to $10^3$. Configurations were generated using the phase reweighting approach, while the value of the phase of the quark determinant were carefully monitored. The $μ$-parameter reweighting technique is employed to precisely locate the point of the phase transition. Among various approximation schemes for calculating the ratio of quark determinants needed for $μ$-reweighting, we found the Taylor expansion of the logarithm of the quark determinant to be the most reliable. Our finite-size analyses show that the transition is first order at $(β, κ, μ/T)=(1.58, 0.1385, 0.584\pm 0.008)$ where $(m_π/m_ρ, T/m_ρ)=(0.822, 0.154)$. It weakens considerably at $(β, κ, μ/T)=(1.60, 0.1371, 0.821\pm 0.008)$ where $(m_π/m_ρ, T/m_ρ)=(0.839, 0.150)$, and a crossover rather than a first order phase transition cannot be ruled out.
Partition function zeros steer the critical behavior of a system. Studying four-flavor lattice QCD at finite temperature and density with the Wilson-clover fermion action and the Iwasaki gauge action using a phase-quenched fermion determinant, we combine statistics from multiple chemical potentials to improve sampling of the configuration space, and aim at unraveling the movement of zeros in finite systems. Preparing for further investigations, we discuss methods and criteria used to sieve through complex parameter space spanned by $(\text{Re}μ, \text{Im}μ)$ and $(\text{Re}μ, \text{Im}β)$, and present statistically robust zeros of the partition function.
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.
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.
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 <N|qbar q|N> 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 σ_{πN} = 38(12) MeV at the physical point and f_{T_s}=σ_s/m_N= 0.012(14)^{+10}_{-3} for the strangeness contribution at our larger than physical sea quark mass.
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.
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.
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.
The type of the vacuum is studied numerically in the maximally Abelian (MA) gauge and in the Landau (LA) gauge of SU(2) gluodynamics. The type of the vacuum is determined by a ratio between the dual coherence and the dual penetration lengths. The dual penetration length is determined from correlations between Wilson loops and electric fields in both gauges. The dual coherence length is found from correlations between Wilson loops and dimension-2 operators both in the MA and the LA gauges. This determination of the coherence length is supported by theoretical and numerical observation that the dimension-2 gluon operators in the studied gauges have a strong correlation with the monopole current determined in the MA gauge. We find numerically that the dual penetration lengths and the dual coherence lengths in the LA and the MA gauges are almost the same. Therefore we conclude, that in both gauges the type of the vacuum in the confinement phase is near to the border between the type 1 and the type 2 dual superconductors.