Monte Carlo simulations molecular systems

This method avoids the convergence and accuracy problems of molecular dynamics or Monte Carlo simulations of systems containing explicit solvent molecules, by evaluating the electrostatic free energy of just one solute conformation surrounded by a dielectric continuum, and by adding the surface term and an estimate of the loss of the configurational entropy upon binding.77... 

In order to extract microscopic quantities such as spatial probabilities and RDFs from collected differential cross scattering data, analysis is typically performed following reverse Monte Carlo (RMC) or empirical potential structure refinement (EPSR) procedures [7]. The latter procedure can be viewed as a Monte Carlo simulation of system utilising a model potential similar to a classical molecular mechanics force field. This model potential is modified in order to bring the total structure factor calculated from the model system as close as possible to the ejq)erimental data. From configurations generated with this refined potential, standard quantities (such as RDFs) may be calculated. 

Two simulation methods—Monte Carlo and molecular dynamics—allow calculation of the density profile and pressure difference of Eq. III-44 across the vapor-liquid interface [64, 65]. In the former method, the initial system consists of N molecules in assumed positions. An intermolecule potential function is chosen, such as the Lennard-Jones potential, and the positions are randomly varied until the energy of the system is at a minimum. The resulting configuration is taken to be the equilibrium one. In the molecular dynamics approach, the N molecules are given initial positions and velocities and the equations of motion are solved to follow the ensuing collisions until the set shows constant time-average thermodynamic properties. Both methods are computer intensive yet widely used. 

Kremer K 1996 Computer simulation methods for polymer physics Monte Carlo and Molecular Dynamics of Condensed Matter Systems vol 49, ed K Binder and G Ciccotti (Bologna Italian Physical Society) pp 669-723... 

Monte Carlo simulations generate a large number of confonnations of tire microscopic model under study that confonn to tire probability distribution dictated by macroscopic constrains imposed on tire systems. For example, a Monte Carlo simulation of a melt at a given temperature T produces an ensemble of confonnations in which confonnation with energy E. occurs witli a probability proportional to exp (- Ej / kT). An advantage of tire Monte Carlo metliod is tliat, by judicious choice of tire elementary moves, one can circumvent tire limitations of molecular dynamics techniques and effect rapid equilibration of multiple chain systems [65]. Flowever, Monte Carlo... 

The input to a minimisation program consists of a set of initial coordinates for the system. The initial coordinates may come from a variety of sources. They may be obtained from an experimental technique, such as X-ray crystallography or NMR. In other cases a theoretical method is employed, such as a conformational search algorithm. A combination of experimenfal and theoretical approaches may also be used. For example, to study the behaviour of a protein in water one may take an X-ray structure of the protein and immerse it in a solvent bath, where the coordinates of the solvent molecules have been obtained from a Monte Carlo or molecular dynamics simulation. 

The correct treatment of boundaries and boundary effects is crucial to simulation methods because it enables macroscopic properties to be calculated from simulations using relatively small numbers of particles. The importance of boundary effects can be illustrated by considering the following simple example. Suppose we have a cube of volume 1 litre which is filled with water at room temperature. The cube contains approximately 3.3 X 10 molecules. Interactions with the walls can extend up to 10 molecular diameters into the fluid. The diameter of the water molecule is approximately 2.8 A and so the number of water molecules that are interacting with the boundary is about 2 x 10. So only about one in 1.5 million water molecules is influenced by interactions with the walls of the container. The number of particles in a Monte Carlo or molecular dynamics simulation is far fewer than 10 -10 and is frequently less than 1000. In a system of 1000 water molecules most, if not all of them, would be within the influence of the walls of the boundary. Clecirly, a simulation of 1000 water molecules in a vessel would not be an appropriate way to derive bulk properties. The alternative is to dispense with the container altogether. Now, approximately three-quarters of the molecules would be at the surface of the sample rather than being in the bulk. Such a situation would be relevcUit to studies of liquid drops, but not to studies of bulk phenomena. 

The Monte Carlo and molecular dynamics simulation methods can be used to explor the conformational space of molecules. During such a simulation the system is able t... 

Ihe allure of methods for calculating free energies and their associated thermod)mai values such as equilibrium constants has resulted in considerable interest in free ene calculations. A number of decisions must be made about the way that the calculatior performed. One obvious choice concerns the simulation method. In principle, eit Monte Carlo or molecular dynamics can be used in practice, molecular dynamics almost always used for systems where there is a significant degree of conformatio flexibility, whereas Monte Carlo can give very good results for small molecules which either rigid or have limited conformational freedom. 

Monte Carlo simulations require less computer time to execute each iteration than a molecular dynamics simulation on the same system. However, Monte Carlo simulations are more limited in that they cannot yield time-dependent information, such as diffusion coefficients or viscosity. As with molecular dynamics, constant NVT simulations are most common, but constant NPT simulations are possible using a coordinate scaling step. Calculations that are not constant N can be constructed by including probabilities for particle creation and annihilation. These calculations present technical difficulties due to having very low probabilities for creation and annihilation, thus requiring very large collections of molecules and long simulation times. 

---