Cartoon of the exchange of interfaces and state swaps for a 4-state system, where each state has two interfaces. The initial path in the left top corner crosses not only the current (black) interface, but also the outermost interface so an exchange is allowed. A shooting move creates a path toward another state, so that state swap is possible. The current interface is then the outermost interface of the new state. Another interface exchange and a subsequent shot creates the red path in the lower left corner cartoon.
Concept of the adaptive Wang-Landau MSTIS scheme. Sampling the outermost interface of the initial state (left) creates an activated path that finds another basin, which results in a very long trial trajectory (middle). This path is analyzed, and if the new state is stable it is treated as a new state with its own set of interfaces. This process is repeated, and new states are added to the database as needed (right).
The 2D potential with the 6 minima plotted in the x,y plane. Samples of the Wang-Landau MSTIS paths at β = 4 are drawn on top. The dynamically unbiased paths connect all states.
Decorrelation in the Wang-Landau MSTIS sampling. Top: The state index as a function of simulation time measured in number of shooting moves. The dotted line shows the results when the DOP is initialized to zero and solid line shows the result for a converged DOP. Note that the decorrelation is almost equal. Bottom: The decorrelation between stable state minima. The main window shows the histogram of interval lengths between stable state visits (replica 0). The inset shows the histogram of interval lengths between stable state visits for the whole ensemble. Note that this interval is naturally much shorter than the one between replica 0, as the replica has to diffuse through the whole interface space.
Left: Convergence of the DOP in the course of the Wang-Landau MSTIS simulation. The replica index is defined as Im + i. Middle: Comparison of the logarithmic crossing probability ln P(λ) for the six states as a function of λ between the replica exchange MSTIS (dotted curves) and Wang-Landau MSTIS (solid curves) simulations. For each state the two approaches result in an identical crossing probability within the error bar (shown for 2 points). Right: The density of paths in the Wang-Landau algorithm converges to crossing probability.
The RPE contains both information about the free energy and the kinetics. Left: FE projections of the RPE in the x, y plane. Note that this projection is naturally equal to the potential V(x, y). Right: Projection of the committor p 3 in the x,y plane. Purple denotes p 3 = 0 and white denotes p 3 = 1.
Convergence of the equilibrium populations as a function of simulation time. Left: Fixed bias. Right: Wang-Landau MSTIS scheme.
The configurations corresponding to the stable states in the LJ7 cluster. From left to right: the pentagonal bipiramid, the capped octahedron, bicapped trigonal bipyramid, and tricapped tetrahedron. 38,39
Buildup of the density of paths for the LJ7 cluster using the adaptive Wang-Landau MSTIS scheme. The index is defined as Im + i.
Converged crossing probability (curves) and density of paths (symbols) for the LJ7 cluster.
Tetrahedral projection of the trajectories consisting of the 4 square distance sum order parameters. The projection is such that the corners represent exact correspondence to the target structure, whereas the center is infinitely away from all structures.
The molecular structure of alanine dipeptide rendered in licorice. Carbons in cyan, oxygen in red, nitrogen in blue, and hydrogen in white. The two order parameters describing the metastable states are the dihedral angles ϕ and ψ.
The free energy surface from a replica exchange MD simulation obtained by projecting the logarithm of the probability to find a (ϕ, ψ) pair. 43
Logarithm of crossing probability against interface for states β, α, E, and F. The red, green, and blue lines are the last three DOP biasing functions, while the black line with circles is the final logarithm of the crossing probability.
Sorted squared distance sum for the four minima.
(a) Outermost interface transition matrix for the four states. Subscripts denote the error in the last digits. Columns are leaving states and rows are arriving states. (b) Flux through the first interface for the four states.
Rate matrix for the four states. Subscripts denote the error in the last digits. Columns are leaving states and rows are arriving states.
Definition of interfaces for states β, α, E, and F.
Rate constant data. (a) Flux through first interface, crossing probability from the first to the last interface, and flux at the outermost interface for each state. (b) Transition probability for states β, α, E, and F at their respective outer interfaces. The rows denote the leaving state and the columns denote the arriving state. (c) Rate constant matrix from single replica MSTIS for states β, α, E, and F.
Article metrics loading...
Full text loading...