Abstract
Nickel has several industrial uses and is a valuable metal, making its selective separation and recycling a priority goal. A novel adsorbent, a Schiff base organically modified silica (ORMOSIL) aerogel, was prepared, for selective nickel removal from wastewater with other metal ions, by including a salen ionophore in the silica-based network. The newly developed adsorbent takes advantage of the salen’s selectivity and of the high porosity of silica aerogels. The aerogel-like adsorbent was prepared via sol-gel chemistry, using a coprecursor approach and ambient pressure drying. The inclusion of the Schiff base in the silica network was accomplished by reacting an amine-containing silica precursor with an aldehyde and confirmed by nuclear magnetic resonance (NMR) analysis. The adsorbent shrunk only 10% after evaporative drying, which resulted in a highly porous material (85% porosity, 4 cm3 g−1 specific pore volume). The low surface area of 28 m2 g-1 was due to the predominantly macroporous structure of the material (mean pore diameter of 563 nm). Adsorption isotherms and kinetic curves with single and binary mixtures of cations at room temperature were used to assess the selectivity of the adsorbent. The adsorption follows a BET (Brunauer-Emmett-Teller) trend. Due to the proximity of the oxygen and nitrogen atoms in the salen and steric hindrance from their neighboring atoms, it is likely that only the smallest hydrated cations can act as a coordination center and interact with both donor atoms. Thus, nickel was fairly removed (50 mg g-1), while other cations barely interacted with the adsorbent (cadmium adsorption maximum of 5 mg g-1). The estimated selectivity coefficient for nickel ranges from 1.8, in relation to copper, to 9.4 relatively to cadmium, which can be relevant for the separation of nickel in several industrial contexts, for instance, from electroplating sludge.
1. Introduction
Heavy metals (HMs) are pollutants that accumulate in the environment, and due to high emissions from human activities, create severe levels of pollution in water courses, sediments, and soils [1, 2]. They are mobile hence leach and are absorbed by organisms, further spreading through ecosystems [3–7]. This issue can be overcome by stopping the pollutant emissions and/or deploying advanced waste treatment techniques. Moreover, many heavy metals are used as raw materials in different industries [8, 9] and have commercial value [9, 10]. Among them, we can cite copper, lead, cadmium, and nickel. They feature high toxicity and are prone to human-made emissions and accumulation in the environment, and concomitantly, they are relevant resources in today’s societies [8, 11, 12]. Copper is the third most important metal by weight [13]; its yearly production is estimated at 16 million tons [9, 14] and is used in electrical equipment (60% of its consumption), plumbing/construction, machinery, transport, coins, art, music instruments, kitchen utensils [11, 13], and to catalyze many reactions [15–17]. Copper is also a micronutrient in animals and plants. Lead and its related compounds are very toxic [18]. Most recent uses of lead [18] have been banned, such as gasoline additives, in solders, paint pigments, and the EU plans to ban lead ammunition and fishing gear [19]. It is still used for lead-acid batteries and as radiation shields [11]. No biological functions are known [18]. Cadmium is highly toxic. It can be applied as a stabilizer in plastics, pigment in dyes, coating agents, accessories/jewelry, Ni-Cd rechargeable batteries, and solar cells [20, 21], with the latter two contributing to impulse its use [21]. The EU has banned it from being used in plastics, paints, and jewelry [22]. In this paper, we focus our particular interest in nickel due to its increasing market value. Nickel is ferromagnetic at room temperature [11] and is known to catalyze reactions [23]. It is used for the production of metal alloys (stainless steel totals 75% of its consumption [24]), nickel cast iron (jewelry, coins), corrosion-resistant coatings by electroplating, and the production of batteries (nickel-metal hydride and Li-ion) [11, 25, 26]. Nickel has biological functions in plants and cyanobacteria [26]. Due to the aforementioned reasons, nickel is a valuable metal, so the minimization of its pollutant discharges, through separation and recovery, has economic impact on businesses. Furthermore, these also contribute to a more circular and sustainable economy, which many countries are promoting.
Several treatment techniques can be applied to remove heavy metals from wastewaters. These were fully described in a recent review of the authors [1]. However, only a few can be used for the simultaneous goal of recovering the cation being removed from water: ion exchange [27–29], adsorption [30–36], and electroplating [27, 28, 37] are the most effective ones. As such, operating these methods with high efficiency and low-cost at large scale is crucial to promote the recovery of metals as valuable resources. Both adsorption processes rely on the development of high performing materials (the ion exchanger and adsorbent, respectively) that are modified to be selective, while electrolysis can be employed for multiple species by assembling an electrolytic cell and adjusting electric potential on demand. However, a major drawback relies on the fact that electrochemical processes are expensive and only viable in cases where wastewater are not highly contaminated. On the other hand, adsorption processes are simple in operation and design, and their cost is mostly associated with the sorbent materials [27, 28, 37]. Natural materials as well as engineered ones have been used as adsorbents for heavy metals, and the number of published research on the topic has been increasing [1].
Schiff bases are the most common family of molecules used as nickel ionophores [38]. They are considered a subclass of imines and are obtained by a condensation reaction of a primary amine and a carbonyl precursor, e.g., an aldehyde [39]. These molecules can be multidentate and are generally based on nitrogen or oxygen donors but can also include phosphor and sulfur [39]. The formation of the ligand-nickel (II) complex was studied with different Schiff bases by several authors developing colorimetric nickel sensors (in which the response to the cation is measured by absorbance in the ultraviolet-visible region or by fluorescence), and it was found that the stability constant is very high (log
A particular case of Schiff base, salens (symmetrical tetradentate molecules with two nitrogen and two oxygen donor atoms [39]) were studied by Gupta and coauthors [44, 45]. They have successfully used salens as ionophores in ion sensing electrodes with poly(vinyl chloride) membranes [44, 45]. In one of the works, the authors have proposed two distinct molecules as nickel ionophores (Figure 1), both with very high affinity for nickel (log

Representation of the ionophore molecules developed by Gupta et al. [44]. Proposed groups responsible for cation sorption are presented in blue.
In this work, we report the development of a selective adsorbent for nickel, obtained by modifying the silica matrix of a silica aerogel with hemi-salen functional groups. Thus, a highly porous and selective adsorbent is obtained. To the best of our knowledge, it is the first time a silica aerogel is modified using such a functional group. The modified aerogel-like monolith is tested in batch tests to assess its suitability and selectivity against common relevant cations.
2. Materials and Methods
2.1. Materials
As silica sources, methyltriethoxysilane (MTES, ≥99%), tetraethyl orthosilicate (TEOS, 98%), and (3-aminopropyl) trimethoxysilane (APTMS, ≥97%) were purchased from Sigma-Aldrich and used as received. Salicylaldehyde (p.a., Merck) was used to obtain the hemi-salen modifying groups on the silica matrix. Anhydrous oxalic acid (≥99%, Sigma-Aldrich) and ammonium hydroxide (25% NH3 in H2O, Sigma-Aldrich) were the sol-gel chemistry catalysts. Methanol (MeOH, 99.8%, VWR) and ethanol (EtOH, ≥99.8%, Fisher) were used as solvents. Deuterated methanol (CD3OD, 99.8%, Eurisotop) and 3-(trimethylsilyl) propionic acid sodium salt (TMSP, 98%, Eurisotop) were used as solvent and internal reference for NMR spectroscopy. Heavy metal solutions were prepared using copper (II) nitrate hemipentahydrate (p.a., Chem-Lab), lead (II) nitrate (≥99.0%, Sigma-Aldrich), cadmium (II) nitrate tetrahydrate (≥99.0%, Sigma-Aldrich), and nickel (II) nitrate hexahydrate (crystals, Sigma-Aldrich). High purity water was used whenever needed. Nitric acid (65%, Fisher) was used to adjust the pH of solutions.
2.2. Modification of APTMS Co-precursor
The addition of a Schiff base functional group to the silica matrix was achieved by modifying the organic group in a silicon alkoxide precursor. The reaction between the primary amine of APTMS and salicylaldehyde was promoted [46, 47], resulting in a hemi-salen modified silica precursor (HSPTMS, (E)-2-(((3-(trimethoxysilyl)propyl)imino)methyl)phenol), as shown in Figure 2.

Schematic representation of the HSPTMS formation reaction. Proposed groups responsible for cation adsorption are presented in blue, and the groups that take part in the reaction are shown in green.
APTMS and salicylaldehyde were mixed in methanol with a 10 mol% excess of salicylaldehyde. This mixture, labelled as A, was stirred until homogenization and left to react in an oven for 24 hours at 27°C. The mixture turned yellow instantly, indicating the formation of the hemi-salen.
2.3. Synthesis of Schiff Base-Modified ORMOSIL Aerogel-Like Monolith
To prepare silica aerogels with the HSPTMS precursor, MTES was diluted in methanol and hydrolyzed (0.1 M oxalic acid). After one day, the hydrolyzed MTES, mixture B, was mixed with the so-called mixture A and stirred. After a couple of minutes, TEOS and then the basic catalyst (1 M ammonium hydroxide) were added. All sol-gel steps were carried out at 27°C. The Si: solvent: water molar ratios in the sol were kept at 1 : 12 : 8, and the silica co-precursors were mixed in an 50/35/15 molar percentage for MTES, TEOS, and HSPTMS, respectively. Gelation and aging occurred in an oven at 27°C for six days. Aged gels were demolded and washed by being immersed in hot ethanol for a total of six days, with ethanol being changed every day. Finally, the gel monoliths were dried in an oven for three days at 60°C.
2.4. Characterization
The formation of HSPTMS was investigated via NMR (Bruker Biospin GmbH, 1H spectra collected at 400 MHz and 13C spectra collected at 100 MHz). 1H and 13C spectra for both reactants and HSPTMS were obtained by dissolution in deuterated methanol. The analysis of the 13C spectrum did not contribute with additional information relatively to the 1H spectrum analysis, and thus, it is not shown. HSPTMS was obtained directly from mixture A via evaporation of the solvent.
The characterization routines of aerogels are described in previous works by the authors [57, 58]. Bulk and skeletal densities were used to obtain adsorbent porosity. Bulk density of the adsorbents was obtained by weighting the samples and measuring their dimensions on the three axes. Skeletal density was obtained with ground samples by He pycnometry (Accupyc 1330, Micrometrics). Linear shrinkage was calculated from the change in diameter of the dried sample comparatively to the gelation mold. The BET specific surface area was obtained through nitrogen adsorption (ASAP 2000, Micrometrics). Pore volume and average pore size were calculated in accordance with a previous work [59]. Microstructure observation was performed with field-emission scanning electron microscopy (FE-SEM) (Merlin Compact/VPCompact FESEM, Carl Zeiss Microscopy GmbH) [60, 61]. Fourier transform infrared spectroscopy (FTIR) (FT/IR 4200, Jasco) was performed, being the spectra obtained with KBr pellets in the wavenumber range of 4000 to 400 cm−1, with 128 scans and a resolution of 4 cm−1. C, H, and N content of powdered samples was determined by elemental analysis (EA 1108 CHNS-O, Fisons). Zeta potential measurements of an aqueous suspension of the sample nanoparticles (0.1%,
Heavy metal concentration in solutions was determined by flame atomic absorption spectroscopy (AAS) with an acetylene-air flame (939 AAS, Unicam).
To provide insight on the interactions with the different ions, the geometry of hydrolyzed HSPTMS in water was computed. The system was composed of a box, with the precursor molecule surrounded by 632 water molecules, of dimensions
2.5. Adsorption Tests
Adsorption experiments followed the procedures used previously [57, 58, 65]. The aerogel-like monolith was prepared for adsorption by milling and selecting particles within a size range of 75 to 250 μm. Batch adsorption tests were conducted by mixing the powdered adsorbent and the cation solution in a test flask, shaken in a rotational stirrer at speed setting 16 rpm (REAX 20, Heidolph Instruments) with an adsorbent dose of 2 g L−1 at pH 4 and 25°C. When the test ended, the solution was filtered, and the concentration of the filtrate was determined. For comparison, a silica aerogel without the Schiff base modification was tested in the same conditions.
Isotherm studies were performed by varying the adsorbate concentration from 20 to 500 mg L−1 and conducted for 24 h. The selectivity of the Schiff base ORMOSIL aerogel-like monolith was studied in batch kinetic tests, with contact times ranging from 5 minutes to 24 h, with an adsorbate concentration of 100 mg L−1 per cation. Binary mixtures containing copper and mixtures containing nickel plus copper and nickel ions isolated were tested.
The viability of recovering the adsorbed nickel was assessed by a desorption test performed with HCl and HNO3 on nickel loaded particles (from the isotherm studies with a starting concentration of 100 mg L−1). The loaded adsorbent was shaken with a 1 M solution of the desorption agent for three hours.
A minimum of two replicates were obtained for all reported data points.
2.6. Analysis of Adsorption Data
All physical quantities and adsorption model parameters here described are defined with mass units. While discussing binary mixtures, these are expressed in molar units when required. Adsorption capacity (
Freundlich, BET, and Langmuir-Henry models were fitted to the equilibrium data to interpret the metal ion-adsorbent interactions. The pseudo-first order and pseudo-second order models, as well as a modified pseudo-second order model that considers an initial adsorption, were tested to describe the kinetic data. The fits were performed using nonlinear algorithms, and their quality was assessed using Akaike and Bayesian information criteria [66].
The Freundlich model [67], Equation (3), describes the adsorption on heterogeneous surfaces [68]. It has two parameters: the Freundlich constant,
In the kinetic models used, the equilibrium uptake or adsorption capacity,
The intraparticle diffusion and Boyd models were further applied to have a deeper analysis of the adsorption of nickel on the aerogels (a similar analysis was not possible for copper due to its fast adsorption). The intraparticle diffusion model (presented by Weber and Morris) is based on Fick’s 2nd law equation [74, 75] and can be described through Equation (8)
If
3. Results and Discussion
3.1. Schiff Base Reaction
The formation of HSPTMS, via the reaction illustrated in Figure 2, was studied, and the NMR spectra are presented in Figure 3. The formation of HSPTMS can be observed, namely, by the disappearance of the carbonyl group from salicylaldehyde (peak a, δ ~10 ppm) and the amine group in APTMS (peak b, δ ~2.7 ppm) in the product, and the appearance of a new proton associated with the C=N bond (peak a, δ ~8.4 ppm). Protons associated with the methoxy groups and propyl chain did not suffer changes in their chemical environment, while the ones associated with the aromatic ring (peaks b-e in the salicylaldehyde and HSPTMS spectra) suffered some alterations, due to a new chemical environment.

1H NMR spectra for salicylaldehyde, APTMS, and HSPTMS.
3.2. Modified Aerogel-Like Monolith Characterization
The obtained modified silica aerogel-like monolith and its porous morphology are depicted in Figure 4. The physical/structural properties of the adsorbent are summarized in Table 1.

Photograph (a) and micrograph at ×30 k magnification (b) of the Schiff base-modified aerogel.
Physical/structural properties of the functional aerogel-like monolith adsorbent.
The ORMOSIL aerogel-like monolith samples feature a distinctive bright yellow color, due to the presence of the hemi-salen. Despite being dried by evaporation of the solvent, the gel remained a monolith and exhibits a small radial shrinkage, resulting in a bulk density and porosity similar to that of amine modified silica aerogels and contrasting to their xerogel counterparts [58]. The skeletal density of this adsorbent is lower than that of silica, as it is observed with ORMOSIL aerogels. As expected, the modification of the silica matrix with the Lewis base group impacted the structural properties of the adsorbent, as this sample is less porous (and hence, denser) than other ORMOSIL aerogels (without Lewis base groups) or native silica aerogels [58, 59, 78, 79]. For the Schiff base-modified aerogel-like monolith, hydrogen bonding with solvent molecules does not counter the effect of nonhydrolysable apolar groups in the network, resulting in low capillary stresses during solvent evaporation and low shrinkages in the final material. This result can be due to: the presence of the azomethine group (imine) instead of a primary or secondary amine in amine modified gels; and/or the high pore radius. Thiol modified silica xerogels also retained monolithity after evaporative drying, and their properties are similar to the aerogel counterparts [57, 80]. The specific surface area of this adsorbent is small, even when compared to other ORMOSIL aerogels, and is due to the predominant macroporous nature of the material, as indicated by the average pore size in Table 1. Previous works show that, in some situations, the specific surface area does not make a significant impact in adsorption performance [58, 80].
The FE-SEM image depicted in Figure 4(b) confirms the high porosity of the sample, revealing numerous voids with a wide distribution of pore sizes. In fact, many pores smaller than the average pore size estimated (563 nm, Table 1) are observable in this image. The pearl necklace type of microstructure expected from silica aerogels is visible. Additionally, the secondary silica particles are very small and seem to be poorly individualized, surely due to the extended aging period. This microstructure is similar to that of amine modified aerogels recently reported by the authors [58]. The evaluated structural properties allow to conclude that the prepared sample can be considered an aerogel, according to the different proposals for this definition [79, 81–83].
The FTIR spectra and the C, H, and N compositions of the aerogel-like monolith are presented in Figure 5 and Table 2, respectively. The interpretation of the experimental results obtained by elemental analysis relies on theoretical assumptions regarding the condensation of hydrolyzed precursor molecules in the sol, previously detailed by the authors [65].

Infrared spectrum of the adsorbent.
C, H, and N contents of the Schiff base-modified aerogel-like monolith.
The FTIR spectrum of the adsorbent reveals the formation of the silica matrix and the incorporation of the Schiff base in its structure. Bands associated with siloxane bonds are visible at 577, 1099, and 1162sh cm−1; bands associated with the ortho substituted benzene ring are visible at 1464 and 3100 cm−1, and the bonded hydroxyl groups generate a broad band at 3431 cm−1. Some bands are resultant from the overlap of different functional groups that vibrate at similar frequencies. The bands at 465 and 791 cm−1 are the result of overlapping vibrations from siloxane bonds and the benzene ring, while the band at ~1632 cm−1 results from the overlapping of multiple bond stretching of the benzene group, the azomethine bond, and the hydroxyl groups. These bonds limit the conclusions gained from this analysis. The remaining bands are due to the methyl and methylene groups in the MTMS and HSPTMS precursors.
Analyzing the carbon content of the sample, the elemental analysis reveals that hydrolysis is not complete, suggesting an incomplete condensation of one hydroxyl group per Si. This result can be justified by the fact that HSPTMS co-precursor may undergo incomplete hydrolysis, in comparison to the remaining co-precursors, due to its slower reactivity and to it being added just before the base. Furthermore, the large organic group in this co-precursor is likely to hinder the formation of siloxane bridges. The nitrogen content in the sample is superior to the predicted, which can be related with the presence of ammonia residues even after the washing stage, as found before in samples without nitrogen-containing precursors [58]. The experimental value obtained indicates that there are 1.8 mmol of hemi-salen groups per gram of adsorbent or
The addiction of the hemi-salen groups caused the silica surface to feature a positive charge (
3.3. Adsorption Isotherms and Mechanisms
The experimental data and the best isotherm fit are presented in Figure 6. The fit parameters for the studied models are reported in Table 3.

Adsorption isotherms in single-metal solutions.
Fit parameters for the isotherm models and maximum uptake by Schiff base-modified aerogel-like monolith.
(a)The model did not fit to the data.
Nickel removal, cation for which some salens were shown to be ionophores, was the highest observed out of the four cations studied. This was expected, given the high formation constants of nickel-Schiff base complexes. From the statistically point of view, the BET model represented the data the best for the adsorption of nickel, copper, and lead, as the data points seem to follow a type II isotherm (gas isotherm classification) [88] or an L3 isotherm (solid-liquid isotherm classification) [89], which is less defined in nickel adsorption. Thus, the adsorption by the Schiff base-modified ORMOSIL aerogel-like monolith is suggested to occur in multilayers of adsorbate deposited on the surface of the adsorbent. However, understanding how metal ions can form multilayers by adsorption onto an ionophore-containing surface is not straightforward. In fact, the adsorption of metal ions from aqueous media is generally controlled by adsorbent surface chemistry and surface area [90] and occurs in the form of coordination bonds [91–98]. In fact, in other materials modified with Lewis bases, the multiple neighboring functional groups participate in the complexation of one cation, and the complexing amine groups can even replace water molecules from the primary hydration shell (which is known as inner-sphere complexation) [91, 95, 96, 99, 100]. The interaction of Lewis base-modified adsorbents with the heavy metal cations may be explained by the Pearson acid-base concept. Thus, the results reported in this work follow a different trend compared to amine modified silica aerogels [58] or similar materials for which the Langmuir model fitted to the data the best in most situations [52–54]. However, by the way the data is presented in the literature, we cannot conclude if testing higher cation concentrations would produce an isotherm shaped like a BET curve. Even so, it should be highlighted that a multilayer-like adsorption for metal ions has been previously reported [101–104]. Furthermore, in a previous work, the existence of multiple types of surface groups available for adsorption generated a L4 isotherm for copper [58]. Different hypotheses for this observation can be raised and are discussed as follows.
The adsorption isotherms can be due to a combination of different factors: the macroporous nature of the adsorbent [88] and the affinity between adsorbate layers, as relevant as cation-adsorbent interactions [89], assuming that multilayers are present. The latter means that adsorption is not an exclusive result of the interaction with active surface sites. Multilayer adsorption is commonplace in the adsorption of organic pollutants due to van der Waals or π-π interactions; however, cations would repel each other. Some anions may interact with the adsorbed cation layer [105]; based on that, Jorgetto et al. [106] proposed that the BET isotherm reflects multiple adsorbate layers due to the electrostatic interactions between alternating layers of cations and counter-ions in solution. However, this situation does not seem so likely with the nitrate anions used in this work [105]. To evaluate the possibility of multilayer of adsorbed cations at the surface of the adsorbent, energy dispersive X-ray analysis (EDS) was conducted with the loaded adsorbent particles. The presence of cations at the surface was not detected (data not shown), with one exception for copper. In a previous work, SEM and EDS analyses showed that the loaded adsorbent particles were coated by metal salts, possibly due to nucleation induced by adsorbed species or to chemical precipitation at the surface of the adsorbent [80]. Such behavior was not observed in this work, and the microstructure of the loaded adsorbent simply seemed more closed (data not shown). This was expected as the loaded particles were dried by evaporating water, which increases capillary pressure in the pores by forming hydrogen bonds with the solid. Only when analyzing a large sample area could the presence of some copper (as high as
To obtain further insight on this matter, the hybrid Langmuir-Henry model (Equation (5)) was analyzed (Table 3). This considers the possibility of the structure of water molecules inside the adsorbent’s pores to be different than in the liquid water (bulk) allowing the dissolution of different amount of electrolyte, in addition to adsorption at active surface sites. If that is occurring, the hybrid Langmuir-Henry model would fit better to experimental data, which is not the case. Lead is the only cation for which the plateau corresponding to the monolayer adsorption capacity is well-defined, resembling a Langmuir curve, suggesting that the cation-adsorbent interactions are strong in this case, and AIC and BIC do not agree on the preferred model; however, we can observe that the data resembles the BET curve.
Based on the previous observations, we must conclude that in spite of the BET model represented the data the best for nickel, copper, and lead, the adsorption of these cations does not occur in multilayers. As such, the increase in adsorption capacity after a plateau (saturation of active sites) may simply reflect changes in the adsorbent’s surface/active site availability (e.g., by rearranging of adsorbed species) or the occurrence of metal nucleation induced by the salen-metal initial interaction [107–109]. EDS spectra support this conclusion, since some copper was quantified using a large area of sample.
A different but related question that shall be addressed is the interaction mechanism between metal ions and the Schiff base. At the highest equilibrium concentration, the Schiff base: Ni ratio is 2 : 1. It could be said that not all potential adsorption sites are occupied; however, as mentioned before, it was proposed that neighboring Lewis base groups are involved in the complexation of one cation. Hence, it is possible that all Schiff base groups in the aerogel are involved in the adsorption process. This hypothesis is explored again below. Nonspecific interactions, like electrostatic interactions, cannot occur with the positively charged surface of the modified silica aerogel-like monolith. The adsorption must occur due to strong interactions with the adsorbent’s surface. Furthermore, it was observed (data not shown) that loaded adsorbent particles exhibit a different color than pristine particles, which was shown to be associated with metal complexation [91]. Thus, we must conclude that the adsorption is due to complexation of cations by the Schiff base, which most likely occurs with partial loss of hydration water (chemisorption) [87].
We propose that the hemi-salen groups interact differently with the cations, due to the latter’s hydrated radius (or its size relative to the distance between the electron donor atoms in the Schiff base) and the proximity of adjacent Schiff base groups. The imine and the hydroxyl group in the Schiff base are relatively close (distance estimated at 2.4 Å, Figure 7), compared to the cation’s hydrated radius (3.9 Å in the equatorial direction of the Jahn-Teller distorted octahedron for copper, 4.1 Å for nickel, 4.6 Å for cadmium, and 5.1 Å for lead [110]). So, considering a single hemi-salen group, it is possible that the bigger cations are not complexed by both atoms, as they would retain a large radius even after partial loss of the hydration shell, and could be repelled by atoms neighboring the imine and hydroxyl groups. The Schiff base could then act as two different active sites, for cadmium and lead, resulting in a weaker interaction between the phases. This justifies the cadmium isotherm and the reduced removal of cadmium and lead. The Freundlich model fitted to the cadmium adsorption data the best, as a linear trend is observed. This shows that cadmium adsorption is due mainly to pore diffusion rather than interactions between solid and liquid phases. A different hypothesis, already mentioned, considers two adjacent hemi-salen groups complexing the metal together, in what would be a tetradentate salen ligand [46, 47, 49]. It is easily understandable that the distance between these groups will determine the cation size that acts as the coordination center. In this situation, potentially all the cations tested can be adsorbed, but because the distribution of the active sites in the adsorbent is not uniform, the distance between surface groups cannot be determined.

Representation of the optimized structure of the hydrolyzed HSPTMS coprecursor.
The removal observed is very reduced when compared to previous amine-modified aerogels [58] or mercapto-amine aerogels [80] but is higher than unmodified aerogels [58, 65]. The cumulative effect of the ionophore ligand, for which cation radius is important, and the impossibility of electrostatic interactions help to justify this observation. A silica aerogel not modified with Lewis bases is not a good heavy metal sorbent and only interacts with lead (the most labile) under the same conditions (uptake for copper of 2.5 mg g-1, 15.7 mg g-1 for lead, 1.6 mg g-1 for cadmium, and 1.4 mg g-1 for nickel). The here presented Schiff base ORMOSIL aerogel-like monolith is much more efficient in nickel removal than the salen modified silica gel (uptake capacity of 4 mg g-1) reported by Kursunlu et al. [53]. When modifying mesoporous silica with salen, Enache and coauthors found that the adsorption of lead increased a little (20% for MCM and up to 5% for HMS) up to 99 mg g-1 [52]. The adsorption of lead increased 49% with the incorporation of the salen reported in this work. The nickel uptake reported in this work is also superior to that of some amine-modified silicas (28 mg g-1) [111] but is inferior to that of EDTA-modified silica (74 mg g-1) [112] or poly (sodium 4-styrenesulfonate)-modified zeolitic imidazolate framework-8 (330 mg g-1) [113]. The Schiff base silica aerogel-like monolith is a better nickel sorbent than natural minerals, zeolites, sawdust and biochar adsorbents [1, 114, 115], and other materials [116, 117].
3.4. Adsorption Kinetics and Selectivity
The HM adsorption by the prepared adsorbent was studied in binary mixtures, as a kinetic test. The total molar uptake results are presented in Table 4 and Figure 8. The kinetic curves for copper and nickel adsorption isolated and in binary mixtures are reported in Table 5 and Figure 9.
Fitting parameters of the kinetic models for the total adsorption in binary mixtures.
(a)Negligible value.

Adsorption kinetics for the total uptake in binary mixtures of (a) nickel and (b) copper.
Fitting parameters of the kinetic models for the adsorption of nickel and copper isolated and in the presence of interfering ions.
(a)Negligible value. (b)System is always at equilibrium—Figure 8(b).

Adsorption kinetics for (a) nickel and (b) copper, isolated and in binary mixtures.
The plots from Figures 8 and 9 illustrate that a significant amount of adsorbates is adsorbed in a short period of time. From there on, the adsorption process continues slowly, which is more pronounced when nickel is present. On the majority of datasets, the modified pseudo-second order kinetic model is reduced to pseudo-second order model, as there is no relevant initial adsorption. The fit quality criteria, AIC and BIC, agree that this model is clearly better than the pseudo-first order model, which indicates that the adsorption of the different sorbates is controlled by the surface reaction, i.e., the interaction with active surface sites. When an initial uptake is significant, AIC penalizes the modified pseudo-second order model for having more parameters (Tables 4 and 5) while BIC considers this model more adequate than the pseudo-first order one. Figures 8 and 9 clearly show that the kinetic data follows the pseudo-second order trend, and so this model was chosen as the most adequate. Figure 6 elucidates that, on the concentrations tested for the kinetic tests, we are near the BET curve knee, and thus, the insight provided by the kinetic tests does not allow to disprove the existence of multilayers since the BET theory assumes that the first layer is Langmuir-like, and so it is reduced to the Langmuir equation when there is only one adsorbed layer [70, 118].
Regarding the adsorption of nickel by the Schiff base-modified aerogel-like adsorbent, the analysis of Boyd model, considered for short adsorption times, shows that intraparticle diffusion is not the rate limiting step (Table 6). This was expected as the pseudo-second order fitted to the data, indicating that chemisorption is the limiting step. It should be noted that the pseudo-first order model, which was shown to be mathematically equivalent to the Boyd model derived for long adsorption times [119], did not fit well to the data. It can also be seen that, apart from nickel adsorption kinetics in the presence of lead, the adsorption process shows two linear regions, as described by Equation (9). The first stage involves a diffusion effect of the boundary layer, followed by a faster process which might be related with the interaction with the ionophore. For the case of Ni(II), in the presence of Pb(II), only one step is detected. Taking into account the magnitude of
Kinetic parameters obtained for intraparticle diffusion and Boyd models for adsorption of Ni(II) in the absence and presence of other metal ions.
When copper and nickel are mixed, the adsorption of copper (Figure 9(b)) is almost instantaneous and remains somewhat constant, while nickel adsorption increases over time and seems that it has not reached equilibrium after 24 hours (Figure 9(a)). Figure 8 shows that the mixtures of cadmium and lead with nickel have a higher uptake than those with copper. Figure 9 clarifies this observation by revealing that nickel adsorption is not affected by the presence of these interfering cations (considering the equilibrium uptakes and their experimental error in Table 5); hence, the equilibrium found has no statistically significant change for all situations presented, and its uptake is higher than that of copper. With copper (Figure 9(b)), interfering cations have a synergistic effect on the former’s adsorption, enhancing it, and its equilibrium uptake is affected by the interferents.
The observation that nickel adsorption is unaffected by the presence of interferents is a clear indication of selectivity of the material towards this cation. Furthermore, this selectivity could not be achieved by mere nonspecific interactions between phases, meaning that the hemi-salen ligands in the aerogel-like monolith are strongly binding this cation, corroborating the idea of strong interactions for some cations. This is most likely the reason for the different shape of the nickel isotherm.
The recovery of nickel, achieved with the desorption tests with both desorption agents (HCl and HNO3), was of 7%, meaning that the recovery of the adsorbed cation is not likely to be achieved without the degradation of the adsorbent as a consequence of strong adsorbent-adsorbate interactions. Thus, this is also an indication of the strong adsorption interaction between the cation and the adsorbent, in agreement with the chemisorption hypothesis drawn in the previous discussion.
The selectivity coefficients in the binary mixtures are reported in Table 7. The selectivity coefficients confirm that nickel is adsorbed preferentially, in the presence of the other tested cations, as expected. In the absence of nickel, a selectivity towards copper is verified as expected from salen ligands.
Selectivity coefficients of the Schiff base-modified aerogel-like monolith.
4. Conclusion
The preparation of a selective adsorbent for nickel was achieved by incorporating a salen ionophore in a silica network, resulting in a bright yellow adsorbent. The ORMOSIL aerogel-like monolith was dried via ambient pressure drying but remained a highly porous (85%) monolith, with a reduced linear shrinkage (10%). The small specific surface area obtained (28 m2 g-1) is typical for some modified silica aerogels and is a result of the macroporous structure of the material. The isotherm studies revealed that the BET curve models the adsorption of the majority of the tested cations, although different interactions with surface sites at higher sorbate concentrations seem more likely than the formation of multilayers. Nickel was the most adsorbed cation (experimental maximum of 50 mg g-1), an improvement towards different reported materials. Kinetic studies evidenced a very fast initial uptake, which is sometimes reflected by what can be considered as an initial boost to the adsorption process that slows down as active sites become occupied. The adsorption of binary mixtures, not common in newly reported adsorbents, revealed a high affinity for both copper and nickel by the Schiff base-modified aerogel-like monolith. However, it is clear that copper is only the preferred sorbate in the absence of nickel, as copper uptake decreases slightly (in the presence of nickel) as time increases, and nickel uptake is not affected by the presence of other species in solution. The selectivity coefficients calculated indicate that the prepared adsorbent removes nickel ~9 times more than cadmium, ~5 times more than lead, and~ 2 times more than copper, confirming the fitness of this material for the separation of nickel from solutions with other metals. Further interfering ions should be tested in the future to assess the viability of the material here presented in separating nickel from complex mixtures/mixtures containing other common ions. The desorption of nickel was not very extensive, so further tests must be conducted in different conditions for the improvement of nickel release. Because of the high stability of Schiff base-nickel complex that was formed, desorption might only be significant if ligands are used as desorption agents, and several ligands should be tested for this purpose. The economic value of nickel may trigger the interest in the new sorbent here developed, in particular due to its high selectivity, which can contribute to follow the circular economy concept. Additionally, these results pave the way for using silica-based aerogels as potential materials for acting as sensors of metal ions.
Footnotes
Data Availability
Data and material may be available upon request to the corresponding authors.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this article.
Acknowledgments
The authors acknowledge Pedro Maximiano for his cooperation with the computational calculus. J.P. Vareda acknowledges the PhD grant SFRH/BD/131280/2017 by Fundação para a Ciência e a Tecnologia, I.P. (FCT, Portugal) funded by national funds from MCTES (Ministério da Ciência, Tecnologia e Ensino Superior) and, when appropriate, cofunded by the European Commission through the European Social Fund. Consumables for the syntheses and characterizations performed at CIEPQPF and CQC research units were funded by national funds through the FCT and when appropriate cofunded by FEDER under the PT2020 Partnership Agreement under the projects POCI-01-0145-FEDER-006910 and POCI-01-0145-FEDER-007630 (FCT Refs. UIDB/EQU/00102/2020 and UIDB/QUI/00313/2020, respectively). Moreover, tests were also performed under the scope of the project “BIOSHELL” (Ref. “BLUEBIO/0003/2019–Recycling Crustaceans Shell Wastes for Developing Biodegradable Wastewater Cleaning Composites”), financed by FCT within program ERA-NET cofund on Blue Bioeconomy (BlueBio)—Unlocking the Potential of Aquatic Bioresources.
