Kinetic, Isotherm and Thermodynamic Studies of the Adsorption of Thymol Blue onto Powdered Activated Carbons from Garcinia cola Nut Shells Impregnated with H 3 PO 4 and KOH: Non-Linear Regression Analysis

In this work, activated carbons (ACs) prepared by chemical activation of garcinia cola nut shell impregnated with H 3 PO 4 (CBH 2/1 ) and KOH (CBK 1/1 ) were used to study the kinetics, equilibrium and thermodynamics of the adsorption of thymol blue from aqueous solution. The characterization of ACs showed the BET measurements gave surface area and total pore volume re-spectively of 328.407 m 2 ∙g −1 and 0.1032 cm 3 ∙g −1 for CBH 2/1 and 25.962 m 2 ∙g −1 and 0.03 cm 3 ∙g −1 for CBK 1/1 ; elemental analysis showed a high percentage of carbon in both ACs. Influence of parameters such as initial pH, contact time, adsorbent mass, initial concentration, obtained using the three-parameter Fritz-Schlunder equation were 32.147 mg∙g −1 for CBH 2/1 and 67.494 mg∙g −1 for CBK 1/1 . The values of the model parameters g and m FS respectively, obtained using the Redlich-Peterson and Fritz-Schlunder III isotherms below 1, showed that the adsorption of thymol blue by the ACs occurred on heterogeneous surfaces. Thermodynamic analyses of the data of the adsorption of thymol blue onto ACs revealed that the adsorption process was temperature dependent, endothermic and spontaneous.


Introduction
The intensive industrialization considered as key to economic development, which started in the last century is at the origin of the proliferation of various pollutants in the environment, and most of them are difficult to eliminate [1].
The wastes resulting from industrial activities represent a danger to man and his environment, due to their very low biodegradability and the slowness of their elimination [1]. The global increase in water pollution has increased awareness among the public, and more attention is being focused on combating the prob-and teratogenic and they are also toxic to human beings, fish species, and microorganisms [9]. The majority of these pollutants are not biodegradable and can constitute risk factors to our health, as well as being a nuisance to the environment. The presence of these pollutants in water, even at very low concentrations is very undesirable because it is highly visible. In fact, the presence of dyes in flowing streams causes a high reduction in the activity of photosynthesis [10].
Thymol blue is an example of such a dye. It is used as a pH indicator in analytical chemistry. When inhaled it leads to irritation of the respiratory tract, or when it gets in contact with the skin, it is very painful [11]. These effects have become a major concern to the society. For this reason, many studies have focused on the use of depollution techniques to eliminate pollutants from wastewaters. The majority of these techniques are very expensive and difficult to implement especially when they are applied to effluents of high flow rates and having low concentrations of pollutants.
Adsorption has proved to be a more effective and cheaper process, and its equipment is easy to implement [1]. However, the choice of adsorbent remains a major worry because most adsorbent materials such as alumina, zeolites, activated clays, and silica gels are also used as catalyst support materials [12].
Therefore, activated carbons are the most commonly used adsorbents because they contain large surface areas, well-developed porosity and contains various functional groups on their surfaces. Activated carbons prepared by chemical activation of lignocellulosic materials have been used for the adsorption of dyes, various organic compounds and antipyretic/analgesic drugs [13] [14] [15]. The advantage of chemical activation is that the pyrolysis stage of the synthesis is achieved at low temperature, thereby greatly reducing the cost of fabrication. Over the last few years, much interest has been shown by researchers on the technique of adsorption, and they have produced low cost activated carbons from lignocellulosic materials which are abundant agricultural wastes, to remove organic compounds, metals and dyes from solutions of waste streams. The costs of activated carbons prepared from lignocellulosic materials are very low compared to the cost of commercial activated carbons [9]. Many activated carbons have been produced from various carbonaceous materials resulting from agricultural wastes such as coconut shells [9], cola nut shells [13], rice and coffee husks [14], Ayous sawdust and Cucurbitaceous peelings [15], hibiscus sabdariffa and stalk biomass [16]. In this work, we will prepare ACs using shells of Garcinia cola (bitter cola) nuts [17], due to the availability and cost free of the material as well as its high carbon content. However, little or no work has been reported on the production of activated carbons using the garcinia cola nut shells.
The aim of this work was to evaluate the adsorption properties of activated carbons prepared from garcinia cola nut shells impregnated by H 3 PO 4

and KOH
for the removal of thymol blue from aqueous solution. Thermodynamics studies were carried out, and the experimental data were analyzed using the two-parameter

Chemicals: Adsorbate (Thymol Blue)
The adsorbent used in this study was thymol blue, which is an anionic dye.
Thymol blue, also called thymol sulphonephthalein is a brownish-green or reddish-brown crystalline powder that is used as a pH indicator and for biological activities. It is insoluble in water, but soluble in alcohol and dilute alkali solutions. Its colour changes from red to yellow at pH 1.2 -2.8 and from yellow to blue at pH 8.0 -9.6. A stock of 500 mg•L −1 of thymol blue was prepared by dissolving 125 mg in 100 mL of dilute alkali solution (Na + + HO − ) using distilled water and then completing to 250 mL. Working solutions of various concentrations were freshly prepared from this stock solution. HCl or NaOH of concentration 0.1 mol•L −1 was used to adjust the initial pH of the solution. NaCl was used to study the influence of ionic strength on the adsorption process

Absorbents
The ACs used had been prepared using optimized experimental conditions by chemical activation of garcinia cola (bitter cola) nut shells. The AC CBH 2/1 was prepared by activating garcinia cola nut shells with H 3 PO 4 in a ratio 2/1, while CBK 1/1 was prepared by activating with KOH in a ratio 1/1. The chemical activation was followed by calcination in a furnace set at 400˚C for 1 hour and a heating rate of 5 ˚C•min −1 as described by Idris-Hermann et al. [17].

Characterizations
Ultimate analyses of raw material and ACs were performed for determining the C, H, N, and S content using CHNSO elemental analyses from HEKAtech CHNS Analyzer. The oxygen content was calculated by the difference on 100%. Textural analysis was performed by N 2 -sorption measurements at 77 K using multipoint porosimeter automatic analyzer (Sorptomatic 1990, Porotec, GmbH), system ADP2005 Extended Report software using to calculated BET-surface area, total pore volume and pore size…, with a relative pressure range of 0.05 < P/P 0 < 0.25, with an overnight pretreatment of 150˚C for 5 hours. Total pore volume (V tot ) was estimated to be the liquid volume of N 2 at the high relative pressure P/P 0 = 0.95. The micropore volume (V micro ) was obtained using the D-R method and average pore diameter (D P ) was calculated from D P = 4V T /S BET [18].
The micropore size distribution parameter was calculated according to the simplified equation of the Horvath-Kawazone (HK) [19].
The BET specific surface area (S BET ) was calculated using the equation below: where V m is the volume occupied by a monolayer, N A is the Avogadro constant 6.023 × 10 23 mol −1 , σ is the surface occupied by the N 2 molecule (16.2 Å 2 ); m is the mass of the sample (g) and V M is the molar volume of N 2 (22414 cm 3 •mol −1 ). In our earlier work, Idris-Hermann et al. [17], the different ACs had been characterized by proximate analysis, bulk density, pH of zero point charge (pH zpc ), methylene blue number, iodine number, surface functional groups, Fourier transform infra-red spectroscopy (FTIR), scanning electron microscopy (SEM), and power X-ray diffraction (XRD).

Batch Adsorption Experiments
C 0 is the initial concentration of the dye, C e is the concentration of the dye at the time t, V is the volume of the adsorbate solution, and m is the mass of the adsorbent.

Effect of Initial pH
To determine the effect of pH, the adsorption of thymol blue by the ACs were investigated over a pH range of 2 to 7 at ambient temperature. For each adsorbent, 50 mg was treated with 20 mL of aqueous solution of 50 mg•L −1 of thymol blue dye. The initial pH of the solution was adjusted by adding HCl or NaOH solution.

Effect of Contact Time
To determine the effect of agitation time on the adsorption process, 50 mg of ground adsorbent was agitated in a 20 mL solution of thymol blue of initial concentration 50 mg•L −1 for different contact times, varying between 0 and 105 minutes. After each contact time t, the solution was rapidly filtered and the residual concentration of thymol blue, C e determined by UV/Vis spectrophotome-ter. The amount (Q e ) of thymol blue adsorbed was calculated using Equation (3).

Effect of Amount of Adsorbent
In this set of experiments, different masses of the adsorbents 30, 50, 70, 90, 110, 130, 150 and 170 mg were treated with 20 mL solution of thymol blue of initial concentration 50 mg•L −1 at the optimum time obtained from the effect of contact time experiments and the optimum pH obtained from the effect of pH experiments.

Effect of Ionic Strength
To determine the effect of ionic strength on the adsorption process, 50 mg of each adsorbent was agitated in a 20 mL solution of thymol blue at pH 2 of initial concentration 50 mg•L −1 containing the salt, NaCl of concentrations ranging between 0.01 and 0.06 mol•L −1 . At equilibrium time (1 hour), the solution was rapidly filtered and the residual concentration of thymol blue determined by spectrophotometer. The amount (Q e ) of thymol blue adsorbed was calculated.

Kinetics of Adsorption Studies
To deduce the controlling mechanisms for the adsorption of thymol blue onto the ACs, data from the kinetics experiments were analyzed using the different non-linear kinetic models presented in Table 1 below. [20] ( ) Intraparticle Diffusion:

Adsorption Isotherm Studies
To study the equilibrium of the adsorption of thymol blue onto the ACs, data from the equilibrium experiments were analyzed using the following two-to five-parameter, non-linear equilibrium models as presented in Table 2 below.

Non-Linear Regression Analysis
In this study, non-linear regression was applied using Microsoft Excel Solver function and using Origin Pro 9, 64 bit for fitting the curve. The best fit for experimental data was determined from the coefficient of determination (R 2 ), residual root mean square error (RMSE), Chi-square test (χ 2 ), Sum Square of Er-  Table 3.

Thermodynamic Studies
The effect of temperature on the adsorption of thymol blue dye onto the ACs was evaluated from thermodynamic studies. 20 mL of thymol blue solution of concentration 50 mg•L −1 (pH = 2) was placed in contact with 50 mg each of AC in a polypropylene flask and placed on a magnetic stirrer/heater and heated with stirring for the equilibrium time 1 hour, at temperatures of 35˚C, 45˚C, 55˚C and 65˚C (308, 318, 328 and 338 K). After the equilibrium time, the concentration of the remaining thymol blue was determined by spectrophotometry, and the maximum amount adsorbed calculated from Equation (3). The thermodynamic quantities, change in Gibbs free energy (∆G 0 ), enthalpy change (∆H 0 ) and entropy change (∆S 0 ) were calculated.

Textural Characterization of the Adsorbents
Textural characteristics of all ACs were studied. As it can be observed, CBK 1/1 25.962 m 2 •g −1 gave the lower value of S BET compared to CBH 2/1 328.407 m 2 •g −1 calculated by Equation (1). The lower porosity of CBK 1/1 is confirmed by the higher value of average pore diameter 4.622 nm, while CBH 2/1 shows D p value of 1.257 nm smaller than 2.0 nm. However, the V micropores of CBH 2/1 is higher than that of CBK 1/1 , which corroborates the microporous characteristic of CBH 2/1 .
The values of total pore volume and BET surface area upon H 3 PO 4 activation are larger because H 3 PO 4 prevents the accumulation of tars on the surface of the carbon. Also, according to pore evolution theory, hydrogen and oxygen compounds in the lignocellulosic material are removed in the form of volatiles materials during decomposition. These volatiles diffuse out of the carbon matrix, leading to the development of pores.
The results obtained by scanning electron microscopy (SEM), power X-ray diffraction (XRD) had been presented in our former work, Idris-Hermann et al. [17].

Elemental Analysis
The elemental composition (CHNS) of the raw material (CB) and the two ACs were studied and the results are presented in Table 4 below. It was obvious from the table that H 3 PO 4 and KOH used as activating agents have significant influence in producing activating carbons with high ratios of C/H and C/N. As expected, the ACs have much higher contents of carbon and lower contents of oxygen and hydrogen compared to the initial raw material. As a matter of fact, the chemical activation during the heat treatment may accelerate the chemical changes in carbonaceous material and enable the removal of volatile compounds. This could be attributed to the fact that H 3 PO 4 and KOH selectively stripe H and O away from the garcinia cola nutshells as H 2 O and H 2 rather than CO, CO 2 or hydrocarbons [36] [37]. Therefore, the carbonaceous material becomes richer in carbon. However, elemental analyses of the activated carbons indicate the absence of sulphur for all of them.

Effect of Initial pH
The effect of initial pH on the adsorption of thymol blue dye by the two ACs was studied by varying the initial pH of the thymol blue solution. The result is shown in Figure 1. According to this figure, the removal of dye from aqueous solution is highly dependent on the pH of the solution. In Figure 1, the results show that, the maximum capacity of adsorption is obtained at pH = 2 and favoured for pH < pH PZC . The pH at the point of zero chage, pH PZC was obtained for the two ACs as 7.751 for CBH 2/1 and 8.324 for CBK 1/1 [17]. When the pH increases, the capacity of adsorption decreases from 14.809 to 3.684 mg•g −1 and 15.476 to 2.514 mg•g −1 for CBH 2/1 and CBK 1/1 respectively. This can be explained by the fact that, when the pH of solution is lower than the pH PZC, more protons will be available for protonating the surface of the ACs, thereby increasing the electrostatic at- and -OH groups of the ACs surface due to the increased concentration of the hydroxide ions. The electrostatic force of repulsion between thymol blue and adsorbent is at the origin of the decreasing capacity of adsorption. A similar result was obtained from the adsorption of amaranth red by pineapple peelings and coconut shells 1 .

Effect of Contact Time
The effect of contact time on the adsorption of thymol blue is given in Figure 2 below. The figure shows the evolution of the adsorption capacity of thymol blue with time of contact at ambient temperature. It can be observed that, the adsorption process took place in three phases. Rapid adsorption of thymol blue took place in the first 5 minutes for the two activated carbons. This is explained as follows: at the start of the adsorption process, all the vacant sites on the adsorbent were completely available and the thymol blue molecules were fast attached. The same tendency had been obtained by Ndifor-Angwa for and others when they adsorbed amaranth red onto treated and untreated lignocellulosic materials [1]. Thereafter, the rate of adsorption deceased gradually with the progress of adsorption within 5 to 60 minutes for the two adsorbents, and reached equilibrium in about 60 minutes. The adsorption capacity at these equilibrium points are 14.974 mg•g −1 for CBH 2/1 and 15.839 mg•g −1 for CBK 1/1 . This can be explained by the fact that the available sites for adsorption have been greatly reduced by the thymol blue molecules already attached there. Also, the already adsorbed molecules tend to repel any other approaching molecules as they are similarly charged, thereby reducing the rate of adsorption. The adsorption process can be attributed to aromatic condensation due to π-π interactions between thymol blue molecules and the surface of the adsorbent. During the aromatic condensation, it is possible that surface functional groups of adsorbent do not occupy the entire nucleophilic sites of thymol blue. Figure 2 also shows that there is no significant change in the rate of adsorption of thymol blue after equilibrium time. This corresponds to the third phase, which lies between 60 and 105 minutes, in which for CBH 2/1 and CBK 1/1 the adsorption capacities are almost constant with time.

Effect of Adsorbent Dosage
The effect of adsorbent dosage on the removal of thymol blue was studied and the results of this study are shown in Figure 3 below. In this figure, the percentage removal of thymol blue by CBH 2/1 and CBK 1/1 increased as the adsorbent dosage increased. The removal percentage of thymol blue increased from 63.519 to 76.623 and 66.445 to 83.058 for CBH 2/1 and CBK 1/1 respectively. This is because an increase in adsorbent mass leads to an increase in surface area, and hence an increase in the number of available adsorption sites for thymol blue to occupy during the adsorption process [38]. A similar trend was reported in Hameed et al. [39] for the adsorption of chromotrope dye onto activated carbons obtained from the seeds of various plants.

Effect of Initial Concentration
The effect of initial concentration on the adsorption of thymol blue was investigated and the results are show in Figure 4 below. In Figure 4, we observe that the quantity of thymol blue adsorbed at equilibrium increases from 10.519 to 34.400 mg•g −1 and 10.893 to 34.681 mg•g −1 for CBH 2/1 and CBK 1/1 respectively.  The increasing adsorption capacities of the ACs with increase in concentration of thymol blue can be due to π-π interactions between the organic compounds and functional groups of the carbon surface. The π-π interactions in most cases are responsible for the mechanism of adsorption of aromatic compounds [14]. This result can also be explained by saying that, increasing the initial concentration of thymol blue leads to an increase in the driving force for mass transfer, and hence an increase in the rate at which thymol blue molecules pass from solution to the particle surface [38]. The quantity of thymol blue adsorption onto all samples of ACs did not show a plateau, suggesting that there was no monolayer formation on the surfaces of the adsorbents. The increase in adsorption capacity with increase in initial concentration indicates that, the ACs have a high potential for the removal of thymol blue from its solutions.

Effect of Ionic Strength
The ionic strength can be viewed as one of the factors that controls nonelectrostatic and electrostatic interactions between adsorbate and the surface of the adsorbent. The effect of ionic strength is also important in the study of adsorption of dye onto adsorbents. Ionic strength effect on the adsorption of thymol blue was investigated and the results are shown in Figure 5. In this figure, the adsorption capacities of dye at equilibrium decreased from 13.347 to 12.072 mg•g −1 for CBH 2/1 and 13.336 to 12.271 mg•g −1 for CBK 1/1 when the ionic strength increased from 0.01 to 0.06 mol•L −1 . These results can be explained by the fact that, the surface of ACs becomes less accessible to thymol blue molecules when the quantity of NaCl salt in solution increases. For this reason, the amount of thymol blue adsorbed decreases. This result is in agreement with the literature, Journal of Encapsulation and Adsorption Sciences which reports that when the electrostatic attraction forces between the adsorbent surface and adsorbate ions are attractive, an increase in ionic strength will decrease the adsorption capacity [40].

Batch Kinetic Studies
Pseudo-first order, pseudo-second order, intraparticle diffusion, Elovich, fractional power, Avrami and Tobin kinetics models are used to investigate and describe the mechanism of the adsorption process (see equation in Table 1). These kinetic models have been presented in different non-linear forms and Table 5 presents the calculated parameters.
According to Table 5, the non-linear models of pseudo-second order, Avrami and Tobin seem to be the most appropriate in adequately describing the adsorption phenomenon for the two ACs, CBH 2/1 and CBK 1/1 studied. As shown in Table 6, the coefficient of determination, R 2 of these models are very good (R 2 > 0.95). However, the relatively higher R 2 values and the lower χ 2 , ERRSQ values and other relatively lower error functions for the Avrami model compared to the pseudo-second order and Tobinmodels confirm that the Avrami non-linear kinetic expression is best suited to describe the adsorption kinetics of thymol blue on the two ACs. The Avrami exponent (n AV ) is a fractional number that is related to the possible changes in the adsorption mechanism that takes place during the adsorption process [41]. Therefore, instead of following only an integer-kinetic order, the adsorption mechanism could follow multiple kinetic orders. That is, during adsorption, the order of reaction can change from one value to another [41]. In this study, n AV exponent showed that the adsorption process was of the multiple kinetic order. The low values of χ 2 and other error functions Table 5. Kinetic data obtained by non-linear fitting analysis (see equations in Table 1 and Table 3 Table 6. Adsorption of thymol blue: non-linear fitting analysis for two-parameter isotherms (see Equations in Table 2 and Table  3). of the Tobin and pseudo-second order models suggest that these two models equally describe the kinetics of thymol blue uptake onto the ACs. The values of V lower than unity obtained for the fractional power kinetic model suggest that fractional power may also describe the adsorption kinetics of the thymol blue-AC system. Adsorption data obtained from the pseudo-second order kinetic model thus indicate that the adsorption of thymol blue onto the ACs CBH 2/1 and CBK 1/1 may also be of a chemical nature.

Adsorption Isotherms
Adsorption isotherms describe the relationship between the quantity of adsorbate adsorbed by the adsorbent (Q e ) and the concentration of adsorbate remaining at equilibrium (C e ). This part aims at finding the models that can accurately describe the data obtained from the experiments of the adsorption thymol blue onto CBH 2/1 and CBK 1/1 , and then to specify the parameters that can be determined. Eight models of two-parameter isotherms, seven models of three-parameter isotherms, two models of four-parameter isotherms and one model of five-parameter isotherms were studied (see equation in Table 2).

Two Parameter Models
The capabilities of the equations of the two-parameter isotherms to model equilibrium adsorption data were examined and Table 6 also presents the calculated parameters of the adsorption isotherms and their characterizations obtained using the non-linear model analysis. It can be observed from the table that for the two ACs, the coefficient of determination for all the models is very good (R 2 > 0.96) except for the Dubinin-Radushkevich model. Thus, in this table the lower values of chi-square (χ 2 ), RMSE, ERRSQ, ARE and other error functions for each isotherm model will determine the best model. Among the two-parameter equations tested, and based on the low values of the error functions and the high determination coefficient, the best representation of the experimental data of the adsorption of thymol blue is obtained using the isotherms of Freundlich, Halsey and Kiselev for CBH 2/1 ; and those of Freundlich, Halsey and Langmuir for CBK 1/1 . In this table, the Halsey adsorption isotherm confirms the multilayer adsorption and the heteroporous nature of the adsorbents. The values of 1/n for the Freundlich model are equal and lower than 1 (1.000 for CBH 2/1 and 0.871 for CBK 1/1 ), thereby indicating the high affinity between these ACs and thymol blue solution. This also indicates that adsorption took place through a physical process, and that there was a greater heterogeneity of the adsorption sites on the adsorbent. The 1/n = 1 for CBH 2/1 indicates that the isotherm is linear, and is of the type C, while for CBK 1/1 , the value of 1/n < 1 indicates the isotherm to be concave, and hence of type L.
For CBH 2/1 the negative value of k n for the Kiselev isotherm shows that there is no complex formed between the adsorbed thymol blue molecules and the AC.
On the other hand, the positive value of k n with CBK 1/1 shows the formation of a complex between the adsorbed thymol blue molecules and the AC. The Jovanovic isotherm describes an approximation of a monolayer localized adsorption without lateral interactions. The results of the Temkin model show that the change in the adsorption energy for the two ACs is higher than zero, thus indicating that the process of adsorption is exothermic. For the Langmuir isotherm, the values of the separation factor falling between 0 -1 ( Figure 6) show a favorable adsorption of thymol blue by the two ACs. This result is in agreement with what was already observed with the Freundlich isotherm.

Three Parameter Models
The capabilities of the equations of the three-parameter isotherms to model equilibrium adsorption data were examined. Table 7 presents the calculated parameters of the adsorption isotherms and their characterizations obtained using the non-linear model analysis.
It can be observed in Table 8 that for the two ACs, the coefficient of determination for all models is very good (R 2 > the 0.96). Thus, in this table, the lower values of the chi-square (χ 2 ), RMSE, ERRSQ, ARE and other error functions for each isotherm model and high value of determination coefficient will determine the best model. According to these values, the best-fitted adsorption isotherm models were found in the order: Fritz-Schlunder III > Redlich-peterson > Radke-Prausnitz > Jossens > Toth > Langmuir-Freundlich > Hill for CBH 2/1 predicting the adsorption to be on heterogeneous surfaces; Fritz-Schlunder III > Redlich-peterson > Radke-Prausnitz > Hill > Langmuir-Freundlich > Jossens > Toth for CBK 1/1 predicting the adsorption process as a cooperative phenomenon Figure 6. Evolution of separation factor against the initial concentration of thymol blue. Journal of Encapsulation and Adsorption Sciences Table 7. Adsorption of thymol blue: Non-linear fitting analysis for three-parameter isotherms (see Equations in Table 2 and Table 3).  Table 8. Adsorption of thymol blue: Non-linear fitting analysis of four-and five-parameter isotherms (see Equations in Table 2 and Table 3). Redlich-Peterson and Fritz-Schlunder III isotherms are lower than 1, which means that the adsorption of thymol blue by the ACs cannot be reduced to the Langmuir isotherm. This implies that the adsorption already predicted to be on heterogeneous surfaces is also confirmed by the Jossens model. For the two ACs, it was observed that n was greater than unity. This suggests that the interactions between the adsorbed molecules are stronger than the adsorptive potential.

Four and Five Parameter Models
Isotherms of four-parameters (Fritz-Schlunder IV and Baudu) and five-parameters (Fritz-Schlunder V) were used to model the equilibrium data of adsorption of thymol blue and Table 8 presents the calculated parameters of the adsorption isotherms and their characterizations obtained using non-linear model analysis.
It can be observed from Table 8 that the coefficients of determination for all models are very good (R 2 > 0.99). Based on the low values of ERRSQ and other error functions, high coefficients of determination, the two four-parameter isotherms gave a good fitting of the experimental results for both CBH 2/1 and CBK 1/1 . The maximum adsorption capacity obtained using the Baudu isotherm were higher than those calculated by the Fritz-Schlunder-III model. However, according to Fritz-Schlunder-III and Baudu isotherm models, CBK 1/1 gave a higher maximum adsorption capacity than CBH 2/1 . Examining the exponent values in the Fritz-Schlunder-IV model, α approached unity while β did not, indicating that the adsorption data could not preferably be fitted with the Langmuir isotherm.
The adsorption data were analyzed according to the non-linear form of the Comparing the selected isotherm models used, we can say from the values of the determination coefficient and the values of error functions, that the models with 3, 4 and 5 parameters better described the adsorption process than the models with 2 parameters for the two ACs.

Effect of Temperature and Thermodynamic Study
The results of the effect of temperature on the adsorption of thymol blue are shown in Figure 7 below. The adsorption of thymol blue appears to be significantly affected by the temperature which influences the quantity adsorbed. We found that as the temperature increased from 35˚C to 65˚C, the adsorption capacity also increased from 15.08 mg•g −1 to 18.062 mg•g −1 for CBH 2/1 and 15.711 mg•g −1 to 18.659 mg•g −1 for CBK 1/1 . Hence, a change in temperature has an impact on the adsorption efficiency and the adsorption capacity of the adsorbent. Journal of Encapsulation and Adsorption Sciences Higher temperatures could also lead to stronger interactions during the formation of intermolecular hydrogen bonding between thymol blue and the two ACs, which is an important contribution to the adsorption process. Also, the rise in solution temperature increased the surface activity and kinetic energy of the thymol blue molecules which in turn increased mass transfer across the liquid film surrounding the different particles of the activated carbons.
The thermodynamic behavior for the adsorption on the two ACs was studied in more detail. The thermodynamic parameters such as the changes in Gibbs free energy (∆G°), enthalpy (∆H°), and entropy (∆S°) were calculated using the following equations: where, K D is the distribution coefficient for adsorption, Q e is the adsorption capacity at different temperature (mg•g −1 ) and C e is the equilibrium dye concentration in solution (mg•L −1 ), R is the gas constant (8.314 J•mol −1 •K −1 ), and T is absolute temperature (K). The Van't Hoff Equation (33) was plotted and the plot was used to calculate the change in enthalpy (∆H°) and entropy (∆S°) of adsorption for the adsorption of thymol blue. The results were used to determine the change in Gibbs free energy (32). These thermodynamic quantities are presented in Table 9 below. The negative values of ∆G° for the temperature range 308 K to 338 K indicate that the adsorption of thymol blue onto the ACs is a spontaneous process. The decrease of ∆G° values with the increase of temperature indicates that the adsorption became less favorable at higher temperatures [42]. The positive values of ∆H° show that the adsorption process is endothermic. The change in entropy (∆S°) was positive, indicating that the entropy of the system increased during the adsorption process. However, it should also be noted that the entropy of the universe (including the system and the surroundings) might increase because the adsorption process was not an isolated one [9]. A positive value of ∆S° reflects the affinity of the adsorbent towards the adsorbate species.

Conclusions
The ACs prepared from garcinia cola nut shell were characterized to obtain their physic-chemical properties. They were used to adsorb thymol blue from aqueous solution by batch process. Textural characteristics obtained by N 2 adsorptiondesorption analysis showed that the surfaces of the ACs were heterogeneous, and contained micropores and mesopores.
Batch adsorption studies showed that increasing initial dye concentration enhanced the interaction between thymol blue and ACs which resulted in the increase of adsorption capacity. Electrostatic attraction and π-π interactions were shown to dominate the adsorption mechanism. The quantity of thymol blue adsorbed increased with increasing temperature. According to seven different error functions and coefficient of determination, kinetic studies shows that, the adsorption data best fitted the Avrami non-linear equation and the pseudo-second order kinetic model thus indicating that the adsorption of thymol blue on the ACs may also be of a chemical nature. Isotherms models confirmed that the system is heterogeneous in nature and the adsorption is of a physical nature.
In future studies, we plan to use these activated carbons for the purification of biogas by adsorption of the H 2 S contained in the synthesized mixture. We also intend to use the activated carbons in studies of column adsorption.