Monte Carlo method, mixing

In the last few years we have witnessed the successful development of several methods for the numerical solution of multi-dimensional quantum Hamiltonians Monte Carlo methods centroid methods,mixed quantum-classical methods, and recently a revival of semiclassical methods. We have developed another approach to this problem, the exponential resummation of the evolution operator. - The rest of this Section will explain briefly this method. 

Equation (3.21) shows that the potential of the mean force is an effective potential energy surface created by the solute-solvent interaction. The PMF may be calculated by an explicit treatment of the entire solute-solvent system by molecular dynamics or Monte Carlo methods, or it may be calculated by an implicit treatment of the solvent, such as by a continuum model, which is the subject of this book. A third possibility (discussed at length in Section 3.3.3) is that some solvent molecules are explicit or discrete and others are implicit and represented as a continuous medium. Such a mixed discrete-continuum model may be considered as a special case of a continuum model in which the solute and explicit solvent molecules form a supermolecule or cluster that is embedded in a continuum. In this contribution we will emphasize continuum models (including cluster-continuum models). 

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. 

The Monte Carlo method is a well-established approach used in characterizing uncertainty that can be used to incorporate ranges of data (distributions) into calculations. The Monte Carlo method involves choosing values from a random selection scheme drawn from probability density functions based on a range of data that characterize the parameters of interest. Monte Carlo analysis can be selectively used to generate input parameters and mixed with point estimates, as appropriate, to calculate risk. 

Finally, we stress that the quantum chemical method presented here has the advantage over DFT-based techniques that it also furnishes wavefunctions that can be used to perform computations of spectra, and therefore have a better contact with the experiment. Another advantage of this approach is that, unlike the diffusion Monte-Carlo method, it can coherently be applied to studies of fermion and mixed boson/fermion doped clusters. An example can be found in our recent work on the Raman spectra of (He)w-Br2(X) clusters [27,28]. 

Mixed effects models under a Bayesian framework have been widely studied and used with the use of Markov chain Monte Carlo methods (10). These methods have gained particular popularity as complex problems became easily formulated using the WinBUGS software (11). See Congdon (12) for an extensive coverage of topics and examples and implementation in WinBUGS. 

Simulation of a mixing process by using a computer provides a better alternative. Monte Carlo simulation techniques are often utilized. Monte Carlo techniques are numerical methods that involve sampling from statistical distributions, either theoretical or empirical, to approximate the real physical phenomena without reference to the actual physical systems. For the general discussion of Monte Carlo methods, readers are referred to Hammersley and Handscomb [19] and Tocher [20]. 

Wang MR, Li ZX (2006) Gas mixing in microchannels using the direct simulation Monte Carlo method. Int J Heat Mass Transf 49(9-10) 1696-1702... 

In the harmonic method and its extensions it is always assumed that the amplitudes of the molecular vibrations about their equilibrium positions and orientations remain small. It will be illustrated by several examples in Sect. 2.3 that this is often not realistic for molecular crystals. In the plastic phases there is not even a well-defined equilibrium orientation of the molecules, but also in the ordered phases the librational amplitudes may become substantial. The motions in such systems have been studied by classical methods, in particular the Molecular Dynamics method and the Monte Carlo method [73]. The advantage of these methods is that they can also be applied to study liquids, and the melting of solids, and other systems (glasses, solutions, mixed crystals) which have lost translational periodicity. Large amplitude motions in molecular crystals can also be studied quantum mechanically, however, by the methods described below. 

A general transport phenomenon in the intercalation electrode with a fractal surface under the constraint of diffusion mixed with interfadal charge transfer has been modelled by using the kinetic Monte Carlo method based upon random walk approach (Lee Pyim, 2005). Go and Pyun (Go Pyun, 2007) reviewed anomalous diffusion towards and from fractal interface. They have explained both the diffusion-controlled and non-diffusion-controlled transfer processes. For the diffusion coupled with facile charge-transfer reaction the... 

There are basically two different computer simulation techniques known as molecular dynamics (MD) and Monte Carlo (MC) simulation. In MD molecular trajectories are computed by solving an equation of motion for equilibrium or nonequilibrium situations. Since the MD time scale is a physical one, this method permits investigations of time-dependent phenomena like, for example, transport processes [25,61-63]. In MC, on the other hand, trajectories are generated by a (biased) random walk in configuration space and, therefore, do not per se permit investigations of processes on a physical time scale (with the dynamics of spin lattices as an exception [64]). However, MC has the advantage that it can easily be applied to virtually all statistical-physical ensembles, which is of particular interest in the context of this chapter. On account of limitations of space and because excellent texts exist for the MD method [25,61-63,65], the present discussion will be restricted to the MC technique with particular emphasis on mixed stress-strain ensembles. 

The differences in the rate constant for the water reaction and the catalyzed reactions reside in the mole fraction of substrate present as near attack conformers (NACs).171 These results and knowledge of the importance of transition-state stabilization in other cases support a proposal that enzymes utilize both NAC and transition-state stabilization in the mix required for the most efficient catalysis. Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82%, 57%, and 1% of chorismate conformers were found to be NAC structures (NACs) in water, methanol, and the gas phase, respectively.172 The fact that the reaction occurred faster in water than in methanol was attributed to greater stabilization of the TS in water by specific interactions with first-shell solvent molecules. The Claisen rearrangements of chorismate in water and at the active site of E. coli chorismate mutase have been compared.173 It follows that the efficiency of formation of NAC (7.8 kcal/mol) at the active site provides approximately 90% of the kinetic advantage of the enzymatic reaction as compared with the water reaction. 

---