Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics.

Langevin dynamics has become a popular tool to simulate the Boltzmann equilibrium distribution. When the repartition of the Langevin equation involves the exact realization of the Ornstein-Uhlenbeck noise, in addition to the conventional density evolution, there exists another type of discrete evolution that may not correspond to a continuous, real dynamical counterpart. This virtual dynamics case is also able to produce the desired stationary distribution. Different types of repartition lead to different numerical schemes, of which the accuracy and efficiency are investigated through studying the harmonic oscillator potential, an analytical solvable model. By analyzing the asymptotic distribution and characteristic correlation time that are derived by either directly solving the discrete equations of motion or using the related phase space propagators, it is shown that the optimal friction coefficient resulting in the minimum characteristic correlation time depends on the time interval chosen in the numerical implementation. When the recommended "middle" scheme is employed, both analytical and numerical results demonstrate that, for good numerical performance in efficiency as well as accuracy, one may choose a friction coefficient in a wide range from around the optimal value to the high friction limit.

A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat.

We introduce a novel simple algorithm for thermostatting path integral molecular dynamics (PIMD) with the Langevin equation. The staging transformation of path integral beads is employed for demonstration. The optimum friction coefficients for the staging modes in the free particle limit are used for all systems. In comparison to the path integral Langevin equation thermostat, the new algorithm exploits a different order of splitting for the phase space propagator associated to the Langevin equation. While the error analysis is made for both algorithms, they are also employed in the PIMD simulations of three realistic systems (the H2O molecule, liquid para-hydrogen, and liquid water) for comparison. It is shown that the new thermostat increases the time interval of PIMD by a factor of 4-6 or more for achieving the same accuracy. In addition, the supplementary material shows the error analysis made for the algorithms when the normal-mode transformation of path integral beads is used.

Exact results for the residual entropy of ice hexagonal monolayer.

Since the problem of the residual entropy of square ice was exactly solved, exact solutions for two-dimensional realistic ice models have been of interest. In this work, we study the exact residual entropy of ice hexagonal monolayer in two cases. In the case that the external electric field along the z-axis exists, we map the hydrogen configurations into the spin configurations of the Ising model on the kagome lattice. By taking the low temperature limit of the Ising model, we derive the exact residual entropy, which agrees with the result determined previously from the dimer model on the honeycomb lattice. In another case that the ice hexagonal monolayer is under the periodic boundary conditions in the cubic ice lattice, the residual entropy has not been studied exactly. For this case, we employ the six-vertex model on the square lattice to represent the hydrogen configurations obeying the ice rules. The exact residual entropy is obtained from the solution of the equivalent six-vertex model. Our work provides more examples of the exactly soluble two-dimensional statistical models.