Molecular dynamics periodic boundary condition

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

A typical molecular dynamics simulation comprises an equflibration and a production phase. The former is necessary, as the name imphes, to ensure that the system is in equilibrium before data acquisition starts. It is useful to check the time evolution of several simulation parameters such as temperature (which is directly connected to the kinetic energy), potential energy, total energy, density (when periodic boundary conditions with constant pressure are apphed), and their root-mean-square deviations. Having these and other variables constant at the end of the equilibration phase is the prerequisite for the statistically meaningful sampling of data in the following production phase. 

Often yon need to add solvent molecules to a solute before running a molecular dynamics simiilatmn (see also Solvation and Periodic Boundary Conditions" on page 62). In HyperChem, choose Periodic Box on the Setup m en ii to enclose a soln te in a periodic box filled appropriately with TIP3P models of water inole-cii les. 

Note MM-i- is derived from the public domain code developed by Dr. Norm an Allinger, referred to as M.M2( 1977), and distributed by the Quantum Chemistry Program Exchange (QCPE). The code for MM-t is not derived from Dr. Allin ger s present version of code, which IS trademarked MM2 . Specifically. QCMPOlO was used as a starting point Ibr HyperChem MM-t code. The code was extensively modified and extended over several years to include molecular dynamics, switching functuins for cubic stretch terms, periodic boundary conditions, superimposed restraints, a default (additional) parameter scheme, and so on. 

The first molecular dynamics simulations of a lipid bilayer which used an explicit representation of all the molecules was performed by van der Ploeg and Berendsen in 1982 [van dei Ploeg and Berendsen 1982]. Their simulation contained 32 decanoate molecules arranged in two layers of sixteen molecules each. Periodic boundary conditions were employed and a xmited atom force potential was used to model the interactions. The head groups were restrained using a harmonic potential of the form ... 

The tests in the two previous paragraphs are often used because they are easy to perform. They are, however, limited due to their neglect of intermolecular interactions. Testing the effect of intennolecular interactions requires much more intensive simulations. These would be simulations of the bulk materials, which include many polymer strands and often periodic boundary conditions. Such a bulk system can then be simulated with molecular dynamics, Monte Carlo, or simulated annealing methods to examine the tendency to form crystalline phases. 

HyperChem uses the TIP3P water model for solvation.You can place the solute in a box of TIP3P water molecules and impose periodic boundary conditions. You may then turn off the boundary conditions for specific geometry optimization or molecular dynamics calculations. However, this produces undesirable edge effects at the solvent-vacuum interface. 

In a normal molecular dynamics simulation with repeating boundary conditions (i.e., periodic boundary condition), the volume is held fixed, whereas at constant pressure the volume of the system must fluemate. In some simulation cases, such as simulations dealing with membranes, it is more advantageous to use the constant-pressure MD than the regular MD. Various schemes for prescribing the pressure of a molecular dynamics simulation have also been proposed and applied [23,24,28,29]. In all of these approaches it is inevitable that the system box must change its volume. 

Periodic boundary conditions, Monte Carlo heat flow simulation, nonequilibrium molecular dynamics, 79—81 Periodic-orbit dividing surface (PODS) geometric transition state theory, 196-201 transition state trajectory, 202-213 Perturbation theory, transition state trajectory, deterministically moving manifolds, 224-228... 

As electric fields and potential of molecules can be generated upon distributed p, the second order energies schemes of the SIBFA approach can be directly fueled by the density fitted coefficients. To conclude, an important asset of the GEM approach is the possibility of generating a general framework to perform Periodic Boundary Conditions (PBC) simulations. Indeed, such process can be used for second generation APMM such as SIBFA since PBC methodology has been shown to be a key issue in polarizable molecular dynamics with the efficient PBC implementation [60] of the multipole based AMOEBA force field [61]. 

This potential was developed to ensure that the molecules inside the sphere never escape and maintain a fully solvated system during molecular dynamics. Here, es, Rs, ew and Rw are the van der Waals constants for the solvent and the wall and rj is the distance between the molecule i and the center of the water sphere, Ro is the radius of the sphere. The quantities A, B and Rb are determined by imposing the condition that W and dW/dr, vanish at r, = Ro. The restraining potential W is set to zero for r, < R0. The van der Waals parameters Es, ew, Rs and Rw can also be specifically defined for different solvents. The constants Awaii and Cwan are computed using a well depth of es = ew = 0.1 kcal and the radius of Rs = Rw = 1.25 A. For the other set of simulations, especially for the hydride ion transfer, we applied periodic boundary conditions by using a spherical boundary shell of 10.0 A of TIP3P40 water to cover the edges of the protein. 

The study of the enantioselective hydrosilylation reaction was performed with a series of combined quantum mechanics/molecular mechanics (QM/MM) calculations [26, 30] within the computational scheme of ab initio (AIMD) (Car-Parrinello) [62] molecular dynamics. The AIMD approach has been described in a number of excellent reviews [63-66], AIMD as well as hybrid QM/MM-AIMD calculations [26, 47] were performed with the ab initio molecular dynamics program CPMD [67] based on a pseudopotential framework, a plane wave basis set, and periodic boundary conditions. We have recently developed an interface to the CPMD package in which the coupling with a molecular mechanics force field has been implemented [26, 68],... 

Once the boundary conditions have been implemented, the calculation of solution molecular dynamics proceeds in essentially the same manner as do vacuum calculations. While the total energy and volume in a microcanonical ensemble calculation remain constant, the temperature and pressure need not remain fixed. A variant of the periodic boundary condition calculation method keeps the system pressure constant by adjusting the box length of the primary box at each step by the amount necessary to keep the pressure calculated from the system second virial at a fixed value (46). Such a procedure may be necessary in simulations of processes which involve large volume changes or fluctuations. Techniques are also available, by coupling the system to a Brownian heat bath, for performing simulations directly in the canonical, or constant T,N, and V, ensemble (2,46). 

The computer simulations employed the molecular dynamics technique, in which particles are moved deterministically by integrating their equations of motion. The system size was 864 Lennard-Jones atoms, of which one was the solute (see Table II for potential parameters). There were no solute-solute interactions. Periodic boundary conditions and the minimum image criterion were used (76). The cutoff radius for binary interactions was 3.5 G (see Table II). Potentials were truncated beyond the cutoff. 

---