Polynomial filtered HMC - An algorithm for lattice QCD with dynamical quarks

dc.contributor.authorKamleh, W.
dc.contributor.authorPeardon, M.
dc.date.issued2012
dc.description.abstractPolynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in Hybrid Monte Carlo (HMC) simulations. The use of a multiple time-scale integration scheme then allows the filtered pseudofermions to be evolved using a coarse step size. The algorithm is applied to <sup>Nf</sup>=2 flavour simulations at two different quark masses. Polynomial filtering yields a reduction in the computational expense of the molecular dynamics integration of between three- and fivefold when compared to standard HMC. Significantly, the observed speedup is better at the lower quark mass. © 2012 Elsevier B.V. All rights reserved.
dc.description.statementofresponsibilityWaseem Kamleh, Mike Peardon
dc.identifier.citationComputer Physics Communications, 2012; 183(9):1993-2000
dc.identifier.doi10.1016/j.cpc.2012.05.002
dc.identifier.issn0010-4655
dc.identifier.issn1879-2944
dc.identifier.orcidKamleh, W. [0000-0002-6177-5366]
dc.identifier.urihttp://hdl.handle.net/2440/74371
dc.language.isoen
dc.publisherElsevier Science BV
dc.rightsCrown Copyright © 2012 Published by Elsevier B.V. All rights reserved.
dc.source.urihttps://doi.org/10.1016/j.cpc.2012.05.002
dc.subjectLattice QCD
dc.subjectHybrid Monte Carlo
dc.subjectAlgorithms
dc.titlePolynomial filtered HMC - An algorithm for lattice QCD with dynamical quarks
dc.typeJournal article
pubs.publication-statusPublished

Files