Monte Carlo/reactive dynamics

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.)... 

The process of adsorption and interaction of probe molecules such as ammonia, carbon monoxide as well as the whole spectrum of organic reactant molecules with zeolite catalysts has been the subject of numerous experimental and computational studies. These interaction processes are studied using several computational methods involving force fields (Monte Carlo, molecular dynamics emd energy minimization) or quantum chemical methods. Another paper [1] discusses the application of force field methods for studying several problems in zeolite chemistry. Theoretical quantum chemical studies on cluster models of zeolites help us to understand the electronic and catalytic properties of zeolite catalysts. Here we present a brief summary of the application of quantum chemical methods to understand the structure and reactivity of zeolites. 

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... 

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. ... 

The need to reliably describe liquid systems for practical purposes as condensed matter with high mobility at a given finite temperature initiated attempts, therefore, to make use of statistical mechanical procedures in combination with molecular models taking into account structure and reactivity of all species present in a liquid and a solution, respectively. The two approaches to such a description, namely Monte Carlo (MC) simulations and molecular dynamics (MD), are still the basis for all common theoretical methods to deal with liquid systems. While MC simulations can provide mainly structural and thermodynamical data, MD simulations give also access to time-dependent processes, such as reaction dynamics and vibrational spectra, thus supplying — connected with a higher computational effort — much more insight into the properties of liquids and solutions. 

Bifurcation and global stability in surface catalyzed reactions using the Monte Carlo method (with D.G. Vlachos and L.D. Schmidt). In H. Swinney, R. Aris, and D. Aronson (eds.), Patterns and Dynamics in Reactive Media, (Vol. 37, pp. 187-206). New York Springer-Verlag, 1991. 

Solvent effects can significantly influence the function and reactivity of organic molecules.1 Because of the complexity and size of the molecular system, it presents a great challenge in theoretical chemistry to accurately calculate the rates for complex reactions in solution. Although continuum solvation models that treat the solvent as a structureless medium with a characteristic dielectric constant have been successfully used for studying solvent effects,2,3 these methods do not provide detailed information on specific intermolecular interactions. An alternative approach is to use statistical mechanical Monte Carlo and molecular dynamics simulation to model solute-solvent interactions explicitly.4 8 In this article, we review a combined quantum mechanical and molecular mechanical (QM/MM) method that couples molecular orbital and valence bond theories, called the MOVB method, to determine the free energy reaction profiles, or potentials of mean force (PMF), for chemical reactions in solution. We apply the combined QM-MOVB/MM method to... 

Recent advances have resulted from the development of more powerful experimental methods and because the classical collision dynamics can now be calculated fully using high-speed computers. By applying Monte Carlo techniques to the selection of starting conditions for trajectory calculations, a reaction can be simulated with a sample very much smaller than the number of reactive encounters that must necessarily occur in any kinetic experiment, and models for reaction can therefore be tested. The remainder of this introduction is devoted to a simple explanation of the classical dynamics of collisions, a description of the parameters needed to define them, and the relationship between these and the rate coefficient for a reaction [9]. 

Whereas selective diffusion can be better investigated using classical dynamic or Monte Carlo simulations, or experimental techniques, quantum chemical calculations are required to analyze molecular reactivity. Quantum chemical dynamic simulations provide with information with a too limited time scale range (of the order of several himdreds of ps) to be of use in diffusion studies which require time scale of the order of ns to s. However, they constitute good tools to study the behavior of reactants and products adsorbed in the proximity of the active site, prior to the reaction. Concerning reaction pathways analysis, static quantum chemistry calculations with molecular cluster models, allowing estimates of transition states geometries and properties, have been used for years. The application to solids is more recent. 

Finally, dynamic Monte Carlo simulations are very useful in assessing the overall reactivity of a catalytic surface, which must include the effects of lateral interactions between adsorbates and the mobility of adsorbates on the surface in reaching the active sites. The importance of treating lateral interactions was demonstrated in detailed ab initio-based dynamic Monte Carlo simulations of ethylene hydrogenation on palladium and PdAu alloys. Surface diffusion of CO on PtRu alloy surfaces was shown to be essential to explain the qualititative features of the experimental CO stripping voltammetry. Without adsorbate mobility, these bifunctional surfaces do not show any catalytic enhancement with respect to the pure metals. 

Recent work (Altmann, 2002) has focussed on combining a dynamic Monte Carlo percolation-grid simulation for reaction kinetics and an enthalpy-based group-interaction viscoelastic model to develop a model for the chemorheological and network properties of reactive systems. More emphasis will be placed on this model in Chapter 6. 

---