Acessibilidade / Reportar erro

A theoretical study of the inversion and rotation barriers in Methyl-substituted Amines

Abstracts

Barrier heights of the internal rotation and inversion motions of methylamine, dimethylamine and trimethylamine molecules were calculated with the CCSD(T)//B3LYP methodology in combination with the cc-pVTZ, cc-pVQZ, and cc-pCVTZ basis sets of Dunning. The complete basis set (CBS) extrapolation procedure and core-valence (CV) correlation effects are also examined to the barrier heights. Our best estimate results (CCSD(T)/CBS+CV//B3LYP/cc-pVQZ) are in very good agreement with the experimental data, indicating the use of this methodology to provide accurate predictions of structures and barrier heights for other systems

energy barriers; methylamine; dimethylamine; trimethylamine; CCSD//B3LYP


Barreiras de energia para os movimentos de inversão e rotação das moléculas metilamina, dimetilamina e trimetilamina foram calculadas usando a metodologia CCSD(T)//B3LYP com os conjuntos de funções de base de Dunning, cc-pVTZ, cc-pVQZ e cc-pCVTZ. O procedimento de extrapolação para o conjunto de base completa (complete basis set, CBS) e os efeitos da correlação dos elétrons do caroço também foram incluídos nos cálculos das barreiras de energia. Nossos melhores resultados estão em excelente concordâncias com os dados experimentais, indicando que a metodologia utilizada é recomendada para a predição exata de propriedades estruturais e de barreiras de energia de outros sistemas moleculares


ARTICLE

A theoretical study of the inversion and rotation barriers in Methyl-substituted Amines

Júnior NascimentoI; Marina PelegriniIII; Luiz F. A. FerrãoI,II; Orlando Roberto-NetoIV; Francisco B. C. MachadoI,* * e-mail: fmachado@ita.br

IDepartamento de Química, Instituto Tecnológico de Aeronáutica, 12228-900 São José dos Campos-SP, Brazil

IIDepartamento de Física, Instituto Tecnológico de Aeronáutica, 12228-900 São José dos Campos-SP, Brazil

IIIDivisão de Ensino, Academia da Força Aérea, 13643-000 Pirassununga-SP, Brazil

IVDivisão de Aerotermodinâmica e Hipersônica, Instituto de Estudos Avançados, 12228-840 São José dos Campos-SP, Brazil

ABSTRACT

Barrier heights of the internal rotation and inversion motions of methylamine, dimethylamine and trimethylamine molecules were calculated with the CCSD(T)//B3LYP methodology in combination with the cc-pVTZ, cc-pVQZ, and cc-pCVTZ basis sets of Dunning. The complete basis set (CBS) extrapolation procedure and core-valence (CV) correlation effects are also examined to the barrier heights. Our best estimate results (CCSD(T)/CBS+CV//B3LYP/cc-pVQZ) are in very good agreement with the experimental data, indicating the use of this methodology to provide accurate predictions of structures and barrier heights for other systems.

Keywords: energy barriers, methylamine, dimethylamine, trimethylamine, CCSD//B3LYP

RESUMO

Barreiras de energia para os movimentos de inversão e rotação das moléculas metilamina, dimetilamina e trimetilamina foram calculadas usando a metodologia CCSD(T)//B3LYP com os conjuntos de funções de base de Dunning, cc-pVTZ, cc-pVQZ e cc-pCVTZ. O procedimento de extrapolação para o conjunto de base completa (complete basis set, CBS) e os efeitos da correlação dos elétrons do caroço também foram incluídos nos cálculos das barreiras de energia. Nossos melhores resultados estão em excelente concordâncias com os dados experimentais, indicando que a metodologia utilizada é recomendada para a predição exata de propriedades estruturais e de barreiras de energia de outros sistemas moleculares.

Introduction

The methyl-substituted amines, methylamine (CH3NH2), dimethylamine [(CH3)2NH], and trimethylamine [(CH3)3N] form an interesting group of molecules which are important in organic syntheses and biological processes, as well as they are efficient corrosion inhibitors of aluminum. Beyond their technological importance, these species are rich in structural features presenting large amplitude internal motions, rotation of the methyl groups and internal inversion of the amino group.

Methylamine has been extensively studied by experimentalists in which its internal rotation and inversion barrier heights were determined using microwave (MW), infrared (IR) and electron diffraction (ED) experiments.1-7 Methylamine is considered a small molecule from the theoretical point of view, therefore it is appropriate to be studied using high level theoretical quantum chemical methodologies, having received many contributions.7-14 Lee et al.8 investigated the origin of the structural stability of methylamine concluding that the stereo electronic effect is the major interaction force affecting its stability. Csonka and Sztraka9 using density functional and post-Hartree-Fock methodologies showed that the density functional BP86/6-311G(d) method provides reliable results for geometry and vibrational frequencies, as compared with the experimental results. Smeyers and Villa10 discussed the influence of the zero-point energy correction to the rotation and inversion motions based on the MP2 results, and During and Zheng7 calculated these motions barriers using B3LYP and MP2 methods. Recently, we have examined the convergence of the geometries, harmonic vibrational frequencies, and the barrier heights of the internal rotation and inversion motions of methylamine using the CCSD(T) theory with a systematic improvement via cc-pVDZ, cc-pVTZ, cc-pVQZ, and cc-pCVTZ basis sets and also including the contribution of the inner-shell and valence-shell electrons.11 More recently, Levi et al.12 carried out accurate calculations of the vibrational frequencies, Kim and Zeroka13 presented new and important results concerning the inversion and rotational barriers and Ducati et al.14 calculated the rotational barrier.

Geometries, vibrational frequencies and the inversion and rotation barriers of dimethylamine have been determined both experimentally and theoretically.15-23 The experimental investigations were obtained using electron diffraction,15 microwave,16,17 and infrared and Raman18,19 spectra. Theoretically, the calculations have been carried out using ab initio methods up to MP4 level.20-23

The trimethylamine molecule has also received various experimental and theoretical studies.24-31 The experimental data of the geometry, vibrational frequencies and internal-inversion and rotational barriers were obtained using microwave,24,26 infrared and Raman,25,27 techniques and employing analysis of the ultraviolet absorption and fluorescence spectra.28 Theoretical results were also obtained using ab initio methods up to MP4 level.29-31

It has been adequately demonstrated that a simultaneous consideration of basis set saturation and core correlation effects using coupled cluster method and extended basis set provide very accurate values of many molecular properties.32,33 Recently, using this methodology, we have investigated the geometries, harmonic vibrational frequencies, and the barrier heights of methylamine.11 In this work, we are concerned with the calculations of the internal and rotational barrier heights of the methyl-substituted amines, methylamine (CH3NH2), dimethylamine [(CH3)2NH], and trimethylamine [(CH3)3N].

Methodology

For all conformers, geometry optimizations were performed using the B3LYP exchange-correlation functional at the density functional theory level34,35 and the cc-pVQZ basis set of Dunning.36 Intrinsic reaction coordinate (IRC)37 calculations were carried out in order to verify the connections of the transition states to the equilibrium geometries, as already done for methylamine by Kim and Zeroka13 and for trimethylamine by Tanaka and Aida.31 Next, using these optimized geometries, the barrier heights of the internal rotation and inversion motion of the methyl-substituted amines were calculated using coupled cluster theory, including a perturbative treatment of the triple excitations (CCSD(T)).38 The correlation-consistent polarized-valence triple and quadruple basis sets of Dunning,36 named as cc-pVTZ, and cc-pVQZ, were employed. As it is well known, these basis set are built for frozen core electron calculations. The CCSD(T) total energies were extrapolated to the complete basis set (CBS) limit by using the extrapolation procedure of Halkier et al.39

where n equals 4 for cc-pVQZ basis set. In order to provide the effect of core-valence correlation (CV) we have also used the correlation consistent core and valence polarization triple zeta basis set (cc-pCVTZ) of Woon and Dunning.40 Thus, these calculations were carried out correlating both the valence electrons and all electrons. Therefore, the barrier heights (ΔE) were calculated using the expression,

As stated above, for methylamine, we have recently used CCSD(T) to optimize the geometries for all conformers and to calculate inversion and rotational barrier heights.11 Kim and Zeroka have shown recently two transition state structures for the inversion mode, with one and two imaginary frequencies, respectively.13 In order to clarify the previous structures of the inversion motion,11 we have also re-optimized the structures of all methylamine conformers using the CCSD(T) methodology.

All the calculations were carried out using the Gaussian 03 code41 and Molpro 2002.6.42

Results and Discussion

The B3LYP/cc-pVQZ optimized geometries of all methylamine conformers are presented in Figure 1, and in Figure 2 are showed the dimethylamine and trimethylamine conformers.



Methylamine has its ground state equilibrium geometry represented by a staggered conformation, and the eclipsed conformation corresponds to the transition state of the internal rotation due to the torsional motion of a methyl moiety along the CN axis. The wagging motions of the amine hydrogen atoms present two planar conformations, which are first-order and second-order inversion transition states, as also determined by Kim and Zeroka.13 All of them belong to the Cs symmetry and are shown in Figure 1. Recently, we have optimized the geometry of methylamine using the CCSD(T)/CBS methodology11 and, in this work, as discussed before, we have re-optimized the structures of all conformers using CCSD(T) methodology, in order to clarify the previous structures of the inversion motion. As one can see, the present results show that, in fact, the previous structures11 refer to the first-order inversion motion with one imaginary frequency. The imaginary frequencies calculated at the CCSD(T)/cc-pVTZ level are 749i cm-1 for the first-order inversion structure, 166i and 740i cm-1 for the second-order, and 317i cm-1for the rotation transition state. These imaginary frequencies calculated at B3LYP/cc-pVQZ level, are respectively 673i, 112i and 660i, and 297i cm-1 as reported in Figure 1. Note that the B3LYP imaginary frequencies differ no more than 80 cm-1 from the CCSD(T)/cc-pVTZ results. All geometric parameters calculated using B3LYP/cc-pVQZ are, in general, in good agreement with the CCSD(T) results and also with previous theoretical calculations7-10,13,14 and experimental data.2-4 We have also carried out the IRC calculations using the B3LYP/cc-pVTZ method and the results confirm that the first-order inversion transition state and the rotation transition state connect to the ground state, as already shown by Kim and Zeroka13 using a more complete two-dimensional potential energy surface. In this work, the IRC paths for the rotation motion were calculated by changing the dihedral angle (d) (H-C-N lone-pair). The changes in the energy relative to the ground state structure (d = ± 60º) are presented in Figure 3. The IRC calculations for the inversion motion are presented in Figure 4 that shows the potential energy as a function of the degree of inversion d, which is defined as:



where, α is the H-N-C-H dihedral angle, which is ± 58.8º at the ground state conformation (αgs) and equal to zero at the first-order inversion transition state geometry.

The origin of the structural stability was confirmed to be essentially due to the stereo electronic effect. As obtained by Lee et al.,8 the change of rotational conformation energy (ΔE) and the change of the C-N distance (ΔrCN) varies almost linearly with either the cosine of the dihedral angle d (H-C-N-lone-pair) or the cosine of 3d.

The two CCSD(T) results of the methylamine barrier heights, the single point calculations with the geometries optimized at the B3LYP/cc-pVQZ level and those obtained with the geometries optimized with the CCSD(T), are presented in Table 1. The results obtained with B3LYP/cc-pVQZ and the experimental data1-3,5 are also presented. All the CCSD(T) results for barrier heights using the B3LYP/cc-pVQZ optimized geometries are similar to the results calculated with the CCSD(T) geometries, differing no more than 0.08 kcal mol-1. It is important to mention that calculations of T1 diagnostic were previously carried out11 at the CCSD/cc-pVQZ level with values equal to 0.008 for all conformations, which are small enough to allow the use of the single reference CCSD(T) wave functions in the study of the present molecular conformations. As already noticed,11 the increase of the basis set and the inclusion of core correlation reduce both barrier height values, as also observed if one compares the CCSD(T)/cc-pVTZ result for the first-order rotation barrier with the one calculated by Ducati et al.14 with CCSD(T)/aug-cc-pVTZ. The inclusion of ZPE correction decreases the values of barrier heights underestimating the CBS-DT (CBS extrapolation using double and triple basis set) results compared to the experimental data. The experimental rotation barriers range from 1.962,3 to 2.055 kcal mol-1, and for inversion, from 4.831 to 5.565 kcal mol-1. The CCSD(T) results with the CBS extrapolation and core-valence correlation corrections (CV) (equation 2), but without ZPE correction, are 1.83 and 4.99 kcal mol-1 for the rotation and inversion barriers. However, the CCSD(T)/cc-pVQZ (valence) results, without ZPE correction are 1.95 and 5.61 kcal mol-1, respectively. As one can see, the inclusion of the ZPE correction from the B3LYP/cc-pVQZ or from the CCSD(T)/CBS-DT calculations provide barrier height values lower than the experimental data, which are accurate enough to encourage us in using this methodology in the calculation of larger systems as dimethyl and trimethyl amines. Also, as affirmed above, the CCSD(T) and CCSD(T)//B3LYP barrier height results differ no more than 0.08 kcal mol-1 without ZPE. Including the B3LYP/cc-pVQZ ZPE correction in the single point CCSD(T) calculations, the difference from the optimized CCSD(T) result is not higher than 0.14 kcal mol-1. Therefore, the CCSD(T) single point calculations are reliable enough to calculate barrier heights for similar structures, and in this work, were carried out for the larger systems.

The B3LYP/cc-pVQZ dimethylamine conformation structures are presented in Figure 2. The ground state equilibrium geometry has a pyramidal structure with Cs point group symmetry, and it will be denoted as the eclipsed conformation. The internal rotation of the methyl groups with C1 point group symmetry has a staggered conformation, and the inversion motion of the amine hydrogen has a planar conformation within C2v point group symmetry. As obtained to methylamine, all geometrical parameters are in good agreement with previous experimental15,16 and theoretical results.21-23 In the planar structure, the H-C-N-C-H atoms are coplanar to the amine hydrogen atom, and in the ground state pyramidal structure, this amine hydrogen atom forms an angle of 53.4 degrees to the H-C-N-C-H plane. These results are in good agreement with the MP2/6-31G(d,p) value, equal to 55.2 degrees, and also with the experimental data of 54.6 degrees.16 The IRC paths using the B3LYP/cc-pVTZ method for the rotation motion were calculated by changing the dihedral angle (d) (H-C-N-H). As calculated for methylamine, the changes in energy relative to the ground state structure (d = ± 60º) are presented in Figure 3. For the inversion the degree of inversion d is defined as shown in Figure 5 and the IRC path is presented in Figure 4.


As analyzed for methylamine, the change of rotational conformation energy (ΔE) and the change of the C-N distance (ΔrCN) varies almost linearly with either the cosine of the dihedral angle d (H-C-N-lone-pair) or the cosine of 3d. Thus, we also concluded that the origin of the structural stability is due to the stereoelectronic effect.

The dimethylamine barrier heights are calculated at the CCSD(T) level using the geometries optimized with B3LYP/cc-pVQZ and are presented in Table 2 together with the experimental data.16-19 As already discussed for the methylamine, the increase of the basis set and the inclusion of the core correlation reduce both barrier height values. The inclusion of core correlation is more important for the inversion barrier. For the internal rotation barrier, our best result obtained at the CCSD(T)/CBS+CV//B3LYP/cc-pVQZ level (equation 2), is 3.11 kcal mol-1, which include the CBS and core-valence correlation effects. The previous theoretical result15 at MP2/6-31G(df,p) level is 3.43 kcal mol-1and our B3LYP/cc-pVQZ result is 2.99 kcal mol-1. The experimental infrared and Raman data are respectively 3.01318 and 3.02919 kcal mol-1. Using ZPE correction, the B3LYP/cc-pVQZ result is 2.67, which is relatively small compared to the experimental data. Using B3LYP ZPE corrections in the CCSD(T)/equation 2 and CCSD(T)/cc-pVQZ results, the rotation barrier are 2.79 and 2.83 kcal mol-1, respectively. For the inversion, our best results using equation 2 and cc-pVQZ basis set are 5.11 and 5.62 kcal mol-1, respectively. The B3LYP/cc-pVQZ is 4.23 kcal mol-1 and previous theoretical results using MP2/6-31G (d,p) and MP4/6-31G(d,p)//MP2/6-31G (d,p) are respectively 6.60 and 6.29 kcal mol-1.22,23 The experimental result obtained from microwave spectra is 4.4 ± 1 kcal mol-1,16,17 which has an error limit somewhat large. A better comparison with the experimental data should include the ZPE correction obtained from the B3LYP/cc-pVQZ calculations, which are 3.24, 4.12 and 4.63 kcal mol-1 calculated using B3LYP/cc-pVQZ, CCSD(T)/equation 2 and CCSD(T)/cc-pVQZ, respectively.

Trimethylamine has its ground state equilibrium geometry well determined24,26,29-31 as a pyramidal structure with C3v point group symmetry with an eclipsed conformation. Figure 2 gives the C3v equilibrium geometry and the two other conformations with Cs point group symmetry, the staggered and planar structures, which correspond respectively to the internal rotation of the methyl groups and to the skeletal inversion. All geometries were optimized using the B3LYP/cc-pVQZ method. The geometrical parameters, in general, are in good agreement with previous experimental24,26,27 and theoretical results.29,30 The IRC paths using the B3LYP/cc-pVTZ method for the rotation motion were calculated by changing the dihedral angle d (H-C-N-lone pair) and the relative energy from the ground state structure (d = ± 60º), which are presented in Figure 3. For the inversion, the IRC path is presented in Figure 4. The calculations were carried out using the degree of inversion d, as defined by Tanaka and Aida.31 As analyzed for methylamine and dimethylamine, the structural stability is essentially due to the stereo electronic effect.

Calculated trimethylamine barrier heights are presented in Table 3 together with the experimental data.24,25,27,28 The inclusion of core correlation is more important for the inversion barrier, as already discussed for the methylamine and dimethylamine. The calculated barrier height results for the torsion of one methyl group are consistent with the experimental and previous theoretical results. The B3LYP/cc-pVQZ result without ZPE correction is 3.91 kcal mol-1 and the result with ZPE is 3.54 kcal mol-1, which seems underestimated compared to the experimental data, i.e., 4.40,24 4.0625 or 3.6327 kcal mol-1, as also observed to methylamine and dimethylamine. Our best result (equation 2) is 4.27 kcal mol-1, but if one uses the ZPE correction from the B3LYP/cc-pVQZ calculations, this value becomes 3.90 kcal mol-1. The previous MP4(SDQ)/aug-cc-pVTZ result without ZPE is 4.04 kcal mol-1.31 For the inversion, the B3LYP/cc-pVQZ without and with ZPE are 7.45 and 6.29 kcal mol-1, respectively. Our best result without ZPE is 9.51 kcal mol-1, which is similar to the previous theoretical study at MP4(SDQ)/aug-cc-pVTZ level, equal to 9.42 kcal mol-1.31 Using the ZPE correction from the B3LYP/cc-pVQZ calculations, our best result is 8.35 kcal mol-1, which is in very close agreement with the experimental value of 8.29 kcal mol-1.28

Conclusions

In this work, we have used the cc-pVTZ, cc-pVQZ, and cc-pCVTZ basis sets of Dunning and the CCSD(T) theory in order to obtain accurate values of inversion and rotation barrier heights of the methylamine, dimethylamine and trimethylamine molecules using the optimized geometries and frequencies calculated at B3LYP/cc-pVQZ level. CBS extrapolation procedure and core-valence correlation effects are also employed in the calculations of the barrier heights using CCSD(T) single point calculations. Our best estimate results (equation 2) are in very good agreement with the experimental data. Therefore, the methodology, CCSD(T)/CBS+CV//B3LYP/cc-pVQZ, used in this study provides accurate and consistent predictions of structures and barrier height values of internal motions in methyl substituted amines, and it is also recommended to be applied to other systems.

Supplementary Information

Supplementary information is available free of charge at http://jbcs.sbq.org.br as a PDF file.

Acknowledgments

The authors acknowledge the continuous research and the fellowship support of the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP).

Submitted: July 13, 2010

Published online: February 10, 2011

FAPESP has sponsored the publication of this article.

Supplementary Information

Table S1 - Click to enlarge

Table S2 - Click to enlarge

Table S3 - Click to enlarge

  • 1. Tsuboi, M.; Hirakawa, A. Y.; Tamagake, K.; J. Mol. Spectrosc. 1967, 22, 272.
  • 2. Takagi, K.; Kojima, T.; J. Phys. Soc. Jpn. 1971, 30, 1145.
  • 3. Belorgeot, C.; Stern, V.; Goff, N.; Kachmarsky, J.; Möller, K. D.; J. Mol. Spectrosc. 1982, 92, 91.
  • 4. Ijima, T.; Bull. Chem. Soc. Jpn. 1986, 59, 853.
  • 5. Kreglewski, M.; J. Mol. Spectrosc. 1989, 133, 10.
  • 6. Oda, M.; Ohashi, N.; Hougen, J. T.; J. Mol. Spectrosc. 1990, 142, 57.
  • 7. Durig, J. R.; Zheng, C.; Struct. Chem. 2001, 12, 137.
  • 8. Lee, S. J.; Mhin, B.J.; Cho, S. J.; Lee, J. Y.; Kim, K. S.; J. Phys. Chem. 1994, 98, 1129.
  • 9. Csonka, G. I.; Sztraka, L.; Chem. Phys. Lett. 1995, 233, 611.
  • 10. Smeyers, Y. G.; Villa, M.; Chem. Phys. Lett. 2000, 324, 273.
  • 11. Pelegrini, M.; Roberto-Neto, O.; Machado, F. B. C.; Chem. Phys. Lett. 2005, 414, 495.
  • 12. Levi, C.; Martin, J. M. L.; Bar, I.; J. Comput. Chem. 2008, 29, 1268.
  • 13. Kim, H.-W.; Zeroka, D.; Int. J. Quantum Chem. 2008, 108, 974.
  • 14. Ducati, L. C.; Custodio, R.; Rittner, R.; Int. J. Quantum Chem. 2010, 110, 2006.
  • 15. Beagley, B.; Hewitt, T. G.; Trans. Faraday Soc. 1968, 64, 2561.
  • 16. Wollrab, J. E.; Laurie, V. W.; J. Chem. Phys. 1968, 48, 5058.
  • 17. Wollrab, J. E.; Laurie, V. W.; J. Chem. Phys. 1971, 54, 532.
  • 18. Durig, J. R.; Griffin, M. G.; Groner, P.; J. Phys. Chem. 1977, 81, 554.
  • 19. Consalvo, D.; van Bladel, J. W. I.; Engeln, R.; Reuss, J.; Chem. Phys. 1993, 171, 221.
  • 20. Skaarup, S.; Griffin, L. L.; Boggs, J. E.; J. Am. Chem. Soc. 1976, 98, 3140.
  • 21. Senent, M. L.; Smeyers, Y. G.; J. Chem. Phys. 1996, 105, 2789.
  • 22. Senent, M. L.; Smeyers, Y. G.; Moule, D. C.; Mol. Phys. 1998, 94, 949.
  • 23. Senent, M. L.; Smeyers, Y. G.; J. Mol. Struct. (Theochem) 1999, 464, 137.
  • 24. Lide Jr., D. R.; Mann, D. E.; J. Chem. Phys. 1958, 28, 572.
  • 25. Möller, K. D.; DeMeo, A. R.; Smith, D. R.; London, L. H.; J. Chem. Phys. 1967, 47, 2609.
  • 26. Wollrab, J. E.; Laurie, V. W.; J. Chem. Phys. 1969, 51, 1580.
  • 27. Durig, J. R.; Craven, S. M.; Bragin, J.; J. Chem. Phys. 1970, 53, 38.
  • 28. Halpern, A. M.; Ondrechen, M. J.; Ziegler, L. D.; J. Am. Chem. Soc. 1986, 108, 3907.
  • 29. Kölmel, C.; Ochsenfeld, C.; Ahlrichs, R.; Theor. Chim. Acta 1992, 82, 271.
  • 30. Marsh, C. M. B.; Schaefer III, H. F.; J. Phys. Chem. 1995, 99, 195.
  • 31. Tanaka, M.; Aida, M.; Chem. Phys. Lett. 2006, 417, 316.
  • 32. Martin, J. M. L.; Taylor, P. R.; Chem. Phys. Lett. 1995, 242, 343.
  • 33. Martin, J. M. L.; Taylor, P. R.; Chem. Phys. Lett. 1996, 248, 336.
  • 34. Becke, A. D.; J. Chem. Phys. 1993, 98, 5648.
  • 35. Lee, C.; Yang, W.; Parr, R. G.; Phys. Rev. B 1988, 37, 785.
  • 36. Dunning Jr., T. H.; J. Chem. Phys. 1989, 90, 1007.
  • 37. González, C.; Schlegel, H. B.; J. Phys. Chem. 1990, 94, 5523.
  • 38. Raghavachari, K.; Trucks, G. W.; Pople, J. A.; Head-Gordon, M.; Chem. Phys. Lett. 1989, 157, 479.
  • 39. Halkier, A.; Helgaker, T.; Jørgensen, P.; Klopper, W.; Koch, H.; Olsen, J.; Wilson, A. K.; Chem. Phys. Lett. 1998, 286, 243.
  • 40. Woon, D. E.; Dunning Jr., T. H.; J. Chem. Phys. 1995, 103, 4572.
  • 41. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A.; Jr.; Vreven, T; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. P.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; González, C.; Pople, J. A.; Gaussian 03, Revision D.01, Gaussian, Inc.: Wallingford CT, 2004.
  • 42. Werner, H. -J.; Knowles, P. J.; Almlöf, J.; Amos, R. D.; Berning, A.; Cooper, D. L.; Deegan, M. J. O.; Dobbyn, A. J.; Eckert, F.; Elbert, S. T.; Hampel, C.; Lindh, R.; Lloyd, A. W.; Meyer, W.; Nicklass, A.; Peterson, K.; Pitzer, R.; Stone, A. J.; Taylor, P. R.; Mura, M. E.; Pulay, P.; Schütz, M.; Stoll, H.; Thorsteinsson, T.; MOLPRO, version 2002.6 a package of ab initio programs. 2002.
  • *
    e-mail:
  • Publication Dates

    • Publication in this collection
      16 May 2011
    • Date of issue
      May 2011

    History

    • Received
      03 July 2010
    • Accepted
      10 Feb 2011
    Sociedade Brasileira de Química Instituto de Química - UNICAMP, Caixa Postal 6154, 13083-970 Campinas SP - Brazil, Tel./FAX.: +55 19 3521-3151 - São Paulo - SP - Brazil
    E-mail: office@jbcs.sbq.org.br