Source author record

T. Kaneko

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

128works
9topics
4close collaborators

Published work

128 publications

preprint2020arXiv

FLAG Review 2019

We review lattice results related to pion, kaon, DD-meson, BB-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor f+(0)f_+(0) arising in the semileptonic K→πK \to π transition at zero momentum transfer, as well as the decay constant ratio fK/fπf_K/f_π and its consequences for the CKM matrix elements VusV_{us} and VudV_{ud}. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)L×SU(2)RSU(2)_L\times SU(2)_R and SU(3)L×SU(3)RSU(3)_L\times SU(3)_R Chiral Perturbation Theory. We review the determination of the BKB_K parameter of neutral kaon mixing as well as the additional four BB parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for mcm_c and mbm_b as well as those for DD- and BB-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant αsα_s. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.

preprint2020arXiv

Study of intermediate states in the inclusive semileptonic B→XclνB \rightarrow X_c l ν decay structure function

We analyze the inclusive semileptonic B→XcℓνB \to X_c \ellν structure functions in 2+1-flavor lattice QCD. The Möbius domain-wall fermion action is used for light, strange, charm and bottom quarks. The structure function receives contributions from various exclusive modes, including the dominant S-wave states Ds(∗)D^{(*)}_s as well as the P-wave states Ds∗∗D_s^{**}. We can identify them in the lattice data, from which we put some constraints on the Bs→Ds∗∗ℓνB_s \to D_s^{**}\ellν form factors.

preprint2016arXiv

Chiral behavior of light meson form factors in 2+1 flavor QCD with exact chiral symmetry

We present a study of chiral behavior of light meson form factors in QCD with three flavors of overlap quarks. Gauge ensembles are generated at single lattice spacing 0.12 fm with pion masses down to 300 MeV. The pion and kaon electromagnetic form factors and the kaon semileptonic form factors are precisely calculated using the all-to-all quark propagator. We discuss their chiral behavior using the next-to-next-to-leading order chiral perturbation theory.

preprint2016arXiv

Renormalization of domain-wall bilinear operators with short-distance current correlators

We determine the renormalization constants for flavor non-singlet fermion bilinear operators of Möbius domain-wall fermions. The renormalization condition is imposed on the correlation functions in the coordinate space, such that the non-perturbative lattice calculation reproduces the perturbatively calculated counterpart at short distances. The perturbative expansion is precise as the coefficients are available up to O(αs4)O(α_s^4). We employ 2+12+1-flavor lattice ensembles at three lattice spacings in the range 0.044--0.080~fm.

preprint2016arXiv

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

preprint2016arXiv

Stochastic calculation of the QCD Dirac operator spectrum with Mobius domain-wall fermion

We calculate the spectral function of the QCD Dirac operator using the four-dimensional effective operator constructed from the Mobius domain-wall implementation. We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number. The spectrum in the entire eigenvalue range is obtained with a single set of measurements. Results on 2+1-flavor ensembles with Mobius domain-wall sea quarks at lattice spacing ~ 0.08 fm are shown.

preprint2015arXiv

η′η^\prime meson mass from topological charge density correlator in QCD

The flavor-singlet component of the eta prime meson is related to the topological structure of the SU(3) gauge field through the chiral anomaly. We perform a 2+1-flavor lattice QCD calculation and demonstrate that the two-point function of a gluonically defined topological charge density after a short Yang-Mills gradient flow contains the propagation of the eta prime meson, by showing that its mass in the chiral and continuum limit is consistent with the experimental value. The gluonic correlator does not suffer from the contamination of the pion contribution, and the clean signal is obtained at significantly lower numerical cost compared to the conventional method with the quark bilinear operators.

preprint2015arXiv

Apparent cross-field superslow propagation of magnetohydrodynamic waves in solar plasmas

In this paper we show that the phase mixing of continuum Alfvén waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfvén and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfvén waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfvén/slow frequency across the magnetic surfaces and is slower than the Alfvén/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

preprint2015arXiv

Decay constants and spectroscopy of mesons in lattice QCD using domain-wall fermions

We report results of masses and decay constants of light and charmed pseudo-scalar mesons using lattice QCD with Möbius domain-wall fermions. Using this formulation we are able to compute pseudo-scalar decay constants through the pseudo-scalar density operator as well as with the axial-vector current. Results are shown from several lattice spacings and pion masses between 230 MeV and 500 MeV. We present an analysis of these results at different quark masses to show the chiral properties of the light mesons masses and decay constants.

preprint2015arXiv

Extracting the η′η^\prime meson mass from gluonic correlators in lattice QCD

Calculation of the eta-prime meson mass is a notoriously difficult problem, as it requires evaluation of the disconnected diagram which is costly and noisy. In this work, we use a gluonic operator to extract the eta-prime state after smearing the link variables through the Wilson flow. With this choice, one can avoid a large cancellation of pion contribution between the connected and disconnected diagrams. We obtain the eta-prime meson mass on lattices with three different lattice spacings and two physical volumes, which allow us to estimate its continuum and large volume limits.

preprint2015arXiv

Light meson electromagnetic form factors from three-flavor lattice QCD with exact chiral symmetry

We study the chiral behavior of the electromagnetic (EM) form factors of pion and kaon in three-flavor lattice QCD. In order to make a direct comparison of the lattice data with chiral perturbation theory (ChPT), we employ the overlap quark action that has exact chiral symmetry. Gauge ensembles are generated at a lattice spacing of 0.11 fm with four pion masses ranging between M_pi \simeq 290 MeV and 540 MeV and with a strange quark mass m_s close to its physical value. We utilize the all-to-all quark propagator technique to calculate the EM form factors with high precision. Their dependence on m_s and on the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields, respectively. A detailed comparison with SU(2) and SU(3) ChPT reveals that the next-to-next-to-leading order terms in the chiral expansion are important to describe the chiral behavior of the form factors in the pion mass range studied in this work. We estimate the relevant low-energy constants and the charge radii, and find reasonable agreement with phenomenological and experimental results.

preprint2015arXiv

Nucleon axial and tensor charges with dynamical overlap quarks

We report on our calculation of the nucleon axial and tensor charges in 2+1-flavor QCD with dynamical overlap quarks. Gauge ensembles are generated at a single lattice spacing 0.12 fm and at a strange quark mass close to its physical value. We employ the all-mode-averaging technique to calculate the relevant nucleon correlation functions, and the disconnected quark loop is efficiently calculated by using the all-to-all quark propagator. We present our preliminary results for the isoscalar and isovector charges obtained at pion masses mπm_π = 450 and 540 MeV.

preprint2014arXiv

Computation of the electromagnetic pion form factor from lattice QCD in the epsilon regime

We calculate the electromagnetic pion form factor in lattice QCD with 2+1 flavors of the dynamical overlap quarks. Up and down quark masses are set below their physical values so that the system is in the so-called epsilon regime with the small size of our lattice ~ 1.8 fm. The finite volume corrections are generally expected to be ~ 100% in the epsilon regime. We, however, find a way to automatically cancel the dominant part of them. Inserting non-zero momenta and taking appropriate ratios of the two and three point functions, we can eliminate the contribution from the zero-momentum pion mode. Then the remaining finite volume effect is a small perturbation from the non-zero modes. Our lattice data agree with this theoretical prediction and the extracted pion charge radius is consistent with the experiment.

preprint2014arXiv

Excitonic Bose-Einstein condensation in Ta2NiSe5 above room temperature

We show that finite temperature variational cluster approximation (VCA) calculations on an extended Falicov-Kimball model can reproduce angle-resolved photoemission spectroscopy (ARPES) results on Ta2NiSe5 across a semiconductor-to-semiconductor structural phase transition at 325 K. We demonstrate that the characteristic temperature dependence of the flat-top valence band observed by ARPES is reproduced by the VCA calculation on the realistic model for an excitonic insulator only when the strong excitonic fluctuation is taken into account. The present calculations indicate that Ta2NiSe5 falls in the Bose-Einstein condensation regime of the excitonic insulator state.

preprint2014arXiv

Full O(α)\mathcal{O}(α) electroweak radiative corrections to e+e−→e+e−γe^+e^- \rightarrow e^+e^- γ at the ILC with GRACE-Loop

By using the GRACE-Loop system, we calculate the full O(α)\mathcal{O}(α) electroweak radiative corrections to the process e+e−→e+e−γe^+e^- \rightarrow e^+e^- γ, which is important for future investigations at the International Linear Collider (ILC). With the GRACE-Loop system, the calculations are checked numerically by three consistency tests: ultraviolet finiteness, infrared finiteness, and gauge-parameter independence. The results show good numerical stability when quadruple precision is used. In the phenomenological results, we find that the electroweak corrections to the total cross section range from ∼−4%\sim -4\% to ∼−21%\sim -21\% when s\sqrt{s} varies from 250250 GeV to 11 TeV. The corrections also significantly affect the differential cross sections, which are a function of the invariant masses and angles and the final-particle energies. Such corrections will play an important role for the high-precision program at the ILC.

preprint2014arXiv

Full O(α)\mathcal{O}(α) electroweak radiative corrections to ttˉγt \bar{t} γ and e−e+γe^-e^+γ productions at ILC with GRACE-Loop

The full O(α)\mathcal{O}(α) electroweak radiative corrections to ttˉγt \bar{t} γ and e−e+γe^-e^+γ productions at the International Linear Collider (ILC) are presented in this paper. The computation is performed with the help of GRACE-Loop system. In the physical results, we discuss on the cross section, electroweak corrections, and the top quark forward-backward asymmetry (AFBA_{FB}) which are the function of the center-of-mass energy.

preprint2014arXiv

Order, criticality and excitations in the extended Falicov-Kimball model

Using exact numerical techniques we investigate the nature of excitonic (electron-hole) bound states and the development of exciton coherence in the one-dimensional half-filled extended Falicov-Kimball model. The ground-state phase diagram of the model exhibits, besides band insulator and staggered orbital ordered phases, an excitonic insulator (EI) with power-law correlations. The criticality of the EI state shows up in the von Neumann entropy. The anomalous spectral function and condensation amplitude provide the binding energy and coherence length of the electron-hole pairs which, on their part, point towards a Coulomb interaction driven crossover from BCS-like electron-hole pairing fluctuations to tightly bound excitons. We show that while a mass imbalance between electrons and holes does not affect the location of the BCS-BEC crossover regime it favors staggered orbital ordering to the disadvantage of the EI. Within the BEC regime the quasiparticle dispersion develops a flat valence-band top in accord with the experimental finding for Ta2_2NiSe5_5.

preprint2014arXiv

Simulation Study of Solar Plasma Eruptions Caused by Interactions between Emerging Flux and Coronal Arcade Fields

We investigate the triggering mechanisms of plasma eruptions in the solar atmosphere due to interactions between emerging flux and coronal arcade fields by using two-dimensional MHD simulations. We perform parameter surveys with respect to arcade field height, magnetic field strength, and emerging flux location. Our results show that two possible mechanisms exist, and which mechanism is dominant depends mostly on emerging flux location. One mechanism appears when the location of emerging flux is close to the polarity inversion line (PIL) of an arcade field. This mechanism requires reconnection between the emerging flux and the arcade field, as pointed out by previous studies. The other mechanism appears when the location of emerging flux is around the edge of an arcade field. This mechanism does not require reconnection between the emerging flux and the arcade field but does demand reconnection in the arcade field above the PIL. Furthermore, we found that the eruptive condition for this mechanism can be represented by a simple formula.

preprint2014arXiv

Strong coupling constant from vacuum polarization functions in three-flavor lattice QCD with dynamical overlap fermions

We determine the strong coupling constant αsα_s from a lattice calculation of vacuum polarization functions (VPF) in three-flavor QCD with dynamical overlap fermions. Fitting lattice data of VPF to the continuum perturbative formula including the operator product expansion, we extract the QCD scale parameter ΛMS‾(3)Λ_{\overline{MS}}^{(3)}. At the ZZ boson mass scale, we obtain αs(5)(MZ)=0.1181(3)(−12+14)α_s^{(5)}(M_Z)=0.1181(3)(^{+14}_{-12}), where the first error is statistical and the second is our estimate of various systematic uncertainties.

preprint2014arXiv

Topology density correlator on dynamical domain-wall ensembles with nearly frozen topological charge

Global topological charge decorrelates very slowly or even freezes in fine lattice simulations. On the other hand, its local fluctuations are expected to survive and lead to the correct physical results as long as the volume is large enough. We investigate this issue on recently generated configurations including dynamical domain-wall fermions at lattice spacings a = 0.08 fm and finer. We utilize the Yang-Mills gradient flow to define the topological charge density operator and calculate its long-distance correlation, through which we propose a new method for extracting the topological susceptibility in a sub-volume. This method takes care of the finite volume correction, which reduces the bias caused by the global topological charge. Our lattice data clearly show a shorter auto-correlation time than that of the naive definition using the whole lattice, and are less sensitive to the global topological history. Numerical results show a clear sea-quark mass dependence, which agrees well with the prediction of chiral perturbation theory.

preprint2013arXiv

Exact-diagonalization study of exciton condensation in electron bilayers

We report on small-cluster exact-diagonalization calculations which prove the formation of electron-hole pairs (excitons) as prerequisite for spontaneous interlayer phase coherence in bilayer systems described by the extended Falicov-Kimball model. Evaluating the anomalous Green's function and momentum distribution function of the pairs, and thereby analyzing the dependence of the exciton binding energy, condensation amplitude, and coherence length on the Coulomb interaction strength, we demonstrate a crossover between a BCS-like electron-hole pairing transition and a Bose-Einstein condensation of tightly bound preformed excitons. We furthermore show that a mass imbalance between electrons and holes tends to suppress the condensation of excitons.

preprint2013arXiv

Full O(α)\mathcal{O}(α) electroweak radiative corrections to e+e−→ttˉγe^+e^- \rightarrow t \bar{t} γ with GRACE-Loop

We present the full O(α)\mathcal{O}(α) electroweak radiative corrections to the process e+e−→ttˉγe^+e^- \rightarrow t \bar{t} γ at the International Linear Collider (ILC). The computation is performed with the help of the GRACE-Loop system. We present the total cross-section and the top quark forward-backward asymmetry (AFBA_{FB}) as a function of the center-of-mass energy and compare them with the process e+e−→ttˉe^+e^- \rightarrow t \bar{t}. We find that the value of AFBA_{FB} in ttˉγt \bar{t} γ production is larger than AFBA_{FB} in ttˉt\bar{t} production. It is an important result for the measurement of the top quark forward-backward asymmetry at the ILC. Applying a structure function method, we also subtract the QED correction to gain the genuine weak correction in both the αα scheme and the GμG_μ scheme (δWGμδ_{W}^{G_μ}). We obtain numerical values for δWGμδ_{W}^{G_μ} which are changing from 2% to -24% when we vary the center-of-mass energy from 360 GeV to 1 TeV.

preprint2013arXiv

Large-scale simulations with chiral symmetry

We carry out a comparative study among five-dimensional formulations of chirally symmetric fermions about the algorithmic performance, chiral symmetry violation and topological tunneling to find a computationally inexpensive formulation with good chiral symmetry. With our choice of the lattice action, we have launched large-scale simulations on fine lattices aiming at a precision study of light and heavy quark physics. We report on the comparative study, current status of the large-scale simulations, and preliminary results on the residual quark mass and auto-correlation.

preprint2013arXiv

Nucleon strange quark content from N_f=2+1 lattice QCD with exact chiral symmetry

We calculate the strange quark content of the nucleon in 2+1-flavor lattice QCD. Chirally symmetric overlap fermion formulation is used to avoid the contamination from up and down quark contents due to an operator mixing between strange and light scalar operators, \bar{s}s and \bar{u}u+\bar{d}d. At a lattice spacing a=0.112(1) fm, we perform calculations at four values of degenerate up and down quark masses, which cover a range of the pion mass M_pi \simeq 300-540 MeV. We employ two different methods: one is a direct method where we calculate the strange quark content by directly inserting the strange scalar operator. The other is an indirect method where the quark content is extracted from a derivative of the nucleon mass in terms of the strange quark mass. With these two methods we obtain consistent results with each other. Our best estimate f_{T_s}=0.009(15)(16) is in good agreement with our previous studies in two-flavor QCD.

preprint2012arXiv

A blind detection of a large, complex, Sunyaev--Zel'dovich structure

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100μJy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological β-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 μK - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

preprint2012arXiv

Chiral behavior of kaon semileptonic form factors in lattice QCD with exact chiral symmetry

We calculate the kaon semileptonic form factors in lattice QCD with three flavors of dynamical overlap quarks. Gauge ensembles are generated at pion masses as low as 290 MeV and at a strange quark mass near its physical value. We precisely calculate relevant meson correlators using the all-to-all quark propagator. Twisted boundary conditions and the reweighting technique are employed to vary the momentum transfer and the strange quark mass. We discuss the chiral behavior of the form factors by comparing with chiral perturbation theory and experiments.

preprint2012arXiv

Intrinsic giant Stark effect of boron-carbon-nitride nanoribbons with zigzag edges

Electronic properties of zigzag boron-carbon-nitride (BCN) nanoribbons, where the outermost C atoms on the edges of graphene nanoribbons are replaced by B or N atoms, are theoretically studied using the first-principles calculations. We show that BCN nanoribbons are metallic, since several bands cross the Fermi level. For BCN nanoribbons in a rich H2_2 environment, the so-called nearly free electron state appears just above the Fermi level because of the intrinsic giant Stark effect due to the internal electric field of a transverse dipole moment. The position of the nearly free electron state can be controlled by applying an electric field parallel to the dipole moment. The hydrogenation of the nitrogen atom is necessary for the appearance of the giant Stark effect in BCN nanoribbons. We also discuss the effect of stacking order on the intrinsic giant Stark effect in bilayer BCN nanoribbons.

preprint2012arXiv

Orthorhombic-to-Monoclinic Phase Transition of Ta2NiSe5 Induced by the Bose-Einstein Condensation of Excitons

Using the band structure calculation and mean-field analysis of the derived three-chain Hubbard model with phonon degrees of freedom, we discuss the origin of the orthorhombic-to-monoclinic phase transition of the layered chalcogenide Ta2_2NiSe5_5. We show that the Bose-Einstein condensation of excitonic electron-hole pairs cooperatively induces the instability of the phonon mode at momentum q→0q\rightarrow 0 in the quasi-one-dimensional Ta-NiSe-Ta chain, resulting in the structural phase transition of the system. The calculated single-particle spectra reproduce the deformation of the band structure observed in the angle-resolved photoemission spectroscopy experiment.

preprint2012arXiv

Spin wave modes in magnetic nanodisks under in-plane magnetic field

The size dependence of spin wave modes in a circular Permalloy (Py) nanodisk under an in-plane magnetic field is systematically studied by using micromagnetics simulations. We show that as the disk diameter is increased, the resonance frequency of the backward mode deceases while that of the uniform mode increases. The avoided crossing of resonance frequencies of the uniform mode and the backward mode appears in the plot of the size dependence of resonance frequencies and the backward mode turns into the so-called "edge mode" for large nanodisks.

preprint2011arXiv

Chiral interpolation in a finite volume

A simulation of lattice QCD at (or even below) the physical pion mass is feasible on a small lattice size of \sim 2 fm. The results are, however, subject to large finite volume effects. In order to precisely understand the chiral behavior in a finite volume, we develop a new computational scheme to interpolate the conventional epsilon and p regimes within chiral perturbation theory. In this new scheme, we calculate the two-point function in the pseudoscalar channel, which is described by a set of Bessel functions in an infra-red finite way as in the epsilon regime, while chiral logarithmic effects are kept manifest as in the p regime. The new ChPT formula is compared to our 2+1- flavor lattice QCD data near the physical up and down quark mass, mud \sim 3 MeV on an L \sim 1.8 fm lattice. We extract the pion mass = 99(4) MeV, from which we attempt a chiral "interpolation" of the observables to the physical point.

preprint2011arXiv

Determination of the chiral condensate from QCD Dirac spectrum on the lattice

We calculate the chiral condensate of QCD with 2, 2+1 and 3 flavors of sea quarks. Lattice QCD simulations are performed employing dynamical overlap fermions with up and down quark masses covering a range between 3 and 100 MeV. On L ~ 1.8-1.9 fm lattices at a lattice spacing ~ 0.11 fm, we calculate the eigenvalue spectrum of the overlap-Dirac operator. By matching the lattice data with the analytical prediction from chiral perturbation theory at the next-to-leading order, the chiral condensate in the massless limit of up and down quarks is determined.

preprint2011arXiv

Kaon semileptonic form factors in QCD with exact chiral symmetry

We report on our calculation of the kaon semileptonic form factors in Nf=2+1 lattice QCD. Chiral symmetry is exactly preserved by using the overlap quark action for a straightforward comparison with chiral perturbation theory (ChPT). We simulate three pion masses down to 290 MeV at a single lattice spacing of 0.11 fm and at a strange quark mass very close to its physical value. The form factors near zero momentum transfer are precisely calculated by using the all-to-all propagator and twisted boundary conditions. We compare the normalizations and slopes of the form factors with ChPT and experiments.

preprint2011arXiv

Nucleon strange quark content from two-flavor lattice QCD with exact chiral symmetry

Strange quark content of the nucleon is calculated in dynamical lattice QCD employing the overlap fermion formulation. For this quantity, exact chiral symmetry guaranteed by the Ginsparg-Wilson relation is crucial to avoid large contamination due to a possible operator mixing with uˉu+dˉd\bar{u}u+\bar{d}d. Gauge configurations are generated with two dynamical flavors on a 16^3 x 32 lattice at a lattice spacing a \simeq 0.12fm. We directly calculate the relevant three-point function on the lattice including a disconnected strange quark loop utilizing the techniques of all-to-all quark propagator and low-mode averaging. Our result f_{T_s} = 0.032(8)(22), is in good agreement with our previous indirect estimate using the Feynman-Hellmann theorem.

preprint2010arXiv

Determination of the chiral condensate from 2+1-flavor lattice QCD

We perform a precise calculation of the chiral condensate in QCD using lattice QCD with 2+1 flavors of dynamical overlap quarks. Up and down quark masses cover a range between 3 and 100 MeV on a 16^3x48 lattice at a lattice spacing around 0.11 fm. At the lightest sea quark mass, the finite volume system on the lattice is in the epsilon-regime. By matching the low-lying eigenvalue spectrum of the Dirac operator with the prediction of chiral perturbation theory at the next-to-leading order, we determine the chiral condensate in 2+1-flavor QCD with strange quark mass fixed at its physical value as Sigma (MS-bar at 2 GeV) = [242(04)(^+19_-18}) MeV}]^3, where the errors are statistical and systematic, respectively.

preprint2010arXiv

Light meson form factors in N_f=2+1 QCD with dynamical overlap quarks

We report on our calculation of pion and kaon form factors in three-flavor QCD using the overlap quark action. Gauge ensembles are generated on a 16^3 \times 48 lattice at a lattice spacing of 0.11 fm with pion masses down to 310 MeV. Connected and disconnected meson correaltors are calculated using the all-to-all quark propagator. We present our preliminary analysis on the chiral behavior of the electromagnetic and scalar form factors as well as a comparison of the shape of the K to pi form factors with experiment.

preprint2010arXiv

Nucleon strange quark content in 2+1-flavor QCD

We calculate the strange quark content of the nucleon directly from its disconnected three-point function in N_f=2+1 QCD. Chiral symmetry is crucial to avoid a possibly large contamination due to operator mixing, and is exactly preserved by employing the overlap quark action. We also use the all-to-all quark propagator and the low-mode averaging technique in order to accurately calculate the relevant nucleon correlator. Our preliminary result extrapolated to the physical point is f_{T_s} = 0.013(12)(16). This is in good agreement with our previous estimate in N_f=2 QCD as well as those from our indirect calculations using the Feynman-Hellmann theorem.

preprint2009arXiv

Calculation of nucleon strange quark content with dynamical overlap quarks

We calculate the nucleon strange quark content directly from disconnected three-point functions. Numerical simulations are carried out in two-flavor QCD using the overlap quark action with up and down quark masses down to a fifth of the physical strange quark mass. To improve the statistical accuracy, we calculate the nucleon two-point functions with the low-mode averaging technique, whereas the all-to-all quark propagator is used for the disconnected quark loop. We obtain the y parameter, which is the ratio of the strange and light quark contents, y = 0.024(45) at the physical point. This is in a good agreement with our earlier calculation from the nucleon spectrum through the Feynman-Hellmann theorem.

preprint2009arXiv

Chiral properties of light mesons with Nf=2+1N_f=2+1 overlap fermions

We present an update of the light meson spectrum with NfN_f=2+1 overlap fermions on a 163×4816^3\times 48 lattice at five different up and down quark masses and two strange quark masses. Based on our experience with the previous simulation with Nf=2N_f=2, we carry out the chiral extrapolation with the prediction of the chiral perturbation theory at the next-to-next-to leading order. We also check the consistency of our analysis by using alternative chiral extrapolation with a reduced theory in which the strange quark mass is integrated out.

preprint2009arXiv

Flavor-singlet mesons in N_f=2+1 QCD with dynamical overlap quarks

We report on our study of flavor-singlet mesons in three-flavor QCD with dynamical overlap quarks. Gauge ensembles are generated on a 16^3 x 48 lattice at a lattice spacing of 0.10 fm with the strange quark masses around its physical value m_s,phys and up and down quark masses down to m_s,phys/5. Connected and disconnected meson correlators are calculated using the all-to-all quark propagator. We present our preliminary results on the spectrum of flavor-singlet pseudoscalar and vector mesons.

preprint2009arXiv

Nucleon sigma term and strange quark content in 2+1-flavor QCD with dynamical overlap fermions

We study the sigma term and the strange quark content of nucleon in 2+1-flavor QCD with dynamical overlap fermions. We analyze the lattice data of nucleon mass taken at two different strange quark masses with five values of up and down quark masses each. Using the reweighting technique, we study the strange quark mass dependence of the nucleon and extract the strange quark content.

preprint2009arXiv

Pion form factors from two-flavor lattice QCD with exact chiral symmetry

We calculate pion vector and scalar form factors in two-flavor lattice QCD and study the chiral behavior of the vector and scalar radii <r^2>_{V,S}. Numerical simulations are carried out on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with quark masses down to \sim m_s/6, where m_s is the physical strange quark mass. Chiral symmetry, which is essential for a direct comparison with chiral perturbation theory (ChPT), is exactly preserved in our calculation at finite lattice spacing by employing the overlap quark action. We utilize the so-called all-to-all quark propagator in order to calculate the scalar form factor including the contributions of disconnected diagrams and to improve statistical accuracy of the form factors. A detailed comparison with ChPT reveals that the next-to-next-to-leading-order contributions to the radii are essential to describe their chiral behavior in the region of quark mass from m_s/6 to m_s/2. Chiral extrapolation based on two-loop ChPT yields <r^2>_V=0.409(23)(37)fm and <r^2>_S=0.617(79)(66)fm, which are consistent with phenomenological analysis. We also present our estimates of relevant low-energy constants.

preprint2008arXiv

Light meson spectrum with Nf=2+1N_f=2+1 dynamical overlap fermions

We report on a numerical simulation with 2+1 dynamical flavors of overlap fermions. We calculate pseudo-scalar masses and decay constants on a 163×48×(0.11fm)416^3\times 48 \times (0.11 {\rm fm})^4 lattice at five different up and down quark masses and two strange quark masses. The lightest pion mass corresponds to ≈310\approx 310 MeV. We also study the validity of the chiral perturbation theory using the results of the numerical simulation with two dynamical flavors and conclude that the one-loop formulae cannot be directly applied in the strange quark mass region. We therefore extrapolate our 2+1-flavor results to the chiral limit by fitting the data to the two-loop formulae of the chiral perturbation theory.

preprint2008arXiv

Light quark masses from unquenched lattice QCD

We calculate the light meson spectrum and the light quark masses by lattice QCD simulation, treating all light quarks dynamically and employing the Iwasaki gluon action and the nonperturbatively O(a)-improved Wilson quark action. The calculations are made at the squared lattice spacings at an equal distance a^2~0.005, 0.01 and 0.015 fm^2, and the continuum limit is taken assuming an O(a^2) discretization error. The light meson spectrum is consistent with experiment. The up, down and strange quark masses in the \bar{MS} scheme at 2 GeV are \bar{m}=(m_{u}+m_{d})/2=3.55^{+0.65}_{-0.28} MeV and m_s=90.1^{+17.2}_{-6.1} MeV where the error includes statistical and all systematic errors added in quadrature. These values contain the previous estimates obtained with the dynamical u and d quarks within the error.

preprint2008arXiv

Pion vector and scalar form factors with dynamical overlap quarks

We calculate the pion vector and scalar form factors in two-flavor QCD. Gauge configurations are generated with dynamical overlap quarks on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with sea quark masses down to a sixth of the physical strange quark mass. Contributions of disconnected diagrams to the scalar form factor is calculated employing the all-to-all quark propagators. We present a detailed comparison of the vector and scalar radii with chiral perturbation theory to two loops.

preprint2008arXiv

Topological susceptibility in (2+1)-flavor lattice QCD with overlap fermion

We determine the topological susceptibility χ_t in the topologically-trivial sector generated by lattice simulations of N_f = 2+1 QCD with overlap Dirac fermion, on a 16^3 x 48 lattice with lattice spacing ~ 0.11 fm, for five sea quark masses m_q ranging from m_s/6 to m_s (where m_s is the physical strange quark mass). The χ_t is extracted from the plateau (at large time separation) of the 2-point and 4-point time-correlation functions of the flavor-singlet pseudoscalar meson η', which arises from the finite size effect due to fixed topology. In the small m_q regime, our result of χ_t agrees with the chiral effective theory. Using the formula χ_t = Σ(m_u^{-1} + m_d^{-1} + m_s^{-1})^{-1} by Leutwyler-Smilga, we obtain the chiral condensate Σ^{MSbar}(2 GeV) = [249(4)(2) MeV]^3.

preprint2007arXiv

2+1 flavor light hadron spectrum and quark masses with the O(a) improved Wilson-clover quark formalism

We present a summary of results of the joint CP-PACS and JLQCD project toward a 2+1 flavor full QCD simulation with the O(a)-improved Wilson quark formalism and the Iwasaki gauge action. Configurations were generated during 2002-2005 at three lattice spacings, a~0.076, 0.100 and 0.122 fm, keeping the physical volume constant at (2.0fm)^3. Up and down quark masses are taken in the range m_{PS}/m_V~0.6-0.78. We have completed the analysis for the light meson spectrum and quark masses in the continuum limit using the full configuration set. The predicted meson masses reproduce experimental values in the continuum limit at a 1% level. The average up and down, and strange quark masses turn out to be m_{ud}^{\bar{MS}}(μ=2 GeV)=3.50(14)({}^{+26}_{-15}) MeV and m_s^{\bar{MS}}(μ=2 GeV)=91.8(3.9)({}^{+6.8}_{-4.1}) MeV. We discuss our future strategy toward definitive results on hadron spectroscopy with the Wilson-clover formalism.

preprint2007arXiv

Implementation of an Optimised Cassegrain System for Radio Telescopes

We present the antenna design for a radio interferometer, the Arcminute Microkelvin Imager, together with its beam pattern measurement. Our aim was to develop a low-cost system with high aperture efficiency and low ground-spill across the frequency range 12-18GHz. We use a modified cassegrain system consisting of a commercially-available paraboloidal primary mirror with a diameter of 3.7m, and a shaped secondary mirror. The secondary mirror is oversized with respect to a ray-optics design and has a surface that is bent towards the primary near its outer edge using a square term for the shaping. The antennas are simple to manufacture and therefore their cost is low. The design increased the antenna gain by approximately 10 per cent compared to a normal Cassegrain system while still maintaining low contamination from ground-spill and using a simple design for the horn.

preprint2007arXiv

Non-perturbative improvement of the axial current with three dynamical flavors and the Iwasaki gauge action

We perform a non-perturbative determination of the improvement coefficient c_A to remove O(a) discretization errors in the axial vector current in three-flavor lattice QCD with the Iwasaki gauge action and the standard O(a)(a)-improved Wilson quark action. An improvement condition with a good sensitivity to c_A is imposed at constant physics. Combining our results with the perturbative expansion, c_A is now known rather precisely for 1/a \gtrsim 1.6 GeV.

preprint2007arXiv

Two-flavor lattice QCD in the epsilon-regime and chiral Random Matrix Theory

The low-lying eigenvalue spectrum of the QCD Dirac operator in the epsilon-regime is expected to match with that of chiral Random Matrix Theory (ChRMT). We study this correspondence for the case including sea quarks by performing two-flavor QCD simulations on the lattice. Using the overlap fermion formulation, which preserves exact chiral symmetry at finite lattice spacings, we push the sea quark mass down to \sim 3 MeV on a 16^3\times 32 lattice at a lattice spacing a \simeq 0.11 fm. We compare the low-lying eigenvalue distributions and find a good agreement with the analytical predictions of ChRMT. By matching the lowest-lying eigenvalue we extract the chiral condensate, Σ(2 GeV)[MSbar] = [251(7)(11) MeV]^3, where errors represent statistical and higher order effects in the epsilon expansion. We also calculate the eigenvalue distributions on the lattices with heavier sea quarks at two lattice spacings. Although the epsilon expansion is not applied for those sea quarks, we find a reasonable agreement of the Dirac operator spectrum with ChRMT. The value of Sigma, after extrapolating to the chiral limit, is consistent with the estimate in the epsilon-regime.

preprint2006arXiv

Nonperturbative O(a)O(a) improvement of the Wilson quark action with the RG-improved gauge action using the Schrödinger functional method

We perform a nonperturbative determination of the O(a)O(a)-improvement coefficient cSWc_{\rm SW} and the critical hopping parameter κcκ_c for NfN_f=3, 2, 0 flavor QCD with the RG-improved gauge action using the Schrödinger functional method. In order to interpolate cSWc_{\rm SW} and κcκ_c as a function of the bare coupling, a wide range of ββ from the weak coupling region to the moderately strong coupling points used in large-scale simulations is studied. Corrections at finite lattice size of O(a/L)O(a/L) turned out to be large for the RG-improved gauge action, and hence we make the determination at a size fixed in physical units using a modified improvement condition. This enables us to avoid O(a)O(a) scaling violations which would remain in physical observables if cSWc_{\rm SW} determined for a fixed lattice size L/aL/a is used in numerical simulations.

preprint2005arXiv

I=2 Pion Scattering Length from Two-Pion Wave Functions

We calculate the two-pion wave function in the ground state of the I=2 SS-wave system and find the interaction range between two pions, which allows us to examine the validity of the necessary condition for the finite-volume method for the scattering length proposed by Lüscher. We work in the quenched approximation employing a renormalization group improved gauge action for gluons and an improved Wilson action for quarks at 1/a=1.207(12)GeV1/a=1.207(12) {\rm GeV} on 163×8016^3 \times 80, 203×8020^3 \times 80 and 243×8024^3 \times 80 lattices. We conclude that the necessary condition is satisfied within the statistical errors for the lattice sizes L≥24L\ge 24 (3.92fm3.92 {\rm fm}) when the quark mass is in the range that corresponds to mπ2=0.273−0.736GeV2m_π^2 = 0.273-0.736 {\rm GeV}^2. We obtain the scattering length with a smaller statistical error from the wave function than from the two-pion time correlator.

preprint2005arXiv

Light hadron spectrum and quark masses in 2+1 flavor QCD

CP-PACS and JLQCD collaborations are carrying out a joint project of the 2+1 flavor full QCD simulation. Gauge configurations are generated for the non-perturbatively O(a)O(a)-improved Wilson quark action and the Iwasaki gauge action using PHMC algorithm at three lattice spacings, a∼0.076a\sim 0.076, 0.010 and 0.122 fm, with a fixed physical volume (2.0fm)3(2.0 fm)^3. We present analysis for the light meson spectrum and quark masses in the continuum limit, which are determined using data obtained from the simulations at the two coarser lattices. Our simulations reproduce experimental values of meson masses. The ud and strange quark masses turn out to be mudMSˉ(μ=2GeV)=3.34(23)MeVm_{ud}^{\bar{MS}}(μ=2 GeV)=3.34(23) MeV and msMSˉ(μ=2GeV)=86.7(5.9)MeVm_s^{\bar{MS}}(μ=2 GeV)=86.7(5.9) MeV. We also show preliminary results at our finest lattice spacing for which simulations are still being continued.

preprint2005arXiv

Non-perturbative O(a)-improvement of Wilson quark action in three-flavor QCD with plaquette gauge action

We perform a non-perturbative determination of the O(a)-improvement coefficient c_SW for the Wilson quark action in three-flavor QCD with the plaquette gauge action. Numerical simulations are carried out in a range of β=12.0-5.2 on a single lattice size of 8^3x16 employing the Schrödinger functional setup of lattice QCD. As our main result, we obtain an interpolation formula for c_SW and the critical hopping parameter K_c as a function of the bare coupling. This enables us to remove O(a) scaling violation from physical observables in future numerical simulation in the wide range of β. Our analysis with a perturbatively modified improvement condition for c_SW suggests that finite volume effects in c_SW are not large on the 8^3x16 lattice. We investigate N_f dependence of c_SW by additional simulations for N_f=4, 2 and 0 at β=9.6. As a preparatory step for this study, we also determine c_SW in two-flavor QCD at β=5.2. At this β, several groups carried out large-scale calculations of the hadron spectrum, while no systematic determination of c_SW has been performed.

preprint2004arXiv

A scaling study of the step scaling function in SU(3) gauge theory with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the renormalization-group improved Iwasaki gauge action and the perturbatively improved Lüscher-Weisz gauge action. We confirm that the step scaling functions from the improved gauge actions agree with that previously obtained from the plaquette action within errors in the continuum limit at both weak and strong coupling regions. We also investigate how different choices of boundary counter terms for the improved gauge actions affect the scaling behavior. In the extrapolation to the continuum limit, we observe that the cut off dependence becomes moderate for the Iwasaki action, if a perturbative reduction of scaling violations is applied to the simulation results. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

A scaling study of the step scaling function of quenched QCD with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the Iwasaki gauge action and the Luescher-Weisz gauge action. In particular, we test the choice of boundary counter terms and apply a perturbative procedure for removal of lattice artifacts for the simulation results in the extrapolation procedure. We confirm the universality of the step scaling functions at both weak and strong coupling regions. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

Bulk first-order phase transition in three-flavor lattice QCD with O(a)O(a)-improved Wilson fermion action at zero temperature

Three-flavor QCD simulation with the O(a)O(a)-improved Wilson fermion action is made employing an exact fermion algorithm developed for odd number of quark flavors. For the plaquette gauge action, an unexpected first-order phase transition is found in the strong coupling regime (β≲β\lesssim 5.0) at relatively heavy quark masses (mPS/mV∼m_{\mathrm{PS}}/m_{\mathrm{V}}\sim 0.74--0.87). Strong metastability persists on a large lattice of size 123×3212^3\times 32, which indicates that the transition has a bulk nature. The phase gap becomes smaller toward weaker couplings and vanishes at β≃5.0β\simeq 5.0, which corresponds to a lattice spacing a≃a\simeq 0.1 fm. The phase transition is not found if the improved gauge actions are employed. Our results imply that realistic simulations of QCD with three flavors of dynamical Wilson-type fermions at lattice spacings in the range a=a= 0.1--0.2 fm require use of improved gauge actions. Possible origins of the phase transition is discussed.

preprint2004arXiv

Chiral extrapolations with small sea quark mass data in two-flavor lattice QCD

We present results on the light hadron spectrum and quark mass in two-flavor QCD calculated with small sea quark masses down to mPS/mV=0.35m_{PS}/m_{V}=0.35. The configurations are generated using the RG improved gauge and tadpole-improved clover quark action at β=1.8β=1.8, where a−1≃1a^{-1} \simeq 1 GeV. Chiral extrapolations are made using not only polynomials and ChPT in the continuum but also formulae of Wilson chiral perturbation theory (WChPT) including O(a2)O(a^2) chiral breaking terms. We examine the viability of WChPT and its influence on quark masses.

preprint2004arXiv

I=2 ππππ Scattering Phase Shift with two Flavors of O(a)O(a) Improved Dynamical Quarks

We present a lattice QCD calculation of phase shift including the chiral and continuum extrapolations in two-flavor QCD. The calculation is carried out for I=2 S-wave ππππ scattering. The phase shift is evaluated for two momentum systems, the center of mass and laboratory systems, by using the finite volume method proposed by Lüscher in the center of mass system and its extension to general systems by Rummukainen and Gottlieb. The measurements are made at three different bare couplings β=1.80β= 1.80, 1.95 and 2.10 using a renormalization group improved gauge and a tadpole improved clover fermion action, and employing a set of configurations generated for hadron spectroscopy in our previous work. The illustrative values we obtain for the phase shift in the continuum limit are δδ(deg.) =−3.50(64)= - 3.50(64), −9.5(30) - 9.5(30) and −16.9(64) - 16.9(64) for s(GeV)\sqrt{s}({\rm GeV}) =0.4=0.4, 0.6 0.6 and 0.8 0.8, which are consistent with experiment.

preprint2004arXiv

I=2 Pion Scattering Length from Two-Pion Wave Function

We present a report on a calculation of scattering length for I=2 SS-wave two-pion system from two-pion wave function. Calculations are made with an RG-improved action for gluons and improved Wilson action for quarks at a−1=1.207(12)GeVa^{-1}=1.207(12) {\rm GeV} on 163×8016^3 \times 80, 203×8020^3 \times 80 and 243×8024^3 \times 80 lattices. We investigate the validity of necessary condition for application of Lüscher's formula through the wave function. We find that the condition is satisfied for lattice volumes L≥3.92fmL\ge 3.92 {\rm fm} for the quark mass range mπ2=0.273−0.736GeV2m_π^2 = 0.273-0.736 {\rm GeV}^2. We also find that the scattering length can be extracted with a smaller statistical error from the wave function than with a time correlation function used in previous studies.

preprint2004arXiv

Lattice QCD calculation of the proton decay matrix element in the continuum limit

We present a quenched lattice QCD calculation of the αand βparameters of the proton decay matrix element. The simulation is carried out using the Wilson quark action at three values of the lattice spacing in the range a\approx 0.1-0.064 fm to study the scaling violation effect. We find only mild scaling violation when the lattice scale is determined by the nucleon mass. We obtain in the continuum limit, |α(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |β(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with αand βin a relatively opposite sign, where the first error is statistical and the second is due to the uncertainty in the determination of the physical scale.

preprint2004arXiv

Light hadron spectroscopy in two-flavor QCD with small sea quark masses

We extend the study of the light hadron spectrum and the quark mass in two-flavor QCD to smaller sea quark mass, corresponding to mPS/mV=0.60m_{PS}/m_{V}=0.60--0.35. Numerical simulations are carried out using the RG-improved gauge action and the meanfield-improved clover quark action at β=1.8β=1.8 (a=0.2a = 0.2 fm from ρρ meson mass). We observe that the light hadron spectrum for small sea quark mass does not follow the expectation from chiral extrapolations with quadratic functions made from the region of mPS/mV=0.80m_{PS}/m_{V}=0.80--0.55. Whereas fits with either polynomial or continuum chiral perturbation theory (ChPT) fails, the Wilson ChPT (WChPT) that includes a2a^2 effects associated with explicit chiral symmetry breaking successfully fits the whole data: In particular, WChPT correctly predicts the light quark mass spectrum from simulations for medium heavy quark mass, such as $m_{PS}/m_V \simgt 0.5$. Reanalyzing the previous data %at mPS/mV=0.80m_{PS}/m_{V}=0.80--0.55 with the use of WChPT, we find the mean up and down quark mass being smaller than the previous result from quadratic chiral extrapolation by approximately 10%, mudMSˉ(μ=2GeV)=3.11(17)m_{ud}^{\bar{\rm MS}}(μ=2 {GeV}) = 3.11(17) [MeV] in the continuum limit.

preprint2004arXiv

Non-perturbative renormalization of meson decay constants in quenched QCD for a renormalization group improved gauge action

Renormalization constants (ZZ-factors) of vector and axial-vector currents are determined non-perturbatively in quenched QCD for a renormalization group improved gauge action and a tadpole improved clover quark action using the Schrödinger functional method. Non-perturbative values of ZZ-factors turn out to be smaller than one-loop perturbative values by O(15O(15%) at lattice spacing of a−1≈a^{-1}\approx 1 GeV. The pseudoscalar and vector meson decay constants calculated with the non-perturbative ZZ-factors show a much better scaling behavior compared to previous results obtained with tadpole improved one-loop ZZ-factors. In particular, the non-perturbative ZZ-factors normalized at infinite physical volume show that scaling violation of the decay constants are within about 10% up to the lattice spacing a−1∼1a^{-1}\sim 1 GeV. The continuum estimates obtained from data in the range a−1=a^{-1}= 1 -- 2 GeV agree with those determined from finer lattices (a−1∼2−4a^{-1}\sim 2-4 GeV) with the standard action.

preprint2003arXiv

An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion

We propose an exact simulation algorithm for lattice QCD with dynamical Kogut-Susskind fermion in which the N_f-flavor fermion operator is defined as the N_f/4-th root of the Kogut-Susskind (KS) fermion operator. The algorithm is an extension of the Polynomial Hybrid Monte Carlo (PHMC) algorithm to KS fermions. The fractional power of the KS fermion operator is approximated with a Hermitian Chebyshev polynomial, with which we can construct an algorithm for any number of flavors. The error which arises from the approximation is corrected by the Kennedy-Kuti noisy Metropolis test. Numerical simulations are performed for the two-flavor case for several lattice parameters in order to confirm the validity and the practical feasibility of the algorithm. In particular tests on a 16^4 lattice with a quark mass corresponding to m_{PS}/m_V ~ 0.68 are successfully accomplished. We conclude that our algorithm provides an attractive exact method for dynamical QCD simulations with KS fermions.

preprint2003arXiv

Calculation of Non-Leptonic Kaon Decay Amplitudes from K→πK\toπ Matrix Elements in Quenched Domain-Wall QCD

We explore application of the domain wall fermion formalism of lattice QCD to calculate the K→ππK\toππ decay amplitudes in terms of the K→πK\toπ and K→0K\to 0 hadronic matrix elements through relations derived in chiral perturbation theory. Numerical simulations are carried out in quenched QCD using domain-wall fermion action for quarks and an RG-improved gauge action for gluons on a 163×32×1616^3\times 32\times 16 and 243×32×1624^3\times 32\times 16 lattice at β=2.6β=2.6 corresponding to the lattice spacing 1/a≈21/a\approx 2GeV. Quark loop contractions which appear in Penguin diagrams are calculated by the random noise method, and the ΔI=1/2ΔI=1/2 matrix elements which require subtractions with the quark loop contractions are obtained with a statistical accuracy of about 10%. We confirm the chiral properties required of the K→πK\toπ matrix elements. Matching the lattice matrix elements to those in the continuum at μ=1/aμ=1/a using the perturbative renormalization factor to one loop order, and running to the scale μ=mc=1.3μ=m_c=1.3 GeV with the renormalization group for Nf=3N_f=3 flavors, we calculate all the matrix elements needed for the decay amplitudes. With these matrix elements, the ΔI=3/2ΔI=3/2 decay amplitude shows a good agreement with experiment in the chiral limit. The ΔI=1/2ΔI=1/2 amplitude, on the other hand, is about 50--60% of the experimental one even after chiral extrapolation. In view ofthe insufficient enhancement of the ΔI=1/2ΔI=1/2 contribution, we employ the experimental values for the real parts of the decay amplitudes in our calculation of ε′/εε'/ε. We find that the ΔI=3/2ΔI=3/2 contribution is larger than the ΔI=1/2ΔI=1/2 contribution so that ε′/εε'/ε is negative and has a magnitude of order 10−410^{-4}. Possible reasons for these unsatisfactory results are discussed.

preprint2003arXiv

Continuum limit of proton decay matrix elements in quenched lattice QCD

We present a lattice QCD calculation of the parameters αand βwhich are necessary in the theoretical estimation of the proton lifetime in grand unified theories (GUTs) using chiral lagrangian approach. The simulation is carried out using the Wilson quark action at three gauge coupling constants in the quenched approximation. We obtain |α(2GeV)|=0.0091(08)(^{+10}_{-19})GeV^3 and |β(2GeV)|=0.0098(08)(^{+10}_{-20})GeV^3 in the continuum limit where the first error is statistical and the second one is due to scale setting.

preprint2003arXiv

Heavy quark expansion parameters from lattice NRQCD

We present a lattice QCD calculation of the heavy quark expansion parameters μπ2μ_π^2 and μG2μ_G^2 for heavy-light mesons and heavy-light-light baryons. The calculation is carried out on a 203×^3\times48 lattice at ββ = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy quarks. We obtain the parameters μπ2μ_π^2 and μG2μ_G^2 in two different methods: a direct calculation of the matrix elements and an indirect calculation through the mass spectrum, and confirm that the both methods give consistent results. We also discuss an application to the lifetime ratios.

preprint2003arXiv

Light hadron spectrum in three-flavor QCD with O(a)-improved Wilson quark action

We report on a calculation of the light hadron spectrum and quark masses in three-flavor dynamical QCD using the non-perturbatively O(a)-improved Wilson quark action and a renormalization-group improved gauge action. Simulations are carried out on a 16^3 \times 32 lattice at β=1.9, where a^{-1} \simeq 2GeV, with 6 ud quark masses corresponding to m_{pi}/m_{rho} \simeq 0.64-0.77 and 2 s quark masses close to the physical value. We observe that the inclusion of dynamical strange quark brings the lattice QCD meson spectrum to good agreement with experiment. Dynamical strange quarks also lead to a reduction of the uds quark masses by about 15%.

preprint2003arXiv

Non-perturbative calculation of ZVZ_V and ZAZ_A in domain-wall QCD on a finite box

We report on a non-perturbative evaluation of the renormalization factors for the vector and axial-vector currents, ZVZ_V and ZAZ_A, in the quenched domain-wall QCD (DWQCD) with plaquette and renormalization group improved gauge actions. We take the Dirichlet boundary condition for both gauge and domain-wall fermion fields on the finite box, and introduce the flavor-chiral Ward-Takahashi identities to calculate the renormalization factors. As a test of the method, we numerically confirm the expected relation that ZV≃ZAZ_V \simeq Z_A in DWQCD. Employing two different box sizes for the numerical simulations at several values of the gauge coupling constant g2g^2 and the domain-wall height MM, we extrapolate ZVZ_V to the infinite volume to remove a/La/L errors. We finally give the interpolation formula of ZVZ_V in the infinite volume as a function of g2g^2 and MM.

preprint2003arXiv

Non-perturbative renormalization of vector and axial vector currents in quenched QCD for a renormalization group improved gauge action

Renormalization constants of vector (ZVZ_V) and axial-vector (ZAZ_A) currents are determined non-perturbatively in quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schrödinger functional method. Meson decay constants fρf_ρ and fπf_π show much better scaling when ZVZ_V and ZAZ_A estimated for infinite physical volume are used instead of ZZ-factors from tadpole-improved one-loop perturbation theory.

preprint2003arXiv

Study of finite volume effects in the non-perturbative determination of $\csw$ with the SF method in full three-flavor lattice QCD

The non-perturbative $\csw$ determined by the Schrödinger functional (SF) method with the RG-improved gauge action in dynamical Nf=3N_f=3 QCD shows a finite volume effect when the numerical simulations are carried out at a constant lattice size L/aL/a. We remove the unwanted finite volume effect by keeping physical lattice extent LL at a constant. The details of the method and the result obtained for non-perturbative $\csw$ with a constant LL are reported.

preprint2002arXiv

Charmonium Spectrum from Quenched Anisotropic Lattice QCD

We present a detailed study of the charmonium spectrum using anisotropic lattice QCD. We first derive a tree-level improved clover quark action on the anisotropic lattice for arbitrary quark mass. The heavy quark mass dependences of the improvement coefficients, i.e. the ratio of the hopping parameters ζ=Kt/Ksζ=K_t/K_s and the clover coefficients cs,tc_{s,t}, are examined at the tree level. We then compute the charmonium spectrum in the quenched approximation employing ξ=as/at=3ξ= a_s/a_t = 3 anisotropic lattices. Simulations are made with the standard anisotropic gauge action and the anisotropic clover quark action at four lattice spacings in the range asa_s=0.07-0.2 fm. The clover coefficients cs,tc_{s,t} are estimated from tree-level tadpole improvement. On the other hand, for the ratio of the hopping parameters ζζ, we adopt both the tree-level tadpole-improved value and a non-perturbative one. We calculate the spectrum of S- and P-states and their excitations. The results largely depend on the scale input even in the continuum limit, showing a quenching effect. When the lattice spacing is determined from the 1P−1S1P-1S splitting, the deviation from the experimental value is estimated to be ∼\sim30% for the S-state hyperfine splitting and ∼\sim20% for the P-state fine structure. Our results are consistent with previous results at ξ=2ξ= 2 obtained by Chen when the lattice spacing is determined from the Sommer scale r0r_0. We also address the problem with the hyperfine splitting that different choices of the clover coefficients lead to disagreeing results in the continuum limit.

preprint2002arXiv

Exploring QCD at small sea quark masses with improved Wilson-type quarks

We explore the region of small sea quark masses below mPS/mV=0.5m_{PS}/m_V=0.5 in two-flavor QCD using a mean-field improved clover quark action and an RG-improved gauge action at a≃0.2a \simeq 0.2 fm on 123×2412^3 \times 24 and 163×2416^3 \times 24 lattices. We find that instability of the standard BiCGStab algorithm at small quark masses can be mostly removed by the BiCGStab(DS-LL) algorithm, which employs LL-th minimal residual polynomials with a dynamical selection of LL. We also find singular spikes of ΔHΔH in the HMC algorithm at moderate values of ΔτΔτ. Nature of the spike is studied. We also study finite-size effects and chiral properties of meson masses.

preprint2002arXiv

Full one-loop electroweak radiative corrections to single Higgs production in e+ e-

We present the full O(α){\cal O}(α) electroweak radiative corrections to single Higgs production in \epemt. This takes into account the full one-loop corrections as well as the effects of hard photon radiation. We include both the fusion and Higgs-strahlung processes. The computation is performed with the help of {\tt GRACE-loop} where we have implemented a generalised non-linear gauge fixing condition. The latter includes 5 gauge parameters that can be used for checks on our results. Besides the UV, IR finiteness and gauge parameter independence checks it proves also powerful to test our implementation of the 5-point function. We find that for a 500GeV machine and a light Higgs of mass 150GeV, the total O(α){\cal O}(α) correction is small when the results are expressed in terms of αQEDα_{\rm QED}. The total correction decreases slightly for higher energies. For moderate centre of mass energies the total O(α){\cal O}(α) decreases as the Higgs mass increases, reaching -10% for MH=350M_H=350GeV and s=500\sqrt{s}=500GeV. In order to quantify the genuine weak corrections we have subtracted the universal virtual and bremsstrahlung correction from the full O(α){\cal O}(α). We find, for MH=150M_H=150GeV, a weak correction slowly decreasing from -2% to -4% as the energy increases from s=300\sqrt{s}=300GeV to s=1\sqrt{s}=1TeV after expressing the tree-level results in terms of GμG_μ

preprint2002arXiv

I=2 Pion Scattering Length with the Wilson Fermion

The calculation of the I=2 pion scattering length in quenched lattice QCD is revisited. The calculation is carried out with the Wilson fermion action employing Lüscher's finite size scaling method at β=5.9β=5.9, 6.1, and 6.3 corresponding to the range of lattice spacing a≃0.12−0.07a\simeq 0.12 - 0.07 fm. We obtain in the continuum limit a0mπ=−0.0410(69)a_0m_π= -0.0410(69), which is consistent with the prediction of chiral perturbation theory a0mπ=−0.0444(10)a_0m_π=-0.0444(10).

preprint2002arXiv

I=2 Pion Scattering Phase Shift with Wilson Fermions

We present a lattice QCD calculation of the scattering phase shift for the I=2 SS-wave two-pion system using the finite size method proposed by Lüscher. We work in the quenched approximation employing the standard plaquette action at β=5.9β=5.9 for gluons and the Wilson fermion action for quarks. The phase shift is extracted from the energy eigenvalues of the two-pion system, which are obtained by a diagonalization of the pion 4-point function evaluated for a set of relative spatial momenta. In order to change momentum of the two-pion system, calculations are carried out on 243×6024^3\times 60, 323×6032^3\times 60, and 483×6048^3\times 60 lattices. The phase shift is successfully calculated over the momentum range 0<p2<0.3GeV20 < p^2 < 0.3 {\rm GeV}^2.

preprint2002arXiv

Light Hadron Spectrum and Quark Masses from Quenched Lattice QCD

We present details of simulations for the light hadron spectrum in quenched QCD carried out on the CP-PACS parallel computer. Simulations are made with the Wilson quark action and the plaquette gauge action on 32^3x56 - 64^3x112 lattices at four lattice spacings (a \approx 0.1-0.05 fm) and the spatial extent of 3 fm. Hadronic observables are calculated at five quark masses (m_{PS}/m_V \approx 0.75 - 0.4), assuming the u and d quarks being degenerate but treating the s quark separately. We find that the presence of quenched chiral singularities is supported from an analysis of the pseudoscalar meson data. We take m_π, m_ρand m_K (or m_ϕ) as input. After chiral and continuum extrapolations, the agreement of the calculated mass spectrum with experiment is at a 10% level. In comparison with the statistical accuracy of 1-3% and systematic errors of at most 1.7% we have achieved, this demonstrates a failure of the quenched approximation for the hadron spectrum: the meson hyperfine splitting is too small, and the octet masses and the decuplet mass splittings are both smaller than experiment. Light quark masses are calculated using two definitions: the conventional one and the one based on the axial-vector Ward identity. The two results converge toward the continuum limit, yielding m_{ud}=4.29(14)^{+0.51}_{-0.79} MeV. The s quark mass depends on the strange hadron mass chosen for input: m_s = 113.8(2.3)^{+5.8}_{-2.9} MeV from m_K and m_s = 142.3(5.8)^{+22.0}_{-0} MeV from m_ϕ, indicating again a failure of the quenched approximation. We obtain Λ_{\bar{MS}}^{(0)}= 219.5(5.4) MeV. An O(10%) deviation from experiment is observed in the pseudoscalar meson decay constants.

preprint2002arXiv

Polynomial Hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors

We present a polynomial hybrid Monte Carlo (PHMC) algorithm for lattice QCD with odd numbers of flavors of O(a)-improved Wilson quark action. The algorithm makes use of the non-Hermitian Chebyshev polynomial to approximate the inverse square root of the fermion matrix required for an odd number of flavors. The systematic error from the polynomial approximation is removed by a noisy Metropolis test for which a new method is developed. Investigating the property of our PHMC algorithm in the N_f=2 QCD case, we find that it is as efficient as the conventional HMC algorithm for a moderately large lattice size (16^3 times 48) with intermediate quark masses (m_{PS}/m_V ~ 0.7-0.8). We test our odd-flavor algorithm through extensive simulations of two-flavor QCD treated as an N_f = 1+1 system, and comparing the results with those of the established algorithms for N_f=2 QCD. These tests establish that our PHMC algorithm works on a moderately large lattice size with intermediate quark masses (16^3 times 48, m_{PS}/m_V ~ 0.7-0.8). Finally we experiment with the (2+1)-flavor QCD simulation on small lattices (4^3 times 8 and 8^3 times 16), and confirm the agreement of our results with those obtained with the R algorithm and extrapolated to a zero molecular dynamics step size.

preprint2001arXiv

Charmonium spectrum from quenched QCD on anisotropic lattices

We present our final results of the charmonium spectrum in quenched QCD on anisotropic lattices. Simulations are made with the plaquette gauge action and a tadpole improved clover quark action employing ξ=as/at=3ξ= a_s/a_t = 3. We calculate the spectrum of S- and P-states and their excitation, and study the scaling behavior of mass splittings. Comparison is made with the experiment and previous lattice results. The issue of hyperfine splitting for different choices of the clover coefficients obtained by Klassen is discussed.

preprint2001arXiv

Chiral property of domain-wall fermion from eigenvalues of 4D Wilson-Dirac Operator

We investigate a chiral property of the domain-wall fermion (DWF) system using the four-dimensional hermitian Wilson-Dirac operator HWH_W. A formula expressing the Ward-Takahashi identity quark mass m5qm_{5q} with eigenvalues of this operator is derived, which well explains the N5N_5 dependence of m5qm_{5q} observed in previous numerical simulations. We further discuss the chiral property of DWF in the large volume in terms of the spectra of HWH_W.

preprint2001arXiv

Equation of state for pure SU(3) gauge theory on anisotropic lattices

We present results for the equation of state for pure SU(3) gauge theory obtained on anisotropic lattices with the anisotropy ξ≡as/at=2ξ\equiv a_s/a_t = 2. The pressure and energy density are calculated on Nt/ξ=4,5N_t / ξ= 4, 5 and 6 lattices with the integral method. They are found to satisfy the leading 1/Nt21/N_t^2 scaling from our coarsest lattice Nt/ξ=4N_t/ξ=4. This enables us to carry out well controlled continuum extrapolations. We find that the pressure and energy density agree with those obtained using the isotropic plaquette action, but have smaller and more reliable errors.

preprint2001arXiv

Equation of state in finite-temperature QCD with two flavors of improved Wilson quarks

We present results of a first study of equation of state in finite-temperature QCD with two flavors of Wilson-type quarks. Simulations are made on lattices with temporal size Nt=4N_t=4 and 6, using an RG-improved action for the gluon sector and a meanfield-improved clover action for the quark sector. The lines of constant physics corresponding to fixed values of the ratio mPS/mVm_{\rm PS}/m_{\rm V} of the pseudo-scalar to vector meson masses at zero temperature are determined, and the beta functions which describe the renormalization-group flow along these lines are calculated. Using these results, the energy density and the pressure are calculated as functions of temperature along the lines of constant physics in the range mPS/mV=0.65m_{\rm PS}/m_{\rm V} = 0.65--0.95. The quark mass dependence in the equation of state is found to be small for $m_{\rm PS}/m_{\rm V} \simlt 0.8$. Comparison of results for Nt=4N_t=4 and Nt=6N_t=6 lattices show significant scaling violation present in the Nt=4N_t=4 results. At high temperatures the results for Nt=6N_t=6 are quite close to the continuum Stefan-Boltzmann limit, suggesting the possibility of a precise continuum extrapolation of thermodynamic quantities from simulations at $N_t\simgt 6$.

preprint2001arXiv

Hybrid Quarkonia with Dynamical Sea Quarks

We present a dynamical lattice calculation with 2 flavours for bottomonium states with an additional gluonic excitation. Using improved actions for the quarks and gauge fields at a lattice spacing of a≈0.1a \approx 0.1 fm, we find 10.977(61)(62) GeV for the energy of the lowest lying bbˉgb\bar bg-hybrid, where the first error is statistical and the second denotes the systematic uncertainty due to the determination of scale. In a parallel quenched simulation we demonstrate explicitly that vacuum polarisation effects are less than 10% of the splitting with the ground state.

preprint2001arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of BKB_K with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at β=2.6β=2.6 and 2.9 corresponding to a−1≈2a^{-1}\approx 2GeV and 3GeV. Effects due to finite fifth dimensional size N5N_5 and finite spatial size NσN_σ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain BK(μ=2GeV)=0.5746(61)(191)B_K(μ= 2 GeV)= 0.5746(61)(191), where the first error is statistical and the second error represents an estimate of scaling violation and O(α2){\cal O}(α^2) errors in the renormalization factor added in quadrature, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is consistent, albeit somewhat small, with BK(μ=2GeV)=0.628(42)B_K(μ= 2 {GeV})= 0.628(42) obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for light quark masses are also reported.

preprint2001arXiv

Light Hadron Spectroscopy with Two Flavors of Dynamical Quarks on the Lattice

We present results of a numerical calculation of lattice QCD with two degenerate flavors of dynamical quarks, identified with up and down quarks, and with a strange quark treated in the quenched approximation. The lattice action and simulation parameters are chosen with a view to carrying out an extrapolation to the continuum limit as well as chiral extrapolations. Gauge configurations are generated with a renormalization-group improved gauge action and a mean field improved clover quark action at three values of ββ and four sea quark masses. The sizes of lattice are chosen so that the physical spatial size is kept constant. Hadron masses, light quark masses and meson decay constants are measured at five valence quark masses. We also carry out complementary quenched simulations with the same improved actions. The quenched spectrum from this analysis agrees well in the continuum limit with the one of our earlier work using the standard action. We find the two-flavor full QCD meson masses in the continuum limit to be much closer to experimental meson masses than those from quenched QCD. We take these results as manifestations of sea quark effects in two-flavor full QCD. For baryon masses full QCD values for strange baryons are in agreement with experiment, while they differ increasingly with decreasing strange quark content, resulting in a nucleon mass higher than experiment. The pattern suggests finite size effects as a possible origin for this deviation. For light quark masses in the continuum limit we obtain values which are reduced by about 25% compared to the values in quenched QCD. We also present results for decay constants where large scaling violations obstruct a continuum extrapolation. Need for a non-perturbative estimate of renormalization factors is discussed.

preprint2001arXiv

Maximum entropy analysis of hadron spectral functions and excited states in quenched lattice QCD

Employing the maximum entropy method we extract the spectral functions from meson correlators at four lattice spacings in quenched QCD with the Wilson quark action. We confirm that the masses and decay constants, obtained from the position and the area of peaks, agree well with the results from the conventional exponential fit. For the first excited state, we obtain mπ1=660(590)m_{π_1} = 660(590) MeV, mρ1=1540(570)m_{ρ_1} = 1540(570) MeV, and fρ1=0.085(36)f_{ρ_1} = 0.085(36) in the continuum limit.

preprint2001arXiv

Spectral function and excited states in lattice QCD with maximum entropy method

We apply the maximum entropy method to extract the spectral functions for pseudoscalar and vector mesons from hadron correlators previously calculated at four different lattice spacings in quenched QCD with the Wilson quark action. We determine masses and decay constants for the ground and excited states of the pseudoscalar and vector channels from position and area of peaks in the spectral functions. We obtain the results, mπ1=660(590)m_{π_1} = 660(590) MeV and mρ1=1540(570)m_{ρ_1} = 1540(570) MeV for the masses of the first excited state masses, in the continuum limit of quenched QCD. We also find unphysical states which have infinite mass in the continuum limit, and argue that they are bound states of two doublers of the Wilson quark action. If the interpretation is correct, this is the first time that the state of doublers is identified in lattice QCD numerical simulations.

preprint2001arXiv

Topological Susceptibility in Lattice QCD with Two Flavors of Dynamical Quarks

We present a study of the topological susceptibility in lattice QCD with two degenerate flavors of dynamical quarks. The topological charge is measured on gauge configurations generated with a renormalization group improved gauge action and a mean field improved clover quark action at three values of β=6/g2β=6/g^2, corresponding to lattice spacings of a≈0.22a \approx 0.22, 0.16 and 0.11 fm, with four sea quark masses at each ββ. The study is supplemented by simulations of pure SU(3) gauge theory with the same gauge action at 5 values of ββ with lattice spacings 0.09 fm$\simlt a \simlt$0.27 fm. We employ a field theoretic definition of the topological charge together with cooling. For the topological susceptibility in the continuum limit of pure SU(3) gauge theory we obtain χt1/4=197−16+13χ_t^{1/4} = 197^{+13}_{-16} MeV where the error shows statistical and systematic ones added in quadrature. In full QCD χtχ_t at heavy sea quark masses is consistent with that of pure SU(3) gauge theory. A decrease of χtχ_t toward light quark masses, as predicted by the anomalous Ward-Takahashi identity for U(1) chiral symmetry, becomes clearer for smaller lattice spacings. The cross-over in the behavior of χtχ_t from heavy to light sea quark masses is discussed.

preprint2000arXiv

Calculation of K→πK\to π matrix elements in quenched domain-wall QCD

We explore the possibility for an evaluation of non-leptonic ΔS=1ΔS=1 KK decay amplitudes through the calculation of K→πK\toπ matrix elements using domain-wall QCD. The relation between the physical K→ππK\toππ matrix elements and K→πK\toπ matrix elements deduced from chiral perturbation theory is recapitulated. Quenched numerical simulations are performed on an 163×32×1616^3\times 32\times 16 lattice at a lattice spacing a−1≈2a^{-1}\approx 2GeV for both the standard plaquette gauge action and a renormalization-group improved gauge action, and reasonable signals for K→πK\toπ matrix elements are obtained. Preliminary results are reported on the K→ππK\toππ matrix elements, and results from two actions are compared.

preprint2000arXiv

Decay constants of B and D mesons from improved relativistic lattice QCD with two flavours of sea quarks

We present a calculation of the B and D meson decay constants in lattice QCD with two (Nf=2) flavours of light dynamical quarks, using an O(a)-improved Wilson action for both light and heavy quarks and a renormalization-group improved gauge action. Simulations are made at three values of lattice spacing a=0.22, 0.16, 0.11 fm and four values of sea quark mass in the range m_PS/m_V \~= 0.8-0.6. Our estimate for the continuum values of the decay constants are fBd = 208(10)(11) MeV, fBs = 250(10)(13)(^{+8}_{-0}) MeV, fDd = 225(14)(14) MeV, fDs = 267(13)(17)(^{+10}_{-0}) MeV for Nf=2 where the statistical and systematic errors are separately listed, and the third error for fBs and fDs show uncertainty of determination of strange quark mass. We also carry out a set of quenched simulations using the same action to make a direct examination of sea quark effects. Taking the ratio of results for Nf=2 and Nf=0, we obtain fb^{Nf=2}/fb^{Nf=0} = 1.11(6), fbs^{Nf=2}/fbs^{Nf=0} = 1.14(5), fd^{Nf=2}/\fd^{Nf=0} = 1.03(6), fds^{Nf=2}/\fds^{Nf=0} = 1.07(5). They show a 10-15% increase in the Nf=2 results over those of Nf=0 for the B meson decay constants, while evidence for such a trend is statistically less clear for the D meson decay constants.

preprint2000arXiv

Dynamical Quark Effects on Light Quark Masses

We present results for light quark masses from a systematic lattice QCD study with two degenerate flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to mPS/mV≈0.8m_{\rm PS}/m_{\rm V} \approx 0.8--0.6 and the lattice spacing a=0.22a=0.22--0.11 fm. In the continuum limit we find mudMSˉ(2GeV)=3.44−0.22+0.14m_{ud}^{\bar{MS}}(2 {\rm GeV})=3.44^{+0.14}_{-0.22} MeV using the ππ and ρρ meson masses as physical input, and msMSˉ(2GeV)=88−6+4m_s^{\bar{MS}}(2 {\rm GeV})=88^{+4}_{-6} MeV or 90−11+590^{+5}_{-11} MeV with the KK or ϕϕ meson mass as additional input. The quoted errors represent statistical and systematic combined, the latter including those from continuum and chiral extrapolations, and from renormalization factors. Compared to quenched results, two flavors of dynamical quarks reduce quark masses by about 25%.

preprint2000arXiv

Eigenvalues of the hermitian Wilson-Dirac operator and chiral properties of the domain-wall fermion

Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous quark mass m_{5q} and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the results for m_{5q}.

preprint2000arXiv

Full QCD Light Hadron Spectrum and Quark Masses: Final Results from CP-PACS

We present the final results of the CP-PACS calculation of the light hadron spectrum and quark masses with two flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to mPS/mV≈0.8m_{\rm PS}/m_{\rm V} \approx 0.8--0.6 and the lattice spacing a=0.22a=0.22--0.11 fm. For the meson spectrum in the continuum limit a clearly improved agreement with experiment is observed compared to the quenched case, demonstrating the importance of sea quark effects. For light quark masses we obtain mudMSˉ(2GeV)=3.44−0.22+0.14m_{ud}^{\bar{MS}}(2GeV)=3.44^{+0.14}_{-0.22} MeV and msMSˉ(2GeV)=88−6+4m_s^{\bar{MS}}(2GeV)=88^{+4}_{-6} MeV (KK-input) and msMSˉ(2GeV)=90−11+5m_s^{\bar{MS}}(2GeV)=90^{+5}_{-11} MeV (ϕϕ-input), which are reduced by about 25% compared to the values in quenched QCD.

preprint2000arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of BKB_K with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at β=2.6β=2.6 and 2.9 corresponding to a−1≈2a^{-1}\approx 2GeV and 3GeV. Effects due to finite fifth dimensional size N5N_5 and finite spatial size NσN_σ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain BK(μ=2GeV)=0.5746(61)B_K(μ= 2 {GeV})= 0.5746(61), where the error is statistical only, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is smaller but consistent with BK(μ=2GeV)=0.628(42)B_K(μ= 2 {GeV})= 0.628(42) obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for strange quark mass are also reported.

preprint2000arXiv

Phase structure and critical temperature of two-flavor QCD with a renormalization group improved gauge action and clover improved Wilson quark action

We study the finite-temperature phase structure and the transition temperature of QCD with two flavors of dynamical quarks on a lattice with the temporal size Nt=4N_t=4, using a renormalization group improved gauge action and the Wilson quark action improved by the clover term. The region of a parity-broken phase is identified, and the finite-temperature transition line is located on a two-dimensional parameter space of the coupling (β=6/g2β=6/g^2) and hopping parameter KK. Near the chiral transition point, defined as the crossing point of the critical line of the vanishing pion mass and the line of finite-temperature transition, the system exhibits behavior well described by the scaling exponents of the three-dimensional O(4) spin model. This indicates a second-order chiral transition in the continuum limit. The transition temperature in the chiral limit is estimated to be Tc=171(4)T_c = 171(4) MeV.

preprint2000arXiv

Quenched charmonium spectrum on anisotropic lattices

We present the results of quenched charmonium spectrum for S- and P-states, obtained by a relativistic heavy quark method on anisotropic lattices. Simulations are carried out using the standard plaquette gauge action and a meanfield-improved clover quark action at at−1=3a_t^{-1} = 3--6 GeV with the renormalized anisotropy fixed to ξ≡as/at=3ξ\equiv a_s/a_t =3. We study the scaling of our fine and hyperfine mass splittings, and compare with previous results.

preprint1999arXiv

Comparative Study of full QCD Hadron Spectrum and Static Quark Potential with Improved Actions

We investigate effects of action improvement on the light hadron spectrum and the static quark potential in two-flavor QCD for a−1≈1a^{-1} \approx 1 GeV and mPS/mV=0.7−0.9m_{PS}/m_V = 0.7-0.9. We compare a renormalization group improved action with the plaquette action for gluons, and the SW-clover action with the Wilson action for quarks. We find a significant improvement in the hadron spectrum by improving the quark action, while the gluon improvement is crucial for a rotationally invariant static potential. We also explore the region of light quark masses corresponding to mPS/mV≥0.4m_{PS}/m_V \geq 0.4 on a 2.7 fm lattice using the improved gauge and quark action. A flattening of the potential is not observed up to 2 fm.

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with renormalization group improved action

A lattice study of the equation of state for pure SU(3) gauge theory using a renormalization-group (RG) improved action is presented. The energy density and pressure are calculated on a 163×416^3\times 4 and a 323×832^3\times 8 lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with RG improved action

We present results for the equation of state for pure SU(3) gauge theory obtained with a renormalization-group (RG) improved action. The energy density and pressure are calculated on a 163×416^3\times 4 and a 323×832^3\times 8 lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.

preprint1999arXiv

Heavy-light decay constants from clover heavy quark action in QCD with two flavors of dynamical quarks

We present results on an analysis of the decay constants f_B and f_Bs with two flavours of sea quark. The calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge coupling, with m_pi/m_rho ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with heavy quarks whose mass is in the range of the b-quark. A detailed comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.

preprint1999arXiv

Heavy-light decay constants from relativistic Nf=2,0 QCD

We present results on an analysis of the decay constants fB and fBs with two flavours of sea quark. The calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge coupling, with m_pi/m_rho ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with heavy quarks whose mass is in the range of the b-quark. A comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.

preprint1999arXiv

Heavy-light spectrum and decay constant from NRQCD with two flavors of dynamical quarks

We report on a study of B mesons on N_f = 2 full QCD configurations using an RG-improved gauge action, NRQCD heavy quark action and tadpole-improved clover light quark action. Results on the heavy-light spectrum and the decay constants from 16^3x32 lattices at a^{-1} ~ 1.5 GeV are presented, and compared with quenched results obtained with the same action combination at matching lattice spacings.

preprint1999arXiv

Hybrid Quarkonia on Asymmetric Lattices

We report on a study of heavy hybrid states using the NRQCD approach on coarse and asymmetric lattices, where we discard vacuum polarisation effects and neglect all spin-correction terms. We find a clear hybrid signal on all our lattices (as=0.15...0.47a_s= 0.15 ... 0.47 fm). We have studied in detail the lattice spacing artefacts, finite volume effects and mass dependence. Within the above approximations we predict the hybrid excitation in Charmonium to be 1.323(13) GeV above its ground state. The bottomonium hybrid was found to be 1.542(8) GeV above its ground state.

preprint1999arXiv

Light hadron spectrum and quark masses in QCD with two flavors of dynamical quarks

We present updated results of the CP-PACS calculation of the light hadron spectrum in Nf=2N_{\rm f}=2 full QCD. Simulations are made with an RG-improved gauge action and a tadpole-improved clover quark action for sea quark masses corresponding to mPS/mV≈0.8m_{\rm PS}/m_{\rm V} \approx 0.8--0.6 and the lattice spacing a=0.22a=0.22--0.09 fm. A comparison of the full QCD spectrum with new quenched results, obtained with the same improved action, shows clearly the existence of sea quark effects in vector meson masses. Results for light quark masses in Nf=2N_{\rm f}=2 QCD are also presented.

preprint1999arXiv

Quenched Light Hadron Spectrum

We present results of a large-scale simulation for the flavor non-singlet light hadron spectrum in quenched lattice QCD with the Wilson quark action. Hadron masses are calculated at four values of lattice spacing in the range a≈a \approx 0.1 - 0.05 fm on lattices with a physical extent of 3 fm at five quark masses corresponding to mπ/mρ≈0.75m_π/m_ρ\approx 0.75 - 0.4. The calculated spectrum in the continuum limit shows a systematic deviation from experiment, though the magnitude of deviation is contained within 11%. Results for decay constants and light quark masses are also reported.

preprint1998arXiv

Full QCD light hadron spectrum from the CP-PACS

We report on an on-going two-flavor full QCD study on CP-PACS using an RG-improved gauge action and a tadpole-improved SW quark action. Runs are made for three lattice spacings a−1≈0.9a^{-1}\approx 0.9, 1.3, and 2.5 GeV on 123×2412^3\times24, 163×3216^3\times32, and 243×4824^3\times48 lattices. Four sea quark masses having mPS/mV≈0.8m_{\rm PS}/m_{\rm V} \approx 0.8--0.6 are simulated, for each of which hadron masses are evaluated for valence quark masses corresponding to mPS/mV≈0.8m_{\rm PS}/m_{\rm V} \approx 0.8--0.5. Results for hadron and light quark masses are presented and compared with those obtained in quenched QCD.

preprint1998arXiv

Quenched Light Hadron Spectrum with the Wilson Quark Action: Final Results from CP-PACS

We report the final results of the CP-PACS calculation for the quenched light hadron spectrum with the Wilson quark action. Our data support the presence of quenched chiral singularities, and this motivates us to use mass formulae based on quenched chiral perturbation theory in order to extrapolate hadron masses to the physical point. Hadron masses and decay constants in the continuum limit show unambiguous systematic deviations from experiment. We also report the results for light quark masses.

preprint1998arXiv

The static quark potential in full QCD

We report results on the static quark potential in two-flavor full QCD. The calculation is performed for three values of lattice spacing a−1≈0.9,1.3a^{-1}\approx 0.9, 1.3 and 2.5 GeV on 123×24,163×3212^3{\times}24, 16^3{\times}32 and 243×4824^3{\times}48 lattices respectively, at sea quark masses corresponding to mπ/mρ≈0.8−0.6m_π/m_ρ\approx 0.8-0.6. An RG-improved gauge action and a tadpole-improved SW clover quark action are employed. We discuss scaling of mρ/σm_ρ/\sqrtσ and effects of dynamical quarks on the potential.

preprint1997arXiv

Hadron spectroscopy and static quark potential in full QCD: A comparison of improved actions on the CP-PACS

We present first results from a full QCD calculation on the CP-PACS, comparing various actions at a−1∼1GeVa^{-1} \sim 1 GeV and mπ/mρ≈0.7m_π/m_ρ\approx 0.7--0.9. We use the plaquette and a renormalization group improved action for the gluons, and the Wilson and the SW-Clover action for quarks. We find that significant improvements in the hadron spectrum results from improving the quarks, while the gluon improvement is required for a rotationally invariant static potential. An ongoing effort towards exploring the chiral limit in full QCD is described.

preprint1997arXiv

Scaling in SU(3) Pure Gauge Theory with a Renormalization Group Improved Action

We study the scaling properties of the static quark potential and the ratio of the critical temperature TcT_c to the square root of the string tension σσ in the SU(3) pure gauge theory using a renormalization group improved action. We first determine the critical coupling βcβ_c on lattices with temporal extension Nt=3N_t=3, 4, and 6, and then calculate the static quark potential at the critical couplings on lattices at zero temperature. We note that the static quark potentials obtained are rotationally invariant with errors of at most 1 - 2 % in all the three cases, and that the potential V(R)V(R) in physical units scales in the whole region of RR investigated. The values of Tc/σT_{c}/\sqrtσ for the three cases in the infinite volume limit are identical within errors. We estimate the value in the continuum limit to be Tc/σ=0.656(4)T_{c}/\sqrtσ = 0.656(4), which is slightly larger than the value in the continuum limit from the one-plaquette action, 0.629(3).

preprint1997arXiv

susy23 v2.0: an Event Generator for Supersymmetric Processes at e+e−e^+e^- Colliders

\susy23 is a Monte-Carlo package for generating supersymmetric (SUSY) processes at e+e−e^+e^- colliders. Twenty-three types of SUSY processes with 2 or 3 final state particles at tree level are included in version 2.0. {\susy23} addresses event simulation requirements at e+e−e^+e^- colliders such as LEP. Matrix elements are generated by GRACE with the helicity amplitude method for processes involving massive fermions. The phase space integration of the matrix element gives the total and differential cross sections, then unweighted events are generated. Sparticle widths and decay branching ratios are calculated. Each final state particle may then decay according to these probabilities. Spin correlations are taken into account in the decays of sparticles. Corrections of initial state radiation (ISR) are implemented in two ways, one is based on the electron structure function formalism and the second uses the parton shower algorithm called QEDPS. Parton shower and hadronization of the final quarks are performed through an interface to JETSET.

preprint1996arXiv

grc4f v1.0: a Four-fermion Event Generator for e+e- Collisions

grc4f is a Monte-Carlo package for generating e+e- to 4-fermion processes in the standard model. All of the 76 LEP-2 allowed fermionic final state processes evaluated at tree level are included in version 1.0. grc4f addresses event simulation requirements at e+e- colliders such as LEP and up-coming linear colliders. Most of the attractive aspects of grc4f come from its link to the GRACE system: a Feynman diagram automatic computation system. The GRACE system has been used to produce the computational code for all final states, giving a higher level of confidence in the calculation correctness. Based on the helicity amplitude calculation technique, all fermion masses can be kept finite and helicity information can be propagated down to the final state particles. The phase space integration of the matrix element gives the total and differential cross sections, then unweighted events are Generated. Initial state radiation (ISR) corrections are implemented in two ways, one is based on the electron structure function formalism and the second uses the parton shower algorithm called QEDPS. The latter can also be applied for final state radiation (FSR) though the interference with the ISR is not yet taken into account. Parton shower and hadronization of the final quarks are performed through an interface to JETSET. Coulomb correction between two intermediate W's, anomalous coupling as well as gluon contributions in the hadronic processes are also included.

preprint1996arXiv

Scaling of the critical temperature and quark potential with a renormalization group improved SU(3) gauge action

We study the scaling property of the ratio of the critical temperature TcT_c to the square root of the string tension σσ in the SU(3) pure gauge theory using a renormalization group improved action. We first determine the critical coupling βcβ_c on lattices with temporal extension Nt=4N_t=4 and 6, and then calculate the static quark potential at the critical couplings on lattices at zero temperature. The values of Tc/σT_{c}/\sqrtσ in the infinite volume limit are identical within errors, while they are slightly larger than the value extrapolated to the continuum limit with the standard action. We also note that the rotational invariance of the static quark potential is remarkably restored in the both cases, and that the potential V(R)V(R) in physical units scales in the whole region of RR investigated.

preprint1995arXiv

Re-entrant Layer-by-Layer Etching of GaAs(001)

We report the first observation of re-entrant layer-by-layer etching based on {\it in situ\/} reflection high-energy electron-diffraction measurements. With AsBr3_3 used to etch GaAs(001), sustained specular-beam intensity oscillations are seen at high substrate temperatures, a decaying intensity with no oscillations at intermediate temperatures, but oscillations reappearing at still lower temperatures. Simulations of an atomistic model for the etching kinetics reproduce the temperature ranges of these three regimes and support an interpretation of the origin of this phenomenon as the site-selectivity of the etching process combined with activation barriers to interlayer adatom migration.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.