The high resolution spectrum of AsH3 (4 0 0) local mode state: Symmetry reduction and rotational re-quantization
Ji-Xin Cheng,Xiao-Gang Wang,Hai Lin,Qing-Shi Zhu
DOI: https://doi.org/10.1016/S1386-1425(98)00157-7
1998-01-01
Abstract:The (4 0 0) stretching overtone band of AsH 3 in the range of 7900–8100 cm −1 was recorded by a Fourier transform spectrometer with a resolution of 0.01 cm −1 . The A 1 or A 2 −E degeneracy of the upper state rotational levels up to J′ =8 was observed. It is interesting to find that this spectrum resembles an A , C hybrid band of an asymmetric top and can be assigned by only one set of quantum numbers J , K a and K c . The rotational levels were determined by combination difference method and fitted in both the symmetric top and the asymmetric top models. In the latter case, the effective rotational constants A , B and C are in good agreement with the values predicted by a simple local mode model, indicating that the dynamical symmetry of AsH 3 in the (4 0 0) local mode state has reduced to that of an asymmetric top. The transition from normal mode to local mode vibration from v =1 to 4 in AsH 3 was discussed. Keywords AsH 3 Symmetry reduction Rotational re-quantization 1 Introduction The rotational energy level structure of the stretching overtone states of small symmetric molecules is of great interest due to its unusual features near the local mode limit. First, it has been demonstrated that there exist simple relations (so called ‘ α -relations’) between the effective vibration–rotation parameters in the normal mode approach [1–3] . The other peculiarity is the reduction of dynamical symmetry [2, 4–7] , that is, T d → C 3 v C 3 v → C s and C 2 v → C s for XY n , ( n =4, 3, 2) molecules respectively. In detail, the spectrum of a T d molecule in the local mode state should resemble a parallel band of a symmetric top molecule and can be assigned by a set of rotational quantum numbers ( J , K ) of symmetric tops. This interesting phenomenon has been observed in GeH 4 ( n 0 0 0), n =3 [8] and 5 [9] , SiH 4 ( n 0 0 0), n =3, 4 and 5 [10, 11] and SnH 4 ( n 0 0 0), n =7 [12] and 6, 8 [13] states. For C 3 v molecules, the rotational analysis on the v =2 [14] , 3 [15] and 4 [16] stretching overtones of SbH 3 provided good examples of the local mode features in XY 3 molecules, i.e. the disappearance of the Coriolis couplings and the simple arithmetic relations between some vibration–rotation parameters. The recent FT study [17] of GeH 3 D ( v =3) and ICLAS study [18] of GeH 3 D ( v =6, 7, 8) together with SbH 3 ( v =4) [16] showed the rotational energy level pattern is close to that of an asymmetric rotor. AsH 3 is a typical pyramid molecule. One expects to be able to observe the local mode behavior in its stretching overtone states. The overtone spectra with medium resolution have been recorded and analyzed with a local mode model [19] . Later Fermi resonance between the stretching modes and the bending modes was investigated [20] . The first [21] and second [22] stretching overtone bands have been rotationally analyzed. In the (3 0 0) state, the A 1(2) –E (A 1(2) denotes A 1 or A 2 ) 3-fold degeneracy of the rotational levels has been observed and an analysis in the asymmetric top model demonstrated the tendency of symmetry reduction from C 3 v to C s [22] . In the present contribution we present a detailed rotational analysis of its third overtone state. It is found that the AsH 3 (4 0 0) state is a good example of symmetry reduction and rotational re-quantization of symmetric top molecules at the local mode limit. 2 Experimental AsH 3 with purity 99.99% was purchased from Nanjing Special Gas Company. A high resolution spectrum was recorded at room temperature on the Bruker 120 HR interferometer, which was equipped with a CaF 2 beam splitter, a InSb detector and a tungsten halogen source. A sample pressure of 33 torr was chosen. No optical filter was applied. A multi-pass absorption cell with effective path length of 105 m was used. Scans (2250) were co-added. The apodized resolution 0.01 cm −1 (MOPD) is comparable with the Doppler width 0.011 cm −1 . The precision of the line position was estimated to be better than 0.001 cm −1 . The spectral range from 5900 to 8900 cm −1 was covered. The spectrum from 7900 to 8100 cm −1 was assigned to the (4 0 0) stretching overtone band according to the local mode model calculation (see below). The other two weaker bands in the range of 6120–6500 cm −1 and 8130–8350 cm −1 were assigned to (2 1 0) and (3 1 0) transitions respectively. The calibration was carried out by measuring the residual H 2 O lines in the range of 7300–7600 cm −1 [23] . The accuracy of the measurements is better than 0.001 cm −1 [24] . 3 Assignment and analysis 3.1 Vibrational analysis Since the third stretching overtone vibration is close to the local mode limit, the local mode theory can be utilized to give a reasonable estimation of the vibration–rotation parameters of the upper states. As has been demonstrated, local mode vibration can be formulated in two equivalent models [25, 26] , harmonically coupled enharmonic oscillator model [27] and normal mode model with Darling–Dennison resonance included (NMDD), while the latter is more commonly used by spectroscopists so far for rotational analysis. In the present case, the three parameters required by the NMDD model were obtained by fitting the known stretching band centers, which are listed in Table 1 . The symmetrized local mode states |400; A 1 〉 and ∣400; E〉 can be written as the linear combination of the normal mode states in the v =4 vibrational manifold, i.e. (1) 400 ; A 1 =0.534 1 ; 3 +3 + 1 ; 3 −3 2 +0.441 0 ; 4 0 +0.687 2 ; 2 0 +0.219 4 ; 0 0 (2) 400 ; E ± =0.180 0 ; 4 ±4 +0.654 1 ; 3 ±1 +0.417 3 ; 1 ±1 +0.360 0 ; 4 ∓2 +0.486 2 ; 2 ∓2 where v 1 ; v 3 l 3 denotes the normal mode basis. The fundamentals v 1 and v 3 have been analyzed in several papers [28–31] . In our previous work [22] , we refitted the fundamental parameters with the published energy levels in Ref. [29] up to J ′=10. The effective vibration–rotation parameters for the ∣400; A 1 ) and ∣400; E〉 states can thus be calculated [32, 33] by using the eigenvectors shown in Eqs. (1) and (2) . The calculated values are listed in Table 2 , which were used to give rise to the initial prediction of the upper state rotational levels. 3.2 Description and assignment of the spectrum The (4 0 0) band in the range of 7900–8100 cm −1 as shown in Fig. 1 corresponds to the transitions to the degenerate local mode pair, i.e. ∣400; A 1 〉 and ∣400; E〉 states. The rotational constants of these two states are nearly equal. On the other hand, as a quasi-spherical molecule, the difference between the constants B and C is rather small so that the rotational levels with the same J ′ are close in energy. As a result the rovibrational levels are strongly coupled by the H 22 terms which are essential for the rotational pattern of local mode state [1–4, 7] . In the spectrum, more transitions become visible due to the mixing of the vibration–rotation wavefunctions. Fig. 2 shows an example of the doubling of some transitions in the R (1) and R (2) branches. The ground state constants were determined by the far infrared spectrum [34] . A refinement has been carried out recently [35] . The ground state (GS) constants used in the present work are taken from Ref. [34] , which are found to be adequately accurate for our study. The calculated effective parameters of the upper states were used to predict the rotational energy levels. Besides, the band origins were estimated from the recorded spectrum. Moreover, the quartic centrifugal distortion constants D J , D K and D JK were initially constrained to the ground state values. All the other high-order parameters used in the calculation were set to zero initially. The selection rules applied for both A 1 and E states were Δ J =0, ±1, Δ K =0, ±1, ±2, ±3, E↔E and A 1 ↔A 2 , which produce a combination of parallel and perpendicular transitions. The ground state combination differences (GSCD) linking levels to a common upper rotational state were used for the assignment of the low J ′ rotational levels. Two interesting features were found and used to verify the assignment. First, as indicated by the prediction, the rotational levels of symmetry A 1(2) and E with the same J ′ and K values are degenerate into one (for K =0) or two (for K ≠0) pairs with splitting <0.01 cm −1 . Second, due to the strong mixing of the vibration–rotation wavefunctions, four to five transitions can be found for many rotational levels of the upper states. The assigned rotational levels were then used in the fit to obtain improved values of the spectroscopic parameters and to predict the levels for higher J ′. The same process was repeated until J ′=8. Altogether 317 transitions were assigned. It should be noted that some high J ′ levels are accidentally perturbed and the A 1(2) –E degeneracy is consequently destroyed. The perturbed levels are marked by * in Table 3 and not used in the fit. For J ′≥9 levels, some lines were observed and assigned but not used in the fit, because most assigned rotational levels were found to deviate from the prediction, moreover, the A 1(2) –E degeneracy can hardly be satisfied. This indicates the existence of perturbation. A preliminary analysis indicates the probable perturber is the (3 0 0; v 2 = v 4 =1) state. Its band center was estimated to be ∼8040.42 cm −1 by an addition of the band centers of the (3 0 0) state [22] and the v 2 = v 4 =1 state [29] . 3.3 Rotational analysis The (4 0 0) A 1 /E states were treated as a pair of fundamentals [14–18, 21, 22] , therefore all the vibration–rotation parameters were effective ones. The rotational levels were calculated by diagonalizing the vibration-rotation matrix. The diagonal element is (3) νl ; Jk H νl ; Jk =ν 0 +B ν J J+1 −k 2 +C ν k 2 −2 Cζ z ν kl−D J v J 2 J+1 2 −D JK v J J+1 k 2 −D K v k 4 +… where the vibrational angular momentum quantum number is l =0, ±1 for the 400 ; A 1 and 400 ; E ± , respectively. The off-diagonal matrix elements include the Coriolis interactions and H 22 terms. The xy -type Coriolis coupling between 400 ; A 1 and 400 ; E takes the following form [36] (4) 400 E + ; Jk+1 H/hc 400 A 1 ; Jk =− 400 A 1 ; Jk+1 H/hc 400 E − ; Jk = 2 B Ω ζ y eff F J, k The H 22 terms include l -, r - and α -resonance operators [37, 32] . The l - and r -resonances connect the rotational levels of different vibrational angular momentum in the E state. The l -resonance matrix element is (5) 400 E + ; Jk+2 H/hc 400 E − Jk =− 1 2 q eff +q eff J J J+1 F J, k F J, k+1 The r -resonance matrix element is (6) 400 E − ; Jk+1 H/ hc 400 E + ; Jk =2r eff 2k+1 F J, k The α -resonances couple the rotational levels between the A 1 and E state with the following matrix elements (7) 400 A 1 ; Jk+2 H/hc 400 E + ; Jk = 400 E − ; Jk+2 H/hc 400 A 1 ; Jk = 1 2 2 α eff xx F J, k F J, k+1 and (8) 400 E + ; Jk+1+ H/hc 400 A 1 ; Jk = 400 A 1 ; Jk+1 H/hc 400 E − ; Jk = 1 2 2 α eff zx F J, k where F ( J , k ) quoted in Eqs. (4)-(8) is J J+1 −k k+1 . All the off-diagonal matrix elements assume the phase convention of Di Lauro and Mills [36] . A least squares fit of the upper state rotational levels up to J ′=8 converged with a standard deviation σ of 1.09×10 −2 cm −1 . Altogether 19 parameters were included in the fit. The effective Coriolis coupling coefficients Cζ z and B Ω ζ y were fixed to the values calculated by the NMDD model because they are too small to be determined. We found r eff has a strong correlation with other H 22 parameters, so it was also fixed to the calculated value. Since the K =±3 splitting can be observed in the spectrum, we tried to include the ε parameter in the fit, however it could not be determined significantly. This indicates that the K =±3 splitting is caused mainly by the H 22 resonances. The parameters of the upper states are listed in Table 2 . It can be seen that the observed parameters are in good agreement with the calculation, as in the case of SbH 3 [14–16] and AsH 3 [21, 22] . The observed and calculated rotational energy levels are listed in the left side of Table 3 . The A 1(2) –E degeneracy is obvious. Besides, it is interesting to note the formation of 6-fold degeneracy of K =0 and 1 levels with the increase of J ′ as shown in Table 3 . This phenomenon can be compared with the formation of 4-fold clusters in H 2 Se [38] . lt should be noted that the quantum numbers v (referred as ∣400; A 1 〉 or ∣400; E〉) and K in Table 3 are only assigned according to the calculation. In fact we found that several pairs of levels exchanged their positions in the fit process, therefore we constrained the weights of these levels to zero and included them only in the last step of the fit. This indicates that K is no longer a ‘good’ quantum number to assign the rotational levels and that a new set of quantum numbers are required. 4 Symmetry reduction and rotational re-quantization Assuming symmetry reduction in the local mode limit, Lehmann established the Hamiltonians in the local mode basis for XY n ( n =2, 3, 4) molecules. Then he obtained the vibration-rotation coupling terms by transforming these Hamiltonians to a symmetrized basis [4] . Alternatively, Child and Zhu derived the effective rotational Hamiltonians of the local mode states by averaging the vibration–rotation Hamiltonians over the corresponding local mode vibrational wavefunctions [6, 7] . For an XY 3 molecule, the effective rotational Hamiltonian at the local mode limit can be transferred into that of an asymmetric top with three different effective rotational constants [39, 7] (9) A= 1 2 B ν +C ν +E + B ν −C ν +E 2 +4G 2 1 2 (10) B=B ν −E (11) C= 1 2 B ν +C ν +E − B ν −C ν +E 2 +4G 2 1 2 where B ν and C ν are the effective rotational constants in the symmetric top Hamiltonian. G and E are related to the α coefficients by the following relations: (12) E= 2 3 α 13 xx − 1 3 q 3 n r ̄ 2 n − 0 r ̄ 2 0 (13) G= 4 3 r 3 + 2 3 α 13 zx n r ̄ 2 n − 0 r ̄ 2 0 On the other hand, the effective off-diagonal H 22 parameters q eff , r eff , α eff zx and α eff xx can be directly related to their corresponding fundamental coefficients [32, 3] . Replacing the symmetrized local mode wavefunctions by linear combination of single bond local-mode wavefunctions and transforming the normal coordinates into the bond displacement coordinates, the integration over the vibrational coordinates yields the following relations [3] : (14) E=−q eff = 2 2 α eff xx (15) G=4r eff = 1 2 α eff zx . lt can be easily proved that Lehmann’s result can be transformed into the same forms as Eqs. (9)-(11) under the following relations (16) q=2E, r=2G where q and r are the rotational coupling coefficients in Lehmann’s model [4] . According to the prediction of Ovchinnikova [39] and Lehmann [4] , the overtone spectrum should resemble an A , C hybrid band of an asymmetric top at the limit of small inter-bond couplings. In order to assign the AsH 3 (4 0 0) band by a set of rotational quantum numbers of asymmetric top molecules, we consider its ground state as an asymmetric top at the oblate symmetric top limit and assign its rotational levels with the quantum numbers J , K a and K c , where K c is equal to K for the symmetric top. In this way, each pair of ground rotational levels with the same J and K c values have identical energies except for the A 1 –A 2 splitting in the case of K= 3 n ( n is an integer). Following the selection rule of an A , C hybrid band, i.e. Δ J =0, ±1 and Δ K a =0, 2±, Δ K c =±1, ±3 or Δ K a =±1, ±3, Δ K c =0, ±2, the rotational levels of the (4 0 0) local mode vibrational state were assigned as those of an asymmetric top. More than five transitions were found for many upper states, as shown in Table 3 . The energy difference of one common upper state obtained from different transitions is <0.01 cm −1 for most levels. Here the quantum number K c corresponds to the quantum number K in the normal mode picture [4] . Fig. 3 shows some details of the assignment. The perturbation causes the splitting of a few rotational levels, corresponding to the destruction of A 1(2) –E degeneracy in the normal mode approach. In these cases the energies were calculated by averaging the split sub-levels. The perturbed levels were given smaller weights in the fit. The obtained rotational levels were fitted as an asymmetric rotor by the following Hamiltonian of A reduction [40] (17) H rot A =B x J x 2 +B y J y 2 +B z J z 2 − Δ JK J 2 J z 2 − Δ K J z 4 − Δ J J 4 − 1 2 δ J J 2 +δ K J z 2 , J + 2 +J − 2 + +… where [ A , B ] + =AB + BA . The initial values of B α ( α = x , y , z ) were constrained to that of the (3 0 0) state [22] while all the quartic and higher centrifugal distortion constants were set to zero initially. The observed and calculated energy levels are listed in the right hand side of Table 3 . Altogether nine spectroscopic parameters were determined by a least-squares fit. The results are shown in Table 4 . The standard deviation of 1.10×10 −2 cm −1 is satisfactory. Although AsH 3 is a quasi-spherical top molecule, the symmetry reduction at the local mode limit does not assume the formation of an accidentally symmetric top. In fact the asymmetry parameter x =2 B −( A − C )/( A − C ) is nearly zero in the present case. Therefore we did not adopt the S reduction. The values of E and G in Eqs. (14) and (15) were obtained by averaging the vibrationally off-diagonal parameters in Table 2 . E= 1 2 2 2 α eff xx −q eff , G= 1 2 4r eff + 1 2 α eff zx . The effective rotational constants A , B and C calculated from Eqs. (9)-(11) are listed in Table 4 . The good agreement between the calculation and observation proves that AsH 3 in its (4 0 0) state has reduced to an asymmetric rotor. Although both normal mode and local mode methods can be applied to assign and analyze AsH 3 (4 0 0) band, however, the rotational states within one J′ block of the local mode pair are so strongly mixed that quantum number K loses its meaning. For example, the eigenfunction of the rotational level of J′=3 and E symmetry with energy 8071.6328 cm −1 can be written as a combination of zero-order rovibrational states ∣ v ; J′K Γ〉, where v denotes A 1 or E states and Γ denotes the rotational symmetry A 1 , A 2 or E.′ (18) Ψ J′=3 ; E =0.45 A 1 ; 31E +0.23 A 1 ; 32E +0.64 E; 30E +0.45 E; 31E +0.23 E; 32E −0.19 E, l +1 ; 33E +0.19 E, l −1 ; 33E On the other hand, the satisfactory assignment as an A , C hybrid band of an asymmetric top demonstrates the rotational levels of a symmetric top should be re-quantized by a set of quantum numbers J , K a , and K c instead of J and K at the local mode limit. The correspondence between the degeneracy and symmetry of the rotational levels in the symmetric and asymmetric top representations is shown in Table 5 . 5 Local mode effect on the spectroscopic parameters It is known that the following simple relations ( α- relations) exist between the effective vibration–rotation parameters of XY 3 molecules when the local mode limit is approached [1–3] . −q eff =α eff BB / 2 , 4r eff =α eff BC / 2 (19) B ν A 1 =B ν E , C ν A 1 =C ν E , Cζ z eff = B Ω ζ y eff =0 From Table 2 it can be seen that the above relations are well satisfied. For example, the observed vibrationally off-diagonal parameters satisfy the ratios q eff a eff xx =0.989 :1.414, r eff : α eff zx =0.990 :5.656, which are very close to the theoretical prediction. Moreover, the Coriolis couplings are nearly zero, in agreement with the theory. It should be noted the α -relations can lead directly to the A 1(2) –E degeneracy. According to the vibration–rotation Hamiltonian defined in Eqs. (5)-(8) , the A and E matrix blocks of J′ =2 levels can be written as A block: 2 ; −1 4r eff 3 α eff xx 0 0 1 ; 1 3 2 α eff zx −3q eff 0 0 ; 0 3 2 α eff zx 3 α eff xx symmetric −1 ; −1 4r eff −2 ; 1 and E block: 2 ; 1 2 2 α eff zx − 6 q eff 0 0 1 ; 0 3 2 α eff zx − 3 2 α eff xx 0 0 ; −1 2 6 r eff 3 α eff xx symmetric −1 ; 1 2 2 α eff zx 2 ; 0 Here the energy levels are identified by ∣k; l 〉. These two blocks are equivalent when the rotational constants of A 1 and E states are equal and the α -relations are satisfied, consequently, the A 1(2) –E degeneracy can be naturally observed. As has been done in the analysis of the spectra of SbH 3 [16] and GeH 3 D [17] , the vibration–rotation Hamiltonian in terms of internal coordinates can be formulated in a simple form by defining a new set of rovibrational parameters [2] as (20) ν + = v s +2v t /3 ν − =ν s −ν t B + = B s +2B t /3 B − =B s −B t C + = C s +2C t /3 C − =C s −C t D J+ = D J s +2D J t /3 D J− =D J s −D J t D JK+ = D JK s +2D JK t /3 D JK− =D JK s − D JK t D K+ = D K s +2D K t /3 D K− =D K s −D K t q + =q eff −2 1/2 α eff xx q − =q eff +2 −1/2 α eff xx r + =−4r eff −2 1/2 α eff zx r − =4r eff −2 1/2 α eff zx where the subscripts s and t stand for A 1 and E states. The new coefficients with plus subscript are diagonal in the local mode basis, while those with minus subscript describe the rotational coupling between different local mode states. We transformed the observed parameters of AsH 3 (2 0 0), (3 0 0) and (4 0 0) bands into the above forms. The results are shown in Table 6 . lt can be seen that those parameters with minus subscript tend to become zero with the increase of the vibrational quantum number from v =2 to 4, however, D J , D JK and D K are exceptions. The probable reason is that, in both (3 0 0) and (4 0 0) state, the r eff parameter was fixed to the theoretical calculation, whose error should be compensated by adjusting the D values in the fit. The above results demonstrate the disappearance of rotational couplings between different local mode states ( n 0 0), (0 n 0) and (0 0 n ) with the increase of vibrational energy, therefore the rotational level structure should resemble that of an asymmetric rotor. This is equivalent to the analysis in Section 4 . 6 Discussion When the vibrational energy is localized in one single bond, the principal axes of the effective rotational moment inertia are rotated by an angle θ in the plane defined by the C 3 axis and the excited bond [4] . As a result, the C 3 axis is no longer the symmetric axis in the local mode vibration. The projection of the angular momentum along the C 3 axis cannot be regarded as a commutator of the molecular Hamiltonian any more. The rotation should be re-quantized in this case because the asymmetric top quantum numbers K a , K c offer a clearer picture to the rotational structure. Using the local mode relation r =( A − C ) sin θ derived by Lehmann [4] and Eq. (16) , the rotated angle of the effective principal axes is estimated to be 16.5° for the (4 0 0) stretching vibration of AsH 3 , which can be compared with θ =13.9° in the (3 0 0) state of AsH 3 [22] . lt has been shown that the high-energy stretching vibrations tend to be bond-localized owing to the anharmonic couplings [27, 41] . The rotational structure of the AsH 3 stretching overtones exemplifies this transition from normal mode to local mode vibration. Fig. 4 shows the J′ =3 levels of the ( n 0 0) states ( n =1, 2, 3 and 4). The (1 0 0) state corresponds to the fundamental normal mode vibrational states v 1 =1 (A 1 ) and v 3 =1 (E), which are separated by 11.26 cm −1 [29] in band center. The level structure shows obvious symmetric top characters and K is a ‘good’ quantum number. With the increase of vibrational quantum number, the lowest pair of A 1 /E normal mode states (local mode pair) become more and more degenerate owing to the Darling–Dennison resonance. Consequently the rotational structure becomes complicated due to the strong rovibrational r -, l - and α -resonances as demonstrated in the first stretching overtone. However, it is interesting that the rotational structure of the higher overtones exhibits simple patterns as shown in the (3 0 0) and (4 0 0) vibrational state. In these cases, the local mode pair become nearly degenerate, i.e. the separation of the band centers is <0.01 cm −1 . The strong resonances due to the vibrationally off-diagonal H 22 operators rearrange the rotational energy levels and lead to A 1(2) –E degeneracy as shown in Fig. 4 . Acknowledgements This work is supported by the National Natural Science Foundation of China and Chinese Academy of Sciences. The authors thank Georges Graner for a careful reading of the manuscript. References 1 L. Halonen A.G. Robiette J. Chem. Phys. 84 1986 6861 2 T. Lukka L. Halonen J. Chem. Phys. 101 1994 8380 3 X.G. Wang Q.S. Zhu J. Chem. Phys. 105 1996 8011 4 K.K. Lehmann J. Chem. Phys. 95 1991 2361 5 K.K. Lehmann J. Chem. Phys. 96 1992 7402 6 M.S. Child Q.S. Zhu Chem. Phys. Lett. 184 1991 41 7 Q.S. Zhu Spectrochim. Acta 48 1992 193 8 Q.S. Zhu B.A. Thrush J. Chem. Phys. 92 1990 2691 9 F.G. Sun X.G. Wang J.L. Liao Q.S. Zhu J. Mol. Spectrosc. 184 1997 12 10 Q.S. Zhu B.S. Zhang Y.R. Ma H.B. Qian Sppectrochim. Acta 46 1990 1217 11 Q.S. Zhu B.S. Zhang Y.R. Ma H.B. Qian Spectrochim. Acta 46 1990 1323 12 M. Halonen X. Zhan J. Chem. Phys. 101 1994 950 13 X Zhan M. Halonen L. Halonen H. Bürger O. Polanz J. Chem. Phys. 102 1995 3911 14 M. Halonen L. Halonen H. Bürger P. Moritz J. Chem. Phys. 95 1991 7099 15 M. Halonen L. Halonen H. Bürger P. Moritz Chem. Phys. Lett. 203 1993 157 16 J. Lummila T. Lukka L. Halonen H. Bürger O. Polanz J. Chem. Phys. 104 1996 488 17 H Bürger L. Halonen Mol. Phys. 87 1996 227 18 H Bürger A. Campargue L. Halonen H. Ruland J. Mol. Spectrosc. 183 1997 99 19 M. Halonen L. Halonen H. Bürger P. Moritz J. Phys. Chem. 96 1992 4225 20 P. Lukka E. Kauppi L. Halonen J. Chem. Phys. 102 1995 5200 21 S.F. Yang X.G. Wang Q.S. Zhu Spectrochim. Acta Part A 53 1997 157 22 O.N. Ulenikov F.G. Sun X.G. Wang Q.S. Zhu J. Chem. Phys. 105 1996 7310 23 R.A. Toth Appl. Opt. 33 1994 4851 24 G. Guelachvili Pure Appl. Chem. 68 1996 193 25 K.K. Lehmann J. Chem. Phys. 79 1983 1098 26 I.M. Mills A.G. Robiene Mol. Phys. 56 1985 743 27 M.S. Child L. Halonen Adv. Chem. Phys. 57 1984 1 28 W.B. Olson A.G. Maki R.L. Sams J. Mol. Spectrosc. 55 1975 52 29 O.N. Ulenikov A.E. Cheglokov G.A. Shevchenko M. Winnewiser B. Winnewiser J. Mol. Spectrosc. 157 1993 141 30 V Dana J.-Y. Mandin G. Tarrago W.B. Olson B. Bezard J. Mol. Spectrosc. 159 1993 468 31 J.-Y. Mandin V. Dana G. Tarrago S. Klee B.P. Winnewisser J. Mol. Spectrosc. 172 1995 319 32 L. Halonen J. Chem. Phys. 86 1987 588 33 X.G. Wang, Ph.D Thesis, University of Science and Technology of China, Hefei, China, 1997. 34 M. Carlotti G. Di. Lonardo L. Fusina J. Mol. Spectrosc. 102 1983 310 35 G. Tarrago V. Dana J.-Y. Mandin S. Klee B.P. Winnewisser J. Mol. Spectrosc. 178 1996 10 36 C. Di. Lauro I.M. Mills J. Mol. Spectrosc. 21 1966 386 37 F.W. Birss Mol. Phys. 31 1976 491 38 J.-M. Flaud C. Camy-Peyret H. Bürger P. Jensen I.N. Kozin J. Mol. Spectrosc. 172 1995 126 39 M.Y. Ovchinnikova Chem. Phys. 120 1988 249 40 J.K.G. Watson J. Durig Vibration Spectra and Structure 6 1977 Elsevier Amsterdam 41 Q.S. Zhu H.F. Li X.G. Wang Chem. Phys. Lett. 212 1993 403