Mass-scaling replica-exchange molecular dynamics optimizes computational resources with simpler algorithm

We develop a novel method of replica-exchange molecular dynamics (REMD) simulation, mass-scaling REMD (MSREMD) method, which improves trajectory accuracy at high temperatures, and thereby contributes to numerical stability. In addition, the MSREMD method can also simplify a replica-exchange routine by eliminating velocity scaling. As a pilot system, a Lennard-Jones fluid is simulated with the new method. The results show that the MSREMD method improves the trajectory accuracy at high temperatures compared with the conventional REMD method. We analytically demonstrate that the MSREMD simulations can reproduce completely the same trajectories of the conventional REMD ones with shorter time steps at high temperatures in case of the Nose-Hoover thermostats. Accordingly, we can easily compare the computational costs of the REMD and MSREMD simulations. We conclude that the MSREMD method decreases the instability and optimizes the computational resources with simpler algorithm under the constant trajectory accuracy at all temperatures.

Application of Simulated Tempering and Magnetizing to a Two-Dimensional Potts Model

We applied the simulated tempering and magnetizing (STM) method to the two-dimensional three-state Potts model in an external magnetic field in order to perform further investigations of the STMs applicability. The temperature as well as the external field are treated as dynamical variables updated during the STM simulations. After we obtained adequate information for several lattice sizes $L$ (up to $160times 160$), we also performed a number of conventional canonical simulations of large lattices, especially in order to illustrate the crossover behavior of the Potts model in external field with increasing $L$. The temperature and external field for larger lattice size simulations were chosen by extrapolation of the detail information obtained by STM. We carefully analyzed the crossover scaling at the phase transitions with respect to the lattice size as well as the temperature and external field. The crossover behavior is clearly observed in the simulations in agreement with theoretical predictions.

Simulated Tempering and Magnetizing: An Application of Two-Dimensional Simulated Tempering to Two-Dimensional Ising Model and Its Crossover

We performed two-dimensional simulated tempering (ST) simulations of the two-dimensional Ising model with different lattice sizes in order to investigate the two-dimensional STs applicability to dealing with phase transitions and to study the crossover of critical scaling behavior. The external field, as well as the temperature, was treated as a dynamical variable updated during the simulations. Thus, this simulation can be referred to as Simulated Tempering and Magnetizing (STM). We also performed the Simulated Magnetizing (SM) simulations, in which the external field was considered as a dynamical variable and temperature was not. As has been discussed by previous studies, the ST method is not always compatible with first-order phase transitions. This is also true in the magnetizing process. Flipping of the entire magnetization did not occur in the SM simulations under $T_mathrm{c}$ in large lattice-size simulations. However, the phase changed through the high temperature region in the STM simulations. Thus, the dimensional extension let us eliminate the difficulty of the first-order phase transitions and study wide area of the phase space. We then discuss how frequently parameter-updating attempts should be made for optimal convergence. The results favor frequent attempts. We finally study the crossover behavior of the phase transitions with respect to the temperature and external field. The crossover behavior was clearly observed in the simulations in agreement with the theoretical implications.