Statistical mechanics Carlo technique

A comprehensive introduction to the field, covering statistical mechanics, basic Monte Carlo, and molecular dynamics methods, plus some advanced techniques, including computer code. 

The partition function Z is given in the large-P limit, Z = limp co Zp, and expectation values of an observable are given as averages of corresponding estimators with the canonical measure in Eq. (19). The variables and R ( ) can be used as classical variables and classical Monte Carlo simulation techniques can be applied for the computation of averages. Note that if we formally put P = 1 in Eq. (19) we recover classical statistical mechanics, of course. 

In this review we put less emphasis on the physics and chemistry of surface processes, for which we refer the reader to recent reviews of adsorption-desorption kinetics which are contained in two books [2,3] with chapters by the present authors where further references to earher work can be found. These articles also discuss relevant experimental techniques employed in the study of surface kinetics and appropriate methods of data analysis. Here we give details of how to set up models under basically two different kinetic conditions, namely (/) when the adsorbate remains in quasi-equihbrium during the relevant processes, in which case nonequilibrium thermodynamics provides the needed framework, and (n) when surface nonequilibrium effects become important and nonequilibrium statistical mechanics becomes the appropriate vehicle. For both approaches we will restrict ourselves to systems for which appropriate lattice gas models can be set up. Further associated theoretical reviews are by Lombardo and Bell [4] with emphasis on Monte Carlo simulations, by Brivio and Grimley [5] on dynamics, and by Persson [6] on the lattice gas model. 

Beyond the clusters, to microscopically model a reaction in solution, we need to include a very big number of solvent molecules in the system to represent the bulk. The problem stems from the fact that it is computationally impossible, with our current capabilities, to locate the transition state structure of the reaction on the complete quantum mechanical potential energy hypersurface, if all the degrees of freedom are explicitly included. Moreover, the effect of thermal statistical averaging should be incorporated. Then, classical mechanical computer simulation techniques (Monte Carlo or Molecular Dynamics) appear to be the most suitable procedures to attack the above problems. In short, and applied to the computer simulation of chemical reactions in solution, the Monte Carlo [18-21] technique is a numerical method in the frame of the classical Statistical Mechanics, which allows to generate a set of system configurations... 

We will delay a more detailed discussion of ensemble thermodynamics until Chapter 10 indeed, in this chapter we will make use of ensembles designed to render the operative equations as transparent as possible without much discussion of extensions to other ensembles. The point to be re-emphasized here is that the vast majority of experimental techniques measure molecular properties as averages - either time averages or ensemble averages or, most typically, both. Thus, we seek computational techniques capable of accurately reproducing these aspects of molecular behavior. In this chapter, we will consider Monte Carlo (MC) and molecular dynamics (MD) techniques for the simulation of real systems. Prior to discussing the details of computational algorithms, however, we need to briefly review some basic concepts from statistical mechanics. 

The Monte Carlo method is a very powerful numerical technique used to evaluate multidimensional integrals in statistical mechanics and other branches of physics and chemistry. It is also used when initial conditions are chosen in classical reaction dynamics calculations, as we have discussed in Chapter 4. It will therefore be appropriate here to give a brief introduction to the method and to the ideas behind the method. 

The same computer revolution that started in the middle of the last century also plays an important, in fact crucial, role in the development of methods and algorithms to study solvation problems. Dealing, for instance, with a liquid system means the inclusion of explicit molecules, in different thermodynamic conditions. The number of possible arrangements of atoms or molecules is enormous, demanding the use of statistical mechanics. Here is where computer simulation, Monte Carlo (MC) or molecular dynamics (MD), makes its entry to treat liquid systems. Computer simulation is now an important, if not central, tool to study solvation phenomena. The last two decades have seen a remarkable development of methods, techniques and algorithms to study solvation problems. Most of the recent developments have focused on combining quantum mechanics and statistical mechanics using MC or... 

Other review articles relating to colloidal suspensions and containing discussions of numerical calculations (Monte Carlo simulations) include contributions by van Megen and Snook (1984) and Castillo et al. (1984). The scope of these articles is not, however, limited to numerical techniques they also provide general reviews of the statistical mechanics of colloidal suspensions. 

The two simulation methods in general use for solving the statistical mechanical equations are Monte Carlo (MC) and molecular dynamics (MD). The two techniques have several common features, but each has certain advantages and limitations. 

Valleau, J. P. and Torrie, G. M., A guide to Monte Carlo for statistical mechanics 2. Byways. In Statistical Mechanics, Part A Equilibrium Techniques, pp. 169-194. New York Plenum (1977). 

Therefore molecular dynamics is a deterministic technique given an initial set of positions and velocities, the subsequent time evolution is in principle completely determined. The computer calculates a trajectory in a 6A-dimensional phase space (3A positions and 3A momenta). However, such trajectory is usually not particularly relevant by itself. Molecular dynamics is a statistical mechanics method. Like Monte Carlo, it is a way to obtain a set of configurations distributed according to some statistical distribution function, or statistical ensemble. 

The next section gives a brief overview of the main computational techniques currently applied to catalytic problems. These techniques include ab initio electronic structure calculations, (ab initio) molecular dynamics, and Monte Carlo methods. The next three sections are devoted to particular applications of these techniques to catalytic and electrocatalytic issues. We focus on the interaction of CO and hydrogen with metal and alloy surfaces, both from quantum-chemical and statistical-mechanical points of view, as these processes play an important role in fuel-cell catalysis. We also demonstrate the role of the solvent in electrocatalytic bondbreaking reactions, using molecular dynamics simulations as well as extensive electronic structure and ab initio molecular dynamics calculations. Monte Carlo simulations illustrate the importance of lateral interactions, mixing, and surface diffusion in obtaining a correct kinetic description of catalytic processes. Finally, we summarize the main conclusions and give an outlook of the role of computational chemistry in catalysis and electrocatalysis. 

There are two main approaches used to simulate polymer materials molecular dynamics and Monte Carlo methods. The molecular dynamics approach is based on numerical integration of Newton s equations of motion for a system of particles (or monomers). Particles follow dctcr-ministic trajectories in space for a well-defined set of interaction potentials between them. In a qualitatively different simulation technique, called Monte Carlo, phase space is sampled randomly. Molecular dynamics and Monte Carlo simulation approaches are analogous to time and ensemble methods of averaging in statistical mechanics. Some modern computer simulation methods use a combination of the two approaches. 

Often we are only interested in the equilibrium structure of a set of molecules. When the system is at high densities or includes a large number of conformations, other methods become computationally unfeasible, and one often uses Monte Carlo (MC) techniques. Monte Carlo methods use random number generators to integrate systems with very high degrees of freedom. These integrals are then used in theories of statistical mechanics to evaluate thermodynamic properties of materials. ... 

---