versión On-line ISSN 0717-9707
J. Chil. Chem. Soc. v.54 n.1 Concepción 2009
J. Chil. Chem. Soc, 54, N° 1 (2009); págs: 43-45
IMPROPER HYDROGEN BONDS - A THEORETICAL STUDY ABOUT THE MOLECULAR STRUCTURE OF INTERMOLECULAR SYSTEMS FORMED BY H3CH+Δ...X AND H3C+Δ...HY WITH X = CL- OR F- AND Y = CL OR F
Departamento de Ciências Farmacéuticas Universidade Federal de Pernambuco 50740-520, Recife, PE-Brazil. e-mail: email@example.com
In this work is presented a theoretical study of the molecular properties of the H3CH...X and H3C+δ...HY intermolecular systems with X = CL- or F- and Y = Cl or F. In the H3CH...X complex, it is formed a traditional hydrogen bond between the CL- or F- anions and the hydrogen atom of the methyl. About the H3C+δ...HY complex, it was observed an improper hydrogen bond because the carbon atom of the methyl cation function as a proton acceptor. In this insight, the capacity of methane to interact with halogen anions (CL- or F-) and molecular acids (HCl or HF) was examined at the B3LYP/6-311++G(3df,3dp) level of calculation. Moreover, the interaction strengths of the H3CH...X and H3C+δ...HY complexes was evaluated by computing their intermolecular distances and binding energies. Finally, QTAIM calculations also were executed with the purpose to examine the intermolecular interactions through the quantification of their electronic densities (p) as well as by the interpretation of the Laplacian operators ().
Keywords: Methane; anions; acids; hydrogen bonds; improper hydrogen bonds.
It is known by the scientific community the existence of several intermolecular interaction types, which can be represented by hydrophobic, dipole-dipole and hydrogen bonding1. Characterized by an interaction between a center containing high electronic density concentration and often an electropositive atom2, the hydrogen bonding plays a vital role in studies of chemical, physical, and biological systems 3. However, few years ago a new class of intermolecular systems known as Dihydrogen-Bonded (DHB)4 complexes has been discovered. These systems are formed through the ability of hydrogen atoms to form an intermolecular contact with protonic centers. In other words, notonecessarily the hydrogen atoms interact with high electronic sites such as lone electron pairs or π bonds 5.
Through the X-rays diffraction analysis, the Crabtree's research group 6 has shown that DHB complexes are formed by proton donors (H+δ) and hydride compounds (H-δ) 7,8. On the other hand, by considering the formation of the NH+δ...H-δY DHB complex, where normally N is an electronegative group and Y is a alkali metal (AM) or alkali earth metals (ΔEM)9-14, some theoretical investigations documented by Grabowski 15 has revealed the existence of H+δ...H-δY DHB complexes formed by Lithium, Sodium, Magnesium, and a different proton acceptor: the methane. Furthermore, it was also verified that the FH+δ...H-δCH3 DHB complex is formed through the interaction of the methane with hydrofluoric acid. In addition, Grabowski and co-workers 16 demonstrated that the methane also interact with alkali metals and alkali earth metals , where the carbon atom is sufficiently electronegative to form electrostatic interactions with Beryllium or Sodium, for instance. However, the methane also presents an opposite behaviorto that observed in DHB complexes, where CL- and I- interact with the hydrogen atoms of the CH bond. Thereby, it is then formed the H3CH+δ...CL- and H3CH+δ...I- complexes. By all these reasons, definitively the methane can be considered a convergence point by taking into account ionic and neutral intermolecular systems, of which we can mention those formed by F- and CL- entities, or else molecules such as hydrochloric acid.
About all aspects cited above, in this work is proposed a theoretical study about the capacity of methane to interact with donors and acceptors of charge density. Such study will be developed through the examination of molecular properties of the H3CH...X and H3C+δ...HY intermolecular systems with X = CL- or F-, and Y = Cl or F. For this proposal, I will evaluate structural parameters and electronic properties with the purpose to describe the intermolecular strength through the measurement of distance, R(H+δ...X), and binding energy, ΔE. Moreover, another electronic evaluation will be effectuated here by means of the Quantum Theory of Atoms in Molecules (QTAIM)17. From QTAIM calculations, the chemical topology is modeled by a numeric integration of the electronic density, allowing then the description of the molecular interactions through the identification of Bond Critical Points (BCP)18 between neighboring atoms. The BCP are the basis used by QTAIM algorithm to identify maxima and minima of electronic density. Thus, covalent bonds containing high electronic density are identified as shared interactions ls and, centers with low electronic densities, in special those identified as being intermolecular and intramolecular contacts are then known as closed-shell interactions, of which the hydrogen bonding is considered a classic example20-25. So, QTAIM analyzes the interaction strength in hydrogen bonded-complexes through the topological parameters, such as the electronic density, ρ, and the Laplacian field, . The Laplacian at the BCP provides a measure of the extent to which the charge density is locally compressed or expanded in the chemical bonds. About this, there are two important criteria often used to characterize the chemical bonds according to QTAIM approach: i) in the case of locally concentrated charge density, <0, it is indicated a shared interaction, e.g., covalent or unsaturated bonds. ii) When the charge density is locally depleted, >0, it is identified a closed-shell interaction, such as hydrogen bonds or else ionic bonds.
Through the application of the B3LYP/6-311++G(3df,3dp) level of theory implemented in GAUSSIAN 98W software26, the calculations of the optimized geometries of the H3CH...X and H3C+δ...HY complexes were performed using default convergence criteria. The hybrid B3LYP 27 of the Density Functional Theory (DFT) 28 was chosen because it has been documented a lot of studies of hydrogen complexes 29 in which their molecular parameters were successfully evaluated by means of hybrid density functionals 30-31. The supermolecule approach 32 was used to determine the values of the binding energies ΔE, which were corrected by the calculations of the Basis Sets Superposition Error (BSSE)33 performed through the MASSAGE keyword in the GAUSSIAN 98W program. Finally, the results of the corrected binding energies ΔEC were computed according to the equation (1), and all QTAIM calculations were executed in the GAUSSIAN 98W program.
RESULTS AND DISCUSSION
Geometry and energy
The Figure 1 illustrates the optimized geometries of the H3CH...Cl- (a), H3CH...F- (b), H3C+δ...HCl (c) and H3C+δ...HF (d) complexes obtained through the B3LYP/6-311++G(3df,3dp) calculations. In Table 1 are listedthe main structural parameters of these complexes, such as their intermolecular distances, R(H+δ...X) and R(C+δ...HY), and deformations on methane molecule, SrCH. Moreover, in Table 1 are also presented the values of the corrected binding energies, ΔEC, as well as the results of the non-corrected binding energies, ΔE, and BSSE amounts.
Firstly, analyzing the R(H+δ...X) intermolecular distances of the H3C H+δ...Cl- and H3CH+δ...F- complexes, the values of 1.9823 Å and 0.9340 Å indicates that the fluorine ion provides a stronger bound system, once H3C+δ...HCl and H3C+δ...HF complexes present longer intermolecular distances, whose values of R(C+δ...HY) are 2.072 Å and 2.1098 Å, respectively. As it is well established that hydrogen bonds often presents typical distance value in range of 2.5 Å5, the results of 0.9340 A of the H3CH+δ...F- is considered a very short interaction. However, an interesting aspect related to the H3CH+δ...F- must be emphasized. Note that the methane structure assumes a planar configuration, or in other words, there is a tendency to form the CH3+ methyl cation and hydrofluoric acid molecular. Indeed, by means of the binding energy of - 116.3 kJ mol-1, the H3CH+δ...F- complex is more stable than H3CH+δ...Cl-, which presented a value for ΔEC of - 4.1 kJ mol-1. As the H3CH+δ...F- complex presented the stronger interaction, it was verified drastic alterations on its structure, e.g., the variation of 0.9955 Å of the CH+δ bond in comparison with results of 1.0881 Å and 2.0837 Å of the methane isolated and when the H3CH+δ...F- complex is formed, respectively.
By giving continuity the discussion of the molecular properties of the H3CH...X and H3C+δ...HY complexes, the Table 2 present the results of the QTAIM calculations. According to the low values of ρ and positive results of , the (H+δ...-X) and (C+δ...HY) contacts were identified as closed-shell interactions 17, in other words, concentration of charge density in separate nuclei. Moreover, I should assume that (H+δ...X) and (C+δ...HY) are considered ionic hydrogen bonds, which are improper types because the carbon of methane is a proton acceptor (c and d complexes), as will be discussed later. However, the interaction of the fluorine ion with the methane still provides special aspects, once the electronic density of 0.365 e.u. for H+δ...F- is higher in comparison with the H+δ...Cl-, C+δ...HCl and C+δ...HF interactions, whose values are in range of 0.021-0.041 e.u. In this sense, the stronger interaction of fluorine ion with the methane suggests another kind of closed-shell interaction, which is formed between the carbon and the H+δ hydrogen atom of HF, as can be seen in the BCP depicted in Figure 2 (b).
Thus, the values of ρ and of 0.021 e.u. and 0.042 e.u. indicates the co-existence of a improper hydrogen bond C...H+δF, although the value of - 3.2 e.u. shows that H+δ and F- interact as a shared contact. In this situation, the electronic density of 0.365 e.u. of the H+δF bond is very similar to 0.369 e.u. of the hydrofluoric acid monomer. At last, corroborating with the results presented above, the formation of H+δ...F- provokes substantial changes in the electronic structure of the intermolecular system. For instance, the weakness of the CH+δ bond has caused the formation of the improper hydrogen bond C...H+δF, whose electronic density of 0.021 e.u. is lower in comparison with 0.029 e.u. of the CHß bond.
The capacity of methane to interact with anions CL- and F- as well as molecular acids HF and HC1 was demonstrated in this work. Due to the accentuated electronegative character of the fluorine, we observe that besides H+δ...F- contact, a C...H+δF improper hydrogen bond was formed. Such aspect was demonstrated by the QTAIM calculations, where positive and negative Laplacian results were quantified between C and H+δ, as well as between H+δ and F. In this last case, the calculated electronic density value of 0.365 e.u. is very close to 0.369 e.u. of the hydrofluoric acid monomer.
The author would like to gratefully to the CAPES and CNPq Brazilian funding agencies. Moreover, I would like also to very gratefully to CENAPAD-SP (Centro Nacional de Processamento de Alto Desempenho - Sao Paulo) by computational facilities to complete this work.
1. T. W. Martín, Z. S. Derewenda, Nat. Struct. Bio. 6, 403, (1999). [ Links ]
2. G. A Jeffrey, W. Saenger, Hydrogen Bonding in Biological Structures, Springer-Verlag, Berlin, 1991. [ Links ]
3. G. R. Desiraju, T. Steiner, The Weak Hydrogen Bond in Structural Chemistry and Biology, Oxford University Press, New York, 1999. [ Links ]
4. R. Bosque, F. Maseras, O. Eisenstein, B. P. P. Patel, W. Yao, R.H. Crabtree, Inorg. Chem. 36, 5505, (1997). [ Links ]
5. I. Rozas, I. Alkorta, J. Elguero, J. Phys. Chem A 101, 9457, (1997). [ Links ]
6. T. B. Richardson, S. de Gala, R. H. Crabtree, P. E. M. Siegbahn, J. Am. Chem. Soc. 117, 12875, (1995). [ Links ]
7. B. G. Oliveira, R.C.M.U. Araújo, M.N. Ramos, Struct. Chem. 19, 185, (2008). [ Links ]
8. B. G. Oliveira, R. C. M. U. Araújo, M. N. Ramos, Struct. Chem. 19, 665, (2008). [ Links ]
9. Q. Liu, R. Hoffman, J. Am. Chem. Soc. 117, 10108, (1995). [ Links ]
11. E. Peris, J. C. Lee, J. R. Rambo, O. Eisenstein, R. H. Crabtree, J. Am. Chem. Soc. 117, 3485 (1995). [ Links ]
12. R. H. Crabtree, Acc. Chem. Res. 23, 95, (1990). [ Links ]
13. R. H. Crabtree, P. E. M. Siegbahn, O. Eisenstein, A. L. Rheingold, T. F. Koetzle, Acc. Chem. Res. 29, 348, (1996). [ Links ]
14. J. Wessel, J. C. Lee, E. Peris, G. P. A. Yap, J. B. Fortín, J. S. Ricci, G. Sini, A. Albinati, T. F. Koetzle, O. Eisenstein, A. L. Rheingold, R. H. Crabtree, Angew. Chem., Int. Ed. Engl. 34, 2507, (1995). [ Links ]
15. S. J. Grabowski, J. Phys. Chem. A 104, 5551, (2000). [ Links ]
16. S. J. Grabowski, W. A. Sokalski, J. Leszczynski, Chem. Phys. Lett. 422, 334, (2006). [ Links ]
17. R. F. W. Bader, Atoms in Molecules. A Quantum Theory, Clarendon, Oxford, 1990. [ Links ]
18. R. F. W. Bader, Chem. Rev. 91, 893, (1991). [ Links ]
19. R. F. W. Bader, P. J. MacDougall, C. D. H. Lau, J. Am. Chem. Soc. 106, 1594,(1984). [ Links ]
20. B. G. Oliveira, F. S. Pereira, R.C.M.U. Araújo, M. N. Ramos, Chem. Phys. Lett. 427, 181, (2006). [ Links ]
21. B. G. Oliveira, M. L. A. A. Vasconcellos, J. Mol. Struct. (THEOCHEM) 114, 83, (2006). [ Links ]
22. B. G. Oliveira, R. C. M. U. Araújo, A. B. Carvalho, E. F. Lima, W. L. V. Silva, M. N. Ramos, A. M. Tavares, J. Mol. Struct. (THEOCHEM) 775, 39, (2006). [ Links ]
23. B. G. Oliveira, R. C. M. U. Araújo, F. F. Chagas, A. B. Carvalho, M. N. Ramos, J. Mol. Model 14, 949, (2008). [ Links ]
24. B. G. Oliveira, R. C. M. U. Araújo, Quim. Nova 30, 791, (2007). [ Links ]
25. B. G. Oliveira, R. C. M. U. Araújo, F. P. Silva, E. F. Lima, E. F. Lima, W. L. V. Silva, A.B. Carvalho, M. N. Ramos, Quim. Nova. 31, 1673, (2008). [ Links ]
26. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G.E. Scuseria, M. A. Robb, J. R. Cheeseman, V. G. Zakrzewski, J. A. Montgomery, Jr., R. E. Stratmann, J. C. Burant, S. Dapprich, J. M. Millam, A. D. Daniels, K. N. Kudin, M. C. Strain, O. Farkas, J. Tomasi, V. Barone, M. Cossi, R. Cammi, B. Mennucci, C. Pomelli, C. Adamo, S. Clifford, J. Ochterski, G. A. Petersson, P. Y. Ayala, Q. Cui, K. Morokuma, N. Rega, P. Salvador, J. J. Dannenberg, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. Cioslowski, J. V. Ortiz, A. G. Baboul, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. Gomperts, R. L. Martín, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gilí, B. Johnson, W. Chen, M.W. Wong, J. L. Andrés, C. González, M. Head-Gordon, E. S. Replogle, J. A. Pople, Gaussian 98W (Revisión A. 1), Gaussian, Inc, Pittsburgh PA, 1998. [ Links ]
27. A.D. Becke, J. Chem. Phys. 107, 8554, (1997). [ Links ]
28. P. Hohenberg, W. Kohn, Phys. Rev. B 136, 864, (1964). [ Links ]
29. J. B. P. da Silva, M. R. Silva Jr., M. N. Ramos, J. Braz. Chem. Soc. 16, 844, (2005). [ Links ]
30. M. A. Muñoz, M. Galán, L. Gómez, C. Carmona, P. Guardado, M. Balón, Chem. Phys. 290, 69, (2003). [ Links ]
31. L. F. Holroyd, T. van Mourik, Chem. Phys. Lett. 442, 42, (2007). [ Links ]
32. F. B. van Duijneveldt, J. N. Murrell, J. Chem. Phys. 46, 1759, (1967). [ Links ]
33. S. F. Boys, F. Bernardi, Mol. Phys. 19, 553, (1970). [ Links ]
(Received 13 March 2008 - Accepted 29 October 2008)