Reactive Monte Carlo

Reactive MC is a technique that permits the incorporation of chemical equilibrium into MC simulations (Johnson et al. 1994 Smith and Triska 1994). It can be combined with GEMC to evaluate the 

---

Shocked Materials Using the Reactive Monte Carlo Method. 

The Langmuir-Hinshelwood picture is essentially that of Fig. XVIII-14. If the process is unimolecular, the species meanders around on the surface until it receives the activation energy to go over to product(s), which then desorb. If the process is bimolecular, two species diffuse around until a reactive encounter occurs. The reaction will be diffusion controlled if it occurs on every encounter (see Ref. 211) the theory of surface diffusional encounters has been treated (see Ref. 212) the subject may also be approached by means of Monte Carlo/molecular dynamics techniques [213]. In the case of activated bimolecular reactions, however, there will in general be many encounters before the reactive one, and the rate law for the surface reaction is generally written by analogy to the mass action law for solutions. That is, for a bimolecular process, the rate is taken to be proportional to the product of the two surface concentrations. It is interesting, however, that essentially the same rate law is obtained if the adsorption is strictly localized and species react only if they happen to adsorb on adjacent sites (note Ref. 214). (The apparent rate law, that is, the rate law in terms of gas pressures, depends on the form of the adsorption isotherm, as discussed in the next section.)... 

To conclude, the introduction of species-selective membranes into the simulation box results in the osmotic equilibrium between a part of the system containing the products of association and a part in which only a one-component Lennard-Jones fluid is present. The density of the fluid in the nonreactive part of the system is lower than in the reactive part, at osmotic equilibrium. This makes the calculations of the chemical potential efficient. The quahty of the results is similar to those from the grand canonical Monte Carlo simulation. The method is neither restricted to dimerization nor to spherically symmetric associative interactions. Even in the presence of higher-order complexes in large amounts, the proposed approach remains successful. 

FIG. 7 Log-log plots of the interface width (w ) versus the Monte Carlo time t, measured at different adsorption probabihties using channels of width L = 30. Data were obtained during the displacement of an A-poisoned phase by the reactive regime. From top to bottom the probabihties are 0.5192, 0.5202, 0.5211, 0.5215, and 0.5238. 

Another interesting version of the MM model considers a variable excluded-volume interaction between same species particles [92]. In the absence of interactions the system is mapped on the standard MM model which has a first-order IPT between A- and B-saturated phases. On increasing the strength of the interaction the first-order transition line, observed for weak interactions, terminates at a tricritical point where two second-order transitions meet. These transitions, which separate the A-saturated, reactive, and B-saturated phases, belong to the same universality class as directed percolation, as follows from the value of critical exponents calculated by means of time-dependent Monte Carlo simulations and series expansions [92]. 

A general method has been developed for the estimation of model parameters from experimental observations when the model relating the parameters and input variables to the output responses is a Monte Carlo simulation. The method provides point estimates as well as joint probability regions of the parameters. In comparison to methods based on analytical models, this approach can prove to be more flexible and gives the investigator a more quantitative insight into the effects of parameter values on the model. The parameter estimation technique has been applied to three examples in polymer science, all of which concern sequence distributions in polymer chains. The first is the estimation of binary reactivity ratios for the terminal or Mayo-Lewis copolymerization model from both composition and sequence distribution data. Next a procedure for discriminating between the penultimate and the terminal copolymerization models on the basis of sequence distribution data is described. Finally, the estimation of a parameter required to model the epimerization of isotactic polystyrene is discussed. 

The model in either its analytic or Monte Carlo form is dependent upon four reactivity ratios defined by ... 

Figure 5.2 Oxygen states present at the ends of -Cu-O-Cu-O- chains are established as the active sites in ammonia oxidation at Cu(110) from a Monte Carlo simulation of the growth of the oxygen adlayer. The reactivity (the experimental curve) is best fitted to the atoms present at chain ends. (Reproduced from Ref. 7).

In a few instances, quantum mechanical calculations on the stability and reactivity of adsorbates have been combined with Monte Carlo simulations of dynamic or kinetic processes. In one example, both the ordering of NO on Rh(lll) during adsorption and its TPD under UHV conditions were reproduced using a dynamic Monte Carlo model involving lateral interactions derived from DFT calculations and different adsorption... 

The independent reaction time (1RT) model was introduced as a shortcut Monte Carlo simulation of pairwise reaction times without explicit reference to diffusive trajectories (Clifford et al, 1982b). At first, the initial positions of the reactive species (any number and kind) are simulated by convolving from a given (usually gaussian) distribution using random numbers. These are examined for immediate reaction—that is, whether any interparticle separation is within the respective reaction radius. If so, such particles are removed and the reactions are recorded as static reactions. 

We have already commented that the master equation method is not suitable, at present, to handle reactive products because, inter alia, the dimensionality of the problem increases with reactions of products. There is no difficulty, in principle, to including reactive products in Monte Carlo simulation, since the time of reaction and the positions of the products can be recorded. In practice, however, this requires a greatly expanded computational effort, which is discouraging. 

It is clear that the different procedures of handling reactive products are based on different approximations therefore, somewhat different results are expected. On the whole, since the IRT methodology is based on the conceived independence of pairwise bimolecular reactions, it needs validation by comparison with well-known examples for which Monte Carlo results are available. Such validations have in fact been made (see Figures 4 and 7 in Clifford et ah, 1986, and Figures 12, 27, and 28 in Pimblott and Green, 1995). 

Sometimes the theoretical or computational approach to description of molecular structure, properties, and reactivity cannot be based on deterministic equations that can be solved by analytical or computational methods. The properties of a molecule or assembly of molecules may be known or describable only in a statistical sense. Molecules and assemblies of molecules exist in distributions of configuration, composition, momentum, and energy. Sometimes, this statistical character is best captured and studied by computer experiments molecular dynamics, Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. Interaction potentials based on quantum mechanics, classical particle mechanics, continuum mechanics, or empiricism are specified and the evolution of the system is then followed in time by simulation of motions resulting from these direct... 

GH Theory was originally developed to describe chemical reactions in solution involving a classical nuclear solute reactive coordinate x. The identity of x will depend of course on the reaction type, i.e., it will be a separation coordinate in an SnI unimolecular ionization and an asymmetric stretch in anSN2 displacement reaction. To begin our considerations, we can picture a reaction free energy profile in the solute reactive coordinate x calculated via the potential of mean force Geq(x) -the system free energy when the system is equilibrated at each fixed value of x, which would be the output of e.g. equilibrium Monte Carlo or Molecular Dynamics calculations [25] or equilibrium integral equation methods [26], Attention then focusses on the barrier top in this profile, located at x. ... 

Xu, J. and S. B. Pope (1999). Assessment of numerical accuracy of PDF/Monte Carlo methods for turbulent reactive flows. Journal of Computational Physics 152, 192-230. 

The choice of other variables R, r, h, 0, and r appropriate for Monte-Carlo averaging is made by pseudo random numbers generated on computer. The reactive cross section can be found by averaging the reaction probability over the impact parameter and rotational state... 

Figure 10. Average number of reactive contacts calculated by the Monte Carlo simulation (X) and by the approximate approach (O) (N = 100)...

Chandler D, Wolynes PG (1981) J Chem Phys 74 4078. Quantum Path Integral Monte Carlo can also provide a way to evaluhte the rate in non-adiabatic reactive flux correlation theory Wolynes PG (1987) J Chem Phys 87 6559... 

Computational efforts using DPT calculations as well as kinetic modeling of reactivities based on Monte Carlo simulations or mean field mefh-ods have been employed to study elementary processes on Pt surfaces. 2 228 Unraveling systematic trends in structure versus reactivity relations remains a formidable challenge due to fhe complex nafure of sfrucfural effects in electrocatalysis. 

It Is to be remarked that the process described by the Infinite set of kinetic (coagulation) equations can be simulated by Monte-Carlo methods ( ). The Information on the number of molecules of the respective size Is stored In the computer memory and weighting for selection of molecules Is applied given by the number and reactivity of groups In the respective molecule. 

The kinetics data of the geminate ion recombination in irradiated liquid hydrocarbons obtained by the subpicosecond pulse radiolysis was analyzed by Monte Carlo simulation based on the diffusion in an electric field [77,81,82], The simulation data were convoluted by the response function and fitted to the experimental data. By transforming the time-dependent behavior of cation radicals to the distribution function of cation radical-electron distance, the time-dependent distribution was obtained. Subsequently, the relationship between the space resolution and the space distribution of ionic species was discussed. The space distribution of reactive intermediates produced by radiation is very important for advanced science and technology using ionizing radiation such as nanolithography and nanotechnology [77,82]. 

---