Next Article in Journal
Marker Substances in the Aroma of Truffles
Previous Article in Journal
Date Fruit (Phoenix dactylifera L.) Cultivar Extracts: Nanoparticle Synthesis, Antimicrobial and Antioxidant Activities
Previous Article in Special Issue
Interactions of Ag Particles Stabilized by 7,7,8,8-Tetracyanoquinodimethane with Olefin Molecules in Poly(ether-block-amide)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mesoporous Silica Modified with 2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde as an Effective Adsorbent for Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study

1
Laboratory of Organic Synthesis, Extraction and Valorization, Faculty of Sciences Ain Chock, Hassan II University, Casablanca 20000, Morocco
2
Laboratory of applied Chemistry and Environment (LCAE), Faculty of Sciences, Mohamed Premier University, Oujda 60000, Morocco
3
Department of Chemistry, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
4
Biochemistry Department, College of Science, King Saud University, P.O. Box 22452, Riyadh 11495, Saudi Arabia
5
Laboratory of Plant Biotechnology Applied to Crop Improvement, Faculty of Science of Sfax, University of Sfax, Sfax 3038, Tunisia
6
Department of Chemistry, Vidya Bharati Mahavidyalaya, Amravati 444 602, India
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(16), 5168; https://doi.org/10.3390/molecules27165168
Submission received: 18 July 2022 / Revised: 11 August 2022 / Accepted: 11 August 2022 / Published: 13 August 2022
(This article belongs to the Special Issue Advances in Organic–Inorganic Hybrid Materials)

Abstract

:
In this study, we will present an efficient and selective adsorbent for the removal of Cu(II) ions from aqueous solutions. The silica-based adsorbent is functionalized by 2-phenylimidazo[1,2-a] pyridine-3-carbaldehyde (SiN-imd-py) and the characterization was carried out by applying various techniques including FT-IR, SEM, TGA and elemental analysis. The SiN-imd-py adsorbent shows a good selectivity and high adsorption capacity towards Cu(II) and reached 100 mg/g at pH = 6 and T = 25 °C. This adsorption capacity is important compared to other similar adsorbents which are currently published. The adsorption mechanism, thermodynamics, reusability and the effect of different experimental conditions, such as contact time, pH and temperature, on the adsorption process, were also investigated. In addition, a theoretical study was carried out to understand the adsorption mechanism and the active sites of the adsorbent, as well as the stability of the complex formed and the nature of the bonds.

1. Introduction

Environmental pollution has become a major concern for mankind and the eco-system. It has impacted human health and it has led to economic and social issues for millions of people worldwide. Furthermore, water pollution is considered one of the most serious problems that requires an immediate and radical solution. It could also lead to serious threats to sustainable environmental development [1]. Rapid industrial development in many countries and human activities result in the contamination of water with heavy metal ions like Cu, Pb, Hg, etc., and this problem has caused several diseases [2]. Cu(II) is toxic at high concentrations and responsible for numerous illnesses due to its bio-accumulable and non-biodegradable nature [3]. Several techniques such as adsorption [4], liquid–liquid extraction [5], membrane filtration [6], ion exchange [7] and electrolysis [8] have been used for the purification of wastewater containing heavy metal ions. Recently, adsorption has emerged as a more promising technique due to: (i) ease of operation, (ii) higher enrichment factor, (iii) reduction of organic solvent usage, low cost and extraction time, (iv) high selectivity, (v) ease of separation and (vi) ability to combine with different modern detection techniques [9]. Various types of materials such as silica [10], activated carbon [11], biopolymers [12], zeolite [13] and clay [14] are successfully employed as effective adsorbents. Compared with other substrates, silica gel has garnered attention in sorption of heavy metals due to its large surface area and high thermal and mechanical stability [15]. Much of the research in recent times has reported the immobilization of organic ligands onto mesoporous silica to decontaminate heavy metal from wastewater [16]. The efficiency of these adsorbents depends mainly on the affinity of the donor atoms (S, O and N) that are deposited on the surface of the materials. Several ligands, bearing donor atoms, have been synthesized and immobilized on the surface of silica for water pollution control [17,18,19,20]. In this context, imidazo[1,2-a] pyridine has attracted attention because it plays a major role in coordination chemistry [21,22]. On the other hand, the condensation of an amine of silica modified with an aldehyde group of imidazo[1,2-a] pyridine leads to the formation of a Schiff base, which is known to have high bonding affinities and to form very stable complexes with metal ions [23,24].
In the present work, imidazo[1,2-a] pyridine-3-carbaldehyde has been used to functionalize silica for the first time and to afford a functional adsorbent. In fact, its structure was fully characterized. The new synthetic material shows a high selectivity towards Cu2+ ions, it also shows a high adsorption capacity compared to other similar adsorbents [9]. The influence of pH value, contact time, metal concentration and temperature for the adsorption of Cu(II) were studied using the batch method. In addition, the selectivity and regeneration capacity were discussed, and the adsorption mechanism was revealed on the basis of theoretical calculations.

2. Experimental Section

2.1. Materials and Methods

All reagents (Aldrich, purity 99.5%) were of analytical grade. Initially, silica gel (60 Å, 70−230 mesh) (E. Merck) was activated at 120 °C for 24 h. The quantification of metal ions in aqueous solutions was determined by atomic absorption (Spectra Varian A.A. 400 spectrophotometer, Shelton, CT, USA). The surfaces were characterized by a CHN analyzer (Microanalysis Center Service, CNRS, Paris, France), Fourier-transform infrared spectroscopy (FTIR, Perkin Elmer System 2000, Waltham, MA, USA) at 25 °C, SEM (FEI-Quanta 200, Hillsboro, OR, USA), TG/DTA (Perkin Elmer Diamond, Waltham, MA, USA) under a 90:10 oxygen/nitrogen atmosphere at 10 °C·min−1, 13C solid state nuclear magnetic resonance (NMR, CP MAX CXP 300 MHz, Lille, France) and BET (ThermoQuest Sorpsomatic 1990 analyzer, Lille, France). Nitrogen adsorption–desorption isotherm plots were obtained on a Thermoquest Sorpsomatic 1990 analyzer after the materials had been purged in a stream of dry nitrogen. The pH determinations were carried out with a pH 2006, J. P. Selecta s. a. pH meter.

2.2. Synthesis of 2-phenylimidazo[1,2-a] pyridine-3-carbaldehyde

2-Phenylimidazo[1,2-a] pyridine (3 mmol) was suspended in CHCl3, the obtained solution was added to a mixture of POCl3 (11 mL) and DMF (8 mL) at 0–5 °C, under stirring for 30 min, followed by refluxing for 6 h and then neutralized using Na2CO3. The residue was extracted by dichloromethane and purified on silica gel to give a white solid.
13C NMR (100 MHz, CDCl3) δ (ppm): 115.10 (C6), 117.22 (C8), 120.62 (C3), 128.36 (C7), 128.92 (2Cph), 129.45 (3Cph), 130.18 (C5), 132.32 (Cph), 147.51 (C8a/C2), 157.90 (C2/C8a), 179.26 (CHO)

2.3. Preparation of Amine-Functionalized Silica SiNH2

A mixture of activated silica (10 g) and dry toluene (50 mL) was stirred under reflux for 2 h. After that, 3-aminopropyltrimethoxysilane (10 mL) was added dropwise to the suspended solution and the refluxing was maintained for 24 h. The SiNH2 formed was filtered and washed by Soxhlet extraction with different organic solvents for 12 h. Finally, the obtained modified silica was dried at 60 °C for 24 h.

2.4. Fabrication of the SiN-imd-py Adsorbent

To prepare the (SiN-imd-py), a mixture of 2-phenylimidazo[1,2-a]pyridine-3-carbaldehyde (0.5 g) and 3- aminopropylsilica (SiNH2) (3.5 g) in 30 mL of methanol was mixed, heated at reflux (140 °C) and stirred for 24 h. At room temperature, the solid residue was filtered. Methanol, dichloromethane, THF and diethyl ether were used in Soxhlet extraction of the product for 10 h. The desired solid product was dried completely (Scheme 1).

2.5. Batch Adsorption Experiment

The parameters influencing the adsorption of metal ions via the material, such as the effect of concentration (10 to 300 mg L1), contact time (5 to 35 min), pH (1 to 7) and temperature (25 to 45 °C), were studied by the batch method. For this purpose, 10 mL of the ionic solution containing the following salts: (Pb(NO3)2,6H2O); (Cd(NO3)2, 6H2O); (Cu(NO3)2,3H2O) and (Zn(NO3)2,6H2O), together with 10 mg of hybrid material, is introduced into a set of test tubes. The pH values are adjusted with diluted HCl and NaOH solutions, and the mixture is stirred at 25 °C. The adsorption capacities are determined as follows [25]:
qe = (C0 − Ce) × V/m
where qe (mg g1) is the adsorption capacity; C0 and Ce (mg L−1) are, respectively, the initial and equilibrium concentrations; m(g) is the mass of the adsorbent; and V (L) is the volume of the solution.

2.6. Computational Methods

Density functional theory (DFT) [26] calculations and QTAIM were used to illustrate the adsorption mechanism [27]. The Gaussian 09 [28] software and gaussview05 [29] were used to optimize the structure. These software and methods were used to investigate the interaction between the material and metal ion Cu2+ [30]. All the structures reported herein were fully optimized in gas phase at the B3LYP/GenECP level of theory [31,32]. The LanL2DZ relativistic pseudo potential [33] basis set was used for the metals, and the 6–311++G (d, p) [34,35] basis set was used for all other atoms. The interaction energy between a single SiN-imd-py hydrogel unit and metal ions (Cu2+) was calculated by:
ΔIE = E(complex) (Ematerial + Eion)
These findings were supported by one of the most well-known DFT indices [36], the nucleophilic Parr functions (P-) [37], which is a useful tool for assessing the potential nucleophilic compound sites. By adding an electron of the Mulliken atomic spin density at the radical cation, we obtain these functions. Bader proposed the technique of topology analysis for analyzing electron density in “atoms in molecules” (AIM) theory, which is also known as “quantum theory of atoms in molecules” (QTAIM) [38]. The electron localization function ELF [39] is used for finding the electron density. Because ELF is a density-based function that can be understood in terms of the relative local excess of kinetic energy density, it has been recognized as a technique that provides a good framework for the investigation of changes in pair electron density [40]. The localized-orbital locator (LOL) [41] is used to further understand the physical characteristics of the chemical bond between two atoms. Schmider and Becke [42] have recently presented a crucial tool, the localized-orbital locator, to explain bonding properties in terms of the local kinetic energy. The latter is roughly defined as a function of electron density and its first and second derivatives.

3. Results and Discussion

3.1. Characterization

3.1.1. Elemental Analysis

Elemental analysis was used to assess the elemental composition of the materials and to confirm the presence or absence of an organic motif on the modified silica SiNH2. (Table 1)
Comparing w.t% of the two materials SiNH2 and Si-imd-py, we notice an increase in the percentages of the elements of C, N and H, which indicates that the functionalization was successfully realized.

3.1.2. FTIR

FTIR characterization technique allowed us to identify the different types of vibration of the functional groups presented in the different materials studied. FTIR spectra of SiG, SiNH2 and SiN-imd-py materials are shown in Figure 1. Indeed, the spectrum of the free silica SiG shows a large band at 3446 cm1 and a weak band at 1652 cm1 attributed to O-H elongation vibrations. The bands at 1086 and 798 cm1 are attributed to Si-O-Si and Si-O vibrations, respectively [43]. For the SiNH2 spectrum, two new adsorption bands were observed around 2932 and 1582 cm1, which are characteristic of C-H and NH2 elongation vibration resettlement [44]. These results indicate that the spacer arm has been immobilized on the silica surface. For the spectrum of the SiN-imd-py material, we also note the appearance of a new band at 1440 cm1 that is attributable to the C=N vibration, which confirms the grafting of the organic motif onto the SiNH2 support.

3.1.3. Scanning Electron Microscope (SEM)

Scanning electron microscopy (SEM) allows the observation of the morphology of material and the distribution of grains on the material surface. The micrographs obtained are shown in Figure 2. As we can see from the SEM images, the morphology of the non-functionalized and functionalized silica particles is different. An irregular and smooth surface was observed for SiG, while the morphology of SiNH2 and SiN-imd-py materials is rough and porous, and agglomerations of molecules appear on the surface of the materials, indicating the presence of organic molecules on the silica surface.

3.1.4. Thermogravimetric Analysis (TGA)

Thermogravimetry makes it possible to follow, as a function of temperature, the evolution of the loss of mass of the sample, mainly caused by the vaporization of water and by the destruction of organic matter. Figure 3 shows the thermogravimetric curves of SiN- imd-py material and both materials, SiNH2 and SiG. For virgin silica, two stages of mass loss can be distinguished. The first stage that corresponds to a weight loss is at 3.15% (between 25 °C and 110 °C). It was assigned to the desorption of water adsorbed on the surface of silica gel. The second stage of mass loss is at 5.85% observed between 110 °C and 800 °C and is attributed to the condensation of silanol groups leading to siloxane groups (Si-O-Si). For the SiNH2, it presented a high mass loss at 6.65% between 200–800 °C. This mass loss is related to the decomposition of the organic chain immobilized on silica. The SiN-imd-py material is also presented in two decomposition stages. The first loss stage, between 25 and 110 °C, is at 4.48%. This is attributable to the removal of adsorbed water. The second mass loss between 200 to 800 °C is 16.68%, due to the degradation of the organic fragment grafted on the silica gel surface.

3.1.5. N2 Physisorption Studies

To investigate the changes in surface area and porosity of SiN-imd-py material, we measured the specific surface area (SBET) according to the Brunauere–Emmett–Teller method [45] and the pore volumes were measured according to the Barret–Joyner–Halenda (BJH) method [46]. Figure 4 shows that the material presents a type IV isotherm, which clearly indicates that the material is mesoporous according to the classification of the International Union of Pure and Applied Chemistry. The figure also reveals the H2-type hysteresis loop for partial pressures P/P0 ≥ 0.4. The results of the physical properties found are reported in Table 1.
From the analysis of the data in Table 2, we observe that the SiN-imd-py material is characterized by a decrease in pore volume and specific surface area compared to the SiNH2 precursor. This decrease is due to the grafting of organic ligands on the silica surface which partially blocks the adsorption of nitrogen molecules.

3.2. Adsorption Studies

3.2.1. Effect of pH

The pH of the solution is one of the most important parameters controlling the adsorption of metal ions. In this study, the effect of various pH values (1–7) at 25 °C was investigated by means of fabricated functionalized material. As shown in Figure 5, when the pH is below 5, the adsorption capacity increases with the pH; when the pH > 5, the adsorption capacity of the adsorbent tends to be stable. The SiN-Imd-Py adsorbent has a higher adsorption capacity for Cu(II) over a wider pH range. Therefore, modification of the imidazopyridine on the SiNH2 surface can significantly improve the adsorption capacity. When the pH value was at 1–2, the adsorbent SiN-imd-py showed a negligible amount of adsorbed Cu(II). At this moment, the grafted imidazopyridine–Schiff base ligand gradually deprotonated, and this implies that the positive charge increases, the electrostatic repulsion with the metal cation increases, and the adsorption performance decreases. The adsorption of Cu(II) on the adsorbent is difficult to quantify at pH above 7 because it can be masked by the precipitation phenomenon in the Cu(OH)2 form. The adsorption capacity of Cu(II) ion reached its maximum at a pH between 6 and 7. Therefore, a pH of 7 was chosen for the extraction of the Cu(II) studied in all subsequent studies.

3.2.2. Effect of Contact Time and Adsorption Mechanism

The effect of time is of considerable practical interest in adsorption and one of the most important characteristics defining the effectiveness of an adsorption. The adsorption kinetics for Cu(II) onto the SiN-imd-py are shown in Figure 6. It can be clearly seen that the adsorption of Cu(II) is very rapid during the first five minutes. This could be attributed to the large number of ligands available for chelation and a high level of free Cu(II) ions at this stage. After 30 min, the sorption rate decreased and finally reached a plateau, thus confirming the rapid adsorption kinetics for Cu(II). Equilibrium would have been reached due to saturation of the adsorption sites. Consequently, the optimal contact time was confirmed at 30 min for the following experiments.
To better comprehend the adsorption procedure, it is necessary to study the adsorption kinetics. This procedure is evaluated by two models, pseudo-first-order and pseudo-second-order, which are generally interpreted by the following equations [47,48]:
Pseudo-first-order model:
ln(qe − qt) = lnqe − k1t
Pseudo-second-order model:
t/qt = 1/K2qe2 + t/qe
where qe (mg g1) and qt (mg g1) are the amounts of adsorbate at equilibrium and at time t (s), respectively. K1(min1) and K2 (g (mg min1)1) are the rate constants of the first-order and the second-order models, respectively. The corresponding parameters for the two kinetic models are listed in Table 2. The theoretical qe value of the pseudo-second-order model fits the experimental values better, indicating that the model is proposed to describe the adsorption kinetics of Cu(II) by the SiN-imd-py adsorbent (Figure 7). The correlation coefficient (R²) of the pseudo-second-order kinetic model is near to 1, much higher than the correlation coefficient of the pseudo-first-order model. The sufficiency of the experimental data with this model indicates that Cu(II) adsorption is controlled by the chemical procedure (chemisorption) due to the existence of chemical interactions between the metal and the ligand grafted on the silica (coordination bonds).

3.2.3. Influence of Initial Concentration

In the adsorption process, the effect of concentration plays an important role because it provides necessary information about adsorbing performance in different concentrations. Accordingly, different initial concentrations of Cu(II) varying from 10 to 300 mg L−1 were employed to investigate the adsorbents further using the batch method. As shown in Figure 8, the amount of Cu(II) adsorbed by the SiN-imd-py gradually increases with increasing initial metal ion concentration. This suggests that the available adsorption sites on the adsorbent surface were handily occupied by the metal ions and that the adsorption capacity was very high. As the initial metal concentration increases, a saturation value is reached. This observation can be explained by the fact that there are fewer active sites on the surface of the adsorbent that are available for further adsorption of metal ions. This result suggests that the chelation of the Imidazopyridine-Schiff base heightens the adsorption capacity of the adsorbent.

3.2.4. Adsorption Isotherms

The equilibrium relationships between the concentration of the adsorbate and the amount of adsorbate accumulated on the adsorbent were explained by the adsorption isotherms. Indeed, the adsorption capacity results were fitted by four classical isotherm models, namely Langmuir, Freundlich, Dubinin-Radushkevich (D-R) and Temkin, which are generally described as the following equations [49,50]:
Langmuir model:
Ce/qe = Ce/qm + 1/qKL
where Ce (mg L1) is the equilibrium concentration, qe and qm (mg g−1) denote the equilibrium and theoretical maximum adsorption of the adsorbent, KL (L. mol1), is the Langmuir affinity.
Freundlich model:
qe = KF Ce1/n
where KF is the Freundlich constant and n is energy or intensity of adsorption.
Dubinin–Radushkevich (D-R) model:
ln (qe) = ln (qm) − βε2
where β (mol J1)2 is a constant related to the mean free energy E of adsorption (E = (2β) − 0.5), and ε (J mol1) is the Polanyi potential related to the equilibrium concentration (ε = RT ln (1 + 1/Ce)).
Temkin model:
qe = RT/bt lnAt + RT/bt ln Ce
where bt (J mol1) and At (L mg1) are Temkin isotherm constant and T defines the Kelvin temperature (K), and R (8.314 J mol1 K1) defines the universal gas constant.
The fitting plots of the Langmuir model are presented in Figure 9, and the theoretical parameters of four adsorption isotherms, as well as the regression coefficients, are listed in Table 3. From the values in Table 4, the Langmuir model provides good correlation coefficient (>0.990), suggesting that the Langmuir model was appropriate to describe the adsorption process of Cu(II) onto SiN-imd-py adsorbent. On the other hand, the calculated equilibrium adsorption capacity q (mg g1) based on the Langmuir model is much closer to experimental data qe (mg g1), which also proved the adequacy of the Langmuir model. The above fact demonstrated that the adsorption process of Cu(II) by the adsorbent was attributed to a homogeneous monolayer adsorption and the maximum adsorption capacity was examined at 103.51 mg g1.

3.2.5. Thermodynamic Studies

The thermodynamic behavior was studied to determine the thermodynamic parameters of Cu(II) adsorption onto SiN-imd-py, including the change in standard enthalpy (ΔH°), standard entropy (ΔS°) and Gibbs energy (ΔG°), which were determined using the following equations [51]:
Kd = (C0 − Ce)/Ce
Ln Kd = ΔS°/R − ΔH°/RT
ΔG° = ΔH° − TΔS°
where Kd implies the adsorption distribution coefficient. R (8.314 J mol1 K) and T(K) are the universal gas constant and the absolute temperature, respectively. The curve Ln(Kd) versus 1/T (Figure 10) was employed to determine the thermodynamic parameters ΔG°, ΔH° and ΔS° which are given in Table 5. The positive value of ΔH° provides clear evidence that the adsorption phenomenon was endothermic. The endothermic nature also involves the absorption of energy in the form of heat from the environment during the adsorption process. The positive values of ΔS° indicate an increase in randomness at the solid–liquid interface during the adsorption of Cu(II) ions, mainly due to the dehydration of metal ions towards the adsorption sites. The negative ΔG° values at all temperatures indicate that the adsorption phenomenon is spontaneous and more favorable at high temperatures. The metal ion adsorption process suggests that a large amount of heat is consumed to transfer the metal ions from the solution to the solid phase. The thermodynamic results show that the elimination of Cu(II) on the adsorbent was endothermic and spontaneous in nature.
It is clear from this table that all the values of ΔG° are negative. The results presented above suggest that the adsorption process involved is spontaneous. On the other hand, the change of the standard free energy decreases with increasing temperatures regardless of the nature of adsorbent. This indicates that a better adsorption is actually obtained at higher temperatures.

3.2.6. Adsorption Selectivity for Cu(II)

The competitive adsorption of heavy metal ions by SiN-imd-py material was studied in a mixed quaternary system (Cu(II), Zn(II), Cd(II) and Pb(II)) by the batch method (Figure 11). The material has a higher selectivity towards Cu(II) compared to other ions. The high selectivity of SiN-imd-py towards Cu(II) is due to the types of ligand immobilized on the silica surface, reflecting the extraordinary ability and efficiency to form more stable complexes with Cu(II).

3.2.7. Desorption and Recycling

Desorption is one of the most important indicators to evaluate whether an adsorbent is practical or not. Adsorbent desorption was examined by adding 2 mol/L HCl to 10 mg of material-Cu(II) with stirring at room temperature for 60 min. Then, the material was filtered and neutralized by a NaOH solution. After washing, the material was dried with the help of a vacuum desiccator and then in the oven for the next adsorption. Table 6 represents the extraction efficiency of the material after five cycles of Cu(II) adsorption–desorption. It is interesting to note that the adsorbents retained more than 95% of their adsorption capacity. Thus, the material has an exceptional recyclability and applicability that could be useful for the purification of copper-contaminated wastewater.

3.2.8. Comparison with Similar Adsorbents

Table 7 shows a comparative study highlighting the interest and efficiency of our adsorbent towards Cu(II) compared to others recently described in the literature.

3.3. Adsorption Mechanism

3.3.1. Active Sites Study

The nucleophilic indexes of Parr (P-) allow us to determine the atomic sites of the donor species. For the material SiN-imd-py, one nitrogen N2 (Figure 12) atom is not considered to be a reactive site because the value of (P-) is negative −0.064 (Table 8). This indicates that it is not involved in coordination with Cu. Nevertheless, the nitrogen N9 and N11 display positives value of p-, revealing that these are suitable sites for an electrophilic attack via Cu2+ ions.
The QTAIM approach is necessary to differentiate between the two nucleophilic sites N9 and N11; as shown in Figure 12 (right), the N9 site (the red arrow) is favored to attack the copper ion Cu2+. By NCI we could confirm the attack of N9 because from the left Figure 12 (the red arrow), no interaction was observed in this region, which further favors the attack of N9 on this side of the molecule [36].

3.3.2. Cu(II) Complexation Study

The optimized geometries of possible complexes are proposed on the basis of NBO analysis, p- indices and QTAIM. In this sense, E was calculated, which is the interaction energy; this was −229,208 for complex1 and −229,199 for complex2 which means that complex1 is more stable in relation to complex2 [61].
To obtain information on the binding nature of the studied complexes we introduced the analysis of the natural binding orbital NBO [62]. The interaction between the N9, O39, O42 and O45 atoms was evaluated by the second order energy E(2) and the electronic configuration (EC) (Table 9). The results of the NBO analysis presented in the table reveal that the four atoms N9, O39, O42 and O45 have almost the same EC, but for Cu2+ this is not the case (EC(Cu2+(1)) < EC(Cu2+(2)).
The secondary energy E(2) values are also listed in the same table, and we see that all energy values found for complex1 are higher than those found for complex2. This shows that the charge transfer is from N9, O39, O42, and O45 to the empty orbitals of Cu2+, and also indicates that complex1cis has undergone an important electron transfer.
Regarding the bond lengths, they are important to determine the type of bond between the atoms as well as the possibility of existence of a bond, because if the bond length is large enough, an adequate interaction between the two atoms cannot take place. The results displayed in Table 10 show that the distance between Cu and other atoms is small (approximately 1.983 Å to 2.249 Å). This fact shows that the interaction between Cu and other atoms is possible. Lower values of the same distance for complex1 show that it is favorable as a stable complex. QTAIM was the most innovative work of the Bader Group, and an interatomic contour line with a (3,–1) cp (Figure 13) implies that electron density is accumulated among the nuclei that are reassembled in this way. The corresponding values of electron density and lagrangian for complex1 (0.832, 0.510) are large compared to complex2 (0.808, 0.129), while the results of the calculation of the laplacian for complex1cis 0.135 are small compared to complex2 0.495; this means that the bonds formed for complex1 are strong compared to complex2. Consequently, the stability of complex1 is important [62].
The ELF localized electron function is recognized as a technique offering a suitable framework for the study of electron density changes. Similarly, Schmider and Becke proposed the localized-orbital locator to define the characteristics of bonding in terms of local kinetic energy. The results (Table 10) show that the ELF and LOL values for complex1cis (respectively 0.102 and 0.252) are higher than the values for complex2 (0.100, 0.250) which reflects the higher stability of complex1.
To highlight the properties of a line, we have chosen to draw this curve (Figure 14), where the dotted line corresponds to the position Y=0 and the red curve represents the position of the two selected atoms Cu and N. For the complex1trans (left), the two values (1.345 and 5.417) indicate the two atoms that we have chosen; also for complex1 (right) the atoms are located at 1.350 and 5.399 according to the Bohr radius.
It is noticeable that there is a small difference between the two curves. For complex1, the peak at position 1.345 rose to 0.009, but the peak at 1.350 of the other complex1 changed to 0.016. This explicitly explains the difference in stability between the two complexes. In Figure 15, the contour maps of electron density on the plane (N9-Cu-O39) for the two complexes show that there are simple bonds between these three atoms. The same remark can be made for the plan maps of the localized-orbital locator in Figure 15 on the right-hand side.

4. Conclusions

In conclusion, we have developed a new hybrid material based on silica and heterocyclic ligand. This material was prepared by immobilization of 2-phenylimidazo[1,2-a] pyridine-3-carbaldehyde on the silica surface. In addition, the synthesized material was characterized by the usual physicochemical techniques (AE, IR-TF, SEM, BET and ATG). The SiN-imd-py material showed a marked adsorption efficiency towards Cu(II) from wastewater (97.11 mg g−1). We also identified the different parameters influencing the adsorption such as pH, contact time, initial concentration and temperature. Imidazopyridine-Schiff base ligand is very sensitive to pH changes as the optimal pH for adsorption is at pH = 6. The equilibrium adsorption time was about 30 min. The kinetic study shows that the adsorption follows the pseudo-second-order model, indicating that the adsorption is chemical. The study of adsorption isotherms is described by the Langmuir model, which shows that the adsorption is monolayer. Actually, the thermodynamic parameters indicate that the metal adsorption process is spontaneous and endothermic. The NCI method and the Parr p-indices help to determine the nucleophilic sites of the material, as in our work. Thanks to these methods, we could find the active sites of SiN-imd-py. To highlight the stability of the complexes and the nature of the bonds, we used the following methods: DFT, NBO, QTAIM, ELF and LOL. Clearly, by means of these methods, we found that complex1 is more stable. It is interesting to note that the adsorbent shows selectivity towards Cu(II). This material is easily regenerated by acid washing. Finally, the adsorbent is very promising for the purification and extraction of Cu(II) from wastewater.

Author Contributions

R.S. and S.T.: conceptualization, supervision, project administration, methodology, resources, data curation, writing—original draft, review and editing; I.H.: elaboration of theoretical aspects, data curation, writing—original draft, review and editing; O.R.: conduct of experiment; V.H.M.: review and editing, methodology; Z.M.A., M.I.A.-Z., A.B.B. and S.R.: project administration, methodology, writing—original draft and paying publication fees. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Saud University for funding this work through research group No (RGP-070).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not available.

References

  1. Töre, G.Y.; Özkoç, B. Recent developments in aquatic macrophytes for environmental pollution control: A case study on heavy metal removal from lake water and agricultural return wastewater with the use of duckweed (Lemnacea). In Phytoremediation Technology for the Removal of Heavy Metals and Other Contaminants from Soil and Water, 1st ed.; Elsevier: Amsterdam, The Netherlands, 2022; pp. 75–127. [Google Scholar] [CrossRef]
  2. Guo, L.-C.; Lv, Z.; Ma, W.; Xiao, J.; Lin, H.; He, G.; Li, X.; Zeng, W.; Hu, J.; Zhou, Y.; et al. Contribution of heavy metals in PM2.5 to cardiovascular disease mortality risk, a case study in Guangzhou, China. Chemosphere 2022, 297, 134102. [Google Scholar] [CrossRef] [PubMed]
  3. Mitra, S.; Chakraborty, A.J.; Tareq, A.M.; Bin Emran, T.; Nainu, F.; Khusro, A.; Idris, A.M.; Khandaker, M.U.; Osman, H.; Alhumaydhi, F.A.; et al. Impact of heavy metals on the environment and human health: Novel therapeutic insights to counter the toxicity. J. King Saud Univ. Sci. 2022, 34, 101865. [Google Scholar] [CrossRef]
  4. Rajendran, S.; Priya, A.; Kumar, P.S.; Hoang, T.K.; Sekar, K.; Chong, K.Y.; Khoo, K.S.; Ng, H.S.; Show, P.L. A critical and recent developments on adsorption technique for removal of heavy metals from wastewater-A review. Chemosphere 2022, 303, 135146. [Google Scholar] [CrossRef] [PubMed]
  5. Sorouraddin, S.M.; Farajzadeh, M.A.; Okhravi, T. Cyclohexylamine as extraction solvent and chelating agent in extraction and preconcentration of some heavy metals in aqueous samples based on heat-induced homogeneous liquid-liquid extraction. Talanta 2017, 175, 359–365. [Google Scholar] [CrossRef]
  6. Fadila, C.; Othman, M.; Adam, M.; Takagi, R.; Yoshioka, T.; Khongnakorn, W.; Rahman, M.; Jaafar, J.; Ismail, A. Adsorptive membrane for heavy metal removal: Material, fabrication, and performance. Mater. Today Proc. 2022; in press. [Google Scholar] [CrossRef]
  7. Pan, S.; Shen, J.; Deng, Z.; Zhang, X.; Pan, B. Metastable nano-zirconium phosphate inside gel-type ion exchanger for enhanced removal of heavy metals. J. Hazard. Mater. 2022, 423, 127158. [Google Scholar] [CrossRef]
  8. Zhou, Q.; Liao, L.; Zhou, H.; Li, D.; Tang, D.; Yu, F. Innovative strategies in design of transition metal-based catalysts for large-current-density alkaline water/seawater electrolysis. Mater. Today Phys. 2022, 26, 100727. [Google Scholar] [CrossRef]
  9. Tighadouini, S.; Radi, S.; Roby, O.; Hammoudan, I.; Saddik, R.; Garcia, Y.; Almarhoon, Z.M.; Mabkhot, Y.N. Kinetics, thermodynamics, equilibrium, surface modelling, and atomic absorption analysis of selective Cu(ii) removal from aqueous solutions and rivers water using silica-2-(pyridin-2-ylmethoxy)ethan-1-ol hybrid material. RSC Adv. 2022, 12, 611–625. [Google Scholar] [CrossRef]
  10. Wieszczycka, K.; Wojciechowska, I.; Filipowiak, K.; Buchwald, T.; Nowicki, M.; Dudzinska, P.; Strzemiecka, B.; Voelkel, A. Novel iminepyridinium -modified silicas as super-adsorbents for metals ions. Appl. Surf. Sci. 2022, 596, 153555. [Google Scholar] [CrossRef]
  11. Mandal, S.; Calderon, J.; Marpu, S.B.; Omary, M.A.; Shi, S.Q. Mesoporous activated carbon as a green adsorbent for the removal of heavy metals and Congo red: Characterization, adsorption kinetics, and isotherm studies. J. Contam. Hydrol. 2021, 243, 103869. [Google Scholar] [CrossRef]
  12. Zhao, C.; Liu, G.; Tan, Q.; Gao, M.; Chen, G.; Huang, X.; Xu, X.; Li, L.; Wang, J.; Zhang, Y.; et al. Polysaccharide-based biopolymer hydrogels for heavy metal detection and adsorption. J. Adv. Res. 2022; in press. [Google Scholar] [CrossRef]
  13. Lv, Y.; Ma, B.; Liu, Y.; Wang, C.; Chen, Y. Adsorption behavior and mechanism of mixed heavy metal ions by zeolite adsorbent prepared from lithium leach residue. Microporous Mesoporous Mater. 2022, 329, 111553. [Google Scholar] [CrossRef]
  14. Es-Sahbany, H.; Hsissou, R.; El Hachimi, M.; Allaoui, M.; Nkhili, S.; Elyoubi, M. Investigation of the adsorption of heavy metals (Cu, Co, Ni and Pb) in treatment synthetic wastewater using natural clay as a potential adsorbent (Sale-Morocco). Mater. Today Proc. 2021, 45, 7290–7298. [Google Scholar] [CrossRef]
  15. Baskaran, K.; Ali, M.; Gingrich, K.; Porter, D.L.; Chong, S.; Riley, B.J.; Peak, C.W.; Naleway, S.E.; Zharov, I.; Carlson, K. Sol-gel derived silica: A review of polymer-tailored properties for energy and environmental applications. Microporous Mesoporous Mater. 2022, 336, 111874. [Google Scholar] [CrossRef]
  16. Tighadouini, S.; Radi, S.; Garcia, Y. Selective chemical adsorption of Cd(ii) on silica covalently decorated with a β-ketoenol-thiophene-furan receptor. Mol. Syst. Des. Eng. 2020, 5, 1037–1047. [Google Scholar] [CrossRef]
  17. Yanovska, E.; Savchenko, I.; Sternik, D.; Kychkiruk, O.; Ol’Khovik, L.; Buriachenko, I. In Situ Immobilization on the Silica Gel Surface and Adsorption Capacity of Poly[N-(4-carboxyphenyl)methacrylamide] on Toxic Metal Ions. Nanoscale Res. Lett. 2017, 12, 313. [Google Scholar] [CrossRef]
  18. Tighadouini, S.; Radi, S.; Ferbinteanu, M.; Garcia, Y. Highly Selective Removal of Pb(II) by a Pyridylpyrazole-β-ketoenol Receptor Covalently Bonded onto the Silica Surface. ACS Omega 2019, 4, 3954–3964. [Google Scholar] [CrossRef]
  19. Rigoletto, M.; Calza, P.; Gaggero, E.; Laurenti, E. Hybrid materials for the removal of emerging pollutants in water: Classification, synthesis, and properties. Chem. Eng. J. Adv. 2022, 10, 100252. [Google Scholar] [CrossRef]
  20. Bilgiç, A.; Çimen, A. Removal of chromium(VI) from polluted wastewater by chemical modification of silica gel with 4-acetyl-3-hydroxyaniline. RSC Adv. 2019, 9, 37403–37414. [Google Scholar] [CrossRef]
  21. Li, B.; Ni, L.-S.; Yong, G.-P. Three new coordination compounds based on a new 3-position substituted imidazo[1,2-a]pyridine ligand: Syntheses, crystal structures and photoluminescent properties. Polyhedron 2018, 154, 21–26. [Google Scholar] [CrossRef]
  22. Kanthecha, D.A.; Bhatt, B.S.; Patel, M.N. Synthesis, characterization and biological activities of imidazo[1,2-a]pyridine based gold(III) metal complexes. Heliyon 2019, 5, e01968. [Google Scholar] [CrossRef]
  23. Rakhtshah, J.; Yaghoobi, F. Catalytic application of new manganese Schiff-base complex immobilized on chitosan-coated magnetic nanoparticles for one-pot synthesis of 3-iminoaryl-imidazo[1,2-a]pyridines. Int. J. Biol. Macromol. 2019, 139, 904–916. [Google Scholar] [CrossRef]
  24. Wang, W.; Wu, G.; Zhu, T.; Yang, Y.; Zhang, Y. Synthesis of -thiazole Schiff base modified SBA-15 mesoporous silica for selective Pb(II) adsorption. J. Taiwan Inst. Chem. Eng. 2021, 125, 349–359. [Google Scholar] [CrossRef]
  25. Xue, X.; Li, F. Removal of Cu(II) from aqueous solution by adsorption onto functionalized SBA-16 mesoporous silica. Microporous Mesoporous Mater. 2008, 116, 116–122. [Google Scholar] [CrossRef]
  26. Rivero, P.; García-Suárez, V.M.; Pereñiguez, D.; Utt, K.; Yang, Y.; Bellaiche, L.; Park, K.; Ferrer, J.; Barraza-Lopez, S. Systematic pseudopotentials from reference eigenvalue sets for DFT calculations. Comput. Mater. Sci. 2015, 98, 372–389. [Google Scholar] [CrossRef]
  27. Yan, X.; Rahman, S.; Rostami, M.; Tabasi, Z.A.; Khan, F.; Alodhayb, A.; Zhang, Y. Carbon Quantum Dot-Incorporated Chitosan Hydrogel for Selective Sensing of Hg2+ Ions: Synthesis, Characterization, and Density Functional Theory Calculation. ACS Omega 2021, 6, 23504–23514. [Google Scholar] [CrossRef]
  28. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision E.01; Gaussian, Inc.: Wallingford, UK, 2013. [Google Scholar]
  29. Viewing Molecular Orbital Calculations with GaussView: A Lab for First- or Second-Year Undergraduate Students. Available online: https://www.ionicviper.org/lab-experiment/viewing-molecular-orbital-calculations-gaussview-lab-first-or-second-year-undergradua (accessed on 25 May 2010).
  30. Tighadouini, S.; Roby, O.; Radi, S.; Lakbaibi, Z.; Saddik, R.; Mabkhot, Y.N.; Almarhoon, Z.M.; Garcia, Y. A Highly Efficient Environmental-Friendly Adsorbent Based on Schiff Base for Removal of Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study. Molecules 2021, 26, 17. [Google Scholar] [CrossRef]
  31. Kruse, H.; Goerigk, L.; Grimme, S. Why the Standard B3LYP/6-31G* Model Chemistry Should Not Be Used in DFT Calculations of Molecular Thermochemistry: Understanding and Correcting the Problem. J. Org. Chem. 2012, 77, 10824–10834. [Google Scholar] [CrossRef]
  32. Ibrahim, S.M.; Halim, S.A. Novel SnZr oxides nanomaterials synthesized by ultrasonic-assisted co-precipitation method: Application in biodiesel production and DFT study. J. Mol. Liq. 2021, 339, 116652. [Google Scholar] [CrossRef]
  33. Check, C.E.; Faust, T.O.; Bailey, J.M.; Wright, B.J.; Gilbert, T.M.; Sunderlin, L.S. Addition of Polarization and Diffuse Functions to the LANL2DZ Basis Set for P-Block Elements. J. Phys. Chem. A 2001, 105, 8111–8116. [Google Scholar] [CrossRef]
  34. Zhou, Y.; Luan, L.; Tang, B.; Niu, Y.; Qu, R.; Liu, Y.; Xu, W. Fabrication of Schiff base decorated PAMAM dendrimer/magnetic Fe3O4 for selective removal of aqueous Hg(II). Chem. Eng. J. 2020, 398, 125651. [Google Scholar] [CrossRef]
  35. Halilu, A.; Hayyan, M.; Aroua, M.K.; Yusoff, R.; Hizaddin, H.F. In Situ Electrosynthesis of Peroxydicarbonate Anion in Ionic Liquid Media Using Carbon Dioxide/Superoxide System. ACS Appl. Mater. Interfaces 2019, 11, 25928–25939. [Google Scholar] [CrossRef]
  36. Tighadouini, S.; Roby, O.; Mortada, S.; Lakbaibi, Z.; Radi, S.; Al-Ali, A.; Faouzi, M.E.A.; Ferbinteanu, M.; Garcia, Y.; Al-Zaqri, N.; et al. Crystal structure, physicochemical, DFT, optical, keto-enol tautomerization, docking, and anti-diabetic studies of (Z)-pyrazol β-keto-enol derivative. J. Mol. Struct. 2022, 1247, 131308. [Google Scholar] [CrossRef]
  37. Chamorro, E.; Pérez, P.; Domingo, L.R. On the nature of Parr functions to predict the most reactive sites along organic polar reactions. Chem. Phys. Lett. 2013, 582, 141–143. [Google Scholar] [CrossRef]
  38. Matta, C.F. On the connections between the quantum theory of atoms in molecules (QTAIM) and density functional theory (DFT): A letter from Richard F. W. Bader to Lou Massa. Struct. Chem. 2017, 28, 1591–1597. [Google Scholar] [CrossRef]
  39. Savin, A.; Silvi, B.; Colonna, F. Topological analysis of the electron localization function applied to delocalized bonds. Can. J. Chem. 1996, 74, 1088–1096. [Google Scholar] [CrossRef]
  40. Savin, P.-D.A.; Nesper, R.; Wengert, S.; Fässler, T.F. ELF: The Electron Localization Function. Angew. Chem. Int. Ed. Engl. 1997, 36, 1808–1832. [Google Scholar] [CrossRef]
  41. Jacobsen, H. Localized-orbital locator (LOL) profiles of transition-metal hydride and dihydrogen complexes. Can. J. Chem. 2009, 87, 965–973. [Google Scholar] [CrossRef]
  42. Schmider, H.; Becke, A. Chemical content of the kinetic energy density. J. Mol. Struct. theochem. 2000, 527, 51–61. [Google Scholar] [CrossRef]
  43. Radi, S.; Tighadouini, S.; Bacquet, M.; Degoutin, S.; Garcia, Y. New Hybrid Material Based on a Silica-Immobilised Conjugated β-Ketoenol–Bipyridine Receptor and Its Excellent Cu(II) Adsorption Capacity. Anal. Methods 2016, 8, 6923–6931. [Google Scholar] [CrossRef]
  44. Radi, S.; Tighadouini, S.; Bacquet, M.; Degoutin, S.; Dacquin, J.-P.; Eddike, D.; Tillard, M.; Mabkhot, Y.N. β-Keto-Enol Tethered Pyridine and Thiophene: Synthesis, Crystal Structure Determination and Its Organic Immobilization on Silica for Efficient Solid-Liquid Extraction of Heavy Metals. Molecules 2016, 21, 888. [Google Scholar] [CrossRef]
  45. Radi, S.; Toubi, Y.; Attayibat, A.; Bacquet, M. New Polysiloxane-Chemically Immobilized C,C-Bipyrazolic Receptor for Heavy Metals Adsorption. J. App. Polym. Sci. 2011, 121, 1393–1399. [Google Scholar] [CrossRef]
  46. Brunauer, S.; Emmett, P.H.; Teller, E. Adsorption of Gases in Multimolecular Layers. J. Am. Chem. Soc. 1938, 60, 309–319. [Google Scholar] [CrossRef]
  47. Largitte, L.; Pasquier, R. A Review of the Kinetics Adsorption Models and Their Application to the Adsorption of Lead by an Activated Carbon. Chem. Eng. Res. Des. 2016, 109, 495–504. [Google Scholar] [CrossRef]
  48. Radi, S.; El Abiad, C.; Moura, N.M.; Faustino, M.A.; Neves, M.G.P. New Hybrid Adsorbent Based on Porphyrin Functionalized Silica for Heavy Metals Removal: Synthesis, Characterization, Isotherms, Kinetics and Thermodynamics Studies. J. Hazard. Mater. 2019, 370, 80–90. [Google Scholar] [CrossRef]
  49. Langmuir, I. The adsorption of gases on plane surface of glass, mica and platinum. J. Am. Chem. Soc. 1918, 40, 1361–1368. [Google Scholar] [CrossRef]
  50. Ali, O.; Mohamed, S. Adsorption of copper ions and alizarin red S from aqueous solutions onto a polymeric nanocomposite in single and binary systems. Turk. J. Chem. 2017, 41, 967–986. [Google Scholar] [CrossRef]
  51. Milonjic, S. A consideration of the correct calculation of thermodynamic parameters of adsorption. J. Serb. Chem. Soc. 2007, 72, 1363–1368. [Google Scholar] [CrossRef]
  52. Chen, Z.; Tang, B.; Niu, Y.; Chen, H.; Liu, Y.; Wang, A.; Bai, L. Synthesis of silica supported thiosemicarbazide for Cu(II) and Zn(II) adsorption from ethanol: A comparison with aqueous solution. Fuel 2021, 286, 119287. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Cao, X.; Sun, J.; Wu, G.; Wang, J.; Zhang, D. Synthesis of pyridyl Schiff base functionalized SBA-15 mesoporous silica for the removal of Cu(II) and Pb(II) from aqueous solution. J. Sol. Gel. Sci. Technol. 2019, 94, 658–670. [Google Scholar] [CrossRef]
  54. G. In situ synthesis of nanocomposite materials based on modified-mesoporous silica MCM-41 and methyl methacrylate for copper(II) adsorption from aqueous solution. J. Iran. Chem. Soc. 2019, 16, 1491–1500. [Google Scholar] [CrossRef]
  55. He, S.; Zhao, C.; Yao, P.; Yang, S. Chemical modification of silica gel with multidentate ligands for heavy metals removal. Desalin. Water Treat. 2016, 57, 1722–1732. [Google Scholar] [CrossRef]
  56. Moftakhar, M.K.; Dousti, Z.; Yaftian, M.R.; Ghorbanloo, M. Investigation of heavy metal ions adsorption behavior of silica-supported Schiff base ligands. Desalin. Water Treat. 2016, 57, 27396–27408. [Google Scholar] [CrossRef]
  57. Radi, S.; Tighadouini, S.; El Massaoudi, M.; Bacquet, M.; Degoutin, S.; Revel, B.; Mabkhot, Y.N. Thermodynamics and kinetics of heavy metals adsorption on silica particles chemically modified by conjugated β-ketoenol furan. J. Chem. Eng. Data 2015, 60, 2915. [Google Scholar] [CrossRef]
  58. Dong, C.; Fu, R.; Sun, C.; Qu, R.; Ji, C.; Niu, Y.; Zhang, Y. Comparison studies of adsorption properties for copper ions in fuel ethanol and aqueous solution using silica-gel functionalized with 3-amino-1,2-propanediol. Fuel 2018, 226, 331–337. [Google Scholar] [CrossRef]
  59. Morcali, M.H.; Zeytuncu, B.; Baysal, A.; Akman, S.; Yucel, O. Adsorption of copper and zinc from sulfate media on a commercial sorbent. J. Environ. Chem. Eng. 2014, 2, 1655–1662. [Google Scholar] [CrossRef]
  60. El Massaoudi, M.; Radi, S.; Bacquet, M.; Degoutin, S.; Adarsh, N.N.; Robeyns, K.; Garcia, Y. A novel environment-friendly hybrid material based on a modified silica gel with a bispyrazole derivative for the removal of Zn(II), Pb(II), Cd(II) and Cu(II) traces from aqueous solutions. Inorg. Chem. Front. 2017, 4, 1821. [Google Scholar]
  61. Weinhold, F.; Landis, C.R. Natural bond orbitals and extensions of localized bonding concepts. Chem. Educ. Res. Pract. 2001, 2, 91–104. [Google Scholar] [CrossRef]
  62. Hammoudan, I.; Chtita, S.; Riffi-Temsamani, D. QTAIM and IRC studies for the evaluation of activation energy on the C=P, C=N and C=O Diels-Alder reaction. Heliyon 2020, 6, e04655. [Google Scholar] [CrossRef]
Scheme 1. Synthesis route of modified chelating material.
Scheme 1. Synthesis route of modified chelating material.
Molecules 27 05168 sch001
Figure 1. FTIR spectra of SiNH2 and Si-Imd-Py (before and after desorption).
Figure 1. FTIR spectra of SiNH2 and Si-Imd-Py (before and after desorption).
Molecules 27 05168 g001
Figure 2. SEM micrographs of SiG, SiNH2 and SiN-imd-py.
Figure 2. SEM micrographs of SiG, SiNH2 and SiN-imd-py.
Molecules 27 05168 g002
Figure 3. Thermogravimetric curves of free silica SiG, SiNH2 and SiN-imd-py.
Figure 3. Thermogravimetric curves of free silica SiG, SiNH2 and SiN-imd-py.
Molecules 27 05168 g003
Figure 4. Nitrogen adsorption–desorption isotherm plots of SiG, SiNH2 and SiN-imd-py.
Figure 4. Nitrogen adsorption–desorption isotherm plots of SiG, SiNH2 and SiN-imd-py.
Molecules 27 05168 g004
Figure 5. The effect of pH on the adsorption capacity of SiN-imd-py towards Cu(II). Adsorption condition: V = 10 mL, m = 10 mg of adsorbent and optimum concentration = 150 mg L−1 for Cu(II) and t = 25 min at 25 °C.
Figure 5. The effect of pH on the adsorption capacity of SiN-imd-py towards Cu(II). Adsorption condition: V = 10 mL, m = 10 mg of adsorbent and optimum concentration = 150 mg L−1 for Cu(II) and t = 25 min at 25 °C.
Molecules 27 05168 g005
Figure 6. Effect of contact time on the adsorption capacity of Cu(II). Adsorption conditions: V = 10 mL, m = 10 mg of adsorbent, pH = 6, and optimum concentrations: 150 mg L−1 for Cu(II) at 25 °C.
Figure 6. Effect of contact time on the adsorption capacity of Cu(II). Adsorption conditions: V = 10 mL, m = 10 mg of adsorbent, pH = 6, and optimum concentrations: 150 mg L−1 for Cu(II) at 25 °C.
Molecules 27 05168 g006
Figure 7. Pseudo-second-order model fit for the adsorption of Cu(II) by SiN-imd-py.
Figure 7. Pseudo-second-order model fit for the adsorption of Cu(II) by SiN-imd-py.
Molecules 27 05168 g007
Figure 8. Effect of concentration on metal ion adsorption onto SiN-imd-py adsorbent. Adsorption conditions: m = 10 mg, V = 10 mL, [Cu(II)] = 10 to 300 mg L−1, pH = 6, time = 25 min at 25 °C.
Figure 8. Effect of concentration on metal ion adsorption onto SiN-imd-py adsorbent. Adsorption conditions: m = 10 mg, V = 10 mL, [Cu(II)] = 10 to 300 mg L−1, pH = 6, time = 25 min at 25 °C.
Molecules 27 05168 g008
Figure 9. Langmuir adsorption model fit of Cu(II) on SiN-imd-py.
Figure 9. Langmuir adsorption model fit of Cu(II) on SiN-imd-py.
Molecules 27 05168 g009
Figure 10. Effect of temperature for the adsorption of metal ions onto SiN-imd-py adsorbent.
Figure 10. Effect of temperature for the adsorption of metal ions onto SiN-imd-py adsorbent.
Molecules 27 05168 g010
Figure 11. Metal ions selectivity effect for SiN-imd-py. Adsorption conditions: contact time = 25 min, V = 10 mL, pH = 6 at 25 °C, m = 10 mg of adsorbent at optimum concentrations: 150 ppm of each studied metal Cu(II), Pb(II), Zn(II) and Cd(II).
Figure 11. Metal ions selectivity effect for SiN-imd-py. Adsorption conditions: contact time = 25 min, V = 10 mL, pH = 6 at 25 °C, m = 10 mg of adsorbent at optimum concentrations: 150 ppm of each studied metal Cu(II), Pb(II), Zn(II) and Cd(II).
Molecules 27 05168 g011
Figure 12. NCI isosurface (left), contour maps of electron density (right) of SiN-Imd-Py.
Figure 12. NCI isosurface (left), contour maps of electron density (right) of SiN-Imd-Py.
Molecules 27 05168 g012
Figure 13. Optimized structures and BCP (bond critical points) for two complexes 1 (left) and 2 (right).
Figure 13. Optimized structures and BCP (bond critical points) for two complexes 1 (left) and 2 (right).
Molecules 27 05168 g013
Figure 14. Curve showing various properties on the line (Cu--N) of the two complexes (left complex, right complex).
Figure 14. Curve showing various properties on the line (Cu--N) of the two complexes (left complex, right complex).
Molecules 27 05168 g014
Figure 15. Contour map of electron density (left) Laplacian and localized-orbital locator (LOL) (right).
Figure 15. Contour map of electron density (left) Laplacian and localized-orbital locator (LOL) (right).
Molecules 27 05168 g015
Table 1. Elemental analysis of SiN-imd-py and SiNH2.
Table 1. Elemental analysis of SiN-imd-py and SiNH2.
ElementSiNH2
%
SiN-imd-py
%
H3.675.55
N9.0216.56
C23.9749.45
Table 2. Physical Properties of Silica Derivatives.
Table 2. Physical Properties of Silica Derivatives.
Silica DerivativesSpecific Surface SBET (m2 g−1)Pore Volume (cm3 g−1)
Free silica305.210.770
SiNH2283.080.690
SiN-Imd-Py231.410.541
Table 3. Kinetics models data of Cu(II) adsorption.
Table 3. Kinetics models data of Cu(II) adsorption.
Metalqe (exp) (mg g−1)Pseudo-First-OrderPseudo-Second-Order
k1 (min−1)qe (mg g−1)k2 (g/mg min)qe (mg g−1)
Cu(II)97.170.160544.94770.93720.005208104.05820.9989
Table 4. Parameters for the Langmuir, Freundlich, Dubini–Radushkevich and Temkin of Cu(II) adsorption.
Table 4. Parameters for the Langmuir, Freundlich, Dubini–Radushkevich and Temkin of Cu(II) adsorption.
MetalLangmuir modelFreundlich model
qe (mg g−1)KL (L mg−1)R2KF (mg g−1)NR2
Cu(II)103.51960.14770.998330.80573.74320.9605
Dubinin–Radushkevich (D-R) modelTemkin model
β (mol2 Kj−2)E (Kj mol−1)R2At (L mg1)bt (J mol−1)R2
Cu(II)4.37 × 10−71069.650.63694.5038152.32620.9829
Table 5. Thermodynamic parameters of SiN-imd-py and other similar adsorbents at various temperatures.
Table 5. Thermodynamic parameters of SiN-imd-py and other similar adsorbents at various temperatures.
MétalAdsorbentsΔH° (kJ mol−1)ΔS° (Jk−1mol−1)ΔG° (kJ mol−1)
298 K308 K318 K
Cu(II)SiN-imd-py
SiPy [9]
3.778
8.1364
14.142
29.163
−0.452
−0.558
−0.593
−0.850
−0.735
−1.141
SG-GPTS-ATS [52]48.99221.86−14.91−17.12−19.32
Table 6. Reusability and recycling of SiN-Imd-Py adsorbent towards Cu(II) in the adsorption–desorption cycles.
Table 6. Reusability and recycling of SiN-Imd-Py adsorbent towards Cu(II) in the adsorption–desorption cycles.
Cycleqe (mg g−1) of Cu(II) Adsorbed on SiN-imd-py
197.17
292.18
390.02
489.07
588.43
Table 7. Comparison of the maximum adsorption capacities of Cu(II) by different adsorbents reported in the literature.
Table 7. Comparison of the maximum adsorption capacities of Cu(II) by different adsorbents reported in the literature.
Silica Gel-LigandReferenceMetal Ion (mg g−1)
2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde
Pyridin-2-ylmethanol
Porphyrin
This work
[9]
[48]
97.17
90.25
19.08
N-propyl-2-pyridylimine[53]35.63
Methyl methacrylate[54]41.36
Dithiocarbamate[55]25.00
(E)-4-(furan-2-ylmethyleneamino) phenol
(E)-2-(furan-2-ylmethyleneamino) phenol
[30]
[30]
36.20
79.36
3-Hydroxysalicylaldiminepropyltriethoxy-silane[56]5.72
Furan ketonenol[57]31.82
3-amino-1,2-propanediol
Commercial Lewatit (L-207)
Bis(pyrazole)butane
[58]
[59]
[60]
31.18
68.09
20.24
Table 8. (p-) of significant atoms in a.u.
Table 8. (p-) of significant atoms in a.u.
2N9N11N
P(-)−0.0640.0840.293
Table 9. NBO analysis calculations for complexes 1 and 2.
Table 9. NBO analysis calculations for complexes 1 and 2.
E(2) Energy in (kcal/mol)
EC of Cu, N9, O39, O42, O45LPN(9)-LP*Cu(38)LPO(39)-RY*Cu(38)LPO(42)-RY*Cu(38)LPO(45)-LP*Cu(38)
C1Cu[core]4S(0.10)3d(4.96)4p(0.11)
N9[core]2S(1.37)2p(4.39)3S(0.01)3p(0.15)3d(0.01)4p(0.02)5p(0.01)
O39[core]2S(1.72)2p(5.21)3p(0.01)
O42[core]2S(1.71)2p(5.21)3p(0.01)
O45[core]2S(1.72)2p(5.21)3p(0.01)
20.7835.0711.92565.01
C2Cu[core]4S(0.20)3d(9.84)4p(0.10)
N9[core]2S(1.36)2p(4.39)3S(0.01)3p(0.10)3d(0.02)4p(0.01)
O39[core]2S(1.72)2p(5.21)3p(0.01)
O42[core]2S(1.72)2p(5.21)3p(0.01)
O45[core]2S(1.72)2p(5.21)3p(0.01)
20.130.210.1555.32
C1: complex1. C2: complex2.
Table 10. Calculation QTAIM, ELF, LOL parameters of the cps and length between atoms in a.u.
Table 10. Calculation QTAIM, ELF, LOL parameters of the cps and length between atoms in a.u.
CPBondingd (Å)ρ∇²ρLagELFLOL
Complex1859N-Cu1.9830.8320.1350.5100.1020.252
8939O-Cu2.1710.3890.2320.5320.5490.194
7842O-Cu2.1550.4610.3160.7100.5440.193
7645O-Cu2.2410.4510.3020.6780.5530.194
Complex2849N-Cu1.9950.8080.4950.1290.1000.250
8839O-Cu2.1690.3830.2310.5170.5500.194
7942O-Cu2.1790.4390.2930.6560.5400.193
7545O-Cu2.2490.4530.3040.6840.5530.194
Lag: Lagrangian, d: bond lengths, ρ: electron density, CP: critical points.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saddik, R.; Hammoudan, I.; Tighadouini, S.; Roby, O.; Radi, S.; Al-Zaben, M.I.; Ben Bacha, A.; Masand, V.H.; Almarhoon, Z.M. Mesoporous Silica Modified with 2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde as an Effective Adsorbent for Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study. Molecules 2022, 27, 5168. https://doi.org/10.3390/molecules27165168

AMA Style

Saddik R, Hammoudan I, Tighadouini S, Roby O, Radi S, Al-Zaben MI, Ben Bacha A, Masand VH, Almarhoon ZM. Mesoporous Silica Modified with 2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde as an Effective Adsorbent for Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study. Molecules. 2022; 27(16):5168. https://doi.org/10.3390/molecules27165168

Chicago/Turabian Style

Saddik, Rafik, Imad Hammoudan, Said Tighadouini, Othmane Roby, Smaail Radi, Maha I. Al-Zaben, Abir Ben Bacha, Vijay H. Masand, and Zainab M. Almarhoon. 2022. "Mesoporous Silica Modified with 2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde as an Effective Adsorbent for Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study" Molecules 27, no. 16: 5168. https://doi.org/10.3390/molecules27165168

Article Metrics

Back to TopTop