Numerical Time Integration in Hybrid Monte Carlo Schemes for Lattice Quantum Chromo Dynamics
|
|
- Cornelia Lynch
- 5 years ago
- Views:
Transcription
1 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
2 Contents Lattice QCD Hybrid Monte Carlo Molecular Dynamics Step Geometric Integration Fermionic fields M. Günther and M. Wandelt, HMC for Lattice QCD 1/30
3 Physical background Quantum Chromo Dynamics (QCD) QCD is the theory of the strong interaction (=color force) between quarks and gluons inside subatomic particles. Quark = elementary particle, fundamental constituent of matter, e.g., protons and neutrons. Gluons = exchange particles for the strong force between quarks. Lattice QCD approach to solve the theory of QCD investigation of elementary particles using computer simulations formulated on a lattice of points in space and time links U µ (x) SU(3) between points x and x + aˆµ M. Günther and M. Wandelt, HMC for Lattice QCD 2/30
4 Task of Lattice QCD Compuation of observables (particle mass, e.g.) observables given as expectation of certain operators O: O = O(U)Π(U)dU with probability space (Ω, σ(ω), F Π ) Ω = (SU(3)) N, N = lattice size Ω probability density Π(U) = exp( S(U))/Z with Z = Ω exp( S(U))dU, action S acting on U Ω. Problem: integral has very high dimension N = , for example M. Günther and M. Wandelt, HMC for Lattice QCD 3/30
5 Importance sampling Approximate integral by O 1 n n O(U i ) i=1 choose configurations U i according to the probability density exp( S(U))/Z prefer configurations that occur with a higher probability sequence is generated by a Markov chain is ergodic satisfies the detailed balance condition = convergence to the unique fixed point distribution given by the density Π(U) M. Günther and M. Wandelt, HMC for Lattice QCD 4/30
6 Hybrid Monte Carlo Method of choice in Lattice QCD uses augmented Markov chain constructs samples of pairs of links U Ω and momenta P ˆΩ = (su(3)) N samples generated according to the probability density ν(u, P ) given by ν(u, P ) = (2π)N/2 1 exp( S(U)) exp( P, P /2) Z } H (2π) N/2 {{}}{{} :=Π(U) :=ϕ(p ) 1 = exp( H(U, P )) with Z H Z H = ν(u, P )d(u, P ), H(U, P ) = P, P /2 + S(U) Ω ˆΩ M. Günther and M. Wandelt, HMC for Lattice QCD 5/30
7 Transition (U 0, P 0 ) (U, P ) Proposal step (U 0, P 0 ) g(u 0, P 0 ) for the moment: arbitrary mapping g : Ω ˆΩ Ω ˆΩ Acceptance step (U, P ) := g(u 0, P 0 ) g(u 0, P 0 ) is accepted with probability ( α((u 0, P 0 )g(u 0, P 0 )) = min 1, ν(g(u ) 0, P 0 )) ν(u 0, P 0 ) = min ( 1, e ) (H(g(U 0,P 0 )) H(U 0,P 0 )) otherwise the old configuration is kept: (U, P ) := (U 0, P 0 ) M. Günther and M. Wandelt, HMC for Lattice QCD 6/30
8 Demand on g Detailed balance condition holds if g is time-reversible: with S = S g(s g(u 0, P 0 ) ) = (U 0, P 0 ) ( ) I 0 flipping the momenta 0 I g is volume-preserving: ( ) det g(u0, P 0 ) = 1. (U 0, P 0 ) Choice of g? M. Günther The and M. Hamiltonian Wandelt, HMC for Lattice flow! QCD 7/30
9 Detailed Balance condition I Transition kernel in HMC K((U 0, P 0 ), (Ω ˆΩ )) = P ((U, P ) (Ω ˆΩ ) (U 0, P 0 )) = α((u 0, P 0 ), g(u 0, P 0 ))δ g(u0,p 0)(Ω ˆΩ ) + (1 α((u 0, P 0 ), g(u 0, P 0 )))δ (U0,P 0)(Ω ˆΩ ). Restricted Transition kernel K L (U 0, Ω ) = P (U 0 Ω U 0 ) = K((U 0, P ), (Ω ˆΩ))ϕ(P )dp ˆΩ Detailed Balance condition for all A, B σ(ω) K L (U, B)Π(U)dU = K L (U, A)Π(U)dU A B M. Günther and M. Wandelt, HMC for Lattice QCD 8/30
10 Detailed Balance condition II Detailed Balance condition for all A, B σ(ω) K L (U, B)Π(U)dU = K L (U, A)Π(U)dU A Left-hand side ( ) ν(g(u, P )) min 1, δ g(u,p ) (B A ˆΩ ν(u, P ) ˆΩ)ϕ(P )Π(U)dP du + ( ( )) ν(g(u, P )) 1 min 1, ϕ(p )Π(U)dP du. ν(u, P ) A B ˆΩ Right-hand side ( ) ν(g(u, P )) min 1, δ g(u,p ) (A B ˆΩ ν(u, P ) ˆΩ)ϕ(P )Π(U)dP du + ( ( )) ν(g(u, P )) 1 min 1, ϕ(p )Π(U)dP du. ν(u, P ) B A ˆΩ M. Günther and M. Wandelt, HMC for Lattice QCD 9/30 B
11 Detailed Balance condition III To be shown min (ν(u, P ), ν(g(u, P ))) δ g(u,p ) (B ˆΩ)d(P, U) = A ˆΩ min (ν(u, P ), ν(g(u, P ))) δ g(u,p ) (A ˆΩ)d(P, U). B ˆΩ With δ g(u,p ) (B ˆΩ) = δ (U,P ) (g 1 )(B ˆΩ)) min (ν(u, P ), ν(g(u, P ))) d(p, U) = (A ˆΩ) g 1 (B ˆΩ) }{{} (B ˆΩ) g 1 (A ˆΩ) F 1= min (ν(u, P ), ν(g(u, P ))) d(p, U) } {{ } F 2= M. Günther and M. Wandelt, HMC for Lattice QCD 10/30
12 Detailed Balance condition IV ϕ(p ) = ϕ( P ), ν(u, P ) = ν(u, P ) & g time reversible ν(g(u, P )) = ν(g 1 (U, P )) & g(u, P ) = g 1 (U, P ) Rewrite F 1 as F 1 = min (ν(u, P ), ν(g(u, P ))) δ (U,P ) (A ˆΩ)d(P, U) = g 1 (B ˆΩ) ( min ν(g 1 ) (U, P )), ν(u, P ) δ g 1 (U,P ) (A ˆΩ) B ˆΩ det dg 1 (U, P ) d(p, U) = d(u, P ) min (ν(g(u, P )), ν(u, P )) δ g 1 (U,P ) (A ˆΩ) B ˆΩ det dg 1 (U, P ) d(p, U) = d(u, P ) min (ν(g(u, P )), ν(u, P )) δ g(u, P ) (A ˆΩ) B ˆΩ det dg 1 (U, P ) d(p, U) = d(u, P ) min (ν(g(u, P )), ν(u, P )) (B ˆΩ) g 1 (A ˆΩ) det dg 1 (U, P ) d(p, U) = F d(u, P ) 2 } {{ } = 1 M. Günther and M. Wandelt, HMC for Lattice QCD 11/30
13 Choice of g Demand on g: time-reversible and volume-preserving mapping Natural choice: use Hamiltonian flow H(U, P ) = P, P /2 + S(U) Problem: analytical flow not available Way-out: use numerical approximation using geometric integration accuracy of less importance (only to get high acceptance rate) Demand on numerical approximation of g time-reversible & volume preserving numerical flow compatible with non-abelian structure of U (SU(3)) N and P (su(3)) N M. Günther and M. Wandelt, HMC for Lattice QCD 12/30
14 Derivation of Equation of Motion Time Derivative of Links U µ (x) SU(3) derived via an infinitesimal rotation on the group manifold U µ (x) = P µ (x)u µ (x) Time Derivative of Momenta P µ (x) su(3) usually: solve Ḣ(U, P ) = 0 for P µ (x) via the link differential operator P µ (x) = x,µ H(U, P ) = x,µ S(U) x,µf(u) = T i i x,µf(u) with i x,µf(u) = df(us) ds s=0 { (U s) ν(y) = e st i U µ(x) if (y, ν) = (x, µ) U ν(y) otherwise. M. Günther and M. Wandelt, HMC for Lattice QCD 13/30
15 Example: Wilson action S g (U) = β Re tr(u N µ(x)σ µ(x)) + terms indep. of U µ (x) Staples Σ µ (x) defined by [ Uν(x)U µ(x + aˆν)u ν(x + aˆµ) 1 + U ν(x aˆν) 1 U µ(x aˆν)u ν(x aˆν + aˆµ) ] ν µ i x,µ Sg(U) = β N Re tr { T i U µ(x)σ µ (x) } x,µs g(u) = β { } N T i Re tr T i U µ(x)σ µ(x) = β 2N {Uµ(x)Σ µ(x)} T A Time Derivative of Momenta P µ (x) su(3) P µ (x) = x,µ S g (U) = β 2N {U µ(x)σ µ(x)} T A M. Günther and M. Wandelt, HMC for Lattice QCD 14/30
16 Numerical integration on Lie groups Structure of ODE-IVP U = diag(p µ (x)u µ (x)), U(t 0 ) = U 0 P = x,µ S(U), P (t 0 ) = P 0 with U (SU(3)) N, P (su(3)) N. Problem for Numerical integration schemes multiplicative structure of Lie group does not fit the additive structure of integrations schemes (RK and others) way-out: bypass via the Lie algebra with its additive structure (Magnus, 1954) M. Günther and M. Wandelt, HMC for Lattice QCD 15/30
17 Bypass via Lie algebra Lie algebra parameterization of U U(t) = Ψ(Ω(t))U 0 with Ψ(t 0 ) = 0 Transformation of ODE ( ) d U = dω Ψ(Ω(t)) Ω(t)U 0 ( )) = (dψ Ω(t) Ω(t) Ψ(Ω(t))U 0 ( )) = (dψ Ω(t) Ω(t) U(t) }{{} =P (t)! ( ) P (t) = dψ Ω(t) Ω(t) Ω(t) = dψ 1 Ω(t) (P (t)). M. Günther and M. Wandelt, HMC for Lattice QCD 16/30
18 Transformed ODE system Solve ODE-IVP on Lie-algebra Ω µ (x) = dψ 1 Ω (P µ(x) µ(x)), Ω µ (x)(t 0 ) = 0, P = x,µ S(Ψ(Ω)U 0 ), P (t 0 ) = P 0 links can be obtained via the transformation U µ (x)(t) = Ψ(Ω µ (x)(t))u 0 Possible choice of parameterization Exponential map: Ψ(Ω) = exp(ω) = 1 k=0 k! Ωk Cayley map: Ψ(Ω) = cay(ω) = (I Ω) 1 (I + Ω) M. Günther and M. Wandelt, HMC for Lattice QCD 17/30
19 Exponential map ψ(ω) = exp(ω) = k=0 1 k! Ωk Ω(t) = d exp 1 Ω (P (t)) = B k k! adk Ω(P (t)) k 0 = P (t) 1 2 [Ω(t), P (t)] + 1 [Ω(t), [Ω(t), P (t)]] Problem: for numerical integration schemes, infinite sum has to be truncated: k 0 q k=0 Compatability with numerical integration schemes: Order p if a) RK scheme has order p and b) q + 2 p holds (Munthe-Kaas). M. Günther and M. Wandelt, HMC for Lattice QCD 18/30
20 Cayley map ψ(ω) = cay(ω) = (I Ω) 1 (I + Ω) Ω(t) = dcay 1 Ω (P (t)) = 1 (I Ω(t))P (t)(i + Ω(t)) 2 No infinite series, no trunaction needed No compatability condition cheap evaluations! M. Günther and M. Wandelt, HMC for Lattice QCD 19/30
21 The basis scheme: leap frog Exponential map: Update (U n, P n ) (U n+1, P n+1 ) P n+ 1 2 = P n h 2 x,µs(u n ), Ω n+1 = Ω n + h d exp 1 Ω (P n+ 1 ) Ω n + h P 2 n+ 1, 2 U n+1 = exp(ω n+1 )U n, P n+1 = P n+ 1 h 2 2 x,µs(u n+1 ). Elimination of the auxiliary Ω variable U n+1 = exp(h (P n h 2 x,µs(u n )))U n, P n+1 = P n h 2 x,µs(u n ) h 2 x,µs(exp(h (P n h 2 x,µs(u n )))) M. Günther and M. Wandelt, HMC for Lattice QCD 20/30
22 The basis scheme: leap frog Cayley map: Update (U n, P n ) (U n+1, P n+1 ) P n+ 1 2 = P n h 2 x,µs(un), Ω n+1 = Ω n + h dcay 1 Ω (P n+ 1 ) 2 = Ω n + h 2 (I Ωn)P n+ 2 1 (I + Ω n), U n+1 = cay(ω n+1 )U n = (I Ω n+1 ) 1 (I + Ω n+1 )U n, P n+1 = P n+ 1 2 h 2 x,µs(u n+1). Elimination of the auxiliary Ω variable U n+1 = cay(h dcay 1 Ω (Pn h 2 x,µs(un)))un, P n+1 = P n h 2 x,µs(un) h 2 x,µs(cay(h dcay 1 Ω (Pn h 2 x,µs(un)))) M. Günther and M. Wandelt, HMC for Lattice QCD 21/30
23 Higher-order schemes Composition Starting point: basis scheme ϕ h t order p, time-reversible & volume-preserving ( ) ( ) Un+1 = ϕ h Un t P n n+1 P n Composition of m basic schemes Volume-preservation ϕ h t n = ϕ γ 1h t n Condition for time-reversibility: ϕ γ 2h t n ϕ γmh t n γ m+1 k = γ k, k = 1,..., m 1 Order p + 1, if the underlying scheme has order p: m m γ j = 1, γ p j = 0 j=1 j=1 For leap-frog to get p = 4: m = 3, γ 1 = γ 3 = (2 3 2) 1/3, γ 2 = 1 2γ 1. M. Günther and M. Wandelt, HMC for Lattice QCD 22/30
24 Higher-order schemes Splitting Starting point: sequence of link and momenta updates Momenta Update U µ (x) = 0 P µ (x) = x,µ S(U) Link Updata U µ (x) = P µ,x U µ,x P µ (x) = 0 in general: order 1 at least order 2, if splitting symmetric important classes: Omelyan and force-gradient schemes M. Günther and M. Wandelt, HMC for Lattice QCD 23/30
25 Higher-order schemes Examples Symmetric five-stage Omelyan scheme P λh U h/2 P (1 2λ)h U h/2 P λh with λ = 1 2 ( ) 1/ ( ) 1/3 Still order 2, but minimal leading error coefficient Five stage force-gradient scheme of order 4 P λh,ξh 3 U h/2 P (1 2λ)h,χh 3 U h/2 P λh,ξh 3 with λ = 1 6, ξ = 0, χ = 1 72 M. Günther and M. Wandelt, HMC for Lattice QCD 24/30
26 Higher-order schemes Discussion Exponential mapping leap-frog scheme used as basic scheme time-reversibility & volume preservation but: truncation order q = 0 Problem: (H(g(U 0, P 0 ) H(U 0, P 0 )) = O(( t) 2 ) only high order only with respect to modified equation of motion Ω µ (x) = P µ (x), U µ,x (t) = exp (Ω µ (x)(t))u µ (x)(0) P µ (x) = x,µ S(U) possible way out: if modified equation of motion based on modified Hamiltonian H(U, P ), use g(u, P ) := H(U, P ) instead of H(U, P ) M. Günther and M. Wandelt, HMC for Lattice QCD 25/30
27 Fermionic fields Additional action: pseudofermionic field Partition function now given by Z = exp ( S g (U))det(DD )du det(dd ) = exp ( S pf (U, φ))dφ, with S pf (U, φ) = φ, (DD ) 1 φ H(U, P ) = 1 2 P, P + S g(u) + S pf (U, φ) M. Günther and M. Wandelt, HMC for Lattice QCD 26/30
28 Computation of pseudofermionic force Additional action: pseudofermionic field x,µ S pf(u,φ) = 2ReT i (DD ) 1 φ, δx,µφ i (δx,µφ)(y) i = 1 δ y,x 2 (1 γ µ)t i U µ (x)d 1 φ(x + ˆµ) +δ y,x+ˆµ 1 2 (1 + γ µ)u 1 µ (x)t i D 1 φ(x) Note: Computation requires two expensive inversions of the Wilson-Dirac operator M. Günther and M. Wandelt, HMC for Lattice QCD 27/30
29 Nested integration - multirate schemes Idea: use different activity levels gauge field S g : fast dynamics, cheap evaluations pseudofermionic field S pf : slow dynamics, expensive evaluations exploit different dynamics by multirate schemes: higher/lower sample rate for fast/slow part introduce additional activity splitting of pseudofermionic field: even-odd preconditioning, domain decomposition and determinant splitting S pf (U, φ ) = M S k [U, φ], m=1 H(U, P ) = P, P /2 + S g (U) + M S k [U, φ]. m=1 M. Günther and M. Wandelt, HMC for Lattice QCD 28/30
30 Sexton and Weingarten First step macro step size h0 micro step size h1 = h 0 /m 1 Leap-frog step ( m1 V H (h 0 ) = V H2 (h 0 /2) V H1 (h 1 )) VH2 (h 0 /2) with M 1 H 1 (U, P ) = P, P /2 + S g[u] + S k [U, φ], H 2 (U, P ) = S M [U, φ] m=1 Nesting next finer step size h2 = h 1 /m 2 further splitting H1 into H 1 (U, P ) = H 11 (U, P ) + H 12 (U, P ) with M 2 H 11 (U, P ) = P, P /2 + S g[u] + S k [U, φ], H 12 (U, P ) = S M 1 [U, φ] m=1 replace VH1 (h 1 ) above by V H1 (h 1 ) = V H12 (h 1 /2) ( V H11 (h 2 )) m2 VH12 (h 1 /2). M. Günther and M. Wandelt, HMC for Lattice QCD 29/30
31 More to read... M. Günther and M. Wandelt, HMC for Lattice QCD 30/30
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 informationDomain Wall Fermion Simulations with the Exact One-Flavor Algorithm
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
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 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 informationAccuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / 11-23 BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal,
More informationAppendix A Notational Conventions
Appendix A Notational Conventions Throughout the book we use Einstein s implicit summation convention: repeated indices in an expression are automatically summed over. We work in natural units where the
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 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 informationNEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING Fernando Casas Universitat Jaume I, Castellón, Spain Fernando.Casas@uji.es (on sabbatical leave at DAMTP, University of Cambridge) Edinburgh,
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 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 informationReal time lattice simulations and heavy quarkonia beyond deconfinement
Real time lattice simulations and heavy quarkonia beyond deconfinement Institute for Theoretical Physics Westfälische Wilhelms-Universität, Münster August 20, 2007 The static potential of the strong interactions
More informationComplex Langevin dynamics for nonabelian gauge theories
Complex Langevin dynamics for nonabelian gauge theories Gert Aarts XQCD 14, June 2014 p. 1 QCD phase diagram QCD partition function Z = DUD ψdψe S YM S F = DU detde S YM at nonzero quark chemical potential
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 informationHamiltonian Monte Carlo with Fewer Momentum Reversals
Hamiltonian Monte Carlo with ewer Momentum Reversals Jascha Sohl-Dickstein December 6, 2 Hamiltonian dynamics with partial momentum refreshment, in the style of Horowitz, Phys. ett. B, 99, explore the
More informationThe symmetries of QCD (and consequences)
The symmetries of QCD (and consequences) Sinéad M. Ryan Trinity College Dublin Quantum Universe Symposium, Groningen, March 2018 Understand nature in terms of fundamental building blocks The Rumsfeld
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 informationAnalytic continuation from an imaginary chemical potential
Analytic continuation from an imaginary chemical potential A numerical study in 2-color QCD (hep-lat/0612018, to appear on JHEP) P. Cea 1,2, L. Cosmai 2, M. D Elia 3 and A. Papa 4 1 Dipartimento di Fisica,
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 informationAnderson átmenet a kvark-gluon plazmában
Anderson átmenet a kvark-gluon plazmában Kovács Tamás György Atomki, Debrecen 25. április 22. Kovács Tamás György Anderson átmenet a kvark-gluon plazmában /9 Együttműködők: Falk Bruckmann (Regensburg U.)
More informationTaylor expansion in chemical potential for 2 flavour QCD with a = 1/4T. Rajiv Gavai and Sourendu Gupta TIFR, Mumbai. April 1, 2004
Taylor expansion in chemical potential for flavour QCD with a = /4T Rajiv Gavai and Sourendu Gupta TIFR, Mumbai April, 004. The conjectured phase diagram, the sign problem and recent solutions. Comparing
More informationAuxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them
Auxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo methods in Fock
More informationDerivation of Electro Weak Unification and Final Form of Standard Model with QCD and Gluons 1 W W W 3
Derivation of Electro Weak Unification and Final Form of Standard Model with QCD and Gluons 1 W 1 + 2 W 2 + 3 W 3 Substitute B = cos W A + sin W Z 0 Sum over first generation particles. up down Left handed
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 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 information6. Auxiliary field continuous time quantum Monte Carlo
6. Auxiliary field continuous time quantum Monte Carlo The purpose of the auxiliary field continuous time quantum Monte Carlo method 1 is to calculate the full Greens function of the Anderson impurity
More information6.1 Quadratic Casimir Invariants
7 Version of May 6, 5 CHAPTER 6. QUANTUM CHROMODYNAMICS Mesons, then are described by a wavefunction and baryons by Φ = q a q a, (6.3) Ψ = ǫ abc q a q b q c. (6.4) This resolves the old paradox that ground
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 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 informationExcited states of the QCD flux tube
Excited states of the QCD flux tube Bastian Brandt Forschungsseminar Quantenfeldtheorie 19.11.2007 Contents 1 Why flux tubes? 2 Flux tubes and Wilson loops 3 Effective stringtheories Nambu-Goto Lüscher-Weisz
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 informationChapter 13. Local Symmetry
Chapter 13 Local Symmetry So far, we have discussed symmetries of the quantum mechanical states. A state is a global (non-local) object describing an amplitude everywhere in space. In relativistic physics,
More informationPhysics 221A Fall 1996 Notes 13 Spins in Magnetic Fields
Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields A nice illustration of rotation operator methods which is also important physically is the problem of spins in magnetic fields. The earliest experiments
More informationA Lattice Study of the Glueball Spectrum
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 288 292 c International Academic Publishers Vol. 35, No. 3, March 15, 2001 A Lattice Study of the Glueball Spectrum LIU Chuan Department of Physics,
More informationLinear Algebra on Lattices: Simit Language Extensions with Applications to Lattice QCD. Gurtej Kanwar
Linear Algebra on Lattices: Simit Language Extensions with Applications to Lattice QCD by Gurtej Kanwar Submitted to the Department of Electrical Engineering and Computer Science in partial fulfillment
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 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 informationarxiv: v1 [hep-lat] 31 Oct 2014
arxiv:4.88v [hep-lat] 3 Oct 24 Zoltán Fodor University of Wuppertal, Department of Physics, Wuppertal D-4297, Germany Jülich Supercomputing Center, Forschungszentrum Jülich, Jülich D-52425, Germany Eötvös
More informationSmall bits of cold, dense matter
Small bits of cold, dense matter Alessandro Roggero (LANL) with: S.Gandolfi & J.Carlson (LANL), J.Lynn (TUD) and S.Reddy (INT) ArXiv:1712.10236 Nuclear ab initio Theories and Neutrino Physics INT - Seattle
More informationFinite-temperature Field Theory
Finite-temperature Field Theory Aleksi Vuorinen CERN Initial Conditions in Heavy Ion Collisions Goa, India, September 2008 Outline Further tools for equilibrium thermodynamics Gauge symmetry Faddeev-Popov
More informationLecture 10. September 28, 2017
Lecture 10 September 28, 2017 The Standard Model s QCD theory Comments on QED calculations Ø The general approach using Feynman diagrams Ø Example of a LO calculation Ø Higher order calculations and running
More informationDimensional reduction near the deconfinement transition
Dimensional reduction near the deconfinement transition Aleksi Kurkela ETH Zürich Wien 27.11.2009 Outline Introduction Dimensional reduction Center symmetry The deconfinement transition: QCD has two remarkable
More informationt Hooft-Polyakov Monopoles on the Lattice
t Hooft-Polyakov Monopoles on the Lattice Davis,Kibble,Rajantie&Shanahan, JHEP11(2000) Davis,Hart,Kibble&Rajantie, PRD65(2002) Rajantie, in progress 19 May 2005 University of Wales, Swansea Introduction
More informationlight-cone (LC) variables
light-cone (LC) variables 4-vector a µ scalar product metric LC basis : transverse metric 24-Apr-13 1 hadron target at rest inclusive DIS target absorbes momentum from γ * ; for example, if q z P z =0
More informationNumerical Stochastic Perturbation Theory and The Gradient Flow
Numerical Stochastic Perturbation Theory and The Gradient Flow Mattia Dalla Brida Trinity College Dublin, Ireland Dirk Hesse Universita degli Studi di Parma, Italia 31st International Symposium on Lattice
More informationRiemann Manifold Methods in Bayesian Statistics
Ricardo Ehlers ehlers@icmc.usp.br Applied Maths and Stats University of São Paulo, Brazil Working Group in Statistical Learning University College Dublin September 2015 Bayesian inference is based on Bayes
More informationarxiv: v1 [hep-lat] 2 May 2012
A CG Method for Multiple Right Hand Sides and Multiple Shifts in Lattice QCD Calculations arxiv:1205.0359v1 [hep-lat] 2 May 2012 Fachbereich C, Mathematik und Naturwissenschaften, Bergische Universität
More informationLecture II: Owe Philipsen. The ideal gas on the lattice. QCD in the static and chiral limit. The strong coupling expansion at finite temperature
Lattice QCD, Hadron Structure and Hadronic Matter Dubna, August/September 2014 Lecture II: Owe Philipsen The ideal gas on the lattice QCD in the static and chiral limit The strong coupling expansion at
More informationDiagrammatic Monte Carlo methods for Fermions
Diagrammatic Monte Carlo methods for Fermions Philipp Werner Department of Physics, Columbia University PRL 97, 7645 (26) PRB 74, 15517 (26) PRB 75, 8518 (27) PRB 76, 235123 (27) PRL 99, 12645 (27) PRL
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 informationClaude Tadonki. MINES ParisTech PSL Research University Centre de Recherche Informatique
Claude Tadonki MINES ParisTech PSL Research University Centre de Recherche Informatique claude.tadonki@mines-paristech.fr Monthly CRI Seminar MINES ParisTech - CRI June 06, 2016, Fontainebleau (France)
More informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationThe Polyakov loop and the Hadron Resonance Gas Model
Issues The Polyakov loop and the Hadron Resonance Gas Model 1, E. Ruiz Arriola 2 and L.L. Salcedo 2 1 Grup de Física Teòrica and IFAE, Departament de Física, Universitat Autònoma de Barcelona, Spain 2
More informationLattice Quantum Chromo Dynamics and the Art of Smearing
Lattice Quantum Chromo Dynamics and the Art of Georg Engel March 25, 2009 KarlFranzensUniversität Graz Advisor: Christian B. Lang Co-workers: M. Limmer and D. Mohler 1 / 29 2 / 29 Continuum Theory Regularization:
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 informationThe Phases of QCD. Thomas Schaefer. North Carolina State University
The Phases of QCD Thomas Schaefer North Carolina State University 1 Plan of the lectures 1. QCD and States of Matter 2. The High Temperature Phase: Theory 3. Exploring QCD at High Temperature: Experiment
More informationParticle Physics WS 2012/13 ( )
Particle Physics WS 01/13 (3.11.01) Stephanie Hansmann-Menzemer Physikalisches Institut, INF 6, 3.101 Content of Today Structure of the proton: Inelastic proton scattering can be described by elastic scattering
More informationProgress in Gauge-Higgs Unification on the Lattice
Progress in Gauge-Higgs Unification on the Lattice Kyoko Yoneyama (Wuppertal University) in collaboration with Francesco Knechtli(Wuppertal University) ikos Irges(ational Technical University of Athens)
More informationRandom Walks on Hyperbolic Groups IV
Random Walks on Hyperbolic Groups IV Steve Lalley University of Chicago January 2014 Automatic Structure Thermodynamic Formalism Local Limit Theorem Final Challenge: Excursions I. Automatic Structure of
More informationQCD Phases with Functional Methods
QCD Phases with Mario PhD-Advisors: Bernd-Jochen Schaefer Reinhard Alkofer Karl-Franzens-Universität Graz Institut für Physik Fachbereich Theoretische Physik Rab, September 2010 QCD Phases with Table of
More informationLorentzian elasticity arxiv:
Lorentzian elasticity arxiv:1805.01303 Matteo Capoferri and Dmitri Vassiliev University College London 14 July 2018 Abstract formulation of elasticity theory Consider a manifold M equipped with non-degenerate
More informationLattice QCD study for relation between quark-confinement and chiral symmetry breaking
Lattice QCD study for relation between quark-confinement and chiral symmetry breaking Quantum Hadron Physics Laboratory, Nishina Center, RIKEN Takahiro M. Doi ( 土居孝寛 ) In collaboration with Hideo Suganuma
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 informationA short overview on strong interaction and quantum chromodynamics
A short overview on strong interaction and quantum chromodynamics Christoph Klein Universität Siegen Doktorandenseminar 08.10.2008 Christoph Klein (Universität Siegen) QCD and strong interaction Doktorandenseminar
More informationBRST and Dirac Cohomology
BRST and Dirac Cohomology Peter Woit Columbia University Dartmouth Math Dept., October 23, 2008 Peter Woit (Columbia University) BRST and Dirac Cohomology October 2008 1 / 23 Outline 1 Introduction 2 Representation
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 informationOn the definition and interpretation of a static quark anti-quark potential in the colour-adjoint channel
On the definition and interpretation of a static quark anti-quark potential in the colour-adjoint channel Effective Field Theory Seminar Technische Universität München, Germany Marc Wagner, Owe Philipsen
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 informationScreening mass of gluons in presence of external Abelian chromomagnetic field
Screening mass of gluons in presence of external Abelian chromomagnetic field N. V. Kolomoyets, V. V. Skalozub Dnepropetrovsk National University Ukraine December 5, 2018 N. V. Kolomoyets Screening mass
More informationLecture 9 Valence Quark Model of Hadrons
Lecture 9 Valence Quark Model of Hadrons Isospin symmetry SU(3) flavour symmetry Meson & Baryon states Hadronic wavefunctions Masses and magnetic moments Heavy quark states 1 Isospin Symmetry Strong interactions
More informationHeisenberg-Euler effective lagrangians
Heisenberg-Euler effective lagrangians Appunti per il corso di Fisica eorica 7/8 3.5.8 Fiorenzo Bastianelli We derive here effective lagrangians for the electromagnetic field induced by a loop of charged
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationExtraction of Quark Distributions on Transverse Spin of the Nucleon at the HERMES Experiment. Hidekazu Tanaka
Extraction of Quark Distributions on Transverse Spin of the Nucleon at the HERMES Experiment A Dissertation By Hidekazu Tanaka February 25 Department of Physics Tokyo Institute of Technology Abstract
More informationAdaptive Multigrid for QCD. Lattice University of Regensburg
Lattice 2007 University of Regensburg Michael Clark Boston University with J. Brannick, R. Brower, J. Osborn and C. Rebbi -1- Lattice 2007, University of Regensburg Talk Outline Introduction to Multigrid
More informationGinsparg-Wilson Fermions and the Chiral Gross-Neveu Model
Ginsparg-Wilson Fermions and the DESY Zeuthen 14th September 2004 Ginsparg-Wilson Fermions and the QCD predictions Perturbative QCD only applicable at high energy ( 1 GeV) At low energies (100 MeV - 1
More informationCanonical partition functions in lattice QCD
Canonical partition functions in lattice QCD (an appetizer) christof.gattringer@uni-graz.at Julia Danzer, Christof Gattringer, Ludovit Liptak, Marina Marinkovic Canonical partition functions Grand canonical
More informationMeron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations
Meron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations Uwe-Jens Wiese Bern University IPAM Workshop QS2009, January 26, 2009 Collaborators: B. B. Beard
More informationNew approaches to strongly interacting Fermi gases
New approaches to strongly interacting Fermi gases Joaquín E. Drut The Ohio State University INT Program Simulations and Symmetries Seattle, March 2010 In collaboration with Timo A. Lähde Aalto University,
More informationA tömeg eredete a látható világegyetemben
A tömeg eredete a látható világegyetemben ELTE Elméleti Fizikai Tanszék MTA-ELTE Lendület rácstérelmélet kutatócsoport (Borsányi Szabolcs, Stephan Dürr,Fodor Zoltán, Christian Hoelbling, Stefan Krieg,
More informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ Bayesian paradigm Consistent use of probability theory
More informationProperties of ground and excited state hadrons from lattice QCD
Properties of ground and excited state hadrons from lattice QCD Daniel Mohler TRIUMF, Theory Group Vancouver, B.C., Canada Newport News, March 22 200 Co-Workers: Christof Gattringer, Christian Lang, Georg
More informationx 3 x 1 ix 2 x 1 + ix 2 x 3
Peeter Joot peeterjoot@pm.me PHY2403H Quantum Field Theory. Lecture 19: Pauli matrices, Weyl spinors, SL2,c, Weyl action, Weyl equation, Dirac matrix, Dirac action, Dirac Lagrangian. Taught by Prof. Erich
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem
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 informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ 1 Bayesian paradigm Consistent use of probability theory
More informationRevisiting the strong coupling limit of lattice QCD
Revisiting the strong coupling limit of lattice QCD Philippe de Forcrand ETH Zürich and CERN with Michael Fromm (ETH) Motivation Intro Algorithm Results Concl. 25 + years of analytic predictions: 80 s:
More informationQuantum simulation with string-bond states: Joining PEPS and Monte Carlo
Quantum simulation with string-bond states: Joining PEPS and Monte Carlo N. Schuch 1, A. Sfondrini 1,2, F. Mezzacapo 1, J. Cerrillo 1,3, M. Wolf 1,4, F. Verstraete 5, I. Cirac 1 1 Max-Planck-Institute
More informationImprovement in Stochastic Noise of Disconnected Diagrams
University of Cyprus Department of Physics Bachelor Thesis Improvement in Stochastic Noise of Disconnected Diagrams Oikonomou Antreas May 207 University of Cyprus Department of Physics Bachelor Thesis
More informationQuantum Field Theory II
Quantum Field Theory II Christof Wetterich April 5, 006 This is a first draft of the script to Prof. Wetterich s QFT II lecture at the university of Heidelberg in the summer semester 006. The latest version
More informationQuantum Field Theory. and the Standard Model. !H Cambridge UNIVERSITY PRESS MATTHEW D. SCHWARTZ. Harvard University
Quantum Field Theory and the Standard Model MATTHEW D. Harvard University SCHWARTZ!H Cambridge UNIVERSITY PRESS t Contents v Preface page xv Part I Field theory 1 1 Microscopic theory of radiation 3 1.1
More informationSU(3) Quantum Spin Ladders as a Regularization of the CP (2) Model at Non-Zero Density: From Classical to Quantum Simulation
SU(3) Quantum Spin Ladders as a Regularization of the CP (2) Model at Non-Zero Density: From Classical to Quantum Simulation W. Evans, U. Gerber, M. Hornung, and U.-J. Wiese arxiv:183.4767v1 [hep-lat]
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationOn local normal forms of completely integrable systems
On local normal forms of completely integrable systems ENCUENTRO DE Răzvan M. Tudoran West University of Timişoara, România Supported by CNCS-UEFISCDI project PN-II-RU-TE-2011-3-0103 GEOMETRÍA DIFERENCIAL,
More informationSpectral Action from Anomalies
Spectral Action from Anomalies Fedele Lizzi Università di Napoli Federico II and Insitut de Ciencies del Cosmo, Universitat de Barcelona Work in collaboration with A.A. Andrianov & M. Kurkov (St. Petersburg)
More informationThe 3-d 3-state Potts model with a non-zero density of static color charges
The 3-d 3-state Potts model with a non-zero density of static color charges Masterarbeit der Philosophisch-naturwissenschaftlichen Fakultät der Universität Bern vorgelegt von Philippe Widmer 2011 Betreut
More informationExcited States of the Nucleon in Lattice QCD
Collaborators: Alan Ó Cais, Waseem Kamleh, Ben G. Lasscock, Derek B. Leinweber and Anthony G. Williams CSSM, University of Adelaide Adelaide Outline Introduction 1 Introduction 2 3 4 2 point Correlation
More informationThe groups SO(3) and SU(2) and their representations
CHAPTER VI The groups SO(3) and SU() and their representations Two continuous groups of transformations that play an important role in physics are the special orthogonal group of order 3, SO(3), and the
More informationSection 11: Review. µ1 x < 0
Physics 14a: Quantum Mechanics I Section 11: Review Spring 015, Harvard Below are some sample problems to help study for the final. The practice final handed out is a better estimate for the actual length
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 informationarxiv:hep-th/ v2 31 Jul 2000
Hopf term induced by fermions Alexander G. Abanov 12-105, Department of Physics, MIT, 77 Massachusetts Ave., Cambridge, MA 02139, U.S.A. arxiv:hep-th/0005150v2 31 Jul 2000 Abstract We derive an effective
More information