Computer simulation Monte Carlo calculations

Now, let us look at Fig. 13. Here, the static structure factor of a three-dimensional homogeneous suspension of polystyrene spheres of diameter 94 nm is shown. The particles volume fraction is 0 = 2.0 x 10 4. Experimental data from static light scattering (closed circles) are compared with computer simulation (Monte Carlo) results (symbol x) and theoretical predictions (lines) obtained from the Ornstein-Zernike equation and different closure relations. The computer simulations and the theoretical calculations where carried out assuming that the interaction between the... 

Recently, molecular dynamics and Monte Carlo calculations with quantum mechanical energy computation methods have begun to appear in the literature. These are probably some of the most computationally intensive simulations being done in the world at this time. 

Molecular dynamics calculations are more time-consuming than Monte Carlo calculations. This is because energy derivatives must be computed and used to solve the equations of motion. Molecular dynamics simulations are capable of yielding all the same properties as are obtained from Monte Carlo calculations. The advantage of molecular dynamics is that it is capable of modeling time-dependent properties, which can not be computed with Monte Carlo simulations. This is how diffusion coefficients must be computed. It is also possible to use shearing boundaries in order to obtain a viscosity. Molec-... 

In a statistical Monte Carlo simulation the pair potentials are introduced by means of analytical functions. In the election of that analytical form for the pair potential, it must be considered that when a Monte Carlo calculation is performed, the more time consuming step is the evaluation of the energy for the different configurations. Given that this calculation must be done millions of times, the chosen analytic functions must be of enough accuracy and flexibility but also they must be fastly computed. In this way it is wise to avoid exponential terms and to minimize the number of interatomic distances to be calculated at each configuration which depends on the quantity of interaction centers chosen for each molecule. A very commonly used function consists of a sum of rn terms, r being the distance between the different interaction centers, usually, situated at the nuclei. In particular, non-bonded interactions are usually represented by an atom-atom centered monopole expression (Coulomb term) plus a Lennard-Jones 6-12 term, as indicated in equation (51). 

During the last few years the progress of computational techniques has made it possible to simulate the dynamic behavior of whole ensembles consisting of several hundred molecules. In this way the limitations of the statistical approach can be at least partly overcome. Two kinds of methods — molecular dynamics and Monte Carlo calculations — were applied to liquids and liquid mixtures and brought new insight into their structure and properties. Even some important characteristics of systems as complicated as associated liquids like water could be... 

Wall, Hiller, and Mazur [300, 301] first used a computer to integrate the classical motion equations for a system of three atoms, and in the 1960s the technique was developed by Blais and Bunker [48, 302-306], and by Karplus [19, 20, 72, 307-311] and Polanyi [71,73, 74, 267, 312, 313] and their coworkers. Recently calculations have been performed on systems simulating abstraction reactions involving more than three atoms [314] and four-center reactions involving four atoms, that is, AB + CD - AC + BD [315-317]. Here we present first a general survey of the Monte Carlo calculations of classical trajectories and then a brief review of some of the results of these calculations. Emphasis is placed on data for reactions that have been studied experimentally and have been mentioned earlier in this chapter. 

The United States is currently a leader in most areas of theoretical/ computational chemistry. In basic theory, Europe has many talented young investigators. Within the next 10 years, given these demographics, the U.S. leadership will be challenged by Europe in electronic structure and basic theory development. This trend does not characterize the entire field of theoretical chemistry. For example, Monte Carlo and molecular dynamics simulation methods were invented in the United States, and to this day, the United States maintains a strong position, especially in quantum Monte Carlo calculations. 

Computer simulations of self-avoiding ring chains have given information on the ratio (5 (r))/(5 (/)) at the limit of large chain length. By suitable extrapolation of data from Monte Carlo simulations of off-lattice chains of N (the number of bonds) up to 99, Bmns and Naghizadeh [78] estimated this ratio to be 0.559. Chen [79] found 0.568 on the basis of Monte Carlo calculations on... 

The two basic sampling methods used for such studies are (1) Monte Carlo calculations, which allow one to obtain, among other properties, adsorption energies and preferred locations for the adsorbed species, and (2) molecular dynamics simulations, which provide dynamical information such as diffusion coefficients. Both methods require repeated calculation of the energy of the system, and accordingly this is still an area dominated by force field calculations. A few studies have been performed by means of density functional calculations, as exemplified by Haase, Sauer, and Hutter, but these calculations require at present substantial computer power. 

The first Monte Carlo calculations were performed on the Eniac by a group of people from Los Alamos, led by N. C. Metropolis. Those calculations were rather primitive, compared with what can be done on today s computing machines but, in one way, they were about as sophisticated as any ever performed, in that they simulated complete chain reactions in critical and supercritical systems, starting with an assumed neutron distribution, in space and velocity, at an initial instant of time and then following all details of the reaction as it develops subsequently. 

---