A separable shadow Hamiltonian hybrid Monte Carlo method.

Hybrid Monte Carlo (HMC) is a rigorous sampling method that uses molecular dynamics (MD) as a global Monte Carlo move. The acceptance rate of HMC decays exponentially with system size. The shadow hybrid Monte Carlo (SHMC) was previously introduced to reduce this performance degradation by sampling instead from the shadow Hamiltonian defined for MD when using a symplectic integrator. SHMC's performance is limited by the need to generate momenta for the MD step from a nonseparable shadow Hamiltonian. We introduce the separable shadow Hamiltonian hybrid Monte Carlo (S2HMC) method based on a formulation of the leapfrog/Verlet integrator that corresponds to a separable shadow Hamiltonian, which allows efficient generation of momenta. S2HMC gives the acceptance rate of a fourth order integrator at the cost of a second-order integrator. Through numerical experiments we show that S2HMC consistently gives a speedup greater than two over HMC for systems with more than 4000 atoms for the same variance. By comparison, SHMC gave a maximum speedup of only 1.6 over HMC. S2HMC has the additional advantage of not requiring any user parameters beyond those of HMC. S2HMC is available in the program PROTOMOL 2.1. A Python version, adequate for didactic purposes, is also in MDL (http://mdlab.sourceforge.net/s2hmc).

MDSIMAID: automatic parameter optimization in fast electrostatic algorithms.

MDSIMAID is a recommender system that optimizes parallel Particle Mesh Ewald (PME) and both sequential and parallel multigrid (MG) summation fast electrostatic solvers. MDSIMAID optimizes the running time or parallel scalability of these methods within a given error tolerance. MDSIMAID performs a run time constrained search on the parameter space of each method starting from semiempirical performance models. Recommended parameters are presented to the user. MDSIMAID's optimization of MG leads to configurations that are up to 14 times faster or 17 times more accurate than published recommendations. Optimization of PME can improve its parallel scalability, making it run twice as fast in parallel in our tests. MDSIMAID and its Python source code are accessible through a Web portal located at http://mdsimaid.cse.nd.edu.