Adaptive importance sampling Monte Carlo simulation of rare transition events.

We develop a general theoretical framework for the recently proposed importance sampling method for enhancing the efficiency of rare-event simulations [W. Cai, M. H. Kalos, M. de Koning, and V. V. Bulatov, Phys. Rev. E 66, 046703 (2002)], and discuss practical aspects of its application. We define the success/fail ensemble of all possible successful and failed transition paths of any duration and demonstrate that in this formulation the rare-event problem can be interpreted as a "hit-or-miss" Monte Carlo quadrature calculation of a path integral. The fact that the integrand contributes significantly only for a very tiny fraction of all possible paths then naturally leads to a "standard" importance sampling approach to Monte Carlo (MC) quadrature and the existence of an optimal importance function. In addition to showing that the approach is general and expected to be applicable beyond the realm of Markovian path simulations, for which the method was originally proposed, the formulation reveals a conceptual analogy with the variational MC (VMC) method. The search for the optimal importance function in the former is analogous to finding the ground-state wave function in the latter. In two model problems we discuss practical aspects of finding a suitable approximation for the optimal importance function. For this purpose we follow the strategy that is typically adopted in VMC calculations: the selection of a trial functional form for the optimal importance function, followed by the optimization of its adjustable parameters. The latter is accomplished by means of an adaptive optimization procedure based on a combination of steepest-descent and genetic algorithms.

[1]  Huajian Gao,et al.  Simulating materials failure by using up to one billion atoms and the world's fastest computer: Work-hardening , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Gerard T. Barkema,et al.  Monte Carlo Methods in Statistical Physics , 1999 .

[3]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[4]  Christoph Dellago,et al.  On the calculation of reaction rate constants in the transition path ensemble , 1999 .

[5]  Christoph Dellago,et al.  Sampling ensembles of deterministic transition pathways , 1998 .

[6]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[7]  Paul R. C. Kent,et al.  Monte Carlo energy and variance-minimization techniques for optimizing many-body wave functions , 1999 .

[8]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[9]  A. Messiah.  Quantum Mechanics , 1961 .

[10]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[11]  Manuel Laso,et al.  Conformational kinetics in liquid n-butane by transition path sampling , 2001 .

[12]  A. Voter.  A method for accelerating the molecular dynamics simulation of infrequent events , 1997 .

[13]  H. Jónsson,et al.  Theoretical studies of atomic-scale processes relevant to crystal growth. , 2000, Annual review of physical chemistry.

[14]  Maurice de Koning,et al.  Importance sampling of rare transition events in Markov processes. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[16]  A. Voter.  Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

[17]  R. Needs,et al.  Quantum Monte Carlo simulations of solids , 2001 .

[18]  Heiko Rieger,et al.  Optimization algorithms in physics , 2001 .

[19]  Bulatov,et al.  Saddle-point configurations for self-interstitial migration in silicon. , 1996, Physical review. B, Condensed matter.

[20]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .