Lattice energy calculation molecular dynamics

Several issues remain to be addressed. The effect of the mutual penetration of the electron distributions should be analyzed, while the use of theoretical densities on isolated molecules does not take into account the induced polarization of the molecular charge distribution in a crystal. In the calculations by Coombes et al. (1996), the effect of electron correlation on the isolated molecule density is approximately accounted for by a scaling of the electrostatic contributions by a factor of 0.9. Some of these effects are in opposite directions and may roughly cancel. As pointed out by Price and coworkers, lattice energy calculations based on the average static structure ignore the dynamical aspects of the molecular crystal. However, the necessity to include electrostatic interactions in lattice energy calculations of molecular crystals is evident and has been established unequivocally. 

Table 1. Comparison of bulk vibrational energies per cation calculated using lattice d5mamics and molecular dynamics.

In principle, we could find the minimum-energy crystal lattice from electronic structure calculations, determine the appropriate A-body interaction potential in the presence of lattice defects, and use molecular dynamics methods to calculate ab initio dynamic macroscale material properties. Some of the problems associated with this approach are considered by Wallace [1]. Because of these problems it is useful to establish a bridge between the micro-... 

These two methods are different and are usually employed to calculate different properties. Molecular dynamics has a time-dependent component, and is better at calculating transport properties, such as viscosity, heat conductivity, and difftisivity. Monte Carlo methods do not contain information on kinetic energy. It is used more in the lattice model of polymers, protein stmcture conformation, and in the Gibbs ensemble for phase equilibrium. 

The lattice energy of a crystal of known structure (atomic positions) is thus calculated by compiling all possible distances between pairs of atoms in different molecules. The method of atom-atom potentials has been employed to investigate phenomena pertaining to static as well as dynamic lattices and the subject has been reviewed by Kitaigorodsky (1973) as well as by Ramdas Thomas (1980). Typical of the problems that have been investigated by this method are defects and planar faults, phase transitions and molecular rotation in crystals. 

What is next Several examples were given of modem experimental electrochemical techniques used to characterize electrode-electrolyte interactions. However, we did not mention theoretical methods used for the same purpose. Computer simulations of the dynamic processes occurring in the double layer are found abundantly in the literature of electrochemistry. Examples of topics explored in this area are investigation of lateral adsorbate-adsorbate interactions by the formulation of lattice-gas models and their solution by analytical and numerical techniques (Monte Carlo simulations) [Fig. 6.107(a)] determination of potential-energy curves for metal-ion and lateral-lateral interaction by quantum-chemical studies [Fig. 6.107(b)] and calculation of the electrostatic field and potential drop across an electric double layer by molecular dynamic simulations [Fig. 6.107(c)]. 

A molecular dynamics simulation in conjunction with experimental evidence was used to elucidate the nature of the interactions between polymer materials and CNTs [239]. Computational time was reduced by representing CNTs as a force field. The calculations indicated an extremely strong noncovalent binding energy. Furthermore, the correlation between the chirality of the nanotubes and mapping of the polymer on to the lattice was discussed [239]. 

As molecular packing calculations involve just simple lattice energy minimizations another set of tests have focused on the finite temperature effects. For this purpose, Sorescu et al. [112] have performed isothermal-isobaric Monte Carlo and molecular dynamics simulations in the temperature range 4.2-325 K, at ambient pressure. It was found that the calculated crystal structures at 300 K were in outstanding agreement with experiment within 2% for lattice dimensions and almost no rotational and translational disorder of the molecules in the unit cell. Moreover, the space group symmetry was maintained throughout the simulations. Finally, the calculated expansion coefficients were determined to be in reasonable accord with experiment. 

---