-
Enhanced path sampling using subtrajectory Monte Carlo moves
Authors:
Daniel T. Zhang,
Enrico Riccardi,
Titus S. van Erp
Abstract:
Path sampling allows the study of rare events like chemical reactions, nucleation and protein folding via a Monte Carlo (MC) exploration in path space. Instead of configuration points, this method samples short molecular dynamics (MD) trajectories with specific start- and end-conditions. As in configuration MC, its efficiency highly depends on the types of MC moves. Since the last two decades, the…
▽ More
Path sampling allows the study of rare events like chemical reactions, nucleation and protein folding via a Monte Carlo (MC) exploration in path space. Instead of configuration points, this method samples short molecular dynamics (MD) trajectories with specific start- and end-conditions. As in configuration MC, its efficiency highly depends on the types of MC moves. Since the last two decades, the central MC move for path sampling has been the so-called shooting move in which a perturbed phase point of the old path is propagated backward and forward in time to generate a new path. Recently, we proposed the subtrajectory moves, stone-skipping (SS) and web-throwing (WT), that are demonstrably more efficient. However, the one-step crossing requirement makes them somewhat more difficult to implement in combination with external MD programs or when the order parameter determination is expensive. In this article, we present strategies to address the issue. The most generic solution is a new member of subtrajectory moves, wire fencing (WF), that is less thrifty than the SS, but more versatile. This makes it easier to link path sampling codes with external MD packages and provides a practical solution for cases where the calculation of the order parameter is expensive or not a simple function of geometry. We demonstrate the WF move in a double well Langevin model, a thin film breaking transition based on classical force fields, and a smaller ruthenium redox reaction at the ab initio level in which the order parameter explicitly depends on the electron density.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
Exchanging replicas with unequal cost, infinitely and permanently
Authors:
Sander Roet,
Daniel T. Zhang,
Titus S. van Erp
Abstract:
We developed a replica exchange method that is effectively parallelizable even if the computational cost of the Monte Carlo moves in the parallel replicas are considerably different, for instance, because the replicas run on different type of processor units or because of the algorithmic complexity. To prove detailed-balance, we make a paradigm shift from the common conceptual viewpoint in which t…
▽ More
We developed a replica exchange method that is effectively parallelizable even if the computational cost of the Monte Carlo moves in the parallel replicas are considerably different, for instance, because the replicas run on different type of processor units or because of the algorithmic complexity. To prove detailed-balance, we make a paradigm shift from the common conceptual viewpoint in which the set of parallel replicas represents a high-dimensional superstate, to an ensemble based criterion in which the other ensembles represent an environment that might or might not participate in the Monte Carlo move. In addition, based on a recent algorithm for computing permanents, we effectively increase the exchange rate to infinite without the steep factorial scaling as function of the number of replicas. We illustrate the effectiveness of the replica exchange methodology by combining it with a quantitative path sampling method, replica exchange transition interface sampling (RETIS), in which the costs for a Monte Carlo move can vary enormously as paths in a RETIS algorithm do not have the same length and the average path lengths tend to vary considerably for the different path ensembles that run in parallel. This combination, coined $\infty$RETIS, was tested on three model systems.
△ Less
Submitted 25 May, 2022;
originally announced May 2022.
-
The dynamics of the DNA denaturation transition
Authors:
Titus S. van Erp,
Michel Peyrard
Abstract:
The dynamics of the DNA denaturation is studied using the Peyrard-Bishop-Dauxois model. The denaturation rate of double stranded polymers decreases exponentially as function of length below the denaturation temperature. Above Tc, the rate shows a minimum, but then increases as function of length. We also examine the influence of sequence and solvent friction. Molecules having the same number of we…
▽ More
The dynamics of the DNA denaturation is studied using the Peyrard-Bishop-Dauxois model. The denaturation rate of double stranded polymers decreases exponentially as function of length below the denaturation temperature. Above Tc, the rate shows a minimum, but then increases as function of length. We also examine the influence of sequence and solvent friction. Molecules having the same number of weak and strong base-pairs can have significantly different opening rates depending on the order of base-pairs.
△ Less
Submitted 25 April, 2012;
originally announced April 2012.
-
Dynamical Rare event simulation techniques for equilibrium and non-equilibrium systems
Authors:
Titus S. van Erp
Abstract:
I give an overview of rare event simulation techniques to generate dynamical pathways across high free energy barriers. The methods on which I will concentrate are the reactive flux approach, transition path sampling, (replica-exchange) transition interface sampling, partial path sampling/milestoning, and forward flux sampling. These methods have in common that they aim to simulate true molecular…
▽ More
I give an overview of rare event simulation techniques to generate dynamical pathways across high free energy barriers. The methods on which I will concentrate are the reactive flux approach, transition path sampling, (replica-exchange) transition interface sampling, partial path sampling/milestoning, and forward flux sampling. These methods have in common that they aim to simulate true molecular dynamics trajectories at a much faster rate than naive brute force molecular dynamics. The advantages and disadvantages of these methods are discussed and compared for a simple one-dimensional test system. These numerical results reveal some important pitfalls of the present non-equilibrium methods that have no easy solution and show that caution is necessary when interpreting their results.
△ Less
Submitted 31 March, 2011; v1 submitted 5 January, 2011;
originally announced January 2011.
-
Kinetics of intermediate-mediated self-assembly in nano-sized materials: a generic model
Authors:
James F. Lutsko,
Vasileios Basios,
Gregoire Nicolis,
Tom P. Caremans,
Alexander Aerts,
Johan A. Martens,
Christine E. A. Kirschhock,
Titus S. van Erp
Abstract:
We propose in this paper a generic model of a non-standard aggregation mechanism for self-assembly processes of a class of materials involving the mediation of intermediates consisting of a polydisperse population of nano-sized particles. The model accounts for a long induction period in the process. The proposed mechanism also gives insight on future experiments aiming at a more comprehensive pic…
▽ More
We propose in this paper a generic model of a non-standard aggregation mechanism for self-assembly processes of a class of materials involving the mediation of intermediates consisting of a polydisperse population of nano-sized particles. The model accounts for a long induction period in the process. The proposed mechanism also gives insight on future experiments aiming at a more comprehensive picture of the role of self-organization in self-assembly processes.
△ Less
Submitted 29 March, 2010;
originally announced March 2010.
-
Comment on "A generalized Langevin formalism of complete DNA melting transition"
Authors:
Titus S. van Erp,
Santiago Cuesta-Lopez,
Johannes-Geert Hagmann,
Michel Peyrard
Abstract:
We show that the calculated DNA denaturation curves for finite (Peyrard-Bishop-Dauxois (PBD) chains are intrinsically undefined.
We show that the calculated DNA denaturation curves for finite (Peyrard-Bishop-Dauxois (PBD) chains are intrinsically undefined.
△ Less
Submitted 26 February, 2009;
originally announced February 2009.
-
Efficient path sampling on multiple reaction channels
Authors:
Titus S. van Erp
Abstract:
Due to the time scale problem, rare events are not accessible by straight forward molecular dynamics. The presence of multiple reaction channels complicates the problem even further. The feasibility of the standard free energy based methods relies strongly on the success in finding a proper reaction coordinate. This can be very difficult task in high-dimensional complex systems and even more if…
▽ More
Due to the time scale problem, rare events are not accessible by straight forward molecular dynamics. The presence of multiple reaction channels complicates the problem even further. The feasibility of the standard free energy based methods relies strongly on the success in finding a proper reaction coordinate. This can be very difficult task in high-dimensional complex systems and even more if several distinct reaction channels exist. Moreover, even if a proper reaction coordinate can be found, ergodic sampling will be a challenge. In this article, we discuss the recent advancements of path sampling methods to tackle this problem. We argue why the path sampling methods, via the transition interface sampling technique, is less sensitive to the choice of reaction coordinate. Moreover, we review a new algorithm, parallel path swapping, that can dramatically improve the ergodic sampling of trajectories for the multiple reaction channel systems.
△ Less
Submitted 21 November, 2007;
originally announced November 2007.
-
Reaction rate calculation by parallel path swapping
Authors:
Titus S. van Erp
Abstract:
The efficiency of path sampling simulations can be improved considerably using the approach of path swapping. For this purpose, we have devised a new algorithmic procedure based on the transition interface sampling technique. In the same spirit of parallel tempering, paths between different ensembles are swapped, but the role of temperature is here played by the interface position. We have teste…
▽ More
The efficiency of path sampling simulations can be improved considerably using the approach of path swapping. For this purpose, we have devised a new algorithmic procedure based on the transition interface sampling technique. In the same spirit of parallel tempering, paths between different ensembles are swapped, but the role of temperature is here played by the interface position. We have tested the method on the denaturation transition of DNA using the Peyrard-Bishop-Dauxois model. We find that the new algorithm gives a reduction of the computational cost by a factor 20.
△ Less
Submitted 25 April, 2007;
originally announced April 2007.
-
Prospects of Transition Interface Sampling simulations for the theoretical study of zeolite synthesis
Authors:
Titus S. van Erp,
Tom P. Caremans,
Christine E. A. Kirschhock,
Johan A. Martens
Abstract:
The transition interface sampling (TIS) technique allows to overcome large free energy barriers within reasonable simulation time, which is impossible for straightforward molecular dynamics. Still, the method does not impose an artificial driving force, but it surmounts the timescale problem by an importance sampling of true dynamical pathways. Recently, it was shown that the efficiency of TIS t…
▽ More
The transition interface sampling (TIS) technique allows to overcome large free energy barriers within reasonable simulation time, which is impossible for straightforward molecular dynamics. Still, the method does not impose an artificial driving force, but it surmounts the timescale problem by an importance sampling of true dynamical pathways. Recently, it was shown that the efficiency of TIS to calculate reaction rates is less sensitive to the choice of reaction coordinate than those of the standard free energy based techniques. This could be an important advantage in complex systems for which a good reaction coordinate is usually very difficult to find. We explain the principles of this method and discuss some of the promising applications related to zeolite formation.
△ Less
Submitted 15 January, 2007;
originally announced January 2007.
-
Efficiency analysis of reaction rate calculation methods using analytical models I: The 2D sharp barrier
Authors:
Titus S. van Erp
Abstract:
We analyze the efficiency of different methods for the calculation of reaction rates in the case of two simple analytical benchmark systems. Two classes of methods are considered: the first are based on the free energy calculation along a reaction coordinate and the calculation of the transmission coefficient, the second on the sampling of dynamical pathways. We give scaling rules for how this e…
▽ More
We analyze the efficiency of different methods for the calculation of reaction rates in the case of two simple analytical benchmark systems. Two classes of methods are considered: the first are based on the free energy calculation along a reaction coordinate and the calculation of the transmission coefficient, the second on the sampling of dynamical pathways. We give scaling rules for how this efficiency depends on barrier height and width, and we hand out simple optimization rules for the method-specific parameters. We show that the path sampling methods, using the transition interface sampling technique, become exceedingly more efficient than the others when the reaction coordinate is not the optimal one.
△ Less
Submitted 8 August, 2006; v1 submitted 10 May, 2006;
originally announced May 2006.
-
Breakdown of Lindstedt Expansion for Chaotic Maps
Authors:
Guido Gentile,
Titus S. van Erp
Abstract:
In a previous paper of one of us [Europhys. Lett. 59 (2002), 330--336] the validity of Greene's method for determining the critical constant of the standard map (SM) was questioned on the basis of some numerical findings. Here we come back to that analysis and we provide an interpretation of the numerical results by showing that no contradiction is found with respect to Greene's method. We show…
▽ More
In a previous paper of one of us [Europhys. Lett. 59 (2002), 330--336] the validity of Greene's method for determining the critical constant of the standard map (SM) was questioned on the basis of some numerical findings. Here we come back to that analysis and we provide an interpretation of the numerical results by showing that no contradiction is found with respect to Greene's method. We show that the previous results based on the expansion in Lindstedt series do correspond to the transition value but for a different map: the semi-standard map (SSM). Moreover, we study the expansion obtained from the SM and SSM by suppressing the small divisors. The first case turns out to be related to Kepler's equation after a proper transformation of variables. In both cases we give an analytical solution for the radius of convergence, that represents the singularity in the complex plane closest to the origin. Also here, the radius of convergence of the SM's analogue turns out to be lower than the one of the SSM. However, despite the absence of small denominators these two radii are lower than the ones of the true maps for golden mean winding numbers. Finally, the analyticity domain and, in particular, the critical constant for the two maps without small divisors are studied analytically and numerically. The analyticity domain appears to be an perfect circle for the SSM analogue, while it is stretched along the real axis for the SM analogue yielding a critical constant that is larger than its radius of convergence.
△ Less
Submitted 30 July, 2004;
originally announced July 2004.
-
Elaborating Transition Interface Sampling Methods
Authors:
Titus S. van Erp,
Peter G. Bolhuis
Abstract:
We review two recently developed efficient methods for calculating rate constants of processes dominated by rare events in high-dimensional complex systems. The first is transition interface sampling (TIS), based on the measurement of effective fluxes through hypersurfaces in phase space. TIS improves efficiency with respect to standard transition path sampling (TPS) rate constant techniques, be…
▽ More
We review two recently developed efficient methods for calculating rate constants of processes dominated by rare events in high-dimensional complex systems. The first is transition interface sampling (TIS), based on the measurement of effective fluxes through hypersurfaces in phase space. TIS improves efficiency with respect to standard transition path sampling (TPS) rate constant techniques, because it allows a variable path length and is less sensitive to recrossings. The second method is the partial path version of TIS. Developed for diffusive processes, it exploits the loss of long time correlation. We discuss the relation between the new techniques and the standard reactive flux methods in detail. Path sampling algorithms can suffer from ergodicity problems, and we introduce several new techniques to alleviate these problems, notably path swapping, stochastic configurational bias Monte Carlo shooting moves and order-parameter free path sampling. In addition, we give algorithms to calculate other interesting properties from path ensembles besides rate constants, such as activation energies and reaction mechanisms.
△ Less
Submitted 6 May, 2004;
originally announced May 2004.
-
Investigating Rare Events by Transition Interface Sampling
Authors:
Daniele Moroni,
Titus S. van Erp,
Peter G. Bolhuis
Abstract:
We briefly review simulation schemes for the investigation of rare transitions and we resume the recently introduced Transition Interface Sampling, a method in which the computation of rate constants is recast into the computation of fluxes through interfaces dividing the reactant and product state.
We briefly review simulation schemes for the investigation of rare transitions and we resume the recently introduced Transition Interface Sampling, a method in which the computation of rate constants is recast into the computation of fluxes through interfaces dividing the reactant and product state.
△ Less
Submitted 25 November, 2003;
originally announced November 2003.
-
Rate constants for diffusive processes by partial path sampling
Authors:
Daniele Moroni,
Peter G. Bolhuis,
Titus S. van Erp
Abstract:
We introduce a path sampling method for the computation of rate constants for systems with a highly diffusive character. Based on the recently developed algorithm of transition interface sampling (TIS) this procedure increases the efficiency by sampling only parts of complete transition trajectories confined within a certain region. The algorithm assumes the loss of memory for highly diffusive p…
▽ More
We introduce a path sampling method for the computation of rate constants for systems with a highly diffusive character. Based on the recently developed algorithm of transition interface sampling (TIS) this procedure increases the efficiency by sampling only parts of complete transition trajectories confined within a certain region. The algorithm assumes the loss of memory for highly diffusive progression along the reaction coordinate. We compare the new technique to the TIS method for a simple diatomic system and show that the computation time of the new method scales linearly, instead of quadraticaly, with the length of the diffusive barrier. The validity of the memory loss assumption is also discussed.
△ Less
Submitted 20 October, 2003;
originally announced October 2003.
-
A Novel Path Sampling Method for the Calculation of Rate Constants
Authors:
Titus S. van Erp,
Daniele Moroni,
Peter G. Bolhuis
Abstract:
We derive a novel efficient scheme to measure the rate constant of transitions between stable states separated by high free energy barriers in a complex environment within the framework of transition path sampling. The method is based on directly and simultaneously measuring the fluxes through many phase space interfaces and increases the efficiency with at least a factor of two with respect to…
▽ More
We derive a novel efficient scheme to measure the rate constant of transitions between stable states separated by high free energy barriers in a complex environment within the framework of transition path sampling. The method is based on directly and simultaneously measuring the fluxes through many phase space interfaces and increases the efficiency with at least a factor of two with respect to existing transition path sampling rate constant algorithms. The new algorithm is illustrated on the isomerization of a diatomic molecule immersed in a simple fluid.
△ Less
Submitted 28 October, 2002;
originally announced October 2002.
-
Aubry transition studied by direct evaluation of the modulation functions of infinite incommensurate systems
Authors:
Titus S. van Erp,
Annalisa Fasolino
Abstract:
Incommensurate structures can be described by the Frenkel Kontorova model. Aubry has shown that, at a critical value K_c of the coupling of the harmonic chain to an incommensurate periodic potential, the system displays the analyticity breaking transition between a sliding and pinned state. The ground state equations coincide with the standard map in non-linear dynamics, with smooth or chaotic o…
▽ More
Incommensurate structures can be described by the Frenkel Kontorova model. Aubry has shown that, at a critical value K_c of the coupling of the harmonic chain to an incommensurate periodic potential, the system displays the analyticity breaking transition between a sliding and pinned state. The ground state equations coincide with the standard map in non-linear dynamics, with smooth or chaotic orbits below and above K_c respectively. For the standard map, Greene and MacKay have calculated the value K_c=.971635. Conversely, evaluations based on the analyticity breaking of the modulation function have been performed for high commensurate approximants. Here we show how the modulation function of the infinite system can be calculated without using approximants but by Taylor expansions of increasing order. This approach leads to a value K_c'=.97978, implying the existence of a golden invariant circle up to K_c' > K_c.
△ Less
Submitted 7 May, 2002; v1 submitted 16 January, 2002;
originally announced January 2002.