A major environmental problem facing agro-food industries is their inability to efficiently manage the residues and by-products generated during the primary processing of raw materials. On the industrial scale, the processing of agricultural commodities generates huge quantities of lignocellulosic residues that are disposed of by land filling or incineration. These practices have a negative effect on man and his environment. For example, open-air incineration generates substances which are toxic to humans and constitute greenhouse gases that have damaging effects on the climate . It is because they are wrongly considered as wastes that these agricultural residues are treated in this way.
A large quantity of these lignocellulose residues can be recycled by fermenting them to generate methane which is a fuel, and ethanol which is both a solvent and a fuel . They can also be composted to obtain an enriched soil conditioner, which can improve soil fertility. Furthermore, a controlled pyrolysis of these lignocellulosic wastes produces activated carbons . Of these alternative treatments of the lignocellulosic residues, the transformation into activated carbons is more flexible; with chemical activation significantly reducing the energy cost of the process.
Activated carbons are carbon-based materials, found in granulated or black powdered form, and having a porous structure with very large surface areas. Amongst other applications, they are used in air filtration systems, purifying contaminated water, deodourization, dechlorination, and in gold mining operations . They are also used as catalysts or catalyst support . A lot of research has been carried out to synthesize activated carbons from lignocellulosic precursors. These precursors include corn shells , sunflower seed shells , Flamboyant (Delonix regia) pods , cotton cakes , rice hulls , coconut shells , garcinia cola pods  and olive stones .
It is in the quest for sustainable development that the recycling of agricultural and industrial waste, especially shea residues (Vitellaria paradoxa), arouses particular interest. Socioeconomically, shea is valuable in that its butter has medicinal, food and cosmetic values , and it is relatively abundant in many African countries. In the Republic of Chad, for example, shea industries produce 556,000 tonnes of almonds per year . While the nuts are used for butter production, nut shells, the more woody part of the fruit could be used as a raw material for the production of activated carbons. This alternative use of the nut shells has the advantage that, on the economic domain, value has been added onto the wastes; on the perspective of environmental pollution, the generation of greenhouse gases is avoided or highly minimised; and finally concerning sanitation the generation of pollutants harmful to man is avoided.
The objective of this research was to optimize the preparation of activated carbons using shea hulls as precursor. The process optimization was carried out using Response Surface Methodology (RSM), applying the Central Composite Design (CCD) in order to obtain maximum information with minimal experiments .
2. Materials and Methods
2.1. Preparation of Activated Carbon
The shea nut shells (CNK) (Figure 1) collected in Baktchoro in the West Tandjilé Division in Chad were washed with tap water to remove the remaining pulp and other forms of attached impurities. They were dehydrated by drying in an oven set at 105˚C for 24 hours.
Finally, they were crushed and sieved to obtain particles of size 0.5 - 2 mm that were preserved for use in this study.
Preparation of activated carbons
20 g of CNK was impregnated by mixing it with 40 mL solutions of phosphoric acid (H3PO4) or sulphuric acid (H2SO4) of concentrations ranging from 1 to 5 mol·L−1. After stirring for 24 hours to ensure reaction between the activating agent and the CNK, the mixture was then filtered and the filtrate dried in an oven at 105˚C for 1 hour to obtain the impregnated material. In each case, the impregnated materials was pyrolysed in a furnace of mark HERAEUS that was set at temperatures ranging from 400˚C to 700˚C for residence times of 30 to 120 minutes, and a constant heating rate of 6˚C·min−1. The residue from filtering the mixture was thoroughly washed with distilled water to remove the residual acid and then dried in a Binder oven set at 105˚C to obtain the activated carbon.
2.2. Determination of Yield and Iodine Number (IN) of Activated Carbons
The yield of activated carbon was determined as the mass of the dried activated carbon relative to the mass of the precursor, and given by the Equation (1):
Figure 1. Photograph of shea (Vitellaria paradoxa): (a) tree and fruit, (b) seeds, (c) nuts and (d) walnut shells.
Iodine Number (IN)
The iodine number of the activated carbons was determined using the procedure established by the American Society for Testing and Materials (ASTM) . This was done by mixing 100 mg of activated carbon with 30 mL of iodine solution of concentration 0.01 mol·L−1 in a sealed bottle. The mixture was stirred away from the sun rays at a speed of 200 rpm for 3 hours. The resulting mixture was quickly filtered using filter paper, and 10 mL of filtrate was titrated with a solution of sodium thiosulfate pentahydrate (Na2S2O3·5H2O), of concentration 0.005 N, using starch solution as indicator. From the chemical reaction given in Equation (2), Equation (3) was deduced and used to determine the amount of iodine adsorbed, thereby giving the iodine number (IN).
(mg·g−1) is the mass of iodine adsorbed per gram of adsorbent; (g) the mass of activated carbon; (mol·L−1) the initial concentration of I2 solution; (mol·L−1) the concentration of sodium thiosulfate; (mL) the volume of sodium thiosulfate; (mL), the volume of iodine solution that was mixed with the activated carbon; (mL) the volume of iodine solution titrated and (g·mol−1) the molar mass of I2.
2.3. Optimizing the Preparation of Activated Carbons
Experimental design: response surface methodology
The preparation of ACs was optimized through Response Surface Methodology. Central Composite Design was chosen to generate the experimental matrix in order to reduce the number of experiments to be performed while simultaneously studying the effects of several factors  . Thus, the total number of experiments performed was calculated from Equation (4) given below .
N is the total number of experiments, 2 indicates the number of levels per factor, k the number of factors studied, 2k the number of factor points, 2k the number of axial points and kc the number of points in the center that determined errors or reproducibility of the data.
Based on literature  and on previously conducted tests in our laboratory , the factors that significantly influence the physical properties of activated carbons are the concentration of activating agents (X1), the carbonization temperature (X2) and residence time (X3). In this work, H3PO4 and H2SO4 were used as activating agents. Table 1 presents the independent experimental variables used as calculated from Equation (5)  . The values of the experimental variables were coded as follows: −1 for the small value, 0 the value at the centre and +1 the high value.
Table 1. Independent experimental variables.
is the value of an independent variable, the actual value of a variable i, the value of the variable at the center point and is the step change.
To assess quantitatively and qualitatively the activated carbons obtained, the mass yield and iodine number were used as responses or performance indicators. The second order polynomial Equation (6) was used to evaluate the interactions between the factors  .
Y is the predicted response, and are the coded values of the factors, and are the regression coefficients for intercept, linear, quadratic and interaction terms respectively and the uncertainty between measured and predicted values.
Validating the statistical model
The empirical model obtained was validated using Statgraphic plus 5.0 software as the evaluation instrument. The model was regarded as good when the experimental values obtained were close to those predicted by the t-test and p-value. When the t-test value was high and the p-value small, then, the correlation coefficient was high. In general, a model was valid when R2 was greater than 75% and the p-value was 5% or less (with 95% confidence level) .
2.4. Characterization of Activated Carbons Prepared at the Optimal Conditions
The activated carbons prepared at the optimal conditions were characterized by various physicochemical methods. These optimal conditions were those provided by the Statgraphic plus 5.0 software; and they were the concentration of activating agents, the carbonization temperature and residence time, that led to high values of iodine number.
2.4.1. EDX Analysis
Elementary analysis using energy dispersion X-ray spectrometry was applied to investigate the presence and percentage of atoms that made up the activated carbons.
2.4.2. Bulk Density
To measure the bulk density, a mass of activated carbon was carefully poured into a graduated and previously weighed 25 cm3 volumetric flask. The difference in mass mAC compared to the occupied volume (Voc) was used to calculate the apparent density as follows:
2.4.3. Moisture Content
To determine the moisture level, an empty and well dried porcelain crucible was weighed (m0) using a precision scale of 1/1000b and 1 g of activated carbon was properly distributed on it and weighed again (m1). The combination consisting of crucible and activated carbon was then placed in the oven set at 100˚C for 1 hour. After cooling to room temperature, the system is again weighed (m2). The humidity, T was calculated by the relationship (8) :
2.4.4. pH and pH at the Point of Zero Charge (pHPZC)
The pH of the ACs was obtained by putting 30 mL of distilled water and 100 mg of AC in contact in a wash-bottle. The mixture was agitated for 24 hours; decanted, and the pH of the solution was measured using an ATC-mark pH-meter.
The pH at the point of zero charge (pHpzc) was determined according to the protocol described by Tchakala et al., (2012) . In washed bottles containing 30 mL of sodium chloride (NaCl) solution of concentration 0.1 mol·L−1, the pH was adjusted from 2 to 10 using 0.1 mol·L−1 solution of hydrochloric acid HCl or 0.1 mol·L−1 sodium hydroxide NaOH to obtain the initial values pHi.
100 mg of activated carbon was introduced into each bottle and the mixture stirred for 24 hours. After leaving the mixtures to settle, the final pH values, pHf were measured. The pHpzc was obtained as the intersection of the curve: and the bisector of the first quadrant.
2.4.5. Specific Surface Area
The porous structure of the ACs was determined by nitrogen adsorption at 77.13 K. This measurement was made using a sorptophotometric porosimeter (Thermo Electron Corporation Advanced Data Processing). The B.E.T method (Brunauer, Emmett and Teller) was used to determine the specific surface of the samples. The total pore volume VT is the volume of liquid nitrogen adsorbed at the relative pressure (P/P0) of 0.005 to 0.4 . Thus, from the B.E.T method coupled with the B.J.H. (Barret-Joyner-Halenda) the distribution of pore size, and consequently the microporous volume, Vmicro, mesoporous volume, Vmeso and the average pore diameter, Dp were calculated .
2.4.6. Quantification of Functional Surface Groups
The quantification of surface functional groups was carried out using the Boehm method (1966) . To do this, 0.1 g of activated carbon in 30 mL of a solution of NaHCO3 (0.1 mol·L−1), or Na2CO3 (0.1 mol·L−1) or NaOH (0.1 mol·L−1) was agitated at room temperature for 48 hours to determine the acid groups, or with 30 mL of HCl (0.1 mol·L−1) to determine the basic groups. The filtrates were titrated with HCl (0.1 mol·L−1) or NaOH (0.1 mol·L−1), depending on whether acid or basics groups were concerned; and in the presence of an indicator. The calculation was based on the principle that NaHCO3 neutralizes carboxylic acid groups, Na2CO3 neutralizes carboxylic and lactone groups and NaOH neutralizes carboxylic, lactone and phenolic groups and the basic groups are neutralized in their entirety by HCl. The number of equivalents (meq·g−1) and the type of functional groups were determined from the relationship (9) :
where, is the number of gram equivalents that reacted, while and are the number of gram equivalents before and after reaction.
2.4.7. Fourier-Transform Infrared Spectroscopy (IR-TF)
A Thermoscientific Nicolet iS5 brand spectrophotometer was used to investigate the chemical functional groups present on the surface of the activated carbons. The IR-TF spectra were recorded between the wavenumbers 3500 and 580 cm−1.
2.4.8. Scanning Electron Microscopy (SEM)
In order to study surface morphology and verify the porosity of the adsorbents, samples of the prepared activated carbons were examined by the PHILIPS XL30 S-FEG scanning electron microscope (SEM) equipped with a spectrometer that uses a cathode luminescence detector. Samples of each activated carbon prepared under the conditions of the best iodine number were selected for the SEM analysis.
3. Results and Discussion
3.1. Experimental Design Using the Central Composite Design
The central composite design was used with the aim of finding the optimal conditions for the preparation of activated carbon in a short time. Mass yield and iodine number were chosen as responses of the factors influencing (X1, X2 and X3) the preparation of the activated carbons. The results of the experimental matrix are presented in Table 2. Repetitive experiments 7, 8 and 11 were used to verify the reproducibility of the data during handling.
3.1.1. Iodine Number
The iodine number of carbon CAK-P (phosphoric acid activated carbon) and CAK-S (sulphuric acid activated carbon) is an important parameter which provides information on their microporosity. The results of the iodine number are presented in Table 2. Experimental values of the iodine number range from 301.95 to 655.81 mg·g−1 and between 419.10 and 781.05 mg·g−1 respectively for CAK-P and CAK-S. The mathematical models of the quadratic regression used to predict the optimal iodine number values are Y2 and Y4 for CAK-P and CAK-S respectively (Equations (10) and (11)):
Table 2. Experimental design matrix for CAK-P and CAK-S preparation using central composite design (CCD).
aExp. val.: experimental values; bPred. val.: predicted values; IN = iodine number.
Table 3 summarizes the analysis of variance (ANOVA) performed to determine the variables that interact significantly during the production of the activated carbons. It appears that the experimental values and predicted values are in agreement, since the determination coefficients R2 of this are close to unity for the two activated carbons (R2 = 91.64% for CAK-P and R2 = 96.57 for CAK-S).
Table 3. Analysis of variance (ANOVA) for the regression mode IN and Yield.
*Significant; Df = degree of freedom; MS = mean square; IN = iodine number; adjR2 = adjusted R2.
Referring to the probable p-values ≤ 0.05, the terms X1, X2, X1X2 and X1X3 for CAK-P and X1, X2, X3 and X1X2 for CAK-S have values less than the said probability p, which shows that all these terms are significant. Therefore, temperature and concentration have a significant effect on the iodine number. This is confirmed by the Pareto diagram analysis shown in Figure 2(a) and Figure 2(b) which clearly leave the corresponding bands cut by a vertical test line. Figure 2(c) and Figure 2(d) show the response surface obtained according to concentration and temperature with a fixed time of 75 minutes. The examination of this figure shows that the iodine number increase with the increase in activation temperature. This is because when the temperature rises, the volatile materials evaporate and some pores are enlarged while more pores are created , thus increasing the sites for adsorption .
In general, the increase in the concentration of activating agents leads to an increase in the iodine number. Tagne et al. (2019)  also reported similar results. For soft lignocellulosic precursors, a high concentration of activating agent leads to less porous activated carbon, which leads to a reduction in the iodine number due to the destruction and/or absence of the appearance of pores. However, the combined action of high concentration and high temperature increase the pores of activated carbon, thereby leading to high iodine number values, as can be seen with experiment 9 in this work.
(a) (b) (c) (d)
Figure 2. Response surfaces: CAK-P (a), CAK-S (b) and Pareto diagrams: CAK-P (c), CAK-S (d) for iodine number.
3.1.2. Activated Carbon Yield
The outcome from the production process of activated carbon is the yield. The performance values shown in Table 2 range from 33.74% - 51.45% and 34.08% - 41.95% respectively for CAK-P and CAK-S. The yield is seen to decrease as the temperature increases. The yields are also seen to be high at low residence time and high concentration of the activating agent, as evidenced by experiment 15 (Table 2). The mathematical models of CAK-P (Y1) and CAK-S (Y3) are given by Equations (12) and (13).
The results of the ANOVA show that the terms X2, X3, X1X3 and X2X3 are significant for the CAK-P and X2, X3, X1X2 and X2X3 are significant for CAK-S. According to Figure 3(c) and Figure 3(d), an increase in carbonization temperature leads to a decrease in yield, while the concentration of activating agents and the residense time have a less significant effect on the response. This decrease in yield can be attributed to the fact that an increase in the calcination temperature increases the amount of carbon combusted and the continued release of volatile material, as is reflected by the loss of material from the precursor. The maximum yield was achieved when the carbonization temperature and residence time were maintained at their lowest values. The results obtained are consistent with those of the literature   and it can be concluded that temperature has a greater impact on yield compared to the concentration of activating agents and the time of calcination on the performance of CAK-P and CAK-S.
3.1.3. Process Optimization and Test Validation
The optimal conditions for the production of the activated carbons were deduced from the maximum yield values and the iodine number. The experimental values of the responses were analyzed by Statgraphic plus 5.0 software and the results are presented in Table 4. We observe that the experimental values are close to the predicted values, implying that the difference between these two values is not significant. However, it is not easy to simultaneously meet the optimal performance and iodine number requirements provided by this software for the preparation of better activated carbons because the zone of interest of the factors is different. The same remark was made by Lekene et al. (2018) . The activated carbons prepared under optimal conditions using iodine number were consequently characterized. At the residence time of 120 min, the optimal conditions obtained were 5 mol·L−1 of H3PO4 and 400˚C carbonization temperature for CAK-P; and 3.5 mol·L−1 of H2SO4 and 700˚C for CAK-S.
(a) (b)(c) (d)
Figure 3. Response Surfaces: CAK-P (a), CAK-S (b) and Pareto diagrams: CAK-P (c), CAK-S (d) for Yield.
Table 4. Validation of the experimental model.
3.2. Physical Characteristics
3.2.1. EDX Analysis
EDX analysis of the activated carbons is shown in Figure 4 while Table 5 gives the weight percentages of the elements present in them. This analysis was coupled with the scanning electron microscopy (SEM/EDX) captured on a selected area 1 section of the SEM image. Carbon appears to be the most representative element in CAK-P and CAK-S, this is explained by the increase in the degree of aromatization during carbonization. The small amounts of the other elements originated from the lignocellulosic precursor, and they become perceptible on the departure of volatile substances during carbonization.
3.2.2. Physico-Chemical Characteristics of the Activated Carbons: pH, pHpzc, Bulk Density, Moisture Content and Boehm Titration
The physico-chemical properties that were determined for the prepared ACs are presented in Table 6. The pH values of 3.34 and 6.91 respectively for CAK-P and CAK-S show that the two activated carbons are acidic. pHpzc values are 2.23 and 6.9 respectively for CAK-P and CAK-S (Table 6, Figure 5). These are the pH values where the two carbons exist at zero electrical potential. This means that for solutions where pH > pHpzc, the electrostatic attraction is favourable for the cationic pollutant because the surface of the material is deprotonated. On the other hand, when pH < pHpzc, the electrostatic attraction is favourable for the anionic pollutant since the surface of the material is protonated .
Figure 4. EDX analysis of activated carbons: CAK-P (a) and CAK-S (b).
Figure 5. Determination of the pH at the point of zero charge (pHpzc) for CAK-P and CAK-S.
Table 5. Elemental content of the activated carbons in percentage.
Table 6. Some physico-chemical characteristics of CAK-P and CAK-S.
The bulk density of an activated carbon provides information on the microcristalline structure and strong inter-linking bonds . The moisture content of these ACs is in line with the standard of ≤5% as stipulated by the American Work and Water Society.
3.2.3. Specific Surface Area
The textural structure of the activated carbons was determined by the BET method, using relative pressure (P/P0) from 0.005 to 0.4. Figure 6 and Figure 7 show the analyses on porosity for CAK-P and CAK-S respectively. According to the classification of IUPAC and taken up by Francisco et al. (2018) , the adsorption and desorption isotherm of N2 is Type I and H4 hysteresis illustrates the coexistence of micropores and mesopores. The result shows that these two activated carbons have more micropores than mesopores. The results are summarized in Table 7. The measurement of specific surfaces according to the Brunauer-Emmett-Teller analysis (BET) yielded 457.38 m2·g−1 and 490.62 m2·g−1 respectively for CAK-P and CAK-S. Thus the microporous surfaces are 522.55 m2·g−1 and 570.56 m2·g−1 and the mesoporous ones are 9.1974 m2·g−1 and 10.779 m2·g−1 respectively for CAK-P and CAK-S. Therefore, the two activated carbons have a similar texture, but the one activated with sulphuric acid has better developed pores. In view of these results it is obvious that these ACs can be used to adsorb small and medium-sized molecules from aqueous solution. Table 7 shows the values obtained.
Table 7. Summary report of micrometrics analysis for CAK-P and CAK-S.
Figure 6. CAK-P: Adsorption-desorption curves of nitrogen at 77.13 K (a), Sample mesopores (b) and micropores (c) Size distribution.
Figure 7. CAK-S: Adsorption-desorption curves of nitrogen at 77.13 K (a), Sample mesopores (b) and micropores (c) Size distribution.
3.2.4. FTIR Spectra
Figure 8 shows the IR-TF spectra recorded in the 3500 to 580 cm−1 region. Low peaks at 2915.36 and 2861.41 cm−1 on the biomass spectrum (CNK) indicate the symmetric and asymmetric stretching vibration of methyl group  . The absence of these peaks on the spectra of the CAK-P and CAK-S activated carbons shows that during the activation with phosphoric and sulphuric acids, there was a reaction resulting in the formation of other oxygenated functional groups. The signal at 2368.65 cm−1 may be due to O-C-O and the peak at 1731.65 cm−1 corresponds to the C-O group of esters and lactones. The appearance of the peak at 1040.36 cm−1 corresponds to the C-O group of primary, secondary, tertiary alcohols and ethers, while that at 1257.96 cm−1 indicates the presence of the same group C-O group of carboxylic acids and esters or the P-O link of ester phosphates or O-C of the P-O-C bridge . Between 800 and 650 cm−1 we have the C-H bond of aromatic polynuclear systems .
3.2.5. Scanning Electron Microscopy (SEM) Analysis
The results obtained from the SEM analyses of the activated carbons are presented in Figure 9. These pictures show that the two activated carbons have
Figure 8. FTIR spectra for biomass (CNK) and activated carbons (CAK-P and CAK-S).
Figure 9. SEM images of the activated carbons (500 um, selected zone 1).
almost identical porosity. It should be remembered that a developed porosity increases the specific surface area of the activated carbons  and therefore increases the number of sites available for the fixing of pollutants during adsorption. It is also observed that both activated carbons have a heterogeneous surface.
In this study, a triacid and diacid were used as activating agents in the preparation of ACs using shea residues (Vitellaria paradoxa) as precursors. Central composite design was used to study the influence of preparation parameters, concentration of activating agents, carbonization temperature and residence time, on expected responses, yield and iodine number. Under the optimal conditions recommended by Statgraphic plus 5.0, thus, 5 mol·L−1, 400˚C and 120 min for CAK-P and 3.5 mol·L−1, 700˚C and 120 min for CAK-S, the iodine numbers were 709.45 mg·g−1 and 817.36 mg·g−1. The results of pHpzc, FTIR and Boehm titration indicate that CAK-P and CAK-S ACs prepared under these conditions are acidic. They have a microporous surface of 522.55 and 570.65 m2·g−1 respectively for CAK-P and CAK-S. Comparative analysis of the physic-chemical characteristics shows that these two ACs are of good quality, the better one being the one activated with sulphuric acid.
The authors greatly thank assistance rendered by M. Kamdem Tamo Arnaud, Doungmo Giscard and Dr. Tagne Tiegam Rufis for their help in the characterization of activated carbons. We also thank the members of our research team for their idea and contribution throughout this manipulation and the writing of the article.
 Balogoun, C.K., Bawa, M.L., Osseni, S. and Aina, M. (2015) Préparation des charbons actifs par voie chimique à l’acide phosphorique à base de coque de noix de coco. International Journal of Biological and Chemical Sciences, 9, 563-580.
 Vintila, T., Ioana, I., Tangne, T.R., Adriana, R.W., Calin, J. and Anagho, S.G. (2019) Residual Biomass from Food Processing Industry in Cameroon as Feedstock for Second-Generation Biofuels. Bioresources, 14, 3731-3745.
 Zhang, T., Walawender, W.P., Fan, L.T., Maohong, F., Daren, D. and Brown, R.C. (2004) Preparation of Activated Carbon from Forest and Agricultural Residues through CO2 Activation. Chemical Engineering Journal, 105, 53-59.
 Mohammed-Khah, A. and Ansari, R. (2009) Activated Charcoal: Preparation, Characterization and Applications: A Review Article. International Journal of Chemical and Technology Research, 1, 859-864.
 Haykiri-Acma, H., Yaman, S. and Kucukbayrak, S. (2006) Gasification of Biomass Chars in Steam-Nitrogen Mixture. Energy Conversion and Management, 47, 1004-1013.
 Vargas, A.M.M., Cazetta, A.L., Garcia, C.A., Moraes, J.C.G., Nogami, E.M., Lenzi, E., Costa, W.F. and Almeida, V.C. (2011) Preparation and Characterization of Activated Carbon from a New Raw Lignocellulosic Material: Flamboyant (Delonix regia) Pods. Journal of Environmental Management, 92, 178-184.
 Tchakala, I., Bawa, M.L., Djaneye-Boundjou, G., Doni, K.S. and Nambo, P. (2012) Optimisation du procédé de préparation des charbons actifs par voie chimique (H3PO4) à partir des tourteaux de Karité et des tourteaux de Coton. International Journal of Biological and Chemical Sciences, 6, 461-478.
 Tchuifon, T.D.R., Anagho, S.G., Nche, G.N. and Ketcha, J.M. (2015) Adsorption of Salicylic and Sulfosalicylic Acid onto Powdered Activated Carbon Prepared from Rice and Coffee Husks. International Journal of Current Engineering and Technology, 5, 1641-1652.
 Idris-Hermann, K.T., Tchuifon, T.R., Doungmo, G. and Anagho, S.G. (2018) Preparation and Characterization of Activated Carbons from Bitter Kola (Garcinia kola) Nut Shells by Chemical Activation Method Using H3PO4; KOH and ZnCl2. Chemical Science International Journal, 23, 1-15.
 Bopda, A., Tchuifon, D.R.T., Nche, G.N., Doungmo, G. and Anagho, S.G. (2019) Non-Linear Equilibrium and Kinetic Study of the Adsorption of 2,4-dinitrophenol from Aqueous Solutionusing Activated Carbon Derived from a Olives Stones and Cotton Cake. African Journal of Environmental Science and Technology, 3, 365-380.
 Djekonbé, P. (2019) Domestication du karité (vitellaria paradoxa c.f. gaertn) au Tchad: Evaluation de vergers et etude de la variabilite morphologique et biochimique des fruits. Thèse Doctorat/PhD, Université de Dschang, Dschang, 180 p.
 Djekota, C., Mouga, M., Djimramadji, A., Djelassem, B., Mbayngoné, E., Maiga, R.D., Rimgoto, K. and Noubady, D. (2014) Potentiel karité au Tchad (Vitellaria paradoxa C.F. Gaertn. Subsp. Paradoxa). Journal of Animal & Plant Sciences, 23, 3646-3656.
 Annadurai, G., Babu, S., Nagarajan, G. and Ragu, K. (2000) Use of Box-Behnken Design of Experiments in the Production of Manganese Peroxidase by Phancrochate chrysosporium (MTCC 767) and Decolorization of Crystal Violet. Bioprocess Engineering, 23, 715.
 Tagne, T.R.F., Ioana, I., Anagho, S.G. and Alin-cristian, M. (2019) Optimization of the Activated Carbon Preparation from Avocado Seeds, Using the Response Surface Methodology. Revista de Chimie (Bucharest), 70, 410-416.
 Armineh, A., Mohammad, A.M. and Mokhtar, A. (2011) Application of Response Surface Methodology for Optimization of Reactive Blue 19 Dye Removal from Aqueous Solutions Using Pulp and Paper Sludge. Fresenius Environmental Bulletin, 20, 1-11.
 Lekene, R.B.N., Ndi, J.N., Rauf, A., Kouotou, D., Placide, D.B.B., Bhanger, M.I. and Ketcha, J.M. (2018) Optimization Conditions of the Preparation of Activated Carbon Based Egusi (Cucumeropsismannii Naudin) Seed Shells for Nitrate Ions Removal from Wastewater. American Journal of Analytical Chemistry, 9, 439-463.
 Atemkeng, D.C., Kamgaing, T., Tchuifon, T.D.R., Doungmo, G., Amola, A.L., Kamdem, T.A. and Anagho, S.G. (2020) Chemical Preparation and Physicochemical Characterization of Powdered Activated Carbons Based on Safou (Dacryodes edulis) Seeds. Journal of Materials and Environmental Science, 11, 896-910.
 Ravikumar, K., Krishnan, S., Ramalingam, S. and Balu, K. (2007) Optimization of Process Variables by the Application of Response Surface Methodology for Dye Removal Using Novel Adsorbent. Dyes and Pigments, 72, 66-74.
 Arunachalam, R. and Annadurai, G. (2011) Optimizred Response Surface Methodology for Adsorption of Dyestuff from Aqueous Solution. Journal of Environmental Science and Technology, 4, 65-72.
 Maazou, S.D.B., Hima, H.I. and Mousbahou, M. (2017) Elimination du chrome par du charbon actif élaboré et caractérisé à partir dela coque du noyau de Balanites aegyptiaca. International Journal of Biological and Chemical Sciences, 11, 3050-3065.
 Guo, Y. and Rockstraw, D.A. (2006) Physical and Chemical Properties of Carbons Synthesized from Xylan, Cellulose, and Kraft Lignin by H3PO4 Activation. Carbon, 44, 1464-1475.
 Danish, M., Hashim, R., Ibrahim, M.N.M. and Sulaiman, O. (2014) Optimized Preparation for Large Surface Area Activated Carbon from Date (Phoenix dactylifera L.) Stone Biomass. Biomass and Bioenergy, 61, 167-178.
 Tan, I., Ahmad, A. and Hameed, B. (2008) Optimization of Preparation Conditions for Activated Carbons from Coconut Husk Using Response Surface Methodology. Chemical Engineering Journal, 137, 462-470.
 Kamgaing, T., Doungmo, G.F., Tchieno, M.M., Kouonang, J.G. and Ketcha, J.M. (2017) Kinetic and Isotherm Studies of Bisphenol a Adsorption onto Orange Albedo (Citrus sinensis) Sorption Mechanisms Based on the Main Albedo Components Vitamin C, Flavones Glycosides and Caortenoids. Journal of Environmental Science and Health, Part A, 52, 757-769.
 Pongener, C., Kibami, D., Rao, K., Goswamee, R.L. and Sinha, D. (2015) Synthesis and Characterization of Activated Carbon from the Biowaste of the Plant Manihot esculenta. Chemical Science Transactions, 4, 59-68.
 Francisco, J.S., Katie, A.C. and Matthias, T. (2018) Characterization of Micro/Mesoporous Materials by Physisorption: Concepts and Case Studies. Accounts of Materials and Surface Research, 3, 34-50.
 Ngakou, C.S., Ngomo, H.M. and Anagho, G.S. (2018) Batch Equilibrium and Effects of Ionic Strength on Kinetic Study of Adsorption of Phenacetin from Aqueous Solution Using Activated Carbon Derived from a Mixture of Ayous Sawdust and Cucurbitaceae Peelings. Current Journal of Applied Science and Technology, 26, 1-24.
 Trachi, M., Bourfis, N., Benamara, S. and Gougam, H. (2014) Préparation et caractérisation d’un charbon actif à partir de la coquille d’amande (Prunus amygdalus) amère. Biotechnology, Agronomy. Society, and Environment, 18, 492-502.