Energy derivatives, calculating

Quantum mechanical calculations are restricted to systems with relatively small numbers of atoms, and so storing the Hessian matrix is not a problem. As the energy calculation is often the most time-consuming part of the calculation, it is desirable that the minimisation method chosen takes as few steps as possible to reach the minimum. For many levels of quantum mechanics theory analytical first derivatives are available. However, analytical second derivatives are only available for a few levels of theory and can be expensive to compute. The quasi-Newton methods are thus particularly popular for quantum mechanical calculations. 

The most extensive calculations of the electronic structure of fullerenes so far have been done for Ceo- Representative results for the energy levels of the free Ceo molecule are shown in Fig. 5(a) [60]. Because of the molecular nature of solid C o, the electronic structure for the solid phase is expected to be closely related to that of the free molecule [61]. An LDA calculation for the crystalline phase is shown in Fig. 5(b) for the energy bands derived from the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for Cgo, and the band gap between the LUMO and HOMO-derived energy bands is shown on the figure. The LDA calculations are one-electron treatments which tend to underestimate the actual bandgap. Nevertheless, such calculations are widely used in the fullerene literature to provide physical insights about many of the physical properties. 

Computing the vibrational frequencies of molecules resulting from interatomic motion within the molecule. Frequencies depend on the second derivative of the energy with respect to atomic structure, and frequency calculations may also predict other properties which depend on second derivatives. Frequency calculations are not possible or practical for all computational chemistry methods. 

Fig. 20. A. Conformation of the Valinomycin-cation complex derived for solution using a combination of proton magnetic resonance data and conformational energy calculations. This structure agrees within tenths of an Angstrom with the crystal structure subsequently determined (100) and shown in Fig. 21. Reproduced with permission from Ref.99).

Fig. 21. Crystal structure of the Valinomycin-K+ complex. Reproduced with permission from Ref.100). This crystal structure confirmed within tenths of an Angstrom the structure derived previously in solution 97 98) and by means of conformational energy calculations...

All of the free energy calculations to this point have involved the standard free energy change, AG°. It is possible, however, to write a general relation for the free energy change, AG, valid under any conditions. This relation is a relatively simple one, but we will not attempt to derive it. It tells us that... 

Values for the partial charges of atoms can be derived from quantum mechanical calculations, from the molecular dipole moments and from rotation-vibration spectra. However, often they are not well known. If the contribution of the Coulomb energy cannot be calculated precisely, no reliable lattice energy calculations are possible. 

So far, we have only considered perturbations in the forward direction. As in conventional free energy calculations, powerful relations can be derived if forward and backward perturbations are combined. With the free energy being a state function, we can reverse the path direction. This leads to several useful relations derived originally by Crooks [16, 18, 19], In particular, we obtain immediately... 

From (9.27), we see that this approach will work nicely if the variance is always small Taylor s theorem with remainder tells us that the error of the first-derivative - mean-field - contribution is proportional to the second derivative evaluated at an intermediate A. That second derivative can be identified with the variance as in (9.27). If that variance is never large, then this approach should be particularly effective. For further discussion, see Chap. 4 on thermodynamic integration, and Chap. 6 on error analysis in free energy calculations. 

A first step toward quantum mechanical approximations for free energy calculations was made by Wigner and Kirkwood. A clear derivation of their method is given by Landau and Lifshitz [43]. They employ a plane-wave expansion to compute approximate canonical partition functions which then generate free energy models. The method produces an expansion of the free energy in powers of h. Here we just quote several of the results of their derivation. 

To exploit the concept of PMF to represent solvent in free energy calculations, practical approximations must be constructed. A common approach is to treat the two components Z H/"P(X) and Z lYelec(X) separately. Approximations for the nonpolar term are usually derived from geometric considerations, as in scaled particle theory, for example [62], The electrostatic contribution is usually derived from continuum electrostatics. We consider these two contributions in turn. 

AMI semi-empirical and B3LYP/6-31G(d)/AMl density functional theory (DFT) computational studies were performed with the purpose of determining which variously substituted 1,3,4-oxadiazoles would participate in Diels-Alder reactions as dienes and under what conditions. Also, bond orders for 1,3,4-oxadiazole and its 2,5-diacetyl, 2,5-dimethyl, 2,5-di(trifluoromethyl), and 2,5-di(methoxycarbonyl) derivatives were calculated <1998JMT153>. The AMI method was also used to evaluate the electronic properties of 2,5-bis[5-(4,5,6,7-tetrahydrobenzo[A thien-2-yl)thien-2-yl]-l,3,4-oxadiazole 8. The experimentally determined redox potentials were compared with the calculated highest occupied molecular orbital/lowest unoccupied molecular orbital (HOMO/LUMO) energies. The performance of the available parameters from AMI was verified with other semi-empirical calculations (PM3, MNDO) as well as by ab initio methods <1998CEJ2211>. 

---