A modified Cayley transform for SU(3) molecular dynamics simulations

Abstract

We propose a modification to the Cayley transform that defines a suitable local parameterization for the special unitary group $\mathrm{SU}(3)$. The new mapping is used to construct splitting methods for separable Hamiltonian systems whose phase space is the cotangent bundle of $\mathrm{SU}(3)$ or, more general, ${\mathrm{SU}(3)}^N$, $N\in \mathbb{N}$. Special attention is given to the hybrid Monte Carlo algorithm for gauge field generation in lattice quantum chromodynamics. We show that the use of the modified Cayley transform instead of the matrix exponential neither affects the time-reversibility nor the volume-preservation of the splitting method. Furthermore, the advantages and disadvantages of the Cayley-based algorithms are discussed and illustrated in pure gauge field simulations.

15 pages, 4 figures

References in corpus (5)

• Testing and tuning symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Solving the matrix exponential function for the groups SU(3), SU(4) and Sp(2)
• Hessian-free force-gradient integrators
• Hessian-free force-gradient integrators and their application to lattice QCD simulations