Analytical expressions energy gradients, calculation

Free energy second derivatives are mainly used to analyse the nature of stationary points on the PES, and to compute harmonic force constants and vibrational frequencies to perform such calculations in solution, one needs analytical expressions for Qa second derivatives with respect to nuclear displacements (the alternative of using numerical differentiation of gradients is far too much expensive except for very small molecules). 

Interatomic forces for molecular dynamics can be calculated from analytical expressions for the gradients of the total energy at the atomic positions. All integral tables and repulsive parameters have been determined for carbon previously. For a review of the DFTB methodology and its application, see ref. 

To show the principles involved in finding an analytical gradient expression consider an HF-LCAO calculation where the electronic energy comes to... 

Analytical calculation of the gradient is the key to efficient geometry optimization. The SCF MO energy expression is Eq. (14.44), and its derivatives with respect to the nuclear coordinates would seem to involve the derivatives of the and (rs tu) integrals... 

---