Modeling and Optimizing the Crystal Violet Dye Adsorption on Kaolinite Mixed with Cellulose Waste Red Bean Peels: Insights into the Kinetic, Isothermal, Thermodynamic, and Mechanistic Study

In this study, a new eco-friendly kaolinite–cellulose (Kaol/Cel) composite was prepared from waste red bean peels (Phaseolus vulgaris) as a source of cellulose to serve as a promising and effective adsorbent for the removal of crystal violet (CV) dye from aqueous solutions. Its characteristics were investigated through the use of X-ray diffraction, Fourier-transform infrared spectroscopy, scanning electron microscopy, energy-dispersive X-ray spectroscopy, and zero-point of charge (pHpzc). The Box–Behnken design was used to improve CV adsorption on the composite by testing its primary affecting factors: loading Cel into the composite matrix of Kaol (A: 0–50%), adsorbent dosage (B: 0.02–0.05 g), pH (C: 4–10), temperature (D: 30–60 °C), and duration (E: 5–60 min). The significant interactions with the greatest CV elimination efficiency (99.86%) are as follows: BC (adsorbent dose vs. pH) and BD (adsorbent dose vs. temperature) at optimum parameters (A: 25%, B: 0.05 g, C: 10, D: 45 °C, and E: 17.5 min) for which the CV’s best adsorption capacity (294.12 mg/g) was recorded. The Freundlich and pseudo-second-order kinetic models were the best isotherm and kinetic models fitting our results. Furthermore, the study investigated the mechanisms responsible for eliminating CV by utilizing Kaol/Cel–25. It detected multiple types of associations, including electrostatic, n-π, dipole–dipole, hydrogen bonding interactions, and Yoshida hydrogen bonding. These findings suggest that Kaol/Cel could be a promising starting material for developing a highly efficient adsorbent that can remove cationic dyes from aqueous environments.


Introduction
Today, the planet faces serious pollution issues, a global concern for the natural ecosystem. Much dye-containing wastewater relating to human activities is discharged without treatment, seriously affecting the environment and drinking water quality [1]. As contaminants in North African nations, including Algeria, dyes are dumped into the clay [25], kaolinite/biopolymer composites [26], cellulose@organically modified montmorillonite [27], cellulose/kaolinite-zeolite composite [28], are employed to remove different pollutants. They have shown promising efficiency for sequestering dyes. Because of the benefits offered by the materials described in the previous paragraph (Kaol, Cel), we decided to synthesize Kaol-Cel biopolymers as a recoverable adsorbent for removing CV dye from aquatic environments. This is the first cellulose-based compound extracted from red bean peels (Phaseolus vulgaris) (RBPs), which gives this study its innovative quality. To achieve maximum CV dye removal by Kaol/Cel, the experimental critical parameters needed to be improved and predicted using response surface techniques, included in Box-Behnken design (BBD). The investigation involved determining the synthesis conditions, including the loading of Cel into the Kaol biopolymeric matrix, as well as other crucial factors related to adsorption (such as adsorbent injection, pH, temperature, and residence period). Moreover, the study explored the kinetics, isotherms, and thermodynamic functions associated with the technique.

Material, Methods, and Instruments
RBPs were acquired from a local farm situated in El-Oued, Southeast Algeria. The samples underwent a thorough cleaning process to eliminate any dust or debris that may have been present on their surfaces. They were subsequently dried at 70 • C for 7-8 h and finely ground using a mixer grinder to ensure the samples were impurities-free. Raw kaolinite was obtained from Aougrout, Adrar, Southwest Algeria. The dye used in the study was CV (C 25 N 3 H 30 Cl, λ max = 590 nm, molecular weight: 407.979 g/mol), purchased from Sigma-Aldrich. Chemical products, such as hydrogen peroxide (H 2 O 2 ), hydrochloric acid (HCl), sodium hydroxide (NaOH), acetic acid (CH 3 COOH), sodium carbonate (Na 2 CO 3 ), sodium hexametaphosphate (NaPO 3 ) 6 , thiourea (SC(NH 2 ) 2 ), and sodium hypochlorite (NaOCl), were furnished by Merck (Darmstadt, Germany). The chemicals employed in this research were pure and of analytical quality. The characterization of Kaol-Cel-25 was carried out using Fourier transform-infrared (FT-IR) analysis (Cary 600 series FT-IR spectrophotometer, Agilent technologies, Santa Clara, CA, USA [29]), X-ray diffraction (XRD) analysis (PANalytical X'Pert PRO diffractometer), scanning electron microscopy (SEM)-energy-dispersive X-ray (EDX) analysis (Hitachi, TM3030Plus, Tabletop Microscope), Brunauer-Emmett-Teller (BET) analysis (Micromeritics ASAP2020 system), and zero-point of charge (pH pzc ) value was derived using the methodology indicated by Kaouah et al. [30].

Pretreatment of Raw Clay
To fraction the clay samples, the method followed here was founded on several processes already applied in previous works [31][32][33]: (i) A series of sieves were used in a cascade to perform preliminary pre-sieving on the unprocessed clay to obtain particle sizes of less than 5 µm; (ii) After that, it was ground into a powder gently and then rinsed with H 2 O 2 (40 mL, 6% w/v) to remove any organic compounds; (iii) After that, the resulting solution was transferred to an Erlenmeyer (1000 mL) was added to the buffer solution of pH 4.8 (80 mL) (16 g sodium acetate and 10 mL of acetic acid); (iv) The resulting solution was transferred to a graduated tube with the addition of a dispersible component, sodium hexametaphosphate (NaPO 3 ) 6 , and allowed to decant for 7 h and 45 min after that floating layer was collected at a depth of 10 cm to produce particles smaller than 2 µm. The granules that were acquired were treated numerous times employing deionized water and centrifuged; finally, (v) One glass slide served as the control and was not subjected to any treatment, while the other slide was subjected to heating at 105 • C for two hours to pick the suitable clay. Glass slides were immersed in a glycerin bath until they became saturated and were then analyzed using XRD to identify the presence of the clay mineral Kaolinite: Al 2 Si 2 O 5 (OH) 4 , which had the following composition: 39.7% alumina, 46.2% silica, and 13.8% water. The BET procedure was utilized for calculating the material's specific area, which was 110.786 m 2 /g.

Extracting Cellulose from Red Bean Peels (RBPs)
The procedure shown in Figure 1 was used to isolate cellulose from the powder obtained from RBPs. The process, detailed in several publications [29,34,35], was slightly modified. Initially, 100 g of dried powder from RBPs was cooked in 3 L of hot water with a pH of 7 for 15-25 min and then filtered. The resulting mass was subjected to acid (HCl) pretreatment by mixing it twice with a 1 M HCl solution (500 mL) at 80-90 • C for 60 min; the residue was filtered and collected. For alkali (NaOH) pretreatment, the agitation of the residue occurred three times. In a 1.5 L solution of 1 M NaOH for 60 min at 80-85 • C and then collected. The filtered residue was then subjected to bleach treatment twice with a 4% (w/v) NaOCl solution (at pH = 5 adjusted with 10% (v/v) CH 3 COOH) for 60 min at 80-90 • C, resulting in white-colored cellulose. The cellulose was rinsed with hot deionized water multiple times until the filtrate pH was neutralized. Finally, the cellulose was freeze-dried for seven hours and ground using a mixer grinder. It was then kept at room temperature for further investigation. treatment, while the other slide was subjected to heating at 105 °C for two hours to pick the suitable clay. Glass slides were immersed in a glycerin bath until they became saturated and were then analyzed using XRD to identify the presence of the clay mineral Kaolinite: Al2Si2O5(OH)4, which had the following composition: 39.7% alumina, 46.2% silica, and 13.8% water. The BET procedure was utilized for calculating the material's specific area, which was 110.786 m 2 /g.

Extracting Cellulose from Red Bean Peels (RBPs)
The procedure shown in Figure 1 was used to isolate cellulose from the powder obtained from RBPs. The process, detailed in several publications [29,34,35], was slightly modified. Initially, 100 g of dried powder from RBPs was cooked in 3 L of hot water with a pH of 7 for 15-25 min and then filtered. The resulting mass was subjected to acid (HCl) pretreatment by mixing it twice with a 1 M HCl solution (500 mL) at 80-90 °C for 60 min; the residue was filtered and collected. For alkali (NaOH) pretreatment, the agitation of the residue occurred three times. In a 1.5 L solution of 1 M NaOH for 60 min at 80-85 °C and then collected. The filtered residue was then subjected to bleach treatment twice with a 4% (w/v) NaOCl solution (at pH = 5 adjusted with 10% (v/v) CH3COOH) for 60 min at 80-90 °C, resulting in white-colored cellulose. The cellulose was rinsed with hot deionized water multiple times until the filtrate pH was neutralized. Finally, the cellulose was freezedried for seven hours and ground using a mixer grinder. It was then kept at room temperature for further investigation.

Preparation of Kaol/Cel Composite
To create Kaol/Cel nanocomposites, the researchers followed a specific procedure [36,37]. Initially, acidification is conducted by dissolving 20 g of the pretreated kaolinite into 200 mL of sulfuric acid (H2SO4) (15% W) in the beaker glass and stirring for 4 h at 80 °C. Then, 10 g of activated kaolinite was mixed with 16 mL of (46% W) NaOH in ice water and subjected to magnetic stirring for six hours. To optimize the process according to BBD, the team loaded different ratios of Cel to the modified Kaol before dissolving it with 0.6 M NaOH and 1 M thiourea solutions. Finally, the resulting mixture was vigorously stirred at 80 °C for four hours and then oven-dried at 60 °C overnight.

Preparation of Kaol/Cel Composite
To create Kaol/Cel nanocomposites, the researchers followed a specific procedure [36,37]. Initially, acidification is conducted by dissolving 20 g of the pretreated kaolinite into 200 mL of sulfuric acid (H 2 SO 4 ) (15% W) in the beaker glass and stirring for 4 h at 80 • C. Then, 10 g of activated kaolinite was mixed with 16 mL of (46% W) NaOH in ice water and subjected to magnetic stirring for six hours. To optimize the process according to BBD, the team loaded different ratios of Cel to the modified Kaol before dissolving it with 0.6 M NaOH and 1 M thiourea solutions. Finally, the resulting mixture was vigorously stirred at 80 • C for four hours and then oven-dried at 60 • C overnight.

Experimental Design
BBD is an effective method for optimizing processes since it conforms to a quadratic surface. As a result, it was included in the overall design of the experiment [38]. A total of forty-six separate tests were conducted to assess the impact of five primary independent variables on the removal of CV. To determine the ideal parameters and specify the experimental domain, preliminary tests were conducted. These studies focused on the following aspects: the loading of Cel into Kaol (A), the adsorbent injection (B), the pH (C), the residence period (D), and the temperature (E). The various experimental levels of independent variables and their respective codes are listed in Table 1. The BBD, as well as the statistical data, were examined with the help of software called Stat-Ease Design-Expert (Version 13.0). For the purposes of assessing and forecasting CV removal, the second-order polynomial model was used. Including all square terms, linear terms, and linear-by-linear interaction items, the quadratic response model can be defined as follows: where Y represents the objective for optimizing the response, while k denotes the number of variables being considered. The indices i and j are used to represent the variable numbers, and β 0 is the constant coefficient, with β i and β ii representing the linear and quadratic coefficients, respectively. The term β ij refers to the interaction coefficient, and ε represents a random error. The values X i and X j correspond to the response of CV dye removal and coded values for the independent factors (−1, 0, and +1). In Equation (1), a positive sign implies that the variables have a synergistic effect, whereas a negative sign indicates that they have an antagonistic effect. To evaluate the precision of the model, the researchers employed Analysis of Variance (ANOVA) to examine the coefficients. The ANOVA provided several parameters, including p-value, F-value, determination coefficient (R 2 ), projected determination coefficient (R 2 pred ), adjusted determination coefficient (R 2 adj ), acceptable precision, degree of freedom (d f ), and standard deviation (SD). The experimental data and model precision were evaluated using these parameters [39]. The researchers used a dependable second-order quadratic model equation to predict the optimal value and describe the interactions between the elements. To determine the factors' best values, we solved the regression equation, evaluated the counter-response surface map, and set limitations for the variable levels. To establish the variables' extreme values, preliminary tests have been performed.

Batch Adsorption Studies
To investigate the adsorption capacity of Kaol/Cel for removal of CV dye, a BBD experimental was employed using varying amounts (0-50%) of the adsorbent with 100 mL of dye solution at concentrations of 50-300 mg/L. The experiments were conducted at different temperature intervals (30-60 • C) and pH values ranging from 4 to 10, with adjustments made using 0.1 N HCl and NaOH solutions at different intervals (5-120 min), as detailed in Table 2. Following the adsorption experiments, the residual concentrations were determined using a Cary Series UV-vis spectrophotometer at λ max = 590 nm after samples were spun at 3400 rpm for 10 min. The adsorption capacity of Kaol/Cel for CV dye removal (q e ; mg/g) and the percentage of dye removal (R%) were calculated using Equations (2) and (3), respectively. where: C 0 (mg/L) and C e (mg/L) are the initial and equilibrium CV dye concentrations, respectively, V (L) is the dye solution's volume, and W is Kaol/Cel weight in grams (g).

Adsorbent Characterization
XRD patterns of initial cellulose are shown in Figure 2a-c. For the XRD pattern of kaolinite (a), the well-known diffraction peaks at 2θ around 6.0 • , 19.6 • , 20.8 • 26.5 • and 29.3 • were attributed to the planes of (110), (002), (020) and (021), respectively, all of which were in agreement with the different characteristic orientation of the layers in the crystal structure of kaolinite (Kaol) [19]. For the XRD pattern of cellulose (b), the peaks observed at 2θ around 5.7 • , 14.1 • , 20.8 • , 27.0 • , and 29.0 • are attributed to the (101), (10-1), (002), (110), and (10-2) planes, respectively, in the cellulose structure. These peak positions represent the interplanar spacing and arrangement of the cellulose chains [40]. The XRD pattern of the Kaol/Cel-25 composite(c) shows a significant decrease in peak intensities compared to the individual kaolin and cellulose patterns. This decrease in density is due to integrating the cellulose molecules into the clay framework. In addition, this work provided further evidence that Kaol particles may be successfully synthesized in the Kaol/Cel-25 formulation [19,40].

Adsorbent Characterization
XRD patterns of initial cellulose are shown in Figure 2a-c. For the XRD pattern of kaolinite (a), the well-known diffraction peaks at 2θ around 6.0°, 19.6°, 20.8° 26.5°and 29.3° were attributed to the planes of (110), (002), (020) and (021), respectively, all of which were in agreement with the different characteristic orientation of the layers in the crystal structure of kaolinite (Kaol) [19]. For the XRD pattern of cellulose (b), the peaks observed at 2θ around 5.7°, 14.1°, 20.8°, 27.0°, and 29.0° are attributed to the (101), (10-1), (002), (110), and (10-2) planes, respectively, in the cellulose structure. These peak positions represent the interplanar spacing and arrangement of the cellulose chains [40]. The XRD pattern of the Kaol/Cel-25 composite(c) shows a significant decrease in peak intensities compared to the individual kaolin and cellulose patterns. This decrease in density is due to integrating the cellulose molecules into the clay framework. In addition, this work provided further evidence that Kaol particles may be successfully synthesized in the Kaol/Cel-25 formulation [19,40].  FT-IR analysis was employed to compare the functional groups on the surface of Kaol/Cel-25 before and after CV adsorption, as illustrated in Figure 3a,b, respectively. The FT-IR spectra of Kaol/Cel-25 before CV adsorption ( Figure 3a) showed bands at 1633 cm −1 , which could be attributed to O-H stretching and coordinated water bending vibrations, respectively [16,41]. The Si-O-Si stretching vibrations of kaolinite or quartz cause the peak at 803 cm −1 , whereas the peaks at 3625 cm −1 are caused by Al-OH-Al and Fe-OH-Al deformation. The spectrum of the utilized cellulose was utilized to determine the various chemical groups of cellulose, such as β-glycosidic linkages at approximately 987 cm −1 , -C-O-C pyranose rings at around 1630 cm −1 , -C-H groups at around 995 cm −1 as well as approximately 2157 cm −1 [42]. The FT-IR spectrum of Kaol/Cel-25 after CV dye adsorption (Figure 3b) displayed the same bands in the spectrum before, with a slight shifting of some bands, indicating the involvement of the functional groups of Kaol in the adsorption of CV dye. Moreover, a new band that appeared at the approximately 1744 cm −1 band is attributed to the elongation of the C=N bond. The 1582 cm −1 band is due to aromatic C-C stretches, while the vibration of aromatic tertiary amine N-C is observed at 1366 cm −1 , indicating that the important functional groups of Kaol/Cel-25 were responsible for the adsorption of CV dye [43].
Moreover, a new band that appeared at the approximately 1744 cm −1 band is attributed to the elongation of the C=N bond. The 1582 cm −1 band is due to aromatic C-C stretches, while the vibration of aromatic tertiary amine N-C is observed at 1366 cm −1 , indicating that the important functional groups of Kaol/Cel-25 were responsible for the adsorption of CV dye [43].  Figure 4a depicts the surface morphology of Kaol, which may be described as an uneven and heterogonous surface with crevices. These fissures can be seen across the surface. The EDX spectrum of Kaol reveals the presence of O, Al, Si, K, Fe, and Ti. These elements are present in various minerals such as kaolinite, quartz, and other clay constituents, as identified through XRD analysis. Before CV dye adsorption, the result of the SEM-EDX analysis for Kaol/Cel can be shown in Figure 4b. The appearance of the surface morphology is a surface with protrusions of varying widths with many big spaces and fissures. The EDX analysis verifies the existence of the components O, C, Al, Si, K, Fe, and Ca. This suggests that the clay matrix effectively incorporated the cellulose particles.
Nevertheless, the morphological structure of Kaol/Cel-25 following CV dye absorption (Figure 4c) was converted into a different state, resulting in a reduced number of many tiny apertures and surface slits of the material. This observation provides evidence that molecules of CV dye have been loaded onto the surface of Kaol/Cel-25. In addition, the EDX examination reveals a spike in the carbonation rate in the related EDX spectrum, which provides more evidence that CV is being adsorbed to the Kaol/Cel-25 surface.

Fitting the Process Models
As was indicated before, BBD was responsible for the design of forty-six (46) different experiments ( Table 2). The BBD method was utilized to explore the individual and interaction impacts that each of the examined factors had on the CV removal efficiency (as a Nevertheless, the morphological structure of Kaol/Cel-25 following CV dye absorption (Figure 4c) was converted into a different state, resulting in a reduced number of many tiny apertures and surface slits of the material. This observation provides evidence that molecules of CV dye have been loaded onto the surface of Kaol/Cel-25. In addition, the EDX examination reveals a spike in the carbonation rate in the related EDX spectrum, which provides more evidence that CV is being adsorbed to the Kaol/Cel-25 surface.

Fitting the Process Models
As was indicated before, BBD was responsible for the design of forty-six (46) different experiments ( Table 2). The BBD method was utilized to explore the individual and interaction impacts that each of the examined factors had on the CV removal efficiency (as a response). The examined factors were the temperature (D), the pH (C), the adsorbent dose (A), and the temperature (D). These four factors were considered to be independent process factors.
The quadratic polynomial model established the mathematical relationship between the process components and the response. The regression model's significance test for the response was conducted, and the findings of the ANOVA are presented in Table 3. The model's F-value is 5.66, and the p-value is less than 0.0001, indicating the significance of the model terms. A p-value of < 0.05 is seen as important under the selected circumstances. The essential model terms for the response, i.e., the elimination of CV, were identified as A, B, BC, BD, and A 2 . However, it was observed that the residence period factor (E) had no effect on the CV elimination, possibly due to the adsorption process's low sensitivity to time changes. The empirical relationship between CV elimination and the significant variables is expressed by the quadratic regression model given in Equation (4).
CV removal (%) = +97.97 + 11.83A + 5.89B − 8.95BC + 6.53BD − 9.66A 2 (4) A positive sign in Equation (4) indicates that the factor has a synergistic impact, while a negative sign means that the factor has an antagonistic effect [44]. As shown in Figure 5, the R-squared (determination coefficient) for the response factors is 0.97, which is a very high number indicating an excellent correlation between the real and projected values.

Interactions Significant for Crystal Violet (CV) Dye Removal
The influence of the adsorbent dose (B) and the solution pH (C) on the CV elimination efficacy showed a significant interaction (as evidenced by the p-value of 0.0085 from Table  3). The remaining independent variables (i.e., Cel loading (Kaol/Cel-25), residence period of 5 min, and temperature of 45 °C), were held constant during the experimentation range. Figure 6a,b depict the 3D surface and 2D contour plots for the BC interaction, respectively. Figure 5 depicts that augmenting the maximum adsorbent injection from 0.02 g to 0.05 g and decreasing the solution pH from 10 to 4 significantly augmented CV elimination efficacy from 61.82% to 99.86%. The greatest CV removal efficiency was achieved at pH 4, with a gradual decrease in dye removal as the pH value increased toward an alkaline environment. Furthermore, Figure 7 shows that the pHpzc of Kaol/Cel-25 was 8.71, establishing that the surface of Kaol/Cel-25 begins to be positively charged at pH values less than pHpzc. Conversely, at pH values above pHpzc, the surface charge of Kaol/Cel-25 becomes negative, which implies that Kaol/Cel-25 may adsorb cationic dyes. Consequently, more essential electrostatic attractions arise between the surface functional groups of negatively charged Kaol/Cel-25 and cationic CV dye, as shown in Equation (5).

Interactions Significant for Crystal Violet (CV) Dye Removal
The influence of the adsorbent dose (B) and the solution pH (C) on the CV elimination efficacy showed a significant interaction (as evidenced by the p-value of 0.0085 from Table 3). The remaining independent variables (i.e., Cel loading (Kaol/Cel-25), residence period of 5 min, and temperature of 45 • C), were held constant during the experimentation range. Figure 6a,b depict the 3D surface and 2D contour plots for the BC interaction, respectively.

Interactions Significant for Crystal Violet (CV) Dye Removal
The influence of the adsorbent dose (B) and the solution pH (C) on the CV elimination efficacy showed a significant interaction (as evidenced by the p-value of 0.0085 from Table  3). The remaining independent variables (i.e., Cel loading (Kaol/Cel-25), residence period of 5 min, and temperature of 45 °C), were held constant during the experimentation range. Figure 6a,b depict the 3D surface and 2D contour plots for the BC interaction, respectively. Figure 5 depicts that augmenting the maximum adsorbent injection from 0.02 g to 0.05 g and decreasing the solution pH from 10 to 4 significantly augmented CV elimination efficacy from 61.82% to 99.86%. The greatest CV removal efficiency was achieved at pH 4, with a gradual decrease in dye removal as the pH value increased toward an alkaline environment. Furthermore, Figure 7 shows that the pHpzc of Kaol/Cel-25 was 8.71, establishing that the surface of Kaol/Cel-25 begins to be positively charged at pH values less than pHpzc. Conversely, at pH values above pHpzc, the surface charge of Kaol/Cel-25 becomes negative, which implies that Kaol/Cel-25 may adsorb cationic dyes. Consequently, more essential electrostatic attractions arise between the surface functional groups of negatively charged Kaol/Cel-25 and cationic CV dye, as shown in Equation (5).   Figure 5 depicts that augmenting the maximum adsorbent injection from 0.02 g to 0.05 g and decreasing the solution pH from 10 to 4 significantly augmented CV elimination efficacy from 61.82% to 99.86%. The greatest CV removal efficiency was achieved at pH 4, with a gradual decrease in dye removal as the pH value increased toward an alkaline environment. Furthermore, Figure 7 shows that the pH pzc of Kaol/Cel-25 was 8.71, establishing that the surface of Kaol/Cel-25 begins to be positively charged at pH values less than pH pzc . Conversely, at pH values above pH pzc , the surface charge of Kaol/Cel-25 becomes negative, which implies that Kaol/Cel-25 may adsorb cationic dyes. Consequently, more essential electrostatic attractions arise between the surface functional groups of negatively charged Kaol/Cel-25 and cationic CV dye, as shown in Equation (5). Figure 5 depicts that augmenting the maximum adsorbent injection from 0.02 g to 0.05 g and decreasing the solution pH from 10 to 4 significantly augmented CV elimination efficacy from 61.82% to 99.86%. The greatest CV removal efficiency was achieved at pH 4, with a gradual decrease in dye removal as the pH value increased toward an alkaline environment. Furthermore, Figure 7 shows that the pHpzc of Kaol/Cel-25 was 8.71, establishing that the surface of Kaol/Cel-25 begins to be positively charged at pH values less than pHpzc. Conversely, at pH values above pHpzc, the surface charge of Kaol/Cel-25 becomes negative, which implies that Kaol/Cel-25 may adsorb cationic dyes. Consequently, more essential electrostatic attractions arise between the surface functional groups of negatively charged Kaol/Cel-25 and cationic CV dye, as shown in Equation (5).   Table 3 was used to assess the significance of the interactions between the independent variables. The results showed that the interaction influence between adsorbent dosage (B) and temperature (D) had a statistically considerable impact on CV removal efficiency, with a p-value of 0.0477. The other variables, namely a 25% Cel loading, a pH of 10, and a 5 min contact time, remained constant throughout the experiment. The interaction between the adsorbent dose and the temperature was further analyzed using 3D and 2D response surface plots, as shown in Figure 8a,b, respectively. The plots indicated that augmenting the adsorbent dose from 0.02 g to 0.05 g led to higher CV dye removal (%), possibly due to increased active adsorption sites or available surface area. Nevertheless, the temperature seemed to possess no significant influence on eliminating CV dye, suggesting an exothermic adsorption process. These results will be discussed further in Section 3.7 on Adsorption Thermodynamics. It was determined whether there was a statistically crucial interaction between each pair of independent variables (as shown in Table 3). The results revealed that the interaction influence between adsorbent dosage (B) and temperature (D) had a significant impact on the CV removal efficiency (p-value = 0.0477). The other independent factors, including Cel loading of 25%, solution pH of 10, and contact time of 5 min, remained constant throughout the experiment. Figure 8a,b illustrate the response surface plots for the interaction between temperature and adsorbent injection. The findings show that augmenting  Table 3 was used to assess the significance of the interactions between the independent variables. The results showed that the interaction influence between adsorbent dosage (B) and temperature (D) had a statistically considerable impact on CV removal efficiency, with a p-value of 0.0477. The other variables, namely a 25% Cel loading, a pH of 10, and a 5 min contact time, remained constant throughout the experiment. The interaction between the adsorbent dose and the temperature was further analyzed using 3D and 2D response surface plots, as shown in Figure 8a,b, respectively. The plots indicated that augmenting the adsorbent dose from 0.02 g to 0.05 g led to higher CV dye removal (%), possibly due to increased active adsorption sites or available surface area. Nevertheless, the temperature seemed to possess no significant influence on eliminating CV dye, suggesting an exothermic adsorption process. These results will be discussed further in Section 3.7 on Adsorption Thermodynamics.   Table 3 was used to assess the significance of the interactions between the independent variables. The results showed that the interaction influence between adsorbent dosage (B) and temperature (D) had a statistically considerable impact on CV removal efficiency, with a p-value of 0.0477. The other variables, namely a 25% Cel loading, a pH of 10, and a 5 min contact time, remained constant throughout the experiment. The interaction between the adsorbent dose and the temperature was further analyzed using 3D and 2D response surface plots, as shown in Figure 8a,b, respectively. The plots indicated that augmenting the adsorbent dose from 0.02 g to 0.05 g led to higher CV dye removal (%), possibly due to increased active adsorption sites or available surface area. Nevertheless, the temperature seemed to possess no significant influence on eliminating CV dye, suggesting an exothermic adsorption process. These results will be discussed further in Section 3.7 on Adsorption Thermodynamics. It was determined whether there was a statistically crucial interaction between each pair of independent variables (as shown in Table 3). The results revealed that the interaction influence between adsorbent dosage (B) and temperature (D) had a significant impact on the CV removal efficiency (p-value = 0.0477). The other independent factors, including Cel loading of 25%, solution pH of 10, and contact time of 5 min, remained constant throughout the experiment. Figure 8a,b illustrate the response surface plots for the interaction between temperature and adsorbent injection. The findings show that augmenting It was determined whether there was a statistically crucial interaction between each pair of independent variables (as shown in Table 3). The results revealed that the interaction influence between adsorbent dosage (B) and temperature (D) had a significant impact on the CV removal efficiency (p-value = 0.0477). The other independent factors, including Cel loading of 25%, solution pH of 10, and contact time of 5 min, remained constant throughout the experiment. Figure 8a,b illustrate the response surface plots for the interaction between temperature and adsorbent injection. The findings show that augmenting the adsorbent injection from 0.02 g to 0.05 g resulted in higher removal of CV, which can be related to an elevation in active adsorption sites and/or surface area. Moreover, the results suggest that the adsorption of dye molecules onto the Kaol/Cel-25 surface was exothermic since the temperature had no discernible influence on CV removal effectiveness. A more detailed discussion on adsorption thermodynamics is included in Section 3.7.

Adsorption Studies
In this work, the influence of residence period vs. initial CV level on the adsorption potential of Kaol/Cel-25 was explored over a range of initial levels from 50 to 300 mg/L. Figure 9 shows the experimental results, while the optimal conditions for CV dye removal were determined based on the highest achieved CV dye removal percentage, observed at a constant Kaol/Cel-25 adsorbent injection (0.035 g), pH 10, and 45 • C, as indicated in Run #24 of Table 2. As the initial CV level was augmented, the equilibrium adsorption capacity was also augmented from 51.7 to 297.7 mg/g, as presented in Figure 9. This increase is a result of a greater collision rate between CV dye and the Kaol/Cel-25 surface, resulting from the elevated initial CV concentration. Additionally, the higher concentration gradient facilitated the diffusion of dye molecules into the internal pores of the adsorbent, leading to their migration towards the active adsorption sites, thus enhancing the adsorption capacity [45]. the adsorbent injection from 0.02 g to 0.05 g resulted in higher removal of CV, which can be related to an elevation in active adsorption sites and/or surface area. Moreover, the results suggest that the adsorption of dye molecules onto the Kaol/Cel-25 surface was exothermic since the temperature had no discernible influence on CV removal effectiveness. A more detailed discussion on adsorption thermodynamics is included in Section 3.7.

Adsorption Studies
In this work, the influence of residence period vs. initial CV level on the adsorption potential of Kaol/Cel-25 was explored over a range of initial levels from 50 to 300 mg/L. Figure 9 shows the experimental results, while the optimal conditions for CV dye removal were determined based on the highest achieved CV dye removal percentage, observed at a constant Kaol/Cel-25 adsorbent injection (0.035 g), pH 10, and 45 °C, as indicated in Run #24 of Table 2. As the initial CV level was augmented, the equilibrium adsorption capacity was also augmented from 51.7 to 297.7 mg/g, as presented in Figure 9. This increase is a result of a greater collision rate between CV dye and the Kaol/Cel-25 surface, resulting from the elevated initial CV concentration. Additionally, the higher concentration gradient facilitated the diffusion of dye molecules into the internal pores of the adsorbent, leading to their migration towards the active adsorption sites, thus enhancing the adsorption capacity [45].

Kinetic Modeling
To comprehend the adsorption pathway and behavior of CV on the surface of Kaol/Cel-25, further investigations were carried out, and understanding the adsorption kinetics is crucial. In this context, two kinetic models, pseudo-first-order (PFO) and pseudo-second-order (PSO), were utilized. These non-linear kinetic models are represented by Equation (6) for PFO [46] and Equation (7) for PSO [47]: = 1 + (7) Figure 9. Influence of CV adsorption contact time vs. initial concentrations on Kaol/Cel-25.

Kinetic Modeling
To comprehend the adsorption pathway and behavior of CV on the surface of Kaol/Cel-25, further investigations were carried out, and understanding the adsorption kinetics is crucial. In this context, two kinetic models, pseudo-first-order (PFO) and pseudo-secondorder (PSO), were utilized. These non-linear kinetic models are represented by Equation (6) for PFO [46] and Equation (7) for PSO [47]: designated by q e (mg/g). In contrast, the quantity of CV dye adsorbed at a specific time t is designated by q t (mg/g). Additionally, the rate constant of the PFO is represented by k 1 (1/min), and the rate constant of the PSO is represented by k 2 (g/mg min). Table 4 displays the correlation coefficients (R 2 ) and model parameters for the PSO and PFO. On the basis of the kinetic adsorption data in Table 4, it can be shown that the PSO model fits the Kaol/Cel-25 dye adsorption on the Kaol/Cel-25 surface better than the PFO model. This is confirmed by the fact that the PSO model had higher R 2 values than the PFO model. Furthermore, the q e (i.e., q e,cal ) values estimated with the PSO model are closer to the experimental q e (i.e., q e,exp ) values than the q e values calculated with the PFO model. These findings and FT-IR analysis suggest that chemical interaction between CV dye and active functional groups on the Kaol/Cel-25 surface is essential in the adsorption process. In the Kaol/Cel-25 surface, the adsorption of CV dye is thus governed by the chemisorption phenomenon.

Isotherms for Adsorption
The adsorption isotherm is a critical piece of information to have when attempting to describe the interaction between the molecules of CV dye and Kaol/Cel-25. As a result, the isotherms that are applied the most often, namely Langmuir [48], Freundlich [49], and Temkin [50], are selected for the purpose of conducting an analysis of the data pertaining to equilibrium adsorption and determining the adsorption potential of Kaol/Cel-25. The non-linear forms of the Langmuir, Freundlich, and Temkin equations are each represented in their own separate Equations (8)- (10): The variables used in the analysis include C e (mg/L) for the equilibrium concentration of CV dye, q max (mg/g) for the maximum quantity of CV dye that can be adsorbed per unit mass of Kaol/Cel-25, and q e (mg/g) for the quantity of CV dye absorbed per unit weight. The constants used include Langmuir (K a ), Freundlich (K f ), and Temkin (K T ) constants (L/mg), as well as n for the dimensionless constant indicating adsorption intensity, b T (J/mol) for the heat of adsorption, T (K) for temperature, and R (8.314 J/mol K) for the universal gas constant.
The absorption of CV dye by Kaol/Cel-25 is shown in Table 5 to be best fit by the Langmuir adsorption isotherm model, which has a higher R 2 (0.99) than the Freundlich and Temkin models. This implies that CV dye form a monolayer coverage on the homogeneous surface of the adsorbent that are energetically equivalent [51]. Moreover, using the Langmuir model, we determined that 294.12 mg/g is the q max for Kaol/Cel-25. Table 6 compares the q max of CV dye adsorption onto Kaol/Cel-25 to that of other adsorbents that have been published for removing CV, highlighting that the cationic dye removal efficiency of Kaol/Cel-25 is highly effective and shows promise for further applications.

Thermodynamic Functions Results
To assess the viability and spontaneity of the CV dye adsorption phenomenon onto the surface of Kaol/Cel-25 and estimate the level of randomness at the interface between the dye and the surface, various adsorption thermodynamic parameters were determined. These parameters include the Gibbs free energy change (∆G • ) (kJ/mol), the entropy change (∆S • ) (kJ/mol. K) and the enthalpy change (∆H • ) (kJ/mol). Such parameters thermodynamic were computed using Equations (11)-(13) [61]: Figure 10 depicts a plot of lnK d against 1/T, from which the thermodynamic parameters (∆H • and ∆S • ) were calculated. The slope of this plot is ∆H • , and the intercept is ∆S • . Figure 10 depicts a plot of lnKd against 1/T, from which the thermodynamic par ters (ΔH° and ΔS°) were calculated. The slope of this plot is ΔH°, and the intercept is The results presented in Table 7 show that the Gibbs free energy change (ΔG°) fo adsorption of CV dye onto the surface of Kaol/Cel-25 is negative, indicating that the nomenon is spontaneous [62]. Furthermore, the negative enthalpy values obtained fo adsorption process of CV by Kaol/Cel-25 suggest that the process is exothermic, wh consistent with the results obtained from the BBD parametric optimization depict Figure 10. Additionally, the adsorption of CV onto Kaol/Cel-25 appears to cause a gr level of disorder at the interface between the solid and solution, as indicated by the tive entropy values [63].

Mechanisms of Adsorption
A proposed adsorption mechanism for CV by Kaol/Cel-25 is illustrated in Figu The various functional groups existing on the surface of Kaol/Cel-25 participate in d ent interactions that facilitate the adsorption pathway of CV. Notably, the most imp interaction is an electrostatic attraction between the CV dye molecules and the surfa the Kaol/Cel-25 adsorbent. Such an adsorption pathway encompasses the electrosta teraction between CV dye cations and the negatively charged sites on the surfa Kaol/Cel-25. Additionally, there is the possibility of n-π interaction, which typical curs when the lone pair electrons on an O are spread out into the orbital of an aro ring of dye [64]. Two kinds of hydrogen bonding take place between the Kaol/Cel-2 the molecular structure of CV. The first and more frequent kind is dipole-dipole hydr bonding interaction between present O on the surface of Kaol/Cel-25 with the O a atoms existing in the CV molecular structure. This interaction is illustrated in Figu The last kind is Yoshida H-bonding, which takes place between -OH groups on the su of Kaol/Cel-25 and the CV aromatic ring [65]. Based on the information presented a it can be inferred that these interactions had a crucial contribution in improving th sorption phenomenon of CV dye on the surface of Kaol/Cel-25. The results presented in Table 7 show that the Gibbs free energy change (∆G • ) for the adsorption of CV dye onto the surface of Kaol/Cel-25 is negative, indicating that the phenomenon is spontaneous [62]. Furthermore, the negative enthalpy values obtained for the adsorption process of CV by Kaol/Cel-25 suggest that the process is exothermic, which is consistent with the results obtained from the BBD parametric optimization depicted in Figure 10. Additionally, the adsorption of CV onto Kaol/Cel-25 appears to cause a greater level of disorder at the interface between the solid and solution, as indicated by the negative entropy values [63].

Mechanisms of Adsorption
A proposed adsorption mechanism for CV by Kaol/Cel-25 is illustrated in Figure 11. The various functional groups existing on the surface of Kaol/Cel-25 participate in different interactions that facilitate the adsorption pathway of CV. Notably, the most impactful interaction is an electrostatic attraction between the CV dye molecules and the surface of the Kaol/Cel-25 adsorbent. Such an adsorption pathway encompasses the electrostatic interaction between CV dye cations and the negatively charged sites on the surface of Kaol/Cel-25. Additionally, there is the possibility of n-π interaction, which typically occurs when the lone pair electrons on an O are spread out into the orbital of an aromatic ring of dye [64]. Two kinds of hydrogen bonding take place between the Kaol/Cel-25 and the molecular structure of CV. The first and more frequent kind is dipole-dipole hydrogen bonding interaction between present O on the surface of Kaol/Cel-25 with the O and N atoms existing in the CV molecular structure. This interaction is illustrated in Figure 11. The last kind is Yoshida H-bonding, which takes place between -OH groups on the surface of Kaol/Cel-25 and the CV aromatic ring [65]. Based on the information presented above, it can be inferred that these interactions had a crucial contribution in improving the adsorption phenomenon of CV dye on the surface of Kaol/Cel-25.

Conclusions
This research aimed to develop an environmentally friendly composite material called kaolinite-cellulose (Kaol/Cel) using waste red bean peels (Phaseolus vulgaris) as a cellulose source. A Box-Behnken design was employed to perform parametric optimization of the synthesis conditions and adsorptive performance of the Kaol/Cel composite to remove crystal violet (CV) dye. The composite's polymeric matrix with uniform Cel distribution was found to contribute to its elevated adsorptive behavior towards CV. The best parameters for CV elimination (99.58%) were identified as 25% Cel loading, 0.035 g adsorbent dosage, pH 10, 45 °C temperature, and 5 min contact time. The equilibrium data were best represented by the Freundlich isotherm model, and the Langmuir model revealed a maximum adsorption capacity (qmax) of 294.12 mg/g at 45 °C. The adsorption mechanism involved several interactions, including electrostatic attractions, n-π stacking interactions, dipole-dipole hydrogen bonding interactions, and Yoshida H-bonding. Additionally, the adsorption process of CV by Kaol/Cel-25 was observed to be exothermic and spontaneous based on adsorption thermodynamic functions. These findings suggest that Kaol/Cel-25 could be a cost-effective composite biosorbent for removing cationic dyes from aquatic environments and could have potential applications in wastewater treatment. This refers to the elimination of organic and inorganic pollutants, as well as the diminution of the chemical oxygen demand.