# Thermodynamic, Nonlinear Kinetic, and Isotherm Studies of Bisphenol A Uptake onto Chemically Activated Carbons Derived from Safou (Dacryodes edulis) Seeds

Authors: Cyrille Donlifack Atemkeng; Arnaud Kamdem Tamo; Giscard Doungmo; Liouna Adoum Amola; Julio Jimmy Kouanang Ngouoko; Théophile Kamgaing
Journal: Journal of Chemistry
 (2022)

---

## Abstract

The interest of this work is to evaluate the possibility of using safou seeds to develop a new low-cost adsorbent and study its application to remove bisphenol A from an aqueous solution for a sustainable and ecological use of this biomass. This was done by optimizing some parameters that influence the adsorption process. The central composite design with four centre points was used to optimize the process variables. The concentration of bisphenol A solution, adsorbent dosage, stirring time, and solution pH on the adsorption capacity were considered, while the response measured was the quantity adsorbed. The activated carbon obtained by treatment with H2SO4 was named NSST and that obtained by treatment with H3PO4 was named NSPT. XRD revealed an amorphous character for the ACs, and EDXS showed they are mainly carbonaceous. Under the optimal adsorption conditions, NSPT showed the best performance. Correlation coefficients R2 and R2adj were of 85.13 and 69.12% for NSPT and 83.71 and 66.17% for NSST. A pseudo-second-order nonlinear kinetic model best described the adsorption kinetic of BPA removal by the ACs. Langmuir’s isotherm best described the adsorption of BPA onto both adsorbents. Thermodynamic studies suggested an exothermic and physisorption process.


---

## Body

## 1. Introduction

Water is the natural resource particularly used in our everyday life and mostly for industrial processes. After use, its release into the environment in more or less highly polluted residual form poses a serious threat to the aquatic fauna because it accumulates in their organisms due to their non- or low biodegradability. These chemicals include endocrine disruptor compounds among which is bisphenol A (BPA). The latter is classified as a priority pollutant by the European Union and the United States Environmental Protection Agency. Also known as 4,4’-(propane-2,2-diyl)diphenol, BPA is an emerging contaminant that is often detected in natural waters at levels of ng/L orμg/L, making it a trace organic contaminant. It is obtained by the reaction between two equivalents of phenol and one equivalent of acetone. This compound is widely applied to produce plastics, resins, and other materials [1]. Bisphenol A, found in tap water, surface water, air, and foodstuffs in considerable quantities [1–3], is now part of the so emerging pollutants. It is capable of mimicking the effects of estrogens that are linked to diabetes, cardiovascular disease, breast and prostate cancer, infertility with changes in sperm quality, and hyperactivity in children [4–6]. In light of the risks, the removal of this chemical from wastewater has attracted attention since the specific migration limit for BPA in food is 0.6 μg/g [7]. Moreover, the maximum concentration in the surface water is 12 μg/L in the USA and 21 μg/L in the Netherlands.Studies have shown that due to the anaerobic or hypoxic state of water, the degradation of bisphenol A in natural water is slow or does not even occur [8]. Therefore, the concentration of bisphenol A in surface water fluctuates greatly depending on the region and discharge time [9]. From the literature, several industrial processes are described to overcome the problem of water pollution by chemicals [10–12]. At a time when the world is facing a major energy challenge, these processes, which are for the most part energy-intensive, lead to the formation of secondary products, which are often toxic and also contribute to global warming. It thus becomes essential to turn to techniques that require renewable resources to the detriment of fossil resources. It is as a result of this that adsorption that can be done on several types of adsorbent supports [13–16] became the most coveted industrial depollution technique. Recently, safou seeds have shown interesting adsorption properties [17]. Adsorption has attracted further attention due to advantages such as ease of use, low energy consumption, high efficiency, cost-effectiveness, flexibility, and a wide variety of existing adsorbent systems [1, 3, 13, 18]. The latter is even more advantageous when the adsorbent support is available and abundant.Nowadays, the use of eco-friendly, low-price, nontoxic, renewable, and biodegradable adsorbents revealed exceptional and compelling sorption efficiency for pollutants [19]. Tropical Africa is full of potential valuable fruit species that are not only abundant but also available almost round the year and tend to generate large amounts of waste residue that can be used for activated carbon preparation. As an example, safou seed residues are widely distributed in Cameroon. Based on this observation, the context of this work is to evaluate the possibility of using safou seeds to develop a new low-cost adsorbent and study its application to remove bisphenol A from aqueous solution for a sustainable and ecological use of this biomass. The cultivation of safou is abundant in Cameroon, and after the consumption of the fruits, the seeds are thrown away as waste. However, they are rich in saturated fatty acids and unsaturated fatty acids [20], which give their molecular surface interesting properties for an efficient adsorption of organic pollutants. To the best of our knowledge based on the literature, safou seeds have not yet been the subject of BPA adsorption study using the response surface methodology (RSM).Since the ability of activated carbon to remove pollutants is severely affected by some factors such as temperature, sorbent dosage, contact time, initial concentration of pollutant, and pH, it is therefore essential to design a suitable process to achieve the maximum uptake capacity of pollutants by the adsorbent [21]. With the traditional methods, the effect of only one variable is studied at a time, and it does not investigate the interactive impact between the variables. Hence, it could not display an overall effect of variables on the appropriate response [22, 23]. The response surface methodology is one of the statistical tools used to rectify this issue. In that case, the RSM has been widely applied in chemical engineering and the environment field. Thus, the central composite design (CCD) under RSM can result in an empirical model due to combined effective parameters and corresponding interactions to optimize the adsorption process [13, 24, 25]. The aims of this research are as follows: (1) to apply ACs characterized by XRD, FTIR, EDXS, pH at the zero-point charge, and BET/BJH methods for an optimal adsorption of BPA in aqueous solution using CCD with four centre points; (2) to reveal the effects of four independent variables (contact time, sorbent dosage, concentration, and pH of pollutant solution) and their interactions on BPA uptake capacity; (3) to model some adsorption kinetics and isotherms by the nonlinear regression method; and (4) to study the effect of temperature on the adsorption process. The interest of this work is the efficient removal of organic pollutants by activated carbon based on agricultural by-products for a sustainable use of local biomass.


## 2. Materials and Methods

### 2.1. Materials

#### 2.1.1. Instrumentation

Determination of the functional groups present on the surfaces the materials was done using Fourier transform infrared (FTIR) spectrophotometry. The IR spectra of the ACs were recorded using a Genesis FTIRTM brand spectrophotometer ATI Mattson equipped with a DTGS (deuterated triglycine sulfate) detector in transmission mode between 4000 and 400 cm−1. The pH measurement was performed with a Mettler Toledo AG pH meter. For the determination of the crystalline or amorphous structure of the activated carbons, the diffraction patterns were recorded using an STOE Stadi-p device (STOE & Co. GmbH, Darmstadt, Germany) with a Cu Kα1 radiation (λ = 1.54056 Å, 40 kV, 30 mA) between 0 and 70°. With regard to energy-dispersive X-ray spectrometry (EDXS), a Euro Vector CHNS-O brand device (Euro EA 3000) made it possible to identify the different elements present in the activated carbons. The Brunauer–Emmett–Teller (BET) coupled with Barrett–Joyner–Halenda (BJH) method gave some textural properties of the ACs by N2 adsorption at 77.13 K using a micrometric sorptometer (Thermo Electron Corporation, Sorptomatic Advanced Data Processing). A UV Genesys 10S UV-Vis spectrophotometer of Thermo Scientific was used at the maximum wavelength of 276 nm to determine the residual concentration of BPA in solution.


#### 2.1.2. Chemicals

All the chemicals used were of analytical grade; bisphenol A with a purity of 97% was purchased from Sigma-Aldrich (Germany) and used as such. BPA stock solutions were prepared by dissolving 0.1 g BPA in 1 L containing less than 10% methanol. Test solutions were made by diluting the stock solutions. 0.10 mol/L NaOH and 0.10 mol/L HCl were used to adjust the pH of the solutions. Distilled water was used throughout this study. The chemical structure and some properties of bisphenol A are presented in Table1.Table 1
Chemical structure and some properties of BPA.
Chemical structureMolecular formulaC15H16O2Molecular weight (g/mol)228.29λmax (nm)276SolubilitySoluble in alcohol and slightly soluble in CCl4


#### 2.1.3. Adsorbents

Two activated carbons (ACs) prepared by chemical means using H3PO4 and H2SO4 were used in this work. They were labeled as follows: NSPT for that activated with H3PO4, and NSST for that activated with H2SO4. The preparation of the ACs was made using CCD, and the variables whose effects were simultaneously observed are the concentration of the oxidizing agent (varying between 1 and 3 mol/L), the calcination temperature (between 400 and 700°C), and the calcination time in the oven (between 30 and 120 minutes) as described by Atemkeng et al. [17].


#### 2.1.4. Characterization of the Adsorbent

Knowledge of some physicochemical properties of an adsorbent material is useful and contributes to the understanding of phenomena, e.g., in the adsorption process. The ACs have undergone few characterizations such as X-ray powder diffraction (XRD), Fourier transform infrared (FTIR) spectrophotometry, energy-dispersive X-ray spectrometry (EDXS), pH at the zero-point charge, and the BET/BJH method.


## 2.1. Materials

### 2.1.1. Instrumentation

Determination of the functional groups present on the surfaces the materials was done using Fourier transform infrared (FTIR) spectrophotometry. The IR spectra of the ACs were recorded using a Genesis FTIRTM brand spectrophotometer ATI Mattson equipped with a DTGS (deuterated triglycine sulfate) detector in transmission mode between 4000 and 400 cm−1. The pH measurement was performed with a Mettler Toledo AG pH meter. For the determination of the crystalline or amorphous structure of the activated carbons, the diffraction patterns were recorded using an STOE Stadi-p device (STOE & Co. GmbH, Darmstadt, Germany) with a Cu Kα1 radiation (λ = 1.54056 Å, 40 kV, 30 mA) between 0 and 70°. With regard to energy-dispersive X-ray spectrometry (EDXS), a Euro Vector CHNS-O brand device (Euro EA 3000) made it possible to identify the different elements present in the activated carbons. The Brunauer–Emmett–Teller (BET) coupled with Barrett–Joyner–Halenda (BJH) method gave some textural properties of the ACs by N2 adsorption at 77.13 K using a micrometric sorptometer (Thermo Electron Corporation, Sorptomatic Advanced Data Processing). A UV Genesys 10S UV-Vis spectrophotometer of Thermo Scientific was used at the maximum wavelength of 276 nm to determine the residual concentration of BPA in solution.


### 2.1.2. Chemicals

All the chemicals used were of analytical grade; bisphenol A with a purity of 97% was purchased from Sigma-Aldrich (Germany) and used as such. BPA stock solutions were prepared by dissolving 0.1 g BPA in 1 L containing less than 10% methanol. Test solutions were made by diluting the stock solutions. 0.10 mol/L NaOH and 0.10 mol/L HCl were used to adjust the pH of the solutions. Distilled water was used throughout this study. The chemical structure and some properties of bisphenol A are presented in Table1.Table 1
Chemical structure and some properties of BPA.
Chemical structureMolecular formulaC15H16O2Molecular weight (g/mol)228.29λmax (nm)276SolubilitySoluble in alcohol and slightly soluble in CCl4


### 2.1.3. Adsorbents

Two activated carbons (ACs) prepared by chemical means using H3PO4 and H2SO4 were used in this work. They were labeled as follows: NSPT for that activated with H3PO4, and NSST for that activated with H2SO4. The preparation of the ACs was made using CCD, and the variables whose effects were simultaneously observed are the concentration of the oxidizing agent (varying between 1 and 3 mol/L), the calcination temperature (between 400 and 700°C), and the calcination time in the oven (between 30 and 120 minutes) as described by Atemkeng et al. [17].


### 2.1.4. Characterization of the Adsorbent

Knowledge of some physicochemical properties of an adsorbent material is useful and contributes to the understanding of phenomena, e.g., in the adsorption process. The ACs have undergone few characterizations such as X-ray powder diffraction (XRD), Fourier transform infrared (FTIR) spectrophotometry, energy-dispersive X-ray spectrometry (EDXS), pH at the zero-point charge, and the BET/BJH method.


## 2.1.1. Instrumentation

Determination of the functional groups present on the surfaces the materials was done using Fourier transform infrared (FTIR) spectrophotometry. The IR spectra of the ACs were recorded using a Genesis FTIRTM brand spectrophotometer ATI Mattson equipped with a DTGS (deuterated triglycine sulfate) detector in transmission mode between 4000 and 400 cm−1. The pH measurement was performed with a Mettler Toledo AG pH meter. For the determination of the crystalline or amorphous structure of the activated carbons, the diffraction patterns were recorded using an STOE Stadi-p device (STOE & Co. GmbH, Darmstadt, Germany) with a Cu Kα1 radiation (λ = 1.54056 Å, 40 kV, 30 mA) between 0 and 70°. With regard to energy-dispersive X-ray spectrometry (EDXS), a Euro Vector CHNS-O brand device (Euro EA 3000) made it possible to identify the different elements present in the activated carbons. The Brunauer–Emmett–Teller (BET) coupled with Barrett–Joyner–Halenda (BJH) method gave some textural properties of the ACs by N2 adsorption at 77.13 K using a micrometric sorptometer (Thermo Electron Corporation, Sorptomatic Advanced Data Processing). A UV Genesys 10S UV-Vis spectrophotometer of Thermo Scientific was used at the maximum wavelength of 276 nm to determine the residual concentration of BPA in solution.


## 2.1.2. Chemicals

All the chemicals used were of analytical grade; bisphenol A with a purity of 97% was purchased from Sigma-Aldrich (Germany) and used as such. BPA stock solutions were prepared by dissolving 0.1 g BPA in 1 L containing less than 10% methanol. Test solutions were made by diluting the stock solutions. 0.10 mol/L NaOH and 0.10 mol/L HCl were used to adjust the pH of the solutions. Distilled water was used throughout this study. The chemical structure and some properties of bisphenol A are presented in Table1.Table 1
Chemical structure and some properties of BPA.
Chemical structureMolecular formulaC15H16O2Molecular weight (g/mol)228.29λmax (nm)276SolubilitySoluble in alcohol and slightly soluble in CCl4


## 2.1.3. Adsorbents

Two activated carbons (ACs) prepared by chemical means using H3PO4 and H2SO4 were used in this work. They were labeled as follows: NSPT for that activated with H3PO4, and NSST for that activated with H2SO4. The preparation of the ACs was made using CCD, and the variables whose effects were simultaneously observed are the concentration of the oxidizing agent (varying between 1 and 3 mol/L), the calcination temperature (between 400 and 700°C), and the calcination time in the oven (between 30 and 120 minutes) as described by Atemkeng et al. [17].


## 2.1.4. Characterization of the Adsorbent

Knowledge of some physicochemical properties of an adsorbent material is useful and contributes to the understanding of phenomena, e.g., in the adsorption process. The ACs have undergone few characterizations such as X-ray powder diffraction (XRD), Fourier transform infrared (FTIR) spectrophotometry, energy-dispersive X-ray spectrometry (EDXS), pH at the zero-point charge, and the BET/BJH method.


## 3. Methods

### 3.1. Batch Procedure

Adsorption experiments were carried out by mechanical agitation at ambient temperature and constant agitation speed of 150 rpm. For each run, the desired volume of BPA solution of known initial concentration was put in contact with the desired mass of each of the ACs at desired pH solution. After stirring, the solution was filtered using Whatman N˚1 filter paper and the filtrates were analyzed in order to obtain the BPA residual concentration. Similar measurements were carried out by considering the optimal operating conditions: pH of the solution, AC doses, stirring duration, and BPA initial concentration. Absorbance values were determined at a maximum wavelength of 276 nm. The amountQe of BPA adsorbed was calculated using the following expression:(1)Qe=Ci−Ce∗Vm,where Ci and Ce in mg/L are, respectively, the initial concentration and the residual concentration of BPA at a precise time, V is the volume of the analyte solution, and m is the AC dosage.


#### 3.1.1. Central Composite Design as Optimization Approach

The experimental data from predesigned tests are analyzed by the response surface methodology (RSM) to create a model relating to the independent variables and the responses [13, 18]. The central composite design reduces the number of error attempts needed to evaluate the effect of interactions between the different factors that govern a process [26]. Herein, the variables whose effects were simultaneously observed are the BPA solution concentration x1 varying between 20 and 180 mol/L, the AC dose x2 varying between 25 and 100 mg, the stirring duration x3 varying between 6 and 120 min, and the pH of BPA solution x4 varying between 2 and 10. The total number of tests N performed for a hybrid composite design consists of a test number of the factorial plane 2k, a test number of the star plane 2k, where k is the number of variant factors, and a centre test number nc, which is used to determine experimental errors and to verify the reproducibility of the results. The number of test performed was determined from the following equation:(2)N=2k+2k+nc=24+24+4=28.These experiments were randomly carried out to minimize the effects of uncontrolled factors using the Statgraphics Plus version 5.0 software. The expected response for each experiment was used to develop an empirical model. During optimization, the relationship of main factors, response, and interactions are fitted to the classical second-order polynomial equation that also includes the linear terms as follows [27]:(3)Y=b0+∑i=1nbixi+∑i=1nbiixi2+∑i=1n∑j>1nbijxixj+ε,where Y is the predicted response, b0 is the constant coefficient, bi are the coefficients of the linear terms, bij are the coefficients of the interaction terms, bii are the coefficients of the quadratic terms, xi and xj are the coded values of factors that have been varied during the removal of BPA onto the activated carbons, and ε is the uncertainty between experimental and predicted values.The expected response was the amount of bisphenol A adsorbed obtained from equation (1). The response was labeled as Y1 and Y2, respectively, for the quantity of BPA adsorbed onto NSPT and NSST AC.


#### 3.1.2. Codification of the Levels of the Variable

Codification of the levels of the variable consists of transforming each studied real value into coordinates inside a scale with dimensionless values, which must be proportional to its localization in the experimental space. Codification is of concern because it enables the investigation of variables of different orders of magnitude without the greater influencing the evaluation of the lesser. The following equation can be applied to transform a real valueZi into a coded value xi according to a determinate experimental design:(4)xi=Zi−Zi0ΔZi,where ΔZi is the distance between the real value in the central point and the real value in the superior or inferior level of a variable, and Zi0 is the real value in the central point. The experimental variables considered in this work are shown in Table 2.Table 2
Experimental variables.
Variable levelLow valueCentre pointHigh valueCoded value−101Concentration of BPA (mg/L)20100180Activated carbon dosage (mg)2562.5100Stirring time (min)663120pH of BPA solution2610


### 3.2. Kinetic and Isotherm Studies

The controlling mechanisms for BPA adsorption onto the ACs were deduced. Data from the kinetic experiments were analyzed using four different nonlinear kinetic models. For the adsorption equilibrium, data were analyzed using four model isotherms. The kinetic and isotherm models studied are presented in Table3.Table 3
Nonlinear form of the kinetic and isotherm model used.
ModelsNonlinear formsReferencesAdsorption kineticsPseudo-first-orderQt=Qe1−e−k1t[28]Pseudo-second-orderqt=k2Qe2t/1+k2Qet[29]Elovichqt=lnαβ/β+lnt/β[30]Intraparticle diffusionqt=Kidt1/2+C[31]Adsorption isothermsLangmuirQe=QmKLCe/1+KLCe[32]FreundlichQe=KfCe1/nQe=Qmexp−K′ε2[32]Dubinin–Radushkevichε=RTln1+1/Ce[33]E=2KD1/2TemkinQe=QmRT/ΔQlnKTCe[34]


#### 3.2.1. Error Function Analysis

For this study, the nonlinear regression method was applied using Microsoft Excel Solver function for fitting the curve. The best fit for experimental data was determined from the coefficient of sum of square of errors (SSE), the average relative error (ARE), the residual root mean square error (RMSE), the nonlinear chi-squared test (χ2), and determination coefficient (R2). The mathematical expression for error function analysis used in this work is presented in Table 4.Table 4
Mathematical expression for error function analysis.
Error functionsAbbreviationsFormulaeReferencesSum of square of the errorsSSESSE=∑i=1nqe,cal,i−qe,exp,i2[35]Average relative errorAREARE=100/n∑i=1nqe,exp,i−qe,cal,i/qe,exp,i[36]Residual root mean square errorRMSERMSE=∑i=1nqe,cal,i−qe,pre,i2/n−2[37]Nonlinear chi-squared testχ2χ2=∑i=1nqe,cal,i−qe,exp,i2/qe,exp,i[38]Coefficient of determinationR2R2=∑qe,cal,i−qe,exp,i2/∑qe,cal,i−qe,exp,i2+qe,cal,i−qm,exp,i2[36]whereqe,calc and qe,meas are, respectively, the calculated adsorbate concentration at equilibrium (mg/g) and the measured adsorbate concentration at equilibrium (mg/g).


### 3.3. Influence of Temperature on the Adsorption Process

The study of the influence of temperature on the adsorption process was carried out with NSPT carbon having presented the best adsorption performance. The temperature of the reaction medium varied from 303 to 353 K during the adsorption of bisphenol A. Under the optimal adsorption conditions obtained above, 30 mL of pollutant solution was put in contact with the required carbon mass. The various residual concentrations of pollutants were determined by dosage with an ultraviolet-visible spectrophotometer at the appropriate wavelength. Some thermodynamic parameters such as the standard free enthalpyΔG0, the standard enthalpy ΔH0, and the standard entropy ΔS0 were determined using the following equations:(5)ΔG0=ΔH0−TΔS0=−RTlnKD,(6)lnKD=ΔS0R−ΔH0R1T,where KD=Qe/Ce is the distribution constant, which represents the ability of the adsorbent to retain adsorbate, R (8.3145 J/mol.K) is the ideal gas constant, and T (K) is the temperature.


## 3.1. Batch Procedure

Adsorption experiments were carried out by mechanical agitation at ambient temperature and constant agitation speed of 150 rpm. For each run, the desired volume of BPA solution of known initial concentration was put in contact with the desired mass of each of the ACs at desired pH solution. After stirring, the solution was filtered using Whatman N˚1 filter paper and the filtrates were analyzed in order to obtain the BPA residual concentration. Similar measurements were carried out by considering the optimal operating conditions: pH of the solution, AC doses, stirring duration, and BPA initial concentration. Absorbance values were determined at a maximum wavelength of 276 nm. The amountQe of BPA adsorbed was calculated using the following expression:(1)Qe=Ci−Ce∗Vm,where Ci and Ce in mg/L are, respectively, the initial concentration and the residual concentration of BPA at a precise time, V is the volume of the analyte solution, and m is the AC dosage.


### 3.1.1. Central Composite Design as Optimization Approach

The experimental data from predesigned tests are analyzed by the response surface methodology (RSM) to create a model relating to the independent variables and the responses [13, 18]. The central composite design reduces the number of error attempts needed to evaluate the effect of interactions between the different factors that govern a process [26]. Herein, the variables whose effects were simultaneously observed are the BPA solution concentration x1 varying between 20 and 180 mol/L, the AC dose x2 varying between 25 and 100 mg, the stirring duration x3 varying between 6 and 120 min, and the pH of BPA solution x4 varying between 2 and 10. The total number of tests N performed for a hybrid composite design consists of a test number of the factorial plane 2k, a test number of the star plane 2k, where k is the number of variant factors, and a centre test number nc, which is used to determine experimental errors and to verify the reproducibility of the results. The number of test performed was determined from the following equation:(2)N=2k+2k+nc=24+24+4=28.These experiments were randomly carried out to minimize the effects of uncontrolled factors using the Statgraphics Plus version 5.0 software. The expected response for each experiment was used to develop an empirical model. During optimization, the relationship of main factors, response, and interactions are fitted to the classical second-order polynomial equation that also includes the linear terms as follows [27]:(3)Y=b0+∑i=1nbixi+∑i=1nbiixi2+∑i=1n∑j>1nbijxixj+ε,where Y is the predicted response, b0 is the constant coefficient, bi are the coefficients of the linear terms, bij are the coefficients of the interaction terms, bii are the coefficients of the quadratic terms, xi and xj are the coded values of factors that have been varied during the removal of BPA onto the activated carbons, and ε is the uncertainty between experimental and predicted values.The expected response was the amount of bisphenol A adsorbed obtained from equation (1). The response was labeled as Y1 and Y2, respectively, for the quantity of BPA adsorbed onto NSPT and NSST AC.


### 3.1.2. Codification of the Levels of the Variable

Codification of the levels of the variable consists of transforming each studied real value into coordinates inside a scale with dimensionless values, which must be proportional to its localization in the experimental space. Codification is of concern because it enables the investigation of variables of different orders of magnitude without the greater influencing the evaluation of the lesser. The following equation can be applied to transform a real valueZi into a coded value xi according to a determinate experimental design:(4)xi=Zi−Zi0ΔZi,where ΔZi is the distance between the real value in the central point and the real value in the superior or inferior level of a variable, and Zi0 is the real value in the central point. The experimental variables considered in this work are shown in Table 2.Table 2
Experimental variables.
Variable levelLow valueCentre pointHigh valueCoded value−101Concentration of BPA (mg/L)20100180Activated carbon dosage (mg)2562.5100Stirring time (min)663120pH of BPA solution2610


## 3.1.1. Central Composite Design as Optimization Approach

The experimental data from predesigned tests are analyzed by the response surface methodology (RSM) to create a model relating to the independent variables and the responses [13, 18]. The central composite design reduces the number of error attempts needed to evaluate the effect of interactions between the different factors that govern a process [26]. Herein, the variables whose effects were simultaneously observed are the BPA solution concentration x1 varying between 20 and 180 mol/L, the AC dose x2 varying between 25 and 100 mg, the stirring duration x3 varying between 6 and 120 min, and the pH of BPA solution x4 varying between 2 and 10. The total number of tests N performed for a hybrid composite design consists of a test number of the factorial plane 2k, a test number of the star plane 2k, where k is the number of variant factors, and a centre test number nc, which is used to determine experimental errors and to verify the reproducibility of the results. The number of test performed was determined from the following equation:(2)N=2k+2k+nc=24+24+4=28.These experiments were randomly carried out to minimize the effects of uncontrolled factors using the Statgraphics Plus version 5.0 software. The expected response for each experiment was used to develop an empirical model. During optimization, the relationship of main factors, response, and interactions are fitted to the classical second-order polynomial equation that also includes the linear terms as follows [27]:(3)Y=b0+∑i=1nbixi+∑i=1nbiixi2+∑i=1n∑j>1nbijxixj+ε,where Y is the predicted response, b0 is the constant coefficient, bi are the coefficients of the linear terms, bij are the coefficients of the interaction terms, bii are the coefficients of the quadratic terms, xi and xj are the coded values of factors that have been varied during the removal of BPA onto the activated carbons, and ε is the uncertainty between experimental and predicted values.The expected response was the amount of bisphenol A adsorbed obtained from equation (1). The response was labeled as Y1 and Y2, respectively, for the quantity of BPA adsorbed onto NSPT and NSST AC.


## 3.1.2. Codification of the Levels of the Variable

Codification of the levels of the variable consists of transforming each studied real value into coordinates inside a scale with dimensionless values, which must be proportional to its localization in the experimental space. Codification is of concern because it enables the investigation of variables of different orders of magnitude without the greater influencing the evaluation of the lesser. The following equation can be applied to transform a real valueZi into a coded value xi according to a determinate experimental design:(4)xi=Zi−Zi0ΔZi,where ΔZi is the distance between the real value in the central point and the real value in the superior or inferior level of a variable, and Zi0 is the real value in the central point. The experimental variables considered in this work are shown in Table 2.Table 2
Experimental variables.
Variable levelLow valueCentre pointHigh valueCoded value−101Concentration of BPA (mg/L)20100180Activated carbon dosage (mg)2562.5100Stirring time (min)663120pH of BPA solution2610


## 3.2. Kinetic and Isotherm Studies

The controlling mechanisms for BPA adsorption onto the ACs were deduced. Data from the kinetic experiments were analyzed using four different nonlinear kinetic models. For the adsorption equilibrium, data were analyzed using four model isotherms. The kinetic and isotherm models studied are presented in Table3.Table 3
Nonlinear form of the kinetic and isotherm model used.
ModelsNonlinear formsReferencesAdsorption kineticsPseudo-first-orderQt=Qe1−e−k1t[28]Pseudo-second-orderqt=k2Qe2t/1+k2Qet[29]Elovichqt=lnαβ/β+lnt/β[30]Intraparticle diffusionqt=Kidt1/2+C[31]Adsorption isothermsLangmuirQe=QmKLCe/1+KLCe[32]FreundlichQe=KfCe1/nQe=Qmexp−K′ε2[32]Dubinin–Radushkevichε=RTln1+1/Ce[33]E=2KD1/2TemkinQe=QmRT/ΔQlnKTCe[34]


### 3.2.1. Error Function Analysis

For this study, the nonlinear regression method was applied using Microsoft Excel Solver function for fitting the curve. The best fit for experimental data was determined from the coefficient of sum of square of errors (SSE), the average relative error (ARE), the residual root mean square error (RMSE), the nonlinear chi-squared test (χ2), and determination coefficient (R2). The mathematical expression for error function analysis used in this work is presented in Table 4.Table 4
Mathematical expression for error function analysis.
Error functionsAbbreviationsFormulaeReferencesSum of square of the errorsSSESSE=∑i=1nqe,cal,i−qe,exp,i2[35]Average relative errorAREARE=100/n∑i=1nqe,exp,i−qe,cal,i/qe,exp,i[36]Residual root mean square errorRMSERMSE=∑i=1nqe,cal,i−qe,pre,i2/n−2[37]Nonlinear chi-squared testχ2χ2=∑i=1nqe,cal,i−qe,exp,i2/qe,exp,i[38]Coefficient of determinationR2R2=∑qe,cal,i−qe,exp,i2/∑qe,cal,i−qe,exp,i2+qe,cal,i−qm,exp,i2[36]whereqe,calc and qe,meas are, respectively, the calculated adsorbate concentration at equilibrium (mg/g) and the measured adsorbate concentration at equilibrium (mg/g).


## 3.2.1. Error Function Analysis

For this study, the nonlinear regression method was applied using Microsoft Excel Solver function for fitting the curve. The best fit for experimental data was determined from the coefficient of sum of square of errors (SSE), the average relative error (ARE), the residual root mean square error (RMSE), the nonlinear chi-squared test (χ2), and determination coefficient (R2). The mathematical expression for error function analysis used in this work is presented in Table 4.Table 4
Mathematical expression for error function analysis.
Error functionsAbbreviationsFormulaeReferencesSum of square of the errorsSSESSE=∑i=1nqe,cal,i−qe,exp,i2[35]Average relative errorAREARE=100/n∑i=1nqe,exp,i−qe,cal,i/qe,exp,i[36]Residual root mean square errorRMSERMSE=∑i=1nqe,cal,i−qe,pre,i2/n−2[37]Nonlinear chi-squared testχ2χ2=∑i=1nqe,cal,i−qe,exp,i2/qe,exp,i[38]Coefficient of determinationR2R2=∑qe,cal,i−qe,exp,i2/∑qe,cal,i−qe,exp,i2+qe,cal,i−qm,exp,i2[36]whereqe,calc and qe,meas are, respectively, the calculated adsorbate concentration at equilibrium (mg/g) and the measured adsorbate concentration at equilibrium (mg/g).


## 3.3. Influence of Temperature on the Adsorption Process

The study of the influence of temperature on the adsorption process was carried out with NSPT carbon having presented the best adsorption performance. The temperature of the reaction medium varied from 303 to 353 K during the adsorption of bisphenol A. Under the optimal adsorption conditions obtained above, 30 mL of pollutant solution was put in contact with the required carbon mass. The various residual concentrations of pollutants were determined by dosage with an ultraviolet-visible spectrophotometer at the appropriate wavelength. Some thermodynamic parameters such as the standard free enthalpyΔG0, the standard enthalpy ΔH0, and the standard entropy ΔS0 were determined using the following equations:(5)ΔG0=ΔH0−TΔS0=−RTlnKD,(6)lnKD=ΔS0R−ΔH0R1T,where KD=Qe/Ce is the distribution constant, which represents the ability of the adsorbent to retain adsorbate, R (8.3145 J/mol.K) is the ideal gas constant, and T (K) is the temperature.


## 4. Results and Discussion

### 4.1. Characterization of the Activated Carbons

#### 4.1.1. Fourier Transform Infrared (FTIR) Spectroscopy of NSPT and NSST Activated Carbons

In order to appreciate the effect of impregnation and calcination, the IR spectra of the NSST and NSPT carbons were superimposed on those of the raw agricultural residue NSB. From these spectra, we observe several absorption bands characteristic of the different functional groups present on their surfaces as shown in Figure1.Figure 1
Superimposed IR spectra of raw safou seeds, and NSPT and NSST activated carbons.The spectral shape of raw safou seeds reveals at 3649.44 cm−1 a vibration band attributable to the O-H bond of polysaccharides, lignin, alcohols, or phenols. At 2923.87 cm−1 and 2854.07 cm−1, bands characteristic of the vibrations of the C-H bonds of the aliphatic CH and CH2 are observed. At 1749.65 cm−1, a band characteristic of the vibrations of the C = C bonds of the aromatic rings appears. Around 1458.04 cm−1, we observe a band that can be attributed to vibrations of the O = C-O bonds of saturated and unsaturated fatty acids and to the C-O bond of the carbonyl groups of lignin [39]. Finally, around 923.91 cm−1, a band is observed that can be attributed to the vibrations of the C-H bond of the polynuclear aromatic rings [40].For the NSPT and NSST activated carbons, no band appears between 3649.44 and 1800 cm−1. This could be explained by a reorganization of the atoms within the molecular matrix of the precursor during the process of oxidation and carbonization. Furthermore, this could be due to the degradation of the material surface during the calcination process, thereby causing the departure of C, H, and O atoms in the form of gases (CO2, CH4, and CO) and the disappearance of the water molecules [37]. For NSPT carbon treated with phosphoric acid, the band appearing around 1397.25 cm−1 and absent on the spectral shape of the raw material and the NSST carbon can be attributed to the vibration of the C-O bond present in the P-O-P chains of polyphosphates [41]. This can be explained by the incorporation of phosphorus within the molecular matrix of the precursor during its impregnation with H3PO4. These observations are proof that activation and carbonization took place but without completely destroying the surface chemical functions of the precursor.


#### 4.1.2. X-Ray Diffraction (XRD) of NSPT and NSST Powder Activated Carbons

In order to have an idea on the amorphous or crystalline character of the activated carbons, the diffractograms of the samples were recorded and are represented in Figure2.Figure 2
Diffractograms obtained by XRD analysis of NSPT and NSST powder samples.These diffractograms show no peak, which may indicate a certain crystallinity of the samples. This observation confirms the amorphous nature of activated carbons, which is an advantageous property for adsorption processes and indicative of a porous structure [42]. Moreover, the observed amorphousness confirms that NSST and NSPT ACs could have high microporosity [43]. Therefore, safou seeds are good precursors for the preparation of powdered activated carbons.


#### 4.1.3. Energy-Dispersive X-Ray Spectrometry (EDXS) of NSPT and NSST Activated Carbons

Results of EDXS analysis of the activated carbons are depicted in Figure3, while Table 5 summarizes the elements present and their respective proportions.Figure 3
Screenshots of images obtained from EDXS analysis of activated carbon samples.Table 5
Elemental surface chemical composition of the activated carbon samples.
ACsElemental composition (%)COMgAlSiPSKCuZnClCaNSPT45.8340.69—0.980.869.380.380.050.900.67—0.12NSST52.0339.80.190.090.340.62.392.750.970.70——From the observation of the data in Table5, we notice that these samples contain mainly carbon. One can thus confirm that the materials are activated carbons because they obey the definition of an activated carbon, which is a highly carbonaceous amorphous material. We also observe that the main heteroatom present is oxygen. This can be explained by the presence of many oxygenated chemical functions on the surface of the NSPT and NSST carbons.


#### 4.1.4. BET/BJH Method and pH at the Zero-Point Charge

From Table6, it can be observed that the use of H2SO4 as an oxidizing agent leads to the development of micropores, while H3PO4 promotes the development of micropores and mesopores. Indeed, the use of phosphoric acid retains carbon and promotes the formation of crosslinks through dehydration, cyclization, and condensation reactions on the one hand. On the other hand, it allows the redistribution of biopolymers, which facilitates their conversion into aromatic rings and forms larger pores. This result can furthermore be used to explain the high removing capacity of NSPT adsorbent compared with that of NSST AC. The higher the porosity of adsorbent material, the higher the specific surface area and the number of sites available for the fixation of pollutants during the adsorption process [1]. The pH at the zero-point charge (pHzpc) was found to be 4.36 for NSPT and 6.94 for NSST as reported by Atemkeng et al. [17].Table 6
Specific surface area and some parameters of the porous structure of the ACs.
CAsSBET (m2/g)MicroporesMesoporesSBET (m2/g)Cv (cm3/g)Ad (nm)SBET (m2/g)Vc (cm3/g)Ad (nm)NSPT326.68443.950.140.3721.250.032.03NSST10.8913.850.011.13—Vc = cumulative volume; Ad = average diameter.


### 4.2. Batch Mode Adsorption Procedure

#### 4.2.1. The Analysis of BPA in Solution

The endocrine disruptor analysis in the solution was performed in the zero-order BPA absorption spectrum. As shown in Figure4, the solution used contained 20 mg/L BPA, whose zero-order absorption spectrum was read at wavelengths between 240 and 400 nm. As indicated, the analyte showed the maximum wavelength (λmax) at 276 nm while preparing the calibration curve of BPA at λmax.Figure 4
Absorption spectrum of BPA 20 mg/L recorded between 240 and 400 nm in aqueous solution and calibration curve.The UV-Vis spectrum in Figure1 recorded between 240 and 400 nm shows an absorption band at 276 nm. This band comes from the π−π∗ transitions characteristic of energy absorption by the C = C ethylenic chromophores of the aromatic ring of BPA [44]. Indeed, after absorption of energy by these chromophores, an excitation of one or more electrons from the π orbitals to the π∗ orbitals occurs.


#### 4.2.2. Experimental Design

The results obtained under the experimental conditions are shown in Table7. By referring to the values of the residues in this table, it can be seen that they are much lower than the experimental and theoretical values. Thus, there is a good correlation between these experimental and theoretical values because they are close. The quantities of bisphenol A adsorbed by the NSPT and NSST activated carbons are, respectively, represented by Y1 and Y2.Table 7
Experimental matrix of the adsorption model of bisphenol A onto NSPT and NSST ACs by the central composite design: experimental response and predicted response.
Testn°Concentration (mg/L)Mass (mg)Time (min)pHQuantity ads (mg/g) by NSPT,Y1Quantity ads (mg/g) by NSST,Y2Exp valPre valResidueExp valPre valResidue1−11116.7937.8301.0376.7477.4730.726200019.96312.4302.46710.07211.3081.2363−11−1−16.0654.5991.4666.1796.0360.1434−11−116.7247.2500.5265.3154.0921.2235−1−1−1118.72717.1331.5944.0795.7721.6936111115.33912.7662.57315.40710.2675.147000012.87212.6700.20213.85012.1118.2618−100010.58110.2230.35811.1599.2891.87911−119.2489.0240.2249.33910.1360.79710000−111.74410.4951.24911.85311.2120.64111001012.46813.4711.00312.61313.3650.75212010011.49712.1050.2547.33410.0742.74131−1−1113.26813.5220.25413.63112.1971.43414111−113.54314.4950.9529.4708.5080.96215−1−11116.89716.1770.7211.07910.3480.269161−11114.90415.7280.82415.44916.5231.0741700−1010.86811.0830.21511.05010.440.611811−1−17.7028.7631.0617.7937.4420.35119−1−1−1−18.62511.5342.9099.0810.1381.05820100010.50512.0801.57510.61313.0792.466211−1−1−111.99210.3131.67911.71911.9250.20622000013.81312.6701.14311.95912.1110.15223000013.80412.6701.13412.35912.1110.24824−1−11−112.98912.5700.41912.71612.8500.13425−111−17.0887.1700.0827.2017.5530.352260−10016.44317.0540.61115.89813.7532.14527000013.85012.6700.9112.06812.1110.043281−11−114.714.5100.1914.24514.3860.141Exp val = experimental value; pre val = predicted value; ads = adsorbed; residue =Exp val−pre val.


#### 4.2.3. Analysis of Variance

The quantities of bisphenol A adsorbed according to each of the variables investigated were evaluated by the following quadratic polynomial equations:(7)Y1=13.4158+0.312A−0.2624B+0.0207C+1.7784D−0.0002A2+0.0004AB+0.0001AC−0.0018AD+0.0013B2+0.0001BC−0.0049BD−0.0001C2−0.0021CD−0.0754D2.(8)Y2=10.3585+0.358A−0.0379B+0.0257C−0.0442D−0.0001A2−0.00003AB−0.00001AC+0.0036AD−0.0001B2+0.0001BC−0.004BD−0.0001C2+0.0020CD−0.0572D2,where the coded variables A, B, C, and D are, respectively, the concentration of the analyte, the activated carbon mass, the duration of agitation, and the pH of the analyte solution. The results of the quadratic regression model of the analysis of variance are presented in Table 8.Table 8
Analysis of variance of the regression model for BPA adsorption onto ACs.
Adsorption of bisphenol A onto NSPT ACAdsorption of bisphenol A onto NSST ACVariableDfY1VariableDfY2MSf-ratiop-valueMSf-ratiop-valueA115.51624.630.0508A164.642219.330.0007∗B1110.24132.880.0001∗B160.907718.220.0009∗C125.68537.660.016∗C139.729711.880.0043∗D116.8495.030.043∗D10.04120.010.913A215.94251.770.2060A212.54260.760.399AB128.98758.650.0115∗AB10.14490.040.838AC110.00142.980.107AC10.06260.020.893AD15.70251.700.214AD121.5186.440.024∗B219.40182.800.117B210.17810.050.821BC12.36080.700.416BC11.42860.430.524BD18.67892.590.131BD15.86481.750.208C210.39800.120.735C210.00570.000.967CD13.97001.180.296CD13.47541.040.326D213.75991.120.308D212.16580.650.435R2 = 85.13%; R2adjusted = 69.12%R2 = 83.71%; R2adjusted = 66.17%∗ = significant factor; Df = degree of freedom; MS = mean square.These tables indicate that equation models (7) and (8) can be used adequately to describe the adsorption of bisphenol A over a wide range of operating conditions. With regard to the adsorption of this pollutant, the quadratic interactions AB, AD, A2, and the effect of D, A, B, and C have a significant influence. They all have probability values ≤0.05. Thus, it can be concluded that this statistical model is adequate to predict the adsorption of bisphenol A by NSPT and NSST ACs in the range of variables studied. For the models presented by the quadratic equations, values of correlation coefficients R2 and R2adjusted of 85.23% and 69.12% for Y1 and of 83.71% and 66.17% for Y2 imply each that the regression model establishes an acceptable correlation between the experimental values. They can thus provide a better relational explanation between the independent variables and the response.


#### 4.2.4. Pareto Diagram and Analysis of Adsorption Response Surfaces of Bisphenol A

Figure5 presents Pareto diagrams and response surfaces that show the variation in the amount of BPA adsorbed. The response surfaces obtained offer the possibility of observing the behavior of the quantity of bisphenol A adsorbed. One can note the synergistic or antagonistic effect of all the variables investigated. From Pareto diagrams, it is possible to establish an interaction link between the variables taken into account during the adsorption process. It can be observed that the quadratic interactions AD, AB, and the effects of the concentration of the pollutant solution, the duration of agitation, and its pH have a significant and synergistic effect on the quantity of bisphenol A adsorbed.Figure 5
Pareto diagrams and response surface plots showing effective factors and their mutual effects on the removal capacity of bisphenol A onto ACs. (a) The effects of pH and concentration for NSST ACs and (b) the effects of adsorbent mass and solution concentration for NSPT.
(a)and (b)From Figure5(a), it is observed that the quantity of bisphenol A adsorbed by the NSST AC increases with the pH. Here, we are witnessing a protonation of the –OH functions of bisphenol A molecules. This phenomenon promotes the establishment of electrostatic interactions of an attractive nature between the bisphenol A molecules positively charged in solution and the NSST activated carbon whose surface is negatively charged. This is due to the fact that the pHZPC (6, 94) of NSST AC is greater than pH 6 of the analyte solution [45].Figure5(b) shows that the quantity of adsorbed pollutant decreases when the mass of carbon increases from 25 to 100 mg. This behavior may be due to the number of adsorption sites, which decreases with the amount of adsorbent [46]. Indeed, the increase in adsorbent dosage makes it possible to increase the adsorption sites, but the bisphenol A molecules have more difficulty approaching these sites because of the congestion. This crowding could be due to the formation of clusters of carbon particles when the mass increases. This thus induces a reduction in the total adsorption surface by preventing the diffusion of analyte molecules inside the pores of the carbon. Consequently, this results in a decrease in the amount of bisphenol A adsorbed per unit mass [47].Finally, the quantity of bisphenol A adsorbed by NSPT and NSST adsorbents increases when the concentration increases from 20 to 180 mg/L. This phenomenon can be explained by the fact that at low concentration, there are still adsorption sites not occupied by the molecules of BPA. Thus, the increase in the concentration of the pollutant in the medium will shift the equilibrium of the system in the direction that favors the accumulation of the analyte on the AC surface. Vacant adsorption sites will be more occupied with the increase in the concentration of the pollutant [48].


#### 4.2.5. Model Validation and Experimental Confirmations

After optimizing the adsorption of bisphenol A, the optimal preparation conditions were obtained. The predicted optimal conditions were 90 mg/L and 180 mg/L of pollutant solution for adsorption onto NSPT and NSST activated carbons, respectively. As for the activated carbon mass, it was 25 mg for each. The stirring time was 90 min for adsorption onto NSPT and 120 min onto NSST, while the pollutant solution pH was 6. Under these predicted optimal conditions, the amounts of bisphenol A adsorbed were 16.706 mg/g for adsorption onto NSST and 17.825 mg/g onto NSPT. The results obtained experimentally from these predicted conditions are compared with those predicted and summarized in Table9.Table 9
Summary of predicted data and experimental data for bisphenol A adsorption.
NSPT activated carbonVariableY1 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre val902590621.19817,8253.373NSST carbon activatedVariableY2 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre.val18025120618.88216.7062.176Exp val = experimental value; pre val = predicted value; residue =Exp val−pre val.In view of the low residual values in this table that are well below the predicted and experimental values, it can be concluded that the responses obtained experimentally are close to those predicted. Under optimal conditions, there is no significant difference between the predicted results and those resulting from experimentation.


### 4.3. Bisphenol A Adsorption Kinetic onto NSPT and NSST Activated Carbons

Figure6 presents the nonlinear regression model of the kinetic models and the variation of the quantity of pollutant adsorbed, respectively, by NSPT and NSST ACs according to the stirring time.Figure 6
Influence of contact time and nonlinear plot of kinetic models. (a) Adsorption onto NSPT (25 mg,V = 150 tr/min, pH = 6, CBPA = 90 mg/L) and (b) adsorption onto NSST (25 mg, V = 150 tr/min, pH = 6, CBPA = 180 mg/L).
(a)and (b)From Figure6, the amounts of bisphenol A adsorbed increase in three phases over time. During the first 15 minutes, the adsorption process is fast with both adsorbents. This observation may be due to the fact that the adsorption sites initially present on the adsorbent surfaces are still free and the mass transfer forces are important. Indeed, when bisphenol A molecules are put together with activated carbon, interactions are established between the adsorption sites and the adsorbate molecules. This could result in the formation of several types of bonds contributing to the retention of the pollutant on the surface of the activated carbons. The adsorption process is less rapid following the second phase, and equilibrium is reached after 80 minutes. The quantities of bisphenol A adsorbed onto NSPT and NSST activated carbons are, respectively, 19.182 and 14.901 mg/g. This observation can be explained by the fact that there is a reduction in the number of adsorption sites after the binding of the first BPA molecules. Thus, those still present in solution compete for the few sites still available. In addition, the pore openings could be clogged by the already adsorbed molecules, thus preventing others from diffusing inside the pores [49]. According to Alvarez et al. [50], the presence of carboxylic and hydroxyl functions on the AC surface would prevent the adsorption of organic compounds. During this phase, the quantities adsorbed can be attributed to π-π interactions. Qualified as dispersive interactions, they take place between the π electrons of the aromatic nucleus (deficient in electrons) of the pollutants and the π electrons of the surface of the activated carbon (rich in electrons) [5, 51]. Figure 7 illustrates the possible interactions that contribute to the bisphenol A adsorption mechanism adapted from Tran et al. [52].Figure 7
Possible interactions contributing to the adsorption mechanism of BPA on activated carbon (except for the filling of pores).The third phase is observed beyond 80 minutes with the gradual formation of a plateau reflecting the saturation of the adsorption sites.In order to decide on the kinetics that better describes the sorption process involved, the different error functions mentioned earlier have been determined. The different values of these errors are shown in Table10 for each of the kinetic models studied.Table 10
Kinetic constants of the models and error functions.
Adsorption of bisphenol A onto NSST activated carbonModelParameterValueR2χ2ARESSERMSELagergrenQe (mg/g)18.070.9840.2650.0324.6812.164K1 (min−1)0.246Pseudo-second-orderQe (mg/g)14.9010.9940.0960.0751.6881.299K2 (g/min.mg)0.025h0 (mg/g.min)9.238t1/2 (min)2.064Elovichα (mg/g.min)6884.3980.8870.3140.2745.2162.284β (g/mg)0.686Intraparticle diffusionKp (mg/g.min0.5)1.7950.9770.4760.8576.9622.639C (mg/g)0.998Adsorption of bisphenol A onto NSPT activated carbonModelParameterValueR2χ2ARESSERMSELagergrenQe (mg/g)14.1810.9660.5100.1056.3852.527K1 (min−1)0.201Pseudo-second-orderQe (mg/g)19.1820.9930.0970.0101.2221.105K2 (g/min.mg)0.024h0 (mg/g.min)5.380t1/2 (min)2.810Elovichα (mg/g.min)377.9430.9470.1450.1571.9521.397β (g/mg)0.671Intraparticle diffusionKp (mg/g.min0.5)1.6870.9740.3293.1824.4832.117C (mg/g)0.700With regard to this table, we note that the values of the coefficients of determinationR2 are close to unity. This would mean that each of the kinetic models is likely to describe the adsorption of bisphenol A molecules onto ACs. However, based on the low values of χ2, it appears that the pseudo-second-order kinetic model better describes the adsorption kinetics of bisphenol A on each of the activated carbons. This conclusion can be confirmed based on the other error functions investigated as mentioned in Table 10. These error functions (ARE, SSE, and RMSE) have the lowest values. In view of the values of the initial adsorption speeds h0 and the half-reaction times t1/2, the adsorption of bisphenol A takes place more quickly onto NSST AC than onto NSPT AC. It could then be concluded that the limiting step of the adsorption process may be a chemical reaction and not a mass transport [29]. Moreover, according to Tran et al. [52], the pseudo-second-order kinetic model describes π-π type interactions for an adsorbate/adsorbent system. With regard to the Elovich model, low values of the β desorption rate constants were observed. These values of the order of 0.6 imply a tendency that there is a multiplicity of sites available for adsorption [53]. As for the Lagergren model, the values of the error functions do not plead in its favor. This observation then indicates that the interactions thus brought into play during the sorption process are not of the van der Waals type. The value of the thickness C of the layers of bisphenol A attached to the surface of the activated carbons is different from zero. This indicates that intraparticle diffusion is also involved in adsorption kinetics as part of the diffusion process. It is still not the only process controlling the scope of the diffusion rate involved in all the steps [54].


### 4.4. Bisphenol A Adsorption Isotherms: Nonlinear Regression Model and Error Functions

The adsorption isotherms have been described by nonlinear regression models, and the plots are presented in Figure8 in the form of quantities adsorbed as a function of the residual concentration.Figure 8
Nonlinear regression models of BPA adsorption isotherm models in aqueous medium: (a) Adsorption of BPA onto NSST and (b) Adsorption of BPA onto NSPT.
(a)(b)Figure8 shows the variation of the quantity of adsorbed pollutant as a function of residual concentration. The different values of the related error functions are shown in Table 11.Table 11
Isotherm constants of the models and error functions.
Adsorption of bisphenol A onto NSST activated carbonModelParameterValueR2χ2ARESSERMSELangmuirQm (mg/g)12,6840.9680.2486.2992.0551.433KL (L/mg)0,632RL0,010Freundlich1/n0.0110.9310.7942.5154.5312.129KF (L/mg)0.138TemkinKT (L/mg)0.0480.8990.83611.4019.3723.061ΔQ (kJ/mol)426.88D-RQm (mg/g)18.5250.8400.9546.87910.1413.184E (kJ/mol)0.002KD-R (L/mg)0.014Adsorption of bisphenol A onto NSPT activated carbonModelParameterValueR2χ2ARESSERMSELangmuirQm (mg/g)26.2470.9980.0222.4270.5120.716KL (L/mg)0.019RL0.248Freundlich1/n0.4780.9800.2010.9224.7312.175KF (L/mg)3.420TemkinKT (L/mg)0.2140.9440.1314.6392.8651.693ΔQ (kJ/mol)117.48D-RQm (mg/g)35.4370.9330.7692.6229.063.615E (kJ/mol)8.452KD-R (L/mg)0.007With regard to Table10, the determination coefficient values of R2 are all close to unity. The low values of χ2 and calculated error functions show that the Langmuir isotherm best describes the adsorption process of bisphenol A onto the activated carbons. In this specific case, the value of the Langmuir separation factor RL is such that 0 < RL < 1. This observation reflects a favorable adsorption [1]. It can thus be assumed that the adsorption is of a homogeneous nature on the surface of the carbons and is probably done by the formation of a monolayer, while the binding sites have equivalent energies. For this phenomenon, the adsorption sites have the same chemical affinity [55] because there is no interaction between the adsorbed pollutant molecules.With regard to the energy values obtained from the Dubinin–Radushkevich model, the process dominating the adsorption of BPA onto NSPT is ion exchange because 8 <E < 16 kJ/mol. For the adsorption onto NSST, the dominant sorption process is physisorption [56]. The heats of adsorption obtained from the Temkin model for the adsorption of bisphenol A on the different ACs are all positive. These involve an exothermic adsorption process [57, 58]. The values of 1⁄n determined from the Freundlich model show that there is an affinity between bisphenol A and activated carbons. Indeed, the more this value tends towards unity, the greater the affinity. They therefore reflect a significant affinity between the various adsorbate/adsorbent systems studied. These values of 0 < 1/n < 1 are suggesting a favorable adsorption of bisphenol A and heterogeneity of the AC adsorption sites [59].


### 4.5. Effect of Temperature on the Adsorption Process of Bisphenol A

The study of the influence of temperature on the adsorption process was made by varying the temperature from 303 to 353 K. It was made using the NSPT AC obtained by impregnation with H3PO4 having shown the best adsorption performance. The variation in the amount of bisphenol A adsorbed as a function of the temperature of the reaction medium is shown in Figure 9.Figure 9
Influence of temperature on bisphenol A adsorption onto NSPT AC (90 mg/L, pH 6, 25 mg,V = 150 tr/min, t = 80 min). (a) Effect of temperature on the quantity adsorbed and (b) linear plot of the distribution coefficient as a function of temperature.
(a)and (b)With regard to Figure9(a), the increase in the temperature of the reaction medium leads to a reduction in the quantity of bisphenol A adsorbed. This phenomenon, in accordance with Svante August Arrhenius law, suggests that the surface reaction is exothermic, and therefore, each increase in temperature disadvantages its progress [60]. In another aspect, the increase in temperature further leads to the solubilization of bisphenol A molecules, thus leading to a decrease in the quantity adsorbed [61]. The plot of ln (KD) = f (1/T) presented in Figure 9(b) made it possible to deduce the entropy, the standard enthalpy, and the standard free enthalpy relating to the adsorbent/adsorbate system studied. The results are listed in Table 12.Table 12
Thermodynamic parameters of bisphenol A adsorption onto NSPT activated carbon.
PollutantT (K)ΔS° (J/mol.K)ΔH° (kJ/mol)ΔG° (kJ/mol)BPA303−140.856−40.1853.4893133.7033234.6233335.27035310.895From Table11, an examination of the negative value of ΔS° shows that the adsorption takes place in a more or less orderly way at the adsorbent/adsorbate interface. The value of the standard adsorption enthalpy of less than 40 kJ/mol shows that physisorption is the dominant process [1, 62]. Furthermore, the ΔH° values are negative, indicating that the adsorption of BPA onto NSPT activated carbon is exothermic. The same observation was reported by Wang and Zhang [5]. One can also note that the values of ΔG° > 0 increase with the increase in temperature. This could be explained by the fact that adsorption becomes difficult and unfavorable when the temperature increases: the adsorption process was nonspontaneous [63]. Moreover, this observation indicates an increase in disorder with energy redistribution during adsorption at the adsorbate/adsorbent interface [62]. However, it is essential to note that the entropy could increase more because it is influenced by that of the universe; during the adsorption process, the system was not isolated.


## 4.1. Characterization of the Activated Carbons

### 4.1.1. Fourier Transform Infrared (FTIR) Spectroscopy of NSPT and NSST Activated Carbons

In order to appreciate the effect of impregnation and calcination, the IR spectra of the NSST and NSPT carbons were superimposed on those of the raw agricultural residue NSB. From these spectra, we observe several absorption bands characteristic of the different functional groups present on their surfaces as shown in Figure1.Figure 1
Superimposed IR spectra of raw safou seeds, and NSPT and NSST activated carbons.The spectral shape of raw safou seeds reveals at 3649.44 cm−1 a vibration band attributable to the O-H bond of polysaccharides, lignin, alcohols, or phenols. At 2923.87 cm−1 and 2854.07 cm−1, bands characteristic of the vibrations of the C-H bonds of the aliphatic CH and CH2 are observed. At 1749.65 cm−1, a band characteristic of the vibrations of the C = C bonds of the aromatic rings appears. Around 1458.04 cm−1, we observe a band that can be attributed to vibrations of the O = C-O bonds of saturated and unsaturated fatty acids and to the C-O bond of the carbonyl groups of lignin [39]. Finally, around 923.91 cm−1, a band is observed that can be attributed to the vibrations of the C-H bond of the polynuclear aromatic rings [40].For the NSPT and NSST activated carbons, no band appears between 3649.44 and 1800 cm−1. This could be explained by a reorganization of the atoms within the molecular matrix of the precursor during the process of oxidation and carbonization. Furthermore, this could be due to the degradation of the material surface during the calcination process, thereby causing the departure of C, H, and O atoms in the form of gases (CO2, CH4, and CO) and the disappearance of the water molecules [37]. For NSPT carbon treated with phosphoric acid, the band appearing around 1397.25 cm−1 and absent on the spectral shape of the raw material and the NSST carbon can be attributed to the vibration of the C-O bond present in the P-O-P chains of polyphosphates [41]. This can be explained by the incorporation of phosphorus within the molecular matrix of the precursor during its impregnation with H3PO4. These observations are proof that activation and carbonization took place but without completely destroying the surface chemical functions of the precursor.


### 4.1.2. X-Ray Diffraction (XRD) of NSPT and NSST Powder Activated Carbons

In order to have an idea on the amorphous or crystalline character of the activated carbons, the diffractograms of the samples were recorded and are represented in Figure2.Figure 2
Diffractograms obtained by XRD analysis of NSPT and NSST powder samples.These diffractograms show no peak, which may indicate a certain crystallinity of the samples. This observation confirms the amorphous nature of activated carbons, which is an advantageous property for adsorption processes and indicative of a porous structure [42]. Moreover, the observed amorphousness confirms that NSST and NSPT ACs could have high microporosity [43]. Therefore, safou seeds are good precursors for the preparation of powdered activated carbons.


### 4.1.3. Energy-Dispersive X-Ray Spectrometry (EDXS) of NSPT and NSST Activated Carbons

Results of EDXS analysis of the activated carbons are depicted in Figure3, while Table 5 summarizes the elements present and their respective proportions.Figure 3
Screenshots of images obtained from EDXS analysis of activated carbon samples.Table 5
Elemental surface chemical composition of the activated carbon samples.
ACsElemental composition (%)COMgAlSiPSKCuZnClCaNSPT45.8340.69—0.980.869.380.380.050.900.67—0.12NSST52.0339.80.190.090.340.62.392.750.970.70——From the observation of the data in Table5, we notice that these samples contain mainly carbon. One can thus confirm that the materials are activated carbons because they obey the definition of an activated carbon, which is a highly carbonaceous amorphous material. We also observe that the main heteroatom present is oxygen. This can be explained by the presence of many oxygenated chemical functions on the surface of the NSPT and NSST carbons.


### 4.1.4. BET/BJH Method and pH at the Zero-Point Charge

From Table6, it can be observed that the use of H2SO4 as an oxidizing agent leads to the development of micropores, while H3PO4 promotes the development of micropores and mesopores. Indeed, the use of phosphoric acid retains carbon and promotes the formation of crosslinks through dehydration, cyclization, and condensation reactions on the one hand. On the other hand, it allows the redistribution of biopolymers, which facilitates their conversion into aromatic rings and forms larger pores. This result can furthermore be used to explain the high removing capacity of NSPT adsorbent compared with that of NSST AC. The higher the porosity of adsorbent material, the higher the specific surface area and the number of sites available for the fixation of pollutants during the adsorption process [1]. The pH at the zero-point charge (pHzpc) was found to be 4.36 for NSPT and 6.94 for NSST as reported by Atemkeng et al. [17].Table 6
Specific surface area and some parameters of the porous structure of the ACs.
CAsSBET (m2/g)MicroporesMesoporesSBET (m2/g)Cv (cm3/g)Ad (nm)SBET (m2/g)Vc (cm3/g)Ad (nm)NSPT326.68443.950.140.3721.250.032.03NSST10.8913.850.011.13—Vc = cumulative volume; Ad = average diameter.


## 4.1.1. Fourier Transform Infrared (FTIR) Spectroscopy of NSPT and NSST Activated Carbons

In order to appreciate the effect of impregnation and calcination, the IR spectra of the NSST and NSPT carbons were superimposed on those of the raw agricultural residue NSB. From these spectra, we observe several absorption bands characteristic of the different functional groups present on their surfaces as shown in Figure1.Figure 1
Superimposed IR spectra of raw safou seeds, and NSPT and NSST activated carbons.The spectral shape of raw safou seeds reveals at 3649.44 cm−1 a vibration band attributable to the O-H bond of polysaccharides, lignin, alcohols, or phenols. At 2923.87 cm−1 and 2854.07 cm−1, bands characteristic of the vibrations of the C-H bonds of the aliphatic CH and CH2 are observed. At 1749.65 cm−1, a band characteristic of the vibrations of the C = C bonds of the aromatic rings appears. Around 1458.04 cm−1, we observe a band that can be attributed to vibrations of the O = C-O bonds of saturated and unsaturated fatty acids and to the C-O bond of the carbonyl groups of lignin [39]. Finally, around 923.91 cm−1, a band is observed that can be attributed to the vibrations of the C-H bond of the polynuclear aromatic rings [40].For the NSPT and NSST activated carbons, no band appears between 3649.44 and 1800 cm−1. This could be explained by a reorganization of the atoms within the molecular matrix of the precursor during the process of oxidation and carbonization. Furthermore, this could be due to the degradation of the material surface during the calcination process, thereby causing the departure of C, H, and O atoms in the form of gases (CO2, CH4, and CO) and the disappearance of the water molecules [37]. For NSPT carbon treated with phosphoric acid, the band appearing around 1397.25 cm−1 and absent on the spectral shape of the raw material and the NSST carbon can be attributed to the vibration of the C-O bond present in the P-O-P chains of polyphosphates [41]. This can be explained by the incorporation of phosphorus within the molecular matrix of the precursor during its impregnation with H3PO4. These observations are proof that activation and carbonization took place but without completely destroying the surface chemical functions of the precursor.


## 4.1.2. X-Ray Diffraction (XRD) of NSPT and NSST Powder Activated Carbons

In order to have an idea on the amorphous or crystalline character of the activated carbons, the diffractograms of the samples were recorded and are represented in Figure2.Figure 2
Diffractograms obtained by XRD analysis of NSPT and NSST powder samples.These diffractograms show no peak, which may indicate a certain crystallinity of the samples. This observation confirms the amorphous nature of activated carbons, which is an advantageous property for adsorption processes and indicative of a porous structure [42]. Moreover, the observed amorphousness confirms that NSST and NSPT ACs could have high microporosity [43]. Therefore, safou seeds are good precursors for the preparation of powdered activated carbons.


## 4.1.3. Energy-Dispersive X-Ray Spectrometry (EDXS) of NSPT and NSST Activated Carbons

Results of EDXS analysis of the activated carbons are depicted in Figure3, while Table 5 summarizes the elements present and their respective proportions.Figure 3
Screenshots of images obtained from EDXS analysis of activated carbon samples.Table 5
Elemental surface chemical composition of the activated carbon samples.
ACsElemental composition (%)COMgAlSiPSKCuZnClCaNSPT45.8340.69—0.980.869.380.380.050.900.67—0.12NSST52.0339.80.190.090.340.62.392.750.970.70——From the observation of the data in Table5, we notice that these samples contain mainly carbon. One can thus confirm that the materials are activated carbons because they obey the definition of an activated carbon, which is a highly carbonaceous amorphous material. We also observe that the main heteroatom present is oxygen. This can be explained by the presence of many oxygenated chemical functions on the surface of the NSPT and NSST carbons.


## 4.1.4. BET/BJH Method and pH at the Zero-Point Charge

From Table6, it can be observed that the use of H2SO4 as an oxidizing agent leads to the development of micropores, while H3PO4 promotes the development of micropores and mesopores. Indeed, the use of phosphoric acid retains carbon and promotes the formation of crosslinks through dehydration, cyclization, and condensation reactions on the one hand. On the other hand, it allows the redistribution of biopolymers, which facilitates their conversion into aromatic rings and forms larger pores. This result can furthermore be used to explain the high removing capacity of NSPT adsorbent compared with that of NSST AC. The higher the porosity of adsorbent material, the higher the specific surface area and the number of sites available for the fixation of pollutants during the adsorption process [1]. The pH at the zero-point charge (pHzpc) was found to be 4.36 for NSPT and 6.94 for NSST as reported by Atemkeng et al. [17].Table 6
Specific surface area and some parameters of the porous structure of the ACs.
CAsSBET (m2/g)MicroporesMesoporesSBET (m2/g)Cv (cm3/g)Ad (nm)SBET (m2/g)Vc (cm3/g)Ad (nm)NSPT326.68443.950.140.3721.250.032.03NSST10.8913.850.011.13—Vc = cumulative volume; Ad = average diameter.


## 4.2. Batch Mode Adsorption Procedure

### 4.2.1. The Analysis of BPA in Solution

The endocrine disruptor analysis in the solution was performed in the zero-order BPA absorption spectrum. As shown in Figure4, the solution used contained 20 mg/L BPA, whose zero-order absorption spectrum was read at wavelengths between 240 and 400 nm. As indicated, the analyte showed the maximum wavelength (λmax) at 276 nm while preparing the calibration curve of BPA at λmax.Figure 4
Absorption spectrum of BPA 20 mg/L recorded between 240 and 400 nm in aqueous solution and calibration curve.The UV-Vis spectrum in Figure1 recorded between 240 and 400 nm shows an absorption band at 276 nm. This band comes from the π−π∗ transitions characteristic of energy absorption by the C = C ethylenic chromophores of the aromatic ring of BPA [44]. Indeed, after absorption of energy by these chromophores, an excitation of one or more electrons from the π orbitals to the π∗ orbitals occurs.


### 4.2.2. Experimental Design

The results obtained under the experimental conditions are shown in Table7. By referring to the values of the residues in this table, it can be seen that they are much lower than the experimental and theoretical values. Thus, there is a good correlation between these experimental and theoretical values because they are close. The quantities of bisphenol A adsorbed by the NSPT and NSST activated carbons are, respectively, represented by Y1 and Y2.Table 7
Experimental matrix of the adsorption model of bisphenol A onto NSPT and NSST ACs by the central composite design: experimental response and predicted response.
Testn°Concentration (mg/L)Mass (mg)Time (min)pHQuantity ads (mg/g) by NSPT,Y1Quantity ads (mg/g) by NSST,Y2Exp valPre valResidueExp valPre valResidue1−11116.7937.8301.0376.7477.4730.726200019.96312.4302.46710.07211.3081.2363−11−1−16.0654.5991.4666.1796.0360.1434−11−116.7247.2500.5265.3154.0921.2235−1−1−1118.72717.1331.5944.0795.7721.6936111115.33912.7662.57315.40710.2675.147000012.87212.6700.20213.85012.1118.2618−100010.58110.2230.35811.1599.2891.87911−119.2489.0240.2249.33910.1360.79710000−111.74410.4951.24911.85311.2120.64111001012.46813.4711.00312.61313.3650.75212010011.49712.1050.2547.33410.0742.74131−1−1113.26813.5220.25413.63112.1971.43414111−113.54314.4950.9529.4708.5080.96215−1−11116.89716.1770.7211.07910.3480.269161−11114.90415.7280.82415.44916.5231.0741700−1010.86811.0830.21511.05010.440.611811−1−17.7028.7631.0617.7937.4420.35119−1−1−1−18.62511.5342.9099.0810.1381.05820100010.50512.0801.57510.61313.0792.466211−1−1−111.99210.3131.67911.71911.9250.20622000013.81312.6701.14311.95912.1110.15223000013.80412.6701.13412.35912.1110.24824−1−11−112.98912.5700.41912.71612.8500.13425−111−17.0887.1700.0827.2017.5530.352260−10016.44317.0540.61115.89813.7532.14527000013.85012.6700.9112.06812.1110.043281−11−114.714.5100.1914.24514.3860.141Exp val = experimental value; pre val = predicted value; ads = adsorbed; residue =Exp val−pre val.


### 4.2.3. Analysis of Variance

The quantities of bisphenol A adsorbed according to each of the variables investigated were evaluated by the following quadratic polynomial equations:(7)Y1=13.4158+0.312A−0.2624B+0.0207C+1.7784D−0.0002A2+0.0004AB+0.0001AC−0.0018AD+0.0013B2+0.0001BC−0.0049BD−0.0001C2−0.0021CD−0.0754D2.(8)Y2=10.3585+0.358A−0.0379B+0.0257C−0.0442D−0.0001A2−0.00003AB−0.00001AC+0.0036AD−0.0001B2+0.0001BC−0.004BD−0.0001C2+0.0020CD−0.0572D2,where the coded variables A, B, C, and D are, respectively, the concentration of the analyte, the activated carbon mass, the duration of agitation, and the pH of the analyte solution. The results of the quadratic regression model of the analysis of variance are presented in Table 8.Table 8
Analysis of variance of the regression model for BPA adsorption onto ACs.
Adsorption of bisphenol A onto NSPT ACAdsorption of bisphenol A onto NSST ACVariableDfY1VariableDfY2MSf-ratiop-valueMSf-ratiop-valueA115.51624.630.0508A164.642219.330.0007∗B1110.24132.880.0001∗B160.907718.220.0009∗C125.68537.660.016∗C139.729711.880.0043∗D116.8495.030.043∗D10.04120.010.913A215.94251.770.2060A212.54260.760.399AB128.98758.650.0115∗AB10.14490.040.838AC110.00142.980.107AC10.06260.020.893AD15.70251.700.214AD121.5186.440.024∗B219.40182.800.117B210.17810.050.821BC12.36080.700.416BC11.42860.430.524BD18.67892.590.131BD15.86481.750.208C210.39800.120.735C210.00570.000.967CD13.97001.180.296CD13.47541.040.326D213.75991.120.308D212.16580.650.435R2 = 85.13%; R2adjusted = 69.12%R2 = 83.71%; R2adjusted = 66.17%∗ = significant factor; Df = degree of freedom; MS = mean square.These tables indicate that equation models (7) and (8) can be used adequately to describe the adsorption of bisphenol A over a wide range of operating conditions. With regard to the adsorption of this pollutant, the quadratic interactions AB, AD, A2, and the effect of D, A, B, and C have a significant influence. They all have probability values ≤0.05. Thus, it can be concluded that this statistical model is adequate to predict the adsorption of bisphenol A by NSPT and NSST ACs in the range of variables studied. For the models presented by the quadratic equations, values of correlation coefficients R2 and R2adjusted of 85.23% and 69.12% for Y1 and of 83.71% and 66.17% for Y2 imply each that the regression model establishes an acceptable correlation between the experimental values. They can thus provide a better relational explanation between the independent variables and the response.


### 4.2.4. Pareto Diagram and Analysis of Adsorption Response Surfaces of Bisphenol A

Figure5 presents Pareto diagrams and response surfaces that show the variation in the amount of BPA adsorbed. The response surfaces obtained offer the possibility of observing the behavior of the quantity of bisphenol A adsorbed. One can note the synergistic or antagonistic effect of all the variables investigated. From Pareto diagrams, it is possible to establish an interaction link between the variables taken into account during the adsorption process. It can be observed that the quadratic interactions AD, AB, and the effects of the concentration of the pollutant solution, the duration of agitation, and its pH have a significant and synergistic effect on the quantity of bisphenol A adsorbed.Figure 5
Pareto diagrams and response surface plots showing effective factors and their mutual effects on the removal capacity of bisphenol A onto ACs. (a) The effects of pH and concentration for NSST ACs and (b) the effects of adsorbent mass and solution concentration for NSPT.
(a)and (b)From Figure5(a), it is observed that the quantity of bisphenol A adsorbed by the NSST AC increases with the pH. Here, we are witnessing a protonation of the –OH functions of bisphenol A molecules. This phenomenon promotes the establishment of electrostatic interactions of an attractive nature between the bisphenol A molecules positively charged in solution and the NSST activated carbon whose surface is negatively charged. This is due to the fact that the pHZPC (6, 94) of NSST AC is greater than pH 6 of the analyte solution [45].Figure5(b) shows that the quantity of adsorbed pollutant decreases when the mass of carbon increases from 25 to 100 mg. This behavior may be due to the number of adsorption sites, which decreases with the amount of adsorbent [46]. Indeed, the increase in adsorbent dosage makes it possible to increase the adsorption sites, but the bisphenol A molecules have more difficulty approaching these sites because of the congestion. This crowding could be due to the formation of clusters of carbon particles when the mass increases. This thus induces a reduction in the total adsorption surface by preventing the diffusion of analyte molecules inside the pores of the carbon. Consequently, this results in a decrease in the amount of bisphenol A adsorbed per unit mass [47].Finally, the quantity of bisphenol A adsorbed by NSPT and NSST adsorbents increases when the concentration increases from 20 to 180 mg/L. This phenomenon can be explained by the fact that at low concentration, there are still adsorption sites not occupied by the molecules of BPA. Thus, the increase in the concentration of the pollutant in the medium will shift the equilibrium of the system in the direction that favors the accumulation of the analyte on the AC surface. Vacant adsorption sites will be more occupied with the increase in the concentration of the pollutant [48].


### 4.2.5. Model Validation and Experimental Confirmations

After optimizing the adsorption of bisphenol A, the optimal preparation conditions were obtained. The predicted optimal conditions were 90 mg/L and 180 mg/L of pollutant solution for adsorption onto NSPT and NSST activated carbons, respectively. As for the activated carbon mass, it was 25 mg for each. The stirring time was 90 min for adsorption onto NSPT and 120 min onto NSST, while the pollutant solution pH was 6. Under these predicted optimal conditions, the amounts of bisphenol A adsorbed were 16.706 mg/g for adsorption onto NSST and 17.825 mg/g onto NSPT. The results obtained experimentally from these predicted conditions are compared with those predicted and summarized in Table9.Table 9
Summary of predicted data and experimental data for bisphenol A adsorption.
NSPT activated carbonVariableY1 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre val902590621.19817,8253.373NSST carbon activatedVariableY2 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre.val18025120618.88216.7062.176Exp val = experimental value; pre val = predicted value; residue =Exp val−pre val.In view of the low residual values in this table that are well below the predicted and experimental values, it can be concluded that the responses obtained experimentally are close to those predicted. Under optimal conditions, there is no significant difference between the predicted results and those resulting from experimentation.


## 4.2.1. The Analysis of BPA in Solution

The endocrine disruptor analysis in the solution was performed in the zero-order BPA absorption spectrum. As shown in Figure4, the solution used contained 20 mg/L BPA, whose zero-order absorption spectrum was read at wavelengths between 240 and 400 nm. As indicated, the analyte showed the maximum wavelength (λmax) at 276 nm while preparing the calibration curve of BPA at λmax.Figure 4
Absorption spectrum of BPA 20 mg/L recorded between 240 and 400 nm in aqueous solution and calibration curve.The UV-Vis spectrum in Figure1 recorded between 240 and 400 nm shows an absorption band at 276 nm. This band comes from the π−π∗ transitions characteristic of energy absorption by the C = C ethylenic chromophores of the aromatic ring of BPA [44]. Indeed, after absorption of energy by these chromophores, an excitation of one or more electrons from the π orbitals to the π∗ orbitals occurs.


## 4.2.2. Experimental Design

The results obtained under the experimental conditions are shown in Table7. By referring to the values of the residues in this table, it can be seen that they are much lower than the experimental and theoretical values. Thus, there is a good correlation between these experimental and theoretical values because they are close. The quantities of bisphenol A adsorbed by the NSPT and NSST activated carbons are, respectively, represented by Y1 and Y2.Table 7
Experimental matrix of the adsorption model of bisphenol A onto NSPT and NSST ACs by the central composite design: experimental response and predicted response.
Testn°Concentration (mg/L)Mass (mg)Time (min)pHQuantity ads (mg/g) by NSPT,Y1Quantity ads (mg/g) by NSST,Y2Exp valPre valResidueExp valPre valResidue1−11116.7937.8301.0376.7477.4730.726200019.96312.4302.46710.07211.3081.2363−11−1−16.0654.5991.4666.1796.0360.1434−11−116.7247.2500.5265.3154.0921.2235−1−1−1118.72717.1331.5944.0795.7721.6936111115.33912.7662.57315.40710.2675.147000012.87212.6700.20213.85012.1118.2618−100010.58110.2230.35811.1599.2891.87911−119.2489.0240.2249.33910.1360.79710000−111.74410.4951.24911.85311.2120.64111001012.46813.4711.00312.61313.3650.75212010011.49712.1050.2547.33410.0742.74131−1−1113.26813.5220.25413.63112.1971.43414111−113.54314.4950.9529.4708.5080.96215−1−11116.89716.1770.7211.07910.3480.269161−11114.90415.7280.82415.44916.5231.0741700−1010.86811.0830.21511.05010.440.611811−1−17.7028.7631.0617.7937.4420.35119−1−1−1−18.62511.5342.9099.0810.1381.05820100010.50512.0801.57510.61313.0792.466211−1−1−111.99210.3131.67911.71911.9250.20622000013.81312.6701.14311.95912.1110.15223000013.80412.6701.13412.35912.1110.24824−1−11−112.98912.5700.41912.71612.8500.13425−111−17.0887.1700.0827.2017.5530.352260−10016.44317.0540.61115.89813.7532.14527000013.85012.6700.9112.06812.1110.043281−11−114.714.5100.1914.24514.3860.141Exp val = experimental value; pre val = predicted value; ads = adsorbed; residue =Exp val−pre val.


## 4.2.3. Analysis of Variance

The quantities of bisphenol A adsorbed according to each of the variables investigated were evaluated by the following quadratic polynomial equations:(7)Y1=13.4158+0.312A−0.2624B+0.0207C+1.7784D−0.0002A2+0.0004AB+0.0001AC−0.0018AD+0.0013B2+0.0001BC−0.0049BD−0.0001C2−0.0021CD−0.0754D2.(8)Y2=10.3585+0.358A−0.0379B+0.0257C−0.0442D−0.0001A2−0.00003AB−0.00001AC+0.0036AD−0.0001B2+0.0001BC−0.004BD−0.0001C2+0.0020CD−0.0572D2,where the coded variables A, B, C, and D are, respectively, the concentration of the analyte, the activated carbon mass, the duration of agitation, and the pH of the analyte solution. The results of the quadratic regression model of the analysis of variance are presented in Table 8.Table 8
Analysis of variance of the regression model for BPA adsorption onto ACs.
Adsorption of bisphenol A onto NSPT ACAdsorption of bisphenol A onto NSST ACVariableDfY1VariableDfY2MSf-ratiop-valueMSf-ratiop-valueA115.51624.630.0508A164.642219.330.0007∗B1110.24132.880.0001∗B160.907718.220.0009∗C125.68537.660.016∗C139.729711.880.0043∗D116.8495.030.043∗D10.04120.010.913A215.94251.770.2060A212.54260.760.399AB128.98758.650.0115∗AB10.14490.040.838AC110.00142.980.107AC10.06260.020.893AD15.70251.700.214AD121.5186.440.024∗B219.40182.800.117B210.17810.050.821BC12.36080.700.416BC11.42860.430.524BD18.67892.590.131BD15.86481.750.208C210.39800.120.735C210.00570.000.967CD13.97001.180.296CD13.47541.040.326D213.75991.120.308D212.16580.650.435R2 = 85.13%; R2adjusted = 69.12%R2 = 83.71%; R2adjusted = 66.17%∗ = significant factor; Df = degree of freedom; MS = mean square.These tables indicate that equation models (7) and (8) can be used adequately to describe the adsorption of bisphenol A over a wide range of operating conditions. With regard to the adsorption of this pollutant, the quadratic interactions AB, AD, A2, and the effect of D, A, B, and C have a significant influence. They all have probability values ≤0.05. Thus, it can be concluded that this statistical model is adequate to predict the adsorption of bisphenol A by NSPT and NSST ACs in the range of variables studied. For the models presented by the quadratic equations, values of correlation coefficients R2 and R2adjusted of 85.23% and 69.12% for Y1 and of 83.71% and 66.17% for Y2 imply each that the regression model establishes an acceptable correlation between the experimental values. They can thus provide a better relational explanation between the independent variables and the response.


## 4.2.4. Pareto Diagram and Analysis of Adsorption Response Surfaces of Bisphenol A

Figure5 presents Pareto diagrams and response surfaces that show the variation in the amount of BPA adsorbed. The response surfaces obtained offer the possibility of observing the behavior of the quantity of bisphenol A adsorbed. One can note the synergistic or antagonistic effect of all the variables investigated. From Pareto diagrams, it is possible to establish an interaction link between the variables taken into account during the adsorption process. It can be observed that the quadratic interactions AD, AB, and the effects of the concentration of the pollutant solution, the duration of agitation, and its pH have a significant and synergistic effect on the quantity of bisphenol A adsorbed.Figure 5
Pareto diagrams and response surface plots showing effective factors and their mutual effects on the removal capacity of bisphenol A onto ACs. (a) The effects of pH and concentration for NSST ACs and (b) the effects of adsorbent mass and solution concentration for NSPT.
(a)and (b)From Figure5(a), it is observed that the quantity of bisphenol A adsorbed by the NSST AC increases with the pH. Here, we are witnessing a protonation of the –OH functions of bisphenol A molecules. This phenomenon promotes the establishment of electrostatic interactions of an attractive nature between the bisphenol A molecules positively charged in solution and the NSST activated carbon whose surface is negatively charged. This is due to the fact that the pHZPC (6, 94) of NSST AC is greater than pH 6 of the analyte solution [45].Figure5(b) shows that the quantity of adsorbed pollutant decreases when the mass of carbon increases from 25 to 100 mg. This behavior may be due to the number of adsorption sites, which decreases with the amount of adsorbent [46]. Indeed, the increase in adsorbent dosage makes it possible to increase the adsorption sites, but the bisphenol A molecules have more difficulty approaching these sites because of the congestion. This crowding could be due to the formation of clusters of carbon particles when the mass increases. This thus induces a reduction in the total adsorption surface by preventing the diffusion of analyte molecules inside the pores of the carbon. Consequently, this results in a decrease in the amount of bisphenol A adsorbed per unit mass [47].Finally, the quantity of bisphenol A adsorbed by NSPT and NSST adsorbents increases when the concentration increases from 20 to 180 mg/L. This phenomenon can be explained by the fact that at low concentration, there are still adsorption sites not occupied by the molecules of BPA. Thus, the increase in the concentration of the pollutant in the medium will shift the equilibrium of the system in the direction that favors the accumulation of the analyte on the AC surface. Vacant adsorption sites will be more occupied with the increase in the concentration of the pollutant [48].


## 4.2.5. Model Validation and Experimental Confirmations

After optimizing the adsorption of bisphenol A, the optimal preparation conditions were obtained. The predicted optimal conditions were 90 mg/L and 180 mg/L of pollutant solution for adsorption onto NSPT and NSST activated carbons, respectively. As for the activated carbon mass, it was 25 mg for each. The stirring time was 90 min for adsorption onto NSPT and 120 min onto NSST, while the pollutant solution pH was 6. Under these predicted optimal conditions, the amounts of bisphenol A adsorbed were 16.706 mg/g for adsorption onto NSST and 17.825 mg/g onto NSPT. The results obtained experimentally from these predicted conditions are compared with those predicted and summarized in Table9.Table 9
Summary of predicted data and experimental data for bisphenol A adsorption.
NSPT activated carbonVariableY1 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre val902590621.19817,8253.373NSST carbon activatedVariableY2 (mg/g)ResidueA (mg/L)B (mg)C (min)DExp valPre.val18025120618.88216.7062.176Exp val = experimental value; pre val = predicted value; residue =Exp val−pre val.In view of the low residual values in this table that are well below the predicted and experimental values, it can be concluded that the responses obtained experimentally are close to those predicted. Under optimal conditions, there is no significant difference between the predicted results and those resulting from experimentation.


## 4.3. Bisphenol A Adsorption Kinetic onto NSPT and NSST Activated Carbons

Figure6 presents the nonlinear regression model of the kinetic models and the variation of the quantity of pollutant adsorbed, respectively, by NSPT and NSST ACs according to the stirring time.Figure 6
Influence of contact time and nonlinear plot of kinetic models. (a) Adsorption onto NSPT (25 mg,V = 150 tr/min, pH = 6, CBPA = 90 mg/L) and (b) adsorption onto NSST (25 mg, V = 150 tr/min, pH = 6, CBPA = 180 mg/L).
(a)and (b)From Figure6, the amounts of bisphenol A adsorbed increase in three phases over time. During the first 15 minutes, the adsorption process is fast with both adsorbents. This observation may be due to the fact that the adsorption sites initially present on the adsorbent surfaces are still free and the mass transfer forces are important. Indeed, when bisphenol A molecules are put together with activated carbon, interactions are established between the adsorption sites and the adsorbate molecules. This could result in the formation of several types of bonds contributing to the retention of the pollutant on the surface of the activated carbons. The adsorption process is less rapid following the second phase, and equilibrium is reached after 80 minutes. The quantities of bisphenol A adsorbed onto NSPT and NSST activated carbons are, respectively, 19.182 and 14.901 mg/g. This observation can be explained by the fact that there is a reduction in the number of adsorption sites after the binding of the first BPA molecules. Thus, those still present in solution compete for the few sites still available. In addition, the pore openings could be clogged by the already adsorbed molecules, thus preventing others from diffusing inside the pores [49]. According to Alvarez et al. [50], the presence of carboxylic and hydroxyl functions on the AC surface would prevent the adsorption of organic compounds. During this phase, the quantities adsorbed can be attributed to π-π interactions. Qualified as dispersive interactions, they take place between the π electrons of the aromatic nucleus (deficient in electrons) of the pollutants and the π electrons of the surface of the activated carbon (rich in electrons) [5, 51]. Figure 7 illustrates the possible interactions that contribute to the bisphenol A adsorption mechanism adapted from Tran et al. [52].Figure 7
Possible interactions contributing to the adsorption mechanism of BPA on activated carbon (except for the filling of pores).The third phase is observed beyond 80 minutes with the gradual formation of a plateau reflecting the saturation of the adsorption sites.In order to decide on the kinetics that better describes the sorption process involved, the different error functions mentioned earlier have been determined. The different values of these errors are shown in Table10 for each of the kinetic models studied.Table 10
Kinetic constants of the models and error functions.
Adsorption of bisphenol A onto NSST activated carbonModelParameterValueR2χ2ARESSERMSELagergrenQe (mg/g)18.070.9840.2650.0324.6812.164K1 (min−1)0.246Pseudo-second-orderQe (mg/g)14.9010.9940.0960.0751.6881.299K2 (g/min.mg)0.025h0 (mg/g.min)9.238t1/2 (min)2.064Elovichα (mg/g.min)6884.3980.8870.3140.2745.2162.284β (g/mg)0.686Intraparticle diffusionKp (mg/g.min0.5)1.7950.9770.4760.8576.9622.639C (mg/g)0.998Adsorption of bisphenol A onto NSPT activated carbonModelParameterValueR2χ2ARESSERMSELagergrenQe (mg/g)14.1810.9660.5100.1056.3852.527K1 (min−1)0.201Pseudo-second-orderQe (mg/g)19.1820.9930.0970.0101.2221.105K2 (g/min.mg)0.024h0 (mg/g.min)5.380t1/2 (min)2.810Elovichα (mg/g.min)377.9430.9470.1450.1571.9521.397β (g/mg)0.671Intraparticle diffusionKp (mg/g.min0.5)1.6870.9740.3293.1824.4832.117C (mg/g)0.700With regard to this table, we note that the values of the coefficients of determinationR2 are close to unity. This would mean that each of the kinetic models is likely to describe the adsorption of bisphenol A molecules onto ACs. However, based on the low values of χ2, it appears that the pseudo-second-order kinetic model better describes the adsorption kinetics of bisphenol A on each of the activated carbons. This conclusion can be confirmed based on the other error functions investigated as mentioned in Table 10. These error functions (ARE, SSE, and RMSE) have the lowest values. In view of the values of the initial adsorption speeds h0 and the half-reaction times t1/2, the adsorption of bisphenol A takes place more quickly onto NSST AC than onto NSPT AC. It could then be concluded that the limiting step of the adsorption process may be a chemical reaction and not a mass transport [29]. Moreover, according to Tran et al. [52], the pseudo-second-order kinetic model describes π-π type interactions for an adsorbate/adsorbent system. With regard to the Elovich model, low values of the β desorption rate constants were observed. These values of the order of 0.6 imply a tendency that there is a multiplicity of sites available for adsorption [53]. As for the Lagergren model, the values of the error functions do not plead in its favor. This observation then indicates that the interactions thus brought into play during the sorption process are not of the van der Waals type. The value of the thickness C of the layers of bisphenol A attached to the surface of the activated carbons is different from zero. This indicates that intraparticle diffusion is also involved in adsorption kinetics as part of the diffusion process. It is still not the only process controlling the scope of the diffusion rate involved in all the steps [54].


## 4.4. Bisphenol A Adsorption Isotherms: Nonlinear Regression Model and Error Functions

The adsorption isotherms have been described by nonlinear regression models, and the plots are presented in Figure8 in the form of quantities adsorbed as a function of the residual concentration.Figure 8
Nonlinear regression models of BPA adsorption isotherm models in aqueous medium: (a) Adsorption of BPA onto NSST and (b) Adsorption of BPA onto NSPT.
(a)(b)Figure8 shows the variation of the quantity of adsorbed pollutant as a function of residual concentration. The different values of the related error functions are shown in Table 11.Table 11
Isotherm constants of the models and error functions.
Adsorption of bisphenol A onto NSST activated carbonModelParameterValueR2χ2ARESSERMSELangmuirQm (mg/g)12,6840.9680.2486.2992.0551.433KL (L/mg)0,632RL0,010Freundlich1/n0.0110.9310.7942.5154.5312.129KF (L/mg)0.138TemkinKT (L/mg)0.0480.8990.83611.4019.3723.061ΔQ (kJ/mol)426.88D-RQm (mg/g)18.5250.8400.9546.87910.1413.184E (kJ/mol)0.002KD-R (L/mg)0.014Adsorption of bisphenol A onto NSPT activated carbonModelParameterValueR2χ2ARESSERMSELangmuirQm (mg/g)26.2470.9980.0222.4270.5120.716KL (L/mg)0.019RL0.248Freundlich1/n0.4780.9800.2010.9224.7312.175KF (L/mg)3.420TemkinKT (L/mg)0.2140.9440.1314.6392.8651.693ΔQ (kJ/mol)117.48D-RQm (mg/g)35.4370.9330.7692.6229.063.615E (kJ/mol)8.452KD-R (L/mg)0.007With regard to Table10, the determination coefficient values of R2 are all close to unity. The low values of χ2 and calculated error functions show that the Langmuir isotherm best describes the adsorption process of bisphenol A onto the activated carbons. In this specific case, the value of the Langmuir separation factor RL is such that 0 < RL < 1. This observation reflects a favorable adsorption [1]. It can thus be assumed that the adsorption is of a homogeneous nature on the surface of the carbons and is probably done by the formation of a monolayer, while the binding sites have equivalent energies. For this phenomenon, the adsorption sites have the same chemical affinity [55] because there is no interaction between the adsorbed pollutant molecules.With regard to the energy values obtained from the Dubinin–Radushkevich model, the process dominating the adsorption of BPA onto NSPT is ion exchange because 8 <E < 16 kJ/mol. For the adsorption onto NSST, the dominant sorption process is physisorption [56]. The heats of adsorption obtained from the Temkin model for the adsorption of bisphenol A on the different ACs are all positive. These involve an exothermic adsorption process [57, 58]. The values of 1⁄n determined from the Freundlich model show that there is an affinity between bisphenol A and activated carbons. Indeed, the more this value tends towards unity, the greater the affinity. They therefore reflect a significant affinity between the various adsorbate/adsorbent systems studied. These values of 0 < 1/n < 1 are suggesting a favorable adsorption of bisphenol A and heterogeneity of the AC adsorption sites [59].


## 4.5. Effect of Temperature on the Adsorption Process of Bisphenol A

The study of the influence of temperature on the adsorption process was made by varying the temperature from 303 to 353 K. It was made using the NSPT AC obtained by impregnation with H3PO4 having shown the best adsorption performance. The variation in the amount of bisphenol A adsorbed as a function of the temperature of the reaction medium is shown in Figure 9.Figure 9
Influence of temperature on bisphenol A adsorption onto NSPT AC (90 mg/L, pH 6, 25 mg,V = 150 tr/min, t = 80 min). (a) Effect of temperature on the quantity adsorbed and (b) linear plot of the distribution coefficient as a function of temperature.
(a)and (b)With regard to Figure9(a), the increase in the temperature of the reaction medium leads to a reduction in the quantity of bisphenol A adsorbed. This phenomenon, in accordance with Svante August Arrhenius law, suggests that the surface reaction is exothermic, and therefore, each increase in temperature disadvantages its progress [60]. In another aspect, the increase in temperature further leads to the solubilization of bisphenol A molecules, thus leading to a decrease in the quantity adsorbed [61]. The plot of ln (KD) = f (1/T) presented in Figure 9(b) made it possible to deduce the entropy, the standard enthalpy, and the standard free enthalpy relating to the adsorbent/adsorbate system studied. The results are listed in Table 12.Table 12
Thermodynamic parameters of bisphenol A adsorption onto NSPT activated carbon.
PollutantT (K)ΔS° (J/mol.K)ΔH° (kJ/mol)ΔG° (kJ/mol)BPA303−140.856−40.1853.4893133.7033234.6233335.27035310.895From Table11, an examination of the negative value of ΔS° shows that the adsorption takes place in a more or less orderly way at the adsorbent/adsorbate interface. The value of the standard adsorption enthalpy of less than 40 kJ/mol shows that physisorption is the dominant process [1, 62]. Furthermore, the ΔH° values are negative, indicating that the adsorption of BPA onto NSPT activated carbon is exothermic. The same observation was reported by Wang and Zhang [5]. One can also note that the values of ΔG° > 0 increase with the increase in temperature. This could be explained by the fact that adsorption becomes difficult and unfavorable when the temperature increases: the adsorption process was nonspontaneous [63]. Moreover, this observation indicates an increase in disorder with energy redistribution during adsorption at the adsorbate/adsorbent interface [62]. However, it is essential to note that the entropy could increase more because it is influenced by that of the universe; during the adsorption process, the system was not isolated.


## 5. Conclusion

The response surface methodology was successfully used to investigate the effects of initial concentration, adsorbent dosage, stirring time, and solution pH on the adsorption capacity of bisphenol A onto two activated carbons. The optimal conditions were obtained for 180 mg/L and 90 mg/L BPA solution, respectively, for NSST and NSPT materials. 25 mg adsorbent dosage and the solution pH 6 were all the same for adsorption on both adsorbents. Stirring time was, respectively, 90 min and 120 min for adsorption onto NSPT and NSST. These optimal conditions resulted in 18.882 and 21.198 mg/g removal for NSST and NSPT, respectively. From the analysis of the response surfaces, all the investigated parameters were found to have significant effects on BPA adsorption. The effect of the analyte concentration, the pH, and the adsorbent dosage were found to have the greatest synergistic impact on the quantity adsorbed. The sorption process was rapid and followed a pseudo-second-order model. The experimental data agreed well with Langmuir’s isotherm model. This study has shown that the response surface methodology based on a central composite design is a suitable method to optimize the operating conditions and maximize removal efficiency of bisphenol A in aqueous solution. A thermodynamic study revealed an exothermic and physisorption process.



---

*Source: 7717148-2022-10-07.xml*