Sushil Kumar‡
ae,
Mahira Bashri‡a,
Safa Gaber
a,
José I. Martínez
b,
Matthew J. O’Connorc,
Sabu Varghesec,
Blaž Belecd,
Gisha Elizabeth Luckachana and
Dinesh Shetty
*ae
aDepartment of Chemistry, Khalifa University of Science & Technology, Abu Dhabi, 127788, United Arab Emirates. E-mail: dinesh.shetty@ku.ac.ae
bInstituto de Ciencia de Materiales de Madrid (ICMM-CSIC), C/Sor Juana Inés de la Cruz 3, 28049, Madrid, Spain
cCore Technology Platform, New York University Abu Dhabi (NYUAD), Saadiyat Island, 129188, Abu Dhabi, United Arab Emirates
dMaterials Research Laboratory, University of Nova Gorica, Vipavska 11c, Slovenia
eCenter for Catalysis & Separations (CeCaS), Khalifa University of Science & Technology, Abu Dhabi, 127788, United Arab Emirates
First published on 21st July 2025
Efficient and selective gold recycling from e-waste mixtures is crucial for advancing sustainability, the circular economy, and waste management objectives. However, the complex composition of e-waste and the underdeveloped practical potential of current adsorbents highlight the need for innovative sorbent development. Here, we present quasi-planar, benzoxazine-core-based two-dimensional covalent organic frameworks (COFs) with synergistic hetero-atom control across the network, achieved via a multi-component synthetic approach. These COFs, designed with specific gold interaction sites, were explored for the highly selective recovery of gold from real e-waste mixtures. By tuning the hetero-atom distribution, we achieved an impressive gold uptake capacity of 3467 mg g−1 and over 92% selectivity in presence of 12 competing metal ions from e-waste leachate. These findings underscore the critical role of molecular-level engineering in COF sorbents, paving the way for the practical recovery of noble metals from e-waste.
New conceptsSelective and efficient gold reclamation from e-waste mixtures is paramount for fostering sustainability, circular economy principles, and optimized waste management strategies. The heterogeneous nature of e-waste and the limited efficacy of existing adsorbents underscore the demand for advanced sorbent innovation. To address this, we introduce two new quasi-planar benzoxazine-core-derived two-dimensional covalent organic frameworks (COFs) with integrated heteroatoms, crafted through a multi-component synthetic methodology. Tailored with precise gold-binding domains, these COFs—namely PGBpy and PGBD—incorporate Namine, NBpy, and Oether linkages, enabling highly selective gold extraction from authentic e-waste samples. This strategic design yielded remarkable gold uptake capacities of 3467 mg g−1 (PGBpy) and 2590 mg g−1 (PGBD), driven by the critical influence of embedded Namine and Oether heteroatoms within the COF structure, which preferentially sequester gold ions amidst competing metal ions in acidic leachates. Both COFs exhibited over 92% gold selectivity, retaining stability and porosity, with DFT analyses confirming the pivotal role of heteroatoms in the benzoxazine core and auxiliary Bpy unit in amplifying adsorption selectivity and efficiency. |
Covalent organic frameworks (COFs) are organized structures of organic units linked by covalent bonds,21 offering great adsorption potential due to their intrinsic porosity, exceptional durability, and regenerative nature.13,22 COFs synthesis allows precise control over pore structure and chemical characteristics through reticular and dynamic covalent chemistry.21 However, few studies explore COFs for gold recovery from e-waste, where the pore interior functionalized with groups such as –SH, –NH2, and –OH to enhance molecular interactions.23,24 Multicomponent synthesis is a useful strategy where we can introduce multiple metal coordination sites decorated pore interiors to improve interactions with the adsorbed species.25
The role of the quasi-planar core containing COF skeleton in a heteroatomic pore environment is an interesting avenue of research. Benzoxazine-containing polymers exhibit high thermal stability,26,27 excellent electronic properties, and strong chemical stability.28,29 While one-dimensional COF featuring quasi-conjugated chains with a benzoxazine-core has recently been reported.30 To the best of our knowledge, a quasi-planar benzoxazine-core-based two-dimensional (2D) COFs have not yet been reported in the literature. We envision that 2D benzoxazine-core COFs, combined with a synergistically heteroatom-controlled backbone and enhanced stability, could be promising candidates for recovering precious metals from e-waste. Considering these facts, we synthesized two such COFs i.e. PGBpy and PGBD, with active gold-interaction pockets having Namine, NBpy, and Oether linkages via multi-component synthesis. The rational design resulted in efficient recovery of gold from e-waste, contributing to a high uptake value of 3467 mg g−1 for gold ions with PGBpy and 2590 mg g−1 with PGBD. The obtained results highlight the importance of inherited heteroatoms i.e. Namine atoms and Oether atoms within the COF backbone in facilitating the selective capture of gold ions in the presence of other competing metal ions present in e-waste leach solution. Both COFs displayed >92% selectivity for gold under acidic conditions, even in real e-waste, without compromising its stability and porosity. DFT studies corroborated the crucial role of heteroatoms in the benzoxazine-core and supplementary Bpy unit in enhancing selectivity and adsorption efficiency.
The crystallinity of PGBpy and PGBD were elucidated using powder X-ray diffraction (PXRD) analysis (Fig. 1c). PXRD pattern displayed a first peak at 2θ = 5.78° (PGBpy) and 5.92° (PGBD), corresponding to diffraction from 110 facets, with a second peak at 2θ = 14° for diffraction from 101 facets. However, in the case of PGBD, the second peak is broad and, thus, difficult to locate precisely. For both COFs, a broad peak in the range 2θ = 18–25° corresponds to diffraction from 002 facets. Three stacking models were computationally simulated to delineate the framework structure of PGBpy and PGBD, and the corresponding model's PXRD patterns were generated and compared with the experimental PXRD pattern (Fig. S3, ESI†). According to our analysis, the theoretical PXRD patterns for the eclipsed AA and staggered AB stacking models exhibit their first diffraction peaks at approximately 2θ = 3.56°, which significantly deviate from the experimental first peak positions observed for PGBpy and PGBD. This mismatch indicates that neither the AA nor AB stacking models accurately represent the structural arrangement of our COFs. Attempting Pawley refinement against these models would indeed yield high R and Rwp values, rendering such refinements uninformative due to the poor fit between theoretical and experimental data. In contrast, the PXRD pattern simulated for the ABC stacking model matches well with experimental PXRD pattern, as illustrated in Fig. 1c. This suggests that the ABC stacking configuration better approximates the actual framework structure. Pawley refinement performed for both COFs resulted in the following fitting parameters: for PGBpy: Rwp = 3.98% and Rp = 4.22%, and for PGBD: Rwp = 2.94% and Rp = 2.17%. Following Pawley refinement, the unit cell dimensions were derived as follows: for PGBpy: a = 30.10 Å, b = 30.04 Å, and c = 8.15 Å; α = 92.57°, β = 92.44°, and γ = 120.67°; and for PGBD: a = 30.18 Å, b = 29.67 Å, and c = 8.08 Å; α = 81.43°, β = 84.38°, and γ = 119.79° (Sections S5 and S6, ESI†). Notably, the closeness in magnitude between parameters a and b, alongside the proximity of the angle between a and b to 120°, suggests a two-dimensional hexagonal lattice arrangement within each layer. This observation suggests the formation of an ordered hexagonal porous structure within PGBpy and PGBD. The high value of Rwp (w/o bck) observed in the PXRD pattern generated after performing the Pawley refinement for the ABC model indicates moderate crystallinity (Fig. S4 and S5, ESI†). Nevertheless, we would like to emphasize that the primary goal of our design is not to achieve maximal crystallinity or porosity but rather to exploit the unique properties of the benzoxazine core—namely its flexibility and heteroatom-rich composition—for enhanced metal ion adsorption.
To determine the morphology of these COFs, scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images were collected. The SEM images of PGBpy and PGBD show spherical morphology with average sizes of 2.5 and 2.0 μm (Fig. 1d and e, respectively). TEM images further confirm the spherical morphology (Fig. 1f and g). In the case of PGBD, some sheet-like structure is attached to the sphere surface, as observed in its TEM (Fig. 1g). These sheets were most probably generated via exfoliation of the surface of the sphere or breakdown of the small sphere during the sonication process employed in sample preparation for SEM/TEM studies. A similar observation is also noticed in PGBpy, where small spheres and sheets are present on a very small scale (Fig. S12, Section S7, ESI†). The chemical structure was further surveyed by XPS analysis (Section S8, ESI†). Both COFs displayed three intense peaks in the spectra, representing C 1s, N 1s, and O 1s signals. In the deconvoluted XPS profile of C 1s, the two signals present at 283.28 and 284.26 eV (for PGBpy), 283.38 and 284.47 eV (for PGBD) are assigned to CC and C–C bonds, respectively, present in the framework (Section S8, ESI†). The broad signal at 285.64 eV (for PGBpy) and 286.34 eV (for PGBD) was assigned to the C–O or C–N bond. The high-resolution N 1s spectra displayed two signals at 398.77 eV (for PGBpy) and 398.82 eV (PGBD), corresponding to the C–N bond formed during the condensation reaction (Section S8, ESI†). In contrast, the signal at 397.81 eV corresponded to the C
N bond of the Bpy unit. Moreover, in the high-resolution spectra of O 1s, the signals present at 531.80 and 533.18 eV (PGBpy); and 530.84 and 532.71 eV (PGBD) are attributed to the C–O–C and adventitious C–OH bond, respectively (Section S8, ESI†). Based on the above analysis, we conclude the formation of the benzoxazine core in the PGBpy and PGBD, along with some unreacted hydroxyl groups within the framework structures.
The porous nature of PGBpy and PGBD, crucial for metal ion adsorption from e-waste was analyzed via nitrogen gas sorption isotherm, revealing a type II isotherm. The BET surface area (SBET) values were 1557 m2 g−1 (PGBpy) and 931 m2 g−1 (PGBD), with mesoporous pore sizes of 28–58 Å (PGBpy) and 28–60 Å (PGBD) (Section S9, ESI†). The pore sizes align with the simulated COF, but larger-than-expected values suggest structural defects or partial layer offset, possibly due to the benzoxazine core's flexibility, which may disrupt ideal ABC stacking and from mesoporous channel.31 The BET isotherms and PXRD patterns, showing a decent surface area and broader peaks respectively, suggest reduced crystallinity and deviate from the conventional advantages of planar COFs, where strong interlayer π–π stacking typically enhances porosity and structural order. The chemical stability of both COFs was investigated under acidic and basic conditions. FTIR spectra of PGBpy after 3 M hydrochloric acid treatment showed no significant shift in the characteristic stretching bands (Section S11, ESI†), while treatment with 3 M NaOH led to the disappearance of characteristic CC and C–N peaks, suggesting framework disintegration. Similar observation was noticed in PXRD studies (Fig. S21, ESI†). It is important to note that the intended application of the developed COFs is to selectively recover gold from e-waste leachate, which is generally acidic in nature, and, therefore, the lack of stability in the basic conditions is less of a concern.
The equilibrium gold adsorption capacities (Qe) of benzoxazine-core based-COFs were obtained by adsorbing aqueous AuCl3 solutions with varied initial concentrations (Ci) ranging from 1 to ∼2000 mg L−1 (Fig. 2b). The adsorption capacities obtained were fitted with the Langmuir and Freundlich adsorption isotherm models (Fig. S27 and S28, ESI†). The adsorption isotherms exhibited a strong correlation with the Langmuir model (R2 >0.998, for PGBpy; and >0.995, for PGBD) compared to the Freundlich model (R2 >0.876, for PGBpy; and >0.914, for PGBD), suggesting the strong monolayer adsorption of Au ions over COF skeleton. In our study, the adsorption isotherms for Au (likely as [AuCl4]− ions before reduction) strongly fit the Langmuir model (R2 >0.998 for PGBpy; >0.995 for PGBD), indicating uniform, single-layer adsorption on specific binding sites within the COF. This behavior is attributed to the heteroatom-rich pore environment, where protonated amine (–NH+) and oxygen (–O– or –OH) sites under acidic conditions selectively bind [AuCl4]− via electrostatic interactions and coordination bonding, as observed in adsorption experiments. The maximum adsorption capacity (Qm) of PGBpy and PGBD was found to be 3467 mg g−1 and 2590 mg g−1, respectively. The large Qm value observed for PGBpy signifies the crucial role played by secondary coordination sites present in the COF backbone, inherited from Bpy unit in addition to the benzoxazine core. On the other hand, this was not the case for PGBD, due to the absence of an additional coordination site, which completely contributes to the heteroatom skeleton in the backbone. This controlled study has helped us to evaluate the role of additional nitrogen atoms present within the PGBpy in enhanced gold adsorption performance.
Since practical gold capture applications are often performed under acidic conditions, the adsorbent chemical stability while maintaining high uptake efficiencies becomes crucial. Considering this fact, we investigated the adsorption performance of COFs across a range of pH conditions, specifically between pH 2 and 12 (Fig. 2c). The findings demonstrated ∼100% gold adsorption within the acidic pH range (2 to 6), where Au3+ primarily exists as planar [AuCl4]− ions within the solution. Under acidic conditions, the protonated Namine and Oether functional groups within the COFs can efficiently adsorb [AuCl4]− ions through electrostatic interactions.32
The selective uptake of gold ions over other competitive metal ions within the complex leaching solution is crucial for a real-life application. In this direction, for both COFs, we investigated selectivity by mimicking the metal composition of the real waste solution, including Au3+, Cu2+, Cr3+, Al3+, Zn2+, Co2+, and Cd2+ ions. According to our findings, PGBpy displayed excellent selectivity where the Au ions had a maximum uptake of 96.3% followed by the Cu2+ (6.2%), Cr3+ (4.2%), Al3+ (3.4%), Zn2+ (3.1%), Co2+ (2.1%), and Cd2+ (1.2%) (Fig. 2d). On the other hand, PGBD exhibited an uptake efficiency of 92.7% compared to Cu2+ (5.7%), Cr3+ (9.9%), Al3+ (2.2%), Zn2+ (3.7%), Co2+ (4.7%), and Cd2+ (2.0%). These observations suggest that the additional coordination sites (NBpy) in PGBpy played a crucial role in improving adsorption efficiency and selectivity. Meanwhile, PGBpy and PGBD demonstrated excellent distribution coefficient (Kd) values for gold ions compared to all the competing ions in the simulated e-waste (Table S9 and S11, ESI†). This outstanding performance can be primarily ascribed to the strategic incorporation of the benzoxazine core, which under acidic conditions (e.g., pH 4–6), promotes the protonation of the amine (–N → –NH+) and partial protonation or hydrogen-bonding activation of oxygen sites which ultimately enhance the electrostatic attraction and coordination capability toward the anionic [AuCl4]−. This dual-protonation effect creates a highly polarized binding environment, amplifying the COF's affinity for gold ions beyond what single-site interactions typically achieve.
For the continuous adsorption of Au ions, the efficiency of the adsorbent for consecutive adsorption/desorption studies is critical.33,34 In this direction, gold-adsorbed COFs were treated with the stripping solution (1 M thiourea and 1 M HCl) to regenerate and reuse the COFs for the next cycle of gold recovery. According to our investigation, both PGBpy and PGBD exhibit excellent gold uptake efficiency of >98% for ten and five consecutive adsorption/desorption cycles respectively (Fig. 2e and Fig. S29, ESI†). Meanwhile, there were no significant changes in the FTIR spectra of regenerated COFs compared to the pristine, suggesting the structural integrity remained intact even after multiple adsorption/desorption cycles (Fig. S31, ESI†).
Furthermore, energy-dispersive X-ray (EDX) mapping performed on Au@PGBpy and Au@PGBD showed the uniform distribution of gold nanoparticles within the COF matrix (Fig. S34 and S35, Section S13, ESI†). After gold adsorption studies, the PXRD pattern of gold adsorbed COFs i.e. Au@COFs were recorded. In the PXRD pattern, the characteristic diffraction peaks of gold nanoparticles in Au@PGBpy and Au@PGBD samples dominate our as-synthesized COFs characteristic peaks, while simultaneously confirming the gold uptake (Fig. S24, ESI†). Whereas, the diffraction peaks present at 2θ = 38.4° and 44.4° correspond to the diffraction from the 111 and 200 planes of Au0 nanoparticle, respectively. To understand the interactions of gold nanoparticles with COFs, a full survey XPS profile analysis of Au@COFs was performed and compared with pristine samples (Fig. S15 and S16, Section S8, ESI†). For Au@PGBpy, the deconvolution of the high-resolution C 1s profile shows two discernible signals at 283.42 and 284.55 eV, signifying the CC and C–C bonds of the PGBpy framework. In contrast, the broader signal appeared at 286.42 eV, which corresponds to the C–O/C–N bond. An observation of a slight shift in the binding energy value by ∼0.78 eV for these C–O/C–N bonds to the higher value compared to pristine PGBpy indicates benzoxazine core interactions with [AuCl4]−. High-resolution N 1s spectra show three distinct peaks at 397.76, 398.78, and 399.68 eV, corresponding to C
N (Bpy), C–N, and Au–N, respectively, suggesting that the nitrogen sites present in the benzoxazine core and Bpy unit of the framework were participated for the Au adsorption (Fig. 3a). Two signals at 529.60 and 531.30 eV in the high-resolution O 1s spectra were attributed to C–O⋯Au and adventitious C–OH⋯Au bonds, respectively (Section S8, ESI†). Interestingly, compared to O 1s profile of pristine PGBpy, a significant shift of ∼1.2–1.4 eV was observed for C–O–C and C–OH bond peaks. This observation suggests the presence of a strong interaction between the Ohydroxyl and Oether with the absorbed [AuCl4]− ions.
In a high-resolution Au 4f XPS profile, peaks observed at 83.21 and 87.48 eV were assigned to the Au0 4f7/2 and Au0 4f5/2 states, respectively (Fig. 3c). Additionally, peaks at 83.99 and 89.37 eV correspond to the Au1+ 4f7/2 and Au1+ 4f5/2 states, while the peaks located at 85.43 and 91.39 eV are attributed to the Au3+ 4f7/2 and Au3+ 4f5/2 states. Au XPS profile analysis confirms the presence of Au3+ and its reduced form, Au1+ and Au0, in the Au@PGBpy material. Notably, the high intensity of the Au0 peaks, relative to those of Au3+ and Au1+ peaks, suggests that the reduction of Au3+ to Au1+ and Au0 occurs instinctively upon the interaction of [AuCl4]− with the benzoxazine core and Bpy unit of the framework. A similar analysis of Au@PGBD further validated the role of the benzoxazine core in providing binding sites for gold ions. A high-resolution C 1s profile displayed three peaks at 283.29, 284.51, and 286.94 eV corresponding to the CC, C–C, and C–O/C–N bonds, respectively. Compared with the C 1s profile of PGBD, a slight shift of ∼0.6 eV in the binding energy of the C–O bond of Au@PGBD suggests that the O atom of the benzoxazine core is involved in interaction with gold ions. Meanwhile, a new peak at the binding energy value of 399.58 eV corresponds to Au–Namine interaction as observed in the high resolution of N 1s for the Au@PGBD, compared with PGBD (Fig. 3b). The high-resolution O 1s profile displays two peaks at 530.64 and 532.49 eV, respectively. The peak at 530.64 eV corresponds to the binding energy of Au–Oether interactions. A shift of ∼0.2 eV for the C–Oether bond w.r.t. the O 1s profile of pristine PGBD suggests the involvement of the O atom of the benzoxazine core in gold adsorption. Furthermore, the 82.99 and 88.24 eV peaks were assigned to Au0 4f7/2 and Au0 4f5/2 states (Fig. 3d). On the other hand, the peaks present at 84.23 and 89.37 eV correspond to Au1+ 4f7/2 and Au1+ 4f5/2 states, while the peaks located at 86.70 and 90.03 eV correspond to Au3+ 4f7/2 and Au3+ 4f5/2 states. Au XPS profile suggests the majority of the gold absorbed from the feed solution is present in Au3+ form compared to its reduced form Au1+ and Au0 in the Au@PGBD. In summary, we can conclude from XPS analysis that in the case of PGBpy, the Bpy serves as an additional site for reducing adsorbed Au3+ to the Au0 state. Conversely, most of the gold ions adsorbed by the PGBD matrix are retained in the Au3+ state due to the absence of supplementary sites, which ultimately causes to lower the adsorption efficiency, compared to PGBpy. It is important to note that the Namine and NBpy atoms present in these COFs can undergo protonation under acidic conditions and thus act as a reduction center for reducing Au3+ species to Au1+ and Au0 species.
This is also validated by the fact that in nitrogen-rich PGBpy, we have observed a greater extent of Au3+ reduction into Au0. To further support our proposed hypothesis and understand the mechanism of gold adsorption, the localization of electronic density in the COF skeleton was investigated using the DFT-based electron localization function (ELF) maps (Fig. 3e). The ELF map shows that Oether, Namine, and NBpy atoms are potential electron-donating sites for proton binding. Considering the presence of abundant O and N atoms within the COFs, we have estimated the binding energies for [AuCl4]− adsorption on O and N atoms (before and after protonation) in the representative PGBpy (Fig. 3f and h). The binding energies on protonated Namine, NBpy, and Oether were found to be −5.31 (Fig. 3g), −3.93 (Fig. 3i), and −4.23 eV (Fig. S37a, Section S13, ESI†). On the other hand, the direct adsorption of [AuCl4]− on the neutral COF structure showed unfavorable binding energy values of −1.44 for Namine (Fig. 3f), −1.39 for NBpy (Fig. 3g), and −1.63 eV for Oether (Fig. S37b, Section S13, ESI†). These findings suggest that the adsorption process begins with the protonation of active sites, which are subsequently attacked by H+ ions from [AuCl4]−. These protonated active sites act as reducing centers, thus facilitating the formation of an intermediate [AuCl3]− species by removing Cl− from [AuCl4]− accompanied by the generation of HCl molecules. Subsequently, the [AuCl3]− binds to another protonated active site and is further reduced to [AuCl2]−. This reduction process sequentially continues step-by-step till Au3+ ions are reduced to Au1+ and Au0 during adsorption. Importantly, the key step for Au3+ adsorption is the protonation of the active sites in our prepared COFs. This is further supported by the pH-dependent studies where the uptake efficiency of our prepared COFs is >99% while only <10% uptake efficiency in alkaline solution was observed (Fig. 2b). We want to highlight here that under acidic conditions, the protonation of the Oether atom can induce the partial cleavage of the ring, resulting in the formation of a hydroxyl group that can interact with [AuCl4]− more effectively, as suggested by the shift in the BE value observed in the high-resolution O 1s XPS profile of Au@PGBpy (Fig. S37b, Section S13, ESI†). However, in our case, FTIR studies show the benzoxazine core remains intact, which is confirmed by the absence of –OH stretch for Au@PGBpy (Fig. S23, Section S11, ESI†).
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5mh00698h |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2025 |