# Pineapple Bark Performance in Dyes Adsorption: Optimization by the Central Composite Design

Authors: Ahlam Fegousse; Abdelali El Gaidoumi; Youssef Miyah; Rabea El Mountassir; Anissa Lahrichi
Journal: Journal of Chemistry
 (2019)

---

## Abstract

This work is concerned with the study of the adsorption in aqueous medium of a three-dye mixture which contains Methylene Blue, Brilliant Green, and Congo Red on the pineapple bark. This adsorbent material has been characterized by scanning electron microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR). The experimental design methodology, based on the response surface methodology (RSM) by the central composite design (CCD), has been applied for the optimization of the parameters, namely, the temperature, dose of the adsorbent, and pH. The yield reached 98.91% under optimal conditions (T = 30°C; adsorbent dose = 2.5 g·L−1; pH = 9.8) at an initial concentration of 20 mg·L−1.


---

## Body

## 1. Introduction

Water is an essential element for the survival of all living organisms. Nowadays, the water sources contaminated by a wide variety of pollutants coming from industrial effluents are a subject of several researchers [1, 2]. Dyes have a large amount of pollutants and can be identified even by the human eye [3]. In this respect, dyes are used in various industries such as textiles, food, paper, rubber, plastics, and cosmetics [4]. Hence, the discharge of wastewater from these industries to water resources poses unavoidable problems. The presence of dyes in trace water is undesirable because most of them are toxic, mutagenic, and carcinogenic [5]. Furthermore, dyes not only prevent the penetration of light and reduce the photosynthetic activities of rivers but also disturb the aquatic balance. Thus, the elimination of dyes from wastewater before discharge is an indispensable task. In this case, several techniques such as flocculation, adsorption, oxidation, electrolysis, biodegradation, ion exchange, and photocatalysis have been used for the removal of dyes in wastewater [6]. Among the several techniques, adsorption has received considerable attention because of its many advantages in terms of cost, ease of use, flexibility and simplicity of design, and insensitivity to toxic pollutants [7, 8]. Similarly, of the various adsorbents, activated carbon may be the most effective adsorbent for the removal of dyes due to its excellent adsorption capacity [9]. However, its use is limited on account of its high cost [10]. Thus, attention has shifted to finding cheaper and more effective alternatives. For this reason, natural materials, agricultural and industrial waste, and biosorbents represent potential alternatives. Additionally, a large number of unconventional and low-cost adsorbents have been proposed by many researchers for the removal of dyes [11, 12] including agricultural waste such as sawdust [13], bark [14], and orange peel [15] and industrial waste, namely, metal hydroxide sludge [16], red mud [17], ashes flying [18], clays [19–21], diatomites [22, 23], zeolites [24, 25], siliceous materials [26, 27], biosorbents [28–31], and others (cyclodextrin [32, 33], starch [34], cotton [35], etc.).The purpose of this work was to study the adsorption of a three-dye mixture that contained Methylene Blue, Brilliant Green, and Congo Red on the pineapple bark according to the surface response method using the central composite design that consists of optimizing physicochemical parameters such as pH, temperature, and mass of the adsorbent.


## 2. Materials and Methods

### 2.1. Preparation of the Adsorbate and Adsorbent

The studied dyes have been supplied by Sigma-Aldrich and used without prior treatment. Congo Red (CR) is an anionic dye, part of the class of azo, with chemical formula C32H22N6O6S2Na2 and molar mass equal to 696.66 g·mol−1. Its molecular structure is illustrated in Figure 1(a). Methylene Blue (MB) and Brilliant Green (BG) which are cationic dyes, with chemical formulas of C16H18ClN3S and C27H34N2O4S and molar masses equal to 319.85 g·mol−1 and 482.63 g·mol−1, respectively. Their molecular structures are shown schematically in Figures 1(b) and 1(c). The colored solution (DM) has been precisely prepared by mixing the three dyes with same proportion in ultrapure water.Figure 1
Chemical structure of dyes: (a) Congo Red, (b) Methylene Blue, and (c) Brilliant Green.The pineapple bark (PB) was washed with ultrapure water followed by free-air drying. After that, it was crushed and sieved to a fraction less than 125 nm. The resulting powder was then washed with ultrapure water to remove impurities and then oven-dried at a temperature of 60°C.


### 2.2. Adsorption Procedure

The experiments have been performed in a batch system with constant stirring in Pyrex beakers using a volume of 50 mL of colored solution. Moreover, 4 mL samples have been taken using a 0.45μm diameter syringe filter and then analyzed by a Jasco V-750 UV-Vis spectrophotometer at a wavelength λmax = 624 nm.The capacity of adsorption was calculated according to(1)qt=C0−Ct∗Vm,where qt is the amount of adsorbed dye per unit adsorbent at instant t (mg/g), m is the used mass of adsorbent (g), C0 is the initial concentration of dye in solution (mg/L), Ct is the concentration of dye in solution at instant t (mg/L), and V is the used volume of solution (L).The concentrations and removal percentages were calculated taking into account the dyes’ purity: ≥35%, ≥82%, and ≥90% for CR, MB, and BG, respectively.


### 2.3. Characterization Techniques

Morphological analysis was observed by scanning electron microscopy (SEM) using a Quanta 200 FEI microscope. Fourier transform infrared spectroscopy (FTIR) analysis was performed in the 400–4000 cm−1 interval using a Bruker Vertex 70 spectrometer to determine the functional analysis of the PB powder.The point zero of charge (pHPZC) was determined by a classical method [36], which consists in preparing 50 mL solutions of NaCl (0.01 M) and adjusting their pH to precise values ranging from 2 to 12 by addition of NaOH or HCl (0.1 M). Next, 0.5 g of adsorbent was introduced into each solution. The suspensions were stirred at room temperature for 24 h before determining the pHfinal. The pHPZC was calculated based on the curve pHfinal − pHinitial = f(pHinitial).


## 2.1. Preparation of the Adsorbate and Adsorbent

The studied dyes have been supplied by Sigma-Aldrich and used without prior treatment. Congo Red (CR) is an anionic dye, part of the class of azo, with chemical formula C32H22N6O6S2Na2 and molar mass equal to 696.66 g·mol−1. Its molecular structure is illustrated in Figure 1(a). Methylene Blue (MB) and Brilliant Green (BG) which are cationic dyes, with chemical formulas of C16H18ClN3S and C27H34N2O4S and molar masses equal to 319.85 g·mol−1 and 482.63 g·mol−1, respectively. Their molecular structures are shown schematically in Figures 1(b) and 1(c). The colored solution (DM) has been precisely prepared by mixing the three dyes with same proportion in ultrapure water.Figure 1
Chemical structure of dyes: (a) Congo Red, (b) Methylene Blue, and (c) Brilliant Green.The pineapple bark (PB) was washed with ultrapure water followed by free-air drying. After that, it was crushed and sieved to a fraction less than 125 nm. The resulting powder was then washed with ultrapure water to remove impurities and then oven-dried at a temperature of 60°C.


## 2.2. Adsorption Procedure

The experiments have been performed in a batch system with constant stirring in Pyrex beakers using a volume of 50 mL of colored solution. Moreover, 4 mL samples have been taken using a 0.45μm diameter syringe filter and then analyzed by a Jasco V-750 UV-Vis spectrophotometer at a wavelength λmax = 624 nm.The capacity of adsorption was calculated according to(1)qt=C0−Ct∗Vm,where qt is the amount of adsorbed dye per unit adsorbent at instant t (mg/g), m is the used mass of adsorbent (g), C0 is the initial concentration of dye in solution (mg/L), Ct is the concentration of dye in solution at instant t (mg/L), and V is the used volume of solution (L).The concentrations and removal percentages were calculated taking into account the dyes’ purity: ≥35%, ≥82%, and ≥90% for CR, MB, and BG, respectively.


## 2.3. Characterization Techniques

Morphological analysis was observed by scanning electron microscopy (SEM) using a Quanta 200 FEI microscope. Fourier transform infrared spectroscopy (FTIR) analysis was performed in the 400–4000 cm−1 interval using a Bruker Vertex 70 spectrometer to determine the functional analysis of the PB powder.The point zero of charge (pHPZC) was determined by a classical method [36], which consists in preparing 50 mL solutions of NaCl (0.01 M) and adjusting their pH to precise values ranging from 2 to 12 by addition of NaOH or HCl (0.1 M). Next, 0.5 g of adsorbent was introduced into each solution. The suspensions were stirred at room temperature for 24 h before determining the pHfinal. The pHPZC was calculated based on the curve pHfinal − pHinitial = f(pHinitial).


## 3. Results and Discussions

### 3.1. Characterizations of the Adsorbent

The use of the SEM technique makes it possible to visualize the morphology of the surface of the adsorbent. The SEM micrographs of PB are shown in Figure2. The obtained images display not only the amorphous nature and the heterogeneous morphology of PB but also indicate the presence of both a grainy structure (Figure 2(b)) close to the appearance of the grapes and a filamentous structure and cellulosic fibers of various sizes and shapes (Figure 2(c)), which promote the fixation of dyes. Consequently, these results are similar to those obtained by Miyah et al. [37].Figure 2
SEM micrographs of PB: magnification (a) ×2000; (b, c) ×4000 (reproduced from Fegousse et al. [38]).FTIR spectrum analysis (Figure3) shows a band around 3444.03 cm−1 which is characteristic of the presence of a band O-H group, a small band around 2925.38 cm−1 which attributes to C-H elongation vibrations, a band at 2361.24 cm−1 which could be due to the nitrile triple bond C≡N, a band at 1633.21 cm−1 which originates from the vibration of the C=O groups, a band at 1398.56 cm−1 which coincides with the OH band aromatic alcohol, a 1255.45 cm−1 band which corresponds to the CO bond, a 1037.19 cm−1 band which is assigned to the primary alcohol function (C-OH), and a band at 668.09 cm−1 which characterizes the function ≡CH [39].Figure 3
FTIR spectra of PB (reproduced from Fegousse et al. [38]).The pHPZC corresponds to the pH value for which the net electric charge of the surface of the material is neutral [40]. pHPZC is very important in adsorption phenomena, especially when electrostatic forces are involved in the mechanisms. Figure 4 shows that the pHPZC of PB is 6.7. This point illustrates an amphoteric behavior of PB. In fact, the surface of PB is positively charged at pH < pHPZC and negatively charged at pH > pHPZC. As pH decreases, the number of negatively charged sites decreases, and the number of positively charged sites increases and vice versa when pH increases [41–43].Figure 4
The point zero of charge of PB.


### 3.2. Effect of Contact Time

As illustrated in Figure5, the equilibrium time is reached in the first five minutes for a concentration range varying between 20 mg·L−1 and 50 mg·L−1. This can be explained by the fact that the available number of active sites on the surface of adsorbent becomes saturated after five minutes [44]. In this case, the adsorption capacity (Qe) remains constant after equilibrium time, which denotes that the adsorption is rapid and confirms strong chemical interactions between the dye molecules and adsorbent. In this respect, these outcomes allow us to have a quick overview on the effectiveness of PB.Figure 5
Effect of contact time on (a) adsorption efficiency and (b) adsorption capacity at different concentrations.m = 2.g·L−1; [DM] = 20–50 mg·L−1; T = 20°C, pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)


### 3.3. Effect of Adsorbent Dose

Figure6(a) shows the monitoring of the adsorption capacity as a function of time at an initial concentration of 20 mg·L−1; a large decrease in the DM adsorbed amount with the increase of the adsorbent dose was observed. This phenomenon is due to the decrease of number of DM molecules (decrease in concentration) and the difficulty of adsorbent-adsorbate interaction, which explains the saturation of the active sites of the adsorbent. Additionally, while the curve which is a function of the adsorption efficiency (Figure 6(b)) has a proportionally inverse appearance to that of the adsorbed quantity, it is obviously noticed that the adsorption efficiency increases with the increase in the adsorbent dose. Hence, 2 g·L−1 of the adsorbent constitutes an optimal dose to remove 97.93% of DM which corresponds to an adsorption capacity of 14.6 mg·g−1.Figure 6
Effect of PB dose on (a) adsorption capacity and (b) adsorption efficiency.m = 0.5–3 g·L−1; [DM] = 20 mg·L−1; T = 20°C; pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)


### 3.4. Effect of Temperature

The outcomes of the evolution of the adsorption capacity as a function of temperature at an initial concentration of 20 mg·L−1 are displayed in Figure 7. In this respect, the analysis of the curve shows that the temperature has a negative effect on the adsorption capacity; it affects the chemical potential of the material: the mobility of the system increases but the interactions decrease. This suggests that the interaction of PB and DM may be exothermic. Thus, the adsorption capacity reaches 9.87 mg·g−1 for a temperature of 20°C which corresponds to a removal efficiency of 98.97%.Figure 7
Effect of temperature on (a) adsorption efficiency and (b) adsorption capacity.m = 2 g·L−1; [DM] = 20 mg·L−1; T = 20–50°C; pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)The adsorption may result in decreasing or increasing the surface energy of the adsorbent materials. The parameters that can describe the adsorption thermodynamics due to the transfer of the solute to the solid-liquid interface are the standard enthalpy (ΔH°), the standard entropy (ΔS°), and the free standard enthalpy (ΔG°). The values of ΔH° and ΔS° were calculated using equation (2) [4] by plotting ln Kd as a function of 1/T (Figure 8).(2)lnKd=ΔS°R−ΔH°RT,where R is the ideal gas constant (R = 8.314 J·mol−1·K−1); T is the absolute temperature of solution (K); and Kd is the distribution coefficient equal to Qe/Ce (Qe and Ce represent the capacity of adsorption and concentration at equilibrium, respectively).Figure 8
lnKd versus 1/T for the adsorption of DM on PB.The values of ΔG° were calculated using the following formula [4]:(3)ΔG°=−RTlnKd.The obtained results show that the standard enthalpy value is negative (Table1), which confirms that the adsorption of the dye molecules at the sites of the PB adsorbent is exothermic. The negative values of the free standard enthalpy and standard entropy indicate that the process of adsorption of the DM on PB is spontaneous and follows a disorder, respectively.Table 1
Thermodynamic parameters of the adsorption process of DM on PB.
Temperature (°C)
ΔH° (KJ·mol−1)
ΔS° (J·mol−1·K−1)
ΔG° (J·mol−1)
R
2
20
−22.44
−51.68
−37591.20
0.9993
30
−38108.04
40
−38624.88
50
−39141.72


### 3.5. Effect of pH

The evolution of DM adsorption as a function of pH at an initial concentration of 20 mg·L−1and dose of 2 g·L−1 is illustrated in Figure 9. The adsorption is favorable at acid pH due to the positive charge of the adsorbent which promotes the adsorption of the negatively charged anionic dye (RC). On the contrary, at basic pH, the adsorption is stable even in the presence of cationic dyes (MB and BG). This indicates that the adsorbent is selective for the anionic dye and only at acidic pH. Furthermore, the adsorption of the cationic dyes is owing to the activation of the surface of the adsorbent by the anionic dye which makes the surface negative and, therefore, active in front of the cationic dyes. Indeed, for acid pH values, the cationic dyes are in competition with the H+ ions presented in high concentration in the solution [45,46]. These ions, positively charged, are more adsorbed than cationic dyes and then beneficial for adsorption of anionic dye with negative charge. Additionally, when the pH increases, there is a decrease in the H+ ions and increase in OH−, which explains the favorability of the adsorbent to absorb the cationic dyes. Based on the obtained results, pH = 6 is the optimum value for the adsorption of DM.Figure 9
pH effect on PB adsorption capacity.m = 2 g·L−1; [DM] = 20 mg·L−1; T = 20°C; pH = 2–12; Vag = 400 rpm; t = 30 min.


### 3.6. Optimization of Adsorption

The aim of this study is to optimize the adsorption of the mixture of dyes by the response surface methodology (RSM). The design used for the optimization, belonging to RSM, is the central composite design (CCD); its equation is polynomial of second order [47]:(4)Yi=b0+b1X1+b2X2+b3X3+b12X1X2+b13X1X3+b23X2X3+b11X12+b22X22+b33X32,where Y represents the percentage of adsorption, Xi represents the independent variable (temperature, pH or dose of the adsorbent), XiXj demonstrates the interactions between the variables, and bi represents the coefficient of the effect resulted by Xi.The number of experiments to be performed is determined by the following formula:(5)N=2K+2K+f,where K is the number of independent variables and f is the number of experiments performed at the center of the experimental scope.To realize 17 experiments according to the computation of the formula (3), each variable is represented by 5 levels: −1.68; −1; 0; +1; 1.68.


### 3.7. Table of Variables

Table2 represents the three studied variables and their corresponding levels. The reaction time is 5 min, and the concentration of DM was fixed at 20 mg·L−1.Table 2
Description of the variables and their levels.
Parameters
Variables
Levels
−1.68
−1
0
+1
+1.68
T (°C)
X
1
20
30
40
50
60
pH
X
2
2
4
6
8
10
m (adsorbent dose) (g·L−1)
X
3
0.5
1
1.5
2
2.5


### 3.8. Experimental Matrix and Postulated Mathematical Model

Table3 summarizes the 17 performed experiments according to the CCD design. Experiments 15, 16, and 17 represent the experiments at the center. According to the results, the percentage of adsorption varies between 75.55% and 97.18%. The mathematical model generated by the Design-Expert software version 11 is(6)Y=97.22+1.83X1+7.22X2+1.77X3−2.95X12−3.35X22−1.43X32+0.52X1X2+1.77X1X3−3.56X2X3.Table 3
Experimental matrix.
Experience
Variables
Response
X
1
X
2
X
3
Y
1
−1
−1
−1
77.01
2
1
−1
−1
76.81
3
−1
1
−1
98.75
4
1
1
−1
98.45
5
−1
−1
1
85.42
6
1
−1
1
90.13
7
−1
1
1
90.75
8
1
1
1
99.70
9
−1.68
0
0
85.15
10
1.68
0
0
92.18
11
0
−1.68
0
75.55
12
0
1.68
0
99.52
13
0
0
−1.68
90.25
14
0
0
1.68
95.70
15
0
0
0
97.30
16
0
0
0
97.18
17
0
0
0
97.25


### 3.9. ANOVA Analysis

Analysis of variance (ANOVA) has been used, and the results are given in Table4. Generally, the adaptation of the postulated model is confirmed by a high Fisher value (F) with as low a probability (P) as possible. From the ANOVA test, the obtained P value of 0.0001 is less than 0.05 value evaluated by the Fischer relation which confirms that the model terms are significant [48,49]. The resulting F value of 230.5 is greater than the Fisher value of F9.7 = 3.8 for a 95% confidence level. So, the model is well suited to experimental data [50–52].Table 4
Results of ANOVA for the fitted model.
Source
Degree of freedom
Sum of squares
Mean square
F value
P value
Model
9
1103.07
122.56
230.5
<0.0001
Residual
7
3.72
0.53
Total
16
1106.7916The fit quality of the model can be justified by the correlation coefficient (R2) which indicates the proportion of the response variability taken into account by the model [53]. The best correlation is obtained for which R2 = 0.80 is suggested [54]. The obtained R value of 0.9966 (Table 5) implied that a large proportion of the variability of the response is considered by the model; 0.34% of the variations are not taken into account.Table 5
Correlation of model.
R
2
R
adj
2
Standard deviation
Mean
Coefficient of variation (%)
0.9966
0.9923
0.73
91.01
0.8


### 3.10. Estimation of Coefficients

The determination of the coefficients makes it possible to judge the significance of the effects of the variables and their interactions. Thus, Table6 summarizes the results found.Table 6
Estimated regression coefficients and their significance.
Coefficients
Estimated values
Standard deviation
F value
P value
b
0
97.22
0.42
230.5
<0.0001
b
1
1.83
0.2
85.95
<0.0001
b
2
7.22
0.2
1338.61
<0.0001
b
3
1.77
0.2
80.29
<0.0004
b
11
−2.95
0.22
184.92
0.0847
b
22
−3.35
0.22
238.34
0.0002
b
33
−1.43
0.22
43.32
<0.0001
b
12
0.52
0.26
4.03
<0.0001
b
13
1.77
0.26
47.14
<0.0001
b
23
−3.56
0.26
190.68
<0.0001Generally, interaction coefficients that have aP value less than 0.05 consider themselves as coefficients having a significant effect on the response. According to the previous table, all P values are less than 0.05; so, all coefficients are considered to have a very crucial effect on the response, except in the case of b11 for which P is equal to 0.0847.


### 3.11. Graphic Study

The graphical representations of three response surfaces are illustrated in Figure10. In this respect, the optimum point of the adsorption yield has been achieved based on the response surface of the X2X3, interaction between pH and adsorbent dose (Figure 10(c)), which have already shown a positive effect on yield (Figures 6 and 9). Besides, the maximum yield that is not acquired according to the X1X2 and X1X3 interactions (Figures 10(a) and 10(b)) is justified by the negative effect of temperature (X1) on the adsorption process (Figure 7). Hence, the optimal conditions are summarized in Table 7.Figure 10
Response surfaces showing adsorption efficiency as function of (a) temperature and pH, (b) temperature and adsorbent dose, and (c) pH and adsorbent dose.
(a)
(b)
(c)Table 7
Established optimal conditions.
Optimal conditions
Adsorption (%)
Temperature (°C)
pH
Adsorbent dose (g·L−1)
30
9.8
2.5
99.96


### 3.12. Evolution of Adsorption under Optimal Conditions

To confirm the validity of the theoretical values of the optimization, a follow-up of the evolution of the adsorption under the optimal conditions has been carried out (Figure11). In this case, we notice that the system leads to a rapid adsorption of 98.91% after 5 min close to the theoretical value.Figure 11
Evolution of the adsorption of DM under optimal conditions.From what have been stated the above, CCD has shown great satisfaction in optimizing the parameters influencing the adsorption of the studied dyes mixture since the statistical tests are well adapted to the experimental results. These outcomes are in total agreement with the various publications, which show the validation of the CCD plane in the study of adsorption of dyes on different types of adsorbents [55–57]. In addition, the yield obtained at optimal conditions is very high, which confirms the performance of the pineapple bark as an adsorbent for the studied mixture.


## 3.1. Characterizations of the Adsorbent

The use of the SEM technique makes it possible to visualize the morphology of the surface of the adsorbent. The SEM micrographs of PB are shown in Figure2. The obtained images display not only the amorphous nature and the heterogeneous morphology of PB but also indicate the presence of both a grainy structure (Figure 2(b)) close to the appearance of the grapes and a filamentous structure and cellulosic fibers of various sizes and shapes (Figure 2(c)), which promote the fixation of dyes. Consequently, these results are similar to those obtained by Miyah et al. [37].Figure 2
SEM micrographs of PB: magnification (a) ×2000; (b, c) ×4000 (reproduced from Fegousse et al. [38]).FTIR spectrum analysis (Figure3) shows a band around 3444.03 cm−1 which is characteristic of the presence of a band O-H group, a small band around 2925.38 cm−1 which attributes to C-H elongation vibrations, a band at 2361.24 cm−1 which could be due to the nitrile triple bond C≡N, a band at 1633.21 cm−1 which originates from the vibration of the C=O groups, a band at 1398.56 cm−1 which coincides with the OH band aromatic alcohol, a 1255.45 cm−1 band which corresponds to the CO bond, a 1037.19 cm−1 band which is assigned to the primary alcohol function (C-OH), and a band at 668.09 cm−1 which characterizes the function ≡CH [39].Figure 3
FTIR spectra of PB (reproduced from Fegousse et al. [38]).The pHPZC corresponds to the pH value for which the net electric charge of the surface of the material is neutral [40]. pHPZC is very important in adsorption phenomena, especially when electrostatic forces are involved in the mechanisms. Figure 4 shows that the pHPZC of PB is 6.7. This point illustrates an amphoteric behavior of PB. In fact, the surface of PB is positively charged at pH < pHPZC and negatively charged at pH > pHPZC. As pH decreases, the number of negatively charged sites decreases, and the number of positively charged sites increases and vice versa when pH increases [41–43].Figure 4
The point zero of charge of PB.


## 3.2. Effect of Contact Time

As illustrated in Figure5, the equilibrium time is reached in the first five minutes for a concentration range varying between 20 mg·L−1 and 50 mg·L−1. This can be explained by the fact that the available number of active sites on the surface of adsorbent becomes saturated after five minutes [44]. In this case, the adsorption capacity (Qe) remains constant after equilibrium time, which denotes that the adsorption is rapid and confirms strong chemical interactions between the dye molecules and adsorbent. In this respect, these outcomes allow us to have a quick overview on the effectiveness of PB.Figure 5
Effect of contact time on (a) adsorption efficiency and (b) adsorption capacity at different concentrations.m = 2.g·L−1; [DM] = 20–50 mg·L−1; T = 20°C, pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)


## 3.3. Effect of Adsorbent Dose

Figure6(a) shows the monitoring of the adsorption capacity as a function of time at an initial concentration of 20 mg·L−1; a large decrease in the DM adsorbed amount with the increase of the adsorbent dose was observed. This phenomenon is due to the decrease of number of DM molecules (decrease in concentration) and the difficulty of adsorbent-adsorbate interaction, which explains the saturation of the active sites of the adsorbent. Additionally, while the curve which is a function of the adsorption efficiency (Figure 6(b)) has a proportionally inverse appearance to that of the adsorbed quantity, it is obviously noticed that the adsorption efficiency increases with the increase in the adsorbent dose. Hence, 2 g·L−1 of the adsorbent constitutes an optimal dose to remove 97.93% of DM which corresponds to an adsorption capacity of 14.6 mg·g−1.Figure 6
Effect of PB dose on (a) adsorption capacity and (b) adsorption efficiency.m = 0.5–3 g·L−1; [DM] = 20 mg·L−1; T = 20°C; pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)


## 3.4. Effect of Temperature

The outcomes of the evolution of the adsorption capacity as a function of temperature at an initial concentration of 20 mg·L−1 are displayed in Figure 7. In this respect, the analysis of the curve shows that the temperature has a negative effect on the adsorption capacity; it affects the chemical potential of the material: the mobility of the system increases but the interactions decrease. This suggests that the interaction of PB and DM may be exothermic. Thus, the adsorption capacity reaches 9.87 mg·g−1 for a temperature of 20°C which corresponds to a removal efficiency of 98.97%.Figure 7
Effect of temperature on (a) adsorption efficiency and (b) adsorption capacity.m = 2 g·L−1; [DM] = 20 mg·L−1; T = 20–50°C; pH = 6; Vag = 400 rpm; t = 60 min.
(a)
(b)The adsorption may result in decreasing or increasing the surface energy of the adsorbent materials. The parameters that can describe the adsorption thermodynamics due to the transfer of the solute to the solid-liquid interface are the standard enthalpy (ΔH°), the standard entropy (ΔS°), and the free standard enthalpy (ΔG°). The values of ΔH° and ΔS° were calculated using equation (2) [4] by plotting ln Kd as a function of 1/T (Figure 8).(2)lnKd=ΔS°R−ΔH°RT,where R is the ideal gas constant (R = 8.314 J·mol−1·K−1); T is the absolute temperature of solution (K); and Kd is the distribution coefficient equal to Qe/Ce (Qe and Ce represent the capacity of adsorption and concentration at equilibrium, respectively).Figure 8
lnKd versus 1/T for the adsorption of DM on PB.The values of ΔG° were calculated using the following formula [4]:(3)ΔG°=−RTlnKd.The obtained results show that the standard enthalpy value is negative (Table1), which confirms that the adsorption of the dye molecules at the sites of the PB adsorbent is exothermic. The negative values of the free standard enthalpy and standard entropy indicate that the process of adsorption of the DM on PB is spontaneous and follows a disorder, respectively.Table 1
Thermodynamic parameters of the adsorption process of DM on PB.
Temperature (°C)
ΔH° (KJ·mol−1)
ΔS° (J·mol−1·K−1)
ΔG° (J·mol−1)
R
2
20
−22.44
−51.68
−37591.20
0.9993
30
−38108.04
40
−38624.88
50
−39141.72


## 3.5. Effect of pH

The evolution of DM adsorption as a function of pH at an initial concentration of 20 mg·L−1and dose of 2 g·L−1 is illustrated in Figure 9. The adsorption is favorable at acid pH due to the positive charge of the adsorbent which promotes the adsorption of the negatively charged anionic dye (RC). On the contrary, at basic pH, the adsorption is stable even in the presence of cationic dyes (MB and BG). This indicates that the adsorbent is selective for the anionic dye and only at acidic pH. Furthermore, the adsorption of the cationic dyes is owing to the activation of the surface of the adsorbent by the anionic dye which makes the surface negative and, therefore, active in front of the cationic dyes. Indeed, for acid pH values, the cationic dyes are in competition with the H+ ions presented in high concentration in the solution [45,46]. These ions, positively charged, are more adsorbed than cationic dyes and then beneficial for adsorption of anionic dye with negative charge. Additionally, when the pH increases, there is a decrease in the H+ ions and increase in OH−, which explains the favorability of the adsorbent to absorb the cationic dyes. Based on the obtained results, pH = 6 is the optimum value for the adsorption of DM.Figure 9
pH effect on PB adsorption capacity.m = 2 g·L−1; [DM] = 20 mg·L−1; T = 20°C; pH = 2–12; Vag = 400 rpm; t = 30 min.


## 3.6. Optimization of Adsorption

The aim of this study is to optimize the adsorption of the mixture of dyes by the response surface methodology (RSM). The design used for the optimization, belonging to RSM, is the central composite design (CCD); its equation is polynomial of second order [47]:(4)Yi=b0+b1X1+b2X2+b3X3+b12X1X2+b13X1X3+b23X2X3+b11X12+b22X22+b33X32,where Y represents the percentage of adsorption, Xi represents the independent variable (temperature, pH or dose of the adsorbent), XiXj demonstrates the interactions between the variables, and bi represents the coefficient of the effect resulted by Xi.The number of experiments to be performed is determined by the following formula:(5)N=2K+2K+f,where K is the number of independent variables and f is the number of experiments performed at the center of the experimental scope.To realize 17 experiments according to the computation of the formula (3), each variable is represented by 5 levels: −1.68; −1; 0; +1; 1.68.


## 3.7. Table of Variables

Table2 represents the three studied variables and their corresponding levels. The reaction time is 5 min, and the concentration of DM was fixed at 20 mg·L−1.Table 2
Description of the variables and their levels.
Parameters
Variables
Levels
−1.68
−1
0
+1
+1.68
T (°C)
X
1
20
30
40
50
60
pH
X
2
2
4
6
8
10
m (adsorbent dose) (g·L−1)
X
3
0.5
1
1.5
2
2.5


## 3.8. Experimental Matrix and Postulated Mathematical Model

Table3 summarizes the 17 performed experiments according to the CCD design. Experiments 15, 16, and 17 represent the experiments at the center. According to the results, the percentage of adsorption varies between 75.55% and 97.18%. The mathematical model generated by the Design-Expert software version 11 is(6)Y=97.22+1.83X1+7.22X2+1.77X3−2.95X12−3.35X22−1.43X32+0.52X1X2+1.77X1X3−3.56X2X3.Table 3
Experimental matrix.
Experience
Variables
Response
X
1
X
2
X
3
Y
1
−1
−1
−1
77.01
2
1
−1
−1
76.81
3
−1
1
−1
98.75
4
1
1
−1
98.45
5
−1
−1
1
85.42
6
1
−1
1
90.13
7
−1
1
1
90.75
8
1
1
1
99.70
9
−1.68
0
0
85.15
10
1.68
0
0
92.18
11
0
−1.68
0
75.55
12
0
1.68
0
99.52
13
0
0
−1.68
90.25
14
0
0
1.68
95.70
15
0
0
0
97.30
16
0
0
0
97.18
17
0
0
0
97.25


## 3.9. ANOVA Analysis

Analysis of variance (ANOVA) has been used, and the results are given in Table4. Generally, the adaptation of the postulated model is confirmed by a high Fisher value (F) with as low a probability (P) as possible. From the ANOVA test, the obtained P value of 0.0001 is less than 0.05 value evaluated by the Fischer relation which confirms that the model terms are significant [48,49]. The resulting F value of 230.5 is greater than the Fisher value of F9.7 = 3.8 for a 95% confidence level. So, the model is well suited to experimental data [50–52].Table 4
Results of ANOVA for the fitted model.
Source
Degree of freedom
Sum of squares
Mean square
F value
P value
Model
9
1103.07
122.56
230.5
<0.0001
Residual
7
3.72
0.53
Total
16
1106.7916The fit quality of the model can be justified by the correlation coefficient (R2) which indicates the proportion of the response variability taken into account by the model [53]. The best correlation is obtained for which R2 = 0.80 is suggested [54]. The obtained R value of 0.9966 (Table 5) implied that a large proportion of the variability of the response is considered by the model; 0.34% of the variations are not taken into account.Table 5
Correlation of model.
R
2
R
adj
2
Standard deviation
Mean
Coefficient of variation (%)
0.9966
0.9923
0.73
91.01
0.8


## 3.10. Estimation of Coefficients

The determination of the coefficients makes it possible to judge the significance of the effects of the variables and their interactions. Thus, Table6 summarizes the results found.Table 6
Estimated regression coefficients and their significance.
Coefficients
Estimated values
Standard deviation
F value
P value
b
0
97.22
0.42
230.5
<0.0001
b
1
1.83
0.2
85.95
<0.0001
b
2
7.22
0.2
1338.61
<0.0001
b
3
1.77
0.2
80.29
<0.0004
b
11
−2.95
0.22
184.92
0.0847
b
22
−3.35
0.22
238.34
0.0002
b
33
−1.43
0.22
43.32
<0.0001
b
12
0.52
0.26
4.03
<0.0001
b
13
1.77
0.26
47.14
<0.0001
b
23
−3.56
0.26
190.68
<0.0001Generally, interaction coefficients that have aP value less than 0.05 consider themselves as coefficients having a significant effect on the response. According to the previous table, all P values are less than 0.05; so, all coefficients are considered to have a very crucial effect on the response, except in the case of b11 for which P is equal to 0.0847.


## 3.11. Graphic Study

The graphical representations of three response surfaces are illustrated in Figure10. In this respect, the optimum point of the adsorption yield has been achieved based on the response surface of the X2X3, interaction between pH and adsorbent dose (Figure 10(c)), which have already shown a positive effect on yield (Figures 6 and 9). Besides, the maximum yield that is not acquired according to the X1X2 and X1X3 interactions (Figures 10(a) and 10(b)) is justified by the negative effect of temperature (X1) on the adsorption process (Figure 7). Hence, the optimal conditions are summarized in Table 7.Figure 10
Response surfaces showing adsorption efficiency as function of (a) temperature and pH, (b) temperature and adsorbent dose, and (c) pH and adsorbent dose.
(a)
(b)
(c)Table 7
Established optimal conditions.
Optimal conditions
Adsorption (%)
Temperature (°C)
pH
Adsorbent dose (g·L−1)
30
9.8
2.5
99.96


## 3.12. Evolution of Adsorption under Optimal Conditions

To confirm the validity of the theoretical values of the optimization, a follow-up of the evolution of the adsorption under the optimal conditions has been carried out (Figure11). In this case, we notice that the system leads to a rapid adsorption of 98.91% after 5 min close to the theoretical value.Figure 11
Evolution of the adsorption of DM under optimal conditions.From what have been stated the above, CCD has shown great satisfaction in optimizing the parameters influencing the adsorption of the studied dyes mixture since the statistical tests are well adapted to the experimental results. These outcomes are in total agreement with the various publications, which show the validation of the CCD plane in the study of adsorption of dyes on different types of adsorbents [55–57]. In addition, the yield obtained at optimal conditions is very high, which confirms the performance of the pineapple bark as an adsorbent for the studied mixture.


## 4. Conclusions

The performance of the pineapple bark in the adsorption of the dyes mixture Methylene Blue, Brilliant Green, and Congo Red has been studied in a batch system. In this respect, the parameters the mass of the adsorbent, the concentration of the dyes mixture, the pH, and the temperature have indeed influenced the adsorption process. Besides, adsorption at a concentration of 20 mg·L−1 is rapid during the first five minutes of contact for a mass of 2 g·L−1. Moreover, the adsorption is favorable in acidic medium and is of the exothermic type. Furthermore, through the experimental design methodology, CCD plan, we have optimized the physicochemical parameters such as the mass of the adsorbent, the pH, and the temperature influencing the adsorption of the dyes mixture into pineapple bark. During 5 min of the adsorption, under the optimum conditions pH = 9.8, T = 30°C, and adsorbent dose = 2.5 g·L−1, the output reaches 98.91% at an initial concentration of 20 mg·L−1 and for a volume of 50 mL.



---

*Source: 3017163-2019-02-14.xml*