Equilibration time, molecular dynamics

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. 

In many molecular dynamics simulations, equilibration is a separate step that precedes data collection. Equilibration is generally necessary to avoid introducing artifacts during the heating step and to ensure that the trajectory is actually simulating equilibrium properties. The period required for equilibration depends on the property of interest and the molecular system. It may take about 100 ps for the system to approach equilibrium, but some properties are fairly stable after 10-20 ps. Suggested times range from 5 ps to nearly 100 ps for medium-sized proteins. 

Molecular dynamics simulations of proteins often begin with a known structure (such as an X-ray diffraction structure) that you want to maintain during equilibration. Since the solvent may contain high energy hot spots, equilibration of the protein and solvent at the same time can change the protein conformation. To avoid this, select only the water molecules and run a molecular dynamics equilibration. This relaxes the water while fixing the protein structure. Then deselect the water and equilibrate the whole system. 

Successful molecular dynamics simulations should have a fairly stable trajectory. Instability and lack of equilibration can result from a large time step, treatment of long-range cutoffs, or unrealistic coupling to a temperature bath. 

The simplest method that keeps the temperature of a system constant during an MD simulation is to rescale the velocities at each time step by a factor of (To/T) -, where T is the current instantaneous temperature [defined in Eq. (24)] and Tq is the desired temperamre. This method is commonly used in the equilibration phase of many MD simulations and has also been suggested as a means of performing constant temperature molecular dynamics [22]. A further refinement of the velocity-rescaling approach was proposed by Berendsen et al. [24], who used velocity rescaling to couple the system to a heat bath at a temperature Tq. Since heat coupling has a characteristic relaxation time, each velocity V is scaled by a factor X, defined as... 

Recently, detailed molecular pictures of the interfacial structure on the time and distance scales of the ion-crossing event, as well as of ion transfer dynamics, have been provided by Benjamin s molecular dynamics computer simulations [71, 75, 128, 136]. The system studied [71, 75, 136] included 343 water molecules and 108 1,2-dichloroethane molecules, which were separately equilibrated in two liquid slabs, and then brought into contact to form a box about 4 nm long and of cross-section 2.17 nmx2.17 nm. In a previous study [128], the dynamics of ion transfer were studied in a system including 256 polar and 256 nonpolar diatomic molecules. Solvent-solvent and ion-solvent interactions were described with standard potential functions, comprising coulombic and Lennard-Jones 6-12 pairwise potentials for electrostatic and nonbonded interactions, respectively. While in the first study [128] the intramolecular bond vibration of both polar and nonpolar solvent molecules was modeled as a harmonic oscillator, the next studies [71,75,136] used a more advanced model [137] for water and a four-atom model, with a united atom for each of two... 

---