Polynomial Filtered Hybrid Monte Carlo
|
|
- Christal Ray
- 5 years ago
- Views:
Transcription
1 Polynomial Filtered Hybrid Monte Carlo Waseem Kamleh and Michael J. Peardon CSSM & U NIVERSITY OF A DELAIDE QCDNA VII, July 4th 6th, 2012
2 Introduction Generating large, light quark dynamical gauge field configurations is expensive. Hybrid Monte Carlo (HMC) is still the most used algorithm for generating dynamical configurations. There have been many improvements to the basic HMC algorithm developed over the years. A partial list of multi-scale type algorithmic improvements: Domain Decomposition method (Luscher). Mass Preconditioning (Hasenbuch). Polynomial Filtering (Peardon & Sexton, WK & Peardon). We begin by quickly reviewing the basics of the HMC algorithm.
3 Hybrid Monte Carlo Review The lattice is embedded in a Hamiltonian system by the addition of a fictitious simulation time τ, along with an additional fictitious field P which are the conjugate momenta to U, H(P, U) = x,µ 1 2 Tr P µ(x) 2 + S[U]. The conjugate momenta P µ (x) are drawn from a Gaussian distribution. In this way, the Hamiltonian H is constructed so that after path integration the expectation values of observables are unaltered. Where are the fermions?
4 Fermion Determinant The fermion fields ψ and ψ are Grassmannian, hence to perform simulations we need to integrate them out. We have that det D w = D ψdψ e d 4 x ψ(x)d w ψ(x). Define the effective action as S eff [U] = S g [U] ln det D w [U]. We can write the action for full QCD in terms of bosonic fields φ using the identity det M = 1 det M 1 = Dφ Dφ e d 4 x φ (x)m 1 φ(x).
5 Pseudofermions The convergence of the Gaussian integral in φ is only guaranteed for Hermitian positive definite (Hpd) M. For Wilson-type fermions, D w is a complex matrix, but det D w is real and positive. So, we can define M = D wd w, and M will be Hpd with det M = det D 2 f. Then the pseudofermionic effective action for full QCD with two flavours of degenerate quarks is S eff [U] = S g [U] + d 4 x φ (x)m 1 [U]φ(x). The pseudofermion fields φ(x) are drawn from a Gaussian distribution.
6 Hybrid Monte Carlo The Hybrid Monte Carlo algorithm creates a Markov chain by alternately performing two steps: A Molecular Dynamics (MD) integration to generate a new configuration (U U, P P ). A Metropolis accept/reject step on the proposed configuration (U, P ). The accept/reject step is based upon the change in the Hamiltonian ρ(u U, P P ) e H. The Molecular Dynamics integration takes place along a trajectory which for sufficiently small integration step sizes h approximately conserves H, hence yielding high acceptance rates.
7 Molecular Dynamics Start by requiring that the Hamiltonian be conserved along the trajectory dh dτ = 0. Then derive the discretised equations of motion, U µ (x, τ + h) = U µ (x, τ) exp ( ihp µ (x, τ) ), δs P µ (x, τ + h) = P µ (x, τ) U µ (x, τ) δu µ (x, τ). The derivative of the action with respect to the gauge fields is known as the force term, F µ (x) = δs δu µ (x).
8 Splitting The Action Lets split our action into its gauge field and pseudofermion field components, S = S g + S pf, where 1 S g = β x,µ<ν 3 Re Tr(1 U µν(x)), S pf = φ (x)(d wd w ) 1 φ(x), x and D w is the Wilson fermion matrix. Each of the terms in the action induces a force term. The size (variance?) of the force term is the dominant factor in determining what step size h is need for a given acceptance probability ρ acc. As the quark mass becomes lighter, the size of the pseudofermion force term increases.
9 Leapfrog Integration The MD equations of motion induce corresponding time evolution operators, V T (h) : {U(τ), P(τ)} {U(τ + h), P(τ)}, V S (h) : {U(τ), P(τ)} {U(τ), P(τ + h)}. The simplest MD integration scheme is the leapfrog V(h) = V S ( h 2 )V T(h)V S ( h 2 ). MD Integration trajectories typically have unit length, and hence as the step size h decreases, the number of integration steps increases.
10 Pseudofermion Force Each time we act with V S (h) we need to evaluate the pseudofermion force term, F pf = δs pf δu. This involves inverting the fermion matrix, and hence is expensive! However, for split actions S = S 1 + S 2 we can use a multiple time scale integration scheme (nested leapfrog), V(h) = V 2 ( h 2 ) [V 1 ( h m ) ] m V 2 ( h 2 ), V 1 (h) = V S1 ( h 2 )V T(h)V S1 ( h 2 ), V 2 = V S2 (h).
11 Multiple Time Scales Multiple time scale integration is effective when the force term F 1 due to S 1 is cheap to evaluate compared to F 2 (that of S 2 ). However, as the step-size for S 2 is larger, we also require that the size of the force term for S 2 is relatively small compared to that of S 1. The gauge force F g is cheap compared to the pseudofermion force F pf, and at heavy quark masses F g > F pf, but at light quark masses the UV fluctuations in the pseudo fermion force become too large for multiple time scales to be effective. This is where polynomial filtering steps in.
12 Polynomial Filtering We can use a polynomial filter P = P(M) to separate the ultraviolet and infrared physics in the pseudofermion force, S poly = χ Pχ, S pf = φ (MP) 1 φ. Recall M = D wd w is Hermitian positive definite. As S poly is fast to evaluate compared to S pf we split the action in the following way, S 1 = S g + S poly, S 2 = S pf. If P 1/z then it will capture the short-distance physics in S poly and act as a UV filter in S pf.
13 Fermionic Determinant Note that the fermionic determinant is unchanged by the introduction of the polynomial filter. To see this we note that Dφ DφDχ Dχ e d 4 x χ Pχ+φ (x)(mp) 1 φ(x) = det P det(mp) = 1 = det M det M 1 In the limit as the order of the polynomial n we have MP 1 and P(M) M 1. In this way the polynomial term can reproduce as much or as little of the pseudofermion term as we want.
14 Chebyshev Filter An effective UV filter is the n th order Hermitian Chebyshev polynomial approximation to 1/z, P n (z) = a n n k=1 (z z k ) 1 z, where we set θ k = 2πk n+1 to obtain the roots z k = λ[ 1 2 (1 + ɛ)(1 cos θ k) i ɛ sin θ k ]. The normalisation is defined by z 0 = 1 2 (1 + ɛ), with a n = 1 z 0 n k=1 (z 0 z k ). The approximation is good between [ɛ, 1], so we rescale with λ = 1 + 8κ.
15 Chebyshev Roots n = 24
16 Intermediate Filter Using this property of the Chebyshev polynomials, we can add an intermediate filter to our action as a bridge between the high and low energy scales, e.g. S 1 = S g S 2 = χ P m χ, S 3 = χ P m n χ, S 4 = φ (MP n ) 1 φ. P n denotes the Chebyshev polynomial of order n. P m n is defined so that P n = P m P m n. This allows us to perform fermion matrix inversions even less frequently. This may(?) be more efficient than a single filter algorithm. Note: filter implementation makes use of multi-shift linear solvers for efficiency.
17 Integration Scheme With a second filter, we have four different scales, five (or more) if we are doing a 2+1 flavour simulation. Nested leapfrog is too cumbersome for fine tuning these scales. We make use of the fact that V i = V Si all commute to introduce a generalised integration scheme. If N i is the number of integration steps per trajectory, then nested leapfrog requires N i N i 1 i The new integration scheme requires only that N i N 1 i
18 Generalised Leapfrog Set N = N 1 and set n i = N 1 /N i The generalised leapfrog algorithm is then: 1 Perform an initial half-step V i ( 1 2 h i) updating P for all i. 2 Loop over j = 1 to N 1 Apply V T (h) to update U. If {0 j mod N i } apply V i (h i ) to update P 3 Apply V T (h) to update U. 4 Perform a final half-step V i ( 1 2 h i) updating P for all i. Can show using BCH that it has errors of O[(h) 2 ]. Reduces to nested leapfrog for N i N i 1.
19 Integrator Error Analysis Given a Hamiltonian H the evolution operator with stepsize h for our system is exp hĥ. Ĥ is the linear operator on the vector space of functions f on phase space (p, q) defined by the Poisson bracket Ĥf = {H, f } = i Write the Hamiltonian as ( H f H ) f. p i q i q i p i H = T + S 1 + S 2 + S 3 + S
20 Integrator Error Analysis For each term in the Hamiltonian we can define a linear operator using the Poisson bracket relation above. We make use of the Baker-Campbell-Hausdorff result, e λâ e λ ˆB e λâ = exp (λ(2â + ˆB) + λ3 ) 6 ([[Â, ˆB], Â] + [[Â, ˆB], ˆB]) + O(λ 4 ) Apply this to the generalised leapfrog integrator: Use the simplest non-trivial case, H = T + S 1 + S 2. Number of integration steps N T = 6, N 1 = 3 and N 2 = 2.
21 Integrator Error Analysis The integrator for this simplest non-trivial case is V(h) = e h 4 Ŝ2 e h 6 Ŝ1 e h 3 ˆT e h 3 Ŝ1 e h 6 ˆT e h 2 Ŝ2 e h 6 ˆT e h 3 Ŝ1 e h 3 ˆT e h 6 Ŝ1 e h 4 Ŝ2. Repeatedly apply BCH result to obtain ( V(h) = exp hĥ + h 3( 1 48 [[Ŝ 2, ˆT], ˆT] [[Ŝ 2, ˆT], Ŝ 2 ] [[Ŝ 1, ˆT], Ŝ 1 ] [[Ŝ 1, ˆT], ˆT] [[Ŝ 1, ˆT], Ŝ 2 ] )) The error in the generalised integrator relative to the leading term is O(h 2 ), just as for the regular leapfrog.
22 Integrator Error Analysis Examine the individual leapfrog integrators H 1 = T + S 1, V 1 (h) = e h 6 Ŝ1 e h 3 ˆT e h 3 Ŝ1 e h 3 ˆT e h 3 Ŝ1 e h 3 ˆT e h 6 Ŝ1 H 2 = T + S 2, V 2 (h) = e h 4 Ŝ2 e h 2 ˆT e h 2 Ŝ2 e h 2 ˆT e h 4 Ŝ2. Using BCH we obtain ( V 1 (h) = exp hĥ 1 + h 3 ( 1 V 2 (h) = exp ) 216 [[Ŝ 1, ˆT], Ŝ 1 ]), ) 108 [[Ŝ 1, ˆT], ˆT] + 1 ( hĥ 2 + h 3 ( 1 48 [[Ŝ 2, ˆT], ˆT] [[Ŝ 2, ˆT], Ŝ 2 ]) The only difference is the cross term [[Ŝ 1, ˆT], Ŝ 2 ]..
23 Pseudofermion Force, Single Filter
24 Polynomial Force, Single Filter
25 ρ acc = erfc(h 2 F /τ2 0 ), Single Filter
26 Characteristic Scale, Single Filter
27 Pseudofermion Force, Double Filter
28 Polynomial Force, Double Filter
29 ρ acc = erfc(h 2 F /τ2 0 ), Double Filter
30 Characteristic Scale, Double Filter
31 Results κ = , m π = 665MeV n p n q 1/h g 1/h p 1/h q 1/h cg
32 Results κ = , m π = 665MeV
33 Pseudofermion Force, κ =
34 Polynomial Force, κ =
35 ρ acc = erfc(h 2 F /τ2 0 ), κ =
36 Characteristic Scale, Double Filter
37 Results κ = , m π 400MeV n p n q 1/h g 1/h p 1/h q 1/h cg
38 Results κ = , m π 400MeV
39
40 Single Flavour QCD Single fermion flavours can be simulated using a rational polynomial, R(M) = a i 1 M + b i M (or some other method e.g. polynomial approx. to 1/ M). Recall M = D wd w is Hermitian positive definite. Can we extend our polynomial filtering technique to single flavour simulations? Suppose we have a polynomial Q such that Q(M) 1 M E.g. numerically calculate coefficients for Chebyshev approximation to 1/ z, then calculate the roots...
41 Single Flavour Polynomial Filter Then we can write S poly1f = χ 1f Qχ 1f, S 1pf = φ 1f RQ 1 φ 1f. As before, the determinant is unaffected by the addition of the polynomial filter. Could add an intermediate filter for the single flavour as before but at the strange quark mass probably not worth it. Knowing roots of Q is needed to rewrite RQ 1 as a sum over poles. This allows the use of efficient linear multi-shift system solvers.
42 Alternative Approximation Note that RQ 1 1. May be advantageous to simply use the Remes algorithm to get a rational approximation to R(z) f (z) = 1 zq(z) and use that instead of the product RQ 1. What are the possible advantages? The rational approximation to f might have improved precision for a given order. Smallest shift for R expressed as a sum over poles might be bigger. = Less iterations to solve.
43 Another Variant Recall that we used the Hermitian Chebyshev approximation to 1/z in our two flavour polynomial filter. The non-hermitian Chebyshev approximation K(z) to 1/z has the same normalisation, but slightly different roots, y k = d(1 cos θ k ) + i d 2 c 2 sin θ k. Valid for an elliptical region in the complex plane. So long as the spectrum of the non-hermitian matrix D w is within this ellipse, we can write n/2 K(D w ) = a n (D w y k )(D w yk ) 1 k=1 D w
44 Factoring the Polynomial Now construct a polynomial using only half the roots (say those with positive imaginary parts), K + (D w ) = n/2 a n (D w y k ) k=1 We need the following two properties of determinants, det(ab) = det A det B det A = (det A) Using these we can deduce that det K +(D w )K + (D w ) = det K(D w ) (det D w ) 1 So long as det D w is real and positive, this is the correct weighting for a single fermion flavour.
45 The Action So we can the construct a polynomial filtered one-flavour action using K +, S poly1f = χ 1f K +K + χ 1f, S 1pf = φ 1f W (D w )W(D w )φ 1f. Here, W needs to be a rational polynomial approximation to {zk +(z)k + (z)} 1. Should be able to obtain this by factoring R(z) the Zolotarev approximation to 1/z and setting W(z) = R + (z)k 1 +. At this point, it might be more efficient to take large n polynomial limit for K +, split it into two and do filtered polynomial HMC for the single flavour...
46 Summary The use of a polynomial approximation to the inverse as a filter successfully separates the UV and IR pseudofermion dynamics. We successfully reduce the cost of dynamical simulations. Polynomial filtering is very simple. Minimal modification to existing code. Technique is not necessarily orthogonal to other improvements (e.g. DD).
47 Summary The generalised leapfrog algorithm is applicable to any multiple time scale integration scheme, far more flexible than nested leapfrog. Polynomial filtering can also be extended to single flavour simulations. Light quark mass results to come, plus comparison.
Domain 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 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 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 informationPhases and facets of 2-colour matter
Phases and facets of 2-colour matter Jon-Ivar Skullerud with Tamer Boz, Seamus Cotter, Leonard Fister Pietro Giudice, Simon Hands Maynooth University New Directions in Subatomic Physics, CSSM, 10 March
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 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 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 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 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 informationQCD on the lattice - an introduction
QCD on the lattice - an introduction Mike Peardon School of Mathematics, Trinity College Dublin Currently on sabbatical leave at JLab HUGS 2008 - Jefferson Lab, June 3, 2008 Mike Peardon (TCD) QCD on the
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 informationMonte Carlo studies of the spontaneous rotational symmetry breaking in dimensionally reduced super Yang-Mills models
Monte Carlo studies of the spontaneous rotational symmetry breaking in dimensionally reduced super Yang-Mills models K. N. Anagnostopoulos 1 T. Azuma 2 J. Nishimura 3 1 Department of Physics, National
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 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] 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 informationNucleon excited states on the lattice
Nucleon excited states on the lattice C.B. Lang, Valentina Verduci Graz, 12.12.12 Overview Motivation: Investigate the nucleon spectrum studying the π N scattering. The approach: How do we extract excited
More informationQFT PS1: Bosonic Annihilation and Creation Operators (11/10/17) 1
QFT PS1: Bosonic Annihilation and Creation Operators (11/10/17) 1 Problem Sheet 1: Bosonic Annihilation and Creation Operators Comments on these questions are always welcome. For instance if you spot any
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 informationAutocorrelation studies in two-flavour Wilson Lattice QCD using DD-HMC algorithm
Autocorrelation studies in two-flavour Wilson Lattice QCD using DD-HMC algorithm Abhishek Chowdhury, Asit K. De, Sangita De Sarkar, A. Harindranath, Jyotirmoy Maiti, Santanu Mondal, Anwesa Sarkar June
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 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 informationFinite Temperature Field Theory
Finite Temperature Field Theory Dietrich Bödeker, Universität Bielefeld 1. Thermodynamics (better: thermo-statics) (a) Imaginary time formalism (b) free energy: scalar particles, resummation i. pedestrian
More informationA Noisy Monte Carlo Algorithm
A Noisy Monte Carlo Algorithm UK/99 06 May 1999 hep-lat/9905033 arxiv:hep-lat/9905033v2 12 May 2000 L. Lin, K.F. Liu, and J. Sloan Dept. of Physics and Astronomy, Univ. of Kentucky, Lexington, KY 40506
More informationIntroduction to lattice QCD simulations
Introduction to lattice QCD simulations Hideo Matsufuru High Energy Accelerator Research Organization (KEK), Oho 1-1, Tsukuba 305-0801, Japan Email: hideo.matsufuru@kek.jp WWW: http://www.rcnp.osaka-u.ac.jp/
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 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 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 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 informationThe interplay of flavour- and Polyakov-loop- degrees of freedom
The interplay of flavour- and Polyakov-loopdegrees of freedom A PNJL model analysis Simon Rößner, Nino Bratović, Thomas Hell and Wolfram Weise Physik Department Technische Universität München Thursday,
More informationFinite Chemical Potential in N t = 6 QCD
Finite Chemical Potential in N t = 6 QCD Rajiv Gavai and Sourendu Gupta ILGTI: TIFR Lattice 2008, Williamsburg July 15, 2008 Rajiv Gavai and Sourendu Gupta ILGTI: TIFRLattice Finite Chemical 2008, Williamsburg
More informationComplex Saddle Points in Finite Density QCD
Complex Saddle Points in Finite Density QCD Michael C. Ogilvie Washington University in St. Louis in collaboration with Hiromichi Nishimura (Bielefeld) and Kamal Pangeni (WUSTL) XQCD4 June 9th, 24 Outline
More informationEffective theories for QCD at finite temperature and density from strong coupling
XQCD 2011 San Carlos, July 2011 Effective theories for QCD at finite temperature and density from strong coupling Owe Philipsen Introduction to strong coupling expansions SCE for finite temperature: free
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 informationThe critical end point of QCD: lattice and experiment
The critical end point of QCD: lattice and experiment Sourendu Gupta ILGTI: TIFR Patnitop 2009 January 2, 2010 SG (ILGTI: TIFR) CEP: lattice and experiment Patnitop 09 1 / 28 Outline 1 On lattice 2 In
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 informationQCD Vacuum, Centre Vortices and Flux Tubes
QCD Vacuum, Centre Vortices and Flux Tubes Derek Leinweber Centre for the Subatomic Structure of Matter and Department of Physics University of Adelaide QCD Vacuum, Centre Vortices and Flux Tubes p.1/50
More informationIs the up-quark massless? Hartmut Wittig DESY
Is the up-quark massless? Hartmut Wittig DESY Wuppertal, 5 November 2001 Quark mass ratios in Chiral Perturbation Theory Leutwyler s ellipse: ( mu m d ) 2 + 1 Q 2 ( ms m d ) 2 = 1 25 m s m d 38 R 44 0
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 informationMonte Carlo simulation inspired by computational optimization. Colin Fox Al Parker, John Bardsley MCQMC Feb 2012, Sydney
Monte Carlo simulation inspired by computational optimization Colin Fox fox@physics.otago.ac.nz Al Parker, John Bardsley MCQMC Feb 2012, Sydney Sampling from π(x) Consider : x is high-dimensional (10 4
More informationLattice QCD. QCD 2002, I. I. T. Kanpur, November 19, 2002 R. V. Gavai Top 1
Lattice QCD QCD 2002, I. I. T. Kanpur, November 19, 2002 R. V. Gavai Top 1 Lattice QCD : Some Topics QCD 2002, I. I. T. Kanpur, November 19, 2002 R. V. Gavai Top 1 Lattice QCD : Some Topics Basic Lattice
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
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 informationarxiv:hep-lat/ v1 6 Oct 2000
1 Scalar and Tensor Glueballs on Asymmetric Coarse Lattices C. Liu a, a Department of Physics, Peking University, Beijing 100871, P. R. China arxiv:hep-lat/0010007v1 6 Oct 2000 Scalar and tensor glueball
More informationNumerical methods for lattice field theory
Numerical methods for lattice field theory Mike Peardon Trinity College Dublin August 9, 2007 Mike Peardon (Trinity College Dublin) Numerical methods for lattice field theory August 9, 2007 1 / 24 Numerical
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 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 information2. Formulation of fermion theory, doubling phenomenon. Euclideanize, introduces 4d cubic lattice. On links introduce (for QCD) SU(3) matrices U n1,n
Chapter 11 Lattice Gauge As I have mentioned repeatedly, this is the ultimate definition of QCD. (For electroweak theory, there is no satisfactory non-perturbative definition). I also discussed before
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 informationUnitary Fermi Gas: Quarky Methods
Unitary Fermi Gas: Quarky Methods Matthew Wingate DAMTP, U. of Cambridge Outline Fermion Lagrangian Monte Carlo calculation of Tc Superfluid EFT Random matrix theory Fermion L Dilute Fermi gas, 2 spins
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationWave functions of the Nucleon
Wave functions of the Nucleon Samuel D. Thomas (1) Collaborators: Waseem Kamleh (1), Derek B. Leinweber (1), Dale S. Roberts (1,2) (1) CSSM, University of Adelaide, (2) Australian National University LHPV,
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 informationGeometry and Physics. Amer Iqbal. March 4, 2010
March 4, 2010 Many uses of Mathematics in Physics The language of the physical world is mathematics. Quantitative understanding of the world around us requires the precise language of mathematics. Symmetries
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 informationThe Removal of Critical Slowing Down. Lattice College of William and Mary
The Removal of Critical Slowing Down Lattice 2008 College of William and Mary Michael Clark Boston University James Brannick, Rich Brower, Tom Manteuffel, Steve McCormick, James Osborn, Claudio Rebbi 1
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationLefschetz-thimble path integral and its physical application
Lefschetz-thimble path integral and its physical application Yuya Tanizaki Department of Physics, The University of Tokyo Theoretical Research Division, Nishina Center, RIKEN May 21, 2015 @ KEK Theory
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 informationLectures on Standard Model Effective Field Theory
Lectures on Standard Model Effective Field Theory Yi Liao Nankai Univ SYS Univ, July 24-28, 2017 Page 1 Outline 1 Lecture 4: Standard Model EFT: Dimension-six Operators SYS Univ, July 24-28, 2017 Page
More information2. The interaction between an electron gas and the potential field is given by
68A. Solutions to Exercises II. April 5. Carry out a Hubbard Stratonovich transformation for the following interaction Hamiltonians, factorizing the terms in brackets, and write one or two sentences interpreting
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationcondensates and topology fixing action
condensates and topology fixing action Hidenori Fukaya YITP, Kyoto Univ. hep-lat/0403024 Collaboration with T.Onogi (YITP) 1. Introduction Why topology fixing action? An action proposed by Luscher provide
More informationHamiltonian Monte Carlo
Chapter 7 Hamiltonian Monte Carlo As with the Metropolis Hastings algorithm, Hamiltonian (or hybrid) Monte Carlo (HMC) is an idea that has been knocking around in the physics literature since the 1980s
More informationDetermination of running coupling αs for N f = QCD with Schrödinger functional scheme. Yusuke Taniguchi for PACS-CS collaboration
Determination of running coupling αs for N f = 2 + QCD with Schrödinger functional scheme Yusuke Taniguchi for PACS-CS collaboration Purpose of this project Determine the fundamental parameter of N f =
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationSymmetries, Fields and Particles. Examples 1.
Symmetries, Fields and Particles. Examples 1. 1. O(n) consists of n n real matrices M satisfying M T M = I. Check that O(n) is a group. U(n) consists of n n complex matrices U satisfying U U = I. Check
More informationNucleon Spectroscopy with Multi-Particle Operators
Nucleon Spectroscopy with Multi-Particle Operators Adrian L. Kiratidis Waseem Kamleh, Derek B. Leinweber CSSM, The University of Adelaide 21st of July - LHP V 2015 - Cairns, Australia 1/41 Table of content
More information2.1 Calculation of the ground state energy via path integral
Chapter 2 Instantons in Quantum Mechanics Before describing the instantons and their effects in non-abelian gauge theories, it is instructive to become familiar with the relevant ideas in the context of
More information7.1 Coupling from the Past
Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov
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 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 informationIntroduction to Integrability
Introduction to Integrability Problem Sets ETH Zurich, HS16 Prof. N. Beisert, A. Garus c 2016 Niklas Beisert, ETH Zurich This document as well as its parts is protected by copyright. This work is licensed
More informationMonte Carlo Methods in Quantum Field Theory
DESY 07-076 Monte Carlo Methods in Quantum Field Theory I. Montvay arxiv:0705.4356v1 [hep-lat] 30 May 2007 Deutsches Elektronen-Synchrotron DESY Abstract In these lecture notes some applications of Monte
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 informationQFT 3 : Problem Set 1
QFT 3 : Problem Set.) Peskin & Schroeder 5. (a.) The basis for the fundamental representation of SU(N) is formed by N N traceless Hermitian matrices. The number of such matrices is = N. For SU(3) this
More informationIntroduction to Supersymmetry
Introduction to Supersymmetry I. Antoniadis Albert Einstein Center - ITP Lecture 5 Grand Unification I. Antoniadis (Supersymmetry) 1 / 22 Grand Unification Standard Model: remnant of a larger gauge symmetry
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 informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.4: Krylov Subspace Methods 2 / 51 Krylov Subspace Methods In this chapter we give
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 informationSimple Evaluation of the Chiral Jacobian with the Overlap Dirac Operator
141 Progress of Theoretical Physics, Vol. 102, No. 1, July 1999 Simple Evaluation of the Chiral Jacobian with the Overlap Dirac Operator Hiroshi Suzuki ) Department of Physics, Ibaraki University, Mito
More informationTime Evolving Block Decimation Algorithm
Time Evolving Block Decimation Algorithm Application to bosons on a lattice Jakub Zakrzewski Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, Jagiellonian University,
More informationNon-perturbative Study of Chiral Phase Transition
Non-perturbative Study of Chiral Phase Transition Ana Juričić Advisor: Bernd-Jochen Schaefer University of Graz Graz, January 9, 2013 Table of Contents Chiral Phase Transition in Low Energy QCD Renormalization
More informationA Renormalization Group Primer
A Renormalization Group Primer Physics 295 2010. Independent Study. Topics in Quantum Field Theory Michael Dine Department of Physics University of California, Santa Cruz May 2010 Introduction: Some Simple
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 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 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 informationSoft Collinear Effective Theory: An Overview
Soft Collinear Effective Theory: An Overview Sean Fleming, University of Arizona EFT09, February 1-6, 2009, Valencia Spain Background Before SCET there was QCD Factorization Factorization: separation of
More informationEffective Field Theory and. the Nuclear Many-Body Problem
Effective Field Theory and the Nuclear Many-Body Problem Thomas Schaefer North Carolina State University 1 Nuclear Effective Field Theory Low Energy Nucleons: Nucleons are point particles Interactions
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationIntroduction to spectroscopy on the lattice
Introduction to spectroscopy on the lattice Mike Peardon Trinity College Dublin JLab synergy meeting, November 21, 2008 Mike Peardon (TCD) Introduction to spectroscopy on the lattice November 21, 2008
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 informationIntroduction to Hamiltonian Monte Carlo Method
Introduction to Hamiltonian Monte Carlo Method Mingwei Tang Department of Statistics University of Washington mingwt@uw.edu November 14, 2017 1 Hamiltonian System Notation: q R d : position vector, p R
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 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 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 informationProperties of monopole operators in 3d gauge theories
Properties of monopole operators in 3d gauge theories Silviu S. Pufu Princeton University Based on: arxiv:1303.6125 arxiv:1309.1160 (with Ethan Dyer and Mark Mezei) work in progress with Ethan Dyer, Mark
More information