Simulation stochastic, Monte Carlo

Having specified the interactions (i.e., the model of the system), the actual simulation then constructs a sequence of states (or the system trajectory) in some statistical mechanical ensemble. Simulations can be stochastic (Monte Carlo (MC)) or deterministic (MD), or they can combine elements of both, such as force-biased MC, Brownian dynamics, or generalized Lan-gevin dynamics. It is usually assumed that the laws of classical mechanics (i.e., Newton s second law) may adequately describe the atoms and molecules in the physical system. 

Finally, "data" can be obtained from computer simulations (26), whether deterministic (molecular dynamics) or stochastic (Monte Carlo). This approach provides a level of microscopic detail not available with any of the above experimental techniques. Results from computer simulations, furthermore, can be both qualitative (for example, observation of cavity dynamics in repulsive supercritical systems (12)) as well as quantitative. However, because true intermolecular potentials are not known exactly, simulation results must be interpreted with caution, especially if they are used to study the behavior of real systems. Through simulations, therefore, one obtains exact answers to ideal (as opposed to real) problems. 

Over the past ten years the numerical simulation of the behavior of complex reaction systems has become a fairly routine procedure, and has been widely used in many areas of chemistry, [l] The most intensive application has been in environmental, atmospheric, and combustion science, where mechanisms often consisting of several hundred reactions are involved. Both deterministic (numerical solution of mass-action differential equations) and stochastic (Monte-Carlo) methods have been used. The former approach is by far the most popular, having been made possible by the development of efficient algorithms for the solution of the "stiff" ODE problem. Edelson has briefly reviewed these developments in a symposium volume which includes several papers on the mathematical techniques and their application. [2]... 

Not having enough space to summarize all of the recent work, we will give a couple of examples. A welcome new development is the integration of energy-loss calculations [79,87] with the stochastic Monte Carlo simulations of... 

The basic idea behind an atomistic-level simulation is quite simple. Given an accurate description of the energetic interactions between a collection of atoms and a set of initial atomic coordinates (and in some cases, velocities), the positions (velocities) of these atoms are advanced subject to a set of thermodynamic constraints. If the positions are advanced stochastically, we call the simulation method Monte Carlo or MG [10]. No velocities are required for this technique. If the positions and velocities are advanced deterministically, we call the method molecular dynamics or MD [10]. Other methods exist which are part stochastic and part deterministic, but we need not concern ourselves with these details here. The important point is that statistical mechanics teUs us that the collection of atomic positions that are obtained from such a simulation, subject to certain conditions, is enough to enable aU of the thermophysical properties of the system to be determined. If the velocities are also available (as in an MD simulation), then time-dependent properties may also be computed. If done properly, the numerical method that generates the trajectories... 

Stochastic (Monte Carlo) methods solve for the neutron flux by simulating particle transport rather than by numerically solving the Boltzmann transport equation (although the simulational algorithm can be derived from the equation itself as well as from the physics). 

An analytical solution to the master equation is only possible for very simple systems. The master equation, however, can readily be simulated by using stochastic kinetics or more specifically kinetic Monte Carlo simulation. Several Monte Carlo algorithms exist. More details on kinetic Monte Carlo simulation can be found in the Appendix. 

Makeev, A., D. Maroudas, and LG. Kevrekidis. 2002a. Coarse stability and bifurcation analysis using stochastic simulators kinetic Monte Carlo examples. Journal of Chemical Physics 116(23) 10083-10091. 

Other methods which are applied to conformational analysis and to generating multiple conformations and which can be regarded as random or stochastic techniques, since they explore the conformational space in a non-deterministic fashion, arc genetic algorithms (GA) [137, 1381 simulation methods, such as molecular dynamics (MD) and Monte Carlo (MC) simulations 1139], as well as simulated annealing [140], All of those approaches and their application to generate ensembles of conformations arc discussed in Chapter II, Section 7.2 in the Handbook. 

Monte Carlo search methods are stochastic techniques based on the use of random numbers and probability statistics to sample conformational space. The name Monte Carlo was originally coined by Metropolis and Ulam [4] during the Manhattan Project of World War II because of the similarity of this simulation technique to games of chance. Today a variety of Monte Carlo (MC) simulation methods are routinely used in diverse fields such as atmospheric studies, nuclear physics, traffic flow, and, of course, biochemistry and biophysics. In this section we focus on the application of the Monte Carlo method for... 

By far the most common methods of studying aqueous interfaces by simulations are the Metropolis Monte Carlo (MC) technique and the classical molecular dynamics (MD) techniques. They will not be described here in detail, because several excellent textbooks and proceedings volumes (e.g., [2-8]) on the subject are available. In brief, the stochastic MC technique generates microscopic configurations of the system in the canonical (NYT) ensemble the deterministic MD method solves Newton s equations of motion and generates a time-correlated sequence of configurations in the microcanonical (NVE) ensemble. Structural and thermodynamic properties are accessible by both methods the MD method provides additional information about the microscopic dynamics of the system. 

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. 

They point out that at the heart of technical simulation there must be unreality otherwise, there would not be need for simulation. The essence of the subject linder study may be represented by a model of it that serves a certain purpose, e.g., the use of a wind tunnel to simulate conditions to which an aircraft may be subjected. One uses the Monte Carlo method to study an artificial stochastic model of a physical or mathematical process, e.g., evaluating a definite integral by probability methods (using random numbers) using the graph of the function as an aid. 

Special considerations are required in estimating paraimeters from experimental measurements when the relationship between output responses, input variables and paraimeters is given by a Monte Carlo simulation. These considerations, discussed in our first paper 1), relate to the stochastic nature of the solution and to the fact that the Monte Carlo solution is numerical rather than functional. The motivation for using Monte Carlo methods to model polymer systems stems from the fact that often the solution... 

The method for estimating parameters from Monte Carlo simulation, described in mathematical detail by Reilly and Duever (in preparation), uses a Bayesian approach to establish the posterior distribution for the parameters based on a Monte Carlo model. The numerical nature of the solution requires that the posterior distribution be handled in discretised form as an array in computer storage using the method of Reilly 2). The stochastic nature of Monte Carlo methods implies that output responses are predicted by the model with some amount of uncertainty for which the term "shimmer" as suggested by Andres (D.B. Chambers, SENES Consultants Limited, personal communication, 1985) has been adopted. The model for the uth of n experiments can be expressed by... 

We have presented applications of a parameter estimation technique based on Monte Carlo simulation to problems in polymer science involving sequence distribution data. In comparison to approaches involving analytic functions, Monte Carlo simulation often leads to a simpler solution of a model particularly when the process being modelled involves a prominent stochastic coit onent. 

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. 

In Sect. 7.4.6, we discussed various stochastic simulation techniques that include the kinetics of recombination and free-ion yield in multiple ion-pair spurs. No further details will be presented here, but the results will be compared with available experiments. In so doing, we should remember that in the more comprehensive Monte Carlo simulations of Bartczak and Hummel (1986,1987, 1993,1997) Hummel and Bartczak, (1988) the recombination reaction is taken to be fully diffusion-controlled and that the diffusive free path distribution is frequently assumed to be rectangular, consistent with the diffusion coefficient, instead of a more realistic distribution. While the latter assumption can be justified on the basis of the central limit theorem, which guarantees a gaussian distribution for a large number of scatterings, the first assumption is only valid for low-mobility liquids. 

The WATS model is formulated in deterministic terms. However, an extension to include simple Monte-Carlo stochastic simulation is possible, taking into consideration a measured variability of the process parameters. 

This criterion requires a search through a nonconvex multidimensional conformation space that contains an immense number of minima. Optimization techniques that have been applied to the problem include Monte Carlo methods, simulated annealing, genetic methods, and stochastic search, among others. For reviews of the application of various optimization methods refer to Pardalos et al. (1996), Vasquez et al. (1994), or Schlick et al. (1999). 

---