Search

We determine the neutral kaon mixing matrix element $B_K$ in the continuum limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action at two different lattice spacings. These lattice fermions have near exact chiral symmetry and therefore avoid artificial lattice operator mixing. We introduce a significant improvement to the conventional NPR method in which the bare matrix elements are renormalized non-perturbatively in the RI-MOM scheme and are then converted into the MSbar scheme using continuum perturbation theory. In addition to RI-MOM, we introduce and implement four non-exceptional intermediate momentum schemes that suppress infrared non-perturbative uncertainties in the renormalization procedure. We compute the conversion factors relating the matrix elements in this family of RI-SMOM schemes and MSbar at one-loop order. Comparison of the results obtained using these different intermediate schemes allows for a more reliable estimate of the unknown higher-order contributions and hence for a correspondingly more robust estimate of the systematic error. We also apply a recently proposed approach in which twisted boundary conditions are used to control the Symanzik expansion for off-shell vertex functions leading to a better control of the renormalization in the continuum limit. We control chiral extrapolation errors by considering both the NLO SU(2) chiral effective theory, and an analytic mass expansion. We obtain $B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding all sources of error in quadrature we obtain $\hat{B}_K = 0.749(27)_{combined}$, with an overall combined error of 3.6%., Comment: 65 pages

We present physical results obtained from simulations using 2+1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing $a$, ($a^{-1}$=\,1.73\,(3)\,GeV and $a^{-1}$=\,2.28\,(3)\,GeV). On the coarser lattice, with $24^3\times 64\times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the finer $32^3\times 64\times 16$ lattice are new. We explain how we use lattice data obtained at several values of the lattice spacing and for a range of quark masses in combined continuum-chiral fits in order to obtain results in the continuum limit and at physical quark masses. We implement this procedure at two lattice spacings, with unitary pion masses in the approximate range 290--420\,MeV (225--420\,MeV for partially quenched pions). We use the masses of the $\pi$ and $K$ mesons and the $\Omega$ baryon to determine the physical quark masses and the values of the lattice spacing. While our data are consistent with the predictions of NLO SU(2) chiral perturbation theory, they are also consistent with a simple analytic ansatz leading to an inherent uncertainty in how best to perform the chiral extrapolation that we are reluctant to reduce with model-dependent assumptions about higher order corrections. Our main results include $f_\pi=124(2)_{\rm stat}(5)_{\rm syst}$\,MeV, $f_K/f_\pi=1.204(7)(25)$ where $f_K$ is the kaon decay constant, $m_s^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(96.2\pm 2.7)$\,MeV and $m_{ud}^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(3.59\pm 0.21)$\,MeV\, ($m_s/m_{ud}=26.8\pm 1.4$) where $m_s$ and $m_{ud}$ are the mass of the strange-quark and the average of the up and down quark masses respectively, $[\Sigma^{\msbar}(2 {\rm GeV})]^{1/3} = 256(6)\; {\rm MeV}$, where $\Sigma$ is the chiral condensate, the Sommer scale $r_0=0.487(9)$\,fm and $r_1=0.333(9)$\,fm., Comment: 129 pages, 59 figures, Published version containing an extended discussion of reweighting and including a new appendix (Appendix C)

We report on the results of a numerical simulation concerning the low-lying spectrum of four-dimensional N=1 SU(2) Supersymmetric Yang-Mills (SYM) theory on the lattice with light dynamical gluinos. In the gauge sector the tree-level Symanzik improved gauge action is used, while we use the Wilson formulation in the fermion sector with stout smearing of the gauge links in the Wilson-Dirac operator. The ensembles of gauge configurations were produced with the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) updating algorithm. We performed simulations on large lattices up to a size of 24^3 x 48 at $\beta=1.6$. Using QCD units with the Sommer scale being set to r_0 = 0.5 fm, the lattice spacing is about a ~ 0.09 fm, and the spatial extent of the lattice corresponds to 2.1 fm. At the lightest simulated gluino mass the spin-1/2 gluino-glue bound state appeared to be considerably heavier than its expected super-partner, the pseudoscalar bound state. Whether supermultiplets are formed remains to be studied in upcoming simulations., Comment: 17 pages, 8 figures (eps), LaTeX2e

Results of a numerical simulation concerning the low-lying spectrum of four-dimensional N=1 SU(2) Supersymmetric Yang-Mills (SYM) theory on the lattice with light dynamical gluinos are reported. We use the tree-level Symanzik improved gauge action and Wilson fermions with stout smearing of the gauge links in the Wilson-Dirac operator. The configurations are produced with the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) algorithm. We performed simulations on lattices up to a size of 24^3x48 at \beta=1.6. Using QCD units with the Sommer scale being set to r_0=0.5 fm, the lattice spacing is about a~0.09 fm, and the spatial extent of the lattice corresponds to 2.1 fm to control finite size effects. At the lightest simulated gluino mass our results indicate a mass for the lightest gluino-glue bound state, which is considerably heavier than the values obtained for its possible superpartners. Whether supermultiplets are formed remains to be studied in upcoming simulations., Comment: 7 pages, talk presented at The XXVII International Symposium on Lattice Field Theory, July 26-31 2009, Peking University, Beijing, China

One of the major challenges of particle physics has been to gain an in-depth understanding of the role of quark flavor and measurements and theoretical interpretations of their results have advanced tremendously: apart from masses and quantum numbers of flavor particles, there now exist detailed measurements of the characteristics of their interactions allowing stringent tests of Standard Model predictions. Among the most interesting phenomena of flavor physics is the violation of the CP symmetry that has been subtle and difficult to explore. Till early 1990s observations of CP violation were confined to neutral $K$ mesons, but since then a large number of CP-violating processes have been studied in detail in neutral $B$ mesons. In parallel, measurements of the couplings of the heavy quarks and the dynamics for their decays in large samples of $K, D$, and $B$ mesons have been greatly improved in accuracy and the results are being used as probes in the search for deviations from the Standard Model. In the near future, there will be a transition from the current to a new generation of experiments, thus a review of the status of quark flavor physics is timely. This report summarizes the results of the current generation of experiments that is about to be completed and it confronts these results with the theoretical understanding of the field., Comment: Report of the CKM workshop, Rome 9-13th Sep. 2008, 340 pages, 106 postscript figures, submitted to Phys. Repts

The lattice provides a powerful tool to non-perturbatively investigate strongly coupled supersymmetric Yang-Mills (SYM) theories. The pure SU(2) SYM theory with one supercharge is simulated on large lattices with small Majorana gluino masses down to about $am_{\tilde g}=0.068$ with lattice spacing $a\simeq 0.125$ fm. The gluino dynamics is simulated by the Two-Step Multi-Boson (TSMB) and the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) algorithms. Supersymmetry (SUSY) is broken explicitly by the lattice and the Wilson term and softly by the presence of a non-vanishing gluino mass. However, the recovery of SUSY is expected in the infinite volume continuum limit by tuning the bare parameters to the SUSY point in the parameter space. This scenario is studied by the determination of the low-energy mass spectrum and by means of lattice SUSY Ward-Identities (WIs)., Comment: 5 pages, 8 figures, 1 table, presented at "Quark Confinement and the Hadron Spectrum VIII", September 1 - 6 2008, Mainz, Germany

The latest results of an ongoing project for the lattice simulation of QCD with a single quark flavor are presented. The Symanzik tree-level-improved Wilson action is adopted in the gauge sector and the (unimproved) Wilson action for the fermion. Results from new simulations with one step of Stout-smearing (rho=0.15) in the fermion action are discussed. The one-flavor theory is simulated by a polynomial hybrid Monte Carlo algorithm (PHMC) at beta=4.0 corresponding to a = 0.13fm, on 16^3x32 and 24^3x48 lattices; the box-size is L = 2.1fm and L = 3.1fm, respectively. At the lightest simulated quark mass the (partially quenched) pion mass is ~300 MeV. The masses of the lightest bound states are computed, including the flavor singlet scalar and pseudoscalar mesons sigma_s and eta_s, the scalar glueball 0^++, and the Delta^++ baryon. Relics of SUSY in the mass spectrum, expected from a large N_c orientifold equivalence with the N=1 supersymmetric Yang-Mills theory, are discussed., Comment: 7 pages, 2 figures, talk presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

We perform Monte Carlo investigations of the 4d ${\cal N}=1$ supersymmetric Yang-Mills (SYM) theory on the lattice with dynamical gluinos in the adjoint representation of the SU(2) gauge group. Our aim is to determine the mass spectrum of the low-lying bound states which is expected to be organised in supermultiplets in the infinite volume continuum limit. For this purpose we perform simulations on large lattices, up to an extension $L/r_0 \simeq 6$ where $r_0 \simeq 0.5 \rm fm$ is the Sommer scale parameter. We apply improved lattice actions: tree-level improved Symanzik (tlSym) gauge action and in the later runs a Stout-smeared Wilson fermion action. The gauge configuration samples are prepared by the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) update algorithm., Comment: 7 pages, 1 figure, 1 table, presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

We have simulated QCD using 2+1 flavors of domain wall quarks on a $(2.74 {\rm fm})^3$ volume with an inverse lattice scale of $a^{-1} = 1.729(28)$ GeV. The up and down (light) quarks are degenerate in our calculations and we have used four values for the ratio of light quark masses to the strange (heavy) quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have measured pseudoscalar meson masses and decay constants, the kaon bag parameter $B_K$ and vector meson couplings. We have used SU(2) chiral perturbation theory, which assumes only the up and down quark masses are small, and SU(3) chiral perturbation theory to extrapolate to the physical values for the light quark masses. While next-to-leading order formulae from both approaches fit our data for light quarks, we find the higher order corrections for SU(3) very large, making such fits unreliable. We also find that SU(3) does not fit our data when the quark masses are near the physical strange quark mass. Thus, we rely on SU(2) chiral perturbation theory for accurate results. We use the masses of the $\Omega$ baryon, and the $\pi$ and $K$ mesons to set the lattice scale and determine the quark masses. We then find $f_\pi = 124.1(3.6)_{\rm stat}(6.9)_{\rm syst} {\rm MeV}$, $f_K = 149.6(3.6)_{\rm stat}(6.3)_{\rm syst} {\rm MeV}$ and $f_K/f_\pi = 1.205(0.018)_{\rm stat}(0.062)_{\rm syst}$. Using non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturbation theory to relate these to $\bar{\rm MS}$ we find $ m_{ud}^{\bar{\rm MS}}(2 {\rm GeV}) = 3.72(0.16)_{\rm stat}(0.33)_{\rm ren}(0.18)_{\rm syst} {\rm MeV}$ and $m_{s}^{\bar{\rm MS}}(2 {\rm GeV}) = 107.3(4.4)_{\rm stat}(9.7)_{\rm ren}(4.9)_{\rm syst} {\rm MeV}$., Comment: 133 pages, 25 figures

The hadron spectrum of one flavor QCD is studied by Monte Carlo simulations. The Symanzik tree-level-improved Wilson action is used for the gauge field and the Wilson action for the fermion. The theory is simulated by a polynomial hybrid Monte Carlo algorithm (PHMC). The mass spectrum of hadronic bound states is investigated at two different lattice spacings: a ~ 0.37r_0 and a ~ 0.27r_0, corresponding to ~0.19fm and ~0.13fm in ordinary QCD. The lattice extension is fixed to L ~ 4.4r_0 (~2.2fm). The lightest simulated quark mass corresponds to a pion with mass ~270MeV. Properties of the theory are analyzed by making use of the ideas of partially quenched chiral perturbation theory (PQChPT). The symmetry of the single flavor theory can be artificially enhanced by adding extra valence quarks, which can be interpreted as u and d quarks. Operators in the valence pion sector can be built. Masses and decay constants are analyzed by using PQChPT formulae at next-to-leading order., Comment: 14 pages, 4 figures; Proceedings of the XXV International Symposium on Lattice Field Theory "Lattice 2007", Regensburg, Germany

We present results for meson masses and decay constants measured on 24^3x64 lattices using the domain wall fermion formulation with an extension of the fifth dimension of L_s=16 for N_f=2+1 dynamical quark flavors. The lightest dynamical meson mass in our set-up is around 331 MeV, while partially quenched mesons reach masses as low as 250 MeV. The applicability of SU(3)xSU(3) and SU(2)xSU(2) (partially quenched) chiral perturbation theory will be compared and we quote values for the low-energy constants from both approaches. We will extract the average light quark and strange quark masses and use a non-perturbative renormalization technique (RI/MOM) to quote their physical values. The pion and kaon decay constants are determined at those values from our chiral fits and their ratio is used to obtain the CKM-matrix element |V_us|. The results presented here include statistical errors only., Comment: 14 pages, 6 figures, 3 tables, joint proceedings for talks presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany

The effect of Stout smearing is investigated in numerical simulations with twisted mass Wilson quarks. The phase transition near zero quark mass is studied on 12^3x24, 16^3x32 and 24^3x48 lattices at lattice spacings a = 0.1 - 0.125 fm., Comment: 7 pages, 5 figures, Poster presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany

One-flavour QCD - a gauge theory with SU(3) colour gauge group and a fermion in the fundamental representation - is studied by Monte Carlo simulations. The mass spectrum of hadronic bound states is investigated in a volume with extensions of L ~ 4.4r_0 (~2.2fm) at two different lattice spacings: a ~ 0.37r_0 (~0.19fm) and a ~ 0.27r_0 (~0.13fm). The lattice action is Symanzik tree-level-improved Wilson action for the gauge field and (unimproved) Wilson action for the fermion., Comment: 21 pages, 4 figures; further references added

We present the first results for neutral kaon mixing using 2+1 flavors of domain wall fermions. A new approach is used to extrapolate to the physical up and down quark masses from our numerical studies with pion masses in the range 240 -- 420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and the kaon is not assumed to be light. Our main result is $B_K^{\bar{\rm MS}}(2 \mathrm{GeV}) = 0.524(10)(28)$ where the first error is statistical and the second incorporates estimates for all systematic errors., Comment: 4 pages, 4 figures. RBC and UKQCD collaborations

We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_\pi = 127(4) MeV, f_K = 157(5) MeV and f_K/f_\pi = 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values., Comment: RBC and UKQCD Collaborations. 17 pages, 14 figures. Typeset with ReVTEX4. v2: replaced with the version published in PRD with improved introduction

A first study of numerical Monte Carlo simulations with two quark doublets, a mass-degenerate one and a mass-split one, interpreted as u, d, s and c quarks, is carried out in the framework of the twisted mass Wilson lattice formulation. Tuning the bare parameters of this theory is explored on 12^3x24 and 16^3x32 lattices at lattice spacings a=0.20fm and a=0.15fm, respectively., Comment: 27 pages, 6 figures

Discretisation errors in two-flavour lattice QCD with Wilson-quarks and DBW2 gauge action are investigated by comparing numerical simulation data at two values of the bare gauge coupling. Both non-zero and zero twisted mass values are considered. The results, including also data from simulations using the Wilson plaquette gauge action, are compared to next-to-leading order chiral perturbation theory formulas., Comment: 53 pages, 11 tables, 18 figures; discussion extended

By employing a twisted mass term, we compare recent results from lattice calculations of Nf=2 dynamical Wilson fermions with Wilson Chiral Perturbation Theory (WChPT). The final goal is to determine some combinations of Gasser-Leutwyler Low Energy Constants (LECs). A wide set of data with different lattice spacings (a ~ 0.2-0.12 fm), different gauge actions (Wilson plaquette, DBW2) and different quark masses (down to the lowest pion mass allowed by lattice artifacts and including negative quark masses) provide a strong check of the applicability of WChPT in this regime and the scaling behaviours in the continuum limit., Comment: 8 pages. To appear in proceedings of Computational Hadron Physics, Nicosia, 14-17 Sep 2005. Notation made uniform with other publications

We summarize four contributions about dynamical twisted mass fermions. The resulting report covers results for N_f=2 obtained from three different gauge actions, namely the standard Wilson plaquette gauge action, the DBW2 and the tree-level Symanzik improved gauge action. In addition, first results for N_f=2+1+1 flavours of twisted mass fermions are discussed., Comment: Contribution to the proceedings of Lattice 2005 resulting from three talks and one poster (Hadron spectrum and quark masses), 21 pages, 12 figures

Twisted mass fermions allow light quarks to be explored but with the consequence that there are mass splittings, such as between the neutral and charged pion. Using a direct calculation of the connected neutral pion correlator and stochastic methods to evaluate the disconnected correlations, we determine the neutral pion mass. We explore the dependence on lattice spacing and quark mass in quenched QCD. For dynamical QCD, we determine the sign of the splitting which is linked, via chiral PT, to the nature of the phase transition at small quark mass., Comment: 6 pages, poster (hadron spectrum and quark masses) at Lattice 2005,Dublin, July 25-30

Lattice QCD with Wilson fermions generically shows the phenomenon of a first order phase transition. We study the phase structure of lattice QCD using Wilson twisted mass fermions and the Wilson plaquette gauge action are used in a range of beta values where such a first order phase transition is observed. In particular, we investigate the dependence of the first order phase transition on the value of the lattice spacing. Using only data in one phase and neglecting possible problems arising from the phase transition we are able to perform a first scaling test for physical quantities using this action., Comment: 15 pages, 7 figures, typo corrected, web-list of authors corrected

Nested multi-step stochastic correction offers a possibility to improve updating algorithms for numerical simulations of lattice gauge theories with fermions. The corresponding generalisations of the two-step multi-boson (TSMB) algorithm as well as some applications with hybrid Monte Carlo (HMC) algorithms are considered., Comment: 10 pages; discussion extended

The effect of changing the lattice action for the gluon field on the recently observed [1] first order phase transition near zero quark mass is investigated by replacing the Wilson plaquette action by the DBW2 action. The lattice action for quarks is unchanged: it is in both cases the original Wilson action. It turns out that Wilson fermions with the DBW2 gauge action have a phase structure where the minimal pion mass and the jump of the average plaquette are decreased, when compared to Wilson fermions with Wilson plaquette action at similar values of the lattice spacing. Taking the DBW2 gauge action is advantageous also from the point of view of the computational costs of numerical simulations., Comment: 34 pages, 3 tables, 14 figures; revised version to appear in EPJC

The phase structure of zero temperature twisted mass lattice QCD is investigated. We find strong metastabilities in the plaquette observable when the untwisted quark mass sweeps across zero., Comment: Talks presented at Lattice2004(spectrum), 6 pages, 6 figures

The quark mass dependence of meson masses and decay constants in N_f=2 QCD is studied at a fixed gauge coupling with different dynamical quark masses. Partially Quenched Wilson Chiral Perturbation Theory (PQWchiPT) is applied to obtain the extrapolation to zero quark mass. It is shown, that in our analysis NNLO-terms play a more important role than O(a)-terms, which can be neglected. Also we compare our results with recent CP-PACS results., Comment: Lattice2004(spectrum), 3 pages, 4 figures

The phase structure of zero temperature twisted mass lattice QCD is investigated. We find strong metastabilities in the plaquette observable when the untwisted quark mass assumes positive or negative values. We provide interpretations of this phenomenon in terms of chiral symmetry breaking and the effective potential model of Sharpe and Singleton., Comment: 31 pages, 7 figures; published version with added clarifications

Our previous calculations of the sea- and valence-quark mass dependence of the pseudoscalar meson masses and decay constants is repeated on a 16^3x32 lattice which allows for a better determination of the quantities in question. The conclusions are similar as before on the 16^4 lattice [1]. The two light dynamical quark flavours we simulate have masses in the range 1/4 m_s < m_{u,d} < 2/3 m_s. The sea quark mass dependence of f_pi and m_pi^2/m_q is well described by the next-to-leading order (NLO) Chiral Perturbation Theory (ChPT) formulas and clearly shows the presence of chiral logarithms. The valence quark mass dependence requires the presence of NNLO contributions in Partially Quenched ChPT (PQChPT) -- in addition to the NLO terms. The O(a) lattice artifacts in these quantities turn out to be small., Comment: 17 pages, 5 figures

The dependence of pseudo-scalar masses and decay constants on the sea and valence quark masses is investigated in the pseudo-Goldstone boson sector of QCD with two light quark flavours. The sea quark masses are at present in the range ${1/3}m_s \leq m_{ud} \leq {2/3}m_s$ whereas the valence quark masses satisfy ${1/2}m_{sea} \leq m_{val} \leq 2m_{sea}$. The values of the Gasser-Leutwyler low energy constants $L_4$, $L_5$, $L_6$ and $L_8$ are estimated. The computation is done with the Wilson-quark lattice action at gauge coupling $\beta=5.1$ on $16^4$ lattices. ${\cal O}(a)$ effects are taken into account by applying chiral perturbation theory for the Wilson lattice action as proposed by Rupak and Shoresh., Comment: 3 pages. Lattice2003(Spectrum)

The dependence of the pseudoscalar meson masses and decay constants on sea and valence quark masses is compared to next-to-leading order (NLO) Chiral Perturbation Theory (ChPT). The numerical simulations with two light dynamical quark flavors are performed with the Wilson-quark lattice action at gauge coupling beta=5.1 and hopping parameters kappa=0.176, 0.1765, 0.177 on a 16^4 lattice. O(a) lattice artifacts are taken into account by applying chiral perturbation theory for the Wilson lattice action. The values of the relevant combinations of Gasser-Leutwyler constants L_4, L_5, L_6 and L_8 are estimated., Comment: 21 pages, 4 figures; discussion of the sea quark mass dependence changed, references added

The dependence of the pseudoscalar meson mass and decay constant is compared to one-loop Partially Quenched Chiral Perturbation Theory (PQChPT) in a numerical simulation with two light dynamical quarks. The characteristic behaviour with chiral logarithms is observed. The values of the fitted PQChPT-parameters are in a range close to the expectation in continuum in spite of the fact that the lattice spacing is still large, namely a=0.28 fm., Comment: 11 pages, 3 figures; discussion of the results in section 3 extended

The design and development of cryogenic propellant tanks with reduced weight and production costs is described with reference to applications for the National Launch System. The development program focused on the use of an aluminum-lithium alloy to demonstrate the production capability, manufacturability, and strength inherent in the novel material. Other key parameters for the alloy include fracture toughness, stress-corrosion resistance, and conformance to NASA specifications for cryogenic propellant tanks. The commercially produced aluminum-lithium alloy product forms are shown to operate acceptably in the temperature range for cryogenic propellant tanks. The alloy under consideration and the tank design are important advances in the development of ultralightweight launch-vehicle structures.