Molecular dynamics constant energy

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. 

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. 

For free dynam ics wdierc the constant temperature check box in the Molecular Dynam ics dialog box is not checked, the total energy K l OT should remain constant,, A fhictuation of the total... 

Just as one may wish to specify the temperature in a molecular dynamics simulation, so may be desired to maintain the system at a constant pressure. This enables the behavior of the system to be explored as a function of the pressure, enabling one to study phenomer such as the onset of pressure-induced phase transitions. Many experimental measuremen are made under conditions of constant temperature and pressure, and so simulations in tl isothermal-isobaric ensemble are most directly relevant to experimental data. Certai structural rearrangements may be achieved more easily in an isobaric simulation than i a simulation at constant volume. Constant pressure conditions may also be importai when the number of particles in the system changes (as in some of the test particle methoc for calculating free energies and chemical potentials see Section 8.9). 

The pressure often fluctuates much more than quantities such as the total energy in constant NVE molecular dynamics simulation. This is as expected because the pressure related to the virial, which is obtained as the product of the positions and the derivativ of the potential energy function. This product, rijdf rij)/drij, changes more quickly with than does the internal energy, hence the greater fluctuation in the pressure. 

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. 

Once HyperChem calculates potential energy, it can obtain all of the forces on the nuclei at negligible additional expense. This allows for rapid optimization of equilibrium and transition-state geometries and the possibility of computing force constants, vibrational modes, and molecular dynamics trajectories. 

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. 

Because of limitations in computer power and time, it is frequently impractical to run a constant energy molecular dynamics simulation. Several approximations to the energy (usually to the potential energy) are possible, which require modifying the Hamilto-... 

If there is no external temperature control (using a simulated constant temperature bath), molecular dynamics simulations are constant energy. 

For a conformation in a relatively deep local minimum, a room temperature molecular dynamics simulation may not overcome the barrier and search other regions of conformational space in reasonable computing time. To overcome barriers, many conformational searches use elevated temperatures (600-1200 K) at constant energy. To search conformational space adequately, run simulations of 0.5-1.0 ps each at high temperature and save the molecular structures after each simulation. Alternatively, take a snapshot of a simulation at about one picosecond intervals to store the structure. Run a geometry optimization on each structure and compare structures to determine unique low-energy conformations. 

For constant temperature dynamics where the constant temperature check box in the Molecular Dynamics Options dialog box is checked, the energy will not remain constant but will fluctuate as energy is exchanged with the bath. The temperature, depending on the value set for the relaxation constant, will approach con-stan cy. 

A molecular dynamics force field is a convenient compilation of these data (see Chapter 2). The data may be used in a much simplified fonn (e.g., in the case of metric matrix distance geometry, all data are converted into lower and upper bounds on interatomic distances, which all have the same weight). Similar to the use of energy parameters in X-ray crystallography, the parameters need not reflect the dynamic behavior of the molecule. The force constants are chosen to avoid distortions of the molecule when experimental restraints are applied. Thus, the force constants on bond angle and planarity are a factor of 10-100 higher than in standard molecular dynamics force fields. Likewise, a detailed description of electrostatic and van der Waals interactions is not necessary and may not even be beneficial in calculating NMR strucmres. 

---