Chiral Fermions Algorithms In Lattice QCD

  • Dafina Xhako Department of Physics Engineering, Faculty of Mathematical Engineering and Physical Engineering, Polytechnic University of Tirana, Tirana, Albania
  • Rudina Zeqirllari Department of Physics, Faculty of Natural Sciences, University of Tirana, Tirana, Albania
  • Artan Boriçi Department of Physics, Faculty of Natural Sciences, University of Tirana, Tirana,
Keywords: algorithms, chiral, critical-slowing, fermions, lattice, QCD


The theory that explains the strong interactions of the elementary particles, as part of the standard model, it is the so-called Quantum Chromodynamics (QCD) theory. In regimes of low energy this theory it is formulated and solved in a lattice with four dimensions using numerical simulations. This method it is called the lattice QCD theory. Quark propagator it the most important element that is calculated because it contains the physical information of lattice QCD. Computing quark propagator of chiral fermions in lattice means that we should invert the chiral Dirac operator, which has high complexity. In the standard inversion algorithms of the Krylov subspace methods, that are used in these kinds of simulations, the time of inversion is scaled with the inverse of the quark mass. In lattice QCD simulations with chiral fermions, this phenomenon it is knowing as the critical slowing-down problem. The purpose of this work is to show that the preconditioned GMRESR algorithm, developed in our previous work, solves this problem. The preconditioned GMRESR algorithm it is developed in U(1) group symmetry using QCDLAB 1.0 package, as good “environment” for testing new algorithms. In this paper we study the escalation of the time of inversion with the quark mass for this algorithm. It turned out that it is a fast inversion algorithm for lattice QCD simulations with chiral fermions, that “soothes” the critical slowing-down of standard algorithms. The results are compared with SHUMR algorithm that is optimal algorithm used in these kinds of simulations. The calculations are made for 100 statistically independent configurations on 64 x 64 lattice gauge U(1) field for three coupling constant and for some quark masses. The results showed that for the preconditioned GMRESR algorithm the coefficient k, related to the critical slowing down phenomena, it is approximately - 0.3 compared to the inverse proportional standard law (k = -1) that it is scaled SHUMR algorithm, even for dense lattices. These results make more stable and confirm the efficiency of our algorithm as an algorithm that avoid the critical slowing down phenomenon in lattice QCD simulations. In our future studies we have to develop the preconditioned GMRESR algorithm in four dimensions, in SU (3) lattice gauge theory.


Download data is not yet available.


M. Lüscher, Annales Henri Poincare, 4, 197-210 (2003).

K. G. Wilson, Phys. Rev. D, 10, 2445 (1974), doi: 10.1103/PhysRevD.10.2445.

D. B. Kaplan, Phys. Lett. B, 228, 342 (1992), doi: 10.1016/0370-2693(92)91112-M.

V. Furman and Y. Shamir, Nucl. Phys. B, 439, 54 (1995), doi: 10.1016/0550-3213(95)00031-M.

R. Narayanan and H. Neuberger, Phys. Lett. B 302, 62 (1993), doi: 10.1016/0370-2693(93)90636-V.

R. Narayanan and H. Neuberger, Nucl. Phys. B 443, 305 (1995), doi: 10.1016/0550-3213(95)00111-5.

A. Boriçi, in: QCD and Numerical Analysis III Proceedings, edited by A. Boriçi, A. Frommer, B. Joó, A. Kennedy and B. Pendleton (Springer, Germany, 2005), pp. 25-39.

A. Boriçi, Nucl. Phys. Proc. Suppl. 83, 771 (2000), doi: 10.1016/S0920-5632(00)91802-4.

A. Boriçi, in: Lattice Fermions and Structure of the Vacuum, edited by V. Mitrjushkin and G. Schierholz, (Kluwer Academic Publishers, London, 2006), pp. 41-54.

H. Neuberger, Phys. Lett. B 417, 141 (1998), doi: 10.1016/S0370-2693(97)01368-3.

A. Boriçi, PhD thesis, No. CSCS TR-96-27, ETH Zurich (1996).

A. Boriçi and P. Forcrand, in: Physics Computing ’94 Proceedings, edited by R. Gruber and M. Tomassini (European Physical Society, Switzerland, 1994), pp. 711–714.

D. Xhako and A. Boriçi, American Journal of Physics and Applications, 2, 67 (2014), doi: 10.11648/j.ajpa.20140202.15.

A. Boriçi, High Energy Physics - Lattice (hep-lat), (2006), arXiv:hep-lat/0610054.

A. Boriçi, Invited talk given at the: Domain Wall Fermions at Ten Years, 2007 March 15-17, Brookhaven National Laboratory, arXiv:hep-lat/0703021.

QCDLAB software,

N. Cundy, Comput. Phys. Commun. 165, 22 (1974), doi: 10.1016/j.cpc.2004.10.005.

D. Hyka (Xhako) and R. Osmanaj (Zeqirllari), in: 35th International Symposium on Lattice Field Theory - Lattice 2017 Proceedings, edited by M. D. Morte, E. G. Sanchez, P. Fritzsch and C. P. Ruano, (EPJ Web of Conferences, France, 2018), pp. 14005-14061.

0 article
How to Cite
Xhako, D., Zeqirllari, R., & Boriçi, A. (2019). Chiral Fermions Algorithms In Lattice QCD. East European Journal of Physics, (1), 34-39.