Experimental and computational studies on activated Bambara groundnut (Vigna subterranean) hulls for the adsorptive removal of herbicides from aqueous solution

The excessive usage of herbicides to control herbs by farmers has become an issue of interest to the environmentalist due to the threat posed by such act to the ecosystem, and therefore, there is the need to curb such practices. The efficiency of Bambara groundnut shell for the removal of pendimethalin (PE) and paraquat dichloride (PQ) herbicides from aqueous solution was established in this work. The activated carbon was prepared via chemical activation method using trioxonitrate (v) acid by determining its void volume, moisture content, bulk density and Fourier transform infrared (FTIR) and scanning electron microscopy (SEM) methods. Batch adsorption techniques were set to optimize the adsorption parameters such as solution pH, adsorbate concentration, contact time, adsorbent dosage and temperature in order to depict the best optimal conditions for the adsorption process. The adsorption process was examined in terms of its equilibrium data, kinetics, thermodynamics involved in the adsorption process as well as computational quantum chemical parameters evaluation. The batch adsorption experiments revealed that the amounts of PE and PQ adsorbed were found to vary with the contact time, adsorbent dosage, pH and initial concentration. The adsorption of PE and PQ decreased with increasing adsorbent dose but increases with increasing initial concentration of the PE and PQ solution. Isotherm studies revealed that the equilibrium data fitted to both Langmuir and Freundlich model with R-squared values of 0.976, 0.993 and 0.909, 0.978 for PE and PQ, respectively, which implied that Langmuir isotherm had a better fit. This was also found to be an indication that the uptake of PE and PQ by ACBGNS occurred through monolayer adsorption on identical homogenous sites. Also, kinetic modeling results obtained showed that the pseudo-second-order model explained the adsorption kinetics of PE and PQ by ACBGNS best, which meant that chemisorption was the slowest step and, thus, the rate determining step. The positive value of ΔH and the positive value of ΔG show the endothermic and spontaneous nature of adsorption of PE and PQ ACBGNS. Batch adsorption experiment and characterization of the ACBGNS have indicated that Bambara ground nut shell can be used to produce activated carbon that can be applied effectively for adsorption of PE and PQ from an aqueous solution. Computational studies results obtained from quantum chemical analysis are consistent with the experimental results obtained from this study.


Background
Farmers in developing countries like Nigeria are faced with various threats posed by unregulated and excessive usage of herbicides during farming processes which results in serious contamination and pollution of surface and ground waters (Neera 2009). These acts call for serious concern by environmentalists and their likes to curb such practices of the reckless use of these chemicals in crop farms, orchards, fields and forest lands (Neera 2009;Salman et al. 2011). Scientists have used several methods to control the unguided usage of these chemicals before, during and after farm practices of which adsorption finds an application. However, activated carbon usage in waste water treatments has been established to be more effective and efficient when compared to other adsorbent materials. This is due to its ease in experimental design and usage even though activated carbons which are commercially available are quite expensive (Ali et al. 2016;El-Shamy et al. 2018;Yanyan et al. 2018;Wei et al. 2019;Khan 2020;Alzhan et al. 2020;Liliana et al. 2020;Ahmed et al. 2020;Ahmed and Sayed 2020;Abdelhamid et al. 2020;Khalil and Kenawy 2020;El-Azazy et al. 2021). Reported studies on the use of Bambara groundnut shell do not involve using it to remove herbicides from an aqueous solution which could have led to understanding of the mechanism of the process (Ademola et al. 2019).
Pendimethalin was known as (3,4-dimethyl-2,6-dinitro-N-pentan-3-ylaniline) and that of paraquat dichloride (1,1-dimethyl-4,4-bipyridinium chloride). Pendimethalin and paraquat dichloride for their interaction with adsorbent surface. However, there is paucity in the usage of activated carbon sourced from Bambara groundnut shells (BGNS) for the adsorption of pendimethalin and paraquat dichloride. Bambara groundnut is an edible fruit that is cultivated mainly because of its high nutritional value. In northern Nigeria, large quantities of Bambara groundnut hulls are usually disposed as waste materials (Akinola et al. 2019). BGNS can be easily transformed and utilized as a good adsorbent. Most of the reported work on the use of BGNS product has to do with the ability of these adsorbents in adsorbing dyes (Akinola et al. 2019). There are few reports on the specific mode of action of the adsorbent and the extent of its adsorption (Akinola et al. 2019;Ayuba and Nyijime 2020).
In describing reactivity and illustrating chemical, molecular and electronic structure of compounds, quantum chemical methods have been found to find a lot of application (Kraka and Cremer 2000). This can assess the potentials and efficacy of the interaction between the adsorbent and adsorbate through mathematical intuitions and modeling. The correlation between structure and reactivity could be used to even predict the efficiency of compounds in adsorption which are yet to be synthesized (Kareison and Lebanov 1996). The use of equations of mathematical origin and models could be used to establish and predict chemical phenomena of which adsorption is such an example (Hinchliffe 1994(Hinchliffe , 1999. In this previous study, it was shown that the activated Bambara groundnut shell is established to be an efficient adsorbent which can be applied for the adsorption of both pendimethalin and paraquat dichloride herbicides form aqueous solution. It functions through surface adsorption sites which depend on experimental parameters including solution pH, amount of adsorbent, time of contact and initial concentration of the herbicide. Physical process favors the adsorption of both PE and PQ herbicides on Bambara groundnut shell. Relatively based on the experimental data generated, the isotherm that fitted best was that of Langmuir model. The kinetics that best describes the adsorption process of both herbicides to the ACBGNS was that of second order (pseudo). Evaluated parameters obtained for both adsorption processes (ΔG, ΔH and ΔS) describe the process to be spontaneous and feasible in nature.
Furthermore, this approach can harness the structural and electronic parameters of these molecules in relation to adsorption to propose a suitable mechanism which may not be achieved experimentally (Growcock 1989;Wang et al. 1999). The aim of the present study was to investigate and established the adsorption efficiency of activated Bambara groundnut shell in the removal of pendimethalin (PE) and paraquat dichloride (PQ) herbicides from aqueous solution.

Preparation of adsorbent
The BGNS sourced after removing the seeds from the pods were thoroughly washed to clear off impurities using water. Shells obtained were subjected to drying at 80 0 C to obtain constant weight. Pulverization and subsequent sieving of the dried shells into fine particles were conducted as described by Ayuba and Nyijime (2021). The final product was kept for further usage in an air safe container. 100 g of the fine powder was mixed with 10% phosphoric acid (H 3 PO 4 ) (100 g sample + 100 mg/l of 10% H 3 PO 4 , wt/v). Pyrolysis in an electric muffle furnace of the fine particles was conducted for 1 h at 450 °C and subsequently removed and cooled at room temperature. Before drying at 110 °C at a period of 24 h, the pyrolyzed carbon was subjected to washing using 2% HCl and distilled water to achieve a neutral pH sample (Ayuba and Nyijime 2021).

Preparation of herbicides solutions
The stock herbicides solution was prepared by mixing 3.6 ml of the 276 mg/l of PQ solution and 2 ml of 500 mg/l of PE solution separately each to a 1-L flask and further diluting to the mark. Serial dilution was carried out using distilled water to give solution concentration of 10, 20, 30, 40, 50 and 60 mg/l of each herbicides solution.

Batch adsorption experiments
Adsorption experiments for PE or PQ onto the prepared adsorbent were studied by weighing 0.1 g of the adsorbent (ACBGNS), and 10cm 3 of 50 ml initial concentration (PQ or PE) was introduced into a 250-ml beaker and was tightly covered and placed in an incubator shaker (Innova model) and run at 25 °C for 2 h. After equilibration, samples were withdrawn and filtered through Whatmann filter paper and the absorbance of the residual herbicides was measured using a Uv-visible spectrophotometer (PerkinElemer) at λ max 243 and 257 nm for pendimethalin and paraquat dichloride, respectively. The same method was used while varying the initial concentration, the contact time, the adsorbent dosage and temperature of adsorption during optimization experiments. The optimized values were used to conduct the main batch adsorption experiments at varying temperatures (60 mg/l concentration, 0.2 g dosage, pH 6 and pH 8, 80 min contact time for PE and PQ, respectively). The amount of the herbicide adsorbed for each gram was obtained using Eq. (1) (Fouodouo et al. 2017).
The initial and equilibrium concentration (mg/l) are defined as C o and C e , volume of the herbicide solution (L) is V, and the mass (g) of ACBGNS is m.

Characterization of the adsorbent
The adsorbent (ACBGNS) was characterized using Cary 630, Agilent Technology model of Fourier Transform Infrared Spectrophotometer and Scanning Electron Microscope (Phenom World Eindhoven). The conducted FTIR experiment used a range of wave numbers (650-4000 cm −1 ) with 32 scans at 8 cm −1 resolution, while the SEM analysis was performed by scanning the micrographs of ACBGNS before and after adsorption using 15.00 kV voltage of acceleration with an × 500 magnification.

Computational details
All the geometries of the molecules have been fully optimized using the B3LYP method within the framework of density functional theory (DFT) housed in Materials Studio 7.0 software (Accelrys, Inc.) in conjunction with the double numerical polarization (DNP) basis set (Lee et al. 1988;Becke 1993). After obtaining the geometry optimized herbicide molecules, quantum chemical parameters were computed which were later used to predict the adsorption mechanisms (Tao et al. 2010;Mohammed et al. 2017;Zohdy et al. 2021). This method has been used to provide useful information about useful quantum chemical parameters (Khaled et al. 2005;Emreguland and Hayvali 2006;Zohdy et al. 2021), through density functional theory (DFT) to predict the structural nature of the adsorbate and adsorbent interaction processes (Lashkari and Arshadi 2004;Mohan and Joseph 2018). The electron affinity (EA) and ionization potential (IP) of the neutral herbicide molecules were calculated in relation to the lowest energy of the unoccupied molecular orbital (E LUMO ) and the highest occupied molecular orbital energy (E HOMO ), respectively, using Koopman's theory (1993).
Energy gas ∆E was calculated using Eq. (4): Calculated ionization potential and electron affinity values were used to determine the values of the absolute electronegativity χ, the absolute hardness (η) and the softness (σ) using expressions in Eqs. (5) and (6): (Pearson 1963).
Value of the global softness (σ) is the reciprocal of the global hardness (Pearson 1963). (2) Page 4 of 13 Nyijime et al. Bull Natl Res Cent (2021) 45:189 For a reaction made of adsorbent and adsorbate having varying electro negativities, the flow of electrons is expected from the adsorbate molecule having lower electronegativity toward that with a higher negativity (adsorbent surface) till equilibrium of the chemical potentials of both systems is attained (Parr and Szentpaly 1999).
Fukui indices were used to establish the reactive sites on the herbicide (adsorbate) which shows chemically reactive regions of the electrophilic and nucleophilic character of the studied herbicide molecules (Verma et al. 2018). Chemically, where the function of the Fukui is found to have large values predicts a region of higher reactivity than where the value is small. To further predict the sites which can serve as hard and soft acid and bases, the HSAB principle can be invoked in relation to the evaluated Fukui indices. To define Fukui function f(r) for nucleophilic attack (Eq. 8) and electrophilic attack (Eq. 9) depending on the direction of electron transfer, the scheme of finite difference approximations obtained from Mulliken population analysis of atoms in the herbicide molecules was used (Alibakhshi et al. 2018;Verma et al. 2018;Zohdy et al. 2021).
Value of q k is the molecule's gross atom charge, the number of electrons (N) in the molecule, the anion is N + 1 describing the LUMO of the neutral molecule with an added electron, and the cation is defined with respect to N − 1 showing a new HOMO after electron removal from the neutral molecule. The ground state geometry of the herbicide molecules resulting from local minima was used for all calculations.

Results
The spectroscopic and micrograph results of this study are reported in Figs. 1 and 2, respectively. The optimized adsorption density of various adsorption parameters including solution pH, contact time, adsorbent dosage, initial herbicide concentration is reported in Figs. 3, 4, 5 and 6. The calculated adsorption isotherm, kinetics study, thermodynamic as well as the calculated quantum chemical parameters and Fukui indices of pendimethalin and paraquat dichloride molecule are reported in Tables 1, 2, 3, 4 and 5, respectively. Figure 7 shows the optimized geometry, HOMO, LUMO of PE and PQ molecules.

Spectroscopic characterization of the adsorbent
FT-IR spectra of pendimethalin and paraquat dichloride (Fig. 1a-c) before and after adsorption show five important peaks which include O-H broad peak at 3117 cm −1 , C ≡ C vibration at 2322 cm −1 , 2084 cm −1 and 2113 cm −1 , peaks at 1994 cm −1 and 1871 cm −1 which can be assigned to C=O group; such assignments can be affected by neighboring aromatic C-C ring stretched vibrations. At 1562 and 1063 cm −1 , the following stretch assignments of C=C and C-O were observed . The peaks at 970 cm −1 , 873 cm −1 and 676 cm −1 correspond to C-H in alkenes and C-H in aromatic rings. There was a minor difference that was observed prior to and after adsorption, which resulted in a shift and broadening of adsorption peaks. The shift of the -OH peak from 3117 to 3339 cm −1 for pendimethalin and 3117 cm −1 to 3119 cm −1 for paraquat dichloride indicates the involvement of the hydroxyl groups in the adsorption of the herbicides (Abbas and Ahmed 2016; Jawad et al. 2018).
The micrograph of the surface of ACBGNS before adsorption (Fig. 2a) shows the morphology of the surface adsorbent having pores, cracks, holes which are expected to enhance adsorption. After adsorption of the two herbicides differently (Fig. 2b, c), the surfaces showed the formation of clusters, patches, with filled cracks on the ACBGNS surfaces which shows the participation of ACBGNS in the adsorption of PE and PQ.
Some of the obtained physical properties of the activated carbon produced 11.2 moisture content (%), 0.242 density (g/cm 3 ) and 1.78 pore volume (cm 3 ). The value of the moisture content pore volume and the bulk density of the produced activated carbon revealed that it had good adsorptive properties. It was noticed from the properties of the activated carbon that though it might not give up to 100% adsorption, it will be good for adsorption of organic and inorganic materials to a large extent (Fito et al. 2019). Maximum weight loss was observed in the present study; 51.82% of the activated carbon was obtained after activation. Figure 3 shows the variation in PE and PQ amount adsorbed onto ACBGNS with solution pH. Removal of PE and PQ by ACBGNS is observed to be at optimum pH of 6 and 8. It can be observed from the plots that ACB-GNS posses high removal of PE at lower pH, while PQ was more adsorbed at higher pH condition. The presence of pectin on the surface of ACBGNS favored the adsorption of PQ which is cationic in nature; this is as a result of negatively charged ions increase in alkaline medium Page 5 of 13 Nyijime et al. Bull Natl Res Cent (2021) (Bello et al. 2017), while the surface of ACBGNS which also contains phenolic or carboxylic functional groups may be the reason for increased adsorption of PE at low pH (Bedin et al. 2015).

Effect of PE and PQ adsorption duration
Variation in the amount of PE and PQ molecules adsorbed with respect to duration is presented in Fig. 4. It was found that the adsorption rate of PE and PQ by ACBGNS was initially fast reaching equilibrium within 80 min. PE and PQ molecules were firmly attached onto the ACBGNS resulting in a higher rate of adsorption. At Page 7 of 13 Nyijime et al. Bull Natl Res Cent (2021) 45:189 a later stage, the rate of PE and PQ removal was found to be constant. The available surfaces for binding the molecules became almost exhausted, and therefore, repulsive forces either between PE or PQ with the ACBGNS surface made the herbicide molecules difficult to occupy the vacant sites (Mondal and Kar 2018). Figure 5 reveals the variation in equilibrium amount of PE and PQ adsorbed with adsorbent dosage (0.2-1.2 g). From the plots, the amount of PE and PQ adsorbed increased with an increase in the adsorbent amount used, while the amount adsorbed for a given amount of the adsorbent decreased considerably for both PE and PQ. This decrease observed for the herbicides in terms of adsorption per unit mass with respect to increasing adsorbent dosage could be attributed to possible overlapping of adsorption sites with an increase in adsorbent dosage. This is an indication that ACBGNS has a high adsorption capacity at low dosage. Figure 6 presents an increase in amount per gram of PE and PQ adsorbed as their initial concentration increases. Driving forces of molecules in solution are low at a lower concentration, therefore resulting in a lower adsorption rate, while when concentration increases, it is expected that the PE and PQ would have higher driving forces which enable them to easily occupy the active sites on the ACBGNS surface from the aqueous phase. This was reported by other researchers (Mondal and Kar 2018;Lee et al. 2019).

Adsorption isotherm
Three different isotherm models, Temkin, Freundlich and Langmuir, were used to treat the generated isotherm data with the intention of proposing adsorption process mechanism. The Langmuir's model has been used to describe the adsorption of monolayer which occurs on the outer surface of an adsorbent with infinite number of active adsorption sites.
C e is the concentration of PE or PQ in mg/L at equilibrium, q e , the amount of herbicide adsorbed per gram of the adsorbent at equilibrium (mg/g), Q o is maximum monolayer coverage capacity (mg/g), K L is the isotherm constant (L/mg) for Langmuir. From the slope and intercept, the values of q max and K L were evaluated through extrapolation of the Langmuir plot of 1 qe versus 1 Ce . The  K F in the Freundlich's equation describes adsorption capacity constant. The adsorption intensity's empirical parameter can be defined by 1/n. From this study, it is clear that the capacity adsorption can depend on major factors including the nature of adsorbent and the type of herbicide used. The calculated values of R L for the two herbicides are given in Table 1 which shows PE and PQ adsorption onto ACBNGS surface is favorable. As tabulated in table, values of n ranged between 0 and 1 qualifying a favorable adsorption.
When the surface coverage is related to the adsorption free energy during the interaction of the herbicides with the adsorbent surface, we consider the isotherm described by the Temkin model. In addition, the heat of adsorption associated with the herbicide molecules on the layer is found to decrease proportionately because of adsorbent-adsorbate interaction. This type of adsorbate-adsorbent interaction is associated with uniformly associated binding energies, to reach a maximum value (Al-Saad et al. 2019). This statement is as presented by the following relation in Eq. (12): The Temkin constant B = RT/b is related to the heat of adsorption, the molar gas constant is R (J/mol K), the temperature is T(K), the adsorption energy variation is     Page 10 of 13 Nyijime et al. Bull Natl Res Cent (2021) 45:189 regression equation and R 2 value observed for this model favorably described the adsorption process of the herbicides onto ACBGNS, across the experimental variables applied. When the values of the correlation coefficient, R 2 , were compared for the three tested isotherms, relatively, the Langmuir model fitted better among others for the adsorption data of herbicide onto ACBGNS.

Kinetic study
The kinetic studies of systems involved in adsorption describe how fast the adsorbate uptake onto the adsorbent is which controls the amount of equilibration time.
The adsorption kinetic curves of the PE and PQ molecules on ACBGNS were obtained by subjecting the experimental data obtained to three models of kinetics comprising pseudo-first and second-order and intraparticle diffusion.

(a) Kinetic model of the pseudo-first-order
Equation (13) can be used to define the model (Benkaddour et al. 2018).
k 1 is the rate constant and q e is the amounts (mg/g) of PQ while q t is that of PE adsorbed at equilibrium and time t, respectively. log (q e − q t ) was plotted against t, and values associated with k 1 were evaluated from the graph's slope. R 2 values, calculated and experimental q e values obtained for PE and PQ adsorption onto ACBGNS are presented in Table 2. R 2 was observed to be very low coupled with a reasonable difference established between the adsorption capacity (q e ) values for experimental and calculated, showing failure in experimental data interpretation.

(b) Kinetic model of the pseudo-second-order
This type of second order has been widely to predict kinetics of adsorption processes and the expression is as in Eq. (14) (Hameed et al. 2009).
Adsorption rate constant is given by k 2 (g/mg min), while t/q t values were plotted against t (figure not shown).
According to the kinetic model of the pseudo-second order, parameters of adsorption of herbicides onto ACB-GNS are given in Table 2. From the results, parameters obtained were found to fit the equation with the R 2 values close to unity. There was good correlation between the calculated and experimental values. The use of this model in adsorption process describing the adsorption of PE and PQ herbicides onto ACBGNS was proved based on such comparison.

(iii) Intraparticle diffusion model
Weber-Moris Eq. (15) predicts the possibility of the sole mechanism through the intra particle diffusion model (Raul et al. 2011).
The constant k int (mg/g min 0.5) describes rate of diffusion, while the boundary layer thickness is defined by C. If the rate-limiting step is only through intra particle diffusion, then a linear plot would be obtained when q t versus t1/2 is extrapolated to the origin. If not, there is involvement of other mechanism. The values of the constant C from Table 2 were found to 15.74 and 11.28 for PE and PQ, respectively. It can be concluded that PQ and PE adsorption by the ACBGNS cannot be predicted by this model.

Thermodynamic adsorption study
Notably parameters obtained from thermodynamic evaluations give vital information about the adsorption nature with respect to feasibility and spontaneity of the process. Temperature effects associated with PQ and PE adsorption onto ACBGNS were studied by temperature variation from 303 to 333 K while keeping all other parameters (pH = 6 and 8, respectively, for PE and PE, adsorbent dosage 0.2 g, initial PQ concentration = 60 mg/l, t = 80 min) at optimized values. The thermodynamic constants for adsorption which can be used to describe the spontaneity and feasibility of an adsorption process are described entropy (∆S), free energy (∆G) and enthalpy (∆H) changes. Equations (16) and (17) describe them thus: the temperature being T(K), the gas constant, R, 8.314 J/mol K and k c the constant of equilibrium which can be obtained from the relation in Eq. (18) (Calvete et al. 2010).
where C a is the adsorbed herbicide in mg/L and the equilibrium herbicide concentration is C e (mg/L). The slope from the van't Hoff 's plot and intercept and the values of ∆S and ∆H were evaluated, and Table 3 contains all the (15) q e = C + K t1/2 int 13 Nyijime et al. Bull Natl Res Cent (2021) 45:189 calculated values. ΔH values which were found to be positive show that an endothermic process is involved fort the adsorption of the herbicides onto ACBGNS. The positive value of ΔS indicates high randomness at PQ and PE molecules adsorption interphase onto ACBGNS (Gunay et al. 2007). Physical adsorption has values in the range of 2.1-20.9 kJ/mol, which describes the type of adsorption for both herbicides involved in this research, while higher range of 80-200 kJ/mol can be used to qualify chemisorption mechanism. ΔG values which are negative show the adsorption is thermodynamically feasible and spontaneous (Mustapha et al. 2019).

Computational study
Quantum chemical parameters of PE and PQ molecules calculated are presented in Table 4 Table 4 describe that PQ herbicide adsorption is higher than that of PE which may be due to the presence of chloride ion in solution and lone pairs of nitrogen increase its activity which is in good agreement with that earlier reported here experimentally. Relationship between adsorption and dipole moment is yet to be fully established in literature with most adsorption processes reported to favor lower μ as shown in Table 4. The μ for PQ is lower than that of PE and therefore better adsorbed. The adsorption ability of a chemical compound can be controlled by its binding energy (E B ), surface volume (SV), surface area (SA). The higher the values of SV, SA, E B obtained, then the process of protonation, hydration and subsequent adsorption on a surface for a given molecule becomes more difficult. Therefore, lower values of these parameters associated with PQ herbicide show that it is better protonated, hydrated and adsorbed on ACBGNS than PE herbicide.
In the presence of charge atmosphere, polarizability can be defined as a measure of the behavior of molecular species. Tendency of the ease of a molecule to be adsorbed is related to better response. This explains why PQ herbicide is better adsorbed on ACBGNS than PE herbicide. Better adsorption of PQ is also favored by their SV and SA values which are low. It is established that molecules compete with difficulty to fit into the adsorption matrix of the adsorbent when the values of surface area and volume of the adsorbate are large.
An index, the function of Fukui can be used to assess the local reactivity (hence adsorption) of a molecule on the basis of the constituent atoms which are electronic in nature (Alokdut et al. 2017). Based on electrophilic and nucleophilic nature of the chemical reactivity of a molecule, Fukui indices provide useful information about the reactive centers (Mohan and Joseph 2018;Zohdy et al. 2021). Threshold values of f + k andf − k control nucleophilic and electrophilic centers of in a given molecule. PE and PQ herbicide's Fukui functions are shown in Table 5. F − for PE and PQ molecule have their highest Mulliken and Hirshfeld charges on N(1) and H(13), C(5) while for the nucleophilic (F + ), PE herbicide on O(17), O(18) and PQ herbicide are on the aromatic ring N(5) shown in Table 5. Therefore, the adsorption of PE on ACBGNS occurs through the nitrogen atom while that of PQ occurs via the aromatic ring of the molecule. Optimized geometry, Frontier orbital diagrams of PE and PQ herbicides are reported in Fig. 7. Within the diagrams, red loop shows positive sites while negative charges are displayed in blue loops on the molecules.

Conclusions
The present study has demonstrated that activated Bambara groundnut shell can be used as a natural, promising, economic and environmentally friendly adsorbent for PE and PQ removal from aqueous solutions. FTIR, SEM, moisture content, pore volume and bulky density were used to analyze for the functional groups and surface area that are responsible for the herbicide ion removal from aqueous solutions. The batch adsorption method is dependent on the contact time, dosage, solution pH and initial herbicide ion concentrations. The experimental data in the adsorption process indicated good correlations with the pseudo-second-order kinetic model and Langmuir isotherm. The thermodynamic study showed that the process is highly endothermic in nature, with positive value of ∆H and change in ∆S. Based on this study, the readily available Bambara groundnut shell which is cheaply available to the environment could be used as a natural adsorbent, due to its high removal efficiencies of the herbicide ion involved in this study. Studied computational parameters and that of the experimental supported each other with high degree of similarities.