## Journal of the Brazilian Chemical Society

##
*Print version* ISSN 0103-5053*On-line version* ISSN 1678-4790

### J. Braz. Chem. Soc. vol.13 no.3 São Paulo June 2002

#### http://dx.doi.org/10.1590/S0103-50532002000300006

__Article__

**The Raman Spectrum of the Squarate (C _{4}O_{4}^{-2} ) Anion: An Ab Initio Basis Set Dependence Study**

*Sandro G. de Miranda and Pedro A. M. Vazquez* *

*Instituto de Química, Universidade Estadual de Campinas, CP 6154, 13083-970 Campinas - SP*

O perfil de excitação Raman do anion esquarato, C

_{4}O_{4}^{-2}, foi calculado por métodosab initiono nível de teoria Hartree-Fock utilizando a Teoria de Resposta Linear (TRL) para seis frequências de excitação: 632,5, 514,5, 488,0, 457,9, 363,8 e 337,1 nm. Um total de cinco conjuntos de funções de base (6-31G*, 6-31+G*, cc-pVDZ, aug-cc-pVDZ e bases de Sadlej) foi investigado visando comparar o desempenho da base 6-31G*, em termos de convergência numérica dos resultados e custo computacional, com os dos demais conjuntos de base. Todos os conjuntos de base testados reproduzem as principais características do espectro Raman do anion esquarato no intervalo de excitação estudado. A base 6-31G* apresentou, em média, o mesmo desempenho nos resultados numéricos dos conjuntos de base maiores mas a uma fração do custo computacional mostrando-se apropriada para o estudo teórico dos derivados e complexos deste anion.

The Raman excitation profile of the squarate anion, C

_{4}O_{4}^{-2}, was calculated usingab initiomethods at the Hartree-Fock using Linear Response Theory (LRT) for six excitation frequencies: 632.5, 514.5, 488.0, 457.9, 363.8 and 337.1 nm. Five basis set functions (6-31G*, 6-31+G*, cc-pVDZ, aug-cc-pVDZ and Sadlej's polarizability basis set) were investigated aiming to evaluate the performance of the 6-31G* set for numerical convergence and computational cost in relation to the larger basis sets. All basis sets reproduce the main spectroscopic features of the Raman spectrum of this anion for the excitation interval investigated. The 6-31G* basis set presented, on average, the same accuracy of numerical results as the larger sets but at a fraction of the computational cost showing that it is suitable for the theoretical investigation of the squarate dianion and its complexes and derivatives.

Keywords: Raman spectra, squarate anion, vibrational spectroscopy,ab initioLinear Response Theory

Introduction

The advances in computer technology and theoretical methods over the last decades has made possible the *ab initio* calculation of molecular properties such as equilibrium geometries, infrared absorption intensities, vibrational frequencies and force fields for a large number of molecules and chemical systems of spectroscopic interest.^{1} These calculated properties proved to be a valuable aid in the interpretation and assignment of the experimental spectroscopic data for a large number of molecules. However, in many cases, the peculiar electronic structure of the system being studied presents Raman spectroscopic features for which no routine computational methods are widely available. One of these cases is the squarate dianion (C_{4}O_{4}^{-2} ), which belongs to the monocyclic oxocarbons (C_{n}O_{n}^{-2}, 3 £ n £ 6 ) class of organic compounds with a wide range of applications. Since the first vibrational spectroscopic study by West and Ito,^{2} the electronic, infrared and Raman spectra of the squarate anion as well as several of its complexes and derivatives have been the subject of several experimental and theoretical studies over the last decades (see, for example, the reviews by Santos^{3} and Schleyer^{4}). The Raman resonance and pre-resonance spectra of aqueous solutions of C_{4}O_{4}^{-2} were studied by Ito *et al.*^{5,6 }in the 632-337 nm spectral range. The experimental excitation profiles obtained by these authors were interpreted theoretically in terms of the Jahn-Teller effect associated to the p,p^{*} electronic transition with the lowest doubly degenerated excited electronic state, and this effect was assumed to be the responsible for the intensity exhibited by the non-totally symmetric modes. As far as we know, no detailed *ab initio* calculations of the excitation profiles for this compound have been performed to date. While the calculation of static, or frequency-independent, Raman intensities has been made routine by almost all current *ab initio* packages, this is not the case for the computation of dynamic, or frequencydependent, Raman intensities. For resonant or near resonance Raman spectra, depending on the nature of the problem and following Hellers' methodology,^{7,8} a series of techniques have been reported using single reference correlated calculations for the ground state mixed with *ab initio* or semi-empirical excited state calculations followed by numerical procedures to estimate the gradient of the excited electronic state energy surface with respect to each of the ground state's normal modes.^{9-12} As pointed out by Hudson,^{10} one of the caveats of this method, is that it is necessary that there be only one electronic transition in the region of study in order to use the excited state calculation. A systematic approach for the computation of dynamic Raman intensities, nonetheless, has been recently proposed by Helgaker and coworkers.^{13} Using the Linear Response Theory (LRT)^{1,14,15 }to compute the frequency-dependent polarizabilities, these authors developed a Hartree-Fock and MCSCF formalism, and a computational implementation where the calculation of Raman intensities, within the Placzek polarizability theory,^{16} involves the differentiation of the frequency-dependent polarizabilities, computed at the laser excitation frequency with respect to nuclear displacements. Although, basis set requirements for the computation of static and dynamic polarizabilities have been the subject of recent studies, only few reports have dealt with Raman intensities calculations. Fleicher and Pulay, using B3LYP/6-31G* level of theory, obtained semiquantitative agreement for the gas phase Raman intensities and depolarization ratios of benzene and coronene.^{17} In a study carried out by Halls and Schlegel,^{18} it was found that the Sadlej^{19,20} polarizability basis set yields the best static intensities quality with respect to the computational cost. These results agree with those previously reported by Van Caillie and Amos,^{21} and by Dulskov and Sauer^{22 }for the calculation of static and dynamic polarizabilities of a range of small molecules using correlated and non correlated methods. On the other hand, earlier vibrational dynamic Raman optical activity (VROA) calculations conducted by Helgaker *et al*.^{13} argued in favor of the aug-cc-pVDZ basis set in order to achieve numerical stability on the calculated Raman intensities and, at least, on the use of the aug-cc-pVTZ basis set to achieve Hartree-Fock limit quality on the vibrational frequencies. However, it has been pointed out that, these bases sets are far too large to be used in the calculation of dynamic Raman intensities for most of the systems of chemical and spectroscopic interest. In the methodology developed by Helgaker, the basis set chosen may impose a severe impact on the time needed to solve the iterative frequency response equations. In this method, the computation of the gradient of the polarizability requires a numerical differentiation with respect to the N nuclear positions demanding 2*3*N computations which results in a very high computational cost. Together, these factors can make the calculation of dynamic Raman intensities by this method unfeasible for medium sized molecules. With respect to the C_{4}O_{4}^{-2} anion, a former *ab initio* study of the vibrational force field of the squarate anion conducted by Torii and Tasumi, using the 6-31G*, 6-31+G* and 6-31+G(df) basis sets, found that the basis set has little effect on the correct description of the geometric and vibrational parameters of this anion, and that these parameters were more sensitive to the selection of theoretical methods for taking electron correlation into account.^{23} Although the squarate anion alone may be studied using theory levels higher than Hartree-Fock and using large bases sets, this is not the case for the majority of its derivatives and complexes for which experimental spectroscopic data is available. In most of these cases, the 6-31G*/Hartree-Fock level of theory represents the best compromise in terms of the computational cost and the quality of the results. Therefore, a survey of the performance of the 6-31G* basis set in terms of computational costs and numerical convergence of the intensities is needed.

Our primary goal in this work is to compare the performance of the 6-31G* basis set to the larger ones, at the Hartree-Fock *ab initio* level of theory, in the calculation of the gas phase excitation profiles of the squarate anion in order to investigate factors such as the computational cost and the quality of the results in terms of its numerical convergence. The inclusion of electronic correlation effects, albeit important for accurate semiquantitative results, will not be considered initially as we are concerned with how well Hartree-Fock calculations can predict dynamic Raman intensities over a wide spectral range of exciting frequencies. These results will serve as a guide to gauge future calculations of the excitation profiles of squarate anion derivatives and complexes.

This paper is organized as follows; in the next section, we detail the computational and theoretical methods as well the bases sets used. Then, we present an analysis of the observed computational demand for each of the basis sets. Next, a summary of the previous studies and a brief theoretical background of the Linear Response Theory as it relates to the computation of Raman intensities are presented, followed by a breakdown of the numerical results. An analysis in terms of the numerical convergence of the excitation profiles and the ability of the bases sets studied to qualitatively reproduce the main spectroscopic features of the squarate anion is presented. We close our discussion by inspecting how the theoretical calculations in the gas-phase can be related to the experimental data in condensed phases.

Computational Details

All computations were carried out with the Dalton electronic structure program.^{24} For the sake of timing and computational resource comparisons, a cluster of FreeBSD/Athlon computers running at 700 MHz and equiped with 128Mbytes of RAM, 64 kbytes of cache level 1, 512 kbytes of level 2 cache, 200Mbytes of virtual memory and standard UDMA 66Mhz disk subsystem was used in all cases. The program was compiled with the GNU Fortran77 2.95.2 compiler and linked against the linear algebra optimized ATLAS library. The following criteria were applied for the basis set selection. To evaluate the importance of the diffuse functions in the calculation of dynamic Raman intensities, the Pople's 6-31+G* basis set was included in the study. Following Helgaker's VROA results,^{13} we also selected two Dunning's double zeta correlation consistent bases sets, cc-pVDZ and aug-cc-pVDZ. In addition, due to its high success rate on predicting static polarizabilities,^{18,21,22 }the Sadlej polarizability basis set was also selected. All the basis sets were obtained from the latest versions available at the Extensible Computational Chemistry Environment Basis Set Database.^{25} For each basis set, the molecular geometry was optimized without symmetry constraints and an analytical Hessian was computed at the stationary equilibrium point for later use in the Raman intensities computations. In order to investigate the behavior of each of the basis sets in this methodology at the Hartree-Fock level over a wide spectral range of excitation frequencies, the polarizabilities were computed at six experimental laser excitation frequencies, in the interval from 632.8 to 337.1 nm, the same as used by Ito and co-workers for the measurement of the Raman excitation profiles of this anion in aqueous solution.^{5,6} A numerical differentiation using symmetric nuclear displacements of ± 0.0001 a.u. about the equilibrium nuclear positions was done in order to obtain the gradients of the polarizabilities at each exciting frequency. The wave function energy convergence threshold was set to 10 ^{-8} a.u. and the iterative response equations convergence threshold was set to 10 ^{-7} a.u. The resulting absolute intensities were normalized, following Ito methodology, with respect to the intensity of the *v*_{2} (a_{1g}) band at 632.8 nm excitation which was taken as unity.

Results and Discussion

*Computational demand of the basis sets*

We commence our discussion by focusing our attention on the elapsed wall times for the calculation of the Raman intensities for each of the basis sets. As reported in Table 1, it is clear that this method may impose severe limitations on the basis set choices for the calculation of Raman intensities even for medium sized molecules like the squarate anion and at the Hartree-Fock level of theory. If we take into account the relative sizes of the bases sets versus the time required to compute the excitation profiles, it seems that there is no direct relation between them. The explanation for this apparent discrepancy resides in the numerical iterative methods employed by the current computational implementation. Most of the computing time is associated with the resolution of the iterative linear response equations until a given numerical convergence threshold is reached. In addition to the usual bottlenecks and overheads associated with the use of large basis sets such as the disk access, this step is highly dependent on the number of simultaneous excitation frequency responses requested, on the basis sets employed for a given theory level and on the convergence thresholds applied both to the wave function and to the response equations. If we take these results as the only criteria for basis set selection, the Sadlej's and the aug-cc-pVDZ basis sets should be discarded without further considerations due to their very slow convergence speed of the response equations. These results shouldn't be extrapolated, though. In another study of the Raman excitation profiles of nitrosyl chloride at MCSCF level of theory, we found that the Sadlej's basis set has the fastest convergence rate, surpassing even the smaller cc-pVDZ basis set^{26}. From the standing point of current experimental research being done with this molecule, however, these results show that survey computational studies at the 6-31G*/HF level of theory will provide results at a computational cost that is less than 50% that of the other basis sets considered. Given the reliability and success of this basis set in several studies involving geometry and frequency calculations we need to investigate, now, how it compares with the other sets with respect to the Raman intensities and excitation profiles.

*Summary of results from previous experimental studies*

The observed spectrum of the squarate anion in aqueous solution is dominated by three intense bands, *v*_{10} at 647 cm ^{-1}, *v*_{2} at 723 cm ^{-1} and *v*_{5} at 1123 cm ^{-1} . These bands were unambiguously assigned, respectively, to the b_{2g} ring deformation, to the a _{1g} totally symmetric CC stretching and to the b _{1g} asymmetric CC stretching.^{ 2,9 }The main features of the spectra reported by Ito *et al* can be summarized as follows^{5 }. Even at the longer excitation wavelengths, n_{10} and *v*_{5} have intensities comparable to the most intense band *v*_{2}. When the excitation wavelength is changed from red towards blue, these three bands show great intensity enhancement, with predominance of *v*_{5} as the strongest band above 488.0 nm and with *v*_{10} and *v*_{2} having approximately the same intensity in the near ultraviolet region. The strong enhancements observed for the three bands at the 337.1 and 363.8 nm excitation wavelengths were attributed to pre-resonance effects at the edge of the p ® p^{*} (*A*_{1g}®^{ 1}*E _{u}* ) electronic transition. The possibility of resonance due to the forbidden, very weak,

*n*® p

^{*}(

*A*

_{1g}®

^{ 1}E

_{g}) absorption occurring in the same spectral interval was discarded based on symmetry considerations that predict, also, the symmetry species b

_{1g}and b

_{2g}as the effective species in the Jahn-Teller effect. These results show resemblance with those reported by Santos

*et al*for the solid sodium salt,

^{27}but, in Santos'

*et al*work, even if no qualitative or semiquantitative numerical estimates of the relative intensities for these three bands in the solid state sodium were reported, their results show that in the Raman spectra of this salt the totally symmetric

*v*band shows a weaker dependence on the exciting frequency in relation to the observed enhancements of the

_{2}*v*and

_{10}*v*bands, being the weakest of the three bands for exciting lines in the 514.5 to 406.7 nm interval. An analysis of the salt spectra shows another feature that attracts our attention and contrasts with the spectroscopic behavior of this anion in different phases of aggregation; the asymmetric

_{5}*v*

_{9}(b

_{2g}) CO stretching, that shows very little intensity enhancement in aqueous solution, presents, in the solid salt spectrum with, response to the exciting frequency which is comparable to that of

*v*

_{5}.

*Dependence of the calculated excitation profiles on the basis set *

In this section, the computed intensities are analyzed in terms of basis set dependence of the numerical convergence of the results. How these results compare with the experimental data is discussed in the next section. The real part of the Linear Response Theory expression for the evaluation of the a_{a,b}*-th* element of the frequency-dependent polarizability tensor for an exact state^{14} is given by:

where w* _{b}* is the exciting frequency of the incident homogeneous periodic electric field and w

*denotes the energy difference between the exact eigenstate |0ñ of H*

_{k}= E_{0}¾ E_{k}_{0}at t = -¥:

and those of the k -th residual exact eigenstates:

of the set {|0ñ, |*k*ñ} of exact eigenstates of the non interacting unperturbed system. Although for the Hartree-Fock method the ground state is the best variational approximation to the eigenstate |0ñ, the remaining excited states are not accurate representations of the {|*k*ñ} set of excited eigenstates of the system for the evaluation of the expected values of the dipole transition matrix elements á0*|*m* ^{i}*|

*k*ñ and á

*k|*m

*|0ñ between the reference state |0ñ and the excited state |*

^{i}*k*ñ . At the same time, different basis sets will produce different sets {

*E*} of eigenvalues affecting the denominators (w

_{k}_{b }± w

_{k}) in the summation in equation 1. In the particular cases where the basis set yields virtual states for which w

_{b}@ w

_{k}, the terms:

will have a very significant weight in the summation in equation 1 enhancing the frequency dependency of the molecular polarizability. These two factors characterize the basis set dependence of the calculation of the frequency-dependent polarizability tensor from the standpoint of Linear Response Theory. Raman bands, on the other hand, derive part of their intensities from the magnitude of the variation of the polarizability tensor with each normal mode of vibration, (¶a* **/ *¶*Q _{i}*) . Unless the {E

*} set of eigenvalues has a very strong dependence on the nuclear positions for very small nuclear displacements, like those used in this work, we should not expect significant changes in the (w*

_{k}_{b}± w

_{k}) denominators due to numerical differentiation. For a given exciting frequency we should expect, instead, that the behavior of the transition dipole elements with respect to the nuclear displacements, due to symmetry changes, will determine the basis set dependence of the intensities.

After these considerations, we turn our attention to the results. The calculated normalized gas phase intensities as well as the aqueous solution experimental results of the two measurements carried out by Ito are summarized in Table . In overall, and in accordance with the experimental evidence, all basis sets correctly predict an enhancement of the intensities as the exciting line wavelength varies from red to blue with a dramatic intensification effect for wavelengths below 457.9 nm. In the UV region, the spectra are dominated by the *v*_{5} (b_{1g} ) band at 1123 cm ^{-1} while the *v*_{10} (b_{2g} ) band at 723 cm ^{-1} shows weaker intensity ratios as observed experimentally in aqueous solution. The calculated intensity ratios exhibit larger values than the experimental ones by factors that vary from 1.3 up to approximately 6.0 in the UV region of excitation. In the 632-457.9 nm excitation spectral interval all basis sets produce similar numerical results, with the exception of Sadlej's basis. The calculated intensity ratios for the totally symmetric *v*_{2} (a_{1g} ) band shows an average deviation of 0.06 among all the basis sets, while for *v*_{10} (b_{2g} ) band this value reduces to 0.04 showing fast basis set convergence for these two vibrational modes. For the *v*_{5} (b_{1g} ) band, however, the basis sets with diffuse functions show an increase in the response to the exciting frequency. To inspect in more detail the observed behavior, the vertical energies of the first low lying E* _{u}* excited states and the corresponding oscillator strengths were computed at the equilibrium geometry using the Linear Response Theory within the Random Phase Approach (RPA) for each of the basis sets. These results are summarized in Table 3. An analysis of the energies displayed in this table shows that the basis sets without diffuse functions produce excited states 1

^{1}

*E*

_{u}and 2

^{1}

*E*

_{u}that are far from the energy of the higher exciting frequencies used in the calculations (see Table 2 ) reducing the contribution of these states to the summation in equation 1. For the basis sets with diffuse functions, on the other hand, this situation changes drastically; the three low lying excited states are very near to each other and nearer to the energies of the exciting frequencies than the states generated by the other basis sets and this results in higher sensitivities of these sets to the exciting frequencies used in the computation of the Raman intensities. This condition is more dramatic for Sadlej's basis sets where all three states have significant weights as compared to the remaining basis sets. Of the three normal modes investigated, the nuclear displacements associated with the

*v*

_{5}(b

_{1g}) vibration produce, for all the basis sets, the strongest variation of the transition dipole matrix elements resulting in stronger enhancements as compared to the remaining vibrational modes.

The results can be summarized as follows: all basis sets reproduce the main features of the observed spectra of the squarate anion in aqueous solution at different excitation wavelengths. The excitation profiles associated with the *v*_{2} (a_{1g}) and *v*_{10} (b_{2g} ) vibrational modes show very little dependence on the basis set if compared to the *v*_{5} (b_{1g} ) mode. RPA calculations of the vertical excitation energies at equilibrium geometry with different basis may precede Raman calculations to define their relationship with respect on the excitation frequencies to be employed in the study of the system of interest. The 6-31G* basis set yields intensity ratios of quality comparable to those returned by larger basis sets.

*Theoretical versus experimental Raman excitation profiles*

In the previous section we investigated the Raman excitation profiles of the three most intense bands of the squarate anion in gas-phase with respect their dependence to the basis set and, now, it is worth to query how these results relates to the observed profiles in aqueous solution. Despite the large amount of experimental data and the body of theoretical investigations of this anion and its derivatives there are very few structural studies of the squarate anion in aqueous solution. Raman bandshape analyses conducted by Ribeiro *et al* strongly suggest the existence of structurally well-defined water cages such as those formed via hydrogen bonds with the oxygens of the solute.^{28} Molecular dynamics and quantum *ab initio* studies carried out by Martins^{29,30} to investigate the structure of the first solvation shell found a well-defined structure consisting of approximately 18 water molecules around the anion with about 12 of these molecules tightly bonded to the oxygens forming a highly symmetric solute-solvent structure. The *ab initio* results pointed out in the direction of a solute whose electronic structure is not considerable affected in the presence of the solvent, making, in principle, conceivable a correlation between the theoretical calculations in the gas phase and the results in aqueous solution. Indeed, if we compare the trends of the theoretical excitation profiles for the squarate anion with those reported for aqueous solution by Ito and those measured in the solid state by Santos we find a better resemblance between the gas phase results with those obtained in water solutions. Figure 1 depicts the simulated Raman spectra using the experimental Raman intensity ratios from Ito^{6} and the computed intensity ratios at 6-31G*/HF level of theory. The curves were produced with the Spectrum^{31 }program using a half-width of 10 cm ^{-1} for the three bands. A visual inspection of Figure 1 shows that the theoretical results reproduce, approximately, the pattern of the observed spectra in aqueous solution. This comparison is valid inasmuch as we take into consideration a series of limitations associated with both the experimental and theoretical methods employed. According to Myers, Raman intensity measurements are accurate to no more than ± 10-15% even in the best cases.^{32} The difficulties associated with the measurement of the excitation profiles of the squarate anion in aqueous solution can be estimated by contrasting, in Table 2, the two measurements made by Ito *et al* in 1975 and 1978 where deviations of the intensity ratios of up to 60% are observed. At the same time, the theoretical methods employed in this study, as discussed, are very dependent on an appropriate description on the excited states. A better description can be achieved by the inclusion of static electronic correlation via MCSCF wave functions but this implies the use of appropriate correlation consistent basis sets and the accurate definition of the active space^{33} in order to take into account all states associated with Raman enhancement with a significant increase in the computational costs. Small active spaces for high symmetry compounds will result in relatively fast calculations but only those states included will be correctly described. The inclusion of dynamic correlation via a single configuration reference state such as MP2/SOPPA^{34} or coupled cluster LRT^{35} are out question for current computational resources for most of the derivatives of this molecule.

Conclusions

In this work, the gas phase *ab initio* pre-resonance profiles for the squarate anion computed at the Hartree-Fock level of theory are reported for the first time. The basis set dependence of the intensities is related both to the sensitivity of the transition dipole matrix elements to the normal modes and to the spectrum of vertical excitation energies produced by each of the basis sets. The overall features of the calculated spectra are in excellent agreement with the experimental data. Our results show that *ab initio* calculations conducted at the 6-31G*/HF level of theory produce numerical results comparable to those obtained with larger basis sets at a fraction of the computational cost in the spectral range of laser exciting lines investigated. The squarate anion has a peculiar electronic structure that requires the inclusion of electronic correlation effects in the description of the excited states that contribute to the Raman intensities. We have no doubt that further theoretical studies where these effects are taken into account in the computation should bear better results. Unfortunately, such kind of investigation is currently out of reach for our computational resources. We hope this work will stimulate further experimental and theoretical investigations of the Raman excitation profiles of the squarate anion, and its derivatives.

Acknowledgments

The authors thank FAEP/Unicamp (grants 0692/00 and 0901/01) and FAPESP (grant 1998/06871-1) for financial support and NIC/Br for computer time. PAMV dedicates this work to the memory of professor Euclydes Francisco da Rocha Fraga.

References

1. Gauss, J. In *Modern Methods and Algorithms of Quantum Chemistry*, Vol. 3, 2nd. ed.; Grotendorst, J., ed.; John von Neumann Institute for Computing: Jülich, 2000, p 541. [ Links ]

2. Ito, M.; West, R.; *J. Am. Chem. Soc*. **1963**, *85*, 2580. [ Links ]

3. de Oliveira, L. F. C.; Mutarelli, S. R.; Gonçalves, N. S.; Santos, P. S.; *Quim. Nova* **1992**, *15*, 55. [ Links ]

4. v. R. Schleyer, P.; Najafian, K.; Kiran, B.; Jiao, H.; *J. Org. Chem.* **2000**, *65*, 426. [ Links ]

5. Iijima, M.; Udagawa, Y.; Ito, M.; *Chem. Phys.* **1975**, 9, 229. [ Links ]

6. Takahashi, M.; Kaya, K.; Ito, M.; *Chem. Phys.* **1978**, *35*, 293. [ Links ]

7. Lee, S. Y.; Heller, E. J.; *J. Chem. Phys.* **1979**, *71*, 4777. [ Links ]

8. Heller, E. J.; Sundberg, R. L.; Tannor, D.; *J. Phys. Chem.* **1982**, *86*, 1822. [ Links ]

9. Torii, H.; Tasumi, M.; Bell, I. M.; Clark, R. J. H.; *Chem. Phys.* **1997**, *216*, 67. [ Links ]

10. Hudson, B. S.; Markham, L. M.; *J. Raman Spectrosc.* **1998**, *29*, 489. [ Links ]

11. Myers, A. B. *J. Raman Spectrosc.*; **1998**, *29*, 437. [ Links ]

12. Waterland, M. R.; Stockwell, D.; Kelley, A. M.; *J. Chem. Phys.* **2001**, *114:14*, 6249. [ Links ]

13. Helgaker, T.; Ruud, K.; Bak, K. L.; Jorgensen, P.; Olsen, J.; *Faraday Disc.* **1994**, *99*, 165. [ Links ]

14. Olsen, J.; Jorgensen, P.; *J. Chem. Phys.* **1985**, *82*, 3235.; [ Links ]Levine, I.; *Quantum Chemistry*, 4th. ed, Prentice-Hall: Englewood Cliffs; 1991. [ Links ]

15. Jorgensen, P.; Jensen, H. J. A.; *J. Chem. Phys.* **1988**, *89*, 3654. [ Links ]

16. Placzek, G. *The Rayleigh and Raman Scattering*; United States Atomic Energy Commision: Lawrence Radiation Laboratory, University of California, Livermore, California, UCRL Translation No. 526 (L), Physics, 1959. [ Links ]

17. Fleicher, U.; Pulay, P.; *J. Raman Spectrosc.* **1998**, *29*, 473. [ Links ]

18. Halls, M. D.; Schlegel, H. B.; *J. Chem. Phys.* **1999**, *111*, 8819. [ Links ]

19. Sadlej, A. J.; *Collect. Czech. Chem. Commun.* **1995**, *53*, 1988. [ Links ]

20. Sadlej, A.; *Theor. Chim. Acta.* **1991**, *79*, 123. [ Links ]

21. Caillie, C. V.; Amos, R. D.; *Chem. Phys. Letters* **1998**, *291*, 71. [ Links ]

22. Dalskov, E. K.; Sauer, S. P. A.; *J. Phys. Chem. A* **1998**, *102*, 5269. [ Links ]

23. Torii, H.; Tasumi, M.; *J. Mol. Struct. (Theochem)* **1995**, *334*, 15. [ Links ]

24. Helgaker, T.; Jensen, H. J. A.; Joergensen, P.; Olsen, J.; Ruud, K.; Aagren, H.; Bak, K. L.; Bakken, V.; Christiansen, O.; Dahle, P.; Dalskov, E.; Enevoldsen, T.; Fernandez, B.; Heiberg, H.; Hettema, H.; Jonsson, D.; Kirpekar, S.; Kobayashi, R.; Koch, H.; Mikkelsen, K.; Norman, P.; Packer, M. J.; Ruden, T. A.; Saue, T.; Sauer, S. P. A.; Sylvester-Hvid, K. O.; Taylor, P. R.; Vahtras, O. *Dalton release 1.1 (2000), an *electronic structure program University of Oslo, Norway. [ Links ]

25. Basis sets were obtained from the Extensible Computational Chemistry Environment Basis Set Database, Version 4/22/01, as developed and distributed by the Molecular Science Computing Facility, Environmental and Molecular Sciences Laboratory which is part of the Pacific Northwest Laboratory, P. O. Box 999, Richland, Washington 99352, USA, and funded by the U. S. Department of Energy. The Pacific Northwest Laboratory is a multi-program laboratory operated by Battelle Memorial Institute for the U. S. Department of Energy under contract DE-AC06-76RLO 1830. Contact David Feller or Karen Schuchardt for further information.

26. Vazquez, P. A. M.; *Abstracts of the XI Simpósio Brasileiro de Química Teórica*, Caxambu, Brasil, 2001. [ Links ]

27. Santos, P. S.; Sala, O.; Noda, L. K.; Goncalves, N. S.; *Spectrochim. Acta* **2000**, *56A*, 1553. [ Links ]

28. Ribeiro, M. C. C.; de Oliveira, L. F. C.; Santos, P. S.; *Chem. Phys.* **1997**, *217*, 71. [ Links ]

29. Martins, L. R. *M. Sci. Thesis*; Universidade Estadual de Campinas: Campinas, Brazil, 2000. http://biq.iqm.unicamp.br/arquivos/teses/ficha43137.htm, accessed in March 2002. [ Links ]

30. Martins, L. R.; Vazquez, P. A. M.; Skaf, M. S.; *J. Mol. Struct. (Theochem)* **2001**, *580, *137*.* [ Links ]

31. Soper, P. D.; E. I. *du Pont de Nemours and Company* 1997 .

32. Myers, A. B.; *J. Raman Spectrosc.* **1997**, *28*, 389. [ Links ]

33. Knowles, P.; Schütz, M.; Werner, H.-J. In *Modern Methods and Algorithms of Quantum Chemistry*, Vol. 1, 1st. ed.; Grotendorst, J., ed.; John von Neumann Institute for Computing: Jülich, 2000, p 97. [ Links ]

34. Packer, M. J.; Dalskov, E. K.; Enevoldsen, T.; Jensen, H. J. A.; Oddershed, J.; *J. Chem. Phys.* **1996**, *105*, 5886. [ Links ]

35. Christiansen, O.; Halkier, A.; Koch, H.; Jorgensen, P.; *J. Chem. Phys.* **1998**, *108*, 2801. [ Links ]

*Received: August 13, 2001 Published on the web: April 5, 2002 *

**FAPESP helped in meeting the publication costs of this article.**

* e-mail: vazquez@iqm.unicamp.br