10.9 Anharmonic Vibrational Frequencies 10.9.2 Vibration Configuration Interaction Theory 10.9.4 Transition-Optimized Shifted Hermite Theory

10.9.3 Vibrational Perturbation Theory

(February 4, 2022)

Vibrational perturbation theory has been historically popular for calculating molecular spectroscopy. Nevertheless, it is notorious for the inability of dealing with resonance cases. In addition, the non-standard formulas for various symmetries of molecules forces the users to modify inputs on a case-by-case basis,^{35, 784, 209} which narrows the accessibility of this method. VPT applies perturbation treatments on the same Hamiltonian as in Eq. (4.1), but divides it into an unperturbed part, $\hat{U}$ ,

\hat{U}=\sum_{i}^{m}\left(-\frac{1}{2}\frac{\partial^{2}}{\partial q_{i}^{2}}+% \frac{{{{\omega}_{i}}}^{2}}{2}{{q_{i}}}^{2}\right)

(10.49)

and a perturbed part, $\hat{V}$ :

\hat{V}=\frac{1}{6}\sum_{{ijk}=1}^{m}{{\eta}_{{ijk}}}{q_{i}}{q_{j}}{q_{k}}+% \frac{1}{24}\sum_{{ijkl}=1}^{m}{{\eta}_{{ijkl}}}{q_{i}}{q_{j}}{q_{k}}{q_{l}}

(10.50)

One can then apply second-order perturbation theory to get the $i$ th excited state energy:

E^{(i)}=\hat{U}^{(i)}+\langle\Psi^{(i)}|\hat{V}|\Psi^{(i)}\rangle+\sum_{j\neq i% }\frac{|\langle\Psi^{(i)}|\hat{V}|\Psi^{(j)}\rangle|^{2}}{\hat{U}^{(i)}-\hat{U% }^{(j)}}

(10.51)

The denominator in Eq. (10.51) can be zero either because of symmetry or accidental degeneracy. Various solutions, which depend on the type of degeneracy that occurs, have been developed which ignore the zero-denominator elements from the Hamiltonian.^{820, 35, 784, 209} An alternative solution has been proposed by Barone,⁶⁵ which can be applied to all molecules by changing the masses of one or more nuclei in degenerate cases. The disadvantage of this method is that it will break the degeneracy which results in fundamental frequencies no longer retaining their correct symmetry. He proposed

E_{i}^{\mathrm{VPT2}}=\sum_{j}\omega_{j}(n_{j}+1/2)+\sum_{i\leq j}x_{ij}(n_{i}% +1/2)(n_{j}+1/2)

(10.52)

where, if rotational coupling is ignored, the anharmonic constants $x_{ij}$ are given by

x_{ij}=\frac{1}{4\omega_{i}\omega_{j}}\left(\eta_{iijj}-\sum_{k}^{m}\frac{\eta% _{iik}\eta_{jjk}}{\omega_{k}^{2}}+\sum_{k}^{m}\frac{2(\omega_{i}^{2}+\omega_{j% }^{2}-\omega_{k}^{2})\eta_{ijk}^{2}}{\left[(\omega_{i}+\omega_{j})^{2}-\omega_% {k}^{2}\right]\left[(\omega_{i}-\omega_{j})^{2}-\omega_{k}^{2}\right]}\right)

(10.53)

Search Results

10.9.3 Vibrational Perturbation Theory