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,
38
Chem. Phys.
(1990),
145,
pp. 427.
Link
,
235
Chem. Phys.
(1988),
123,
pp. 187.
Link
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,
ˆU,
ˆU=m∑i(-12∂2∂q2i+ωi22qi2) | (10.71) |
and a perturbed part, ˆV:
ˆV=16m∑ijk=1ηijkqiqjqk+124m∑ijkl=1ηijklqiqjqkql | (10.72) |
One can then apply second-order perturbation theory to get the ith excited state energy:
E(i)=ˆU(i)+⟨Ψ(i)|ˆV|Ψ(i)⟩+∑j≠i|⟨Ψ(i)|ˆV|Ψ(j)⟩|2ˆU(i)-ˆU(j) | (10.73) |
The denominator in Eq. (10.73) 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.
929
Rev. Mod. Phys.
(1951),
23,
pp. 90.
Link
,
38
Chem. Phys.
(1990),
145,
pp. 427.
Link
,
235
Chem. Phys.
(1988),
123,
pp. 187.
Link
An
alternative solution has been proposed by Barone,
73
J. Chem. Phys.
(2005),
122,
pp. 014108.
Link
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
EVPT2i=∑jωj(nj+1/2)+∑i≤jxij(ni+1/2)(nj+1/2) | (10.74) |
where, if rotational coupling is ignored, the anharmonic constants xij are given by
xij=14ωiωj(ηiijj-m∑kηiikηjjkω2k+m∑k2(ω2i+ω2j-ω2k)η2ijk[(ωi+ωj)2-ω2k][(ωi-ωj)2-ω2k]) | (10.75) |