Next Article in Journal
Clickbait Convolutional Neural Network
Previous Article in Journal
Reinforcement Learning Based Data Self-Destruction Scheme for Secured Data Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetry Adaptation of the Rotation-Vibration Theory for Linear Molecules

1
Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
2
Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Bergische Universität Wuppertal, D-42097 Wuppertal, Germany
*
Author to whom correspondence should be addressed.
Symmetry 2018, 10(5), 137; https://doi.org/10.3390/sym10050137
Submission received: 10 April 2018 / Revised: 23 April 2018 / Accepted: 25 April 2018 / Published: 29 April 2018

Abstract

:
A numerical application of linear-molecule symmetry properties, described by the D h point group, is formulated in terms of lower-order symmetry groups D n h with finite n. Character tables and irreducible representation transformation matrices are presented for D n h groups with arbitrary n-values. These groups can subsequently be used in the construction of symmetry-adapted ro-vibrational basis functions for solving the Schrödinger equations of linear molecules. Their implementation into the symmetrisation procedure based on a set of “reduced” vibrational eigenvalue problems with simplified Hamiltonians is used as a practical example. It is shown how the solutions of these eigenvalue problems can also be extended to include the classification of basis-set functions using , the eigenvalue (in units of ) of the vibrational angular momentum operator L ^ z . This facilitates the symmetry adaptation of the basis set functions in terms of the irreducible representations of D n h . 12 C 2 H 2 is used as an example of a linear molecule of D h point group symmetry to illustrate the symmetrisation procedure of the variational nuclear motion program Theoretical ROVibrational Energies (TROVE).

1. Introduction

The geometrical symmetry of a centrosymmetric linear molecule in its equilibrium geometry is described by the D h point group (see Table 1). While the molecular vibrational states (assuming a totally symmetric singlet electronic state) span the representations of this point group of infinite order, the symmetry properties of the combined rotation-vibration states must satisfy the nuclear-statistics requirements and transform according to the irreducible representations (irreps) of the finite molecular symmetry group
D h ( M ) = { E , ( p ) , E * , ( p ) * }
where, for the centrosymmetric linear molecule A–B–C–…–C–B–A, the permutation operation ( p ) is the simultaneous interchange of the two A nuclei, the two B nuclei, the two C nuclei, etc., E is the identity operation, E * is the spatial inversion operation, which inverts the positions of all particles through the molecular centre of mass, and ( p ) * = ( p ) E * is the permutation-inversion operation [1]. The irreps of D h (M) are given in Table 2 (see also Table A-18 of Ref. [1]). Table 2 presents several alternative notations for the irreducible representations. These alternative notations and multiple names for the same concept are perhaps not quite in agreement with time-honored principles such as Occam’s Razor, but they represent a nice example of the development of spectroscopic notation. We note already here that D h (M) is isomorphic (and, for a triatomic linear molecule A–B–A like CO 2 , identical) to the group customarily called C 2 v (M) (Table A-5 of Ref. [1]), the molecular symmetry group of, for example, the H 2 O molecule, whose equilibrium structure is bent. The molecular symmetry (MS) groups D h (M) and C 2 v (M) are determined by applying the principle of feasibility first introduced by Longuet-Higgins [2] (see also Ref. [1]), and one obviously obtains isomorphic MS groups for all chain molecules A–B–C–…–C–B–A, irrespective of these molecules having linear or bent equilibrium structures. Longuet-Higgins [2] (see also Ref. [1]) further showed that for a so-called rigid non-linear molecule (in this context, a rigid molecule is defined as one whose vibration can be described as oscillations around a single potential energy minimum), the MS group is isomorphic to the point group describing the geometrical symmetry at the equilibrium geometry. H 2 O is a rigid non-linear molecule, whose geometrical symmetry at equilibrium is described by the C 2 v point group which is indeed isomorphic to the MS group C 2 v (M). For the rigid linear molecule CO 2 , however, as already mentioned, the geometrical symmetry at equilibrium is described by the infinite-order point group D h which obviously is not isomorphic to the MS group D h (M) = C 2 v (M) of order four.
One can argue that the MS group as defined by Longuet-Higgins [2] (see also Ref. [1]) provides the simplest symmetry description of a molecule required for understanding its energy level pattern and the properties deriving from this pattern. For rigid non-linear molecules, this symmetry description is identical to that arising from the molecular point group at equilibrium, and this explains the many successful, traditional applications of point group symmetry, especially in chemical contexts. For a rigid linear molecule, the infinite-order point group obviously provides a much more detailed symmetry description than the finite MS group. Again, one can argue that the MS group provides the symmetry operations relevant for describing the ‘fully coupled’ rovibronic (rotation-vibration-electronic) wavefunctions of a molecule and that the additional point group symmetry is redundant and unnecessary. In practice, however, the point group symmetry gives rise to useful information, in particular for the electronic, vibrational, and rotational basis functions used to express the fully coupled wavefunctions, and so it is advantageous to also employ the point group symmetry. The particular problems associated with the symmetry description of linear molecules were described early on by Hougen [3], and by Bunker and Papoušek [4]. The latter authors introduced the so-called Extended Molecular Symmetry (EMS) Group which, for a centrosymmetric linear molecule, is isomorphic to the D h point group. We discuss the EMS group in more detail below.
The aim of the present work is to illustrate how D h symmetry can be implemented into general nuclear-motion programs. As an example we use TROVE [5,6], a numerical variational method to solve for the ro-vibrational (rotational-vibrational) spectra of (small to medium) general polyatomic molecules, which has been used to simulate the hot spectra of various polyatomic molecules [7,8,9,10,11,12,13,14,15,16,17,18,19] as part of the ExoMol project [20,21], a database of ab initio spectra of molecules of astrophysical importance. Since calculations of this type are based on electronic properties of the molecule, primarily the potential energy and dipole moment surfaces, obtained by solving the electronic Schrödinger equation, one can claim that such calculations are also rovibronic. However, as the calculations are directly concerned with the rotation and vibration, we characterize them as ro-vibrational. The use of molecular symmetry has applications in diverse fields, including molecular spectroscopy and the construction of molecular wavefunctions, ligand-field theory, material science, and electronic structure calculations [1,22,23,24,25]. While the use of a symmetry-adapted basis set has been shown to make calculations of ro-vibrational energies far more efficient by reducing the size of the Hamiltonian matrix blocks to be diagonalized [22], it is not strictly necessary. This is in contrast to intensity calculations, which would hardly be practicable without knowledge about the symmetry of the ro-vibrational states, mainly due to the selection rules imposed by the nuclear spin statistics associated with different irreducible representations [17,22]. A symmetry-free intensity calculation would involve the actual numerical computation of a colossal number of vanishing intensities for transitions that do not satisfy the symmetry selection rules. In this context, it is also important that some energy levels have zero weights and do not exist in nature. Without knowledge of how the eigenvectors transform under the symmetry operations, it is therefore impossible—or at least immensely complicated—to describe the molecular spectrum correctly. We show how a symmetrization procedure similar to that of TROVE can be extended to enable symmetry classification, in particular of vibrational basis functions, in the D h point group and thereby introduce the possibility of labelling these basis functions by the value of the vibrational angular momentum quantum number (see, for example, Ref. [26]). In practice, it turns out that the infinitely many elements in D h represent a problem in numerical calculations; we circumvent this by employing, instead of D h , one of its subgroups D n h with a finite value of n. We discuss below how to choose an adequate n-value.
In numerical calculations, the vibrational and rotational basis functions are initially symmetry classified before they are combined to form the ro-vibrational basis. For a linear molecule, only combinations with k = are physically meaningful, where k is the z-axis-projection of the rotational angular momentum quantum number and is the vibrational angular momentum quantum number (see, for example, Refs. [1,18,26,27,28]). With the extended symmetrization procedure of the present work, the vibrational basis functions can be labelled by their -values and it becomes straightforward to construct the meaningful combinations. In a given ro-vibrational calculation, the required extent of the rotational excitation is defined by the maximum value J max of the angular momentum quantum number J. The maximum values of | k | and | | , K max and L max , respectively, are then K max = L max = J max . However, in practise the numerical calculations are computationally limited by the total number of quanta representing vibrational bending modes, which controls the maximum value for | | , and thus | k | . We find that the group D n h suitable for symmetry classification in numerical calculations has an n-value determined by K max = L max .
The TROVE symmetrization approach makes use of a set of simplified, ‘reduced’ vibrational Hamiltonians, each one describing one vibrational mode of the molecule. The symmetrization is achieved by utilizing the fact that each of these Hamiltonian operators commutes with the operations in the symmetry group of the molecule in question [22], so that eigenfunctions of a reduced vibrational Hamiltonian generate irreducible representations of the symmetry group. Consequently, we obtain a symmetry-adapted ro-vibrational basis set numerically by solving the eigenvalue problems for the reduced Hamiltonians; the vibrational basis functions are products of the eigenfunctions thus obtained. In the course of the present work, we have implemented in TROVE a general subroutine to generate automatically all transformation matrices associated with the irreducible representations of a given symmetry group D n h ; this can be applied to general nuclear motion routines. The matrices are chosen to describe the transformation of vibrational basis functions that are eigenfunctions of the operator L ^ z (with eigenvalues ) representing the vibrational angular momentum.
To the best of our knowledge, no general transformation matrices for D h have been reported in the literature although the corresponding character tables have been published many times (see, for example [29]). Hegelund et al. [30] have described the transformation properties of the customary rigid-rotor/harmonic-oscillator basis functions (see, for example, Refs. [1,26,31]) for D n h point groups with arbitrary n ⩾ 3 (see also Section 12.4 of Ref. [1]). The basis functions span the irreducible representations of D n h and the coefficients obtained, defining the transformation properties, can straightforwardly be organized as transformation matrices. The present paper aims at providing the missing information for D h . As an illustration, we present how this symmetry information is implemented in TROVE as part of the automatic symmetry adaptation technique [22].
The paper is structured as follows. Section 2 gives an overview of the rotational and vibrational symmetry classifications and groups for a centrosymmetric linear molecule, and Section 3 presents the corresponding irreducible-representation transformation matrices and character tables. The symmetrization approach implemented in TROVE is outlined in Section 4, followed by some numerical examples in Section 5. Our conclusions are given in Section 6.

2. Rotational and Vibrational Symmetry

2.1. The Groups D h (M), D h (EM), and D h

Our general aim is to construct a symmetry adapted basis set for centrosymmetric linear molecules (such as, for example, CO 2 or C 2 H 2 ) to be used in variational solutions of the ro-vibrational Schrödinger equation [22]. We employ basis functions that are products of rotational and vibrational factor functions:
Φ J , k , v , l = Φ J , k rot Φ v , vib ,
where J is the rotational angular momentum quantum number, k is the projection of the angular momentum on the molecule-fixed z axis, v is a generic vibrational quantum number and is the vibrational angular momentum quantum number. We use the physically meaningful basis functions having k = [1,26]. A common choice for the rotational basis functions are the rigid-symmetric-rotor functions [1,26] Φ J , k rot = | J , k , m , from which we will omit the rotational quantum number m (the projection of the rotational angular momentum on the space-fixed z-axis) from Φ J , k rot since nothing depends on it in the situation of no external electric or magnetic fields, as considered here.
The complete internal wavefunction Φ int (see Chapter 8 of Ref. [1]) is subject to Fermi-Dirac and Bose-Einstein statistics [1]. In the present work, we neglect the dependence of the energy on the nuclear spin, and so we take Φ int = Φ elec Φ rv Φ ns where Φ elec is the electronic wavefunction (which, as mentioned above, we assume to describe a totally symmetric singlet electronic state), Φ rv is the rotation-vibration wavefunction represented in the variational calculation by a linear combination of the basis functions in Equation (2), and Φ ns is a nuclear-spin wavefunction. Nuclear spin statistics requires Φ int to change sign under the operation ( p ) in Equation (1) if ( p ) involves an odd number of odd permutations of fermions [1], and to be invariant under ( p ) in all other cases. If ( p ) Φ int = + Φ int ( Φ int ), Φ int has Σ g + or Σ u ( Σ u + or Σ g ) symmetry in the group D h (M) of Equation (1) (Table 2), depending on whether the parity p is + 1 ( 1 ). The parity is the character under E * : E * Φ int = p Φ int .
The nuclear-spin wavefunction Φ ns does not depend on the spatial coordinates of the nuclei and so it is invariant under the “geometrical” symmetry operations of the point group D h . It is also invariant under E * but it may have its sign changed by ( p ) . Thus, it can have Σ g + or Σ u + symmetry in D h (M) (Table 2).
We note that only the operation ( p ) D h (M) [Equation (1)] is relevant for the discussion of Fermi-Dirac and Bose-Einstein statistics in connection with the complete internal wavefunction Φ int . E * D h (M) is also a “true” symmetry operation, but the operations in the point group D h do not occur naturally in this context. However, as mentioned above it is advantageous also to make use of the D h symmetry, and for this purpose Bunker and Papoušek [4] defined the EMS group D h (EM) which is isomorphic to D h . The operations in D h (EM) can be written as (see also Section 17.4.2 of Ref. [1])
E ( = E 0 ) , E ε , ( p ) ε , E ε * , ( p ) ε *
where the angle ε satisfies 0 ε < 2 π and is chosen independently for the operations in Equation (3). A general element O ε of the EMS group is defined as follows: (i) The effect of O ε on the spatial coordinates of the nuclei and electrons in the molecule is the same as that of the element O of the MS group; (ii) The effect of O ε on the Euler angles [1] θ and ϕ is the same as the effect of O of the MS group; (iii) The effect of O ε on the Euler angle [1] χ is defined by Eqs. (17-101)–(17-104) of Ref. [1], so as to mimic a rotation by ε about the molecular axis.
The irreducible representations of D h and D h (EM) are listed in Table 1. Four of them are one-dimensional (1D): Σ g + , Σ g , Σ u + , and Σ u and an infinite number are two-dimensional (2D): Π g / u , Δ g / u , Φ g / u , Γ g / u , H g / u , I g / u , …. As indicated in the table, an equivalent notation is A 1 g , A 2 g , A 2 u , A 1 u for 1D irreps and E n g and E n u , where n = 1 , 2 , , for 2D irreps, see Table 1. The rotational basis functions Φ J , k rot = | J , k span the irreducible representations of D h (EM) as given in Table 3.
For a linear centrosymmetric molecule, both the rotational and vibrational basis functions can be classified according to the irreps of the infinite-order D h (EM). This group is defined such that the effect of the operations on the vibronic coordinates are identical to those of the point group D h . It follows from the discussion given above, however, that only the operations in the MS group D h (M) (corresponding to ε = 0 for the operations in D h (EM)) are relevant for determining the requirements of Fermi-Dirac and Bose-Einstein statistics. D h (EM) operations with ε > 0 are artificial (in the sense that the complete ro-vibrational Hamiltonian does not commute with them [1]) and therefore the basis function Φ J , k , v , l from Equation (2) must be invariant to them—we can view this as a “reality check” of Φ J , k , v , l , which turns out to be invariant to the artificial operations for k = . It is seen from Table 1 that consequently, Φ J , k , v , l can only span one of the four irreducible representations Σ g + , Σ g , Σ u + , and Σ u of the EMS group D h (EM). In the Introduction, we already gave examples of the weird and wonderful universe of spectroscopic notation. We now extend this universe by pointing out that according to the labelling scheme of Ref. [32], the four irreps Σ g + , Σ g , Σ u + , and Σ u are also denoted e ortho, e para, f ortho and f para. The correspondence depends on whether J is even or odd and is given in Table 4.
The rotational and vibrational factor wavefunctions Φ J , k rot and Φ v , vib , respectively, in Equation (2) are symmetry classified in D h (EM) and there are no restrictions as to their possible symmetries. However, the fact that the product function Φ J , k , v , l must transform according to a 1D irrep introduces restrictions on the possible combinations of Φ J , k rot and Φ v , vib ; these restrictions limit the physically useful combinations to those with k = . For example, the vibrational state ν 5 (with vibrational basis functions Φ v 5 = 1 , = ± 1 vib of Π u symmetry in D h (EM)) of acetylene C 2 H 2 can be combined with the Φ J , k rot rotational wavefunctions having ( J , k ) = ( 1 , ± 1 ) (and Π g symmetry) to produce three ro-vibrational combinations with symmetries Σ u + , Σ u and Π u in D h (EM). However, only the Σ u + and Σ u states can be used in practice and the Π u state must be discarded.

2.2. The Point Groups D n h and Their Correlation with D h

In numerical, variational ro-vibrational calculations, the symmetrization and symmetry classification of rotational and vibrational basis functions facilitate the actual calculations, since the matrix representation of the ro-vibrational Hamiltonian, which is diagonalized numerically in a variational calculation, becomes block diagonal according to the symmetries of the basis functions [1]. In addition, the resulting eigenfunctions are automatically symmetrized and can be labelled by the irrep that they generate. Without this, the calculations would produce redundant energies, there would be no way to determine the appropriate nuclear-spin statistics to be applied to a given state, and it would be impossible to identify the rotation-vibration transitions allowed by symmetry selection rules [1]. In particular, one could not determine the nuclear spin-statistical weight factors g ns entering into intensity calculations (for 12 C 2 H 2 the spin-statistical weight factors are 1 for Σ g + and Σ u (para) ro-vibrational states and 3 for Σ g and Σ u + (ortho) states). There are no allowed electric dipole transitions between ortho and para states [17,32].
D h is the geometrical symmetry group of a (horizontal, say) circular disc whose upper and lower surfaces are equivalent so that one can turn the disc upside-down without any observable change resulting from this. Similarly, D n h is the geometrical symmetry group of a (horizontal, say) regular polygon with n vertices (i.e., a regular n-gon) whose upper and lower surfaces are equivalent. That is, we can think of D h as the limiting case of a progression of D n h groups: D h = lim n D n h . As mentioned above, we aim at implementing D h symmetry for the ro-vibrational basis functions employed in variational calculations. However, owing to the infinitely many operations of D h and the corresponding infinitely many irreps, this is impracticable. Consequently, we resort to the strategy often used in numerical calculations and approximate ∞ by a large, finite number n or, in other words, we approximate D h by D n h with a suitably large n-value. In order to do this, we must discuss the correlation between D n h and D h .
In Table 5 we list the symmetry operations in D n h . It is seen that we must distinguish between n even and odd. The difference in group structure—and an accompanying difference in the labelling of the irreps—are caused by the fact that for n even, the point group inversion i (as explained in connection with Equations (4)–(7) of Ref. [1], the point group inversion operation i is different from the spatial inversion operation E * and should be careful to distinguish between the two) is present in D n h , whereas for n odd it is not. Since iD h , in some sense an even-n D n h is more similar to D h than an odd-n D n h . It could be argued that only even-n D n h groups should be considered in the limit of n ; this is the approach we have taken in the example calculations of Section 4. However, for completeness we also discuss odd-n D n h groups here. The corresponding information is potentially useful for the treatment of polyatomic molecules with an MS group isomorphic to an odd-n D n h group.
It will be shown (see Section 4 and Section 5) that the optimum value for n for the D n h group used in a TROVE calculation to approximate D h depends on the maximum value L max of the vibrational angular momentum number required for a given calculation. We have L max = K max , the maximum value on the z-axis projection of the rotational angular momentum. In practical calculations we are usually limited by L max , as determined by the maximum total value of vibrational bending quanta, rather than by K max , as determined by the maximum quanta of rotational excitation.
The general formulation of the irreducible representations of D n h for arbitrary n is outlined in Section 3 below.

3. General Formulation of the Character Tables and the Irreducible Representation Transformation Matrices of the D n h Groups

3.1. General Structure

Let us consider for a moment the point group C 3 v . It contains the six operations
C 3 v = E , C 3 , C 3 2 , σ ( x z ) , σ ( 2 ) , σ ( 3 ) .
We have chosen a right-handed axis system such that C 3 and C 3 2 are rotations of 2 π 3 and 4 π 3 , respectively, around the z axis. The positive direction of the rotations is defined as the direction in which a right-handed screw will rotate when it advances in the positive direction of the z axis. The x and y axes are perpendicular to the z axis and chosen such that the group operation σ ( x z ) is a reflection in the x z plane. The operations σ ( 2 ) and σ ( 3 ) are then reflections in planes that contain the z axis and form angles of 2 π 3 and 4 π 3 , respectively, with the x z plane.
It is straightforward to verify the following relations
E = C 3 C 3 C 3 = σ ( x z ) 2
C 3 2 = C 3 C 3
σ ( 2 ) = C 3 C 3 σ ( x z )
σ ( 3 ) = C 3 σ ( x z ) .
In Equations (5)–(8) all operations in the group C 3 v have been expressed as products of the two operations C 3 and σ ( x z ) . These operations are called the generating operations for C 3 v . It is clear that in order to symmetry classify an operator (or a function) in C 3 v , it is sufficient to know how the operator (or function) transforms under the generating operations C 3 and σ ( x z ) . With this knowledge, Equations (5)–(8) can be used to construct the transformation properties under all other operations. All point groups can be defined in terms of generating operations. Hegelund et al. [30] have shown that for a general group C n v the generating operations can be chosen as C n and σ ( x z ) by analogy with the choice for C 3 v .
Two simple isomorphic cyclic groups, C s and C i , can now be introduced:
C s = E , σ h
C i = E , i
where σ h is a reflection in a horizontal plane (perpendicular to the n-fold axis) and i is the point group inversion. The irreps of these groups are given in Tables A-2 and A-3 of Ref. [1]. It can be shown [35] that the D n h groups can be written as direct products of these simple groups:
D n h = C n v C s ( n odd )
D n h = C n v C i ( n even ) .
That is, an odd-n D n h contains all elements RC n v together with all elements that can be written as R σ h , and an even-n D n h contains all elements RC n v together with all elements that can be written as R i .
As explained in Section 12.4 of Ref. [1] all operations in a D n h group can be obtained as products involving three generating operations which are typically denoted by R + , R + , and R . The generating operations for the D n h groups are summarized in Table 6. R + is chosen as C n for all n, but for n odd, ( R + , R ) = ( σ h , C 2 ( x ) ) , whereas for n even, ( R + , R ) = ( i , C 2 ( x ) ) , where C 2 ( x ) is a rotation by π about the molecule-fixed x axis.
Owing to the direct product structure of the D n h groups (Equations (11) and (12)) it would in fact have been more logical to choose σ ( x z ) as a generating operation for D n h instead of C 2 ( x ) . However, this does not seem to be the customary choice (see, for example, Hegelund et al. [30]) and we attempt here to following accepted practice as much as possible. With the relations
σ ( x z ) = C 2 ( x ) σ h ( n odd ) , and
σ ( x z ) = C 2 ( x ) C n n / 2 i ( n even ) .
It is straightforward to express the elements of D n h in terms of the chosen generating operations. Here, C n n / 2 is a rotation by π about the z axis.
When the transformation properties of an object under R + , R + , and R are known, the transformation properties under all other operations in a D n h point group can be unambiguously constructed.

3.2. Irreducible Representations

As described above, the structure of the D n h point groups alternates for even and odd n-values. Consequently, so do the transformation matrices generated by the rotation-vibration basis functions (see Refs. [29,30] and Section 5.1.2 of Jensen and Hegelund [35]). The irreducible representations of D n h point groups are easily constructed for arbitrary n as described in Section 5.8.2 of Ref. [1], as listed in Table 7. The irreps are expressed in terms of the characters under the generating operations R + , R + , and R which are also given in Table 7 (The labelling of the irreducible representations by A, B, and E is recommended by the International Union of Pure and Applied Chemistry; see Section 2.3.3 of Ref. [36]. Why these three labels were chosen remains one of the mysteries in the weird and wonderful world of spectroscopic notation).
Comparison of Table 1 and Table 7 shows that an even-n D n h group has four 1D irreps called A 1 g , A 2 g , A 1 u , and A 2 u and ( n 2 ) 2D irreps, of which half are called E r g and the other half E r u (r = 1,2, …, n / 2 1 ). All of these irreps correlate with irreps of D h denoted by the same names in Table 1. In addition, the even-n D n h group has another four 1D irreps called B 1 g , B 1 u , B 2 g , B 2 u associated with a sign change of the generating function under the C n rotation (Table 7). These B-type irreps have no counterparts in D h and so basis functions of these symmetries are useless, if not nonphysical, in the context of approximating D h by D n h . We noted above that the point group inversion operation i is contained in D h and in even-n D n h , but not in odd-n D n h . Therefore the labelling of the irreps of odd-n D n h differs from that used for D h and in even-n D n h . However, Table 8 gives the correspondence between the irreps of odd-n D n h and those of even-n D n h and D h , and so we have established the correlation between the D n h and the D h irreps for all n-values.

3.3. Transformation Matrices

In practical applications of representation theory, such as the symmetry adaptation and description of basis functions that are the subject of the present work, it is not sufficient to have the irreducible-representation characters of Table 7 only. We also need groups of matrices that constitute irreducible representations of the D n h group with an arbitrary finite n-value. For the 1D irreps (of type A and B, in the notation of Table 1 and Table 7) the 1 × 1 transformation matrix is simply equal to the character in Table 7. For the 2D irreps (of type E, in the notation of Table 1) we require 2 × 2 matrices whose traces are the characters in Table 7. Once a set of irreducible-representation matrices are known, symmetrized basis functions (with transformation properties defined by the irreducible-representation matrices) can in principle be determined by the projection-operator technique described in Section 6.3 of Ref. [1].
Representation matrices are not uniquely determined. Having determined one set of, say, 2 × 2 matrices M R that constitute an E-type irreducible representation of a D n h group, for any 2 × 2 matrix V with a non-vanishing determinant we can construct an equivalent representation consisting of the matrices V M R V 1 as explained in Section 5.4.1 of Ref. [1]. We normally consider representation matrices effecting the transformation under the group operations of particular wavefunctions, coordinates or operators.
We consider here the transformation/representation matrices generated by the rotational basis functions | J , k . The relative phases of these functions are chosen in the customary manner as given in Section 11.2.3 of Ref. [1] so that the matrix elements of the “molecule-fixed” angular momentum ladder operators are real and positive. To determine all transformation matrices, it is sufficient initially to know the transformation properties of these functions under the generating operations R + , R + , and R (see Table 7). When these are known, the transformation matrix for any group operation R is uniquely determined; one determines the product involving R + , R + , and R that equals R and the desired transformation matrix for R is the analogous matrix product of the representation matrices of R + , R + , and R , respectively.
The transformation properties of the | J , k functions under the generating operations R + , R + , and R are straightforwardly determined from the results of Hegelund et al. [30] which are reproduced in Section 12.4 of Ref. [1]. Table 9 gives the 1 × 1 matrices generated by | J , 0 and the 2 × 2 matrices M R generated by ( | J , K , | J , K ) under R + , R + , and R for K> 0. It is advantageous also to generate the matrices M R = V M R V 1 generated by the so-called Wang functions | J , 0 , + = | J , 0 for K = 0 and
| J , K , + | J , K , = V | J , K | J , K
For K> 0, where
V = 1 2 1 2 i 2 i 2 .
The function | J , 0 , + obviously generates the same 1 × 1 matrices as | J , 0 . These, along with the 2 × 2 matrices M R and M R are included in Table 9.
We are now in a situation to generate the 2 × 2 transformation matrices for 2D irreps of the D n h group, spanned by ( | J , K , + , | J , K , ) for K > 0. Towards this end, we use the fact that any element RD n h can be expressed as a product of the generating operations R + , R + , and R , and that the transformation matrix M R generated by R can be expressed as the analogous matrix product of the representation matrices M R + M R + , and M R in Table 9. In forming the matrix products, one can make use of the fact that all M R = V M R V 1 , where the M R matrices are generated by ( | J , K , | J , K ) (Table 9) for K > 0 and the matrix V is defined in Equation (15). For example, all D n h groups contain the operations C n r , where r = 1, 2, …, n 1 . The operation C n r = R + r thus generates the transformation matrix
M C n r = M C n r = V M C n r V 1 = V e + i r K ε 0 0 e i r K ε V 1 = cos ( r K ε ) sin ( r K ε ) sin ( r K ε ) cos ( r K ε ) ,
With C n r = R + r C n v , r = 1, 2, 3, …, n 1 . In general, C n v further contains n reflections in planes that contain the C n axis, customarily chosen as the z axis of the molecule-fixed axis system. As discussed for C 3 v in connection with Equations (5)–(8), we can start with one such reflection, σ ( x z ) say, and then obtain the other n 1 reflections as C n r σ ( x z ) , r = 1, 2, 3, …, n 1 . However, σ ( x z ) is not chosen as a generating operation for D n h (see Table 7), but we can use Equations (13) and (14) to express the n reflections as
σ ( r ) = C n r C 2 ( x ) σ h ( n odd ) , and
σ ( r ) = C n r C 2 ( x ) C n n / 2 i ( n even )
where r = 0, 1, 2, 3,…, n 1 . Again, C n n / 2 is a rotation by π about the z axis. For n odd, all n reflections are of the same type and the reflection planes all contain one vertex of the regular n-gon whose geometrical symmetry we consider. For n even, we obtain two different reflection types: n / 2 reflections obtained for even r = 0, 2, 4, n 2 , and n / 2 reflections obtained for odd r = 1, 3, 5, n 1 . The even-r type reflection are of the σ v type, with the reflection plane containing two vertices of the regular n-gon, while the reflection planes of the odd-r type, said to be of σ d type, bisect the angle between neighbouring pairs of σ v reflection planes and contain no vertices of the n-gon (the use of the label subscripts v (for vertical) and d (for diagonal) is recommended by the International Union of Pure and Applied Chemistry; see Section 2.1.2 of Ref. [36]. Why v and d were chosen vanishes in the historical fog occupying a large part of the world of spectroscopic notation). We have now constructed all elements of C n v , and we straightforwardly augment this group by 2 n elements; R i for n even and R σ h for n odd, with RC n v (Equations (11) and (12)). The n 1 operations C n r i ( C n r σ h ), r = 1, 2, 3,…, n 1 are improper rotations S n ( r ) for n even(odd). We see from Equations (18)–(19) that the remaining operations of type σ ( r ) i (n even) and σ ( r ) σ h (n odd) can be written as
C 2 ( r ) = σ ( r ) σ h = C n r C 2 ( x ) ( n odd ) , and
C 2 ( r ) = σ ( r ) i = C n r C 2 ( x ) C n n / 2 ( n even ) .
These operations are rotations by π about axes perpendicular to the C n axis which are contained in the plane of the regular n-gon. For n odd, each of these C 2 axes passes through one vertex of the n-gon and all of the C 2 rotations are equivalent. For n even, there are two types of C 2 rotations, depending on whether the σ ( r ) operation in Equation (21) is of type σ v or σ d . If it is σ v then the rotation by π is of type C 2 and the corresponding rotation axis passes through two vertices of the n-gon. If, on the other hand, it is of type σ d then the rotation axis of the corresponding C 2 rotation is contained in a σ d reflection plane and bisects the angle between two neighbouring C 2 axes.
We have now explained how for an arbitrary n-value, each operation in D n h can be expressed as a product of the generating operations given in Table 7. To generate a corresponding set of representation/transformation matrices, we must derive the analogous matrix products of the representation matrices in Table 9. The resulting representation matrices are given in Table 10 for n even and in Table 11 for n odd.
It is seen in Table 12 that the rotational basis functions | J , 0 and ( | J , K , + , | J , K , ) generate g-type (even) symmetries of D n h only for n even. This is because, as explained in Section 4.5 of Ref. [1], the point group inversion i and its MS-group counterpart O ^ i do not change the Euler angles, i.e., the rotational coordinates. Consequently, the rotational functions are invariant to these operations. We have added to Table 10 transformation matrices also for u-type (odd) irreps. These matrices can be thought of as generated by functions ( | J , K , ± | v u = 1 , where | v u = 1 is the vibrational wavefunction for the fundamental level of a (probably hypothetical) vibrational mode ν u of A 1 u symmetry.
If one uses Table 10 and Table 11 to determine a set of representation/transformation matrices for a given 2D irrep, it is important to realize that the matrices given are generated by the ( | J , K , + , | J , K , ) rotational basis functions or, in the case of E κ u for n even, by the ro-vibrational functions | J , K , ± | v u = 1〉 defined above. One must choose J and K values so that they are commensurable with the irrep considered. The most important restriction here is that for n odd, E κ ( E κ ) symmetry requires even(odd) K (see Table 12). For a given K value, one can make the always-physical choice of J = K.

4. Symmetrization Using the TROVE Approach

In this section we use the variational nuclear-motion program TROVE as an example of a practical implementation of the symmetrization of ro-vibrational basis states using D n h representations. TROVE uses a general numerical symmetrization approach to build a symmetry-adapted ro-vibrational basis set [22]. The procedure will be outlined here and extended to include classification based on the vibrational angular momentum quantum number, , as necessary for dealing with linear molecules of D h point group symmetry, using 12 C 2 H 2 as an example. This classification is general and can be implemented into other similarly constructed variational routines.
The use of a symmetry-adapted basis set can considerably reduce the size of the Hamiltonian matrix to be diagonalized. This is due to the useful property that the matrix elements between basis functions of different symmetry are zero by definition:
Ψ μ J , Γ s , α | H r v | Ψ μ J , Γ t , α = H μ , μ δ s , t δ α , α ,
where Γ s and Γ t give the irreducible representations (irreps) of D n h that the basis functions, Ψ μ J , Γ s , α and Ψ μ J , Γ t , α , transform according to, and α and α represent their degenerate components (if present). The block diagonal structure of a Hamiltonian matrix in the D n h irreducible representation is given in Figure 1; the symmetry blocks of non-vanishing matrix elements can be diagonalized separately.
TROVE utilises the concept of a sum-of-product basis set, where the primitive basis functions are
Φ k , v , l J ( θ , ϕ , χ , ξ 1 , ξ 2 ξ N ) = | J , k , m ϕ v 1 ( ξ 1 ) ϕ v 2 ( ξ 2 ) ϕ v N ( ξ N ) ,
With 1D vibrational basis functions ϕ v i ( ξ i ) (where ξ i is a generalised vibrational coordinate) and rigid-rotor (spherical harmonics) rotational basis functions | J , k , m . The 1D vibrational basis functions are either obtained by solving the corresponding reduced 1D Schrödinger equations or are taken as the harmonic or Morse oscillators.

Symmetrization of the Basis Set for 12 C 2 H 2 Using the ( 3 N 5 ) Coordinate TROVE Implementation

As an illustration of the practical application of the finite D n h group being used in place of D h , we show an example of the construction of the vibrational basis set in case of the linear molecule 12 C 2 H 2 . We use the recent implementation of the ( 3 N 5 ) coordinates approach in TROVE (see Ref. [18]) and select a set of seven vibrational coordinates used for 12 C 2 H 2 : Δ R , Δ r 1 , Δ r 2 , Δ x 1 , Δ y 1 , Δ x 2 , Δ y 2 , as illustrated in Figure 2. The transformation matrices defining their symmetry properties are listed in Table 13 (with even n used in this example). These relate to the symmetry operations of Table 5, and the general irrep transformation matrices for D n h of even n given in Table 10.
For the stretching primitive basis functions, ϕ v 1 ( ξ 1 ) , ϕ v 2 ( ξ 2 ) and ϕ v 3 ( ξ 3 ) , we use eigenfunctions of the corresponding 1D reduced stretching Hamiltonian operators H ^ i ( 1 D ) obtained by freezing all other degrees of freedom at their equilibrium values in the J = 0 Hamiltonian. We use the Numerov-Cooley approach [5,37,38] to solve these eigenvalue problems. For the bending basis functions ϕ v 4 ( ξ 4 ) , , ϕ v 7 ( ξ 7 ) , 1D harmonic oscillators are used. Here ξ = { Δ R lin , Δ r 1 lin , Δ r 2 lin , Δ x 1 , Δ y 1 , Δ x 2 , Δ y 2 } , where R lin , Δ r 1 lin , and Δ r 2 lin are linearized (rectilinear) versions of Δ R , Δ r 1 , Δ r 2 , respectively, as described in Ref. [18].
According to the TROVE symmetrization technique [22], the symmetry adapted vibrational basis functions are formed from linear combinations of the products of the 1D vibrational basis functions ϕ v i ( ξ i ) as follows, with 12 C 2 H 2 used as an example. For 12 C 2 H 2 , the vibrational part of the basis set of Equation (23) is divided into three sub-sets:
ϕ v 1 ( 1 D ) ( ξ 1 ) = ϕ v 1 ( ξ 1 ) ,
ϕ v 2 v 3 ( 2 D ) ( ξ 2 , ξ 3 ) = ϕ v 2 ( ξ 2 ) ϕ v 3 ( ξ 3 ) ,
ϕ v 4 v 5 v 6 v 7 ( 4 D ) ( ξ 4 , ξ 5 , ξ 6 , ξ 7 ) = ϕ v 4 ( ξ 4 ) ϕ v 5 ( ξ 5 ) ϕ v 6 ( ξ 6 ) ϕ v 7 ( ξ 7 ) .
These sub-sets are then used as basis sets for the corresponding reduced Hamiltonian operators: stretching H ^ ( 1 D ) and H ^ ( 2 D ) , and bending H ^ ( 4 D ) . The reduced Hamiltonians H ^ ( N D ) ( N = 1 , 2 , 4 ) are constructed by averaging the total vibrational Hamiltonian operator H ^ ( J = 0 ) over the other ground vibrational basis functions. For example, H ^ ( 4 D ) is given by:
H ^ ( 4 D ) = 0 | 0 , 0 | H ^ ( J = 0 ) | 0 , 0 | 0 ,
where | 0 = ϕ 0 ( 1 D ) ( ξ 1 ) and | 0 , 0 = ϕ 0 , 0 ( 2 D ) ( ξ 2 , ξ 3 ) .
According to the idea of the so-called complete set of commuting operators (CSCO) [22] which the TROVE symmetrization approach is based on, the eigenfunctions of the reduced operator H ^ ( N D ) must transform according to one of the irreps of the symmetry group of the system, since H ^ ( N D ) commutes with the symmetry operators of this group. Thus the symmetrization of the basis set is generated automatically by solving the appropriate eigenvalue problem, provided that the corresponding irreps have been determined. To this end, the symmetry operators of the appropriate group can be applied to the eigenfunctions and their transformation properties on a set of sampled geometries (usually 40–60 in the case of TROVE) analysed. Some states of the same energy (either with accidental or actual degeneracy) may appear as random mixtures of each other, and have to be processed simultaneously and even further reduced to irreps, if necessary (see Section 5 for an example).
Applying this procedure to stretching functions gives rise to A-type symmetries: e.g., for D n h (even n), the eigenfunctions of H ^ ( 1 D ) span the A 1 g irrep, while the eigenfunctions of H ^ ( 2 D ) span the A 1 g and A 2 u irreps.
The 4D bending basis set, based on the 1D harmonic oscillators of Equation (26), has the disadvantage of being extremely degenerate: combinations of ϕ v 4 v 5 v 6 v 7 ( 4 D ) give rise to large clusters of the same energies. According to the TROVE symmetrization approach these combinations must be processed together, which makes this process extremely slow. In order to facilitate this step we first transform the 4D bending sets (Equation (26)) to become eigenfunctions of the vibrational angular momentum operator,
L ^ z = λ , λ ξ λ ζ λ , λ z p ^ λ
where p λ is a vibrational momentum operator and ζ λ , λ z are Coriolis coefficients [39] as described in Ref. [18].
TROVE is equipped to compute matrix elements of quadratic forms, therefore we use L ^ z 2 instead of L ^ z . Using the ϕ v 4 v 5 v 6 v 7 ( 4 D ) basis functions we find eigenfunctions of L ^ z 2 by diagonalizing the matrix formed by combinations of the 4D bending basis set of Equation (26):
ϕ v 4 v 5 v 6 v 7 ( 4 D ) | L ^ z 2 | ϕ v 4 v 5 v 6 v 7 ( 4 D ) .
The eigenfunctions of L ^ z 2 are consequently characterized by their vibrational angular momentum quantum number = | | = 2 and can thus be divided into independent sub-sets with different symmetry properties: the L = 0 sub-set must be a mixture of A-type functions, while the L > 0 sub-sets consist of the E L -type irreps ( E L g and E L u ). These mixtures are then further reduced to irreps using the TROVE symmetrization scheme outlined above, in which the reduced 4D-eigenvalue problem, using the eigenfunctions of L ^ z 2 as the basis set, is solved for a 4D isotropic harmonic oscillator Hamiltonian:
H ^ 4 D = 1 2 p ^ 4 2 + p ^ 5 2 + p ^ 6 2 + p ^ 7 2 + 1 2 λ ξ 4 2 + ξ 5 2 + ξ 6 2 + ξ 7 2 ,
where λ is a related to the harmonic vibrational wavenumber and p ^ i are the vibrational momenta, conjugate to ξ i . Thus we obtain eigenfunctions which can be divided into sub-sets of the same energies and values of . These sub-sets must transform independently, thereby significantly decreasing the time spent on the symmetry sampling step by breaking the symmetry space into small sets and making numerical calculations more computationally viable. Although the L ^ z 2 -diagonalisation step is not strictly necessary for the general symmetrization procedure that follows it, this increase in efficiency is a big advantage.
As mentioned above, in addition to the -quantum number being advantageous in building the vibrational basis sets, it is also required for coupling the basis set functions according to the linear molecule angular momentum rule k = (see, for example, Refs. [18,27,28]). The maximum value for L max = K max is specified as an input into the TROVE numerical routine.
As a result of applying the procedure described above, a symmetry-adapted vibrational basis set Φ v , L Γ vib , α is generated. Here Γ vib is the irrep of the basis function according to D n h , and α indicates a degenerate component in the case of 2D irreps.
The symmetry-adapted rotational basis set in TROVE is represented by:
| J , K , τ Γ rot = i σ 2 | J , K + ( 1 ) J + K + τ | J , K ,
where K = 0 is a special case, given by:
| J , 0 , τ Γ rot = | J , 0 .
Here | J , k is a rigid rotor wavefunction, τ ( = 0 , 1 ) is a parameter used to define the parity of a state, where σ = ( K mod 3 ) for τ = 1 and σ = 0 for τ = 0 (see [22,40,41]). The irreps Γ s of these functions are listed in Table 14, where τ defines their degenerate component. The symmetry properties of | J , K , τ Γ rot can be derived from those of | J , k using the method described in Section 3.3.
The symmetrized rotational and vibrational basis functions are then combined to form a full ro-vibrational symmetry-adapted basis set:
Ψ v , K J , Γ s = α , τ T α , τ ( Γ vib , Γ rot ) Γ s Φ v , K Γ vib , α | J , K , τ Γ rot ,
where T α , τ ( Γ vib , Γ rot ) Γ s are symmetrization coefficients with α indicating a degenerate component in the case of 2D irreps, Γ s is a 1D irreps in D h (see Section 3) and the K = L condition for linear molecules in the ( 3 N 5 ) -approach [18] was applied. Note that the symmetrized basis functions use K and L instead of k and in Equation (23).

5. Numerical Example

Some test calculations were carried out using TROVE [5] for 12 C 2 H 2 using a small basis set. These calculations utilise the symmetrization procedure of Section 4.

5.1. Symmetrization

Here we give an example of building a symmetry adapted basis set for the 4D bending function of Equation (26) using the TROVE symmetrization approach. In this example, the size of the primitive basis sets was controlled by the polyad number
P = 2 v 1 + v 2 + v 3 + v 4 + v 5 + v 6 + v 7 P max ,
With P max = 2. Here, the quantum numbers v k for k = 1 7 correspond to the vibrational primitive functions ϕ v k ( ξ k ) .
Using the 4D reduced Hamiltonian in Equation (27) with this small basis set we obtain the following contracted eigenfunctions (with only the first seven given here):
Ψ 1 L = 0 = 0.9793 | 0000 0.0095 ( | 2000 + | 0200 + | 0020 + | 0002 ) + 0.1425 ( | 1010 + | 0101 ) , Ψ 2 L = 1 = 1 2 ( | 1000 | 0010 ) , Ψ 3 L = 1 = 1 2 ( | 0100 | 0010 ) , Ψ 4 L = 1 = 1 2 ( | 1000 + | 0100 ) , Ψ 5 L = 1 = 1 2 ( | 0100 + | 0010 ) , Ψ 6 L = 2 = 0.3505 ( | 2000 | 0200 + | 0020 | 0002 ) + 0.4957 ( | 1010 | 0101 ) + + 0.0651 ( | 1100 | 1001 | 0110 + | 0011 ) , Ψ 7 L = 2 = 0.0460 ( | 2000 | 0200 + | 0020 | 0002 ) 0.0651 ( | 1010 | 0101 ) + + 0.4957 ( | 1100 | 1001 | 0110 + | 0011 ) ,
One can see that after this step some of the eigenfunctions ( Ψ 1 , Ψ 2 , Ψ 3 , Ψ 4 and Ψ 5 ) are already in the form of an irreducible representation, while Ψ 6 and Ψ 7 need to be further reduced.
In order to define the L-values, the matrix elements of L ^ z are computed as in Equations (28) and (29). In this example, the sets with degenerate eigenvalues and identical L values are: { Ψ 1 L = 0 } ( E ˜ 1 = 0 cm 1 ), { Ψ 2 L = 1 , Ψ 3 L = 1 } ( E ˜ 2 = 636.11 cm 1 ), { Ψ 4 L = 1 , Ψ 5 L = 1 } ( E ˜ 3 = 763.12 cm 1 ), { Ψ 6 L = 2 , Ψ 7 L = 2 } ( E ˜ 4 = 1215.84 cm 1 ). The pair of eigenfunctions Ψ 2 L = 1 and Ψ 3 L = 1 , for example, both have L = 1 and are also degenerate ( ± 10 12 ) with the eigenvalue 636.11 cm 1 . All degenerate states are combined into the same set and are assumed to share the symmetry transformation properties, now with the additional constraint that those states in the same set must also possess the same value of L. For our example, this gives the following symmetries and L-values: { Ψ 1 L = 0 } A 1 g , { Ψ 2 L = 1 , Ψ 3 L = 1 } E 1 g , { Ψ 4 L = 1 , Ψ 5 L = 1 } E 1 u , { Ψ 6 L = 2 , Ψ 7 L = 2 } E 2 g .
The irreducible form of the wavefunctions Ψ 6 , Ψ 7 is now given by:
Ψ 6 L = 2 = 2 4 ( | 2000 | 0200 + | 0020 | 0002 ) 1 2 ( | 1010 | 0101 ) ,
Ψ 7 L = 2 = 1 2 ( | 1100 | 1001 | 0110 + | 0011 ) .

5.2. Even vs. Odd D n Symmetries

For the example calculations using even vs. odd D n h that are outlined below the primitive and contracted basis sets were controlled by the polyad number as given by Equation (34), with P max = 8 for the primitive basis set and reduced to 6 after contraction (see Refs. [5,18] for more details).
In place of using the infinite group D h , we use a finite group D n h , with a value of n large enough to cover all required excitations of the vibrational angular momentum L = | | up to up L max and of the rotational quantum number K up to K max (with the constraint L max = K max ) such that n = 2 L max + 1 or n = 2 L max + 2 (depending on whether n is odd or even, respectively). For example, in order to be able to cover the rotational excitation up to K = 10 ( E 10 g and E 10 u ), it is necessary to use at least the D 21 h symmetry.
Even though odd and even values of n lead to different symmetry operations (see Table 5), both lead to the same resulting eigenvalues energies in the TROVE calculations, with example energies and assignments given in Table 15, on the condition that n 2 L max + 1 (odd n) or n 2 L max + 2 (even n). For a maximum value of the z-projection of the vibrational angular momentum, L max = K max = 4, different values of n were used for D n h in the symmetrization approach described in Section 4.
If a lower value than n = 2 L max + 1 (for odd n) or n = 2 L max + 2 (for even n) is used, then the symmetrization procedure will lead to the wrong classification of states, resulting in, for example, the wrong nuclear statistics in intensity calculations. It should be noted that for practical numerical calculations we are limited by the maximum number of vibrational bending quanta which can be included in calculations, giving a limit on L max . We therefore refer to this as the deciding factor in what n for D n h to use. However, it is also dependent on K max , the maximum value required for the z-projection of rotational angular momentum quantum number J, which would ideally be limited by J max (working under the assumption that K = L for the ( 3 N 5 ) -approach to dealing with linear molecules; see [18] and, for example, [27,28]).

6. Conclusions

We have presented an outline of the method used to treat linear molecules of the D h point group using finite D n h symmetry (with arbitrary user-defined n) in numerical calculations, and given the implementation in nuclear motion routine TROVE as an example application. We have shown how a symmetrization scheme such as that used in TROVE can be extended by including the vibrational angular momentum operator L ^ z into the set of commuting operators, allowing the classification of basis sets based on vibrational angular momentum quantum number, L. Character tables and irreducible representation transformation matrices for D n h of general integer odd or even n have been presented, along with some numerical examples for 12 C 2 H 2 .
The work on 12 C 2 H 2 presented in Ref. [18] utilises this symmetrization procedure for linear molecules, as does work extending the room temperature line list of Ref. [18] to higher temperatures, which is currently in progress.

Author Contributions

Katy L. Chubb, Sergei N. Yurchenko and Per Jensen all contributed to the writing of the manuscript and the associated theoretical work. Katy L. Chubb and Sergei N. Yurchenko ran the example numerical calculations.

Acknowledgments

This work was supported by the STFC projects No. ST/M001334/1 and ST/J002925. We acknowledge the use of the UCL Legion High Performance Computing Facility (Legion@UCL), and associated support services, in the completion of this work. We also acknowledge use of the DiRAC@Darwin HPC cluster, which is the UK HPC facility for particle physics, astrophysics and cosmology and is supported by STFC and BIS.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bunker, P.R.; Jensen, P. Molecular Symmetry and Spectroscopy, 2nd ed.; NRC Research Press: Ottawa, ON, Canada, 1998. [Google Scholar]
  2. Longuet-Higgins, H. The symmetry groups of non-rigid molecules. Mol. Phys. 1963, 6, 445–460. [Google Scholar] [CrossRef]
  3. Hougen, J.T. Classification of Rotational Energy Levels. II. J. Chem. Phys. 1963, 39, 358–365. [Google Scholar] [CrossRef]
  4. Bunker, P.; Papoušek, D. The symmetry groups of linear molecules. J. Mol. Spectrosc. 1969, 32, 419–429. [Google Scholar] [CrossRef]
  5. Yurchenko, S.N.; Thiel, W.; Jensen, P. Theoretical ROVibrational Energies (TROVE): A robust numerical approach to the calculation of rovibrational energies for polyatomic molecules. J. Mol. Spectrosc. 2007, 245, 126–140. [Google Scholar] [CrossRef]
  6. Yachmenev, A.; Yurchenko, S.N. Automatic differentiation method for numerical construction of the rotational-vibrational Hamiltonian as a power series in the curvilinear internal coordinates using the Eckart frame. J. Chem. Phys. 2015, 143, 014105. [Google Scholar] [CrossRef] [PubMed]
  7. Yurchenko, S.N.; Tennyson, J. ExoMol line lists IV: The rotation-vibration spectrum of methane up to 1500 K. Mon. Not. R. Astron. Soc. 2014, 440, 1649–1661. [Google Scholar] [CrossRef]
  8. Yurchenko, S.N.; Barber, R.J.; Yachmenev, A.; Thiel, W.; Jensen, P.; Tennyson, J. A variationally computed T = 300 K line list for NH3. J. Phys. Chem. A 2009, 113, 11845–11855. [Google Scholar] [CrossRef] [PubMed]
  9. Sousa-Silva, C.; Hesketh, N.; Yurchenko, S.N.; Hill, C.; Tennyson, J. High Temperature partition functions and thermodynamic data for ammonia and phosphine. J. Quant. Spectrosc. Radiat. Transf. 2014, 142, 66–74. [Google Scholar] [CrossRef]
  10. Sousa-Silva, C.; Al-Refaie, A.F.; Tennyson, J.; Yurchenko, S.N. ExoMol line lists-VII. The rotation-vibration spectrum of phosphine up to 1500 K. Mon. Not. R. Astron. Soc. 2015, 446, 2337–2347. [Google Scholar] [CrossRef]
  11. Underwood, D.S.; Yurchenko, S.N.; Tennyson, J.; Jensen, P. Rotational spectrum of SO3 and a theoretical evidence for the formation of rotational, energy level clusters in its vibrational ground state. J. Chem. Phys. 2014, 140, 244316. [Google Scholar] [CrossRef] [PubMed]
  12. Underwood, D.S.; Tennyson, J.; Yurchenko, S.N.; Clausen, S.; Fateev, A. ExoMol line lists XVII: A line list for hot SO3. Mon. Not. R. Astron. Soc. 2016, 462, 4300–4313. [Google Scholar] [CrossRef]
  13. Al-Refaie, A.F.; Yurchenko, S.N.; Yachmenev, A.; Tennyson, J. ExoMol line lists-VIII: A variationally computed line list for hot formaldehyde. Mon. Not. R. Astron. Soc. 2015, 448, 1704–1714. [Google Scholar] [CrossRef]
  14. Owens, A.; Yurchenko, S.N.; Yachmenev, A.; Tennyson, J.; Thiel, W. Accurate ab initio vibrational energies of methyl chloride. J. Chem. Phys. 2015, 142, 244306. [Google Scholar] [CrossRef] [PubMed]
  15. Owens, A.; Yurchenko, S.N.; Yachmenev, A.; Thiel, W.; Tennyson, J. ExoMol molecular line lists XXII. The rotation-vibration spectrum of silane up to 1200 K. Mon. Not. R. Astron. Soc. 2017, 471, 5025–5032. [Google Scholar] [CrossRef]
  16. Al-Refaie, A.F.; Ovsyannikov, R.I.; Polyansky, O.L.; Yurchenko, S.N.; Tennyson, J. A variationally calculated room temperature line-list for H2O2. J. Mol. Spectrosc. 2015, 318, 84–90. [Google Scholar] [CrossRef]
  17. Al-Refaie, A.F.; Polyansky, O.L.; Ovsyannikov, R.I.; Tennyson, J.; Yurchenko, S.N. ExoMol line lists XV: A hot line-list for hydrogen peroxide. Mon. Not. R. Astron. Soc. 2016, 461, 1012–1022. [Google Scholar] [CrossRef]
  18. Chubb, K.L.; Yachmenev, A.; Tennyson, J.; Yurchenko, S.N. TROVE: Treating linear molecule HCCH. J. Chem. Phys. 2018. submitted. [Google Scholar]
  19. Mant, B.P.; Yachmenev, A.; Tennyson, J.; Yurchenko, S.N. ExoMol molecular line lists-XXVII: Spectra of C2H4. Mon. Not. R. Astron. Soc. 2018. submitted for publication. [Google Scholar]
  20. Tennyson, J.; Yurchenko, S.N. ExoMol: Molecular line lists for exoplanet and other atmospheres. Mon. Not. R. Astron. Soc. 2012, 425, 21–33. [Google Scholar] [CrossRef]
  21. Tennyson, J.; Yurchenko, S.N.; Al-Refaie, A.F.; Barton, E.J.; Chubb, K.L.; Coles, P.A.; Diamantopoulou, S.; Gorman, M.N.; Hill, C.; Lam, A.Z.; et al. The ExoMol database: molecular line lists for exoplanet and other hot atmospheres. J. Mol. Spectrosc. 2016, 327, 73–94. [Google Scholar] [CrossRef]
  22. Yurchenko, S.N.; Yachmenev, A.; Ovsyannikov, R.I. Symmetry-Adapted Ro-vibrational Basis Functions for Variational Nuclear Motion Calculations: TROVE Approach. J. Chem. Theory Comput. 2017, 13, 4368–4381. [Google Scholar] [CrossRef] [PubMed]
  23. Laane, J.; Ocola, E.J. Applications of Symmetry and Group Theory for the Investigation of Molecular Vibrations. Acta Appl. Math. 2012, 118, 3–24. [Google Scholar] [CrossRef]
  24. Fritzsche, S. Application of point-group symmetries in chemistry and physics: A computer-algebraic approach. Int. J. Quantum Chem. 2005, 106, 98–129. [Google Scholar] [CrossRef]
  25. Häser, M. Molecular point-group symmetry in electronic structure calculations. J. Chem. Phys. 1991, 95, 8259–8265. [Google Scholar] [CrossRef]
  26. Papoušek, D.; Aliev, M.R. Molecular Vibrational-Rotational Spectra: Theory and Applications of High Resolution Infrared; Studies in Physical and Theoretical Chemistry; Elsevier: Amsterdam, The Netherlands, 1982. [Google Scholar]
  27. Watson, J.K.G. Vibration-rotation hamiltonian of linear molecules. Mol. Phys. 1970, 19, 465–487. [Google Scholar] [CrossRef]
  28. Wlodarczak, G. Linear Polyatomic Molecules: Introduction. In Linear Polyatomic Molecules; Springer: Berlin/Heidelberg, Germnay, 2012; pp. 6–24. [Google Scholar]
  29. Chen, C. Symmetry Adapted Analysis of Linear Molecules. J. Chin. Chem. Soc. 1973, 20, 191–202. [Google Scholar] [CrossRef]
  30. Hegelund, F.; Rasmussen, F.; Brodersen, S. The selection rules and the transition moment for rotation–vibrational transitions in axial molecules. J. Raman Spectrosc. 1973, 1, 433–453. [Google Scholar] [CrossRef]
  31. Hirano, T.; Nagashima, U.; Jensen, P. Bending wavefunctions for linear molecules. J. Mol. Spectrosc. 2018, 343, 54–61. [Google Scholar] [CrossRef]
  32. Chubb, K.L.; Joseph, M.; Franklin, J.; Choudhury, N.; Furtenbacher, T.; Császár, A.G.; Gaspard, G.; Oguoko, P.; Kelly, A.; Yurchenko, S.N.; et al. MARVEL analysis of the measured high-resolution spectra of C2H2. J. Quant. Spectrosc. Radiat. Transf. 2018, 204, 42–55. [Google Scholar] [CrossRef]
  33. Brown, J.M.; Hougen, J.T.; Huber, K.P.; Johns, J.W.C.; Kopp, I.; Lefebvre-Brion, H.; Merer, A.J.; Ramsay, D.A.; Rostas, J.; Zare, R.N. The labeling of parity doublet levels in linear molecules. J. Mol. Spectrosc. 1975, 55, 500–503. [Google Scholar] [CrossRef]
  34. Herman, M.; Lievin, J. Acetylene—From intensity alternation in spectra to ortho and para molecule. J. Chem. Educ. 1982, 59, 17. [Google Scholar] [CrossRef]
  35. Jensen, P.; Hegelund, F. Lecture Notes: Molecular Rotation-Vibration Theory; ResearchGate: Gatersleben, Germany, 2014; pp. 1–97. [Google Scholar] [CrossRef]
  36. Schutte, C.J.H.; Bertie, J.E.; Bunker, P.R.; Hougen, J.T.; Mills, I.M.; Watson, J.K.G.; Winnewisser, B.P. Notations and conventions in molecular spectroscopy: Part 2. Symmetry notation. Pure Appl. Chem. 1997, 69, 1641–1650. [Google Scholar] [CrossRef]
  37. Noumerov, B.V. A method of extrapolation of perturbations. Mon. Not. R. Astron. Soc. 1924, 84, 592–602. [Google Scholar] [CrossRef]
  38. Cooley, J.W. An Improved eigenvalue corrector formula for solving the Schrödinger equation for central fields. Math. Comp. 1961, 15, 363–374. [Google Scholar]
  39. Sørensen, G.O. A New Approach to the Hamiltonian of Nonrigid Molecules. In Large Amplitude Motion in Molecules II; Topics in Current Chemistry; Dewar, M.J.S., Ed.; Springer: Berlin/Heidelberg, Germeny, 1979; Volume 82, pp. 97–175. [Google Scholar]
  40. Bunker, P.R.; Jensen, P. Spherical top molecules and the molecular symmetry group. Mol. Phys. 1999, 97, 255–264. [Google Scholar] [CrossRef]
  41. Yurchenko, S.N.; Carvajal, M.; Jensen, P.; Lin, H.; Zheng, J.J.; Thiel, W. Rotation-vibration motion of pyramidal XY3 molecules described in the Eckart frame: Theory and application to NH3. Mol. Phys. 2005, 103, 359–378. [Google Scholar] [CrossRef]
Figure 1. The block diagonal structure of a Hamiltonian matrix in the D n h irreducible representation. The empty (white) cells indicate blocks of vanishing matrix elements. It should be noted that, although B-symmetries will be present for even values of n, they are not physical and do not appear as a block of matrix elements to be diagonalized.
Figure 1. The block diagonal structure of a Hamiltonian matrix in the D n h irreducible representation. The empty (white) cells indicate blocks of vanishing matrix elements. It should be noted that, although B-symmetries will be present for even values of n, they are not physical and do not appear as a block of matrix elements to be diagonalized.
Symmetry 10 00137 g001
Figure 2. HCCH as described using the ( 3 N 5 ) coordinates employed in TROVE. R is the vector (of length R) pointing from the first to the second carbon atom, C 1 to C 2 , while r i are the two CH i bond vectors (of lengths r i ). The Δ x 1 , Δ x 2 , Δ y 1 and Δ y 2 notation of this diagram is to reflect the Cartesian projections of the CH i bond vectors.
Figure 2. HCCH as described using the ( 3 N 5 ) coordinates employed in TROVE. R is the vector (of length R) pointing from the first to the second carbon atom, C 1 to C 2 , while r i are the two CH i bond vectors (of lengths r i ). The Δ x 1 , Δ x 2 , Δ y 1 and Δ y 2 notation of this diagram is to reflect the Cartesian projections of the CH i bond vectors.
Symmetry 10 00137 g002
Table 1. Common character table for the point group D h and the Extended Molecular Symmetry (EMS) group D h (EM) a .
Table 1. Common character table for the point group D h and the Extended Molecular Symmetry (EMS) group D h (EM) a .
D h (EM): E 0 E ε E ε * ( 12 ) π * ( 12 ) π + ε * ( 12 ) ε
1212
D h :E 2 C ε σ v ( ε / 2 ) i 2 S π + ε C 2 ( ε / 2 )
Σ g + , A 1 g :111111
Σ u + , A 2 u :111 1 1 1
Σ g , A 2 g :11 1 11 1
Σ u , A 1 u :11 1 1 1 1
Π g , E 1 g :2 2 cos ε 2 02 2 cos ε 2 0
Π u , E 1 u :2 2 cos ε 2 0 2 2 cos ε 2 0
Δ g , E 2 g :2 2 cos 2 ε 02 2 cos 2 ε 0
Δ u , E 2 u :2 2 cos 2 ε 0 2 2 cos 2 ε 0
Φ g , E 3 g :2 2 cos 3 ε 02 2 cos 3 ε 0
Φ u , E 3 u :2 2 cos 3 ε 0 2 2 cos 3 ε 0
a The elements of D h are defined as follows: C ε is a rotation by ε about the molecular axis, σ v ( ε / 2 ) is a reflection in a plane containing the molecular axis, i is the point group inversion operation, S π + ε is an improper rotation by π + ε about the molecular axis, and C 2 ( ε / 2 ) is a rotation by π about an axis perpendicular to the molecular axis (see also Ref. [1]). Here, ε = 0 2 π . See the text for the definitions of the D h (EM) operations.
Table 2. Character table for the D h (M) Molecular Symmetry (MS) group a .
Table 2. Character table for the D h (M) Molecular Symmetry (MS) group a .
Γ 1 Γ 2 Γ 3 Γ 4 Γ 5 Γ 6 E ( p ) E * ( p ) *
Σ g + + s A 1 A + A g A 1 g 1111
Σ u + + a B 2 B + B u A 2 u 1 1 1 1
Σ g a B 1 B B g A 2 g 1 1 1 1
Σ u s A 2 A A u A 1 u 11 1 1
a Γ 1 Γ 6 are several alternative notations for the irreducible representations of D h (M). Γ 3 is customarily used for C 2 v (M) and Γ 5 is for C 2 h (M) (Table A-8 of Ref. [1]). g and u stand for the German gerade (even) and ungerade (odd), related to the permutation-inversion operation ( p ) * .
Table 3. The irreducible representation Γ of D h (EM) spanned by the rotational wavefunction | J , k of a linear molecule in the absence of external electric and magnetic fields. The irrep depends on k, the z-axis projection in units of , of the rotational angular momentum.
Table 3. The irreducible representation Γ of D h (EM) spanned by the rotational wavefunction | J , k of a linear molecule in the absence of external electric and magnetic fields. The irrep depends on k, the z-axis projection in units of , of the rotational angular momentum.
k Γ
0(J even) Σ g +
(J odd) Σ g
±1 Π g
±2 Δ g
±3 Φ g
Table 4. Symmetry labels for the ro-vibrational states of a linear molecule such as 12 C 2 H 2 . The e / f labels are defined in Ref. [33] and o r t h o / p a r a define the nuclear-spin state [32,34]. Γ a , Γ b and Γ c are alternative notations for the irreducible representations of D h (M) (see Table 2 for an expanded list).
Table 4. Symmetry labels for the ro-vibrational states of a linear molecule such as 12 C 2 H 2 . The e / f labels are defined in Ref. [33] and o r t h o / p a r a define the nuclear-spin state [32,34]. Γ a , Γ b and Γ c are alternative notations for the irreducible representations of D h (M) (see Table 2 for an expanded list).
Γ a Γ b Γ c e / f Ortho / Para
J odd: Σ g + A 1 g + s f p a r a
Σ u A 1 u s e p a r a
Σ g A 2 g a e o r t h o
Σ u + A 2 u + a f o r t h o
J even: Σ g + A 1 g + s e p a r a
Σ u A 1 u s f p a r a
Σ g A 2 g a f o r t h o
Σ u + A 2 u + a e o r t h o
Table 5. Symmetry operations of the D n h groups, for even and odd n. σ h , σ v and σ d represent reflections in planes perpendicular to the molecular axis, containing the molecular axis, and bisecting the angle between a pair of C 2 axes, respectively. An improper rotation S n r is a rotation by r ( 2 π n ) ( r = 1 n 2 ) followed by a reflection in the plane perpendicular to the molecular axis and containing the nuclear center-of-mass. C n r represents rotations by r ( 2 π n ) , where r = 1 n 1 . See Ref. [1] for further details on these symmetry operations.
Table 5. Symmetry operations of the D n h groups, for even and odd n. σ h , σ v and σ d represent reflections in planes perpendicular to the molecular axis, containing the molecular axis, and bisecting the angle between a pair of C 2 axes, respectively. An improper rotation S n r is a rotation by r ( 2 π n ) ( r = 1 n 2 ) followed by a reflection in the plane perpendicular to the molecular axis and containing the nuclear center-of-mass. C n r represents rotations by r ( 2 π n ) , where r = 1 n 1 . See Ref. [1] for further details on these symmetry operations.
Symmetry OperationNumber of OperationsDescription
Even n:
E1Identity
C n r n 1 Rotations about the n-fold molecular axis
C 2 / C 2 nn rotations by π about axes
perpendicular to the molecular axis
i1Point group inversion
S n r n 2 Improper rotation (see caption)
σ h 1Horizontal reflection (see caption)
σ v n / 2 Vertical reflection (see caption)
σ d n / 2 Diagonal reflection (see caption)
Total: 4 n
Odd n:
E1Identity
C n r n 1 Rotations about the n-fold molecular axis
C 2 nn rotations by π about axes
perpendicular to the molecular axis
S n r n 1 Improper rotation (see caption)
σ h 1Horizontal reflection (see caption)
σ v nVertical reflection (see caption)
Total: 4 n
Table 6. Generating operations for the D n h groups (n even and n odd).
Table 6. Generating operations for the D n h groups (n even and n odd).
Point Group R + R + R
D n h , n odd C n σ h C 2 ( x )
D n h , n even C n i C 2 ( x )
Table 7. Irreducible representations for the D n h groups and their characters under the generating operations R + , R + and R .
Table 7. Irreducible representations for the D n h groups and their characters under the generating operations R + , R + and R .
D n h
(n even)
E R +
(= C n )
R +
(=i)
R
(= C 2 ( x ) )
A 1 g 1−111
A 2 g 1−11−1
B 1 g 1−111
B 2 g 1−11−1
E r g a 2 2 cos 2 π r n 20
A 1 u 1−1−11
A 2 u 1−1−1−1
B 1 u 1−1−11
B 2 u 1−1−1−1
E r u a 2 2 cos 2 π r n −20
D n h
(n odd)
E R +
(= C n )
R +
(= σ h )
R
(= C 2 ( x ) )
A 1 1111
A 2 111−1
E r b 2 2 cos 2 π r n 20
A 1 11−11
A 2 11−1−1
E r b 2 2 cos 2 π r n −20
a r = 1, 2, …, n 2 1 ; b r = 1, 2, …, n 1 2 .
Table 8. The correspondence between the g / u (gerade/ungerade) notation of the irreps of D n h (even n) and the / notation of the irreps of D n h (odd n), based on K (the absolute value of the projection, in units of , onto the molecule-fixed z-axis of the rotational angular momentum).
Table 8. The correspondence between the g / u (gerade/ungerade) notation of the irreps of D n h (even n) and the / notation of the irreps of D n h (odd n), based on K (the absolute value of the projection, in units of , onto the molecule-fixed z-axis of the rotational angular momentum).
K Γ (even n) Γ (odd n)D h (EM)
0 A 1 g A 1 Σ g +
A 1 u A 1 Σ u +
A 2 g A 2 Σ g
A 2 u A 2 Σ u
>0, even E k g E k Δ g , Γ g , I g
E k u E k Δ u , Γ u , I u
>0, odd E k g E k Π g , Φ g , H g
E k u E k Π u , Φ u , H u
Table 9. Transformation matrices for the D n h groups generated by the rotational basis functions | J , 0 , + = | J , 0 for K = 0 and ( | J , K , | J , K ) and ( | J , K , + , | J , K , ) for K > 0, with ε = 2 π / n .
Table 9. Transformation matrices for the D n h groups generated by the rotational basis functions | J , 0 , + = | J , 0 for K = 0 and ( | J , K , | J , K ) and ( | J , K , + , | J , K , ) for K > 0, with ε = 2 π / n .
D n h (n even)E R + = C n R + = i R = C 2 ( x )
| J , 0 , + 1111
| J , K | J , K 1 0 0 1 e + i K ε 0 0 e i K ε 1 0 0 1 0 1 1 0
| J , K , + | J , K , 1 0 0 1 cos K ε sin K ε sin K ε cos K ε 1 0 0 1 1 0 0 1
D n h ( n odd) E R + = C n R + = σ h R = C 2 ( x )
| J , 0 , + 111 ( 1 ) J
| J , K | J , K 1 0 0 1 e + i K ε 0 0 e i K ε ( 1 ) K 0 0 ( 1 ) K 0 ( 1 ) J ( 1 ) J 0
| J , K , + | J , K , 1 0 0 1 cos K ε sin K ε sin K ε cos K ε ( 1 ) K 0 0 ( 1 ) K ( 1 ) J 0 0 ( 1 ) J
Table 10. Irreducible-representation transformation matrices of the D n h group for n even, generated by the rotational basis functions ( | J , K , + , | J , K , ) for K> 0. ε = 2 π n , r is an integer used to identify the group operations, and κ = | K + n t | ; the integer t is determined such that 1 ⩽ κ n / 2 1 .
Table 10. Irreducible-representation transformation matrices of the D n h group for n even, generated by the rotational basis functions ( | J , K , + , | J , K , ) for K> 0. ε = 2 π n , r is an integer used to identify the group operations, and κ = | K + n t | ; the integer t is determined such that 1 ⩽ κ n / 2 1 .
ε r r E κ g E κ u
E 1 0 0 1 1 0 0 1
C n r r K ε 1 n 1 cos ε r sin ε r sin ε r cos ε r cos ε r sin ε r sin ε r cos ε r
C 2 2 r K ε 0 n 2 1 1 K cos ε r sin ε r sin ε r cos ε r 1 K cos ε r sin ε r sin ε r cos ε r
C 2 ( 2 r + 1 ) K ε 0 n 2 1 1 K cos ε r sin ε r sin ε r cos ε r 1 K cos ε r sin ε r sin ε r cos ε r
i 1 0 0 1 1 0 0 1
σ h ( 1 ) K 1 0 0 1 ( 1 ) K 1 0 0 1
σ v 2 r K ε 0 n 2 1 1 K cos ε r sin ε r sin ε r cos ε r 1 K cos ε r sin ε r sin ε r cos ε r
σ d ( 2 r + 1 ) K ε 0 t n 2 1 1 K cos ε r sin ε r sin ε r cos ε r 1 K cos ε r sin ε r sin ε r cos ε r
S n ( r ) r K ε 1, 2, 3, …, n / 2 1 ,
n / 2 + 1 , …, n 1 a
1 K cos ε r sin ε r sin ε r cos ε r 1 K cos ε r sin ε r sin ε r cos ε r
a We omit r = 0 and r = n / 2 from this list because S n ( 0 ) = σ h and S n ( n / 2 ) = i.
Table 11. Irreducible-representation transformation matrices of the D n h group for n odd, generated by the rotational basis functions ( | J , K , + , | J , K , ) for K > 0. a ε = 2 π n , r is an integer used to identify the group operations, and κ = | K + n t | ; the integer t is determined such that 1 ⩽ κ ( n 1 ) / 2 .
Table 11. Irreducible-representation transformation matrices of the D n h group for n odd, generated by the rotational basis functions ( | J , K , + , | J , K , ) for K > 0. a ε = 2 π n , r is an integer used to identify the group operations, and κ = | K + n t | ; the integer t is determined such that 1 ⩽ κ ( n 1 ) / 2 .
ε r r E κ E κ
E 1 0 0 1 1 0 0 1
C n r r K ε 1 n 1 cos ε r sin ε r sin ε r cos ε r cos ε r sin ε r sin ε r cos ε r
C 2 r K ε 0 n 1 1 J cos ε r sin ε r sin ε r cos ε r 1 J cos ε r sin ε r sin ε r cos ε r
σ h 1 0 0 1 1 0 0 1
σ v r K ε 0 n 1 1 J cos ε r sin ε r sin ε r cos ε r 1 J cos ε r sin ε r sin ε r cos ε r
S n ( r ) r K ε 1 n 1 cos ε r sin ε r sin ε r cos ε r cos ε r sin ε r sin ε r cos ε r
a E κ ( E κ ) functions have even(odd) K (Table 12). The present table has been simplified accordingly.
Table 12. Irreducible representations of the D n h groups generated by the rotational basis functions | J , 0 for K = 0 and ( | J , K , + , | J , K , ) for K > 0. The integer t is determined such that 0 ⩽ κ n / 2 for n even and 0 ⩽ κ ( n 1 ) / 2 for n odd, with κ = | K + n t | .
Table 12. Irreducible representations of the D n h groups generated by the rotational basis functions | J , 0 for K = 0 and ( | J , K , + , | J , K , ) for K > 0. The integer t is determined such that 0 ⩽ κ n / 2 for n even and 0 ⩽ κ ( n 1 ) / 2 for n odd, with κ = | K + n t | .
K κ D n h (n even)D h (EM)
=00 A 1 g Σ g +
>00 A 1 g A 2 g Σ g + Σ g +
>0 n / 2 B 1 g B 2 g Σ g + Σ g +
>0 κ = 1, 2, …, n / 2 1 E κ g Π g , Δ g , Φ g , Γ g , …
K κ D n h ( n odd)D h (EM)
=00J even A 1 Σ g +
J odd A 2 Σ g
>0, odd0 A 1 A 2 Σ u + Σ u
>0, even0 A 1 A 2 Σ g + Σ g
>0, odd1, 2, …, ( n 1 ) / 2 E κ Π g , Φ g , H g
>0, even1, 2, …, ( n 1 ) / 2 E κ Δ u , Γ u , I u
Table 13. Transformation properties based on those of Table 10 for the D n h group (relating to the symmetry operations of Table 5), where n is even, for transforming the set of 7 vibrational coordinates ( Δ R , Δ r 1 , Δ r 2 , Δ x 1 , Δ y 1 , Δ x 2 , Δ y 2 ) used in the calculations of Ref. [18] for linear molecule 12 C 2 H 2 , as illustrated in Figure 2. The two-component vectors ρ 1 = ( Δ x 1 , Δ y 1 ) T and ρ 2 = ( Δ x 2 , Δ y 2 ) T transform as E 1 u , with the transformation matrices M R E 1 u from Table 10. m is an integer for the bounds given for each operation, used to form ε m , where ε = 2 π n in all cases.
Table 13. Transformation properties based on those of Table 10 for the D n h group (relating to the symmetry operations of Table 5), where n is even, for transforming the set of 7 vibrational coordinates ( Δ R , Δ r 1 , Δ r 2 , Δ x 1 , Δ y 1 , Δ x 2 , Δ y 2 ) used in the calculations of Ref. [18] for linear molecule 12 C 2 H 2 , as illustrated in Figure 2. The two-component vectors ρ 1 = ( Δ x 1 , Δ y 1 ) T and ρ 2 = ( Δ x 2 , Δ y 2 ) T transform as E 1 u , with the transformation matrices M R E 1 u from Table 10. m is an integer for the bounds given for each operation, used to form ε m , where ε = 2 π n in all cases.
Irrep ε m mTransformation
E Δ R
Δ r 1
Δ r 2
ρ 1
ρ 2
C n m m ε 1 n 1 Δ R
Δ r 1
Δ r 2
M C n r E 1 u ( ε m ) · ρ 1
M C n r E 1 u ( ε m ) · ρ 2
C 2 2 m ε 0 n 2 1 Δ R
C 2 ε ( 2 m + 1 ) 0 n 2 1 Δ r 2
Δ r 1
M C 2 / E 1 u ( ε m ) · ρ 2
M C 2 / E 1 u ( ε m ) · ρ 1
i Δ R
Δ r 2
Δ r 1
ρ 2
ρ 1
σ h Δ R
Δ r 2
Δ r 1
ρ 2
ρ 1
σ d ε ( 2 m + 1 ) 0 n 2 1 Δ R
σ v 2 m ε 0 n 2 1 Δ r 1
Δ r 2
M σ d / v E 1 u ( ε m ) · ρ 1
M σ d / v E 1 u ( ε m ) · ρ 2
S n m m ε 1 n 2 Δ R
Δ r 2
Δ r 1
M S n r E 1 u ( ε m ) · ρ 2
M S n r E 1 u ( ε m ) · ρ 1
Table 14. Symmetries of the symmetrized rotational basis set used by TROVE, Equations (31) and (32) for different combinations of J, K and τ (where τ ( = 0 , 1 ) and K = | k | ); each 2D representation E K g state has an a and b component, represented by the different values of τ . See Table 8 for an explanation of the differing notation of Γ r o t for even and odd values of n.
Table 14. Symmetries of the symmetrized rotational basis set used by TROVE, Equations (31) and (32) for different combinations of J, K and τ (where τ ( = 0 , 1 ) and K = | k | ); each 2D representation E K g state has an a and b component, represented by the different values of τ . See Table 8 for an explanation of the differing notation of Γ r o t for even and odd values of n.
K τ Γ rot
Even nOdd n
00 A 1 g A 1
1 A 2 g A 2
>0, odd0 E k g b E k b
1 E k g a E k a
>0, even0 E k g a E k a
1 E k g b E k b
Table 15. An example of some rotational, vibrational and ro-vibrational assignments (see Section 4 for the meaning of the rotational assignments and e.g., [32] for the vibrational assignments) with associated symmetries ( Γ r , Γ v and Γ r v , respectively) from ro-vibrational calculations using TROVE of 12 C 2 H 2 using different (even/odd) values of n for D n h . In each case L max = 4. The energies are identical for symmetries of higher n than those shown here, but converge towards the experimental values as the polyad number Equation (34) is increased; a low value is used here for demonstration purposes. The symmetry assignment will remain unchanged for more accurate calculations; these will be published elsewhere.
Table 15. An example of some rotational, vibrational and ro-vibrational assignments (see Section 4 for the meaning of the rotational assignments and e.g., [32] for the vibrational assignments) with associated symmetries ( Γ r , Γ v and Γ r v , respectively) from ro-vibrational calculations using TROVE of 12 C 2 H 2 using different (even/odd) values of n for D n h . In each case L max = 4. The energies are identical for symmetries of higher n than those shown here, but converge towards the experimental values as the polyad number Equation (34) is increased; a low value is used here for demonstration purposes. The symmetry assignment will remain unchanged for more accurate calculations; these will be published elsewhere.
Exp. Energy (cm 1 ) [32]Energy (cm 1 )JK τ ν 1 ν 2 ν 3 ν 4 l 4 ν 5 l 5 D 12 h D 13 h
This Work Γ r v Γ r Γ v Γ r v Γ r Γ v
2.3532862.356491101 0000 0 0 0 A 2 g A 2 g A 1 g A 2 A 2 A 1
614.044355625.810547111 0001 1 0 0 A 2 g E 1 g E 1 g A 2 E 1 E 1
1232.7491621283.603736101 0002 0 0 0 A 2 g A 2 g A 1 g A 2 A 2 A 1
7.0598227.069433200 0000 0 0 0 A 1 g A 1 g A 1 g A 1 A 1 A 1
618.745653630.518518211 0001 1 0 0 A 1 g E 1 g E 1 g A 1 E 1 E 1
1235.8743921276.518756220 0002 2 0 0 A 1 g E 2 g E 2 g A 1 E 2 E 2

Share and Cite

MDPI and ACS Style

Chubb, K.L.; Jensen, P.; Yurchenko, S.N. Symmetry Adaptation of the Rotation-Vibration Theory for Linear Molecules. Symmetry 2018, 10, 137. https://doi.org/10.3390/sym10050137

AMA Style

Chubb KL, Jensen P, Yurchenko SN. Symmetry Adaptation of the Rotation-Vibration Theory for Linear Molecules. Symmetry. 2018; 10(5):137. https://doi.org/10.3390/sym10050137

Chicago/Turabian Style

Chubb, Katy L., Per Jensen, and Sergei N. Yurchenko. 2018. "Symmetry Adaptation of the Rotation-Vibration Theory for Linear Molecules" Symmetry 10, no. 5: 137. https://doi.org/10.3390/sym10050137

APA Style

Chubb, K. L., Jensen, P., & Yurchenko, S. N. (2018). Symmetry Adaptation of the Rotation-Vibration Theory for Linear Molecules. Symmetry, 10(5), 137. https://doi.org/10.3390/sym10050137

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop