Kinetic simulations stochastic methods

The relative fluctuations in Monte Carlo simulations are of the order of magnitude where N is the total number of molecules in the simulation. The observed error in kinetic simulations is about 1-2% when lO molecules are used. In the computer calculations described by Schaad, the grids of the technique shown here are replaced by computer memory, so the capacity of the memory is one limit on the maximum number of molecules. Other programs for stochastic simulation make use of different routes of calculation, and the number of molecules is not a limitation. Enzyme kinetics and very complex oscillatory reactions have been modeled. These simulations are valuable for establishing whether a postulated kinetic scheme is reasonable, for examining the appearance of extrema or induction periods, applicability of the steady-state approximation, and so on. Even the manual method is useful for such purposes. 

More detailed studies of eleetroeatalytie processes, which incorporate heterogeneous surfaee geometries and finite surface mobilities of reactants, require kinetic Monte Carlo simulations. This stochastic method has been successfully applied in the field of heterogeneous catalysis on nanosized catalyst particles [59,60]. Since these simulations permit atomistic resolution, any level of structural detail may easily be incorporated. Moreover, kinetic Monte Carlo simulations proceed in real time. The simulation of current transients or cyclic voltammograms is, thus, straightforward [61]. 

Stochastic methods simulate the dynamic changes that occur in the structure of the adlayer of catalytic surface and thus model the elementary surface kinetics [82-ioo] j jjg temporal changes of a system can be followed by solving the stochastic master equation which simulates the dynamic changes in the system as it moves from one state i) to another state (j). The master equation, which can written as... 

Numerical integration (sometimes referred to as solving or simulation) of differential equations, ordinary or partial, involves using a computer to obtain an approximate and discrete (in time and/or space) solution. In chemical kinetics, these differential equations are typically the rate laws that describe the time evolution of the system. One obtains results for the mean concentrations, without any information about the (typically very small) fluctuations that are inevitably present. Continuation and sensitivity analysis techniques enable one to extrapolate from a numerically obtained solution at one set of parameters (e.g., rate constants or initial concentrations) to the behavior of the system at other parameter values, without having to carry out a full numerical integration each time the parameters are changed. Other approaches, sometimes referred to collectively as stochastic methods (Gardiner, 1990), can provide data about fluctuations, but these require considerably more computational labor and are often impractical for models that include more than a few variables. 

In the molecular dynamics approach, one represents the system by a set of particles randomly distributed in space and then solves Newton s equations for the motion of each particle in the system. Particles react according to the rules assumed for the kinetics when they approach each other sufficiently closely. While the particle motion is deterministic, a stochastic aspect is present in the choice of the initial particle positions and velocities. Again, an average over replicate runs will reduce the statistical fluctuations. By necessity, only a very small system over a very short time (picoseconds) can be simulated (Kawezynski and Gorecki, 1992). We will not treat stochastic methods in any further detail here, but, instead. 

Different stochastic methods have been considered to overcome the accretion limit problem. The main focus has been on the master equation method [47 9] and macroscopic Monte Carlo simulations [49-51]. Of the two, the master equation method can be more easily coupled to rate equations, which handle the gas phase chemistry in the most straightforward way, since it treats the gas and surface kinetics with equations of similar form. 

Reaction kinetic models can be simulated not only by solving the kinetic system of differential equations but also via simulating the equivalent stochastic models. Computer codes are available that solve the stochastic kinetic equations. One of these is the Chemical Kinetics Simulator (CKS) program that was developed at IBM s Almaden Research Centre. It provides a rapid, interactive method for the accurate simulation of chemical reactions. CKS is a good tool for teaching the principles of stochastic reaction kinetics to students and trainees. 

Although the collision and transition state theories represent two important methods of attacking the theoretical calculation of reaction rates, they are not the only approaches available. Alternative methods include theories based on nonequilibrium statistical mechanics, stochastic theories, and Monte Carlo simulations of chemical dynamics. Consult the texts by Johnson (62), Laidler (60), and Benson (59) and the review by Wayne (63) for a further introduction to the theoretical aspects of reaction kinetics. 

The dWi are Gaussian white noise processes, and their strength a is related to the kinetic friction y through the fluctuation-dissipation relation.72 When deriving integrators for these methods, one has to be careful to take into account the special character of the random forces employed in these simulations.73 A variant of the velocity Verlet method, including a stochastic dynamics treatment of constraints, can be found in Ref. 74. The stochastic... 

A number of different techniques have been developed for studying nonhomogeneous radiolysis kinetics, and they can be broken down into two groups, deterministic and stochastic. The former used conventional macroscopic treatments of concentration, diffusion, and reaction to describe the chemistry of a typical cluster or track of reactants. In contrast, the latter approach considers the chemistry of simulated tracks of realistic clusters using probabilistic methods to model the kinetics. Each treatment has advantages and limitations, and at present, both treatments have a valuable role to play in modeling radiation chemistry. 

The irradiation of water is immediately followed by a period of fast chemistry, whose short-time kinetics reflects the competition between the relaxation of the nonhomogeneous spatial distributions of the radiation-induced reactants and their reactions. A variety of gamma and energetic electron experiments are available in the literature. Stochastic simulation methods have been used to model the observed short-time radiation chemical kinetics of water and the radiation chemistry of aqueous solutions of scavengers for the hydrated electron and the hydroxyl radical to provide fundamental information for use in the elucidation of more complex, complicated chemical, and biological systems found in real-world scenarios. 

Ray Kapral came to Toronto from the United States in 1969. His research interests center on theories of rate processes both in systems close to equilibrium, where the goal is the development of a microscopic theory of condensed phase reaction rates,89 and in systems far from chemical equilibrium, where descriptions of the complex spatial and temporal reactive dynamics that these systems exhibit have been developed.90 He and his collaborators have carried out research on the dynamics of phase transitions and critical phenomena, the dynamics of colloidal suspensions, the kinetic theory of chemical reactions in liquids, nonequilibrium statistical mechanics of liquids and mode coupling theory, mechanisms for the onset of chaos in nonlinear dynamical systems, the stochastic theory of chemical rate processes, studies of pattern formation in chemically reacting systems, and the development of molecular dynamics simulation methods for activated chemical rate processes. His recent research activities center on the theory of quantum and classical rate processes in the condensed phase91 and in clusters, and studies of chemical waves and patterns in reacting systems at both the macroscopic and mesoscopic levels. 

After the kinetic model for the network is defined, a simulation method needs to be chosen, given the systemic phenomenon of interest. The phenomenon might be spatial. Then it has to be decided whether in addition stochasticity plays a role or not. In the former case the kinetic model should be described with a reaction-diffusion master equation [81], whereas in the latter case partial differential equations should suffice. If the phenomenon does not involve a spatial organization, the dynamics can be simulated either using ordinary differential equations [47] or master equations [82-84]. In the latter case but not in the former, stochasticity is considered of importance. A first-order estimate of the magnitude of stochastic fluctuations can be obtained using the linear noise approximation, given only the ordinary differential equation description of the kinetic model [83-85, 87]. 

It is often stated that MC methods lack real time and results are usually reported in MC events or steps. While this is immaterial as far as equilibrium is concerned, following real dynamics is essential for comparison to solutions of partial differential equations and/or experimental data. It turns out that MC simulations follow the stochastic dynamics of a master equation, and with appropriate parameterization of the transition probabilities per unit time, they provide continuous time information as well. For example, Gillespie has laid down the time foundations of MC for chemical reactions in a spatially homogeneous system.f His approach is easily extendable to arbitrarily complex computational systems when individual events have a prescribed transition probability per unit time, and is often referred to as the kinetic Monte Carlo or dynamic Monte Carlo (DMC) method. The microscopic processes along with their corresponding transition probabilities per unit time can be obtained via either experiments such as field emission or fast scanning tunneling microscopy or shorter time scale DFT/MD simulations discussed earlier. The creation of a database/lookup table of transition... 

---