Molecular Dynamics or Monte Carlo simulations

Fig. 5. To generate an ensemble using Molecular Dynamics or Monte-Carlo simulation techniques the interaction between all pairs of atoms within a given cutoff radius must be considered. In contrast, to estimate changes in free energy using a stored trajectory only those interactions which are perturbed need be determined making the approach highly efficient.

In some cases the atomic charges are chosen to reproduce thermodynamic properties calculated using a molecular dynamics or Monte Carlo simulation. A series of simulations is performed and the charge model is modified until satisfactory agreement with experiment is obtained. This approach can be quite powerful despite its apparent simplicity, but it is only really practical for small molecules or simple models. 

TIk experimentally determined dipole moment of a water molecule in the gas phase is 1.85 D. The dipole moment of an individual water molecule calculated with any of thv se simple models is significantly higher for example, the SPC dipole moment is 2.27 D and that for TIP4P is 2.18 D. These values are much closer to the effective dipole moment of liquid water, which is approximately 2.6 D. These models are thus all effective pairwise models. The simple water models are usually parametrised by calculating various pmperties using molecular dynamics or Monte Carlo simulations and then modifying the... 

Surface tension is usually predicted using group additivity methods for neat liquids. It is much more difficult to predict the surface tension of a mixture, especially when surfactants are involved. Very large molecular dynamics or Monte Carlo simulations can also be used. Often, it is easier to measure surface tension in the laboratory than to compute it. 

Conventional molecular dynamics or Monte Carlo simulations. 

The idea of a finite simulation model subsequently transferred into bulk solvent can be applied to a macromolecule, as shown in Figure 5a. The alchemical transformation is introduced with a molecular dynamics or Monte Carlo simulation for the macromolecule, which is solvated by a limited number of explicit water molecules and otherwise surrounded by vacuum. Then the finite model is transferred into a bulk solvent continuum... 

From this short discussion, it is clear that atomistically detailed molecular dynamics or Monte Carlo simulations can provide a wealth of information on systems on a local molecular atomistic level. They can, in particular, address problems where small changes in chemical composition have a drastic effect. Since chemical detail is avoided in mesoscopic models, these can often capture such effects only indirectly. 

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. 

Techniques such as those described in Section III. 1, involving molecular dynamics or Monte Carlo simulations, have been used extensively to study ionic solvation. A classical... 

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 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. 

So far, our discussion has focussed on stationary quantum chemical methods, which yield results for fixed atomic nuclei, i.e. for frozen molecular structures like minimum structures on the Born-Oppenheimer potential energy surface. Processes in supramolecular assemblies usually feature prominent dynamical effects, which can only be captured through explicit molecular dynamics or Monte Carlo simulations [95-98]. Molecular dynamics simulations proved to be a useful tool for studying the detailed microscopic dynamic behavior of many-particle systems as present in physics, chemistry and biology. The aim of molecular dynamics is to study a system by recreating it on the computer as close to nature as possible, i.e. by simulating the dynamics of a system in all microscopic detail over a physical length of time relevant to properties of interest. 

---