Home > Publications database > Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo > print

001     58568
005     20180211163830.0
024 7 _ |2 DOI
|a 10.1016/j.cpc.2008.08.006
024 7 _ |2 WOS
|a WOS:000262065100003
037 _ _ |a PreJuSER-58568
041 _ _ |a eng
082 _ _ |a 004
084 _ _ |2 WoS
|a Computer Science, Interdisciplinary Applications
084 _ _ |2 WoS
|a Physics, Mathematical
100 1 _ |a Cundy, N.
|b 0
|0 P:(DE-HGF)0
245 _ _ |a Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo
260 _ _ |a Amsterdam
|b North Holland Publ. Co.
|c 2009
300 _ _ |a 26 - 54
336 7 _ |a Journal Article
|0 PUB:(DE-HGF)16
|2 PUB:(DE-HGF)
336 7 _ |a Output Types/Journal article
|2 DataCite
336 7 _ |a Journal Article
|0 0
|2 EndNote
336 7 _ |a ARTICLE
|2 BibTeX
336 7 _ |a JOURNAL_ARTICLE
|2 ORCID
336 7 _ |a article
|2 DRIVER
440 _ 0 |a Computer Physics Communications
|x 0010-4655
|0 1439
|y 1
|v 180
500 _ _ |a Record converted from VDB: 12.11.2012
520 _ _ |a The computational costs of calculating the matrix sign function of the overlap operator together with fundamental numerical problems related to the discontinuity of the sign function in the kernel eigenvalues are the major obstacle towards simulations with dynamical overlap fermions using the Hybrid Monte Carlo algorithm. In a previous paper of the present series we introduced optimal numerical approximation of the sign function and have developed highly advanced preconditioning and relaxation techniques which speed up the inversion of the overlap operator by nearly an order of magnitude.In this fourth paper of the series we construct an HMC algorithm for overlap fermions. We approximate the matrix sign function using the Zolotarev rational approximation, treating the smallest eigenvalues of the Wilson operator exactly within the fermionic force. Based on this we derive the fermionic force for the overlap operator. We explicitly solve the problem of the Dirac delta-function terms arising through zero crossings of eigenvalues of the Wilson operator. The main advantage of scheme is that its energy violations scale better than O(Delta tau(2)) and thus are comparable with the violations of the standard leapfrog algorithm over the course of a trajectory. We explicitly prove that our algorithm satisfies reversibility and area conservation. We present test results from our algorithm on 4(4), 6(4), and 8(4) lattices. (C) 2008 Elsevier B.V. All rights reserved.
536 _ _ |a Scientific Computing
|c P41
|2 G:(DE-HGF)
|0 G:(DE-Juel1)FUEK411
|x 0
588 _ _ |a Dataset connected to Web of Science
650 _ 7 |a J
|2 WoSType
653 2 0 |2 Author
|a Lattice quantum chromodynamics
653 2 0 |2 Author
|a Overlap fermions
653 2 0 |2 Author
|a Hybrid Monte Carlo
700 1 _ |a Krieg, S.
|b 1
|u FZJ
|0 P:(DE-Juel1)132171
700 1 _ |a Arnold, G.
|b 2
|u FZJ
|0 P:(DE-Juel1)VDB57262
700 1 _ |a Frommer, A.
|b 3
|0 P:(DE-HGF)0
700 1 _ |a Lippert, T.
|b 4
|u FZJ
|0 P:(DE-Juel1)132179
700 1 _ |a Schilling, K.
|b 5
|0 P:(DE-HGF)0
773 _ _ |a 10.1016/j.cpc.2008.08.006
|g Vol. 180, p. 26 - 54
|p 26 - 54
|q 180<26 - 54
|0 PERI:(DE-600)1466511-6
|t Computer physics communications
|v 180
|y 2009
|x 0010-4655
856 7 _ |u http://dx.doi.org/10.1016/j.cpc.2008.08.006
909 C O |o oai:juser.fz-juelich.de:58568
|p VDB
913 1 _ |k P41
|v Scientific Computing
|l Supercomputing
|b Schlüsseltechnologien
|0 G:(DE-Juel1)FUEK411
|x 0
914 1 _ |y 2009
915 _ _ |0 StatID:(DE-HGF)0010
|a JCR/ISI refereed
920 1 _ |0 I:(DE-Juel1)JSC-20090406
|k JSC
|l Jülich Supercomputing Centre
|g JSC
|x 0
920 1 _ |0 I:(DE-82)080012_20140620
|k JARA-HPC
|l Jülich Aachen Research Alliance - High-Performance Computing
|g JARA
|x 1
970 _ _ |a VDB:(DE-Juel1)92268
980 _ _ |a VDB
980 _ _ |a ConvertedRecord
980 _ _ |a journal
980 _ _ |a I:(DE-Juel1)JSC-20090406
980 _ _ |a I:(DE-82)080012_20140620
980 _ _ |a UNRESTRICTED
981 _ _ |a I:(DE-Juel1)VDB1346

LibraryCollectionCLSMajorCLSMinorLanguageAuthor
Marc 21