Liz Nayibe Martínez Saavedraa,
Ricardo Gonçalves Penidoa,
Lucas de Azevedo Santosb,
Teodorico C. Ramalhob,
Bruno E. Lobo Baetaa,
Márcio C. Pereirac and
Adilson Candido da Silva*a
aDepartamento de Química, Instituto de Ciências Exatas e Biológicas, Universidade Federal de Ouro Preto, 35400-000 Ouro Preto, Minas Gerais, Brazil. E-mail: adilsonqui@iceb.ufop.br; adilsonufla@gmail.com; Tel: +55 31 3559 1934
bDepartamento de Química, Universidade Federal de Lavras, 37200-000 Lavras, Minas Gerais, Brazil
cInstituto de Ciência, Universidade Federal dos Vales do Jequitinhonha e Mucuri, Engenharia e Tecnologia, 39803-371 Teófilo Otoni, Minas Gerais, Brazil
First published on 14th August 2018
The effects of solvent on the synthesis of molecularly imprinted polymers (MIPs) for the selective adsorption of quinoline were evaluated in this work. The MIPs were synthesized by the “bulk” method using the quinoline molecule (IQ) as a template in different solvents, such as toluene (MIPT) and chloroform (MIPC). The adsorbents were characterized by thermogravimetric analysis (TGA), Fourier transform infrared (FT-IR) spectroscopy, scanning electron microscopy (SEM), and N2 adsorption/desorption measurements. The influences of time, adsorbate concentration, and temperature on the adsorption of quinoline by MIPT and MIPC were evaluated. Maximum adsorption capacities (qe) of 35.23 and 24.10 mg g−1 were obtained for MIPT and MIPC, respectively. Thermodynamic studies indicate that occur physisorption and a spontaneous process (ΔadsG° < 0) entropically directed. Finally, the highest selectivity and reusability of MIPC for quinoline adsorption was ascribed to the better interaction between the chloroform and monomer, which favors the formation of porous adsorbents with higher numbers of adsorption sites.
The low removal rates of sulfide compounds in these processes are attributed to the presence of nitrogen compounds, which can strongly inhibit the reactions of the hydroprocessing by competitive adsorption. This occurs because the nitrogen atoms in the contaminants strongly adsorb on the active sites of the catalyst. Quinoline is one of the most powerful inhibitors of the HDT process due to poisoning of the catalysts.4,5 Adsorption treatments using porous materials, such as zeolites, activated carbon, activated alumina, silica gel, kaolinite, and montmorillonite, are the preferred alternatives for the removal of these contaminants; however, they are not selective.4,5 On the other hand, molecularly imprinted polymers (MIPs) are highly selective adsorbents because they are capable of forming complementary sites with the target molecule on their surfaces, which leads to the formation of specific cavities in the polymer network. Because of their benefits, such as selectivity, stability and relatively easy synthesis, the use of MIPs to remove specific compounds from environmental samples has attracted attention in recent years.6,7 Extensive efforts have been directed at the synthesis of these materials, as well as adsorption studies in organic phase8–11 and aqueous phase.12 It is important to note that the selective adsorption of quinoline favors its recovery and, consequently, its use as a molecular platform for the synthesis of other compounds of interest, such as quinolinic acid. The primary parameter studied during the synthesis of MIPs is the solvent. The ideal solvent should ensure selective interaction between the monomer and the target molecule. Thus, careful choice of solvent is required, with adequate polarity to favor the interactions between the monomer and the template molecule. The synthesis of the adsorbent must be performed in solvents with polarities similar to the medium in which the material will be used. Another characteristic that must be evaluated in a solvent for the synthesis of MIPs is its vapor pressure; solvents with higher vapor pressures are more readily volatilized, favoring pore formation during the synthesis and resulting in adsorbents with greater surface areas. Thus, it is fundamental to study the solvent that will be used during the synthesis when the aim is the production of a MIP with high selectivity and adsorption capacity.
In view of the above, the present study aims to evaluate the influence of the solvent in obtaining selective and efficient MIPs for the adsorption of quinoline in an organic matrix.
The infrared spectra were obtained in the 400 to 4000 cm−1 range with an ABB Bomem brand instrument coupled with accessories for attenuated total reflectance (ATR) and diffuse reflectance (DRIFT) analysis. The thermogravimetric analyses were performed on SDT2960 Simultaneous TGA-DTA equipment, model 2960, TA Instruments, with temperature variation from 25 °C to 1000 °C and a heating rate of 10 °C per minute in an inert atmosphere using N2. The micrographs were produced at different magnification levels using a JEOL/EO scanning electron microscope (SEM), model JSM-5500, version 1.0, with a coupled reflex photo capture system. The specific surface areas and porosities of the materials were measured using AutosorbIQ equipment.
Removal of the template after the polymerization was performed in a Soxhlet extractor by washing the material with a methanol/acetic acid (9:1 v/v) mixture. The final time of extraction was considered to be that at which the presence of residual template in the final extraction solution was not observed. After maximum removal of the target molecule, the polymer was dried at 40 °C for 24 hours and stored at room temperature.
Fig. 1 schematically illustrates the two possible molecular imprinting procedures in MIPs, i.e. molecular interaction through hydrogen bonding or electrostatic interaction with ion pair formation.
The quantification of quinoline was performed in triplicate before and after the adsorption assays to determine the capacity of adsorption, qe (mg g−1), of the adsorbent in equilibrium. The Langmuir, Freundlich, Sips, and Temkin models were used to fit the adsorption isotherms.
The results of the kinetic studies were evaluated by adjusting the data to the order n (PON) model, expressing the adsorption rate and the number of active sites per adsorbent particle.13 The pseudo-second order model was also evaluated as described by Ho and Mckay (1999), and the rate constant k2ads (g mg−1 min−1) was determined. In addition, the Elovich model was evaluated according to equations proposed by Roginsky and Zeldovich in 1934; the constant a (mg g−1 min−1) is considered to be the initial rate of the adsorption process, and b (g mg−1) is a parameter related to the covered surface and the energy activation. The intraparticle diffusion model was tested according to Weber and Morris (1963) by determining Kpi (g mg−1 min−1/2), which is the diffusion rate constant between the particles of step i, where Ci is the intersection of step i. The best model was determined by analysis of the error functions X2, Adj. R2 and RMSE.
The Freundlich model suggests that the binding sites are not independent and assumes that the surface is energetically heterogeneous, allowing interactions between the adsorbed molecules20 (eqn (1)):
q = KCf1/n | (1) |
The Sips model is a versatile expression that can simulate behaviors of the Langmuir and Freundlich models. It is a model of three parameters, deduced for the prediction of heterogeneous adsorption systems, and avoids limiting the concentration of adsorbate output in the Freundlich model.15 Eqn (2) represents the mathematical model:
(2) |
(3) |
ΔG° = −RTlnKp | (4) |
(5) |
The Van't Hoff equation shows the dependence between the adsorption equilibrium constant and the temperature.17
Through the parameters outlined in eqn (6)–(8), specific coefficients were determined to evaluate the selectivities of the imprinted adsorbents.
The distribution coefficient Kd (L g−1) for each compound was calculated by the equation
(6) |
The selectivity coefficient k of quinoline concerning the analogous molecule (identified as X) was calculated by eqn (7), as follows:
(7) |
The relative selectivity coefficient k′ in relation to MIP and NIP was calculated as follows:
(8) |
Fig. 2 (a) TGA profiles and (b) FTIR spectra of quinoline and MIPs before extraction, and (c) FTIR spectra of MIPs and NIPs after extraction. |
Fig. 2b and c show the FTIR spectra of the polymers. For all polymers, similarity was observed in the bands. A wide band at 3170 to 3672 cm−1 with a maximum of 3500 cm−1, attributed to the stretching (O–H) mode of the hydroxyl groups of the MAA functional monomer,22 was observed. Bands at 2900 cm−1, resulting from (C–H) stretching, and at 1384 cm−1, resulting from (CH2) stretching/distortion, were observed.23 Strong bands appeared at 1750 (CO), 1230, and 1125 cm−1 (C–O–C), which are characteristic of MAA and ester groups (EGDMA), respectively.24 At 1480 cm−1, a band originating from the stretching (–COO) of the MAA was observed. In Fig. 2b, it can be observed that the band at 1480 cm−1 attributed to the vibration of the COO group of the monomer has a lower intensity than the same band shown in Fig. 2c. This may occur due to the interaction of quinoline with the COO group of the monomer, as suggested in Fig. 1. It is worth emphasizing that variation in the intensity of a band is not enough to confirm the presence or absence of the template. However, for these types of materials, the mass fraction of the template (quinolone) in the MIP is very low because all other compounds were added in excess during the synthesis of the material. Therefore, the specific vibrations of the monomers, radical initiators and cross-link agents are more intense, which makes it difficult to see specific vibrations between the monomer and template. However, a change in the chemical environment of the system with the insertion of a new interaction between the monomer and the template may decrease the intensity of the –COO vibration frequencies because as shown by theoretical calculations, the interaction between the template and the monomer occurs at the carboxylic acid function of methacrylic acid. The bands at 937, 875 and 750 cm−1 are due to out-of-plane deformation of (C–H) on different substituted benzene rings. The spectra initially indicate that the quinoline extraction was successful because the spectra of the MIPs with quinoline molecular imprinting are similar to the spectra of the non-molecular imprinted NIPs.
Table 1 shows the specific areas (BET), pore sizes, and volume data for the adsorbents. According to the results, the MIP that developed the largest specific area was the MIPC (205.90 m2 g−1, which is twice the value obtained for the MIPT). The lower boiling point of chloroform (61.2 °C) compared to toluene (110.6 °C) can explain this result. Solvents with higher vapor pressures are responsible for increasing the surface area of the adsorbent during synthesis, contributing to better performance.
In both MIPs, there was greater development in the amount of mesopores (20 to 500 Å) compared to micropores (<20 Å). The mean micropore size for all the materials was 17 Å, and that of the mesopores was 120 to 130 Å. Note that the specific areas of the NIPs were smaller than those obtained for the corresponding MIPs, indicating that the porogenic solvent in the presence of the template probably develops a better porous structure. Thus, it is possible to infer that the imprinted polymers are mesoporous adsorbents, which have exciting characteristics for use in adsorption processes (Table 1).
Polymer | Specific area – BET (m2 g−1) | Average pore sizes (Å) | Pore volume (cm3 g−1) | |||
---|---|---|---|---|---|---|
Micro | Meso | Micro | Meso | Total | ||
NIPT | 81.8 | 17.7 | 131.15 | 6.5 × 10−2 | 1.12 | 1.18 |
MIPT | 100.0 | 17.2 | 132.35 | 7.0 × 10−3 | 2.14 | 2.13 |
NIPC | 161.1 | 17.7 | 105.59 | 5.9 × 10−2 | 1.17 | 1.23 |
MIPC | 205.9 | 17.8 | 133.29 | 7.3 × 10−2 | 1.06 | 1.13 |
The SEM images obtained for the materials are shown in Fig. 3. The different micrometric scales of 10 μm to 20 μm allow us to determine the external morphology of each of the materials with increases of 800 up to 1000 times, showing the macroporous cavities. The macropores indicate uniformity and similarity between the NIPs and MIPs. The general surfaces of all the polymers are exhibited, showing spherical particles of different sizes. For the materials synthesized with toluene, the surfaces appear to be wider, whereas for the materials synthesized with chloroform, the surfaces appear to be more complex and narrow.
Table 2 shows the parameters obtained for each model. It can be observed that the pseudo-order model n (PON) best fit the adsorption process of MIPT. The parameter n is responsible for measuring the order of the reaction; the value of this variable for MIPT was close to 1, which indicates that the pseudo-first order equation could also fit the experimental data. The adjustment of this model may indicate that the adsorption process for materials synthesized with toluene occurs with the occupation of a single active site on a solid and energetically homogeneous surface.28
Model | Polymer | NIPT | MIPT | NIPC | MIPC |
---|---|---|---|---|---|
qe,exp (mg g−1) | 14.28 | 15.48 | 11.29 | 12.85 | |
Pseudo second order | qe,calc (mg g−1) | 16.62 | 18.52 | 10.70 | 12.24 |
K2 (min−1) | 1.9 × 10−3 | 1.3 × 10−3 | 0.08 | 0.03 | |
X2 | 0.90 | 1.38 | 0.19 | 0.26 | |
Adj. R2 | 0.95 | 0.95 | 0.53 | 0.79 | |
R2 | 0.97 | 0.98 | 0.59 | 0.82 | |
RMSE | 0.31 | 0.39 | 0.14 | 0.17 | |
Pseudo order n | qe,calc (mg g−1) | 14.91 | 15.66 | 11.25 | 13.08 |
kn (min−1) | 0.01 | 0.02 | 0.01 | 3.6 × 10−3 | |
n | 1.27 | 0.98 | 2.85 | 2.51 | |
X2 | 0.51 | 0.45 | 4.8 × 10−3 | 3.8 × 10−3 | |
Adj. R2 | 0.95 | 0.96 | 0.68 | 0.90 | |
R2 | 0.98 | 0.98 | 0.75 | 0.92 | |
RMSE | 0.24 | 0.38 | 0.02 | 0.02 | |
Elovich | a (mg g−1 min−1) | 1.38 | 1.20 | 1.4 × 106 | 989.53 |
b (g mg−1) | 0.32 | 0.28 | 1.83 | 0.97 | |
X2 | 1.05 | 2.34 | 0.03 | 8.3 × 10−3 | |
Adj. R2 | 0.92 | 0.91 | 0.96 | 0.99 | |
R2 | 0.96 | 0.96 | 0.97 | 0.99 | |
RMSE | 0.34 | 0.51 | 0.05 | 0.02 | |
Intraparticle diffusion | K1p | 1.33 | 1.68 | 0.16 | 0.06 |
K2p | 0.001 | 0.11 | 0.31 | 0.13 |
When comparing the values of the kinetic constant kn of NIPT (0.015 min−1) and MIPT (0.023 min−1), it is possible to note that for MIPT, adsorption occurred more rapidly; this behavior may be related to the lower limit of mass transfer of the quinoline molecules to the selective sites compared to non-specific sites generated in the NIPT polymer synthesized without molecular imprinting.
The model that best fit the adsorption process for MIPC was the Elovich model. This model suggests that the active sites of the respective adsorbents are of a heterogeneous nature and, therefore, exhibit different activation energies in the adsorption process.29 The constant a is related to the adsorption rate, and the constant b is related to the surface coverage.30 For NIPC, the reaction rate was faster and the surface coverage was higher, which indicates an increase of the adsorption surface available to the adsorbent with increasing adsorption. The opposite occurred for MIPC, which indicates a decrease in available sites for quinoline with increasing time in the adsorption process.
The values of the Kp diffusion rates for the corresponding adsorption systems were higher for the constant k1 than for k2, indicating that the quinoline adsorption was generally fast at the beginning and decreased with time. This explains why these materials have wide external surfaces.
Fig. 5 Non-linear Sips and Temkin models obtained at different temperatures for each of the MIPT and MIPC adsorption isotherms. |
Polymer | NIPC | MIPC | |||||||
---|---|---|---|---|---|---|---|---|---|
Model | Parameters | Temperature (K) | |||||||
288.15 | 298.15 | 308.15 | 318.15 | 288.15 | 298.15 | 308.15 | 318.15 | ||
qexp (mg g−1) | 7.86 | 11.88 | 11.46 | 13.69 | 16.99 | 22.43 | 24.10 | 21.73 | |
Freundlich | Kf (mg g−1) (L g−1)1/n | 1.01 | 0.70 | 1.11 | 1.48 | 1.44 | 1.08 | 1.91 | 2.00 |
n (L mg−1) | 1.98 | 1.37 | 1.69 | 1.70 | 1.50 | 1.15 | 1.32 | 1.45 | |
X2 | 0.42 | 1.02 | 1.59 | 2.30 | 1.91 | 2.03 | 7.34 | 4.45 | |
Adj. R2 | 0.91 | 0.92 | 0.86 | 0.87 | 0.93 | 0.95 | 0.87 | 0.90 | |
R2 | 0.95 | 0.96 | 0.93 | 0.94 | 0.96 | 0.98 | 0.94 | 0.95 | |
RMSE | 0.64 | 1.01 | 1.26 | 1.51 | 1.38 | 1.42 | 2.70 | 2.11 | |
Langmuir | Qmax (mg g−1) | 10.90 | 27.51 | 18.34 | 22.62 | 32.48 | 44.79 | 29.78 | 44.25 |
b (L mg−1) | 0.01 | 0.02 | 0.03 | 0.04 | 0.01 | 0.01 | 0.02 | 0.02 | |
X2 | 0.15 | 0.61 | 0.92 | 1.23 | 0.93 | 1.04 | 5.35 | 2.62 | |
Adj. R2 | 0.96 | 0.95 | 0.92 | 0.93 | 0.96 | 0.96 | 0.90 | 0.94 | |
R2 | 0.98 | 0.98 | 0.96 | 0.96 | 0.98 | 0.98 | 0.95 | 0.97 | |
RMSE | 0.39 | 0.78 | 0.96 | 1.11 | 0.96 | 1.20 | 2.31 | 1.61 | |
Sips | Qmax (mg g−1) | 8.16 | 13.53 | 11.0 | 14.00 | 19.34 | 28.46 | 24.77 | 23.12 |
b (L mg−1) | 0.06 | 0.04 | 0.07 | 0.07 | 0.06 | 0.05 | 0.09 | 0.08 | |
n | 0.57 | 0.46 | 0.36 | 0.41 | 0.56 | 0.53 | 0.33 | 0.43 | |
X2 | 0.05 | 0.05 | 0.10 | 0.16 | 0.18 | 0.20 | 0.13 | 0.18 | |
Adj. R2 | 0.98 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | |
R2 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | 0.99 | |
RMSE | 0.23 | 0.23 | 0.32 | 0.40 | 0.43 | 0.45 | 0.37 | 0.43 | |
Temkin | bT (kJ mol−1) | 0.98 | 0.53 | 0.61 | 0.50 | 0.35 | 0.25 | 0.23 | 0.29 |
AT (L mg−1) | 0.39 | 0.23 | 0.27 | 0.30 | 0.29 | 0.25 | 0.30 | 0.33 | |
X2 | 0.15 | 0.92 | 0.66 | 0.78 | 0.28 | 1.60 | 3.57 | 1.62 | |
Adj. R2 | 0.96 | 0.93 | 0.94 | 0.95 | 0.98 | 0.96 | 0.93 | 0.96 | |
R2 | 0.97 | 0.97 | 0.97 | 0.98 | 0.99 | 0.98 | 0.97 | 0.98 | |
RMSE | 0.38 | 0.96 | 0.81 | 0.83 | 0.53 | 1.26 | 1.89 | 1.27 |
Polymer | ΔadsG (kJ mol−1) | TΔadsS° (kJ mol−1) | ΔadsH° (kJ mol−1) | |||
---|---|---|---|---|---|---|
Temperature (K) | ||||||
288.15 | 298.15 | 308.15 | 318.15 | |||
NIPT | −18.43 | −19.92 | −22.85 | −26.56 | 80.59 | 59.99 |
MIPT | −17.79 | −19.71 | −23.02 | −27.21 | 93.35 | 72.95 |
NIPC | −18.13 | −20.02 | −21.84 | −23.09 | 50.05 | 30.11 |
MIPC | −17.39 | −18.59 | −20.48 | −21.76 | 44.75 | 25.93 |
Polymer | Selectivity parameters | Target and analogous molecules | ||
---|---|---|---|---|
Quinoline | DBT | Carbazole | ||
NIPT | qexp (mg g−1) | 1.25 | 0.16 | 6.99 |
Kd (L g−1) | 0.02 | 0.003 | 0.19 | |
K | — | 7.90 | 0.13 | |
MIPT | qexp (mg g−1) | 10.29 | 7.09 | 9.13 |
Kd (L g−1) | 0.34 | 0.19 | 0.28 | |
k | — | 1.76 | 1.21 | |
k′ | — | 0.22 | 8.91 | |
NIPC | qexp (mg g−1) | 5.62 | 0.57 | 7.12 |
Kd (L g−1) | 0.14 | 0.01 | 0.19 | |
K | — | 12.36 | 0.72 | |
MIPC | qexp (mg g−1) | 13.61 | 5.26 | 7.65 |
Kd (L g−1) | 0.59 | 0.13 | 0.22 | |
k | — | 4.47 | 2.70 | |
k′ | — | 0.36 | 3.72 |
From Fig. 5c and d, an increase in adsorption for MIPC was observed at 298.15 and 308.15 K, with qe values of 22.43 and 24.10 mg g−1, respectively. The difference in the qe values of MIPC and NIPC was ∼13 mg g−1, indicating better specific recognition at 308.15 K. This shows that the synthesis conditions for these materials favored the formation of a more selective polymer, resulting in greater recognition capacity. Considering the results presented so far, it is worth pointing out that the use of chloroform as a porogenic solvent in the MIPC preparation favored the formation of specific interactions between the template and the monomer during the monomer/template active complex formation step. The model that best fit the MIPC data adsorption isotherm is the Sips model. The values of the parameter n presented in Table 3 are all lower than 1, which shows that the surfaces of these materials have an energetically heterogeneous distribution of active sites. The values of the Sips (n) parameter were compared with those obtained from the Freundlich (1/n) parameter; they were lower than 1, reinforcing the positive cooperation phenomenon between the adsorbent and the adsorbate presented on the heterogeneous surfaces resulting from the corresponding polymers. In addition, if the n values of the Freundlich model are greater than 1, favorable adsorption is predicted.33
ΔΔG | Volume | |||
---|---|---|---|---|
Clo | Tol | Clo | Tol | |
Quinoline | 0.00 | 0.00 | 1335.92 | 1320.61 |
DBT | 68.07 | 66.61 | 1609.87 | 1586.63 |
Carbazole | 36.72 | 46.32 | 1442.77 | 1325.16 |
Adsorption processes can be classified as physical or chemical according to the enthalpy change of adsorption ΔadsH°; if the resulting value is less than 80 kJ mol−1, the phenomenon that occurs in the mechanism of adsorption is physical, but if the value is in the range between 80 and 400 kJ mol−1, the phenomenon that occurs is chemisorption. Thus, the obtained adsorption enthalpy values (ΔadsH°) are in the range in which the adsorption process occurs by the physisorption phenomenon; this indicates the presence of stronger electrostatic interactions, probably ion-dipole, for the materials synthesized with toluene. In addition, the positive value obtained from the parameter ΔadsH° shows that the adsorption of the quinoline was endothermally directed. Due to the results obtained up to the present moment, it is possible to point out that the interaction process between the template/monomer can be better explained by electrostatic interaction. One hypothesis is that during the formation of the active complex, methacrylic acid transfers H+ to the basic functional group of the quinoline molecule. This action results in a positive electrostatic interaction between the carboxylate anion and the protonated amino group of the quinoline.
The positive result of TΔadsS° indicates an increase in the number of degrees of freedom in the system; this phenomenon is a consequence of desolvation between quinoline and n-heptane after adsorption, where for a process to be spontaneous, the following condition must be fulfilled: ΔadsG° = ΔadsH° − TΔadsS° < 0. In this study, the relation TΔadsS° > ΔadsH° was obtained; it can be said that the increase in the degree of freedom at the adsorbent interface is due to an entropically directed process.
Table 5 shows the obtained parameters. The Kd values with quinoline were higher for the MIPs compared to the NIPs. Similar results were obtained with the values of k. For the MIPs, the values of k were >1; these values indicate that the adsorption of quinoline occurred preferentially. MIPC obtained the highest k values for both interferences (carbazole and DBT), showing that the synthesis conditions utilized favored higher specificity between the monomer and the template molecule; this indicates that the synthesis with chloroform provides a greater interaction with quinoline.
Fig. 7 illustrates the interactions of different adsorbates with MAA using theoretical calculations. Initially, it was not possible to distinguish the interactions between ion pairs and hydrogen bonds because all calculations converged to this last scenario. However, the results for the relative free energy of interaction provide interesting evidence. The values of ΔΔG, shown in Table 6, confirm that the adsorption selectivity (Fig. 6) is directly proportional to the free energy interaction, which increases in the order A > C > B. In addition, the volume factor is also relevant and directly proportional to the adsorption selectivity, maintaining the order A > C > B. Furthermore, according to our theoretical results, the differences in volume units between quinoline and its analogues are greater in chloroform than in toluene. For example, DBT is 273.95 bohr mol−3 greater and carbazole is 106.85 bohr mol−3 greater than quinoline in chloroform, but is only 266.02 bohr mol−3 and 4.55 bohr mol−3 greater in toluene, respectively.
Using topological electron density analysis, it is possible to define the covalent or ionic character of the interactions between atoms. In order to understand the specific interaction X⋯MAA, we will analyze the critical points (BCP) 1 and 2 appearing between X and MAA (Fig. 7). In this sense, it is important to infer two magnitudes: ∇2ρ and |G/V|. ∇2ρ indicates whether the electronic density converges (negative values) or diverges (positive values) in the BCP. Therefore, covalent interactions present negative values for this magnitude, while electrostatic interactions have positive values.35–37 In |G/V|, the proportion between the kinetic energy (G) and the potential (V) is given. If |G/V| > 1, in BCP, V has a smaller modulus than G. Thus, the electronic density tends to position itself around the nuclei, which is typical behavior of electrostatic interactions. Meanwhile, for covalent interactions, |G/V| < 1.36–38
Table 7 shows the values of |G/V| and ∇2ρ for the BCPs 1 and 2 in Fig. 7. Note that for all cases, ∇2ρ is positive, which characterizes electrostatic interactions. The values of |G/V| converge to this same conclusion, except for the critical point 1 for dimer A. At this point, |G/V| is 0.754 in chloroform and 0.769 in toluene; this proportion is usually found for hydrogen bonds.38 However, BCPs in hydrogen bonds with covalent character must have ∇2ρ < 0, which is not observed in this case.37 Therefore, the interactions between X and MAA should be predominantly electrostatic, which may explain the selectivity difference.
BCP | G | V | |G/V| | ∇2ρ | ||||
---|---|---|---|---|---|---|---|---|
Clo | Tol | Clo | Tol | Clo | Tol | Clo | Tol | |
a Values are truncated at the third decimal place. | ||||||||
1-A | 0.039 | 0.037 | −0.051 | −0.049 | 0.754 | 0.769 | 0.105 | 0.104 |
2-A | 0.008 | 0.009 | −0.007 | −0.007 | 1.201 | 1.206 | 0.039 | 0.041 |
1-B | 0.004 | 0.005 | −0.003 | −0.003 | 1.446 | 1.422 | 0.021 | 0.024 |
2-B | 0.013 | 0.013 | −0.010 | −0.010 | 1.255 | 1.256 | 0.062 | 0.062 |
1-C | 0.022 | 0.017 | −0.019 | −0.014 | 1.154 | 1.215 | 0.099 | 0.081 |
2-C | 0.003 | 0.001 | −0.002 | −0.001 | 1.314 | 1.477 | 0.013 | 0.005 |
Because it is an electrostatic phenomenon, the interaction between X and MAA will compete for interactions with the solvent. In fact, solvents with high polarity will weaken the electrostatic interactions between X and MAA through new electrostatic interactions between solvent⋯X and solvent⋯MAA; thus, the MAA carboxyl groups will be organized outside the pores to maximize the solvent⋯MAA interactions.38–40 For this reason, the pores of MIPC will have less availability to form hydrogen bonds than the pores of MIPT. In consequence, only the strongest X⋯MAA interactions will prevail, giving more selectivity to the MIP. For our DFT calculations, the hydrogen bond interaction between quinoline⋯MAA was by far the strongest among the tested analogs. Thus, quinoline is preferred over DBT and carbazole.
After five adsorption/desorption cycles, the adsorption capacities of the MIPs were retained. For MIPT, the percentage of regeneration of 70% in the first cycle increased as the number of cycles increased, ending with regeneration of 82% in cycle 5. In the MIPC, the percentage decreased from 100% to 83% in cycle 2; however, as the number of cycles increased, the regeneration percentage also increased, with 100% recovery of the material at the end. These data show that the material can be used for several cycles without losing its capacity; in addition, the quinoline can be recovered after the process using a solid–liquid extraction system. This possibility guarantees greater sustainability for the process; thus, the application of this material is more attractive.
The computational calculations indicate that electrostatic interactions are predominant between the monomer and the template. The produced adsorbent materials presented high thermal stability, which enables their use in petrochemical industries. MIPC demonstrated the capability to regenerate itself and to be reused with good efficiency in different cycles, which facilitates the recovery of quinoline after the process.
This journal is © The Royal Society of Chemistry 2018 |