Potential energy surfaces force-constant matrix

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Bom-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and force constant fit reasonable well those obtained by other approaches. 

The saddle point on a three-dimensional potential-energy surface, characterized by one negative force constant in the harmonic force constant matrix. 

A normal-mode representation of the Hamiltonian for the reduced system involves the diagonalization of the projected force constant matrix, which in turn generates a reduced-dimension potential-energy surface in terms of the mass-weighted coordinates of the reaction path [64] ... 

Because of their importance to nucleation kinetics, there have been a number of attempts to calculate free energies of formation of clusters theoretically. The most important approaches for the current discussion are harmonic models, " Monte Carlo studies, and molecular dynamics calcula-tions. In the harmonic model the cluster is assumed to be composed of constituent atoms with harmonic intermolecular forces. The most recent calculations, which use the harmonic model, have taken the geometries of the clusters to be those determined by the minimum in the two-body additive Lennard-Jones potential surface. The oscillator frequencies have been obtained by diagonalizing the Lennard-Jones force constant matrix. In the harmonic model the translational and rotational modes of the clusters are treated classically, and the vibrational modes are treated quantum mechanically. The harmonic models work best at low temjjeratures where anharmonic-ity effects are least important and the system is dominated by a single structure. 

The curvature of the potential energy surface in the direction Qesk at Res, is measured by the force constant, Kk = Kq + [46], The diagonal matrix ele-... 

In summary, a high diversity of potential models for molecular mechanics calculations of zeolites hitherto exists. From the theoretical point of view, an appropriate force field should be able to predict structures and vibrations with similar accuracy. On the other hand, the structure of a system under study is determined by the energy minimum, whereas normal modes are dependent on the curvature (second derivative) of the potential energy surface. Consequently, force fields obviously successful in predicting structural features might not automatically be appropriate for simulating vibrational spectra. The only way to overcome this difficulty is to include experimental spectroscopic data into the parametrization process [60]. Alternatively, besides structures and energies a matrix of force constants obtained in quantum mechanical calculations can be included into the quantum mechanical data base used to tune the parameters of the potential function [51]. 

The second-derivative matrix of a potential energy (hyperf surface. Synonymous with force constant matrix. 

---