Description and analysis of the F-MOFs-M
The SEM micrographs of the sample are shown in Fig. 1b. The original fiberglass membrane has a glazed surface, excellent randomly oriented fiber, which is the construction of the MOFs vector. To grow MOFs on the fiberglass membrane, the ZrCl4 was doped onto the fiberglass membrane as a precursor of Zr-O bonds. The main reason to employ the fiberglass membrane was that the surface of the fiberglass membrane possessed many hydroxide radicals (Li et al. 2020). The precursor could not only act as a bridge to assist UIO-66-NH2(-COOH) to grow on the surface of the fiberglass membrane but also provide excellent mechanical stability for the final composite membrane F-MOFs-M. The UIO-66-NH2(-COOH) growing on the membrane surface seems like a “chocolate bar” (Fig. 1c). Herein, for the first time, we produced a straightforward two-step hydrothermal synthesis to create a composite membrane loaded with UIO-66-NH2(-COOH). The interior structure of the F-MOFs-M possessed abundant porosities, which favorited the adsorption of heavy metals from the wastewater. Additionally, physical synthesis had also been applied to the preparation of composite materials (Hakimifar &Morsali 2019). Huo et al (Huo &Yan 2012) reported a fast magnetization of MOFs MIL-101 for rapid magnetic solid-phase extraction of polycyclic aromatic hydrocarbons (PAHs) environmental water samples by physical synthesis. The F-MOFs-M with dispersive skeletons and tightly hierarchical structures are superior to the composite membrane prepared by the physical synthesis (Fig. 1d), in terms of providing a large contact area between the water and adsorbent while also enhancing membrane flux.
The XRD patterns of the UIO-66-NH2 and the F-MOFs-M are shown in Fig. 2a. The diffraction peaks of UIO-66-NH2(-COOH) on the surface of the membrane were similar to those for UIO-66-NH2 (Browe et al. 2017), indicating that the modification of UIO-66-NH2 with -COOH groups did not impact the original crystal structure (Katz et al. 2013). The chemical composition of UIO-66-NH2 and the F-MOFs-M are characterized by FT-IR spectra presented in Fig. 2b. The stretching vibration of the C-H, C-C, benzene ring, O-H, and C-N bond corresponded to the peaks at 1244 cm− 1, 1386 cm− 1, 1502 cm− 1, 2996 cm− 1, and 3475 cm− 1, respectively. The peaks at 1572 cm− 1(the C-O bond) and 1442 cm− 1 (the C-C vibration mode ) (Zhao et al. 2019) came from -COOH and C = O bond in the infrared spectra of UIO-66-NH2(-COOH). Two notable peaks at 3475 cm− 1and 3462 cm− 1 were ascribed to the amino groups, which were due to the different stretching modes of the N − H bond (El-Mehalmey et al. 2018). Besides, the vibration of the C = N and C = O bonds caused two significant peaks in the infrared spectra of UIO-66-NH2(-COOH) at 1576 cm− 1 and 1683 cm− 1. These peaks revealed that the –NH2 on UIO-66-NH2(-COOH) had been successfully converted from C − N to N-H using 2-Aminoterephthalic Acid and 1, 3-Benzenedicarboxylic acid.
X-ray photoelectron spectrometry is used to validate the chemical compositions of the F-MOFs-M (Fig. 3). As represented in Fig. 3a, the deconvolution of the Zr3d3/2 and Zr3d5/2 peaks at around 185.26 eV and 182.83 eV are attributed to zirconium in the Zr4+ state, indicating that Zr ions existed as ZrO2 in the tetragonal phase (Zhao et al. 2018). The O1s XPS spectrum is shown in Fig. 3b. The peak at 529.73 eV was attributed to lattice oxygen due to Zr-O bonds. The peak at approximately 531.68 eV was attributed to oxygen components of the carboxylate groups O-C = O due to the unreactive ligand. As shown in Fig. 3c, the binding energy of the N1s2/1 peaks at around 401.52 eV and 399.16 eV were attributed to -NH2/NH3+ and N-H in the 2-aminoterephthalic acid, respectively (Zhao et al. 2019). In a word, the surface of F-MOFs-M had a large number of –NH2 and –COOH, which contributed to improving the hydrophilicity and coordination.
The smallest unit cell of UIO-66-NH2(-COOH) consisted of a cubic structure of cationic Zr6O4(OH)4 nodes, which was indicated by ESI-HR-MS. As shown in Fig. 4, the characteristic peaks at m/z = 675.12, m/z = 681.41, m/z = 687.47, and m/z = 699.72 were ascribed to 90Zr6O8H7−, 91Zr6O8H7− 92Zr6O8H7− and 94Zr6O8H7−, respectively. The isotopic abundance ratios of 90Zr, 91Zr, 92Zr, and 94Zr were consistent with the natural isotopic composition of Zr. These data combined with Zr isotopic fingerprints allow identification of the Zr-MOF. Many researchers studied MOFs with hexazirconium oxo hydroxo (Zr6O4(OH)4) cluster nodes (Audu et al. 2016, Deria et al. 2015, Katz et al. 2013, Mondloch et al. 2014). For example, UIO-66 consisted of a cubic framework of cationic Zr6O4(OH)4 nodes (formed in situ via hydrolysis of ZrCl4) and 1,4-benzene dicarboxylate linkers (BDC). Unfortunately, the best operating condition of the instrument was a mass range of 50 to 1000 m/z using the measurement for negative mode. The theoretical mass of UIO-66-NH2(-COOH) overstepped the range of m/z = 1000. Therefore, the molecular weight of the complete MOFs could not be obtained. In this context, the characteristic peaks at m/z = 675.12, m/z = 681.41, m/z = 687.47, and m/z = 699.72 were ascribed to the smallest unit cell of UIO-66-NH2(-COOH). It could be deduced from the results that the UIO-66-NH2(-COOH) had at least two missing-linker sites per node (Fig. 4).
All the aforementioned investigations demonstrated that the F-MOFs-M had super hydrophilicity, sheet structure, and alive-linker site making it promising for the removal of harmful ions in wastewater.
Adsorption of performances
Effect of initial pH
The pH value should be optimized to improve the adsorption efficiency and reduce interference from the matrix. As shown in Fig. 5, the removal efficiencies of CrO42− by F-MOFs-M were higher than 75% at the pH ranging from 1 to 9, The high removal rate occurred at pH > 2; The removal rates of AsO43− were over 80% within the pH range of 1–10, and the highest removal efficiency occurred at pH 9. Insignificant influences of pH on SbO3−, Ni2+, and Pb2+ adsorptions were observed in the pH range of 1 to 8, and the optimal pH for relatively higher removal efficiencies for these metals was about 7. Cationic metals could be precipitated and negative charges on the surface of MOFs might reject anion groups at the pH values higher than 8. All the results indicated that the F-MOFs-M had an excellent capacity for capturing all the target metal ions in a wide pH range.
Adsorption isotherms and kinetics
Adsorption isotherms: the adsorption isotherms of all the tested ions on the F-MOFs-M are shown Fig. 6. It showed that the adsorption of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ increased rapidly at lower equilibrium concentrations, then a plateau appeared with the further increase of equilibrium concentration. Two often-used empirical equations, the Langmuir and Freundlich isotherm models, were used to describe the characteristics of the adsorption isotherms, which could be expressed as nonlinear and linear forms (Equ. i, ii):
\(qe=qmax*k1\frac{Ce}{1+k1*Ce} \left(\text{i}\right)\) \(\frac{Ce}{qe}=\frac{1}{qm*k1}+\frac{Ce}{qm} \left(\text{i}\text{i}\right)\)
where qe (mg·g− 1) was the adsorption amount of the target ions at equilibrium; qmax (mg·g− 1) was the maximum adsorption capacity; \(qm\) was the mono-layer adsorption capacity (mg·g− 1); the \(k1\) was the Langmuir constant, and \(Ce\) (mg·L− 1) was the equilibrium concentration.
The Freundlich models also had two forms expressed as:
\(qe=k2*{Ce}^{(1/n)}\) (iii)
$$log\left(qe\right)=log\left(k2\right)+\frac{1}{n}logCe \left(vi\right)$$
Where qe (mg·g− 1) was the adsorption amount at equilibrium concentration \(Ce\) (mg·L− 1) n and \(k2\) were the simulating equation constants relevant to adsorption characteristics.
The simulating results are listed in Tables 2 and 3. All the adsorption isotherms are well described by the two equations, with Langmuir equation being better for all the tested ions except for AsO43− as indicated by the higher determination coefficients (R2). The adsorption of AsO43− by F-MOFs-M might prefer multilayer adsorption, in consistent with the previous reports (Zhang et al. 2019b).
The maximum adsorption capacity (qmax) of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ by F-MOFs-M obtained from Langmuir equation were 54.96 mg·g− 1,188.9 mg·g− 1, 324.0 mg·g− 1, 132.5 mg·g− 1, and 15.87 mg·g− 1, respectively, which were accordingly lower than those of the synthesized bifunctional MOFs [UIO-66-NH2(-COOH)] material due to the partly occupations of the active sites used to connect with the membrane (Su et al. 2020). However, the maximum adsorption capacities of most target ions by F-MOFs-M were still higher than many reported MOFs (Table 3). Besides, MOFs loaded on fiber glass (F-MOFs-M) helped to increase the stability of MOFs, and supplied great advantage for the later handling of the absorbent after absorption by simply picking it up.
Adsorption kinetics: the absorption kinetics of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ on F-MOFs-M are shown in Fig. 7. It was clear that the concentration of target elements were decreased by rising time. In the first 30 min, CrO42− and SbO3− contents were rapidly decreased, while AsO43−, Pb2+, and Ni2+ contents were tardily reduced and after that, the adsorption rate of those target ions was remarkably decreased by raising the reaction time. After 120min of reaction, the adsorbed amounts of target elements on F-MOFs-M reached equilibrium, suggesting that the F-MOFs-M had fast adsorption dynamics for the removal of multifarious ions from water. The F-MOFs-M hydrolysis time is about 7 days, much longer than the equilibrium time of target elements, which guaranteed the removal of target elements in wastewater.
The experimental results were fitted with the pseudo-second-order kinetic model using the equation (v) below:
$$\frac{t}{Qt}=\frac{t}{Qe}+\frac{1}{k3*{Qe}^{2 }} \left(v\right)$$
Where k3was characteristic constants of the pseudo-second-order adsorption (g·mg− 1 min− 1), Qt was the amount of analyst adsorbed at the time \(t\) (mg·g− 1), and Qe (mg·g− 1) was the adsorption capacity unit of adsorbent at equilibrium (mg·g− 1).
The high correlation coefficients (R2) of all elements were revealed near 0.995 between time (t) and t/Qt (Fig. 7). Meantime, the adsorption rate coefficient (k3) of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ were calculated to be 1.28×10− 3 g·mg− 1 min− 1, 9.02×10− 2 g·mg− 1 min− 1, 1.78×10− 2 g·mg− 1 min− 1, 8.01×10− 3 g·mg− 1 min− 1, and 3.12×10− 3 g·mg− 1 min− 1, respectively, which were higher than those adsorbed materials under similar conditions (Table S1). The results indicated that the F-MOFs-M possessed a faster adsorption capacity. Moreover, the experimental data was certified to the calculated Qe value, indicating the adsorption process of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ onto the composite membranes fitted pseudo-second-order kinetic model.
Simultaneous capture of cations and anions
To study the simultaneous capture of cations and anions, the removal performances of F-MOFs-M on single solutions of AsO43−, CrO42−, SbO3−, Pb2+, and Ni2+ and those mixed solutions were researched (Fig. 8). The adsorption efficiency for each ion from the mixed solution is shown in Fig. 8a. It shows that oxyanions could be captured more preferentially than cations. The adsorption efficiencies were between 75% and 85% for oxyanions, and 20%-50% for cations. When the single cations separately appeared, each removal efficiency of Pb2+ and Ni2+ are 82% and 66%, and Cr3+, K+, Ca2+, and Mg2+ were just only 12%, 11%, and 9%, respectively (Fig. 8b). Therefore, with the presence of the oxyanions inhibited the capture of cations, and the more remarkable inhibiting effect was shown for Pb2+.
When multiple ions are competitive for adsorption, the adsorption partition-coefficient (Kd) is a relative metric of the selectivity adsorption (Cosgrove et al. 2019, Xu et al. 2020), and it was calculated by the equation below:
Kd5=Qe/Ce
Where Qe is the adsorption capacity after equilibration (mg/g); Ce is the equilibrium concentration (mg/L); Kd5 is the adsorption partition coefficient at the initial concentration of 5 mg/L.
As shown in Table 4, by comparing the Kd5 values the results show a selective adsorption order of SbO3− > CrO42− > AsO43−> Pb2+ > Ni2+. Generally, the electronegativity directly determined the adsorption selectivity of ions, the metal ions with stronger electronegativity easily interacted with the adsorbent (Pan et al. 2022). However, the adsorption capacity of Pb2+ was much lower than that of SbO3− (Liu et al. 2019). To explore this, the results of the Zeta potential of the F-MOFs-M show that the F-MOFs-M carries a positive charge on the surface during the reaction (Fig. 9). A positive surface charge favored the preferential adsorption of anion groups. Therefore, the adsorption selectivity depended on the surface charge of the F-MOFs-M and the self-characteristics of the adsorbate.
As observed, the F-MOFs-M had a better affinity to oxyanions in the presence of cations. The UIO-66-NH2(-COOH), which has two missing-linker sites per node at least (Fig. 4), captured oxyanions significantly faster and in greater quantities than cations. According to the Hofmeister effect (Hua et al. 2018, Nezamzadeh-Ejhieh &Afshari 2012), and many of the studies reported (Xiong et al. 2020, Xu et al. 2020) that the chemisorption ability of some metals obeyed the sequence of Pb2+༞Ni2+, Pb2+༞Cd2+, Ni2+༞Cd2+. Herein, the F-MOFs-M obtained a preferential adsorption order of oxyanions(SbO3− > CrO42− > AsO43−)༞cations(Pb2+༞Ni2+).
Recycling and stability of the F-MOFs-M
The reproducibility of the resulting membrane is crucial for its practical application (Glomstad et al. 2017). And therefore, the elution and reusage of the F-MOFs-M were further studied. The SbO3−, CrO42−, AsO43−, Pb2+, and Ni2+ loaded membranes were stirred continuously in HNO3 (0.5 M) at 25℃ for 360min. the results showed that the cations-loaded (Pb2+, and Ni2+) F-MOFs-M were easily eluted, achieving almost 98.5% and 97.6% regeneration rates. Meanwhile, the anions (SbO3−, CrO42−, and AsO43−)-loaded F-MOFs-M were also easily eluted with over 98% elution rates. To obtain the recycling of the F-MOFs-M, the regenerated membranes were then used for the 5 times regenerated cycles. It is certain that the F-MOFs-M still possesses over 90% of the original removal rates even after 5 times reuse (Fig. 10), which indicated that the F-MOFs-M could be regenerated and recycled without losing its removal efficiency.
The bottleneck of MOFs in practical application as absorbents is their structural instability in solution. Therefore, the hydrolysis process of the fabricated F-MOFs-M in an aqueous solution was investigated using HR-ICP-MS (Table 1) within 3 h to 60 days. The results showed that the UIO-66-NH2(-COOH) could stably exist in an aqueous solution within at least 7 days. Obviously, it possessed an excellent chemical stability, which was closely related to the electrostatic configuration in water.
Adsorption Mechanisms
The previous investigations showed the F-MOFs-M exhibited superior adsorption performance for the tested harmful ions, which might attribute to the abundant functional groups and the exposed active sites in the UIO-66-NH2(-COOH) structure. As shown in Fig. 2a, after-adsorption treatment, the sharp diffraction modes corresponding to the extruded UIO-66-NH2(-COOH) arise three broad peaks at around 14.68°, 19.33°, and 51.70° following the (400), (420), and (137) diffraction peaks of the metal peaks, indicating the successful adsorption of heavy metals. The second evidence comes from Fig. 2b by the FT-IR spectra, the infrared spectrum of UIO-66-NH2(-COOH), the peak of the O-H disappeared and the metal-O peak appeared. A characteristic peak at about 875 cm− 1 proved the existence of a metal-O bond after adsorption, and three peaks at 564 cm− 1, 584 cm− 1and 568 cm− 1 were closely related to the symmetric tensile vibration of Zr-O (Audu et al. 2016, Huo et al. 2019). For instance, the XPS spectra of As3d (Fig. 11) is shown two visible shifts at 46.72 eV to 45.68 eV after AsO43− loading, suggesting that the chemical environment of AsO43− had changed. Furthermore, the peaks resulting from the Zr isotope ratio were used to calculate the smallest unit cell of UIO-66-NH2(-COOH) by ESI-HR-MS. The results show that the Zr atoms contained unsaturated binding sites, and target oxyanions are coordinated to the Zirconium oxide octahedron nodes (Fig. 4) (Katz et al. 2013). The active binding sites on the Zirconium oxide octahedron nodes were excellent binders for the AsO43−.
Additionally, as displayed in N1s and O1s spectrum after-absorption (Fig. 3c and b), the characteristic peaks at 401.52 eV and 399.16 eV are ascribed to -N-H, N+, and C-N bonds, respectively. The -NH and N+ own a higher peak after absorption, suggesting that the cations are being promoted to obtain inert nitrogen. Compared with the O1s XPS spectrum of before-absorption (Fig. 3b), the peak for chemisorbed oxygen is more prominent. The peak at 533.72 eV could consider being –COO−, decreased its binding energy, which might be formed the hydrogen bonds with cations or promoted the transformation of C = O oxygen species as active oxygen. According to reports (Liu et al. 2020, Tang et al. 2021, Zhao et al. 2019), the C = O bond could increase the chemical adsorption capacities of ions such as Hg2+, Pb2+, and Ni2+. As mentioned above, the lone pair of electrons functional group went through coordination and covalence with target cations.