Synthesis of Chemically Modified Acid-Functionalized Multiwall Carbon Nanotubes with Benzimidazole for Removal of Lead and Cadmium Ions from Wastewater

In this work, acid-functionalized multiwalled carbon (MWCNTs–CO2H) nanotube was successfully functionalized with a heterocyclic scaffold, namely benzimidazole, to give novel functionalized multiwalled carbon nanotubes (BI@MWCNTs). Then, FTIR, XRD, TEM, EDX, Raman spectroscopy, DLS, and BET analyses were used to characterize the synthesized BI@MWCNTs. The effectiveness of the adsorption of two heavy metal ions, Cd2+ and Pb2+, in single metal and mixed metal solutions on the prepared material was investigated. Influencing parameters for the adsorption method, for example duration, pH, starting metal concentration, and BI@MWCNT dosage, were examined for both metal ions. Moreover, adsorption equilibrium isotherms fit with the Langmuir and Freundlich models perfectly, while the intra-particle diffusion models provide pseudo-second order adsorption kinetics. The adsorption of Cd2+ and Pb2+ ions onto BI@MWCNTs revealed an endothermic and a spontaneous method with great affinity as a result of the negative values of Gibbs free energy (ΔG) and the positive values of enthalpy (ΔH) and entropy (ΔS). Both Pb2+ and Cd2+ ions were completely eliminated from aqueous solution (100 and 98%, respectively) using the prepared material. Additionally, BI@MWCNTs have a high adsorption capacity and were regenerated in a simple way and reused for six cycles, which make them a cost-effective and efficient absorbent for the removal of such heavy metal ions from wastewater.


Introduction
According to Food and Agriculture Organization (FAO) of the United Nations report in 2020, "3.2 billion people live in agricultural areas with high to very high water shortages or scarcity, of whom 1.2 billion people-roughly one-sixth of the world's population live in severely water-constrained agricultural areas". The scarcity of water, mainly in arid and semi-arid areas of the world, is exerting exceptional pressure on sources and necessitates offering satisfactory water for human and different uses. Water recycle/reuse has been confirmed to be successful and promising for reliable water delivery. In many developing countries, this has led to the reuse of community wastewater in agriculture for irrigation. Nevertheless, it will increase human exposure to the pollutants. Inorganic pollutants form a major part of this problem, such as heavy metal ions. There are abundant metal ions which are significantly toxic to human beings and ecological environments, including lead (Pb), cadmium (Cd), mercury (Hg), zinc (Zn), nickel (Ni), etc. They exist, in various forms, in the soil and in natural and wastewater, and may become contaminants in food and drinking water. It was reported that plants grown and irrigated with untreated wastewater were highly polluted by heavy metals, which led to many health problems [1]. This was explained as a first result of heavy metal accumulation from irrigation water in soil [2]. Therefore, eliminating or decreasing the heavy metal contaminants in wastewater to safe limits is an important way to avoid their accumulation into the soil and, consequently, in was suspended in dry toluene (50 mL), and 2-aminobenzimidazole (0.1 g, 0.7 mmol) was added. The reaction mixture was subjected to sonication for 1 h, then reflux for 24 h. The hot reaction mixture was filtered off and washed with methanol. The resulting black solid was dried to yield BI@MWCNTs in an excellent yield (>95%, 0.1 g). The reaction steps are shown in Scheme 1. benzimidazole (BI@MWCNTs). The synthesis steps are explained hereafter. Firstly, the MWCNT-COOH (0.1 g) was sonicated with an excess of thienyl chloride for 2 h, then refluxed in a water bath for 24 h. The excess thienyl chloride was removed under vacuum (caution) to yield the acid chloride MWCNTs-COCl (0.1 g). The obtained acid chloride was suspended in dry toluene (50 mL), and 2-aminobenzimidazole (0.1 g, 0.7 mmole) was added. The reaction mixture was subjected to sonication for 1 h, then reflux for 24 h. The hot reaction mixture was filtered off and washed with methanol. The resulting black solid was dried to yield BI@MWCNTs in an excellent yield (>95%, 0.1 g). The reaction steps are shown in Scheme 1.

Characterization of Functionalized MWCNT Materials
Data from an X-ray diffractometer (Shimadzu-7000, Torrance, CA, USA) using a CuK radiation beam (=0.154060 nm) was collected. Fourier transform infrared (FTIR) analysis with a range of 400 to 4000 cm −1 was used (A Bruker ALFA FTIR spectrometer, Billerica, MA, USA). Additionally, a scanning electron microscope (a JOEL scanning electron microscope-type electron probe (SEM) JXA-840, Tokyo, Japan) was used to investigate the surface morphology of the synthesized nanomaterials. An XploRA PLUS device (Horiba, Kyoto, Japan) verified the Raman spectra. Furthermore, particle size distribution of prepared samples was characterized using dynamic light scattering (DLS) that was carried out by using Otsuka Electronics DLS-6500 system with the He-Ne laser (λ = 633 nm, 10 mW). The temperature of the sample solution was controlled at ∼25 °C. In addition, the surface area and total pore volume of Brunauer-Emmett-Teller (BET) structures were determined.

Adsorption Tests
Applying a thermostatic shaker set at 150 rpm, a specified weight of BI@MWCNTs was added to 50 mL of worked solutions of metal salt at 25 °C in order to conduct the adsorption assays. By diluting 1 g/L standard solution of Pb(NO3)2 or Cd(NO3)2.4H2O with deionized water, altered starting metal solution concentrations (50, 100, 150, and 200 mg/L) were created. Subsequently, for the adsorption operations, the liquid and solid phases of the BI@MWCNTs were separated using centrifugation at 4000 rpm for 15 min using atomic adsorption spectroscopy; the amount of remaining metal in the solution was calculated. The quantities of adsorbed metal ions were calculated using the following equations:

Characterization of Functionalized MWCNT Materials
Data from an X-ray diffractometer (Shimadzu-7000, Torrance, CA, USA) using a CuK radiation beam (=0.154060 nm) was collected. Fourier transform infrared (FTIR) analysis with a range of 400 to 4000 cm −1 was used (A Bruker ALFA FTIR spectrometer, Billerica, MA, USA). Additionally, a scanning electron microscope (a JOEL scanning electron microscope-type electron probe (SEM) JXA-840, Tokyo, Japan) was used to investigate the surface morphology of the synthesized nanomaterials. An XploRA PLUS device (Horiba, Kyoto, Japan) verified the Raman spectra. Furthermore, particle size distribution of prepared samples was characterized using dynamic light scattering (DLS) that was carried out by using Otsuka Electronics DLS-6500 system with the He-Ne laser (λ = 633 nm, 10 mW). The temperature of the sample solution was controlled at~25 • C. In addition, the surface area and total pore volume of Brunauer-Emmett-Teller (BET) structures were determined.

Adsorption Tests
Applying a thermostatic shaker set at 150 rpm, a specified weight of BI@MWCNTs was added to 50 mL of worked solutions of metal salt at 25 • C in order to conduct the adsorption assays. By diluting 1 g/L standard solution of Pb(NO 3 ) 2 or Cd(NO 3 ) 2 .4H 2 O with deionized water, altered starting metal solution concentrations (50, 100, 150, and 200 mg/L) were created. Subsequently, for the adsorption operations, the liquid and solid phases of the BI@MWCNTs were separated using centrifugation at 4000 rpm for 15 min using atomic adsorption spectroscopy; the amount of remaining metal in the solution was calculated. The quantities of adsorbed metal ions were calculated using the following equations: where the quantities of adsorbed metal per unit weight of the adsorbent at time (t) and at equilibrium are represented by the variables (q t ) and (q e ) (mg g −1 ), respectively. The metal concentrations are (C 0 , C t , and C e ) (mg L −1 ) at the beginning time, time t, and at equilibrium time, respectively. V (L) denotes the volume of the metal solution, and m (g) is the adsorbent's weight. The percent of removal (R%) was determined as follows: The following equation was used to evaluate the distribution coefficient (K d ) of Cd 2+ and Pb 2+ onto the produced adsorbent: When evaluating the actual performance of the adsorbent, the distribution coefficient or partition coefficient (K d ) may be utilized.
The following equation was used to investigate the selectivity coefficient (α) for the obligatory nature of a certain adsorbate in the presence of intrusive ions: where K d (I) is the K d value of the other metal in the multi-metal solutions, which in this case is cadmium, and K d (T) is the K d value of the targeted metal.

Adsorption Kinetics Process
Since they represent the rate at which an adsorbent is absorbed, kinetic studies are crucial. Based on kinetic studies, the rate of adsorption method and the mechanism of the adsorption can be clarified. Based on kinetic measurements with starting Cd 2+ and Pb 2+ ion concentrations ranging from 50 to 200 mg/L, the pH was set at 5.5 (to prevent precipitate of hydroxyl) and the adsorption duration was varied from 1 to 30 min with a dosage of 0.02 mg/mL of adsorbent for our kinetic investigations. The pseudo-first order and pseudo-second order models were used to clarify the kinetics of adsorption [26,27]. Equations (6) and (7) are as follows: q t = k 2 q e 2t 1 + k 2 q e t × 100 (7) where k 1 and k 2 are the adsorption rate constants of the pseudo-first-and pseudo-second order, respectively. For low solute concentrations, the pseudo-first order kinetic model is more appropriate. One way to express it in writing is explained as follows: where k 1 (min −1 ) is the rate constant of the pseudo-first order adsorption, q e is the amount of metal adsorbed at saturation per gram of adsorbent, and q t is the amount of metal adsorbed at time t per gram of adsorbent (mg gm −1 ), respectively. The pseudo-second order kinetic model depends on both the quantity of solute adsorbed at equilibrium and the amount adsorbed on the surface of the adsorbent. An intra-particle diffusion model was applied at various starting concentrations to examine the diffusion mechanism of Cd(II) and Pb(II) onto produced BI@MWCNT adsorbents [27], as follows: where C is the magnitude of the intercept, a constant that indicates the importance of the boundary layer or mass transfer effect, and k i (mg g −1 min −1/2 ) is the intra-particle diffusion rate constant, which is the slope of the straight line of q t vs. t 1/2 .

Adsorption Isotherms
Using the Langmuir and Freundlich isotherms, the adsorption thermodynamics and isotherms were put to the test to confirm the metal uptake behavior of the BI@MWCNTs. The adsorption isotherms was used to show the sorption capacity of the BI@MWCNTs at various starting concentrations at equilibrium. Adsorption isotherms provide a full understanding of the nature of interactions by describing how the adsorbate interacts with adsorbents. For this investigation, a number of isotherm equations have been generated and used, and the two significant isotherms were utilized. By utilizing several adsorption isotherm models, the mechanism of interaction between the adsorbed metal and BI@MWCNTs was investigated. The Langmuir and Freundlich models suit the experimental data well [28,29]. The Langmuir model is stated as follows and takes into account the progression of an adsorbate monolayer on an adsorbent's normal surface: where C e is the equilibrium adsorbate concentration (in mmol L −1 ), q e is the quantity of adsorbate adsorbed per unit weight of adsorbent (in mmol g −1 ), q m is the adsorption capacity (mmol g −1 ), or monolayer capacity, and k L is a constant (L mmol −1 ), as shown in the following equations [9,28]. A key feature of the Langmuir model is the separation factor (R L ), a dimensionless quantity. The following is how the R L equation is written out: where k L is the Langmuir constant (L mg −1 ) and C o is the adsorbate's maximum initial concentration (mg L −1 ). The value of R L determines whether the isotherm's form is irreversible (R L = 0), favorable (0 < R L < 1), linear (R L = 1), or unfavorable (R L > 1). The Freundlich isotherm, on the other hand, takes into account adsorption between adsorbates and adsorbents with a heterogeneous surface. The energy level at the adsorption sites reveals that the rate of adsorption varies. In order to obtain Freundlich isotherms, we use the following equation: ln q e = ln k f + 1 (n ln C e ) (12) where k f is a constant (function of temperature and adsorption energy) and n is a constant related to adsorption intensity, determined by plotting ln q e vs. ln C e , which gives a straight line with slope of 1/n and an intercept of ln k f . The magnitude of the "n" shows an indication of the favorability of adsorption.

Adsorption Thermodynamics
At 25 • C and 35 • C, thermodynamic adsorption studies were conducted with an initial concentration of 100 mg L −1 . The following equations [27,30] were used to compute the thermodynamic parameters: where K D is the equilibrium partition constant, ∆G is the Gibbs free energy change (kJ mol −1 ), R is the universal gas constant (8.314 J mol −1 K −1 ), T (K) is the temperature, ∆H is the enthalpy change (kJ mol −1 ), and ∆S is the entropy change (kJ mol −1 K −1 ). The K D values for each temperature were used to calculate the values of ∆G, while the slope and intercept of the ln K D vs. 1/T plot were used to determine the values of ∆H and ∆S, respectively.

Test for Regeneration and Reuse
Evaluation of an adsorbent's practical viability requires regeneration and reuse tests. In order to test the reversibility of metal ion adsorption with BI@MWCNTs, 4 mg of BI@MWCNTs were combined with 100 mL of 50 mg L −1 Cd(II) and Pb(II). The mixtures were kept for 15 min with shaking at 120 rpm. The adsorbent was separated from the mixture by centrifugation. Separated adsorbents were washed, dried in an oven at 40 • C for two hours, and then reused. A total of 0.1 M HCl was added after the adsorption process, and the solution was then ultrasonically treated at ambient temperature for 30 min. Metal ion concentration was determined using atomic adsorption spectroscopy. To investigate regenerability, three cyclic adsorption-desorption processes were used.

Characterization of BI@MWCNTs
To investigate the successful functionalization of the MWCNTs-CO 2 H with BI, the infrared spectra of BI, MWCNTs-CO 2 H, and BI@MWCNTs were verified. Figure 1 shows the FTIR spectra of BI, MWCNTs-CO 2 H, and BI@MWCNTs. The spectra indicate the disappearance of the primary amino pair peak of the aminobenzimidazole in the prepared BI@MWCNTs with the presence of the -OH absorption band of the remaining unfunctionalized carboxylic group at 3333 cm −1 . A clear peak, which could be attributed to conjugated C=O stretching vibrations, was seen at 1630 cm −1 . The presence of additional peaks in the range of 1600-1500 cm −1 were attributed to the aromatic ring system and C=N of the benzimidazole moiety. Furthermore, the 3020 cm −1 , 1215 cm −1 , and 769 cm −1 peaks in the BI@MWCNTs' spectrum, which are brought on by the aromatic C-H stretching, C-N stretching, and C-H out of plane bending modes, respectively, confirmed the functionalization of MWNCTs-CO 2 H with BI [31,32].

Test for Regeneration and Reuse
Evaluation of an adsorbent's practical viability requires regeneration and reuse te In order to test the reversibility of metal ion adsorption with BI@MWCNTs, 4 mg BI@MWCNTs were combined with 100 mL of 50 mg L −1 Cd (II) and Pb (II). The mixtu were kept for 15 min with shaking at 120 rpm. The adsorbent was separated from mixture by centrifugation. Separated adsorbents were washed, dried in an oven at 40 for two hours, and then reused. A total of 0.1 M HCl was added after the adsorption p cess, and the solution was then ultrasonically treated at ambient temperature for 30 m Metal ion concentration was determined using atomic adsorption spectroscopy. To inv tigate regenerability, three cyclic adsorption-desorption processes were used.

FTIR
To investigate the successful functionalization of the MWCNTs-CO2H with BI, infrared spectra of BI, MWCNTs-CO2H, and BI@MWCNTs were verified. Figure 1 sho the FTIR spectra of BI, MWCNTs-CO2H, and BI@MWCNTs. The spectra indicate the d appearance of the primary amino pair peak of the aminobenzimidazole in the prepa BI@MWCNTs with the presence of the -OH absorption band of the remaining unfuncti alized carboxylic group at 3333 cm −1 . A clear peak, which could be attributed to conjuga C=O stretching vibrations, was seen at 1630 cm −1 . The presence of additional peaks in range of 1600-1500 cm −1 were attributed to the aromatic ring system and C=N of the b zimidazole moiety. Furthermore, the 3020 cm −1 , 1215 cm −1 , and 769 cm −1 peaks in BI@MWCNTs' spectrum, which are brought on by the aromatic C-H stretching, C stretching, and C-H out of plane bending modes, respectively, confirmed the function zation of MWNCTs-CO2H with BI [31,32].   Figure 2A. The newly prepared material, BI@MWCNTs, is shown in Figure 2B. The BI adsorption causes the nanotubes' diameter to grow up to 20 nm and gives them their distinctive shading.  Figure 2A. The newly prepared material, BI@MWCNTs, is shown in Figure 2B. The BI adsorption causes the nanotubes' diameter to grow up to 20 nm and gives them their distinctive shading.

Raman Spectroscopy
Raman spectra provided additional support for MWCNTs functionalization. As shown in Figure 3, both the spectra for MWCNTs-COOH and BI@MWCNTs contain the distinctive G (tangential mode) and D (disorder mode) bands linked to the sp 2 carbon and the carbon nanotube at 1585 cm −1 and 1330 cm −1 , respectively. Furthermore, the intensity of the D band for the MWCNTs-CO2H is more than the intensity of the D band for the BI@MWCNTs, which confirmed that the functionalization with BI had successfully occurred. In addition, the intensity ratio ID/IG for BI@MWCNTs (0.98) was higher than that for MWCNTs-CO2H (0.85), which provided additional proof that the out of the plane rotation of the heterocyclic side chain contributed more to these materials than MWCNTs-CO2H [32].

Raman Spectroscopy
Raman spectra provided additional support for MWCNTs functionalization. As shown in Figure 3, both the spectra for MWCNTs-COOH and BI@MWCNTs contain the distinctive G (tangential mode) and D (disorder mode) bands linked to the sp 2 carbon and the carbon nanotube at 1585 cm −1 and 1330 cm −1 , respectively. Furthermore, the intensity of the D band for the MWCNTs-CO 2 H is more than the intensity of the D band for the BI@MWCNTs, which confirmed that the functionalization with BI had successfully occurred. In addition, the intensity ratio ID/IG for BI@MWCNTs (0.98) was higher than that for MWCNTs-CO 2 H (0.85), which provided additional proof that the out of the plane rotation of the heterocyclic side chain contributed more to these materials than MWCNTs-CO 2 H [32].  Figure 2A. The newly prepared material, BI@MWCNTs, is shown in Figure 2B. The BI adsorption causes the nanotubes' diameter to grow up to 20 nm and gives them their distinctive shading.

Raman Spectroscopy
Raman spectra provided additional support for MWCNTs functionalization. As shown in Figure 3, both the spectra for MWCNTs-COOH and BI@MWCNTs contain the distinctive G (tangential mode) and D (disorder mode) bands linked to the sp 2 carbon and the carbon nanotube at 1585 cm −1 and 1330 cm −1 , respectively. Furthermore, the intensity of the D band for the MWCNTs-CO2H is more than the intensity of the D band for the BI@MWCNTs, which confirmed that the functionalization with BI had successfully occurred. In addition, the intensity ratio ID/IG for BI@MWCNTs (0.98) was higher than that for MWCNTs-CO2H (0.85), which provided additional proof that the out of the plane rotation of the heterocyclic side chain contributed more to these materials than MWCNTs-CO2H [32].

XRD
The recorded MWCNTs-CO 2 H and BI@MWCNTs X-ray diffraction (XRD) patterns are shown in Figure 4 as intensity counts vs. 2θ. One strong diffraction peak at 2θ = 25.4 • , which corresponds to the (002) plane, and one low-intensity peak at 2θ = 43 • , which corresponds to the (100) reflection plane, were identified in the MWCNTs-CO 2 H and BI@MWCNTs' diffractograms (JCPDS card no. 75-1621) and were attributed to the graphitic structure. When MWCNTs-CO 2 H was converted to BI@MWCNTs, no appreciable changes in the 2θ values of the XRD patterns were seen. The enhanced crystallinity of BI@MWCNTs is evidenced by the fact that the diffraction peak intensity of BI@MWCNTs is higher than that of MWCNTs-CO 2 H. This might be because when MWCNTs-CO 2 H underwent treatment, reaction process led to more purification and the removal of amorphous carbon and other impurities in the carbon nanotube.

XRD
The recorded MWCNTs-CO2H and BI@MWCNTs X-ray diffraction (XRD) patterns are shown in Figure 4 as intensity counts vs. 2θ. One strong diffraction peak at 2θ = 25.4°, which corresponds to the (002) plane, and one low-intensity peak at 2θ = 43°, which corresponds to the (100) reflection plane, were identified in the MWCNTs-CO2H and BI@MWCNTs' diffractograms (JCPDS card no. 75-1621) and were attributed to the graphitic structure. When MWCNTs-CO2H was converted to BI@MWCNTs, no appreciable changes in the 2θ values of the XRD patterns were seen. The enhanced crystallinity of BI@MWCNTs is evidenced by the fact that the diffraction peak intensity of BI@MWCNTs is higher than that of MWCNTs-CO2H. This might be because when MWCNTs-CO2H underwent treatment, reaction process led to more purification and the removal of amorphous carbon and other impurities in the carbon nanotube.

Dynamic Light Scattering (DLS)
Dynamic light scattering (DLS) analyses of 2-aminobenzimidazole (BI), MWCNTs-CO2H, and the functionalized carbon nanotubes with benzimidazole (BI@MWCNTs) were studied to affirm the particle size distribution. The samples were analyzed throughout the experimental period, from 0 to 3 h, to study the dispersion stability of the material in glyceline and water, as a reference. All samples were dispersed in glyceline and water to suitable concentrations and were analyzed in triplicate. The obtained size of the prepared particles is bigger than what was previously observed through SEM and ranged from 600 to 2000 nm, because the DLS investigation measures the hydrodynamic diameter and the characteristics of nanoparticles exclusively depending on their size. There are numerous techniques available for describing the surfaces of nanoparticles. Due to the differences in how to measure the particle dimensions, the pore size distribution values of functionalized MWCNTs are different from those of SEM ( Figure 2) compared to DSL, as shown in Figure 5. This is due to the fact that particles are suspended in a solution for particle size distribution (DLS), and measurements are taken while the particles are in solution using

Dynamic Light Scattering (DLS)
Dynamic light scattering (DLS) analyses of 2-aminobenzimidazole (BI), MWCNTs-CO 2 H, and the functionalized carbon nanotubes with benzimidazole (BI@MWCNTs) were studied to affirm the particle size distribution. The samples were analyzed throughout the experimental period, from 0 to 3 h, to study the dispersion stability of the material in glyceline and water, as a reference. All samples were dispersed in glyceline and water to suitable concentrations and were analyzed in triplicate. The obtained size of the prepared particles is bigger than what was previously observed through SEM and ranged from 600 to 2000 nm, because the DLS investigation measures the hydrodynamic diameter and the characteristics of nanoparticles exclusively depending on their size. There are numerous techniques available for describing the surfaces of nanoparticles. Due to the differences in how to measure the particle dimensions, the pore size distribution values of functionalized MWCNTs are different from those of SEM ( Figure 2) compared to DSL, as shown in Figure 5. This is due to the fact that particles are suspended in a solution for particle size distribution (DLS), and measurements are taken while the particles are in solution using Brownian motion [33]. The hydrodynamic diameter may increase as a result of the individual particles aggregating. Brownian motion [33]. The hydrodynamic diameter may increase as a result of the individual particles aggregating.

BET Analysis
Brunauer-Emmett-Teller (BET) analysis was used to measure the specific surface areas of the adsorbents, and an automatic surface analyzer was used to measure the total pore volume and pore size (Quantachrome, Boynton Beach, FL, USA). Figure 6 shows N2 adsorption isotherms, which measured at 77 K, and nitrogen adsorption-desorption isotherms for the MWCNTs-CO2H and BI@MWCNTs. The result is summarized in Table 1 with the BET measurements of the MWCNTs-CO2H and BI@MWCNTs. Due to the intertubular structure, hysteresis was seen for both MWCNTs-CO2H and BI@MWCNTs. The noticed decrease effect may be attributed to the well-known capillary condensation [16]. The isotherms obtained from the various samples were categorized as type IV isotherms with H3 type hysteresis loops by the original IUPAC classification [16]. This finding suggests that samples may contain mesoporous and macroporous structures. The microporous structure with a mean pore size of 14 nm in the sample is responsible for the abrupt increase in the N2 adsorption plot of MWCNTs-CO2H at a relative pressure of > 0.9. Furthermore, the MWCNTs-CO2H have high surface area (252.6 m 2 g −1 ) and high pore volume (0.799 cm 3 g −1 ), implying the presence of additional sites, functionalization with BI, and co-lateral growth of BI. In contrast, the BI@MWCNT surface area and pore volume were decreased to 45.3 m 2 g −1 and 0.3041 cm 3 g −1 , respectively.

BET Analysis
Brunauer-Emmett-Teller (BET) analysis was used to measure the specific surface areas of the adsorbents, and an automatic surface analyzer was used to measure the total pore volume and pore size (Quantachrome, Boynton Beach, FL, USA). Figure 6 shows N 2 adsorption isotherms, which measured at 77 K, and nitrogen adsorption-desorption isotherms for the MWCNTs-CO 2 H and BI@MWCNTs. The result is summarized in Table 1 with the BET measurements of the MWCNTs-CO 2 H and BI@MWCNTs. Due to the intertubular structure, hysteresis was seen for both MWCNTs-CO 2 H and BI@MWCNTs. The noticed decrease effect may be attributed to the well-known capillary condensation [16]. The isotherms obtained from the various samples were categorized as type IV isotherms with H3 type hysteresis loops by the original IUPAC classification [16]. This finding suggests that samples may contain mesoporous and macroporous structures. The microporous structure with a mean pore size of 14 nm in the sample is responsible for the abrupt increase in the N 2 adsorption plot of MWCNTs-CO 2 H at a relative pressure of > 0.9. Furthermore, the MWCNTs-CO 2 H have high surface area (252.6 m 2 g −1 ) and high pore volume (0.799 cm 3 g −1 ), implying the presence of additional sites, functionalization with BI, and co-lateral growth of BI. In contrast, the BI@MWCNT surface area and pore volume were decreased to 45.3 m 2 g −1 and 0.3041 cm 3 g −1 , respectively.

Effect of pH
One of the key influences regulating the method of the adsorption is the solution pH, which has an impact on how the adsorbate interacts with the adsorbent surface charges. In order to find the ideal pH level, the influence of the pH on the adsorption effectiveness of BI@MWCNTs toward Pb 2+ and Cd 2+ was investigated in the pH range of 2.5 to 10. The percentage of removal significantly enhanced after the pH was raised from 3 to 6, as shown in Figure 7. Nevertheless, the removal effectiveness progressed slowly when the solution was in the range of 5-6. The quantity of adsorption reduced after pH surpassed 7. The Pb 2+ and Cd 2+ adsorption behavior on the produced BI@MWCNTs in relation to pH values might be explained by a variety of factors. For instance, at lower pH values (less than 4), the surfaces of the BI@MWCNTs would be surrounded by a large number of H3O + ions that contest with the metal ions. Furthermore, at lower pH values (less than 4), the nitrogen atoms of the benzimidazole scaffold (either the pyridine-or pyrrole-like) might be easily protonated, which would subsequently affect the adsorption efficiency of the BI@MWCNTs. However, increasing the pH of the solution to higher values (the weak acidic condition) to a range of 5-6 increases the electrostatic attraction between the metal ions and the surface of the BI@MWCNTs, progressively improving the removal effectiveness [34]. Furthermore, Pb 2+ and Cd 2+ are present in varied formulae in solution at altered pH levels, while when the pH exceeds 7, the precipitates of hydroxides (Cd(OH)2., Cd(OH)3-Pb(OH)2, and Pb(OH)3) are developed rapidly, and these precipitates are hardly adsorbed onto the BI@MWCNTs [34]. This is because the benzimidazole nucleus has high-

Effect of pH
One of the key influences regulating the method of the adsorption is the solution pH, which has an impact on how the adsorbate interacts with the adsorbent surface charges. In order to find the ideal pH level, the influence of the pH on the adsorption effectiveness of BI@MWCNTs toward Pb 2+ and Cd 2+ was investigated in the pH range of 2.5 to 10. The percentage of removal significantly enhanced after the pH was raised from 3 to 6, as shown in Figure 7. Nevertheless, the removal effectiveness progressed slowly when the solution was in the range of 5-6. The quantity of adsorption reduced after pH surpassed 7. The Pb 2+ and Cd 2+ adsorption behavior on the produced BI@MWCNTs in relation to pH values might be explained by a variety of factors. For instance, at lower pH values (less than 4), the surfaces of the BI@MWCNTs would be surrounded by a large number of H 3 O + ions that contest with the metal ions. Furthermore, at lower pH values (less than 4), the nitrogen atoms of the benzimidazole scaffold (either the pyridine-or pyrrole-like) might be easily protonated, which would subsequently affect the adsorption efficiency of the BI@MWCNTs. However, increasing the pH of the solution to higher values (the weak acidic condition) to a range of 5-6 increases the electrostatic attraction between the metal ions and the surface of the BI@MWCNTs, progressively improving the removal effectiveness [34]. Furthermore, Pb 2+ and Cd 2+ are present in varied formulae in solution at altered pH levels, while when the pH exceeds 7, the precipitates of hydroxides (Cd(OH) 2 ., Cd(OH) 3− Pb(OH) 2 , and Pb(OH) 3 ) are developed rapidly, and these precipitates are hardly adsorbed onto the BI@MWCNTs [34]. This is because the benzimidazole nucleus has high-class structural features and an electron-rich scaffold, which might produce an electrostatic repulsion. At a pH greater than 8, the adsorption of Pb 2+ and Cd 2+ fall significantly and it becomes problematic to eliminated using BI@MWCNTs. Therefore, in the next studies, a pH of 5 will be the pH of choice for eliminating Pb 2+ and Cd 2+ and will be applied in this study. class structural features and an electron-rich scaffold, which might produce an electrostatic repulsion. At a pH greater than 8, the adsorption of Pb 2+ and Cd 2+ fall significantly and it becomes problematic to eliminated using BI@MWCNTs. Therefore, in the next studies, a pH of 5 will be the pH of choice for eliminating Pb 2+ and Cd 2+ and will be applied in this study.

Adsorption Competition
Affording to adsorbent's adsorption capacity ratio to the final concentration of adsorbate in the liquid at equilibrium, the distribution coefficient (Kd) essentially signifies the accurate presentation of an adsorbent [31,35], whereas the adsorption capacity and removal effectiveness can be changed by operational settings. In this investigation, Equation (5) was utilized to assess the effectiveness of the produced adsorbent material using the distribution coefficient (Kd) for Pb 2+ and Cd 2+ . The values of Kd for Pb 2+ and Cd 2+ increased with increasing pH to reach pH 6 and then decreased by further pH increasing, as shown in Figure 8. This is because higher pH activates more adsorption locations and makes extra metal ions accessible to transfer from the solution to the BI@MWCNTs' surface. Additionally, at each pH level, BI@MWCNTs demonstrated a more definite favorite for Pb 2+ over Cd 2+ . These can be attributed to Pb 2+ showed strong affinity for BI@MWCNTs in the solution containing two metals. Table 2 shows the distribution coefficient values of Pb 2+ are more than that of Cd 2+ , and this is showing that BI@MWCNTs exhibited selective adsorption for Pb 2+ over Cd 2+ according to distribution coefficient values onto BI@MWCNTs.

Adsorption Competition
Affording to adsorbent's adsorption capacity ratio to the final concentration of adsorbate in the liquid at equilibrium, the distribution coefficient (K d ) essentially signifies the accurate presentation of an adsorbent [31,35], whereas the adsorption capacity and removal effectiveness can be changed by operational settings. In this investigation, Equation (5) was utilized to assess the effectiveness of the produced adsorbent material using the distribution coefficient (K d ) for Pb 2+ and Cd 2+ . The values of K d for Pb 2+ and Cd 2+ increased with increasing pH to reach pH 6 and then decreased by further pH increasing, as shown in Figure 8. This is because higher pH activates more adsorption locations and makes extra metal ions accessible to transfer from the solution to the BI@MWCNTs' surface. Additionally, at each pH level, BI@MWCNTs demonstrated a more definite favorite for Pb 2+ over Cd 2+ . These can be attributed to Pb 2+ showed strong affinity for BI@MWCNTs in the solution containing two metals. Table 2 shows the distribution coefficient values of Pb 2+ are more than that of Cd 2+ , and this is showing that BI@MWCNTs exhibited selective adsorption for Pb 2+ over Cd 2+ according to distribution coefficient values onto BI@MWCNTs.   The optimization of the adsorption processes of Pb 2+ and Cd 2+ onto BI@MWCNTs, was accomplished by utilizing three factors, namely time (min), metal starting concentration (mg/L), and BI@MWCNT dosage (mg). The pH of the solutions was maintained at 5 during the whole process period. Figure 9 shows that the percentage of removal efficiency was increased with longer contact time before decreasing with initial metal concentrations from 50 to 150 mg/L. The percent removal then remained roughly constant. This observation showed that the metal ions were initially superficially adsorbed, and the rate of adsorption increased quickly. The metal ions began to adsorb into the porous structure of BI@MWCNTs and reached a constant value due to the external surface of the BI@MWCNTs becoming saturated, meaning that no further adsorption took place. The optimum values of the three parameters, as shown in Figure 9A-C, were 20 min contact time, 50 mg/ L of Pb 2+ and Cd 2+ concentrations, and 6 mg of BI@MWCNTs at a pH of 5.0 and temperature of 25 °C in order to completely remove Pb 2+ . However, the optimum removal% of Pb 2+ and Cd 2+ was 100% and 98%, respectively, using 6 mg of BI@MWCNTs, pH 5.0, and an initial metal ion concentration 50 mg/L, and a contact time 20 min was accomplished. Using statistical analysis conditions, a confirmation test (three replicates) was conducted.

Pb 2+ and Cd 2+ Adsorption Optimization
The optimization of the adsorption processes of Pb 2+ and Cd 2+ onto BI@MWCNTs, was accomplished by utilizing three factors, namely time (min), metal starting concentration (mg/L), and BI@MWCNT dosage (mg). The pH of the solutions was maintained at 5 during the whole process period. Figure 9 shows that the percentage of removal efficiency was increased with longer contact time before decreasing with initial metal concentrations from 50 to 150 mg/L. The percent removal then remained roughly constant. This observation showed that the metal ions were initially superficially adsorbed, and the rate of adsorption increased quickly. The metal ions began to adsorb into the porous structure of BI@MWCNTs and reached a constant value due to the external surface of the BI@MWCNTs becoming saturated, meaning that no further adsorption took place. The optimum values of the three parameters, as shown in Figure 9A-C, were 20 min contact time, 50 mg/L of Pb 2+ and Cd 2+ concentrations, and 6 mg of BI@MWCNTs at a pH of 5.0 and temperature of 25 • C in order to completely remove Pb 2+ . However, the optimum removal% of Pb 2+ and Cd 2+ was 100% and 98%, respectively, using 6 mg of BI@MWCNTs, pH 5.0, and an initial metal ion concentration 50 mg/L, and a contact time 20 min was accomplished. Using statistical analysis conditions, a confirmation test (three replicates) was conducted.

Adsorption Kinetics
The influence of varying the adsorption time on the rate of the adsorption process was examined by applying two models of kinetics, namely the pseudo-first order and pseudo-second order models. These models were useful in understanding the Pb 2+ and Cd 2+ adsorption mechanisms using the prepared BI@MWCNTs at the adsorbed initial concentration (50 mg/L). The adsorption capacity was thought to be controlled by a single mechanism performing on a single class of adsorbing sites according to the pseudofirst order theory [28]. The calculated values were found to be very different from their similar experimental values when the calculated values were compared to the experimental adsorption capacities. Consequently, the pseudo-first order reaction kinetic model could not effectively characterize Pb 2+ and Cd 2+ adsorption on BI@MWCNTs. Alternatively, Figure 10 shows that the pseudo-second order kinetic model was reliant on the quantity of adsorbed adsorbate on the surface of the BI@MWCNTs. The agreement between the calculated (q e,cal ) and experimental q e values is shown in Table 3. As demonstrated in Table 3, the calculated q e (q e,cal ) values agreed with the experimental q e values (q e,exp ), and the correlation coefficient (R 2 ) values that were attained were consistently higher than those of the pseudo-first order. According to the obtained results, chemisorption was the rate-determining step in the adsorption procedure, and the adsorption rate for both ions was influenced by the approachability of adsorption locations on the surface of the BI@MWCNTs. Furthermore, k 2 values decrease with the increase in the initial concentrations of Pb 2+ or Cd 2+ , which may be due to higher competition for the adsorption sites at high concentrations compared with low concentrations [36].

Adsorption Kinetics
The influence of varying the adsorption time on the rate of the adsorption process was examined by applying two models of kinetics, namely the pseudo-first order and pseudo-second order models. These models were useful in understanding the Pb 2+ and Cd 2+ adsorption mechanisms using the prepared BI@MWCNTs at the adsorbed initial concentration (50 mg/L). The adsorption capacity was thought to be controlled by a single mechanism performing on a single class of adsorbing sites according to the pseudo-first order theory [28]. The calculated values were found to be very different from their similar experimental values when the calculated values were compared to the experimental adsorption capacities. Consequently, the pseudo-first order reaction kinetic model could not effectively characterize Pb 2+ and Cd 2+ adsorption on BI@MWCNTs. Alternatively, Figure  10 shows that the pseudo-second order kinetic model was reliant on the quantity of adsorbed adsorbate on the surface of the BI@MWCNTs. The agreement between the calculated (qe,cal) and experimental qe values is shown in Table 3. As demonstrated in Table 3, the calculated qe (qe,cal) values agreed with the experimental qe values (qe,exp), and the correlation coefficient (R 2 ) values that were attained were consistently higher than those of the pseudo-first order. According to the obtained results, chemisorption was the rate-determining step in the adsorption procedure, and the adsorption rate for both ions was influenced by the approachability of adsorption locations on the surface of the BI@MWCNTs. Furthermore, k2 values decrease with the increase in the initial concentrations of Pb 2+ or Cd 2+ , which may be due to higher competition for the adsorption sites at high concentrations compared with low concentrations [36].

Adsorption Isotherms
As shown in Table 4, using the Langmuir and Freundlich adsorption isotherm models at dissimilar initial concentrations and two altered temperatures (25 • C and 35 • C), the data revealed that the Langmuir model, which suggests that Pb 2+ and Cd 2+ ions adsorb with monolayer coverage, best represented the adsorption data. The determined value q m of Cd 2+ adsorbed on the BI@MWCNTs was 916 mg/g, whereas the determined value q m of Pb 2+ adsorbed on the BI@MWCNTs was 931 mg/g, which was noticed to be increased at high temperatures. In addition, the Freundlich model's 1/n F values, which indicate surface heterogeneity or adsorption intensity, varied between 0 and 1, revealing a larger heterogeneity as the values move toward 0. The Pb 2+ and Cd 2+ ions adsorbed onto the BI@MWCNTs primarily via the internal surface by regulatory intra-particle diffusion, which is consistent with the conclusion that a heterogeneous chemisorption process took place. The findings show that both metals can effectively adsorb using the Langmuir and Freundlich isotherm models.  Table 5 shows the thermodynamic factors that were investigated at 25 • C, including Gibbs free energy (∆G), enthalpy (∆H), and entropy (∆S). The data show that the negative values of ∆G indicate that the adsorption system on the BI@MWCNTs was physisorption and a spontaneous process [37]. On the other hand, the positive values of ∆H and ∆S elucidate that the adsorption processes on the BI@MWCNTs were random and endothermic processes, respectively [27].

Regeneration Adsorption Tests
After undergoing the adsorption process, 4 mg of the used BI@MWCNTs was spread in 0.1 M HCl solution and sonicated for 30 min with the intention of assessing its ability to be used again. The BI@MWCNTs were then cleaned and reused. Figure 11 shows the percentage of removal as a motionless high after five cycles, and the differences were relatively minor. This is a strong indication that BI@MWCNTs are an excellent and costeffective adsorbent for both lead and cadmium elimination in wastewater applications.  Table 5 shows the thermodynamic factors that were investigated at 25 °C, including Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS). The data show that the negative values of ΔG indicate that the adsorption system on the BI@MWCNTs was physisorption and a spontaneous process [37]. On the other hand, the positive values of ΔH and ΔS elucidate that the adsorption processes on the BI@MWCNTs were random and endothermic processes, respectively [27].

Regeneration Adsorption Tests
After undergoing the adsorption process, 4 mg of the used BI@MWCNTs was spread in 0.1 M HCl solution and sonicated for 30 min with the intention of assessing its ability to be used again. The BI@MWCNTs were then cleaned and reused. Figure 11 shows the percentage of removal as a motionless high after five cycles, and the differences were relatively minor. This is a strong indication that BI@MWCNTs are an excellent and cost-effective adsorbent for both lead and cadmium elimination in wastewater applications.  Figure 11. Pb 2+ and Cd 2+ adsorption-desorption cycles on BI@MWCNTs (adsorbed initial concentration, 50 mg/L; adsorbent dose, 4 mg/L; pH, 5.0; temperature, 25 • C). Table 6 shows comparative values of the previously reported percentages of removal for Pb 2+ and Cd 2+ from solutions using different carbon nanotube-, carbon-, and graphene oxide-based materials with the prepared BI@MWCNTs in this study. As illustrated, the percentage of removal, using the synthesized BI@MWCNTs in this study for (Pb 2+ ) and (Cd 2+ ) are 100 and 98%, respectively, which are excellent values for the elimination of both metal cations compared to the previously reported carbon-based adsorbents. Furthermore, the prepared functionalized MWCNTs in this study (BI@MWCNTs) have been recovered and reused six times, which makes them a cost-effective and efficient absorbent for the removal of (Pb 2+ ) and (Cd 2+ ) from wastewater.

Conclusions
A well-known functionalized benzimidazole was used to modify acid functionalized multiwalled-carbon nanotubes (MWCNTs-CO 2 H) and yield a novel functionalized nanomaterial (BI@MWCNTs). The structure and physical characteristics of the BI@MWCNTs were confirmed and measured by surface area analysis, dynamic light scattering (DLS) analysis, Raman spectroscopy, scanning electron microscopy (SEM), and FT-IR. The main conclusions of this study can be discussed in the following points: Effect of pH-the influence of the pH on the adsorption effectiveness of BI@MWCNTs toward Pb 2+ and Cd 2+ was investigated in the pH range of 2.5 to 10. The optimum efficiency for removal of both ions was achieved at pH 5. BI@MWCNTs exhibited selective adsorption for Pb 2+ over Cd 2+ according to distribution coefficient values. The optimization of the adsorption processes of Pb 2+ and Cd 2+ onto BI@MWCNTs was accomplished by utilizing three factors, namely time (min), metal starting concentration (mg/L), and BI@MWCNT dosage (mg). The optimum values of the three parameters are 20 min contact time, 50 mg/L of Pb 2+ , and Cd 2+ concentrations, and 6 mg of BI@MWCNTs, carried out at a pH of 5.0 and a temperature of 25 • C in order to completely remove Pb 2+ (100%) and 98% of Cd 2+ . In terms of adsorption kinetics, the influence of varying the adsorption time on the rate of the adsorption process was examined. According to the obtained results, chemisorption was the rate determining step in the adsorption procedure, and the adsorption rate for both ions was influenced by the approachability of adsorption locations on the surface of the BI@MWCNTs. The pseudo-second order kinetic model was reliant on the quantity of adsorbed adsorbate on the surface of the BI@MWCNTs. In terms of adsorption isotherms, using the Langmuir and Freundlich adsorption isotherm models at dissimilar initial concentrations and two altered temperatures (25 • C and 35 • C), the data revealed that the Langmuir model, which suggests that Pb 2+ and Cd 2+ ions adsorb with monolayer coverage on the BI@MWCNTs' surface, best represented the adsorption data. The thermodynamic factors that were investigated at 25 • C, including Gibbs free energy (∆G), enthalpy (∆H), and entropy (∆S), gave the following results. The Pb 2+ and Cd 2+ adsorbent systems were shown to be spontaneous by the negative values of ∆G, while the positive values of ∆S show that randomness increases at the solid-liquid interface during the adsorption process, and the positive values of ∆H show that the adsorption route is endothermic. Regeneration adsorption tests have revealed that the percentage of removal as a motionless high after five cycles, and the differences were relatively minor. The comparative study revealed that the percentage of removal using the prepared BI@MWCNTs for (Pb 2+ ) and (Cd 2+ ) are 100 and 98%, respectively, which are excellent values for the elimination of both metal cations compared to the previously reported carbon-based adsorbents. Furthermore, the prepared functionalized MWCNTs in this study (BI@MWCNTs) have been recovered and reused six times, which make them a cost-effective and efficient absorbent for the removal of (Pb 2+ ) and (Cd 2+ ) from wastewater.