Anales de la Asociación Química Argentina
versión impresa ISSN 0365-0375
An. Asoc. Quím. Argent. v.94 n.1-3 Buenos Aires ene./jul. 2006
Electron delocalizaton interactions and NMR spin-spin coupling constants in saturated cage compounds
1Contreras, R. H.; 2Díez, E., 3Esteban, A.L.; 4Della, E.W.; 4Lochert, I. J.
1 Departamento de Física, FCEyN, Universidad de Buenos Aires, Ciudad Universitaria, P. 1
and CONICET, (C1428EHA) Buenos Aires, Argentina. FAX: +54-11-4576-3357, E-Mail: email@example.com
2 Departamento de Química Física Aplicada, Facultad de Ciencias, C2, Universidad Autónoma de Madrid, E-28049 Madrid, Spain.
3 Departamento de Química Física, Universidad de Alicante, Apartado 99, E-03080 Alicante, Spain
4 Department of Chemistry, School of Physical Sciences, Flinders University, Bedford Park, South Australia 5042, Australia.
Received March 1st, 2006. In final form May 10th, 2006
Dedicated to Prof. Imre G. Csizmadia on the occasion of his 75th birthday
Since the cage substrate of 1-X;3-Y-bicyclo[1.1.1]pentanes is notably strained, it is expected that in compounds of this type important trans-cage spin-spin coupling constants, SSCCs, could be observed. Such SSCCs should be dominated to a great extent by the Fermi Contact, FC, contribution since this term, for long-range couplings, is transmitted by electron-transfer interactions, which in strained compounds are particularly enhanced. In this paper, both theoretically and experimentally, 4J(13Ca,1H3) SSCCs in 1-X-bicyclo[1.1.1]pentanes, where the a-atom of group X is 13C, are studied. Results are compared with similar trans-cage SSCCs in analogous compounds. Insights into electron delocalization interactions defining the experimental trends are obtained through the study of s-hyperconjugative interactions of cage bonds. These interactions are studied using the Natural Bond Orbital, NBO, method. Results suggest that multipath additivity should hold in compounds where trans-cage SSCCs are determined by charge transfer interactions.
Dado que la estructura cerrada del sustrato de los pentanos 1-X;3-Y-biciclo[1.1.1]se encuentra tensionada, es de esperar que en compuestos de este tipo existan acoplamientos spin-spin "trans-cage", SSCCs. Tales SSCCs deberían ser dominadas por la contribución del Contacto Fermi, FC, ya que este término, para acoplamientos de largo alcance, se transmite por interacciones de transferencia electrónica, las cuales se potencian en compuestos tensionados. En este trabajo se estudian, teórica y prácticamente, 4J(13Ca,1H3) SSCCs en 1-X-biciclo[1.1.1]pentanos donde el átomo s del grupo X es 13C. Los resultados son comparados con SSCCs "trans-cage"en compuestos análogos. Se obtiene información sobre las interacciones de delocalización electrónica que definen las tendencias experimentales a partir del estudio de interacciones s-hiperconjugadas de los enlaces en la estructura. Estas interacciones son estudiadas usando el método NBO. Los resultados sugieren que la aditividad multimecanística debería ser válida en compuestos donde los "trans-cage" SSCCs son determinados por interacciones de transferencia de carga.
Indirect spin-spin coupling constants, SSCCs, originate in magnetic electron-nucleus interactions . The non-relativistic theory of these spectral parameters was presented by Ramsey more than a half century ago [2,3]. According to such a theory, isotropic SSCCs, commonly known as "scalar couplings", are made up of four different contributions, namely, Fermi contact, FC, spin-dipolar, SD, paramagnetic spin-orbit, PSO, and diamagnetic spin-orbit, DSO, eq. (1).
JNM = FCJNM + SD JNM + PSO JNM + DSO JNM (1)
There is no experimental form to determine separately the influence of each contribution. The accurate theoretical calculation of these spectroscopic parameters is a real challenge for Quantum Chemistry since their calculation is computationally very demanding. For this reason the most evolved post-Hartree-Fock approaches can, at present, only be applied to small compounds containing only light atoms[4,5]. In the mid 1990's the first calculations of SSCCs within the DFT framework were reported, [6,7,8] and a few years later, using the coupled-perturbed (CP-DFT) approach, a large breakthrough in such calculations was achieved [9,10,11]. Within this approach reasonable good results of SSCCs for medium sized compounds can be calculated .
In many instances different approaches [3,14,15,16] were developed in order to study how the nucleus-electron magnetic interaction is transmitted between regions inside a molecule, or from one molecule to another in a complex. With those approaches it is sought to increase the ability of the high resolution NMR spectroscopy to provide invaluable information about fine details of molecular electronic structures. Although in this paper a better understanding of the transmission in long-range SSCCs saturated compounds is sought, a different approach is followed, based in the following considerations. Many SSCCs are now known to be by far dominated by the FC interaction, although there are also known conspicuous exceptions to this statement [17, 18]. Owing to the importance of the FC term special emphasis was made to understand how this interaction is transmitted through the electronic molecular system. It is known that the FC interaction induces a slight electron spin polarization, which is transmitted through the electronic system by exchange interactions (spin correlation), i.e. its transmission is closely related to that of the so called "Fermi hole" . This suggests that the overlap between two occupied localized molecular orbitals, LMOs, is an efficient pathway for transmitting the FC spin information . There are many examples where this condition can be observed, e.g. when two moieties are spacially close. In this particular case it constitutes what is known in the literature as the "through-space" transmission of the FC term . It can also be observed for two LMOs representing two bonds which share one atom; this mechanism should be important for transmission of, for instance, geminal SSCCs. This phenomenon attenuates rapidly when increasing the number of bonds separating the coupling nuclei, it would be very difficult to observe any measurable effect beyond three bonds (vicinal couplings) if this is the only mechanism for transmitting the FC term. Why then long-range couplings, nJNM, with n > 3, are frequently observed?. Although electron delocalization is known to play an important role in transmitting the FC interaction , it is better answering this question by considering separately, a) unsaturated compounds, b) saturated compounds.
For unsaturated compounds this question was answered many years ago [22,23]; for saturated compounds we have recently addressed this topic  which, until few years ago, it was a kind of a puzzle . In saturated compounds delocalization interactions are not so well known as in unsaturated systems; for this reason in general long range couplings in saturated compounds are not so well known as, for instance, in conjugated systems. However, s-hyperconjugative interactions can be quantitatively described by the NBO approach of Weinhold et al.  and, therefore, the ability of different saturated compounds to transmit long-range couplings can be assessed, at least qualitatively, comparing adequately long-range SSCCs with such interactions. In this work we study long-range SSCCs 4J(Ca,H3) transmitted through s-delocalization like those described recently when rationalizing "trans-cage" SSCCs in 1-F,4-X-cubanes , I. It was observed that this system is very efficient for transmitting the FC term since there are six equivalent coupling pathways, and each of them is equally efficient owing mainly to these two factors: the arrangement of bonds is such that they favor electron delocalization interactions; each of these interactions is important since the substrate of these compounds is strained and therefore, electron delocalization interactions are enhanced . In a strained cage hydrocarbon the C—C bonds are very good electron donors, while the corresponding (C—C)* antibonds are rather poor electron acceptors. In scheme 3 the six equivalent coupling pathways for transmitting the FC term of 5J(F1,Ca) SSCC in I are sketched, highlighting the six donor bonds that undergo important interactions of types (C—C) —> (C1—F)* and (C—C)—>(C1—X)*. With this idea in mind it is possible to rationalize, for instance, 7J(F,F) = 2.8 Hz in difluorodiadamantane , where the suggested three-coupling pathways for transmitting such a long-range SSCC are highlighted, scheme 4.
In this paper the "trans-cage" long-range 4J(13Ca,1H3) SSCCs are studied in several members of series 1 X-bicyclo[1.1.1] pentanes, II, where the a atom of substituent X is 13C. Since these experimental SSCCs were not reported , they were also measured for X = CH3, CH2OH, t-Bu, Ph, COCH3, CO2H, CO2CH3, CN as part of this work. For five members of this series, X = CH3, CH2OH, COCH2, CO2H, CN, i.e. for compounds containing the simplest side-chains, the respective four contributions to 4J(13Ca,1H3) SSCCs were calculated at the B3LYP/EPR-III level. The relevant s-hyperconjugative interactions were studied using the Weinhold's Natural Bond Orbitals, NBO, approach .
Figure 1: 1-F,4-X-cubane, I.
Figure 2: 1-X-bicyclo[1.1.1] pentane, II
Figure 3: Coupling pathways for 5J(F1,X) in I
Figure 4: 7J(F,F) = 2.8 Hz in difluorodiadamantane
Experimental and Computational Details
3a. Experimental details
Syntheses of compounds of series II (X= CH3, CH2OH, t-Bu, Ph, COCH3, CO2H, CO2CH3, CN) were described elsewhere . Samples were prepared in CDCl3 at concentrations of ca 0.6 mol -1. 13C NMR spectra were recorded on a Varian Gemini 300BB instrument, operating at 75.462 MHz. 4J(13Ca,1H) SSCCs were determined from proton-coupled spectra, measured at 100-500 Hz spectral width with 30-40K digital resolution. Experimental couplings thus obtained are considered to be accurate to ± 0.1 Hz. These measurements were carried out in the Department of Chemistry, Faculty of Science and Engineering, Flinders University of South Australia.
3b. Computational details
The geometries for five members of series II, were optimized using the hybrid B3LYP functional, which corresponds to Lee et al. correlated functional  and the exchange part is treated according to the Becke's three parameter approach [31,32]. For such optimizations the 6-311G(d,p) basis set was chosen. Calculations of all four terms of SSCCs, i.e. FC, SD, PSO and DSO, were carried out using the B3LYP functional and the EPR-III basis set . This basis set was chosen since SSCCs calculated at the B3LYP/EPR-III level are known to be close to the basis-set converged limit . The EPR-III basis set is of a triple-zeta quality and includes diffuse and polarization functions. The s part of the EPR-III basis set is enhanced to better reproduce the electronic density in the nuclear regions; this point is particularly important when calculating the FC term. The CP-DFT perturbative approach [9,10] as implemented in the Gaussian 03 package of programs  was used for calculating all the three second-order terms of spin-spin couplings, i.e. FC, SD, and PSO; the DSO term is treated as a first-order quantity and therefore it is evaluated as an expectation value over the electronic ground state. All DFT calculations were carried out with the Gaussian 03 package of programs . s-Hyperconjugative interactions were studied using the NBO approach .
Results and Discussion
In table 1 the experimental values measured in this work for the "trans-cage" 4J(13Ca,1H) SSCCs measured in 8 members of series II are displayed. These values are compared with the corresponding 4J(13Ca,19F) SSCCs (taken from Ref. ) in a series similar to II, but where the H atom placed at position 3, has been replaced by an F atom, i.e. 1-X;3-F-bicyclo[1.1.1] pentanes, II-F. This extra series of experimental values was added to table 1 since it was considered convenient to observe that the sensitivity to a given substituent is notably enhanced in the substrate with F at position 3 as compared with the substrate with H at the same position. It is worth noting that there is no exact parallelism in these two series of values, although they follow the same trend. SSCCs in the series II-F were not calculated in this work due to the difficulty to obtain these parameters when at least a fluorine atom is involved . However, some considerations will be made to rationalize the larger efficiency for transmitting the trans-cage couplings. To obtain a better idea about the sensitivity for transmitting the FC term in "trans-cage" SSCCs, it is interesting to consider also the parent hydrocarbon, i.e. bicyclo[1.1.1]pentane where it is known that 4J(H1,H3) = 18 Hz , and compound 1-F-bicyclo[1.1.1] pentane, where 4J(F1,H3) = 71 Hz  ( = 69.71 Hz according to Ref. ). At this point it is important to recall that, when comparing SSCCs involving different isotopic species, it must be taken into account that JAB SSCC is proportional to the gyromagnetic ratios gA and gB of nuclei A and B, respectively. In this way 4J(H,H), 4J(C,H), 4J(C,F) SSCCs can be converted into "equivalent" 4J(H,H)eq SSCCs, multiplying, respectively, by ; gH/gH = 1; gH/gC = 3.9761; and (gH/gF)2 = 1.0624.
Table 1. Comparison of experimental trans-cage 4J(Ca,X3) couplings in series II and II-F.
All SSCCs are given in Hz.
(a) This work
b)Taken from Ref. .
In table 2 calculated 4J(Ca,H3) SSCCs in five members of series II are compared with the respective experimental values. All four contributions to such SSCCs are explicitly shown to convey the idea of the relative importance of each contribution, i.e. they are by far dominated by the FC term. It is also noteworthy that the PSO and DSO terms almost cancel each other, a situation observed for several types of SSCCs, although there are also many cases where such cancellation does not take place. Their values are shown with two decimal figures to stress the low sensitivity of non-contact contributions to the X-substituent. For comparing with experimental values, total calculated SSCCs are rounded off to only 1 decimal figure. For each of the compounds with X = CH2OH; COCH3; and CO2H two conformations were considered, as shown in table 2. In general the agreement between calculated and experimental couplings is very good, observing the largest difference for X = CN; however, the experimental trend is correctly reproduced.
Table 2. Comparison between calculated 4J(Ca,H3) SSCCs (in Hz) with their experimental values, measured as part of this work.
a) This corresponds to the preferential conformation with the dihedral angle OCaC1C2 = 180°.
b) The second conformer corresponds to OCaC1C2 = 0°; it is 2.9 kcal/mol above the preferential conformation.
c) This corresponds to the preferential conformation with the dihedral angle OCaC1C2 = 0°.
d) The second conformer corresponds to OCaC1C2 = 180°; it is 0.7 kcal/mol above the preferential conformation.
e) This corresponds to the preferential conformation with the dihedral angle OCaC1C2 = 0°.
f) The second conformer corresponds to OCaC1C2 = 180°; it is 0.4 kcal/mol above the preferential conformation.
It is now well known that in the bicyclo[1.1.1]pentane cage the six C—C bonds, C1—C2, C1—C4 and C1—C5, i.e. C1—Ci; and C2—C3, C4—C3 and C5—C3, i.e. Ci—C3, are excellent electron donors, and their corresponding antibonds, (C—C)*, are bad electron acceptors for s-hyperconjugative interactions. For this reason, we can expect that the main hyperconjugative interactions transmitting the spin information associated to the FC term from Ca to H3 (or, equivalently, from H3 to Ca) should be (C1—Ci)—>(C3—H)*, and (Ci—C3)—>(C1—Ca)*, interactions called of type 1 and type 2, respectively. It is interesting to note that for an X group with a three-fold symmetry axes the three interactions of type 1 are equivalent, the same holds for the three (Ci—C3)—>(C1—Ca)* interactions. From an experimental point of view an effective equivalence of this type holds since in all cases the barrier for inner rotations around the C1—Ca bond are low enough and all compounds displayed in table 1 show an effective three-fold symmetry axes at room temperature. For this reason, in table 3 the sum of the three interactions of type 1 and the sum of the three interactions of type 2, are given in the first and second columns, respectively.
Table 3. Electron delocalization interactions 1 and 2 (in kcal/mol) expected to be efficient for transmitting the FC term of 4J(Ca,H3) SSCCs in series II
. The Ca s % character of the C1—Ca bond is also shown since the actual value of this coupling also depends on it.
a 1: Sum of the three (C1—Ci)—>(C3—H)* interactions (I = 2, 4 and 5).
b 2: Sum of the three (Ci—C3)—>(C1—Ca)* interactions (i = 2, 4 and 5).
The Ca s % character of the C1—Ca is also shown since the actual value of 4J(Ca,H3) SSCC depends on it . It is observed that the sum of interactions (C1—Ci)—>(C3—H)*, 1, is much less sensitive to substituents placed at position 1 than the sum of interactions (Ci—C3)—>(C1—Ca)*, 2. For comparison purposes in table 4 the corresponding interactions in the parent hydrocarbon, bicyclo[1.1.1]pentane, in II(X = CH3) and in 1-F-bicyclo[1.1.1]pentane are also shown, together with the respective experimental values of trans-cage SSCCs, together with their "equivalent" J(H,H) values. Due to symmetry reasons, in bicyclo[1.1.1]pentane the (C1—Ci)—>(C3—H)* interactions are exactly the same as the (Ci—C3)—>(C1—Ca)* interactions. It is noted that the former are close to (C1—Ci)—>(C3—H)* in compounds of series II; which indicates that the main difference in the ability for transmitting the FC interaction corresponding to 4J(X,H3) can be ascribed to the different (C1—X)* antibond acceptor capability. It is noted that, according to Table 4, the (C1—F)* antibond in 1 F bicyclo[1.1.1]pentane, is notably a better electron acceptor than (C1—Xa)* antibond in series II.
Table 4: Comparison of s-hyperconjugative interactions of types 1 and 2 in bicyclo[1.1.1]pentane; in II(X = CH3) and in 1-F-bicyclo[1.1.1]pentane.
The experimental trans-cage 4J(X,H3) SSCCs (in Hz) are also compared.
a 1: Sum of the three (C1——Ci)—>(C3—H)* interactions (i = 2, 4 and 5).
b 2: Sum of the three (Ci—C3)—>(C1—Ca)* interactions (i = 2, 4 and 5).
c J(H,H)eq = (gH/gX)4J(X,H3).
d Taken from Ref. .
e Taken from Ref. .
In the series of compounds 1-X-bicyclo[1.1.1]pentanes trans-cage 4J(13Ca,1H3) SSCCs studied in this work are new examples of the efficiency for transmitting the FC spin information through a strained cage compound. The capability for these cage compounds for transmitting substituent effects on such couplings is fundamentally determined by the different electron acceptor ability of the (Xa—C1)* antibonding orbital for different substituents. On the other hand, the better ability of compounds 1-X;3-F-bicyclo[1.1.1]pentanes (see table 1) than 1-X-bicyclo[1.1.1]pentanes originate in the much better electron acceptor ability of the (C3—F)* antibond in the former than that of the (C3—H)* antibond in the latter.
Since electron delocalization interactions take place in all three bridges of these bicyclo-compounds, it is observed that SSCCs transmitted through strained cage compounds, should follow a multipath additivity. It is important to recall that, recently, the multipath additiviy in cyclopropane and bicyclobutane has been questioned .
RHC gratefully acknowledges financial support from CONICET (PIP 2140) and UBACyT (X-222). ED and ALE gratefully acknowledge support from the "Dirección General de Enseñanza Superior e Investigación Científica" (BQU2000-0211-CO2).
EWD and IJL are grateful to the Australian Research Council for financial support.
 Emsley, J. W.; Feeney, J.; Sutcliffe, L. H. High Resolution Nuclear Magnetic Resonance Spectroscopy, Pergamon Press, Oxford, (1966). [ Links ]
 Ramsey, N. F.; Purcell, E. M., Phys. Rev. 1952, 85, 143. [ Links ]
 Ramsey, N. F., Phys. Rev. 1953, 91, 303. [ Links ]
 Helgaker, T.; Jaszunski, M.; Ruud, K., Chem. Rev., 1999, 99, 293. [ Links ]
 Fukui, H., Prog. NMR Spectrosc., 1999, 35, 267. [ Links ]
 Malkin, V. G.; Malkina, O. L.; Salahub, D.R., Chem. Phys. Lett., 1994, 221, 91. [ Links ]
 Malkina, O. L.; Salahub, D. R.; Malkin, V.G., J. Chem. Phys., 1996, 105, 8793. [ Links ]
 Dickson, R. M.; Ziegler, T., J. Phys. Chem., 1996, 100, 5286. [ Links ]
 Sychrovský, V.; Gräfenstein, J.; Cremer, D., J. Chem. Phys., 2000, 113, 3530. [ Links ]
 Helgaker, T.; Watson, M.; Handy, N. C., J. Chem. Phys., 2000, 113, 9402. [ Links ]
 Barone, V.; Peralta, J. E.; Contreras, R. H.; Snyder, J. P., J. Phys. Chem. A, 2002, 106, 5607. [ Links ]
 Contreras, R. H.; Barone, V.; Facelli, J. C.; Peralta, J. E., Ann. Reps. NMR Spectrosc., 2003, 51, 167. [ Links ]
 Engelmann, A. R.; Contreras, R. H.; Facelli, J. C., Theoret. Chim. Acta, 1981, 59, 17. [ Links ]
 Wu, A.; Cremer, D., Phys. Chem. Chem. Phys., 2003, 5, 4541. [ Links ]
 Wu, A.; Gräfenstein, J.; Cremer, D., J. Phys. Chem. A, 2003, 107, 7043. [ Links ]
 Gräfenstein, J.; Cremer, J. D., Chem. Phys. Lett., 2004, 387, 415. [ Links ]
 Barone, V.; Provasi, P. F.; Peralta, J. E.; Zinder, J. P.; Sauer, S. P. A.; Contreras, R. H., J. Phys. Chem. A, 2003, 107, 4748. [ Links ]
 Wrackmeyer, B., Estruct. Chem., 2005, 16, 67. [ Links ]
 Soncini, A.; Lazzeretti, P., J. Chem. Phys. 2003, 118, 7165. [ Links ]
 Arnold, W. D.; Oldfield, E., J. Am. Chem. Soc., 2000, 122, 12835. [ Links ]
 Matta, C. F.; Hernández -Trujillo, J.; Bader, R. F. W., J. Phys. Chem. A, 2002, 106, 7369. [ Links ]
 Hoffman, R.A., Molec. Phys., 1958, 1, 326. [ Links ]
 Barfield, M.; Chakrabarti, B., 1969, 57, 757. [ Links ]
 Contreras, R. H.; Esteban, A. L.; Díez, E.; Head, N. J.; Della, E. W. ,Molec. Phys., 2006, 104, 485. [ Links ]
 Gakh, Y. G.; Gakh, A. A.; Gronenborn, A. M., Magn. Reson. Chem., 2000, 38, 551. [ Links ]
 Reed, A. E.; Curtiss, L. A.; Weinhold, F., Chem. Rev., 1988, 88, 899. [ Links ]
 Contreras, R. H.; Peralta, J. E., Prog. NMR Spectrosc., 2000, 37, 321. [ Links ]
 Olah, G. A.; Shih, J. G.; Krishnamurthy, V. V.; Singh, B. P., J. Am. Chem. Soc., 1984, 106, 4492. [ Links ]
 Della, E. W.; Lochert, I. J.; Peruchena, N. M.; Aucar, G. A.; Contreras, R. H. J. Phys. Org. Chem., 1996, 9, 168. [ Links ]
 Lee, C.; Yang, W.; Parr, R. G., Phys. Rev. B, 1988, 37, 785. [ Links ]
 Becke, A. D., Phys. Rev. A, 1988, 38, 3098. [ Links ]
 Becke, A. D., J. Chem. Phys. 1993, 98, 5648. [ Links ]
 Barone, V., J. Chem. Phys., 1994, 101, 6834. [ Links ]
 Peralta, J. E.; Scuseria, G. E.; Cheeseman, J. R.; Frisch, M., J. Chem. Phys. Lett. 2003, 375, 452. [ Links ]
 Gaussian 03, Revision B.01, 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. B.; 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.; Gonzalez, C. and Pople, J. A. Gaussian 03, Revision B.05 2003, Gaussian, Inc., Pittsburgh PA. [ Links ]
 NBO version 3.1; E. D. Glendening, A. E. Reed, J. E. Carpenter, F. Weinhold, (Included in the GAUSSIAN 03 package of programs). [ Links ]
 Adcock, W.; Krstic, A. R., Magn. Reson. Chem., 2000, 38, 115. [ Links ]
 Jaszunski, M.; Chem. Phys. Let., 2004, 385, 122, and references therein. [ Links ]
 Wiberg, K. B., Connor, D. S. J. Am. Chem. Soc., 1966, 88, 4437. [ Links ]
 Barfield, M.; Della, E. W.; Pigou, P. E.; Walter, S. R., J. Am. Chem. Soc., 1982, 104, 3549. [ Links ]
 Contreras, R. H.; Ruiz de Azúa, M. C.; Giribet, C. G.; Aucar, G. A.; Lobayan de Bonczok, R., J. Mol. Struct. (Theochem), 1993, 284, 249. [ Links ]
 Soncini, A.; Lazzeretti, P., ChemPhysChem., 2006, 7, 679. [ Links ]