Open access peer-reviewed chapter

Possibility of Complexation of the Calix[4]Arene Molecule with the Polluting Gases: DFT and NCI-RDG Theory

Written By

Bouzid Gassoumi, Fatma Ezzahra Ben Mohamed, Houcine Ghalla and Rafik Ben Chaabane

Submitted: June 24th, 2020 Reviewed: September 1st, 2020 Published: November 25th, 2020

DOI: 10.5772/intechopen.93838

Chapter metrics overview

366 Chapter Downloads

View Full Metrics


The calix[4]arenes (abbreviated as CX[4]) are characterized by a specific hydrophobic cavity formed by a four cyclically phenol groups to encapsulate a gas or small molecules. Recently, the CX[4] molecule is used in a specific media and in pharmaceutical drug delivery. The pollution problem will be a vital subject in the future because the increase of the explosions of the gaseous pollutants in the environment. In this report, we have encapsulated the polluting gases NO3, NO2, CO2 and N2 by the calix[4]arene molecule. In this work, The binding energies of the CX[4]-gas has been calculated including the BSSE (Basis Set Superposition Error) counterpoise (CP). The red-shift of the O-H bonding interactions obtained by adding the gas in the sensitive area of calix[4]arene is clearly explained by the infrared spectrum analysis. The Molecular electrostatic potential (MEP) of the stable CX[4]-gas complexes have been investigated in the endo-vs. exo-cavity regions. Finally, the non-covalent interactions analyses of the stable host-guests complexes have been estimated by using DFT calculations.


  • complexation
  • specific gases
  • H-bonding
  • molecular electrostatic potential and binding energy

1. Introduction

The recognition of the electrostatic and magnetic properties in the selective guest in microscopic systems facilitates to know several anionic cationic or neutral guests’ complexes [1, 2, 3, 4, 5]. In this work, the CX[4] molecule is specified by its own chemical composition and the hydrophobic cavity form [6]. This molecule is characterized by specific parameters such as the diameter and the height, which facilitates the encapsulation of the cationic, anionic, neutral guests or small molecules [7, 8, 9]. Also, the cavity of the CX[4] molecule have attracted the researchers to test a new guest materials to be functional in the medical [10, 11, 12] or micro-biological field [13, 14]. In literature, we have noted that there are some works which discuss the encapsulation of the CX[4] with gases molecules (CH4, NH3 and C2H2) [15, 16, 17, 18]. In our work, we have studied the photo-physical properties of the CX[4]-gas complexes (CX[4]-NO3, CX[4]-NO2, CX[4]-CO2 and CX[4]-N2). We discuss the possibility of the encapsulation of these gases by the CX[4] molecule to show the sensibility of this molecule to the polluting gases outside or inside the cavity. We have chosen the NO3, NO2, CO2 and N2 gas as guest’s because they can be formed a dipole-dipole interactions and a CH…π hydrogen-bonding with the CX[4] molecule. By using the density functional theory (DFT) calculations, we have described the dynamic stabilities of the endo-vs. exo-cavity of the CX[4]-gas complexes. The vibrational properties of the CX[4]-gas complexes have been studied. The Molecular electrostatic potential studies of these host-guests complexes have been performed. The Non-covalent interaction via RDG function are very important to know the nature of the interactions between the specific guests and the CX[4] molecule.


2. Computational method

The stable structures of the studied systems (CX[4] and CX[4]-gas) have been calculated with the DFT method by using the B3LYP [16, 19, 20, 21] coupled to the D3BJ (empirical Becke and Johnson damping dispersion corrections) in combination with the 6-31 + G(d) basis set, as implemented in GAUSSIAN 09 package [22] and the Gauss View [23] as a visual program. We have been calculated the binding energies of the CX[4]-gas take into account the Basis Set Superposition Error (BSSE) counterpoise (CP) correction energy of Boys and Bernardi [24].

The binding energies are given by the following formula:

Δ E C X 4 gas = E C X 4 gas E C X 4 E gas + BSSE E1

where ECX[4]-gas, ECX[4] and Egas are the total energies of host-guest and host or guest molecules. The electronic parameters of the studied complexes are very effective to show the sensibility of the specific gas inside or outside the cavity. We analyzed then the infrared spectrum of these complexes using the DFT/B3LYP-D3 method. The nature of the interaction between the CX[4] molecule and the pollutant gases is better explained using the Non covalent interaction via RDG analysis [25].


3. Results and discussions

3.1 Geometry optimizations

The CX[4]-gas complexes have been optimized at the DFT/B3LYP-D3 ( Figure 1 ). In this study, we have placed the NO3, NO2, CO2 and N2 gas in the exo or endo-cavity positions. Table 1 presents the Binding energy Eb (in (kcal/mol)) of CX[4]-gas complexes. The binding energy value of the CX[4]-CO2 is equal to 21.33 kcal/mol (see Figure 1 ). We show that, the CX[4]-CO2(endo) complex is more stable than that of CX[4]-CO2(exo). In addition, from the CX[4]-NO3 complex, the NO3 gas has placed in the endo or the exo-cavity position. In this situation, we have obtained a divergence in the case of the interaction of CX[4] with NO3 outside the cavity that is why, we have tested the case where the NO3 gas perpendicular to the 4-fold axis of CX[4] and a parallel position of the NO3 (NO3 parallel to the 4-fold axis). We have been noted a very weak energy of the CX[4]-NO3(perp.) complex in comparison with the CX[4]-NO3(paral.) complex. Moreover, we have noted that the interaction between the CX[4] molecule and the NO3(paral.) is stabilized by a low dipole moment. We have calculated two positions for the guest N2: the N2 gas located on the outside of cavity and this gas located perpendicular to the 4-fold axis of CX[4] molecule. We have shown that the CX[4]-N2(endo) has a very strongest energy, which explain that this geometry is more stable than CX[4]-N2(exo) (see Table 1 ). The binding energy (Eb) of the stable complex is equal to 18.90 kcal/mol. Finally, from the CX[4]-NO2(exo) complex, we note that the NO2 gas is located in the outside of the H-link network of the phenolic hydroxyl groups. The binding energy of the CX[4]-NO2(endo) complex is equal to 20.54 kcal/mol (see Table 1 ). As we see, The binding energy of the CX[4]-NO2(endo)complex is higher than that of CX[4]-CO2(exo/endo). Therefore, The CX[4]-NO2(endo) complex is more stable than others, which is explained by the lowest dipole moment.

Figure 1.

Optimized geometries of the CX[4]-CO2(endo) (a), CX[4]-CO2(exo) (b), CX[4]-NO3(paral.) (c), CX[4]-NO3(perp.) (d), CX[4]-NO2(endo) (e), CX[4]-NO2(exo) (f*), CX[4]-NO2(exo) (g*), CX[4]-N2(endo) (h) and CX[4]-N2(exo) (i) complexes using B3LYP-D3BJ/6-31 + G(d) method (top view).

Complexes Eb BSSE Eb (with BSSE)
CX[4]-NO3 (paral.) 24.62 6.17
CX[4]-NO3(perp.) 16.51 6.05 22.56
CX[4]-N2(exo) 16.89 1.70 18.59
CX[4]-N2(endo) 16.90 1.80 18.70
CX[4]-NO2(exo)f* 18.17 1.82 19.99
CX[4]-NO2(exo)g* 18.11 1.80 19.91
CX[4]-NO2(endo) 17.83 2.71 20.54
CX[4]-CO2(endo) 18.62 2.71 21.33
CX[4]-CO2(exo) 18.10 2.71 20.81

Table 1.

Binding energy Eb (in (kcal/mol)) of CX[4]-gas complexes.

3.2 Red-shifted of O-H bonds

The reorganization of the O…H stretching vibration has been calculated after the system perturbation with the pollutant gases. Moreover, The IR spectrum calculation of the four stable complexes is given in Figure 2 . In literature, Furel et al. have been studied the experimental infrared spectrum of the CX[4] molecule. In this work, we have analyzed the stretching region of the IR spectrum between 2900 and 3500 cm−1, respectively 3254 (OH), 3168(OH asymmetric), 3045(Car.H), 2951 (asym.CH2), 2916 cm−1(sym.CH2). The infrared spectrum of the CX[4]-gas have been compared to the IR spectra of the free CX[4] molecule (See Figure 2 ). In addition, to take into account the an-harmonic effect our calculated frequencies scaled by 0.956. However, we have noted that, the CX[4]-CO2 is characterized by two peaks located in the vicinity of 3177 cm−1 and 3181 cm−1respectively. The CX[4]-N2 complex have the same results. These tow peaks are specified by the O-H asymmetric stretching vibration. We have shown that the frequency band located at 3160 cm−1 is corresponding to the O-H stretching vibration of the phenol O-H groups. Concerning the CX[4]-NO2 complex, we have noted the appearance of the two peaks in the neighborhood of 3170 cm−1 and 3193 cm−1. These peaks are corresponding to the vibrations of the O-H and O-H asymmetric bonds. Finally, we have shown a burst of the OH peak in the CX[4]-NO3 gas what form four peaks located in the vicinity of 2928 cm−1 (asym.CH2), 3100 cm−1(asym.OH), 3204 cm−1(OH) and 3298 cm−1(free OH) respectively. The noted values for the comparing of the red-shifted O-H vibration between the CX[4]-gas and the free CX[4]molecule are 44 cm−1, 24 cm−1 and 9 cm−1 for CX[4]-NO3, CX[4]-NO2, CX[4]-CO2 and CX[4]-N2 successively.

Figure 2.

Infrared spectrum of the stable host-guests complexations (H-bonding region (the unit is cm−1)).

3.3 Molecular electrostatic potential study

In Figure 3(ad) , we have been created the MEP map of the stable host-guests. These graphs indicate the relation between the supra-molecular structure and the physic-cal-chemical properties of the CX[4]-gas complexes. In this work, we have explained the more nucleophile or electrophile sites in these stable host-guests. The color code of the maps varying from −0.005 to 0.005 (isoval = 0.001). From Figure 3 , we show that the electrophilic sites surrounded by N atoms and the nucleophilic sites surrounded by O atoms. In addition, these host-guests complexes are characterized by the existence of the positive charges located at the level of the phenolic branch. This part demonstrates how the MESP can explore the region of the unravel molecular interactions between the CX[4] molecule and these hosts.

Figure 3.

Molecular electrostatic potential analysis of the CX[4]-CO2 (a), CX[4]-N2 (b), CX[4]-NO2 (c) and CX[4]-NO3 (d) complexes calculated by B3LYP-D3/6–31 + G(d) level.

3.4 Non covalent interactions analysis

Figure 4(a) shows the existence of a weak van der Waals (VdW) type interactions between CO2 gas and CX[4]. Concerning the CX[4]-N2 complex ( Figure 4(b) ), we show a green color between the guest and the host that indicates the existence of the weak Van der Waals interactions and the blue color at the lower edge level indicates the existence of the O-H-bonding type interactions. Also, we have shown clearly a red color located in the center of the phenol rings indicates a strong repulsion. The NCI-plots have been confirmed these results (see Figure 4 ). In addition, from CX[4]-NO2 ( Figure 4(c) ) complex, we find the existence of a red color which explain the steric effect interactions, blue color (hydrogen bonds type interactions) and a green color (week van der Waals type interactions). The type of majority bonds of the links between the NO3 gas and the CX[4] molecule is the weak VdW type interactions. The NCI-RDG analysis shows that the VdW type interactions and the hydrogen bonding interactions between the guest and the host are very necessary for the stability of the encapsulated complexes.

Figure 4.

NCI-RDG plots of the electron density and its reduced gradient of the inclusion complexes for CX[4]-gas (CX[4]-CO2(a), CX[4]-N2(b), CX[4]-NO2(c) and CX[4]-NO3(d)). The iso-surfaces were constructed with RGD = 0.5a.u and the colors scaling from −0.01 to −0.01 a.u.


4. Conclusion

The CX[4] and the CX[4]-gas complexes have been optimized using the density functional theory (DFT). Our work has clearly explained the sensibility of the pollutant gas inside the cavity, which is very important in comparison with the gas located outside the cavity. The IR spectrum study has explained the role of the NO3 gas in the red-shifted of the O-H bonds in comparison with the other gases. The MEP results is clearly explained the charge distribution reactivity. The NCI-RDG analysis clearly shows the strong interactions of the gas NO3 and NO2 with the endo-cavity environments of the CX[4] molecule. Finally, the non-covalent interactions analyses show that the calix[4]arene maybe useful for encapsulated the pollutant gas in the future. The sensitivity of the calix[4]arene molecule for these polluting gases opens a way to test the interaction of CX[4] with other types of biological molecules (proteins).



The authors acknowledge financial support from the Tunisian’s Ministry of high education and scientific research. In this work, we were granted access to the HPC resources of the FLMSN, ‘Fédération Lyonnaise de Modélisation et Sciences Numériques’, partner of EQUIPEX EQUIP@MESO and to the ‘Centre de calcul CCIN2P3’ at Villeurbanne, France.


  1. 1. Athar, M., Lone, M.Y., Jha, P.C.: Recognition of anions using urea and thiourea substituted calixarenes: A density functional theory study of non-covalent interactions. Chemical Physics. 501, 68-77 (2018). doi:10.1016/j.chemphys.2017.12.00
  2. 2. Kumagai, S., Hayashi, K., Kameda, T., Morohashi, N., Hattori, T., Yoshioka, T.: Identification of number and type of cations in water-soluble Cs+ and Na+ calix[4]arene-bis-crown-6 complexes by using ESI-TOF-MS. Chemosphere. 197, 181-184 (2018). doi:10.1016/j.chemosphere.2018.01.040
  3. 3. Arena, G., Contino, A., Gulino, F.G., Magrı̀, A., Sciotto, D., Ungaro, R.: Complexation of small neutral organic molecules by water soluble calix[4]arenes. Tetrahedron Letters. 41, 9327-9330 (2000). doi:10.1016/S0040-4039(00)01687-7
  4. 4. Ortolan, A., Oestroem, I., Caramori, G., Parreira, R., Muñoz-Castro, A., Bickelhaupt, F.M.: Anion Recognition by Organometallic Calixarenes: Analysis from Relativistic DFT Calculations. Organometallics. 37, (2018). doi:10.1021/acs.organomet.8b00292
  5. 5. Cabral, B.J.C., Coutinho, K., Canuto, S.: Dynamics of endo- vs. exo-complexation and electronic absorption of calix[4]arene-Ar2. Chemical Physics Letters. 612, 266-272 (2014). doi:10.1016/j.cplett.2014.08.036
  6. 6. Haino, T., Rudkevich, D.M., Shivanyuk, A., Rissanen, K., Jr, J.R.: Induced-Fit Molecular Recognition with Water-Soluble Cavitands. Chemistry – A European Journal. 6, 3797-3805 (2000).doi:10.1002/1521-3765(20001016)6:20<3797::AID-CHEM3797>3.0.CO;2-1
  7. 7. Wei, A.: Calixarene-encapsulated nanoparticles: self-assembly into functional nanomaterials. Chem Commun (Camb). 1581-1591 (2006). doi:10.1039/b515806k
  8. 8. Wang, J.-F., Huang, L.-Y., Bu, J.-H., Li, S.-Y., Qin, S., Xu, Y.-W., Liu, J.-M., Su, C.-Y.: A fluorescent calixarene-based dimeric capsule constructed via a M II –terpyridine interaction: cage structure, inclusion properties and drug release. RSC Advances. 8, 22530-22535 (2018). doi:10.1039/C8RA02146E
  9. 9. Atwood, J.L., Koutsantonis, G.A., Raston, C.L.: Purification of C60 and C70 by selective complexation with calixarenes. Nature. 368, 229-231 (1994). doi:10.1038/368229a0
  10. 10. Da Silva, E., Lazar, A.N., Coleman, A.W.: Biopharmaceutical applications of calixarenes. Journal of Drug Delivery Science and Technology. 14, 3-20 (2004). doi:10.1016/S1773-2247(04)50001-1
  11. 11. Fahmy, S.A., Ponte, F., Abd El-Rahman, M.K., Russo, N., Sicilia, E., Shoeib, T.: Investigation of the host-guest complexation between 4-sulfocalix[4]arene and nedaplatin for potential use in drug delivery. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy. 193, 528-536 (2018). doi:10.1016/j.saa.2017.12.070
  12. 12. Schühle, D.T., Peters, J.A., Schatz, J.: Metal binding calixarenes with potential biomimetic and biomedical applications. Coordination Chemistry Reviews. 255, 2727-2745 (2011)
  13. 13. Balasaheb Nimse, S., Kim, T.: Biological applications of functionalized calixarenes. Chemical Society Reviews. 42, 366-386 (2013). doi:10.1039/C2CS35233H
  14. 14. Hua, B., Shao, L., Zhang, Z., Sun, J., Yang, J.: Pillar[6]arene/acridine orange host–guest complexes as colorimetric and fluorescence sensors for choline compounds and further application in monitoring enzymatic reactions. Sensors and Actuators B: Chemical. 255, 1430-1435 (2018). doi:10.1016/j.snb.2017.08.141
  15. 15. B.J.C. Cabral, K. Coutinho, S. Canuto, Dynamics of endo- vs. exo-complexation and electronic absorption of calix[4]arene-Ar2, Chemical Physics Letters. 612 (2014) 266-272.
  16. 16. B. Gassoumi, H. Ghalla, R.Ben. Chaabane, DFT and TD-DFT investigation of calix[4]arene interactions with TFSI− ion, Heliyon. 5 (2019) e02822.
  17. 17. Y. Inokuchi, K. Soga, K. Hirai, M. Kida, F. Morishima, T. Ebata, Ultraviolet Photodissociation Spectroscopy of the Cold K+·Calix[4]arene Complex in the Gas Phase, J. Phys. Chem. A. 119 (2015) 8512-8518.
  18. 18. S. Kaneko, Y. Inokuchi, T. Ebata, E. Aprà, S.S. Xantheas, Laser spectroscopic and theoretical studies of encapsulation complexes of calix[4]arene, J Phys Chem A. 115 (2011) 10846-10853.
  19. 19. L.I. Shamova, G.A. Shamov, I.S. Antipin, A.I. Konovalov, Modeling K+ and Ag+ Complexation by Thiacalix[4]arene Amides Using DFT: The Role of Intramolecular Hydrogen Bonding, J. Phys. Chem. A. 113 (2009) 5691-5699.
  20. 20. G. Mazzone, M.E. Alberto, F. Ponte, N. Russo, M. Toscano, Anion-π weak interactions in a heteroaromatic calixarene receptor. A theoretical investigation, Inorganica Chimica Acta. 470 (2018) 379-384.
  21. 21. B. Gassoumi, M. Chaabene, H. Ghalla, R.B. Chaabane, Role of hydrogen bonding interactions within of the conformational preferences of calix[n = 4,6,8]arene: DFT and QTAIM analysis, J Mol Model. 26 (2019) 12.
  22. 22. G09 |,
  23. 23. Dennington RI, Keith T, Millam J, GaussView Version 5.0.8. Semichem Inc
  24. 24. S.F. Boys, F. Bernardi, The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors, Molecular Physics. 19 (1970) 553-566.
  25. 25. Johnson ER, Keinan S, Mori-Sánchez P, et al (2010) Revealing Noncovalent Interactions. J Am Chem Soc 132:6498-6506.

Written By

Bouzid Gassoumi, Fatma Ezzahra Ben Mohamed, Houcine Ghalla and Rafik Ben Chaabane

Submitted: June 24th, 2020 Reviewed: September 1st, 2020 Published: November 25th, 2020