Abstract
Concepts and applications of lattice quantum chromodynamics (LQCD) are introduced. After discussing how to define quarks and gluons on the Euclidean hypercubic lattice, the strong coupling expansion and the weak coupling expansions are reviewed to see the vital role played by the quantum fluctuations in QCD. Fundamental techniques for numerical LQCD simulations such as the Markov Chain Monte Carlo method and the Hybrid Monte Carlo method are discussed in some details. As examples of the high precision LQCD simulations, numerical results of the heavy quark potential and the hadron masses are shown. Recent LQCD results on the baryon-baryon interactions are briefly discussed.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Similar content being viewed by others
Notes
- 1.
Rigorous definitions are as follows. (1) The Markov chain is said to be irreducible if one can find a finite positive integer n( < ∞) such that (T n) ϕ ϕ′ > 0 for all ϕ and ϕ′. (2) The period, d(ϕ), is defined by the greatest common divisor of the set of positive integers n( ≥ 1) such that (T n) ϕ ϕ > 0 is satisfied. If d(ϕ) = 1 for all ϕ, the Markov chain is said to be aperiodic [10].
References
K.G. Wilson, Phys. Rev. D 10, 2445 (1974)
M. Creutz, Phys. Rev. D 21, 2308 (1980)
N. Brambilla et al., Eur. Phys. J. C 74, 2981 (2014)
K.G. Wilson, Nucl. Phys. Proc. Suppl. 140, 3 (2005)
M. Creutz, Quarks, Gluons and Lattices (Cambridge University Press, Cambridge, 1985)
H.J. Rothe, World Sci. Lect. Notes Phys. 82, 1 (2012)
C. Hoelbling, Acta Phys. Polon. B 45, 2143 (2014)
A. Ukawa, J. Stat. Phys. 160, 1081 (2015)
J.W. Negele, H. Orland, Quantum Many Particle Systems (Addison-Wesley, Redwood, 1988)
O. Häggström, Finite Markov Chains and Algorithmic Applications (Cambridge University Press, Cambridge, 2002)
H. Suwa, S. Todo, Phys. Rev. Lett. 105, 120603 (2010)
S. Duane, A.D. Kennedy, B.J. Pendleton, D. Roweth, Phys. Lett. B 195, 216 (1987)
N. Metropolis, A.W. Rosenbluth, M.N. Rosenbluth, A.H. Teller, E. Teller, J. Chem. Phys. 21, 1087 (1953)
S. Schaefer, PoS LATTICE 2012, 001 (2012)
G.S. Bali, Phys. Rep. 343, 1 (2001)
S. Durr et al., Science 322, 1224 (2008)
S. Borsanyi et al., Science 347, 1452 (2015)
The Review of Particle Physics (2015). http://pdg.lbl.gov/
M. Asakawa, T. Hatsuda, Y. Nakahara, Prog. Part. Nucl. Phys. 46, 459 (2001)
Z. Fodor, C. Hoelbling, Rev. Mod. Phys. 84, 449 (2012)
R. Machleidt, Int. J. Mod. Phys. E26, 1740018 (2017)
S. Aoki et al., Eur. Phys. J. C 74, 2890 (2014)
M. Lüscher, Nucl. Phys. B 354, 531 (1991)
N. Ishii, S. Aoki, T. Hatsuda, Phys. Rev. Lett. 99, 022001 (2007)
N. Ishii et al., HAL QCD collaboration. Phys. Lett. B 712, 437 (2012)
T. Iritani, HAL QCD collaboration. JHEP 1610, 101 (2016)
S. Okubo, R.E. Marshak, Ann. Phys. 4, 166 (1958)
M. Oka, K. Shimizu, K. Yazaki, Prog. Theor. Phys. Suppl. 137, 1 (2000)
T. Inoue et al., HAL QCD collaboration. Nucl. Phys. A 881, 28 (2012)
T. Inoue et al., HAL QCD collaboration. Phys. Rev. Lett. 111, 112503 (2013)
A. Akmal, V.R. Pandharipande, D.G. Ravenhall, Phys. Rev. C 58, 1804 (1998)
T. Doi et al., HAL QCD collaboration. PoS LATTICE2015, 086 (2016)
Acknowledgements
The author thanks Takumi Doi and Atsushi Nakamura for useful comments and information on various aspects of LQCD simulations. He also thank the members of HAL QCD Collaboration for fruitful discussions on the hadron-hadron interactions on the lattice. This work was supported in part by MEXT SPIRE and JICFuS and also by RIKEN iTHES Project.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Appendix
Appendix
3.1.1 Four Vectors and Dirac Matrices
In the (3+1)-dimensional Minkowski spacetime, coordinates, derivatives and four vectors with μ = 0, 1, 2, 3 are
In the 4-dimensional Euclidean space, we define the corresponding vectors for μ = 4, 1, 2, 3 as
In the (3+1)-dimensional Minkowski spacetime with the metric g μ ν = diag(1, −1, −1, −1), the Dirac matrices satisfy the following relations for μ = 0, 1, 2, 3,
In the standard Dirac representation, we have
where σ j are the Pauli matrices; \(\sigma _{1} = \left (\begin{array}{cc} 0&1\\ 1 &0\\ \end{array} \right ),\sigma _{2} = \left (\begin{array}{cc} 0& - i\\ i & 0\\ \end{array} \right ),\sigma _{3} = \left (\begin{array}{cc} 1& 0\\ 0 & -1\\ \end{array} \right ).\)
In the 4-dimensional Euclidean space with the metric δ μ ν = diag(1, 1, 1, 1), we define the Euclidean Dirac matrices as
which satisfy the relations,
3.1.2 SU(N) Algebra
Let \(\mathcal{T}^{a}\) (a = 1, ⋯ , N 2 − 1) are the Hermitian generators of the SU(N) group. They satisfy the Lie algebra
where f abc is the structure constant being totally anti-symmetric in its indices. \((\mathcal{T}^{b})^{2}\) commutes with every generator \(\mathcal{T}^{a}\) and is called the quadratic Casimir operator.
For N = 2, f abc reduces to the anti-symmetric tensor ε ijk with ε 123 = 1. For N = 3, the non-vanishing components of f abc read \(f_{123} = 1,f_{147} = -f_{156} = f_{246} = f_{257} = f_{345} = -f_{367} = 1/2,f_{458} = f_{678} = \sqrt{3}/2\).
In the fundamental representation, \(\mathcal{T}^{a}\) is written by the N × N matrices t a as
where λ a for N = 2 reduce to the Pauli matrices σ i , while those for N = 3 reduce to the Gell-Mann matrices.
Some useful relations of t a for general N are
In the adjoint representation, \(\mathcal{T}^{a}\) is written by (N 2 − 1) × (N 2 − 1) matrices T a as
which satisfy the relations
3.1.3 Gaussian and Grassmann Integrals
Basic Gaussian and Grassmann integrals are
Here x (z) is a real (complex) number, while \(\bar{\xi }\) and ξ are anti-commuting Grassmann numbers (\(\{\xi,\bar{\xi }\}= 0\), and \(\xi ^{2} =\bar{\xi } ^{2} = 0\)). a and b are assumed to be real and positive numbers, while c is an arbitrary complex number. Equation (3.109) can be shown by rewriting the integral in terms of the real and imaginary parts of z or in terms of the polar coordinates of z. Equation (3.110) can be shown by noting that \(\mathrm{e}^{-c\bar{\xi }\xi } = 1 - c\bar{\xi }\xi\) and ∫ d ξ = ∂∕∂ ξ (integral = derivative) for Grassmann variables.
Generalization of the above results to the case of multiple variables is straightforward. For x = (x 1, ⋯ , x n ), z = (z 1, ⋯ , z n ), ξ = (ξ 1, ⋯ , ξ n ), and \(\bar{\xi }= (\bar{\xi }_{1},\cdots \,,\bar{\xi }_{n})\) with \(\{\xi _{k},\xi _{l}\} =\{\bar{\xi } _{k},\bar{\xi }_{l}\} =\{\xi _{k},\bar{\xi }_{l}\} = 0\), we have
Here A is a non-singular and real-symmetric matrix whose eigenvalues a l satisfy a l > 0 for all l. B is a non-singular complex matrix whose complex eigenvalues b l obtained by the biunitary transformation (UBV †) satisfy Re b l > 0 for all l. C is an arbitrary complex matrix with no conditions. Note that B and C do not have to be Hermitian matrices. In field theories, the label “l” summarizes all possible indices including spin, flavor, color, spacetime points etc. and “Det” denotes the determinant for all these indices.
3.1.4 Method of Characteristics
We need to construct a general solution of the following partial differential equation,
For this purpose, we introduce the running coupling \(\bar{g}(\lambda )\) through \(\lambda d\bar{g}/d\lambda = -\beta (\bar{g})\) whose formal solution reads
Then the solution of Eq. (3.114) can be written as
This can be explicitly checked by applying the partial derivative on both sides,
where we have used the relation \(\partial \bar{g}/\partial g =\beta (\bar{g})/\beta (g)\) obtained from Eq. (3.115).
In general, the first-order partial differential equation (PDE) can be transformed to a set of ordinary differential equations (ODEs) and can be solved by the method of characteristics. As an illustration, let us consider the following PDE,
This is equivalent to the coupled ODEs,
where s parametrizes the “flow” of the coordinates. \((\bar{t}(s),\bar{x}(s))\). This is called the characteristic curve as shown in Fig. 3.15.
The function f can be obtained by integrating the last equation of Eq. (3.119) on the characteristic curve from the initial point (t i, x i) to the final point (t F, x F) ≡ (t, x),
where h stands for an integration of the known function \(c(\bar{t},\bar{x})\) on the characteristic curve. Equation (3.120) implies that the desired function at (t, x) is obtained essential by a “pullback” of the point to (t i, x i) along the characteristic curve. It is a straightforward exercise to generalize the above derivation to the system with more coordinates, (t, x).
3.1.5 Leapfrog Integrator in Molecular Dynamics
Let us start with a Tayler expansion of the field ϕ:
where we have used the equation of motion, \(\dot{\phi }(s) \equiv d\phi (s)/ds =\pi (s)\). To evaluate π(s +ɛ∕2), we take the midpoint prescription which does not have O(ɛ 2) error,
Equations (3.121) and (3.122) give a procedure to move the molecular dynamics one-step forward, (ϕ(s), π(s −ɛ∕2)) → (ϕ(s +ɛ), π(s +ɛ∕2)). The initial and final steps need to receive special care,
which have only O(ε 2) accuracy. An illustration of this leapfrog integrator is shown in Fig. 3.16. Since the initial and final steps introduce O(ɛ 2) error irrespective of the length of the MD trajectory, and the intermediate steps introduce O(ɛ 3) ×ɛ −1 = O(ɛ 2) error as a whole, one finds Δ H = O(ɛ 2) after one MD trajectory before the Metropolis test.
The leapfrog integrator satisfies the reversibility and symplectic property, which can be checked explicitly by using the above definitions (Exercise 3.13).
Rights and permissions
Copyright information
© 2017 Springer International Publishing AG
About this chapter
Cite this chapter
Hatsuda, T. (2017). Lattice Quantum Chromodynamics. In: Hjorth-Jensen, M., Lombardo, M., van Kolck, U. (eds) An Advanced Course in Computational Nuclear Physics. Lecture Notes in Physics, vol 936. Springer, Cham. https://doi.org/10.1007/978-3-319-53336-0_3
Download citation
DOI: https://doi.org/10.1007/978-3-319-53336-0_3
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-53335-3
Online ISBN: 978-3-319-53336-0
eBook Packages: Physics and AstronomyPhysics and Astronomy (R0)