Molecular-dynamics Calculations

The diffusion constant can be calculated through molecular dynamics simulations using the Einstein relation. 

The calculated numbers show an increase with temperature, which at low temperatures may be fitted to an Arrhenius type expression. The barrier depends, however, as we have seen on the direction of the diffusion process and hence if the initial diffusion direction is chosen at random only some average will be obtained. If the temperature is high such that a majority of molecules have energies above the barrier, then the Arrhenius-type behavior cannot be expected. This is also clearly seen from the numbers in the Table 3.2. The diffusion constant approaches a limiting value at high temperatures. 

TABLE 3.2 Diffusion Constant D(T) in cm /sec Obtained in [5 for Physisorbed Hydrogen Diffusion on a Cu(111) Surface 

Molecular mechanics still gives an inherently enthalpic view of the world. Molecular dynamics (MD) calculations should give more reliable and rigorous predictions of thermodynamic properties, as they deal with entropy by averaging over time and temperature. Ideally, a full atomistic model, including solvent, should be used. This can be prohibitively expensive, so solvent-continuum models such as Poisson-Boltzmann or are 

Free energy perturbation (FEP) theoryis based on statistical mechanics, and computes the change going from state A to state B. Although this method has the correct physical basis, the calculations do not converge unless the difference between the two states is small. Calculations are therefore divided into a number of small steps that do not have to be physically realistic, for example mutating a Cl atom into a F atom (equation 4.3). 

It can be seen from Table 2.2 that the predicted Henry s law constants agree well with the high-quality experimental values showing maximum deviations of 20%. The force field used has been constructed by adjustment to experimental data for the pure components. Since thermophysical data published in the literature often shows a high scattering and is frequently not available for the conditions required in process intensification, MD simulations based on high-precision force fields are an alternative to costly and time-consuming measurements. 

---

There is, of course, a mass of rather direct evidence on orientation at the liquid-vapor interface, much of which is at least implicit in this chapter and in Chapter IV. The methods of statistical mechanics are applicable to the calculation of surface orientation of assymmetric molecules, usually by introducing an angular dependence to the inter-molecular potential function (see Refs. 67, 68, 77 as examples). Widom has applied a mean-held approximation to a lattice model to predict the tendency of AB molecules to adsorb and orient perpendicular to the interface between phases of AA and BB [78]. In the case of water, a molecular dynamics calculation concluded that the surface dipole density corresponded to a tendency for surface-OH groups to point toward the vapor phase [79]. 

The theoretical treatments of Section III-2B have been used to calculate interfacial tensions of solutions using suitable interaction potential functions. Thus Gubbins and co-workers [88] report a molecular dynamics calculation of the surface tension of a solution of A and B molecules obeying Eq. III-46 with o,bb/ o,aa = 0.4 and... 

It was noted in connection with Eq. III-56 that molecular dynamics calculations can be made for a liquid mixture of rare gas-like atoms to obtain surface tension versus composition. The same calculation also gives the variation of density for each species across the interface [88], as illustrated in Fig. Ill-13b. The density profiles allow a calculation, of course, of the surface excess quantities. 

Molecular dynamics calculations have been made on atomic crystals using a Lennard-Jones potential. These have to be done near the melting point in order for the iterations not to be too lengthy and have yielded density functioi). as one passes through the solid-vapor interface (see Ref. 45). The calculations showed considerable mobility in the surface region, amounting to the presence of a... 

The alternative simulation approaches are based on molecular dynamics calculations. This is conceptually simpler that the Monte Carlo method the equations of motion are solved for a system of A molecules, and periodic boundary conditions are again imposed. This method pennits both the equilibrium and transport properties of the system to be evaluated, essentially by numerically solvmg the equations of motion... 

The summation of pair-wise potentials is a good approximation for molecular dynamics calculations for simple classical many-body problems [27], It has been widely used to simulate hyperthennal energy (>1 eV) atom-surface scattering ... 

Allen M P, Warren M A, Wilson M R, Sauron A and Wiliam S 1996 Molecular dynamics calculation of elastic constants in Gay-Berne nematic liquid crystals J. Chem. Phys. 105 2850-8... 

Andersen H C 1983 RATTLE a velocity version of the SHAKE algorithm for molecular dynamics calculations J. Comput. Phys. 52 24-34... 

To demonstrate the basic ideas of molecular dynamics calculations, we shall first examine its application to adiabatic systems. The theory of vibronic coupling and non-adiabatic effects will then be discussed to define the sorts of processes in which we are interested. The complications added to dynamics calculations by these effects will then be considered. Some details of the mathematical formalism are included in appendices. Finally, examples will be given of direct dynamics studies that show how well the systems of interest can at present be treated. 

H. C. Andersen. Rattle A velocity version of the Shake algorithm for molecular dynamics calculations. J. Comp. Phys., 52 24-34, 1983. 

Prepare a molecii le for a molecii lar dynamics sim illation. If the forces on atoms are too large, th e in legralion algorithm may-fail during a molecular dynamics calculation. ... 

HyperChem uses th e ril 31 water m odel for solvation. You can place th e solute in a box of T1P3P water m oleeules an d impose periodic boun dary eon dition s. You may then turn off the boundary conditions for specific geometry optimi/.aiion or molecular dynamics calculations. However, th is produces undesirable edge effects at the solvent-vacuum interface. 

In a molecular dynamics calculation, you can add a term to adjust the velocities, keeping the molecular system near a desired temperature. During a constant temperature simulation, velocities are scaled at each time step. This couples the system to a simulated heat bath at Tq, with a temperature relaxation time of "r. The velocities arc scaled bv a factor X. where... 

You can include geometric restraints—for interatomic distances, bond angles, and torsion angles—in any molecular dynamics calculation or geometry optim i/.ation. Here are some applications of restrain ts ... 

Molecular dynamics calculations can automatically average and save these values ... 

Toxvaerd S 1990. Molecular Dynamics Calculation of the Equation of State of Alkanes. Journal of Chemical Physics 93 4290-4295. 

There are many algorithms for integrating the equations of motion using finite difference methods, several of which are commonly used in molecular dynamics calculations. All algorithms assume that the positions and dynamic properties (velocities, accelerations, etc.) can be approximated as Taylor series expansions ... 

Woodcock L V 1971. Isothermal Molecular Dynamics Calculations for Liquid Salts. Chemical Phi, Letters 10 257-261. 

For fluids, this is computed by a statistical sampling technique, such as Monte Carlo or molecular dynamics calculations. There are a number of concerns that must be addressed in setting up these calculations, such as... 

For ah initio, semiempirical, or molecular dynamics calculations, the amount of CPU time necessary is generally the factor of greatest concern to researchers. For very large molecules, memory use is of concern for molecular mechanics... 

Mesoscale simulations model a material as a collection of units, called beads. Each bead might represent a substructure, molecule, monomer, micelle, micro-crystalline domain, solid particle, or an arbitrary region of a fluid. Multiple beads might be connected, typically by a harmonic potential, in order to model a polymer. A simulation is then conducted in which there is an interaction potential between beads and sometimes dynamical equations of motion. This is very hard to do with extremely large molecular dynamics calculations because they would have to be very accurate to correctly reflect the small free energy differences between microstates. There are algorithms for determining an appropriate bead size from molecular dynamics and Monte Carlo simulations. 

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

Chem3D uses a MM2 force field that has been extended to cover the full periodic table with the exception of the /block elements. Unknown parameters will be estimated by the program and a message generated to inform the user of this. MM2 can be used for both energy minimization and molecular dynamics calculations. The user can add custom atom types or alter the parameters used... 

NWChem (we tested Version 3.2.1) is a program for ah initio, band-structure, molecular mechanics, and molecular dynamics calculations. The DFT band-structure capability is still under development and was not included in the Linux version tested. NWChem is unique in that it was designed from scratch for efficient parallel execution. The user agreement is more restrictive than most, apparently because the code is still under active development. At the time of this book s publication, limited support was available for users outside of the EMSL facility. 

PLS (partial least-squares) algorithm used for 3D QSAR calculations PM3 (parameterization method three) a semiempirical method PMF (potential of mean force) a solvation method for molecular dynamics calculations... 

Molecular dynamics calculations use equations 25-27. HyperChem integrates equations 26 and 27 to describe the motions of atoms. In the absence of temperature regulation, there are no external sources or depositories of energy. That is, no other energy terms exist in the Hamiltonian, and the total energy of the system is constant. 

You can often use experimental data, such as Nuclear Overhauser Effect (NOE) signals from 2D NMR studies, as restraints. NOE signals give distances between pairs of hydrogens in a molecule. Use these distances to limit distances during a molecular mechanics geometry optimization or molecular dynamics calculation. Information on dihedral angles, deduced from NMR, can also limit a conformational search. 

The setup of these calculations is very similar for both quantum and molecular mechanics. In practice, molecular dynamics calculation s using the nl) initio and semi-empirical quantum mechanical SCFmethods are limited to relatively small systems. Each time step requires a complete calculation of the wave function and the forces. 

HyperChem always computes the electronic properties for the molecule as the last step of a geometry optimization or molecular dynamics calculation. However, if you would like to perform a configuration interaction calculation at the optimized geometry, an additional single point calculation is required with the Cl option being turned on. 

In order to conserve the total energy in molecular dynamics calculations using semi-empirical methods, the gradient needs to be very accurate. Although the gradient is calculated analytically, it is a function of wavefunction, so its accuracy depends on that of the wavefunction. Tests for CH4 show that the convergence limit needs to be at most le-6 for CNDO and INDO and le-7 for MINDO/3, MNDO, AMI, and PM3 for accurate energy conservation. ZINDO/S is not suitable for molecular dynamics calculations. 

---