We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow Hamiltonian $til
de{H}$. This property allows for a new tuning method of the molecular dynamics integrator and also allows for a new class of integrators (force-gradient integrators) which is expected to reduce significantly the computational cost of future large-scale gauge field ensemble generation.
We extract to high statistical precision an excited spectrum of single-particle isoscalar mesons using lattice QCD, including states of high spin and, for the first time, light exotic JPC isoscalars. The use of a novel quark field construction has en
abled us to overcome the long-standing challenge of efficiently including quark-annihilation contributions. Hidden-flavor mixing angles are extracted and while most states are found to be close to ideally flavor mixed, there are examples of large mixing in the pseudoscalar and axial sectors in line with experiment. The exotic JPC isoscalar states appear at a mass scale comparable to the exotic isovector states.
We show how to improve the molecular dynamics step of Hybrid Monte Carlo, both by tuning the integrator using Poisson brackets measurements and by the use of force gradient integrators. We present results for moderate lattice sizes.
