1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Adaptive single replica multiple state transition interface sampling
Rent:
Rent this article for
USD
10.1063/1.4813777
/content/aip/journal/jcp/139/4/10.1063/1.4813777
http://aip.metastore.ingenta.com/content/aip/journal/jcp/139/4/10.1063/1.4813777

Figures

Image of FIG. 1.
FIG. 1.

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.

Image of FIG. 2.
FIG. 2.

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).

Image of FIG. 3.
FIG. 3.

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.

Image of FIG. 4.
FIG. 4.

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.

Image of FIG. 5.
FIG. 5.

Left: Convergence of the DOP in the course of the Wang-Landau MSTIS simulation. The replica index is defined as + . Middle: Comparison of the logarithmic crossing probability ln (λ) 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.

Image of FIG. 6.
FIG. 6.

The RPE contains both information about the free energy and the kinetics. Left: FE projections of the RPE in the , plane. Note that this projection is naturally equal to the potential (, ). Right: Projection of the committor in the x,y plane. Purple denotes = 0 and white denotes = 1.

Image of FIG. 7.
FIG. 7.

Convergence of the equilibrium populations as a function of simulation time. Left: Fixed bias. Right: Wang-Landau MSTIS scheme.

Image of FIG. 8.
FIG. 8.

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.

Image of FIG. 9.
FIG. 9.

Buildup of the density of paths for the LJ7 cluster using the adaptive Wang-Landau MSTIS scheme. The index is defined as + .

Image of FIG. 10.
FIG. 10.

Converged crossing probability (curves) and density of paths (symbols) for the LJ7 cluster.

Image of FIG. 11.
FIG. 11.

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.

Image of FIG. 12.
FIG. 12.

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 ψ.

Image of FIG. 13.
FIG. 13.

The free energy surface from a replica exchange MD simulation obtained by projecting the logarithm of the probability to find a (ϕ, ψ) pair.

Image of FIG. 14.
FIG. 14.

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.

Tables

Generic image for table
Table I.

Sorted squared distance sum for the four minima.

Generic image for table
Table II.

(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.

Generic image for table
Table III.

Rate matrix for the four states. Subscripts denote the error in the last digits. Columns are leaving states and rows are arriving states.

Generic image for table
Table IV.

Definition of interfaces for states β, α, E, and F.

Generic image for table
Table V.

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.

Loading

Article metrics loading...

/content/aip/journal/jcp/139/4/10.1063/1.4813777
2013-07-24
2014-04-18
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Adaptive single replica multiple state transition interface sampling
http://aip.metastore.ingenta.com/content/aip/journal/jcp/139/4/10.1063/1.4813777
10.1063/1.4813777
SEARCH_EXPAND_ITEM