Page 271 - Strategies and Applications in Quantum Chemistry From Molecular Astrophysics to Molecular Engineer
P. 271

254                                                         V. BARONE ET AL.
                            equation formally equivalent to that governing the motion of a particle with a unit mass in a
                            one dimensional space. This model is intimately linked to the use of local vibrational basis
                            functions centered at different points along the path. In our approach  [16,26,27], cubic
                            splines are used to interpolate the potential along the path and to generate a larger set of
                            equispaced points on which cubic splines are also used as basis functions. This kind of
                            treatment avoids any modelling of the ab-initio data and involves only analytical integrals.
                            Although the size of the basis set is larger than the one necessary when employing Hermite,
                            Morse or Gaussian functions, the spline approach remains competitive since the matrices to
                            be diagonalized are banded with a constant width of 7. Furthermore, no new integrals are
                            introduced by the computation of expectation values of observables (also represented by
                            spline  fittings), and  the  additional  computational  effort  depends on  the number  of
                            eigenstates to be taken into account, rather than on the dimension of the primitive spline
                            basis  set. The  expectation  value  of  a  given observable in the eigenstate j
                            corresponding to the eigenvalue  is given by



                            The temperature dependence of the observable is obtained by assuming a Boltzmann
                            population of the vibrational levels, so that









                            3. Results

                            Full geometry optimizations and calculations of harmonic force constants were performed at
                            the UMP2/6-311G** level. Although  this  is  not the main  concern of this  study, it  is
                            noteworthy  that the  relatively  unexpansive  theoretical  treatment we  have  developped
                            provides structural and spectroscopic parameters in close agreement with experiment (see
                            Table 1).  More precisely, the  harmonic  approximation  seems  quite adequate for
                            whereas strong anharmonicities affect the CH stretchings and the out-of-plane motion of
                                 The wave number of this latter vibration is increased to   (in much better
                            agreement with the experimental value of    by our one-dimensional anharmonic
                            treatment. Such a strong positive correction is in agreement with experimental estimates
                            [29]. From another point of view, the two radicals are well suited to point out the influence
                            of the shape of the potential well on vibrational effects: a simple well  and  a  double-
                            well with an high inversion  barrier
                            The  influence of out-of-plane bending on geometrical parameters, electronic energy and
                            coupling  constants is  shown in  Figures  2-4. The  linear relationship  between the s
                            coordinate and the  angle is well evidenced in Figure 2a. We recall that in our approach,
                            although the geometries in internal coordinates used to build the path are mass independent,
                            the arc length s varies with the atomic masses, whereas the reduced mass governing the
                            motion always remains unitary. The larger mass of fluorine versus hydrogen then explains
                            the lower  slope of the curve versus s for   than for   Also noteworthy is the
                            increase of the CH and CF bond lengths upon inversion (Figure 2b).
   266   267   268   269   270   271   272   273   274   275   276