Monte Carlo simulations principles

In principle, simulation teclmiques can be used, and Monte Carlo simulations of the primitive model of electrolyte solutions have appeared since the 1960s. Results for the osmotic coefficients are given for comparison in table A2.4.4 together with results from the MSA, PY and HNC approaches. The primitive model is clearly deficient for values of r. close to the closest distance of approach of the ions. Many years ago, Gurney [H] noted that when two ions are close enough together for their solvation sheaths to overlap, some solvent molecules become freed from ionic attraction and are effectively returned to the bulk [12]. 

Predictive modeling [38], Tachugi design principles [2], Monte Carlo simulations to simulate impacts of different product and process conditions on Q attribute level [40]... 

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. 

In principle, the ideal description of a solution would be a quantum mechanical treatment of the supermolecule consisting of representative numbers of molecules of solute and solvent. In practice this is not presently feasible, even if only a single solute molecule is included. In recent years, however, with the advances in processor technology that have occurred, it has become possible to carry out increasingly detailed molecular dynamics or Monte Carlo simulations of solutions, involving hundreds or perhaps even thousands of solvent molecules. In these, all solute-solvent and solvent-solvent interactions are taken into account, at some level of sophistication. 

Eurthermore, uncertainties in the exposure assessment should also be taken into account. However, no generally, internationally accepted principles for addressing these uncertainties have been developed. For predicted exposure estimates, an uncertainty analysis involving the determination of the uncertainty in the model output value, based on the collective uncertainty of the model input parameters, can be performed. The usual approach for assessing this uncertainty is the Monte Carlo simulation. This method starts with an analysis of the probability distribution of each of the variables in the uncertainty analysis. In the simulation, one random value from each distribution curve is drawn to produce an output value. This process is repeated many times to produce a complete distribution curve for the output parameter. 

K. Reuter and M. Scheffler, First-Principles Kinetic Monte Carlo Simulations for Heterogeneous Catalysis Application to the CO Oxidation at RuO2(110), Phys. Rev. B 73 (2006), 045433. 

Burmaster and Anderson (1994) have compiled a list of principles of good practice, which were originally aimed at Monte Carlo simulations, but are valuable also for other techniques in uncertainty analysis. These recommendations later appeared in modified and supplemented form in various handbooks and other publications, e.g., the USEPA Guiding Principles for Monte Carlo Analysis (USEPA 1997). 

Molecular Dynamics and Monte Carlo simulations have been used to predict the adsorption isotherms and transport diffusivities of Xe adsorbed in A1P04-31. The results of these calculations can be used to predict the properties of Xe diffusing through membranes made from A1P04-31 crystals directly from atomic-scale principles. 

The EAM and MEAM potentials once determined from electronics principles calculations [178] have been used to reproduce physical properties of many metals, defects, and impurities. For example, EAM molecular statics, molecular dynamics, and Monte Carlo simulations were performed on hydrogen embrittlement effects on dislocation motion and plasticity [46,179-181]. These potentials have been used to analyze plasticity [74,144,145,148-150,182,183], cracks and fracture [117,184], and fatigue [119, 120, 185, 186]. 

The Monte-Carlo principle uses random numbers in the region 0 < random < 1. The main randomized values in the Monte-Carlo simulation, see Figure 1, of a gas-phase adsorption process for radioactive species are ... 

Reuter K, Scheffler M. First-principles kinetic Monte Carlo simulations for heterogeneous catalysis application to the CO oxidation at RuC>2(110). Phys Rev B. 2006 73(4). 

The polarizable point dipole model has also been used in Monte Carlo simulations with single particle moves.When using the iterative method, a whole new set of dipoles must be computed after each molecule is moved. These updates can be made more efficient by storing the distances between all the particles, since most of them are unchanged, but this requires a lot of memory. The many-body nature of polarization makes it more amenable to molecular dynamics techniques, in which all particles move at once, compared to Monte Carlo methods where typically only one particle moves at a time. For nonpolarizable, pairwise-additive models, MC methods can be efficient because only the interactions involving the moved particle need to be recalculated [while the other (N - 1) x (]V - 1) interactions are unchanged]. For polarizable models, all N x N interactions are, in principle, altered when one particle moves. Consequently, exact polarizable MC calculations can be... 

In principle, the diffusion steps (a) and (e) could be studied through molecular dynamics simulations as long as rehable forces fields are available to describe the zeolite structure and its interaction with the substrates. Also, if the adsorption takes place without charge transfer between the reagents/products and the zeolite, steps (b) and (d) could also be investigated either by molecular dynamics or Monte Carlo simulations. Step (c) however can only be followed by quantum mechanical techniques because the available force fields cannot yet describe the breaking and formation of chemical bonds. 

In addition to the exposnre model documentation components noted above, the American Indnstrial Health Council (AIHC) (1994) and USEPA (1997a) have recommended data, particularly those based on Monte Carlo simulation, that are also relevant for simulation models being used as part of the overall assessment process. The USEPA has also issued guidelines for data quality assessment (USEPA, 1996b) relevant to model documentation. Some of these principles are listed below more details are provided in USEPA (1992a, 1997b), AIHC (1994) and Burmaster and Anderson (1994). 

Evidence for the Time-Temperature Superposition Principle from Monte Carlo Simulations of the Glass Transition in Two-Dimensional Polymer Melts. 

Snapshot of a Pd(lOO) surface during a simulation of ethylene hydrogenation at 298 K, 25 torr of ethylene and 100 torr of dihydrogen. [Figure from First-Principles-Based Monte Carlo Simulation of Ethylene Hydrogenation Kinetics on Pd, by E. W. Hansen and M. Neurock, m Journal of Catalysis, Volume 196 241, copyright 2000 by Academic Press, reproduced by permission of the publisher and authors.]... 

The aim of this chapter is to provide the reader with an overview of the potential of modern computational chemistry in studying catalytic and electro-catalytic reactions. This will take us from state-of-the-art electronic structure calculations of metal-adsorbate interactions, through (ab initio) molecular dynamics simulations of solvent effects in electrode reactions, to lattice-gas-based Monte Carlo simulations of surface reactions taking place on catalyst surfaces. Rather than extensively discussing all the different types of studies that have been carried out, we focus on what we believe to be a few representative examples. We also point out the more general theory principles to be drawn from these studies, as well as refer to some of the relevant experimental literature that supports these conclusions. Examples are primarily taken from our own work other recent review papers, mainly focused on gas-phase catalysis, can be found in [1-3]. 

Hansen, E., Neurock, M. First-principles-based Monte Carlo simulation of ethylene hydrogenation kinetics on Pd. J. Catal. 2000,1%, 241-52. 

---