The potential energy surface for the model system. The potential energy function is given in units of and takes a functionThe separation of contour lines (different colors) is . The grids are given for illustration. In the calculations the size of the small areas is .
The convergence of with different initial guesses.
The calculated free energy surface for the model system. The free energy (in units of and the separation between contour lines is ) was calculated using the non-Boltzmann distribution approach using 40 different temperatures ranging between and .
The ratio of the successful over total trajectories increases with the adaptive kinetic sampling steps. The trajectory length is .
The natural logarithm of visiting probability by the successful trajectories (calculated after reweighting of their initial configurations according to a Boltzmann distribution at the temperature ). The separation between contour lines is .
An example of the sampled spatial points (black) and those that lead to successful trajectories (red). The sampling thus not allows calculation of thermodynamic properties at the desired temperature (the points at the center) but also covers with a high probability that reaction portion of the configuration space.
The calculated rate as a function of the length of trajectories (the trajectories are obtained with a constant energy using Newtonian equations).
The convergence of for a simulation for maltose. It can be seen that they converge in a short simulation of a few picoseconds. In comparison, their convergence took about when the previous updating method was used (Ref. 26).
(a) The sampled energy for a large system (kinase in explicit solvent with 19 179 atoms) using the enhanced sampling method (dashed line) and normal MD simulations (solid line). (b) The energy distribution is significantly broadened in the enhanced sampling simulations (dashed line) compared to the normal MD (solid line).
The potential of mean force (dashed line) as a function of the coordinate, which shows a lower barrier compared to that obtained along the pathway with the lower barrier (calculated using configurations with ).
Article metrics loading...
Full text loading...