Domain Wall Fermion Simulations with the Exact One-Flavor Algorithm
|
|
- Jewel George
- 5 years ago
- Views:
Transcription
1 Domain Wall Fermion Simulations with the Exact One-Flavor Algorithm David Murphy Columbia University The 34th International Symposium on Lattice Field Theory (Southampton, UK) July 27th, 2016 D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
2 Motivation Recent advances in algorithms used by RBC/UKQCD have driven down the cost of simulating degenerate light quark pairs: 1 Extensive tuning of forces via Hasenbusch mass splitting [Phys. Lett. B ] 2 Shamir DWF Möbius DWF zmöbius DWF [Izubuchi et al., Lattice 2015] 3 Reduced L s approximation during molecular dynamics [proposed by Brower et al., arxiv: ; implemented by G. McGlynn, Lattice 2015] 4 Mixed precision CG: single precision inner solves, double precision defect correction RHMC strange and charm quark determinants are now the most expensive part of our evolution strategy Gauge 5800 s 11.6% Light quarks (L s = 14) s 36.1% L s = 14/L s = 32 correction det s 3.3% Strange and charm quarks (L s = 32) s 48.5% DED 170 s 0.3% Total s Table: Timings for one HMC trajectory of N f = physical mass a 1 3 GeV ensemble on a 12,288-node BG/Q partition [G. McGlynn, Ph.D. thesis] D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
3 Motivation (cont d.) These techniques are of limited use for RHMC: Unclear how to restart inner single-precision solver in the context of multishift CG Cost of multishift D D inversions at light quark masses largely negates potential gain from Hasenbusch splitting For I = 0 K ππ calculations with G-parity boundary conditions D D describes 4 flavors, and RHMC is needed for light quarks as well [C. Kelly, 15:50] Goal: Explore the exact one-flavor algorithm as an alternative to RHMC for simulating single quark flavors or degenerate G-parity quark pairs D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
4 The Exact One-Flavor Algorithm Introduced by TWQCD for efficient one-flavor simulations on GPU clusters History: 2009: Introduction [TWQCD, arxiv: ] 2014: Detailed derivation [Chen and Chiu, Phys. Lett. B ] 2014: Demonstration of 20% overall speed-up in N f = simulation with EOFA vs. RHMC for heavy quark [Chen and Chiu, arxiv: ] Main idea: Use block manipulations in spin space to write [ ] D(m1) 1 det = D(m 2) det (M 1 L) det (M R) with M L and M R Hermitian and positive-definite Contrast with RHMC, where we instead compute [ ] { [ D(m1) D D(m 1) det = det D(m 2) D D(m 2) using a rational approximation to the square root ]} 1/2 Exact in the sense that one does not need to take a fractional power of the fermion determinant to remove unwanted flavors D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
5 Derivation Sketch I [Chen and Chiu, Phys. Lett. B ] We begin by factoring the DWF Dirac operator: D w: Wilson Dirac operator c,d: Möbius parameters L ss : 5D hopping term (D DWF ) xx,ss = ( (c + d) D w + 1 ) xx δ ss + ( (c d) D w 1 ) xx L ss { ( = (D w) xx δ ss + δ ) } { xx d + c (1 + L) (1 L) 1 1 ss }{{ } d (1 L) + c (1 } + L) D EOFA (D w) xx δ ss + δ xx (P +M + + P M ) ss M ± does not involve the gauge field and can be worked out explicitly Shamir kernel: (c, d) = (1/2, 1/2) D DWF and D EOFA are the same operator Möbius kernel: (c, d) = (α/2, 1/2) D EOFA is dense but has a somewhat better condition number than D DWF H γ 5R 5D EOFA is Hermitian, but not positive-definite D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
6 Derivation Sketch II [Chen and Chiu, Phys. Lett. B ] Consider a slight generalization of D EOFA, in the chiral representation of γ µ : ( ) W M 5 + M +(m 1) (σ t) D(m 1, m 2) = (σ t) W M 5 + M (m 2) Applying the Schur det. identity to D EOFA (m) = D(m, m), and taking a ratio: det (D EOFA (m 1)) det (W M5 + M+(m1)) det (H (m1)) = Schur complements det (D EOFA (m 2)) det (W M 5 + M (m 2)) det (H +(m 2)) Defining ± R 5 (M ±(m 2) M ±(m 1)) = kω ±Ω ± ( m2 m1), also have: det ( D(m1, m ) 2) det ( det (W M5 + M+(m1)) det (H+(m2) +) ) = 1 = = D(m1, m 2) det (W M 5 + M (m 2)) det (H (m 1) + ) Finally, using Sylvester s determinant identity (det(1 + AB) = det(1 + BA)): [ ( )] 1 [ ( det (D EOFA(m 1)) det (D = det 1 + kω 1 EOFA(m 2)) H (m 1) Ω det 1 + kω 1 + Ω + H +(m 2) +P + }{{}}{{} M L Hermitian and pos.-def.! M R )] 1 D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
7 Derivation Sketch III [Chen and Chiu, Phys. Lett. B ] To simulate, we express the determinant as a path integral: 1 det (M 1 L) det (M = Dφ LDφ L R) DφRDφ R e φ L M Lφ L φ R M Rφ R φ L and φ R are bosonic pseudofermion fields with two spin components M L and M R contain nested matrix inverses expensive to compute with CG Easier to work with block form: M EOFA = 1 kp Ω 1 H(m 1 ) Ω P + kp + Ω 1 + Ω + P + H(m 2 ) + P + for which we have, in terms of an ordinary four-component spinor: det (D EOFA (m 1)) det (D EOFA (m = DφDφ e φ M EOFA φ 2)) EOFA evolution requires inversions of the general form ( H(m) + α ±P ± ) ψ = φ D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
8 The Hybrid Monte Carlo Algorithm We perform dynamical QCD simulations using the HMC algorithm Idea: Generate a Markov chain of gauge field configurations by evolving a Hamiltonian system in (unphysical) molecular dynamics time Generalized coordinate: Gauge field (U x,µ) Conjugate momentum: π x,µ su(3) Hamiltonian: H = 1 2 π2 + S g[u] + S f [U] Equations of motion: { τ U x,µ = π x,µu x,µ ( ) τ π x,µ = T a x,µ a Sg[U] + S f [U] Metropolis accept/reject step corrects for finite precision integration Accept new gauge field U x,µ with probability Paccept = min ( 1, e H) Implementing this for EOFA requires: 1 Action S f [U] 2 Pseudofermion heatbath 3 Pseudofermion contribution to momentum update π x,µ π x,µ T a ( a x,µs f [U] ) τ Tests performed on two N f = Shamir DWF ensembles: 16I: , β = 2.13, (am l, am h ) = (0.01, 0.032) [Phys. Rev. D76, ] 32I: , β = 2.25, (am l, am h ) = (0.004, 0.03) [Phys. Rev. D93, ] D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
9 I. Action: Equivalence of EOFA and RHMC for Möbius DWF Relation between RHMC determinant ratios and EOFA determinant ratios: { [ ]} D 1/2 ( D DWF (m 1) (c + d) L s ) + m 1 (c d) Ls 12L 3 T [ ] DEOFA (m 1) det = det D D DWF (m 2) (c + d) Ls + m 2 (c d) Ls D EOFA (m 2) ] [ D(mq) D(mq+ mq) det Exact RHMC EOFA m q log(rw) RHMC EOFA N rw (a) Free 4 5 lattice, m q = 0.01, α = 4.0 (Möbius scale), 20 stochastic hits. (b) 16I. Reweight am h = in N rw steps with 10 hits per step. Observe significantly smaller systematic and statistical errors for EOFA Cheaper to perform quark mass reweighting with EOFA! D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
10 II. Heatbath At start of HMC trajectory: draw η with P (η) e η2 /2, compute φ = M 1/2 EOFA η Still requires rational approximation x 1/2 α 0 + N p α l=1 l/(β l + x) Defining γ l = (1 + β l ) 1, can show Np ] M 1/2 EOFA α01 + α lγ l [1 + kγ lp Ω 1 Ω P kγ lp +Ω 1 + Ω +P + H(m 1) γ l P H(m 2) γ lβ l +P + l=1 ±P ± has large number of zero modes not invertible! (no multishift) λ(meofa) 10 3 λ max 10 2 λ min am 1 (am 2 1.0) Figure: Spectral range of M EOFA. η η φ M φ / η η N p Figure: Heatbath relative error using N p poles. D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
11 III. Momentum Update EOFA pseudofermion force is { a x,µs f [U] = kχ L γ5r5 ( a x,µ D w ) χl kχ R γ5r5 ( a x,µ D w ) χr χ L = [H(m 1)] 1 Ω P φ, χ R = [H(m 2) +P +] 1 Ω +P +φ Compared to RHMC: Cheaper to evaluate Somewhat smaller total force, with F L F R Percentage EOFA RHMC F EOFA = F RHMC = Percentage EOFA (L) EOFA (R) F L = F R = F Figure: Lattice-wide distribution of EOFA and RHMC total forces by link F Figure: Lattice-wide distribution of EOFA force contributions from L and R terms by link D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
12 Test: Reproduce 16I Ensemble [Phys. Rev. D76, ] Parallel N f = evolutions: compare RHMC heavy quark to EOFA heavy quark P vs. MD traj. 1 H vs. MD traj. 0.2 Q (%) EOFA RHMC β 2.13 am l 0.01 am h Table: Parameters EOFA RHMC EOFA RHMC am π 0.243(2) 0.244(1) af π (8) (6) am K 0.326(2) 0.326(1) af K (6) (4) am Ω 0.990(9) 0.995(10) am res(m l) (4) (4) Table: Spectrum D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
13 Optimizations I: Accelerating EOFA Inversions Ref.: RHMC multishift EOFA (unpreconditioned) EOFA (preconditioned) EOFA (preconditioned + mixed CG) am phys s = am phys l = Time (s) EOFA solves accelerated with: Even-odd preconditioning Optimized BG/Q assembly generated by BAGEL [P. Boyle] Mixed-precision CG Test on 32I ensemble at physical am l and am s determined by chiral fits [Phys. Rev. D93, ] Reference (dashed line): even-odd preconditioned multishift solve of D D with same quark mass am phys s am phys l EOFA 7.6 s s RHMC 37.7 s s Ratio Table: Comparison of total CG inversion time between optimized EOFA and RHMC D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
14 Optimizations II: Accelerating EOFA Heatbath Following TWQCD, we use the chronological inversion method introduced by Brower et al. to forecast CG guesses [Nucl. Phys. B484, ] Forecast solution to (H + α l ±P ±) ψ = φ by minimizing Ψ[x] = x (H + α l ±P ±) x φ x x φ over the space of accumulated solutions for previous shifts α l CG Iterations L Zero Last Solution Forecasted R Pole Test: 16I, 10 poles in rational approximation Find 30% reduction in total iteration count No additional gain from sharing solutions between L and R terms D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
15 Conclusions and Next Steps We have: Independently implemented TWQCD s EOFA, and performed basic algorithmic tests Reproduced 16I ensemble using EOFA to evolve the heavy flavor Partially implemented and tested a number of optimizations: 1 Even-odd preconditioning 2 Optimized BG/Q assembly for sparse matrix applications (BAGEL) 3 Mixed-precision CG 4 Forecasted CG guesses for heatbath 5 Hasenbusch mass splitting Next steps: Better ideas to ameliorate the cost of the EOFA heatbath? Tuning Sexton-Weingarten integration and/or Hasenbusch mass splitting? G-parity simulations with EOFA light quarks Port to P. Boyle s Grid library: is a further speed-up possible if we perform deflated sparse matrix inversions with HDCR? [Poster by A. Yamaguchi and P. Boyle] We expect a gain over RHMC with a fully optimized and tuned implementation of EOFA, but the exact gain is difficult to predict! Thank you! D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
16 Extra Slides D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
17 Spectrum of the EOFA and RHMC actions Can think of EOFA as a preconditioning step applied to RHMC action Mapping of eigenvalue spectrum: (λ min, λ max) RHMC (1, λ max) EOFA Similar range, but M EOFA spectrum more densely concentrated near low end Computes same determinant ratio, but easier to stochastically estimate M EOFA! Notation: action is S = φ Mφ MEOFA = 1 kp Ω 1 H(m1) Ω P + 1 kp+ω + Ω+P+ H(m2) +P+ MRHMC = ( D D(m2) ) 1/4 ( D D(m1) ) 1/2 ( D D(m2) ) 1/ λ(m EOFA ) λ(m RHMC ) Figure: 4 5 free lattice, am 1 = 0.1, am 2 = 1.0, am 5 = 1.8. Lines mark 95% of density. D. Murphy (Lattice 2016) EOFA Domain Wall Fermion Simulations July 27th, / 17
Polynomial Filtered Hybrid Monte Carlo
Polynomial Filtered Hybrid Monte Carlo Waseem Kamleh and Michael J. Peardon CSSM & U NIVERSITY OF A DELAIDE QCDNA VII, July 4th 6th, 2012 Introduction Generating large, light quark dynamical gauge field
More informationLattice simulation of 2+1 flavors of overlap light quarks
Lattice simulation of 2+1 flavors of overlap light quarks JLQCD collaboration: S. Hashimoto,a,b,, S. Aoki c, H. Fukaya d, T. Kaneko a,b, H. Matsufuru a, J. Noaki a, T. Onogi e, N. Yamada a,b a High Energy
More informationThe Rational Hybrid Monte Carlo algorithm. M. A. Clark Center for Computational Science, Boston University, Boston, MA 02215, USA
Center for Computational Science, Boston University, Boston, MA 02215, USA E-mail: mikec@bu.edu The past few years have seen considerable progress in algorithmic development for the generation of gauge
More informationA construction of the Schrödinger Functional for Möbius Domain Wall Fermions
A construction of the Schrödinger Functional for Möbius Domain Wall Fermions Graduate school of Science, Hiroshima University, Higashi-Hiroshima, Japan E-mail: m132626@hiroshima-u.ac.jp Ken-Ichi Ishikawa
More informationPrecision Lattice Calculation of Kaon Decays with Möbius Domain Wall Fermions
Precision Lattice Calculation of Kaon Decays with Möbius Domain Wall Fermions Hantao Yin Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate School
More informationLattice QCD at non-zero temperature and density
Lattice QCD at non-zero temperature and density Frithjof Karsch Bielefeld University & Brookhaven National Laboratory QCD in a nutshell, non-perturbative physics, lattice-regularized QCD, Monte Carlo simulations
More informationThe phase diagram of QCD from imaginary chemical potentials
The phase diagram of QCD from imaginary chemical potentials Massimo D Elia Genoa University & INFN Quarks, Hadrons, and the Phase Diagram of QCD, St. Goar, september 3, 2009 In collaboration with Francesco
More informationChiral Fermions on the Lattice:
Chiral Fermions on the Lattice: A Flatlander s Ascent into Five Dimensions (ETH Zürich) with R. Edwards (JLab), B. Joó (JLab), A. D. Kennedy (Edinburgh), K. Orginos (JLab) LHP06, Jefferson Lab, Newport
More informationNumerical Time Integration in Hybrid Monte Carlo Schemes for Lattice Quantum Chromo Dynamics
Numerical Time Integration in Hybrid Monte Carlo Schemes for Lattice Quantum Chromo Dynamics M. Günther and M. Wandelt Bergische Universität Wuppertal IMACM Dortmund, March 12, 2018 Contents Lattice QCD
More informationopenq*d simulation code for QCD+QED
LATTICE 2017 Granada openq*d simulation code for QCD+QED Agostino Patella CERN & Plymouth University on behalf of the software-development team of the RC collaboration: Isabel Campos (IFCA, IFT), Patrick
More informationPoS(LATTICE 2013)487. Vacuum polarization function in N f = 2+1 domain-wall fermion. Eigo Shintani. Hyung-Jin Kim
Vacuum polarization function in N f = 2+1 domain-wall fermion PRISMA Cluster of Excellence, Institut für Kernphysik and Helmholtz Institute Mainz, Johannes Gutenberg-Universität Mainz, D-55099 Mainz, Germany
More informationarxiv:hep-lat/ v1 5 Oct 2006
KEK-CP-180 RIKEN-TH-77 UTHEP-57 JLQCD s dynamical overlap project arxiv:hep-lat/0610036v1 5 Oct 006 JLQCD Collaboration: a,b, S. Aoki c, H. Fukaya d, S. Hashimoto a,b, K-I. Ishikawa e, K. Kanaya c, H.
More informationSimulation of lattice QCD with domain-wall fermions
Journal of Physics: Conference Series OPEN ACCESS Simulation of lattice QCD with domain-wall fermions To cite this article: Ting-Wai Chiu and the Twqcd Collaboration 2013 J. Phys.: Conf. Ser. 454 012044
More informationarxiv: v1 [hep-lat] 7 Oct 2007
Charm and bottom heavy baryon mass spectrum from lattice QCD with 2+1 flavors arxiv:0710.1422v1 [hep-lat] 7 Oct 2007 and Steven Gottlieb Department of Physics, Indiana University, Bloomington, Indiana
More informationarxiv: v1 [hep-lat] 10 Jul 2012
Hybrid Monte Carlo with Wilson Dirac operator on the Fermi GPU Abhijit Chakrabarty Electra Design Automation, SDF Building, SaltLake Sec-V, Kolkata - 700091. Pushan Majumdar Dept. of Theoretical Physics,
More informationarxiv: v1 [hep-lat] 28 Aug 2014 Karthee Sivalingam
Clover Action for Blue Gene-Q and Iterative solvers for DWF arxiv:18.687v1 [hep-lat] 8 Aug 1 Department of Meteorology, University of Reading, Reading, UK E-mail: K.Sivalingam@reading.ac.uk P.A. Boyle
More informationarxiv: v1 [hep-lat] 23 Dec 2010
arxiv:2.568v [hep-lat] 23 Dec 2 C. Alexandrou Department of Physics, University of Cyprus, P.O. Box 2537, 678 Nicosia, Cyprus and Computation-based Science and Technology Research Center, Cyprus Institute,
More informationProbing the Chiral Limit in 2+1 flavor Domain Wall Fermion QCD
Probing the Chiral Limit in 2+1 flavor Domain Wall Fermion QCD Meifeng Lin for the RBC and UKQCD Collaborations Department of Physics Columbia University July 29 - August 4, 2007 / Lattice 2007 @ Regensburg
More informationT.W. Chiu, Chung-Yuan Christian Univ, May 13, 2008 p.1/34. The Topology in QCD. Ting-Wai Chiu Physics Department, National Taiwan University
T.W. Chiu, Chung-Yuan Christian Univ, May 13, 2008 p.1/34 The Topology in QCD Ting-Wai Chiu Physics Department, National Taiwan University The vacuum of QCD has a non-trivial topological structure. T.W.
More informationarxiv: v1 [hep-lat] 12 Nov 2014
Exploring the phase structure of 12-flavor SU(3) arxiv:111.336v1 [hep-lat] 12 Nov 21 Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa, USA E-mail: zechariah-gelzer@uiowa.edu Yuzhi
More informationNucleon structure from 2+1-flavor dynamical DWF ensembles
Nucleon structure from 2+1-flavor dynamical DWF ensembles Michael Abramczyk Department of Physics, University of Connecticut, Storrs, CT 06269, USA E-mail: michael.abramczyk@uconn.edu Meifeng Lin Computational
More informationTopological susceptibility in (2+1)-flavor lattice QCD with overlap fermion
T.W. Chiu, Lattice 2008, July 15, 2008 p.1/30 Topological susceptibility in (2+1)-flavor lattice QCD with overlap fermion Ting-Wai Chiu Physics Department, National Taiwan University Collaborators: S.
More informationarxiv: v2 [hep-lat] 23 Dec 2008
arxiv:8.964v2 [hep-lat] 23 Dec 28, F. Farchioni, A. Ferling, G. Münster, J. Wuilloud University of Münster, Institute for Theoretical Physics Wilhelm-Klemm-Strasse 9, D-4849 Münster, Germany E-mail: k_demm@uni-muenster.de
More informationCurrent Physics Projects by JLQCD
Current Physics Projects by JLQCD Jun Noaki 野秋 淳一 for JLQCD Collaboration Lattice Hadron Physics V, Cairns July 20 24, 2015 Geography JLQCD members (as of Jul. 2015) Tsukuba (KEK) : G. Cossu, B. Fahy,
More informationAn application of the UV-filtering preconditioner to the Polynomial Hybrid Monte Carlo algorithm
An application of the UV-filtering preconditioner to the Polynomial Hybrid Monte Carlo algorithm PACS-CS Collaboration: a, S. Aoki b,c, T. Ishikawa d, N. Ishizuka b,d, K. Kanaya b, Y. Kuramashi b,d, M.
More informationIsospin and Electromagnetism
Extreme Scale Computing Workshop, December 9 11, 2008 p. 1/11 Isospin and Electromagnetism Steven Gottlieb Extreme Scale Computing Workshop, December 9 11, 2008 p. 2/11 Questions In the exascale era, for
More informationMass Components of Mesons from Lattice QCD
Mass Components of Mesons from Lattice QCD Ying Chen In collaborating with: Y.-B. Yang, M. Gong, K.-F. Liu, T. Draper, Z. Liu, J.-P. Ma, etc. Peking University, Nov. 28, 2013 Outline I. Motivation II.
More informationThe Conformal Window in SU(3) Yang-Mills
The Conformal Window in SU(3) Yang-Mills Ethan T. Neil ethan.neil@yale.edu Department of Physics Yale University Lattice 2008 Williamsburg, VA July 15, 2008 Ethan Neil (Yale) Conformal Window in Yang-Mills
More informationThe heavy-light sector of N f = twisted mass lattice QCD
The heavy-light sector of N f = 2 + 1 + 1 twisted mass lattice QCD Marc Wagner Humboldt-Universität zu Berlin, Institut für Physik mcwagner@physik.hu-berlin.de http://people.physik.hu-berlin.de/ mcwagner/
More informationThe Polyakov Loop and the Eigenvalues of the Dirac Operator
The Polyakov Loop and the Eigenvalues of the Dirac Operator Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA E-mail: soeldner@bnl.gov Aiming at the link between confinement and
More informationMobility edge and locality of the overlap-dirac operator with and without dynamical overlap fermions
Mobility edge and locality of the overlap-dirac operator with and without dynamical overlap fermions JLQCD Collaboration: a,b, S. Aoki c,d, H. Fukaya e, S. Hashimoto a,b, K-I. Ishikawa f, K. Kanaya c,
More informationChiral fermions and chemical potential
Chiral fermions and chemical potential Some thoughts Rajamani Narayanan Department of Physics Florida International University NCSU, July 21 Introduction of chemical potential on the lattice for chiral
More informationTesting a Fourier Accelerated Hybrid Monte Carlo Algorithm
Syracuse University SURFACE Physics College of Arts and Sciences 12-17-2001 Testing a Fourier Accelerated Hybrid Monte Carlo Algorithm Simon Catterall Syracuse University Sergey Karamov Syracuse University
More informationThermodynamics using p4-improved staggered fermion action on QCDOC
Thermodynamics using p4-improved staggered fermion action on QCDOC for the RBC-Bielefeld Collaboration Brookhaven National Laboratory and Columbia University, USA E-mail: chulwoo@bnl.gov We present an
More informationPoS(EPS-HEP 2009)057. Bottomonium Studies at BaBar. Veronique Ziegler. SLAC National Accelerator Laboratory
SLAC National Accelerator Laboratory E-mail: vziegler@slac.stanford.edu Selected studies in bottomonium physics carried out by the BaBar experiment at the SLAC PEP-II e + e collider are presented. They
More informationCold QCD. Meeting on Computational Nuclear Physics. Washington, DC July Thomas Luu Lawrence Livermore National Laboratory
Cold QCD Meeting on Computational Nuclear Physics Washington, DC July 2012 Thomas Luu Lawrence Livermore National Laboratory This work was performed under the auspices of the U.S. Department of Energy
More informationCase Study: Quantum Chromodynamics
Case Study: Quantum Chromodynamics Michael Clark Harvard University with R. Babich, K. Barros, R. Brower, J. Chen and C. Rebbi Outline Primer to QCD QCD on a GPU Mixed Precision Solvers Multigrid solver
More informationEfficient implementation of the overlap operator on multi-gpus
Efficient implementation of the overlap operator on multi-gpus Andrei Alexandru Mike Lujan, Craig Pelissier, Ben Gamari, Frank Lee SAAHPC 2011 - University of Tennessee Outline Motivation Overlap operator
More informationAccelerating Quantum Chromodynamics Calculations with GPUs
Accelerating Quantum Chromodynamics Calculations with GPUs Guochun Shi, Steven Gottlieb, Aaron Torok, Volodymyr Kindratenko NCSA & Indiana University National Center for Supercomputing Applications University
More informationFun with the S parameter on the lattice
Fun with the S parameter on the lattice David Schaich (Boston Colorado Syracuse) from work with the LSD Collaboration and USQCD BSM community arxiv:1009.5967 & arxiv:1111.4993 Origin of Mass 2013 Lattice
More informationParallelization of the Dirac operator. Pushan Majumdar. Indian Association for the Cultivation of Sciences, Jadavpur, Kolkata
Parallelization of the Dirac operator Pushan Majumdar Indian Association for the Cultivation of Sciences, Jadavpur, Kolkata Outline Introduction Algorithms Parallelization Comparison of performances Conclusions
More informationBaryon Resonance Determination using LQCD. Robert Edwards Jefferson Lab. Baryons 2013
Baryon Resonance Determination using LQCD Robert Edwards Jefferson Lab Baryons 2013 Where are the Missing Baryon Resonances? What are collective modes? Is there freezing of degrees of freedom? What is
More informationChiral magnetic effect in 2+1 flavor QCD+QED
M. Abramczyk E-mail: mabramc@gmail.com E-mail: tblum@phys.uconn.edu G. Petropoulos E-mail: gregpetrop@gmail.com R. Zhou, E-mail: zhouran13@gmail.com Physics Department, University of Connecticut, 15 Hillside
More informationThe electric dipole moment of the nucleon from lattice QCD with imaginary vacuum angle theta
The electric dipole moment of the nucleon from lattice QCD with imaginary vacuum angle theta Yoshifumi Nakamura(NIC/DESY) for the theta collaboration S. Aoki(RBRC/Tsukuba), R. Horsley(Edinburgh), YN, D.
More informationMixed action simulations: approaching physical quark masses
Mixed action simulations: approaching physical quark masses Stefan Krieg NIC Forschungszentrum Jülich, Wuppertal University in collaboration with S. Durr, Z. Fodor, C. Hoelbling, S. Katz, T. Kurth, L.
More informationPROTON DECAY MATRIX ELEMENTS FROM LATTICE QCD. Yasumichi Aoki RIKEN BNL Research Center. 9/23/09 LBV09 at Madison
PROTON DECAY MATRIX ELEMENTS FROM LATTICE QCD Yasumichi Aoki RIKEN BNL Research Center 9/23/09 LBV09 at Madison Plan low energy matrix elements for N PS,l GUT QCD relation what is exactly needed to calculate
More informationAxial symmetry in the chiral symmetric phase
Axial symmetry in the chiral symmetric phase Swagato Mukherjee June 2014, Stoney Brook, USA Axial symmetry in QCD massless QCD Lagrangian is invariant under U A (1) : ψ (x) e i α ( x) γ 5 ψ(x) μ J 5 μ
More informationPoS(LAT2006)208. Diseases with rooted staggered quarks. Michael Creutz Brookhaven National Laboratory, Upton, NY 11973, USA
Brookhaven National Laboratory, Upton, NY 11973, USA E-mail: creutz@bnl.gov Calculations using staggered quarks augmented with a root of the fermion determinant to reduce doubling give a qualitatively
More informationThe Λ(1405) is an anti-kaon nucleon molecule. Jonathan Hall, Waseem Kamleh, Derek Leinweber, Ben Menadue, Ben Owen, Tony Thomas, Ross Young
The Λ(1405) is an anti-kaon nucleon molecule Jonathan Hall, Waseem Kamleh, Derek Leinweber, Ben Menadue, Ben Owen, Tony Thomas, Ross Young The Λ(1405) The Λ(1405) is the lowest-lying odd-parity state of
More informationCascades on the Lattice
Cascade Physics - Jlab 2005 Cascades on the Lattice Kostas Orginos College of William and Mary - JLab LHP Collaboration LHPC collaborators R. Edwards (Jlab) G. Fleming (Yale) P. Hagler (Vrije Universiteit)
More informationThe Λ(1405) is an anti-kaon nucleon molecule. Jonathan Hall, Waseem Kamleh, Derek Leinweber, Ben Menadue, Ben Owen, Tony Thomas, Ross Young
The Λ(1405) is an anti-kaon nucleon molecule Jonathan Hall, Waseem Kamleh, Derek Leinweber, Ben Menadue, Ben Owen, Tony Thomas, Ross Young The Λ(1405) The Λ(1405) is the lowest-lying odd-parity state of
More informationHierarchically Deflated Conjugate Gradient
Hierarchically Deflated Conjugate Gradient Peter Boyle, University of Edinburgh July 29, 2013 Eigenvector Deflation Krylov solvers convergence controlled by the condition number κ λmax λ min Lattice chiral
More informationImproving many flavor QCD simulations using multiple GPUs
Improving many flavor QCD simulations using multiple GPUs M. Hayakawa a, b, Y. Osaki b, S. Takeda c, S. Uno a, N. Yamada de a Department of Physics, Nagoya University, Nagoya 464-8602, Japan b Department
More informationCAA: Covariant Approximation Averaging. a new class of error reduction techniques
CAA: Covariant Approximation Averaging a new class of error reduction techniques Taku Izubuchi, with T. Blum E. Shintani, C. Lehner, T. Kawanai, T. Ishikawa, J.Yu, R. Arthur, P. Bolye,. RBC/UKQCD in preparation
More informationA study of chiral symmetry in quenched QCD using the. Overlap-Dirac operator
FSU-SCRI-98-128 A study of chiral symmetry in quenched QCD using the Overlap-Dirac operator Robert G. Edwards, Urs M. Heller, Rajamani Narayanan SCRI, Florida State University, Tallahassee, FL 32306-4130,
More informationRG scaling at chiral phase transition in two-flavor QCD
RG scaling at chiral phase transition in two-flavor QCD Graduate School of Science, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-826, Japan E-mail: ishikawa@theo.phys.sci.hiroshima-u.ac.jp Yoichi
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Takustrasse 7, Berlin, Germany
DECORRELATION OF THE TOPOLOGICAL CHARGE IN TEMPERED SIMULATIONS OF FULL QCD H. STÜBEN Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustrasse 7, 14195 Berlin, Germany Abstract. The improvement of
More informationPoS(LATTICE 2013)500. Charmonium, D s and D s from overlap fermion on domain wall fermion configurations
Charmonium, D s and D s from overlap fermion on domain wall fermion configurations,, Y. Chen, A. Alexandru, S.J. Dong, T. Draper, M. Gong,, F.X. Lee, A. Li, 4 K.F. Liu, Z. Liu, M. Lujan, and N. Mathur
More informationLattice Monte Carlo for carbon nanostructures. Timo A. Lähde. In collaboration with Thomas Luu (FZ Jülich)
Lattice Monte Carlo for carbon nanostructures Timo A. Lähde In collaboration with Thomas Luu (FZ Jülich) Institute for Advanced Simulation and Institut für Kernphysik Forschungszentrum Jülich GmbH, D-52425
More informationof the gauge eld U which one obtains after integrating through K = (8 + 2m) 1. It is easy to show that the lattice Dirac operator
A new simulation algorithm for lattice QCD with dynamical quarks B. Bunk a, K. Jansen b, B. Jegerlehner c, M. Luscher b, H. Simma b and R. Sommer d a Humboldt Universitat, Institut fur Physik Invalidenstrasse
More informationarxiv: v1 [hep-lat] 4 Nov 2014
Meson Mass Decomposition,2, Ying Chen, Terrence Draper 2, Ming Gong,2, Keh-Fei Liu 2, Zhaofeng Liu, and Jian-Ping Ma 3,4 arxiv:4.927v [hep-lat] 4 Nov 24 (χqcd Collaboration) Institute of High Energy Physics,
More informationThe Fermion Bag Approach
The Fermion Bag Approach Anyi Li Duke University In collaboration with Shailesh Chandrasekharan 1 Motivation Monte Carlo simulation Sign problem Fermion sign problem Solutions to the sign problem Fermion
More informationCHARMED BOTTOM BARYON SPECTROSCOPY. Zachary S. Brown, William Detmold, Stefan Meinel, Konstantinos Orginos
CHARMED BOTTOM BARYON SPECTROSCOPY Zachary S. Brown, William Detmold, Stefan Meinel, Konstantinos Orginos 1 OUTLINE Landscape of heavy baryon spectroscopy Details of our calculation Extrapolations Results
More informationFirst results from dynamical chirally improved fermions
First results from dynamical chirally improved fermions arxiv:hep-lat/595v1 1 Sep 25 C. B.Lang Karl-Franzens-Universität Graz, Austria E-mail: christian.lang@uni-graz.at Karl-Franzens-Universität Graz,
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationLattice QCD calculation of direct CP violation and long distance effects in kaon mixing and rare decays
Lattice QCD calculation of direct CP violation and long distance effects in kaon mixing and rare decays Department of Physics, Columbia University, New York, NY 10027, USA E-mail: nhc@phys.columbia.edu
More informationDDalphaAMG for Twisted Mass Fermions
DDalphaAMG for Twisted Mass Fermions, a, b, Constantia Alexandrou a,c, Jacob Finkenrath c, Andreas Frommer b, Karsten Kahl b and Matthias Rottmann b a Department of Physics, University of Cyprus, PO Box
More informationMultigrid Methods for Linear Systems with Stochastic Entries Arising in Lattice QCD. Andreas Frommer
Methods for Linear Systems with Stochastic Entries Arising in Lattice QCD Andreas Frommer Collaborators The Dirac operator James Brannick, Penn State University Björn Leder, Humboldt Universität Berlin
More informationMonte Carlo Algorithms
Last Modified Wednesday, 15 June 011 Monte Carlo Algorithms Tait Institute, School of Physics & Astronomy, The University of Edinburgh Tuesday, 14 June 011 STRONGnet Summer School 011, Bielefeld Outline
More informationη and η mesons from N f = flavour lattice QCD
η and η mesons from N f = 2++ flavour lattice QCD for the ETM collaboration C. Michael, K. Ottnad 2, C. Urbach 2, F. Zimmermann 2 arxiv:26.679 Theoretical Physics Division, Department of Mathematical Sciences,
More informationA Lattice Path Integral for Supersymmetric Quantum Mechanics
Syracuse University SURFACE Physics College of Arts and Sciences 8-3-2000 A Lattice Path Integral for Supersymmetric Quantum Mechanics Simon Catterall Syracuse University Eric B. Gregory Zhongshan University
More informationLattice Quantum Chromodynamics on the MIC architectures
Lattice Quantum Chromodynamics on the MIC architectures Piotr Korcyl Universität Regensburg Intel MIC Programming Workshop @ LRZ 28 June 2017 Piotr Korcyl Lattice Quantum Chromodynamics on the MIC 1/ 25
More informationlattice QCD and the hadron spectrum Jozef Dudek ODU/JLab
lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab a black box? QCD lattice QCD observables (scattering amplitudes?) in these lectures, hope to give you a look inside the box 2 these lectures how
More informationThe light hadron spectrum from lattice QCD
Ch. Hoelbling with S. Durr Z. Fodor J. Frison S. Katz S. Krieg T. Kurth L. Lellouch Th. Lippert K. Szabo G. Vulvert The Budapest-Marseille-Wuppertal collaboration Rab 2008 Outline 1 Motivation 2 Setup
More informationTowards QCD Thermodynamics using Exact Chiral Symmetry on Lattice
Towards QCD Thermodynamics using Exact Chiral Symmetry on Lattice Debasish Banerjee, Rajiv V. Gavai & Sayantan Sharma T. I. F. R., Mumbai arxiv : 0803.3925, to appear in Phys. Rev. D, & in preparation.
More informationHadronic physics from the lattice
Hadronic physics from the lattice Chris Michael c.michael@liv.ac.uk University of Liverpool Hadronic physics from the lattice p.1/24 Hadronic Structure - Introduction What are hadrons made of? Is a meson
More informationLATTICE 4 BEGINNERS. Guillermo Breto Rangel May 14th, Monday, May 14, 12
LATTICE 4 BEGINNERS Guillermo Breto Rangel May 14th, 2012 1 QCD GENERAL 2 QCD GENERAL 3 QCD vs QED QCD versus QED Quantum Electrodynamics (QED): The interaction is due to the exchange of photons. Every
More informationSpectroscopy of charmed baryons from lattice QCD
Institute of Physics, University of Graz, 8010 Graz, Austria. E-mail: padmanath.madanagopalan@uni-graz.at Robert G. Edwards Jefferson Laboratory, 12000 Jefferson Avenue, Newport News, VA 23606, USA E-mail:
More informationarxiv: v1 [hep-lat] 19 Jul 2009
arxiv:0907.3261v1 [hep-lat] 19 Jul 2009 Application of preconditioned block BiCGGR to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD Abstract H. Tadano a,b, Y. Kuramashi c,b, T.
More informationPseudoscalar Flavor-Singlet Physics with Staggered Fermions
Pseudoscalar Flavor-Singlet Physics with Staggered Fermions UKQCD Collaboration, Alan Irving, Chris M. Richards Department of Mathematical Sciences, University of Liverpool, Liverpool, L69-7ZL, UK E-mail:
More informationLattice Simulations with Chiral Nuclear Forces
Lattice Simulations with Chiral Nuclear Forces Hermann Krebs FZ Jülich & Universität Bonn July 23, 2008, XQCD 2008, NCSU In collaboration with B. Borasoy, E. Epelbaum, D. Lee, U. Meißner Outline EFT and
More informationQCD thermodynamics with two-flavours of Wilson fermions on large lattices
QCD thermodynamics with two-flavours of Wilson fermions on large lattices Bastian Brandt Institute for nuclear physics In collaboration with A. Francis, H.B. Meyer, O. Philipsen (Frankfurt) and H. Wittig
More informationLattice QCD + Hydro/Cascade Model of Heavy Ion Collisions
Lattice QCD + Hydro/Cascade Model of Heavy Ion Collisions Michael Cheng Lawrence Livermore National Laboratory 2010 Winter Workshop on Nuclear Dynamics Ocho Rios, Jamaica, January 2-9, 2010 Outline Calculation
More informationarxiv: v1 [hep-ph] 13 Nov 2013
euniversityofmanchester November 14, 2013 arxiv:1311.3203v1 [hep-ph] 13 Nov 2013 Odd- and even-parity charmed mesons revisited in heavy hadron chiral perturbation theory Mohammad Alhakami 1 School of Physics
More informationlattice QCD and the hadron spectrum Jozef Dudek ODU/JLab
lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab the light meson spectrum relatively simple models of hadrons: bound states of constituent quarks and antiquarks the quark model empirical meson
More informationDefining Chiral Gauge Theories Beyond Perturbation Theory
Defining Chiral Gauge Theories Beyond Perturbation Theory Lattice Regulating Chiral Gauge Theories Dorota M Grabowska UC Berkeley Work done with David B. Kaplan: Phys. Rev. Lett. 116 (2016), no. 21 211602
More informationFermions in higher representations. Some results about SU(2) with adjoint fermions.
Fermions in higher representations. Some results about SU(2) with adjoint fermions. Agostino Patella Swansea University with L. Del Debbio, C. Pica arxiv:0805.2058 [hep-lat] Lattice08, Williamsburg 18/7/2008
More informationLow-energy QCD II Status of Lattice Calculations
Low-energy QCD II Status of Lattice Calculations Hartmut Wittig Institute for Nuclear Physics and Helmholtz Institute Mainz Determination of the Fundamental Parameters of QCD Nanyang Technical University
More informationCritical end point of Nf=3 QCD at finite temperature and density
Critical end point of Nf=3 QCD at finite temperature and density a,b, Xiao-Yong Jin b, Yoshinobu Kuramashi b,c,d, Yoshifumi Nakamura b, and Akira Ukawa b a Institute of Physics, Kanazawa University, Kanazawa
More informationB-meson decay constants with domain-wall light quarks and nonperturbatively tuned relativistic b-quarks
B-meson decay constants with domain-wall light quarks and nonperturbatively tuned relativistic b-quarks RBC and UKQCD collaborations Oliver Witzel Center for Computational Science Lattice 2013, Mainz,
More informationLattice QCD with Eight Degenerate Quark Flavors
Lattice QCD with Eight Degenerate Quark Flavors Xiao-Yong Jin, Robert D. Mawhinney Columbia University Lattice 2008 Outline Introduction Simulations and results Preparations Results Conclusion and outlook
More informationKaon Physics in the Standard Model
Kaon Physics in the Standard Model RBC and UKQCD Collaborations USQCD All Hands Meeting April 28, 2017 Norman H. Christ Columbia University Five Topics 2 nd order weak with charm Long distance contribution
More informationSpectral action, scale anomaly. and the Higgs-Dilaton potential
Spectral action, scale anomaly and the Higgs-Dilaton potential Fedele Lizzi Università di Napoli Federico II Work in collaboration with A.A. Andrianov (St. Petersburg) and M.A. Kurkov (Napoli) JHEP 1005:057,2010
More informationSpectral Properties of Quarks in the Quark-Gluon Plasma
Lattice27 : 2, Aug., 27 Spectral Properties of Quarks in the Quark-Gluon Plasma Masakiyo Kitazawa (Osaka Univ.) F. Karsch and M.K., arxiv:78.299 Why Quark? Because there are quarks. in the deconfined phase
More informationMass of Heavy Mesons from Lattice QCD
Mass of Heavy Mesons from Lattice QCD David Richards Jefferson Laboratory/Hadron Spectrum Collaboration Temple, March 2016 Outline Heavy Mesons Lattice QCD Spectroscopy Recipe Book Results and insight
More informationLQCD at non-zero temperature : strongly interacting matter at high temperatures and densities Péter Petreczky
LQCD at non-zero temperature : strongly interacting matter at high temperatures and densities Péter Petreczky QCD and hot and dense matter Lattice formulation of QCD Deconfinement transition in QCD : EoS
More informationRG scaling at chiral phase transition! in two-flavor QCD
Lattice 2016 July 25-30 2016 RG scaling at chiral phase transition in two-flavor QCD K.-I.Ishikawa(Hiroshima U.) Y. Iwasaki(U.Tsukuba) Yu Nakayama(CalTech) T. Yoshie(U. Tsukuba) 1 Objective Clarify the
More informationHeavy-quark hybrid mesons and the Born-Oppenheimer approximation
Heavy-quark hybrid mesons and the Born-Oppenheimer approximation Colin Morningstar Carnegie Mellon University Quarkonium Workshop, Fermilab Sept 20, 2003 9/20/2003 Hybrid mesons (C. Morningstar) 1 Outline!
More informationExploring composite Higgs scenarios with mass-split models
Exploring composite Higgs scenarios with mass-split models Oliver Witzel Lattice for BSM Physics 2018 Boulder, CO, USA, April 05, 2018 1 / 22 motivation Experimental observations Discovery of the Higgs
More informationNeutrinos and flavour in a brane model
Neutrinos and flavour in a brane model Raymond R. Volkas The University of Melbourne Neutrinos in Cosmology, in Astro-, Particle- and Nuclear Physics September 2009 Raymond R. Volkas (U Melbourne) Neutrinos
More information