Services on Demand
Article
Indicators
 Cited by SciELO
Related links
 Cited by Google
 Similars in SciELO
 Similars in Google
Share
Boletín de la Sociedad Chilena de Química
Print version ISSN 03661644
Bol. Soc. Chil. Quím. vol.47 no.4 Concepción Dec. 2002
http://dx.doi.org/10.4067/S036616442002000400003
TORSIONAL VIBRATION AND INTERNAL ROTATION IN (XYn)2
MOLECULES. THE HÜCKEL  MÖBIUS LCVW.
Jorge Ricardo Letelier D.*a , Alejandro ToroLabbéb
and
YingNan Chiuc
aDepartamento de Química, Facultad de Ciencias Físicas y Matemáticas
Universidad de Chile, Ave. Tupper 2069, Casilla 2777, Santiago, Chile.
Email jletelie@tamarugo.cec.uchile.cl
bFacultad de Química, Universidad Católica de Chile.
Vicuña Mackenna 4860, Santiago, Chile.
Email atola@puc.cl
cDepartment of Chemistry, The Catholic University ofAmerica
Washington, DC 20064, USA
(Received: March 06, 2002  Accepted: June 10, 2002)
ABSTRACT
The geometrical foundation of the HückelMöbius concept is revisited and used to represent approximate torsional wavefunctions of the hindered internal rotation of a coaxial (XYn)2 type molecule, constructing linear combinations of vibrational wavefunctions (LCVW) centered at the periodic potential minima.
Making use of the sign variation of the "resonance" interaction integral with different vibrational levels, we construct energy levels that form the cluster of properly alternating nondegenerate and degenerate levels for the highbarrier case.
Keywords: Torsional Vibration, Hindered rotation eigenvalues
RESUMEN
Se revisan los fundamentos geométricos del concepto HückelMöbius utilizado para representar funciones de onda torsionales aproximadas para la rotación interna impedida de moléculas tipo (XYn)2 coaxiales, construyendo combinaciones lineales de funciones vibracionales (LCVW) centradas en lo mínimos del potencial periódico.
Haciendo uso de la variación de signo de la integral de interacción "resonante", construimos el patrón de niveles de energía que conforman el cluster de niveles de energía alternantes, degenerados y no degenerados, para el caso de una barrera alta.
PALABRAS CLAVES: Vibración torcional, autovalores, totación impedida.
INTRODUCTION
Hindered internal rotation has been of interest for a long time since its discovery in ethane by Kemp and Pitzer^{1} and its early treatment by Pitzer and by Wilson^{2}. It has been thoroughly reviewed from the point of view of microwave spectroscopy^{3}. Comprehensive symmetry treatments generally make use of LonguetHiggins' symmetry group of nonrigid molecules. These treatments were motivated by the work of LonguetHiggins^{4}, Bunker ^{5} and Hougen^{6} .These treatments use a molecular double point group and arrive at symmetry classifications of vibrational and rotational levels. Specifically, for hindered internal rotation of ethane, hydrazine, ethylene and (XY_{2})_{2} coaxial rotors^{7}, these treatments lead to the cluster of energy levels with appropriate alternation of nondegenerate (A) and degenerate (E) levels required by symmetry as shown in the work of Harter and Patterson^{8}.
The general pattern and energy clusters for two coaxial rotors of equal groups (XYn)_{2}, that makes use of the HückelMöbius concept has been derived and thoroughly discussed by Chiu^{9}. In this work, we construct torsionalvibrational wavefunctions by taking appropriate linear combination vibrational wavefunctions (LCVW) centered at potential minima of the hindered (torsional) rotator with the purpose of approximating these torsional wavefunctions at high barrier. The process is analogous to taking linear combinations of atomic orbitals (LCAO) to construct molecular orbitals. Chiu went further to show that these LCVW have the desired transformation properties under the appropriate operations of the double group of Bunker and Hougen.
METHOD
We briefly review here the basis of the method; detailed analysis can be found elsewhere^{9}. Following Hougen6 and Bunker^{5} we consider the combined rotational wavefunctions with respect to the common axis of the two groups c_{a} and c_{b}, with rotational angles c a and c b respectively.
(1) 
Where K = ka + kb is the angular momentum for the overall rotation of the whole molecule around the common axis, and k = (ka  kb) is for the free relative internal rotation of the two groups during which the overall momentofinertia axes are not changed. When the barrier is high this free relative internal rotation becomes torsional vibration. Because of the 1/2 factor in the definition of the overall rotational angles and the nfold rotation of either the a or b group of (XYn), Cn, will produce half of the expected (2p /n) change in c and g . This is the reason why the potential for internal rotation has 2nfold rather than nfold symmetry. Furthermore the ntimes repeated nfold rotation,, of a or b, while regenerating the same physical configuration of the molecule, may or may not lead to the same F rot or F tor depending on whether K or k is even or odd, viz.
(2) 
The doublevalued nature or , is the basis for the doublegroup treatment. However, because the rotational wavefunction as a whole must be invariant under either or we see immediately that the even must occur with even (and odd K^{o} with odd K^{o}). This is so that the change c and g to c + p and g + p will leave the total rotational wavefunction Y invariant which is a physical basis6 of the double group. In other words, two orientations are equivalent, and the number of symmetry operators is doubled,
viz.
(3a) 
(3b) 
where the combinations and will always be even, guaranteeing invariance of Y in (3b) under two orientations.
With the above as background, the LCVW torsional wavefunctions that we construct must have those with odd property combining with odd K (= 1, 3 ...) and those with even combining with even K (= 0, 2 ...). It is found that the Hückel combination will always have even property and Möbius combination will always have odd property and there is just the right number of each. To show this we recall that because of the mathematical definition of g = 1/2(c a  c b), there are 2n minima for the nfold symmetry group XYn, in coaxial(XYn)_{2}. By simple cyclic group theory, there are 2n linear combinations, notoriously of the Hückel type, of the vibrational wavefunctions with appropriate energies similar to LCAOMO for periodic cyclic systems, as follows
Where k = 0, 1, 2, ... (2n 1). To reach equation (4b), it has been assumed only nearest neighbor interaction, where b is related to the transmission coefficient for tunneling from one potential minimum to the next. It is the analog of the resonance integral in LCAO theory. We then divide the linear combinations into two sets: one with even k e(= 2^{L} ) and one with odd k o(= 2^{L '} + 1). We retrieve the nfold symmetry in each set before we apply the or , operator to test their Hückel or Möbius property, viz.
(Hückel) k ^{e}= 2L  (5) 
and where L' are the integers:
(Möbius) k ^{o}= 2L ' + 1  (6) 
Where the superscripts o and e denote odd and even numbers respectively and where the nfold symmetry is apparent, based on summation index from 0 to n1 in the wavefunctions and based on the doublegroup characters .
APPLICATION TO (XYn)2 ¾ TYPE MOLECULES
The concepts and methods of calculation, developed in the previous section, are now applied specifically to the case of the (XYn)2 type molecules, where n = 2, 3, 4, 5.
In the foregoing analysis, the overall rotational angles and , as well as the angles c _{a} and c _{b} are referred to laboratory coordinates. In this context, g truly gives rise to a doublevalued problem. We simplify the analysis by choosing instead the internal (referred to molecular axes) relative torsional angle q and the problem becomes then a nfold symmetry instead of a 2nfold one. The torsional wavefunctions and energies are now:
The potential used has the general form
(8) 
Where , is scaling constant used to express energy in units of the rotational energy of a rotor of whose moment of inertia is , a_{0} being the unit of atomic length (ie. first Bohr radius). In this manner, the handling of potential functions and energy values are dimensionless. As an example, a typical cyclic potential of 4fold symmetry is shown in Figure 1.
Fig 1. The nfold potential energy function, as given by equation (8), for the relative torsional oscillation. The basis functions for the LCVW are taken as the solutions of the portion AB. 
In this calculation, the wavefunctions F _{s} of equation (7a) are taken as the solutions of a onedimensional oscillator, whose potential corresponds to the portion around one of the minimum of the cyclic potential.
Fig. 2. Wave functions, as given by equation (7a) for k =1, 3, 5 (odd) for a torsional oscillator of 2fold symmetry. Potential minima are found at q = p /2 and q = 3p /2. 
With reference to Figure 1, we solved the eigenvalue problem within the portion AB indicated there and because the potential function V(q ) in eq. (8) is far from harmonic, the solutions are found numerically using other numerical methods developed by us^{11}. The eigenfunctions are then of numerical nature, therefore, the linear combinations of used in eq. (7a) are also formed numerically. Within this approximation, only the lowest clusters of eigenvalues are considered to be reasonably well described.
Fig. 3.Wave functions, as given by equation (7a) for k =2, 4, 6 (even) for a torsional oscillator of 2fold symmetry. Potential minima are found at q = p /2 and q = 3p /2. 
Fig. 4.Wave functions, as given by equation (7a) for k =1, 3, 5 (odd) for a torsional oscillator of 3fold symmetry. Potential minima are found at q = p /4, p and q = 3p /4. 
Fig. 5.Wave functions, as given by equation (7a) for k =2, 4, 6 (even) for a torsional oscillator of 3fold symmetry. Potential minima are found at q = p /4, p and q = 3p /4 
Within a given cluster of eigenvalues, the parameter a is set equal to the corresponding eigenvalue (v = 1, 2, ...) of the local potential (as shown schematically in Figure 1). The cluster of eigenvalues then in eq. (7b) corresponds to the splitting due to interaction between the displaced torsional oscillators.
The "Resonance" Integral.
Due to the analogy of the LCVW method to LCAOMO it seems natural and tempting to use a similar approximation for the resonance integral b by making it proportional to the overlap. There is a difference however, which is that these are FranckCondon overlap integrals between displaced oscillators at neighboring minima (if nearest neighbor approximation is considered). This "resonance integral", which determines the splitting of torsional levels, is also related to the tunneling probability, for this reason Chiu^{9}, has coined them "tunneling interaction" integral.
Since the value of b depends on the width, as well as the height and shape of the barrier, there is no simple way of knowing nor estimating it a priori. The effect of overlap can be incorporated (and also the sign alternation that comes with it) by using tables of FranckCondon overlap integrals R_{vv} for displaced harmonic oscillator^{12} (HO), if effectively harmonic oscillators wavefunctions are used. In our case, since we used a potential of the form give in equation (8), the wavefunctions in eq. (7a) are not HO. Again, the FranckCondon overlaps were carried out numerically and b is made proportional to the corresponding R_{vv. }In our calculations, we have considered every quantity equal to unity, that is, torsional moment of inertia =1.0, etc. to make it system independent.
Results for the lowest eigenvalue cluster for several nfold symmetry potentials are given in Table I. In this table, the cluster pattern is readily observed. Their corresponding wavefunctions are depicted in figures 2 to 5 (for n = 2 and 3fold symmetry only). In these figures, the alternation pattern of doublesingle level is apparent.
Clusters of torsional vibration eigenvalues for the nfold multiplewell potential
for (XY_{n})_{2} type molecules ( n =2,3,4,5). W_{n} refers to the number of well minima.
(Energy is given in arbitrary (dimensionless) units.)
 
Eigenvalue  W_{2}  W_{3}  W_{4}  W_{5} 
 
1  4.27917  3.18521  2.29056  1.64666 
2  4.27842  3.17238  2.25023  1.58105 
3  0.32276  3.14330  2.13236  1.36303 
4  0.34513  2.82333  2.00432  1.05760 
5  4.01714  2.98497  4.50004  0.76804 
6  4.22409  3.40708  4.76873  5.97405 
7  6.56528  6.79841  5.55685  7.34790 

These results, although contain many approximations, qualitatively describe the energy level pattern and are in reasonably good agreement with those calculated by direct numerical integration^{11}.
Improvements of the method.
The energy values computed using this method can be improved by adding non nearest neighbor interactions, much in the same way as it is found in the treatment of linear crystals with periodic boundary conditions^{10}. This requires the inclusion of several types of "resonance integrals" b _{1}, b_{2}, b_{3}, appropriate to first, second, thirdneighbor interaction, and so on. Energies and wave functions would be as follow:
(9) 
As an extension, the method can also be adapted to include nonsymmetric cyclic potentials. As an example, consider the case of periodic cyclic potentials where local potential of different depth alternate, we have then
(10) 
Here F ^{(1)} and F ^{(2)} represent the vibrational wave functions of the two local potentials and only nearestneighbor interaction has been considered.
ACKNOWLEDGMENTS
Funding for this work by FondecytChile, Grant N° 1000971 and the ChiuFeng Chia Research Fund are gratefully acknowledged.
REFERENCES
[1] J. D. Kemp and K. S. Pitzer, J. Chem. Phys., 4, 749 (1936); J. D. Kemp and K. S. Pitzer, J. Am. Chem. Sec., 59, 276 (1937); K. S. Pitzer, Chem. Rev., 27, 39 (1940). [ Links ] [ Links ] [ Links ]
[2] E. B. Wilson, Jr., J. Chem. Phys., 6, 740 (1938). [ Links ]
[3] a) C. H. Townes and A. K. Schawlow, Microwave Spectroscopy, McGrawHill, New York, 1955; b) H. C. Alien, Jr. and P. C. Cross, Molecular VibRotors, Wiley, New York, 1963; c) T, M. Sugden and C. N. Keuney, Microwave Spectroscopy of Gases, Van Nostrand, London, 1965; d) J. E. Wollrab, Rotational Spectra and Molecular Structure, Academic Press, New York, 1967; e) W. Cordy, William V. Smith and R. F. Trambarulo, Microwave Spectroscopy, Wiley, New York, 1953; f) W. Cordy and R. L. Cook, Microwave Molecular Spectre, Interscience, New York, 1970. [ Links ] [ Links ] [ Links ] [ Links ] [ Links ] [ Links ]
[4] H. C. LonguetHiggins, Mol. Phys., 6, 445 (1963). [ Links ]
[5] a) P. R. Bunker and H. C. LonguetHiggins, Proc. Roy. Soc. London, Ser.A 280, 340 (1964); b ) P. R. Bunker, Mol. Phys., 8, 81 (1964). [ Links ] [ Links ]
[6] J. T. Hougen, Can. J. Phys., 42, 1920 (1964).; J. T. Hougen, Can. J. Phys., 43, 935 (1965); J. T. Hougen, Pure Appl. Chem., 11, 481 (1965); J. T. Hougen, J. Mol. Spectrosc., 82, 92 (1980); J. T. Hougen, J. Mol. Spectrosc., 89, 296 (1981). [ Links ] [ Links ] [ Links ] [ Links ] [ Links ]
[7] A. J. Merer and J. K. G. Watson, J. Mol. Spectrosc., 47, 499 (1973). [ Links ]
[8] W. G. Harter and C. W. Patterson, J. Chem. Phys. 66, 4872, 4886 (1977). [ Links ]
[9] Y.N. Chiu, J. Molec. Struct. 108, 211 (1984). [ Links ]
[10] A.P. Sutton, Electronic Structure of Materials, Oxford Univ. Press, London, 1994 [ Links ]
[11] a) C.A. UtretasDíaz, J. R. Letelier, Computers and Chem., 19, 39 (1995); b) J. R. Letelier, A.ToroLabbé, C.A. UtrerasDíaz, Spectrochim. Acta, 47A, 29 (1991); (c) J.R. Letelier, C.A. UtrerasDíaz, Spectrochim. Acta, 53A, 247 (1997). [ Links ] [ Links ] [ Links ]
[12] J. R. Henderson, R. A. Willett, M. Muramoto and D. C. Richardson, Tables of Harmonic FranckCondon Overlap Integrals Including Displacement of Normal Coordinates, Douglas Rep. SM45807, Douglas Missle and Space System Division, 1964.Sec. London, Ser. A, 280 (1964) [ Links ]