The Monte Carlo MC method

In molecular dynamics, the computational pull that drives the system in its exploration of phase space is the passing of time consequent to the solution of Newton s equation of motion. This repetitive task is already very dull, but another computational technique for molecular simulation, the Monte Carlo method [4], rests upon an even 

---

The Monte Carlo (MC) method can be used to efficiently calculate thermal equilibrium properhes (see Fig. 3.2). However, since it is an energy-barrier-based method, it will fail to generate dynamic features such as the precession of the spins, and will be able to generate the dynamic magnetizahon in the overdamped limit (X —> oo) only if an appropriate algorithm is used [35]. 

Polymerization rate represents the instantaneous status of reaction locus, but the whole history of polymerization is engraved within the molecular weight distribution (MWD). Recently, a new simulation tool that uses the Monte Carlo (MC) method to estimate the whole reaction history, for both hnear [263-265] and nonlinear polymerization [266-273], has been proposed. So far, this technique has been applied to investigate the kinetic behavior after the nucleation period, where the overall picture of the kinetics is well imderstood. However, the versatility of the MC method could be used to solve the complex problems of nucleation kinetics. 

The concept of liquidlike clusters seems to have first emerged as a result of the explosion of simulations by the molecular dynamics (MD) method of McGinty, ° Cotterill et al., Damgaard Kristensen et al., and Briant and Burton, and by the Monte Carlo (MC) method of Lee et al. These were followed quickly by further MD simulations and MC simula-... 

In classical mechanics there exist, apart from the mean field theory, two popular methods to describe the dynamics of molecular systems, viz., the molecular dynamics (MD) method and the Monte Carlo (MC) method (Hansen and McDonald, 1976). In both methods the system is represented by a finite number, usually about 100 to 300, of molecules. In order to reduce boundary effects, this finite system is periodically repeated in all directions. 

The first molecular simulations were performed almost five decades ago by Metropolis et al. (1953) on a system of hard disks by the Monte Carlo (MC) method. Soon after, hard spheres (Rosenbluth and Rosenbluth, 1954) and Lennard-Jones (Wood and Parker, 1957) particles were also studied by both MC and molecular dynamics (MD). Over the years, the simulation techniques have evolved to deal with more complex systems by introducing different sampling or computational algorithms. Molecular simulation studies have been made of molecules ranging from simple systems (e.g., noble gases, small organic molecules) to complex molecules (e.g., polymers, biomolecules). 

The Monte Carlo (MC) method gives a qnalitative description of collective diffusive motion in terms of probabilities. MC calculations are simple, yet powerfirl enough to describe non-equihbrium adatom-surface interactions. 

The Monte Carlo (MC) method, used to simulate the properties of liquids, was developed by Metropolis et al. (1953). Without going into any detail, it should be pointed out that there are two important features of this MC method that make it extremely useful for the study of the liquid state. One is the use of periodic boundary conditions, a feature that helps to minimize the surface effects that are likely to be substantial in such a small sample of particles. The second involves the way the sample of configurations are selected. In the authors words Instead of choosing configurations randomly, then weighing them with exp[—/i ], we choose configurations with probability exp[—/6 ] and weight them evenly. ... 

Previous works dealing with disordered sur ces have been dedicated mainly to random, or correlated topographies. In the latter case, the combination of heterogeneity and ad-ad interactions effects produce complex behaviors on the equilibrium properties. An exact statistical mechanical treatment is unfortunately not yet available and, therefore, the theoretical description of adsorption has relied on simplified models. One way of circumventing this complication is the Monte Carlo (MC) method, which has demonstrated to be a valuable tool to study surface processes [3,4]. 

Two computer simulation methods are widely used for ionic solutions the Monte Carlo (MC) method and molecular dynamics (MD). 

Molecular simulations Computer modeling of the motion of an assembly of atoms or molecules. In molecular simulations only the motions of nuclei are considered, (i.e one assumes that the electronic Schrodinger equation has been solved providing intermolecular interaction potentials.) In practice empirical interaction potentials are utilized in most cases. Two main approaches are used the Monte Carlo (MC) method and molecular dynamics (MD). The former relies on statistical sampling of the configuration space of the systems, whereas the latter solves classical mechanics (Newton) equations to find trajectories of molecules. 

The specific method devised by Metropolis et aL (1953) to compute the properties of liquids is now known as the Monte Carlo (MC) method. In fact, this is a special procedure to compute multidimensional integrals numerically. 

The Monte Carlo (MC) method starts with a particular arrangement of all the particles (solute and solvent molecules) in the system—a configuration. Then, a three-step procedure is applied. 

The above discussion on the two components of should lead to a better understanding of physical adsorption. Theoretically, polymer adsorption(so) can be treated by the Scheutjens-Fleer (SF)(si) mean-field theory, the Monte Carlo (MC) method,( 2) or the scaling approach. (83) In Figure 10, two profiles are given for the cases of adsorption (x = 1) and depletion (x = 0) using the SF theory, where x is the Flory-Huggins interaction parameter(84) between a polymer and a solvent with respect to pure components. The polymer coil expands if X < 0.5 and contracts if x These two cases are referred to as good and poor solvents, respectively. From the volume fraction profile c )(z), we can calculate other adsorption parameters, such as F, the adsorbed amount ... 

A correct calculation of solvation thermodynamics and solution structure is conceivable only in terms of the methods of statistical physics, in particular, the computer experiment schemes, including, in the first place, the molecular dynamics (MD) and the Monte-Carlo (MC) methods [10]. By means of the MD method Newton s classic equations of motion are solved numerically with the aid of a computer assuming that the potential energy of molecular interaction is known. In this manner, the motion of molecules of the liquid may be observed , the phase trajectories found and then the values of the necessary functions are averaged over time and determined. This method permits both the equilibrium and the kinetical properties of the system to be calculated. 

The most rigorous calculations of the mechanisms of 8 2 reactions in solution were conducted by Jorgensen and co-workers [55,56]. They examined as an example the process Cl" -h CH3CI using the ab initio (6-3IG basis set) and the Monte Carlo (MC) methods. So as to carry out the MC calculations of the reaction in question in water and dimethyl formamide (DMF), the following problems had to be solved. 

The Monte Carlo (MC) method for predicting isotherms in a range of pore sizes at different pressures at a given temperature consists of randomly moving and rotating molecules in a small section of a model pore. The equations for describing the sorbate-sorbate and sorbate-sorbent interactions are usually the same as those described for the DFT method (LJ potential and Steele s 10-4-3 potential respectively). 

One of the ways of circumventing the problem of finding multiple energy minima of complex molecules is to turn to more sophisticated techniques that are capable of sampling phase space efficiently without the need to home in on particular minimum energy conformations. The two most useful techniques are molecular dynamics (MD) and the Monte Carlo (MC) method. Both approaches make use of the same types of potential functions used in molecular mechanics, but are designed to sample conformation space such that a Boltzmann distribution of states is generated. MC and MD techniques for molecular systems have been widely reviewed [11-14], and only the basics of the two methods are described below. 

Laureau, A., Aufiero, M., Rubiolo, P., Merle-Lucotte, E., Heuer, D., 2015a. Coupled neutronics and thermal-hydraubcs transient calculations based on a fission matrix approach application to the molten salt fast reactor. In Proceedings of the Joint International Conference on Mathematics and Computation (M C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method, Nashville, USA. 

To perform high-dimensional integrals numerically, the Monte Carlo (MC) method is often used. The MC method turns a multidimensional integral into a sum over a stochastic sequence of points called a trajectory, so that in the limit of an infinitely long trajectory, the value of the integral is numerically reproduced. Explicitly, the MC algorithm yields the ratio of two integrals... 

If one is interested in equilibrium canonical (fixed temperature) properties of liquid interfaces, an approach to sample phase space is the Monte Carlo (MC) method. Here, only the potential energy function l/(ri,r2,. .., rjy) is required to calculate the probability of accepting random particle displacement moves (and additional moves depending on the ensemble type ). All of the discussion above regarding the boundary conditions, treatment of long-range interactions, and ensembles applies to MC simulations as well. Because the MC method does not require derivatives of the potential energy function, it is simpler to implement and faster to run, so early simulations of liquid interfaces used However, dynamical information is not available with... 

---