Sampling rare events in nonequilibrium and nonstationary systems
Abstract
Although many computational methods for rare event sampling exist, this type of calculation is not usually practical for general nonequilibrium conditions, with macroscopically irreversible dynamics and away from both stationary and metastable states. A novel method for calculating the time-series of the probability of a rare event is presented which is designed for these conditions. The method is validated for the cases of the Glauber-Ising model under time-varying shear flow, the Kawasaki-Ising model after a quench into the region between nucleation dominated and spinodal decomposition dominated phase change dynamics, and the parallel open asymmetric exclusion process. The method requires a subdivision of the phase space of the system: it is benchmarked and found to scale well for increasingly fine subdivisions, meaning that it can be applied without detailed foreknowledge of the physically important reaction pathways.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- December 2010
- DOI:
- 10.1063/1.3525099
- arXiv:
- arXiv:1001.2456
- Bibcode:
- 2010JChPh.133x4101B
- Keywords:
-
- computational fluid dynamics;
- drops;
- Ising model;
- metastable states;
- nonequilibrium flow;
- nucleation;
- probability;
- sampling methods;
- shear flow;
- spinodal decomposition;
- time series;
- two-phase flow;
- 47.55.db;
- 47.11.-j;
- 05.50.+q;
- 02.50.Cw;
- 47.70.Nd;
- Drop and bubble formation;
- Computational methods in fluid dynamics;
- Lattice theory and statistics;
- Probability theory;
- Nonequilibrium gas dynamics;
- Condensed Matter - Mesoscale and Nanoscale Physics;
- Condensed Matter - Statistical Mechanics;
- Physics - Atomic and Molecular Clusters;
- Physics - Computational Physics
- E-Print:
- 11 pages, 6 figures