It is common to treat the vibration-rotation problem in methane using the Nielsen transformed-Hamiltonian approach [10, 14, 60, 65]. Some of the terms in such a transformed Hamiltonian are obtained by considering explicitly various orders of the appropriate contact transformations. Other terms are introduced phenomenologically. Because of the variety of notation used and the different methods of introducing new interaction terms, relationships between the parameters determined in various recent vibration-rotation spectral analyses [43-45] of CH4 are not all known with certainty. But in all cases, the form of the interaction terms in the Hamiltonian operator is constrained by the group theoretical requirement that they be of species A1 in the point group Td.
Consider first the form of the centrifugal distortion operator, which consists of terms quartic in the molecule-fixed components of the total angular momentum operator. Since Jx , Jy , Jz transform as F1, we are interested in the species of F14. However, we need actually only consider species belonging to the symmetrized fourth power [F14]. Operators belonging to other species contained in F1 × F1 × F1 × F1 can be written in terms of products cubic and lower in the angular momentum components, by making use of the commutation relations
$$J_i J_j - J_j J_i = -i\hbar \, e_{ijk} J_k ~ . $$
(eq. 62)
Table 9 indicates that [F14] contains the representation A1 only twice. Thus, there are for methane only two linearly independent fourth power expressions in the total angular momentum components, associated with two independent quartic centrifugal distortion constants. The quartic centrifugal distortion operator is expressed in many different notations in the literature. One which parallels symmetric top notation somewhat is the following (suppressing ℏ's)
$$\begin{eqnarray*} {\cal H}_{\rm cent. \, dist.} = -D_J \mbox{$J$}^4 &-& {\textstyle{1\over2}} \, D_t \left[ 3\mbox{$J$}^4 - 30\mbox{$J$}^2 J_z^2 + 35 J_z^4 - 6\mbox{$J$}^2 + 25J_z^2 \right .\\ &+& \left .{\textstyle{5\over2}} (J_x+iJ_y)^4 + {\textstyle{5\over2}} (J_x-iJ_y)^4 \right] ~ . \end{eqnarray*}$$
(eq. 63)
The construction of this operator from first principles is beyond the scope of this work, but it can be verified using (eq. 5) and Table 3, that it is of species A1 for any values of DJ and Dt . We note also that the second term contains the operator J2Jz2, associated in a symmetric top with the constant DJK, and the operator Jz4, associated in a symmetric top with the constant DK.
Jx , Jy , Jz can also be considered to transform according to the irreducible representation Dg(1) of the continuous three-dimensional rotation-reflection group. The symmetrized fourth power [(Dg(1))4] contains the species Dg(0) + Dg(2) + Dg(4). It can be shown [14] that the first term in (eq. 63) transforms as Dg(0) under the three-dimensional rotation-reflection group, while the second term corresponds to a linear combination of terms transforming according to Dg(4). If we represent functions transforming according to Dg(4) by fm(4), where -4 ≤ m ≤ +4, the second term in (eq. 63) is proportional to [14]
$$(2 \cdot 5 \cdot 7)^{1/2} \, f_0^{(4)} + 5 \left[f_{+4}^{(4)}+f_{-4}^{(4)} \right] ~ . $$
(eq. 64)
It is common [10, 14, 65] to handle higher-order vibration-rotation operators for υs = 1 states, s = 3 or 4, rather formally as follows. Consider a set of (2k1 + 1) vibrational operators fm1(k1) which transform according to the irreducible representation D(k1) of the continuous three-dimensional rotation group. These vibrational operators are taken to be functions of the vibrational coordinates Qsx , Qsy , Qsz , their conjugate linear momenta, and the components Lsx , Lsy , Lsz of the vibrational angular momentum. Consider also a set of (2k2 + 1) rotational operators gm2(k2) which transform according to D(k2) and are taken to be functions of Jx , Jy , Jz . It is then possible to couple these operators using the vector coupling techniques to obtain vibration-rotation operators
$$h_m^{(k)} = \sum_{m_1 m_2} \, f_{m_1}^{(k_1)} \, g_{m_2}^{(k_2)} (k_1 k_2 m_1 m_2 \, | \, k_1 k_2 k \, m ) $$
(eq. 65)
which transform according to D(k). Matrix elements of these coupled vibration-rotation operators can be obtained formally using the results of irreducible tensor theory.
It is possible to construct all operators involving components of the nuclear spin angular momentum raised to a given power by considering only molecule-fixed components of these spin operators, since the molecule-fixed components are automatically invariant with respect to overall rotations in the laboratory-fixed axis system. However, the spin basis set functions are characterized by projections quantized along the laboratory-fixed Z axis, so that it is necessary in calculating matrix elements to use the analog of (eq. 35) to write the spin operators in terms of laboratory-fixed components and elements of the direction cosine matrix.
The application of symmetry considerations to the resulting expressions, which involve tensors in the rotational operators and tensors in the nuclear spin operators, is particularly complicated because care must be taken to insure that final coupled tensor operators transform as Dg(0) under overall rotations in the laboratory, and according to A1 under the point-group symmetry operations. Loosely speaking, the laboratory-fixed "projections" of the rotational tensors must be correctly coupled to the laboratory-fixed "projections" of the spin tensors, and the molecule-fixed "projections" of the rotational tensors must be correctly coupled to the molecule-fixed "projections" of the spin tensors.
The formalism of irreducible tensorial sets displays its greatest power in systematizing the evaluation of matrix elements. If we consider two wave functions Ψm1(j1) and Ψm2(j2), and an operator Om3(j3), which transform, respectively, according to the irreducible representation D(j1), D(j2), and D(j3) of the continuous three-dimensional rotation group, then matrix elements can be written [22, 62, 64] in the form
$$\langle\Psi_{m_1}^{(j_1)}| \, O_{m_3}^{(j_3)} \, |\Psi_{m_2}^{(j_2)}\rangle = \langle j_1| |O^{(j_3)}| |j_2\rangle \, (j_2 j_3 m_2 m_3 \, | \, j_2 j_3 j_1 m_1) ~ . $$
(eq. 66)
The first quantity on the right is a so-called reduced matrix element; its value does not depend on m1, m2, or m3 and must in general be calculated for the specific molecular problem under consideration. The second factor is a vector coupling coefficient; its value is determined by symmetry and can be tabulated once and for all, or conveniently manipulated algebraically. Because of (eq. 66), (2j1 + 1) (2j2 + 1) (2j3 + 1) matrix elements can be related to a single unknown parameter, achieving thereby an enormous simplification and reduction of labor.
Equation (66) gives an example of a matrix element in an uncoupled basis set. It is also written in terms of 3 - j symbols in the literature. Extensive formalism has been developed to yield matrix elements of coupled operators in coupled basis sets, which requires the introduction of 6 - j and 9 - j symbols, etc. There are numerous different phase conventions used in the literature in the definitions of 3 - j symbols and related quantities. Users of irreducible tensor formalism must be absolutely certain that their own applications are internally self-consistent.