Concept of the Energy-Dependent Temperature for Direct Langevin Dynamics Simulations of Rare Events in Systems with Arbitrary High Energy Barriers

We suggest an algorithm which allows single-stage direct Langevin dynamics simulations of transitions over arbitrarily high energy barriers. For this purpose, we propose a concept of the energy-dependent temperature (EDT): near the energy minima this temperature is high, but it tends toward room temperature for energies approaching the barrier value. In the resulting algorithm simulation time required for the computation of the escape rate over the barrier does not increase with barrier height. Switching times computed via our EDT algorithm agree very well with those obtained with the forward flux sampling (FFS). As the simulation time required by our method does not increase with the energy barrier, we achieve a very large speed-up compared even to our strongly optimized version of FFS (and all other multistage algorithms). In addition, our approach is free from the instability occurring in all multistage "climbing" methods where a product of a large number of transition probabilities between the interfaces must be computed.