Next Article in Journal
Carbonic Anhydrase 2 Deletion Delays the Growth of Kidney Cysts Whereas Foxi1 Deletion Completely Abrogates Cystogenesis in TSC
Next Article in Special Issue
Towards Enhanced Tunability of Aqueous Biphasic Systems: Furthering the Grasp of Fluorinated Ionic Liquids in the Purification of Proteins
Previous Article in Journal
Studies on the Thermal Decomposition Course of Nitrogen-Rich Heterocyclic Esters as Potential Drug Candidates and Evaluation of Their Thermal Stability and Properties
Previous Article in Special Issue
Comprehensive Evaluation of Polyaniline-Doped Lignosulfonate in Adsorbing Dye and Heavy Metal Ions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sorption Behavior of Azo Dye Congo Red onto Activated Biochar from Haematoxylum campechianum Waste: Gradient Boosting Machine Learning-Assisted Bayesian Optimization for Improved Adsorption Process

by
Diego Melchor Polanco Gamboa
1,
Mohamed Abatal
1,*,
Eder Lima
2,
Francisco Anguebes Franseschi
3,
Claudia Aguilar Ucán
3,
Rasikh Tariq
4,
Miguel Angel Ramírez Elías
3 and
Joel Vargas
5
1
Facultad de Ingeniería, Universidad Autónoma del Carmen, Ciudad del Carmen 24115, Campeche, Mexico
2
Institute of Chemistry, Federal University of Rio Grande do Sul (UFRGS), Av. Bento Goncalves 9500, P.O. Box 15003, Porto Alegre 91501-970, RS, Brazil
3
Facultad de Química, Universidad Autónoma del Carmen, Calle 56 No. 4 Av. Concordia, Ciudad del Carmen 24180, Campeche, Mexico
4
Tecnologico de Monterrey, Institute for the Future of Education, Ave. Eugenio Garza Sada 2501, Monterrey 64849, Nuevo León, Mexico
5
Instituto de Investigaciones en Materiales, Unidad Morelia, Universidad Nacional Autónoma de México, Antigua Carretera a Pátzcuaro No. 8701, Col. Ex Hacienda de San José de la Huerta, Morelia 58190, Michoacán, Mexico
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(9), 4771; https://doi.org/10.3390/ijms25094771
Submission received: 9 April 2024 / Revised: 24 April 2024 / Accepted: 24 April 2024 / Published: 27 April 2024

Abstract

:
This work aimed to describe the adsorption behavior of Congo red (CR) onto activated biochar material prepared from Haematoxylum campechianum waste (ABHC). The carbon precursor was soaked with phosphoric acid, followed by pyrolysis to convert the precursor into activated biochar. The surface morphology of the adsorbent (before and after dye adsorption) was characterized by scanning electron microscopy (SEM/EDS), BET method, X-ray powder diffraction (XRD), and Fourier-transform infrared spectroscopy (FTIR) and, lastly, pHpzc was also determined. Batch studies were carried out in the following intervals of pH = 4–10, temperature = 300.15–330.15 K, the dose of adsorbent = 1–10 g/L, and isotherms evaluated the adsorption process to determine the maximum adsorption capacity (Qmax, mg/g). Kinetic studies were performed starting from two different initial concentrations (25 and 50 mg/L) and at a maximum contact time of 48 h. The reusability potential of activated biochar was evaluated by adsorption–desorption cycles. The maximum adsorption capacity obtained with the Langmuir adsorption isotherm model was 114.8 mg/g at 300.15 K, pH = 5.4, and a dose of activated biochar of 1.0 g/L. This study also highlights the application of advanced machine learning techniques to optimize a chemical removal process. Leveraging a comprehensive dataset, a Gradient Boosting regression model was developed and fine-tuned using Bayesian optimization within a Python programming environment. The optimization algorithm efficiently navigated the input space to maximize the removal percentage, resulting in a predicted efficiency of approximately 90.47% under optimal conditions. These findings offer promising insights for enhancing efficiency in similar removal processes, showcasing the potential of machine learning in process optimization and environmental remediation.

1. Introduction

Water pollution is a serious concern because almost all industrial activities generate pollutants that are discharged into water resources without treatment or with a deficiency treatment. Dyes are an example of stable organic compounds used in many industries, such as textile, paper, plastic, printing, cosmetics, and others. When dyes reach aquatic environments, they risk the aquatic ecosystem and even human health because some are toxic, carcinogenic, and mutagenic [1]. In addition, small amounts of dye cause a change in the aesthetic perception of the water and also affect the penetration of sunlight and the amount of oxygen dissolved in water [2].
Many types of dyes are on the market, and their classification depends mainly on how they are applied to the substrate and their chemical structure. The first one includes direct dyes, reactive dyes, vat dyes, and disperse dyes, and the latter includes azo dyes, indigo dyes, acid dyes, basic dyes, and anthraquinone dyes, among others [3]. About 50% of the dyes on the market are azo dyes, which is why azo dyes are considered the most important widely used dyes [4].
Congo red (CR) is a di-azo dye containing two azo groups (-N=N-) which are attached to two aromatic radicals. CR is widely used in the textile industry due to the fact that it can easily form complexes with polysaccharides such as cellulose [5,6]. Another application of Congo red is in histology, identifying amyloid proteins in tissues to diagnose amyloidosis [7].
In the aqueous solution, the molecular structure of CR is affected by pH variations, which give different colors. Congo red solutions remain red at neutral and alkaline pH values but turn blue at acid pH values [6]. Thus, CR is also used as an acid–base indicator in several industries. As mentioned above, CR dye has different uses in various industries that generate waste that can be discharged into water resources. The problem with azo dyes is their toxicity; when they degrade under anaerobic conditions, the benzidine formed is a known carcinogen. Therefore, benzidine is a byproduct of CR degradation that may cause bladder cancer in humans [8]. In addition to their carcinogenic effect, azo dyes have other toxic effects that can affect human health by damaging organs such as skin, kidneys, liver, and even central nervous system [9].
According to its adverse effects on humans and plants [10,11], Congo red dye presents a significant risk to human health and the environment, so it is vital to achieve its elimination from aqueous media. There are many technologies for removing dyes, particularly for Congo red dye in aqueous solution. The most common methods applied are biological degradation, photocatalysis, and adsorption. Microorganisms have developed resistance to high dye concentrations and may degrade them into carbon dioxide, inorganic salts, and water by enzymatic action. Different species of genus Aspergillus fungi, such as Aspergillus flavus [12], Aspergillus niger [13], and Aspergillus fumigatus [14], have been successfully used to biodegrade CR. Photocatalysis is an advanced removal technology that has received more attention recently due to its main advantages, such as rapid dye degradation and no generation of secondary waste [15]. Recent and novel investigations to synthesize new photocatalysts and apply them to Congo red removal have been carried out, such as orange peel extract biosynthesized zinc oxide nanoparticles being used to remove CR via photocatalysis [16] and novel rare earth metal doped ZnO photocatalysts for degradation of CR [17].
Adsorption has also been used for the removal of CR. Some biosorbents such as Antigonon leptopus leaf powder [18], mango leaves powder [19], Moringa oleifera seed coat [20], and peel waste of Hylocereus undatushave (white dragon fruit) [21] have been used to remove CR dye efficiently. On the other hand, the most used adsorbent is commercial activated carbon due to its high uptake capacity, high surface area, and environmental friendliness. Non-toxicity is one of its advantages. However, its main disadvantage is its high operating costs [22]. So, a need for new materials for synthesizing activated biochars has risen. Different biomass, such as Guar gum [23] and walnut shell [24], have been used to produce activated carbon.
In this work, activated biochar was prepared from Haematoxylum campechianum waste, taking advantage of its abundance in the state of Campeche in Mexico. This material was used to remove the azo dye CR in aqueous solution. ABHC was characterized before and after adsorption of the dye by SEM/EDS, BET method, X-ray powder diffraction (XRD), and Fourier-transform infrared spectroscopy (FTIR), and, lastly, pHpzc was also determined. The effects of contact time, dye concentration, temperature, pH of the aqueous medium, and the activated biochar dose were evaluated to compare and evaluate the removal efficiency. Cycles of adsorption–desorption were performed to determine the reusability of the activated biochar.
Another research contribution of this work is to present a novel approach to optimize chemical removal processes using advanced machine learning techniques. The efficient removal of contaminants from various mediums is a critical challenge in environmental engineering and remediation. Traditional optimization methods often rely on labor-intensive experimentation and trial-and-error approaches. Here, we contribute to the field by proposing a data-driven approach that leverages the power of Gradient Boosting regression models and Bayesian optimization. By harnessing a rich dataset containing information on initial concentration, time, temperature, pH, and dose, we aim to maximize the removal percentage through precise tuning of input variables. This research not only offers a pragmatic solution to enhance removal efficiency but also underscores the potential of machine learning in addressing complex environmental challenges.

2. Results and Discussion

2.1. Characterization of ABHC

As shown in Figure 1, the point of zero charge of ABHC was 6.5. This result is similar to those obtained from other biomass precursor materials such as Dipterocarpus alatus (pHpzc = 6.3) [25] and rice husk (pHpzc = 6.8) [26]. Therefore, when the solution pH is above the pHpzc (pH > pHpzc), the surface of ABHC is negatively charged and the cationic species will be preferentially removed. In contrast, if values of pH are below pHpzc (pH < pHpzc), the charge of the surface of ABHC will become positive, and then anionic species are preferentially attracted via electrostatic interactions [27].
XRD patterns of ABHC before and after Congo red adsorption are presented in Figure 2. For ABHC, there is a highest intensity diffraction peak around 2θ = 26.5, which can be assigned to the crystalline hexagonal phase of graphite [28]. There are other diffraction peaks, with minor intensity, around 2θ = 22.3 and 42.9, which are attributed to crystalline graphite [29]. After CR adsorption, the diffraction pattern did not present significant changes, indicating that there is no change in the crystallinity of ABHC as a result of the adsorption process [24]. Similar results were observed in previous studies. Homagai et al. (2022) reported that the XRD patterns did not show remarkable changes after adsorption of crystal violet on modified rice husk. This result was due to the amorphous nature attributed to the lignin, cellulose, and hemicellulose present in the biosorbents [30]. Furthermore, Li et al. (2020) reported that there was no change in XRD patterns after adsorption of Congo red and methylene blue dye on walnut shell-based activated carbon and also attributed their result to the amorphous nature of the adsorbent [24]. Therefore, it can be assumed that the similarity on XRD diffraction patterns of ABHC before and after Congo red adsorption could be attributed to the presence of lignin, cellulose, and hemicellulose, which are confirmed by FT-IR.
The grain size of the ABHC and ABHC-CR samples was calculated from the XRD patterns using the Debye–Scherrer Equation (1).
D = K λ B C o s ( θ B )
where λ is the wavelength (λ(CuKα) = 0.15405 nm), K is an empirical constant related to the crystalline shape (for no spherical shapes, K = 0.9), B is the full width at half maximum of the peak (FWHM) in radian on the 2θ scale, and θ is the Bragg’s diffracting angle corresponding to the maximum of the peak. FWHM were determined using the Origin Pro version 2021 based on the Gaussian line.
According to Equation (1), the average particle size of the ABHC and ABHC-CR samples was 25.56 nm and 25.69 nm, respectively. This result implies that the particle size of AC was not affected after adsorption of CR.
Scanning electron microscope (SEM) analysis was carried out to investigate the physical surface morphology of ABHC before and after dye adsorption.
The SEM micrographs of ABHC (Figure 3a) show that the particles of the synthesized biochar material before dye adsorption have a rough surface and an irregular shape with a variety of randomly distributed cavities, which can provide easy access transport toward the adsorption sites [31]. After Congo red adsorption, ABHC has slight changes on its surface, and some cavities were filled with CR molecules. Additionally, SEM micrographs after CR adsorption (Figure 3b) show new flake-like deposits formed on the activated biochar surface due to the interactions between the functional groups on the activated biochar and the dye molecule. The elemental composition of the biochar material was performed by energy dispersive X-ray spectroscopy (EDS), as shown in Figure 3c. In ABHC, the material consists predominantly of carbon and oxygen, and the sum of these two elements is 98.0% per weight. The rest of the composition (2.0%) corresponds to metallic fractions (Ca, Al, and P). Lastly, the elemental analysis presents 1.7% per weight of N and 0.2% of S (Figure 3d). Nitrogen is found in the Congo red molecule (azo dye). Additionally, sulfur is also found in the molecule of the dye (Figure 16), which confirms the adsorption of dye onto the adsorbent.
Figure 4 shows the adsorption–desorption isotherm of N2 on the ABHC sample. As observed in this figure, the isotherm belonged to type IV according to International Union of Pure and Applied Chemistry (IUPAC classification), indicating that the ABHC had mesoporous structure [32].
The mean pore diameter of ABHC calculated by Barrett, Joyner, and Halenda (BJH) equations was 2.14 nm, with a surface area of 124.15 m2/g. The difference in the surface properties can be attributed to the type of biomass precursor. The mean pore diameter value of biochar allows it to be classified as a mesoporous type, and several activated carbons with mesoporous type have been used to remove Congo red dye effectively [33]. Additionally, the width of the CR molecule is approximately 2.62 nm [34], requiring a mesoporous type for being adsorbed. Lafi et al. (2019) [35] prepared activated carbon using coffee waste as a precursor and reported a mean pore diameter of 4.04 nm and surface area of 219.69 m2/g and removed CR efficiently. Values of pore diameter and surface area of ABHC are close to those mentioned above. Therefore, it is suspected that ABHC can be used to remove CR dye from the aqueous solution.
Figure 5 shows the FTIR spectrum for ABHC and ABHC-CR, and the assigned bands are summarized in Table 1. For activated biochar from Haematoxylum campechianum (ABHC) broadband at 3358 cm−1 attributed to OH groups present in cellulose [33], the band at 1705 cm−1 could be attributed to C=O present in lignin [5], and the band at 1224 cm−1 can represent C-O stretching of aryl, carboxylates, or ether groups presents in lignin [33]. The groups mentioned above are known to participate in pollutant adsorption processes. After Congo red adsorption, new bands appeared. At 1421 cm−1, this band can be attributed to C=C stretching in the aromatic ring [36], and the band at 1045 cm−1 corresponds to S=O [37]. These functional groups are present in the Congo red dye molecule (Figure 16) and confirm the presence of the adsorbate in the adsorbent. The change in the bands from 3358, 779, and 1224 cm−1 to 3368, 769, and 1196 cm−1 indicates that OH, N-H, and C-O groups are mainly involved in the adsorption of CR onto ABHC and suggests hydrogen bonding, n–π, and π–π interactions [38] as part of the dye adsorption interactions.

2.2. Kinetic Study

The behavior of the experimental data and their fit to the kinetics models are shown in Figure 6. For both initial concentrations (25 and 50 mg/L), equilibrium seems to be reached slowly, around 24 h of contact, and the experimental data are closer to the pseudo-second-order model (blue curve) (Figure 6a,b). The results of the fit to the kinetic models and the corresponding parameter values are shown in Table 2.
For CiCR = 25 mg/L, the three models presented a good fit to the experimental data with values of R2 > 0.900. PFO, PSO, and Elovich models show values of R2 = 0.957, R2 = 0.988, and R2 = 0.928, respectively. However, error functions show better results for PSO (Table 3). Additionally, the adsorption capacity at equilibrium (qe,calc) obtained with the pseudo-second-order model was 15.09 mg/g and is closer to the experimental adsorption capacity at equilibrium (qe,exp = 15.21 mg/g) than qe obtained with the pseudo-first-order model (13.87 mg/g). Therefore, the best model to describe the kinetics of the adsorption of Congo red at CiCR = 25 mg/L is the pseudo-second-order.
Figure 6b shows a good fit to the PSO and Elovich adsorption kinetic models at an initial concentration (CiCR = 50 mg/L) with values for R2 > 0.900 and small values for the error functions (Table 3). The adsorption capacity at equilibrium (qe,cal) obtained with the pseudo-second-order (PSO) model was 23.81 mg/g and is closer to qe,exp = 23.74 mg/g than qe,cal obtained with the pseudo-first-order PFO model. Furthermore, it is important to mention the good fit with the Elovich model (R2 = 0.968), which assumes that the activation energy increases with adsorption time and that the surface of the adsorbent is heterogeneous [39]. In summary, the adsorption kinetics of CR dye at the two initial concentrations (25 and 50 mg/L) follow pseudo-second-order (PSO). A tendency to fit the pseudo-second-order model at low initial dye concentrations is observed. Sabarinathan et al. (2019) [40] reported the same trend in the adsorption kinetics of methylene blue onto molecular polyoxometalate, and Ho & Mckay (1998) [41] studied the adsorption of two dyes, Basic Blue 69 and Acid Blue 25, onto peat. According to the results of the kinetic study, the capacity of adsorption at equilibrium (qe,exp) does not increase considerably between 24 and 48 h. Hence, the next equilibrium experiments were carried out at the equilibrium time of 24 h.

2.3. Adsorption Isotherms of CR at Different Solution pH

Congo red is used as a pH indicator. At acidic pH values (2.18–3.16), CR solutions turn blue, and there is a shift in the maximum absorbance wavelength [42], and at pH ≥ 3.86, there is no significant shift in the maximum absorbance wavelength. Hence, the study of the effect of pH was carried out in the range of 4 to 10 (Figure 7). According to the best values of the determination coefficient R2 (Table 4) and error functions (Table 5), the Redlich–Peterson and Langmuir models fit best over the entire studied pH range. β values from the Redlich–Peterson model are close to 1. Therefore, the fit to the Langmuir model is further supported, and it can be assumed that the adsorption of Congo red onto activated biochar prepared from Haematoxylum campechianum waste (ABHC) implies the following considerations: all the adsorption sites are equal; each adsorption site only keeps one molecule of the adsorbate and all sites do not sterically and energetically depend on the adsorbed amount of the adsorbate [43]. Table 4 shows the maximum adsorption capacity, Qmax = 53.93, 114.8, 68.33, 30.83, and 10.51 mg/g at pH = 4.0, 5.4, 7.0,8.4, and 10.2, respectively. At high pH values (7.0, 8.4, and 10.2), Qmax values decrease. This behavior can be attributed to the pKa value of Congo red = 4.5 [44], which means that CR would be negatively charged at basic values of pH solution and also at the pHpzc of ABHC, where, at pH > 6.5, the surface of ABHC is negatively charged. Therefore, there will be repulsion between the adsorbent surface’s negative charge and the adsorbate’s negative charge. A decrease in Qmax, from 114.8 to 53.93 mg/g at pH from 4.0 to 5.4 was also observed. This decrease is due to the positive charge of ABHC (pH < pHpzc = 6.5) and the protonation of CR, resulting in repulsion. The highest value of Qmax= 114.831 mg/g was obtained at a pH solution of 5.4, suggesting the best adsorption of CR at this pH value. Other works have reported an optimal pH value for the removal of CR around 5 [27,45,46].

2.4. Adsorption Isotherms of CR at Different Adsorbent Dose

Figure 8 depicts the effect of the dose of ABHC on the adsorption of Congo red dye. It is observed that for all doses (1 g/L to 10 g/L), the experimental data are close to Freundlich, Langmuir, and Redlich–Peterson nonlinear isotherms; however, for low doses (1 g/L and 2 g/L), the experimental data are also close to mode. In addition, a decrease in the adsorption capacity at equilibrium (qe) is observed when the dose of activated biochar increases (Figure 8a–d).
The parameters of each adjusted nonlinear isotherm model at a different dose of ABHC are listed in Table 6. At low dose values (1 g/L and 2 g/L), the values of R2 were superior to 0.900 for all isotherm models. For Langmuir isotherm, the R2 values were 0.940 and 0.988, while for Freundlich, R2 = 0.976 and 0.993 for doses of 1 g/L and 2 g/L, respectively. The good fit to both models suggests that the adsorption of CR dye onto activated biochar involves several interactions, which had already been mentioned previously in FTIR studies. At high doses of ABHC (5 g/L and 10 g/L), the best adjustments are obtained with the Langmuir (R2 = 0.970 and R2 = 0.966), Freundlich (R2 = 0.985 and R2 = 0.982), and Redlich–Peterson (R2 = 0.985 and R2 = 0.982).
In the whole range of doses (1 g/L to 10 g/L), the best-fitted model was Redlich–Peterson with the highest values of R2 and the lowest values of error functions (Table 7); this is a hybrid model of the Langmuir and Freundlich models, and when its dimensionless parameter β is close to one, the model approaches Langmuir, and when β is close to zero, the model approaches Freundlich [47]. At low dose, dose = 1 g/L, β = 0.742 and at dose = 2 g/L, β = 0.669 values are close to 1, so adsorption is close to the features of the Langmuir model. By increasing the dose to 5 g/L (β = 0.367 is closer to zero), adsorption is close to the features of the Freundlich model. Finally, Table 6 shows that Qmax (from the Langmuir isotherm model) decreases as the dose of activated biochar increases (from Qmax = 92.86 mg/g at dose = 1 g/L to Qmax = 12.22 mg/g at dose = 10 g/L); this behavior is due to the fact that as the dose of ABHC increases, there is a decrease in the amount of free active sites because of the agglomeration of sorption sites at high dose [9] or the solute concentration gradient between the solution of dye and the adsorbent surface, so the amount of CR adsorbed per unit weight of ABHC is reduced with increasing ABHC dose [48]. The highest value of the parameters Qmax= 92.86 mg/g and n = 3.079 were achieved with the lowest dose used in this work; therefore, the optimum dose of ABHC is 1 g/L.

2.5. Adsorption Isotherms of CR at Different Temperatures

The effect of the temperature on the adsorption of Congo red is shown in Figure 9. A change in the shape of the isotherm is observed when temperature increases from T = 300.15 to T = 330.15. In the whole studied temperature range, experimental data best fit the Redlich–Peterson model.
At low temperatures (Figure 9a), experimental data are closer to the Freundlich nonlinear model than the Langmuir nonlinear model. However, at higher temperatures (Figure 9c), experimental data and the shape of the isotherm are closer to the nonlinear Langmuir model than the Freundlich nonlinear model.
Table 8 summarizes the parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear isotherms. The maximum adsorption capacity, Qmax (mg/g) (from the Langmuir model), decreased significantly from 92.86 to 29.20 mg/g at 300.15 and 330.15 K, respectively. In other words, the adsorption of CR onto ABHC is less favorable at high temperatures. In addition, KF (from the Freundlich model) values decreased when increasing the temperature from 300.15 to 330.15 K, which also indicates unfavorable adsorption at high temperatures. This behavior has been commonly reported in the adsorption of azo dyes, suggesting a better adsorption of these types of dyes at low temperatures [49].
According to R2 and error function values (Table 9), at 300.15 and 313.15 K, the best fit to the isotherm models corresponds to the following order: Redlich–Peterson > Freundlich >Langmuir. While, at 330.15 K, Redlich–Peterson > Langmuir> Freundlich. That is, the best isotherm model to describe the adsorption of Congo red dye onto ABHC in the full range of studied temperatures is the Redlich–Peterson, which assumes the considerations of Freundlich and Langmuir models [47]. However, values of the parameter β presented an increasing trend as the temperature increased (β = 0.742 at T = 300.15 K and β = 0.899 at T = 330.15 K). Therefore, the dye adsorption tends to be in the Langmuir model at high temperatures.

2.6. Comparison of the Maximum Adsorption Capacity Qmax onto Different Activated Carbons

As mentioned in the introduction, CR dye is difficult to remove from the aqueous medium; this proves the low maximum adsorption capacity obtained with other activated carbons (Table 10). However, [50,51] and this work have reported high values of Qmax (Qmax > 100 mg/g), which indicates these activated carbons are good sorbents to remove CR. In Table 10, the experimental conditions are also shown. It is observed that most of the works reported Qmax at temperatures that are not too high, in the range of 25–30 °C, like this work. With respect to the pH value, high pH values are not observed, i.e., high adsorption capacities are not obtained at very high pH (as mentioned previously in the pH effect). However, Qmax values reported at lower pH values (2–3) are sometimes questioned due to the change in the wavelength of maximum absorbance of Congo red at acidic pH values that can lead to errors at the moment of quantifying CR solutions [42]. Activated biochar from Haematoxylum campechianum (ABHC) presented a high adsorption capacity of CR dye at not-so-extreme experimental conditions compared to the other activated carbons.

2.7. Adsorption and Desorption Cycles of CR

Two different NaOH concentrations (0.01 and 0.2 M) were used in two adsorption/desorption cycles. For 0.01M NaOH (Figure 10a), Qads decreased from 44.70 mg/g to 5.50 mg/g for cycles 1 and 2, respectively. Regarding the Qdes, after cycle number 1, only 8.58 mg/g of CR was desorbed, and after cycle number 2, no desorption of CR was achieved. While, for 0.2M NaOH (Figure 10b), the adsorption capacity (Qads) resulted in 7.69 mg/g for cycle number 2; comparing these results with those mentioned above for 0.01M NaOH, there is only an increase of approximately 2 mg/g in the adsorption capacity using 0.2M NaOH. Whereas, for the desorption capacity, after cycle number 1, Qdes increased from 8.58 mg/g to 18.01 mg/g, and for cycle 2, a desorption capacity of 3.02 mg/g was achieved, suggesting better desorption of Congo red with increasing NaOH concentration. However, these desorption results are low, which means that CR is well sorbed onto the surface of the ABHC and that the interactions are strong. The low adsorption removal after two cycles indicates the low possibility of reuse, which is a disadvantage of activated carbon that has already been reported [57]. In order to improve the reusability of ABHC, it is suggested to apply any other regeneration method such as another eluent or even temperature.

2.8. Machine Learning-Assisted Optimization for Improving the Adsorption Process

Machine learning (ML) techniques have become pivotal in enhancing the predictability [58] and understanding of adsorption processes across various domains, including environmental and materials science. This review synthesizes recent advancements in the integration of ML with adsorption modeling, highlighting its application in predicting the behavior of diverse adsorbents and adsorbates, from heavy metals to organic pollutants and gases. Guanwei Yin et al. (2021) explored different ML models—Multi-Layer Perceptron (MLP), Passive Aggressive Regression, and Decision Tree Regressor—to predict dye adsorption from aqueous solutions. Among these, the Decision Tree Regressor emerged as the most effective, suggesting its suitability for correlating adsorption equilibrium data due to high accuracy and low error rates (R2 = 0.99) [59]. Chen Zhao et al. (2024) demonstrated how multivalent ions like Ca2+, K+, Na+, and Mg2+ impact the adsorption of azo dyes using ML models coupled with Density Functional Theory (DFT). Their study revealed that the Gradient Boosting Decision Tree model provided the best fit, suggesting that mixed background ions significantly influence adsorption behavior, as validated through experimental and theoretical analyses [60]. Shoushi Zhao et al. (2024) applied ML to predict adsorption of various metal cations by clay minerals, identifying key descriptors through extensive feature engineering. The Extreme Gradient Boosting model was highlighted for its exceptional predictive accuracy (R2 = 0.977), underscoring the model’s capability in geochemical applications [61]. Qing-Yun Cai et al. (2024) focused on predicting protein adsorption capacities using QSAR models derived from ML. Utilizing Random Forest and Gradient Boosting methods, they achieved high predictive performances across various datasets, showing that ML can significantly contribute to the design of more efficient separation processes in biochemical engineering [62]. Hyeonmin Lee and Yongju Choi (2023) developed ML models to predict the adsorption capacity of sediment-amended activated carbon for hydrophobic organic contaminants. The top-performing model achieved an R2 of 0.94, illustrating the utility of ML in enhancing the predictability of adsorption processes in complex environmental matrices [63]. Lisheng Guo et al. (2024) employed multiple ML algorithms to model heavy metal adsorption by bentonite, with the extreme Gradient Boosting Regression (XGB) model showing the best performance. This study also provided insights into the influence of various factors on adsorption capacity, supported by a web-based graphical interface for model access [64]. Kai Chen et al. (2023) developed ML models to predict the adsorption percentage of oxyanions on goethite, assessing the influence of specific surface area and other descriptors. The study provided detailed insights into the non-linear relationships and crystal face compositions affecting adsorption capacities [65]. Raja Selvaraj et al. (2024) utilized ANN and ANFIS to model tetracycline adsorption on activated carbon derived from fruit biomass. This approach not only yielded high correlation coefficients but also provided a detailed mechanism of adsorption through statistical physics models, highlighting the method’s applicability in pharmaceutical waste treatment [66]. The reviewed studies collectively underscore the transformative potential of machine learning in adsorption science. By effectively predicting adsorption capacities and understanding the underlying mechanisms through models like Decision Trees, Gradient Boosting, and Extreme Gradient Boosting, ML facilitates more precise control and optimization of adsorption processes across varied applications. This integration not only enhances model accuracy but also contributes significantly to environmental sustainability and material efficiency.
The scatter matrix provided in Figure 11 presents a visual representation of pairwise relationships across multiple variables, with dot colors corresponding to different ranges of the “Removal percentage”. From the distribution of points, it can discern patterns and correlations within the data. The descriptive statistics complement the scatter matrix by quantifying the central tendency and dispersion of the same variables. With a count of 105 for each variable, it is a sizeable dataset. The “Initial concentration” has a mean of approximately 201.55, but its standard deviation is relatively high, at 293.31, indicating a wide range of values, which is also visible in the scatter plot with a spread across the axis. The “Time” variable shows less variation relative to its mean, suggesting more consistent data. The “Temperature” data are less varied, with a small standard deviation around its mean of 303.02. The “pH” values have a mean near neutral, at 7.01, and “Dosis” averages around 1.93 with a larger standard deviation, showing some variability in the dataset. The “Removal percentage” has a mean of 35.71, with a broad standard deviation of 22.47, which indicates variability in the removal efficiency. In the scatter matrix, there does not seem to be a clear linear pattern or strong correlation between most variables and the “Removal percentage”, as indicated by the wide scatter of points. However, a higher concentration of points along the lower range of “Initial concentration” suggests that lower concentrations may be associated with a broader range of “Removal percentage” outcomes. “Time” seems to show clusters at specific intervals, perhaps indicative of experimental time points or specific durations at which data were collected. Temperature appears fairly constant except for a few higher measurements, which might represent a different experimental condition. The pH levels are centered around neutral, with fewer instances of more acidic or basic conditions. The dosages show a strong clustering at the lower end, suggesting most experiments used a smaller dose, with fewer instances using higher dosages. Overall, the combination of scatter matrix and descriptive statistics provides a comprehensive overview of the dataset, indicating a wide variance in some variables and suggesting potential trends that might warrant further investigation, such as the influence of “Initial concentration” on the “Removal percentage”. These insights could be crucial for optimizing conditions to enhance the removal process’s efficiency.
The heatmap visualization of this matrix, as constructed with seaborn and matplotlib, implementing Pearson Coefficient of Correlation, enhance the comprehensibility of these relationships, presenting a clear and concise graphical representation of the data, where the ‘coolwarm’ color palette effectively differentiates between positive and negative associations, displayed in Figure 12. In examining the correlation matrix for the dataset, there is a notable linear relationship between the variables and the removal percentage, which serves as the target variable in the context of the removal process’s efficacy. The initial concentration is significantly inversely correlated with the removal percentage (r = −0.61), indicating a tendency for the removal efficiency to decrease as the initial concentration of the substance increases. On the other hand, time exhibits a positive but weak correlation with the removal percentage (r = 0.14), suggesting only a minor influence on removal efficiency. The temperature’s negative and moderate correlation (r = −0.28) with removal percentage may reflect a thermally dependent process where increased temperatures could potentially impede the removal efficiency. In the case of pH, there exists a moderately negative correlation with removal percentage (r = −0.40), hinting at a higher pH potentially reducing the efficacy of removal. Dosis, however, shows a negligible negative correlation with removal percentage (r = −0.06), indicating a very slight decrease in removal efficiency with increased dosis, though this relationship is so weak that it might not be significant. This correlation matrix not only highlights the strength and direction of relationships between multiple variables but also underscores the intricate dynamics within the removal process. The negative correlations observed with initial concentration and pH notably suggest that these variables are crucial for understanding the process in depth. While correlations provide insights into potential patterns within the dataset, they do not confirm causation and should prompt further investigative analysis, such as regression modeling, to unravel the complex interactions between these factors and optimize the conditions for maximum removal efficiency.
For the optimization [67] of removal percentage, predictive models of removal percentage, using a variety of machine learning algorithms, were evaluated for their effectiveness. The Gradient Boosting algorithm outshined its counterparts, exhibiting superior performance in terms of accuracy and computational efficiency. Henceforth, this discussion delineates the development process of a Gradient Boosting regression model, which is presented as the most efficacious algorithm for establishing a regression model between the input variables—’Initial concentration’, ‘Time’, ‘Temperature’, ‘pH’, and ‘Dosis’—and the target variable, ‘Removal percentage’. The Gradient Boosting Regressor was constructed using Python’s renowned Scikit-learn library, a tool of choice for its comprehensive suite of machine learning algorithms. The model was tailored with hyperparameters set to 100 estimators, ensuring a robust learning process without overfitting. The data division was configured with a traditional 80–20 split between the training and testing datasets, facilitated by Scikit-learn’s train_test_split function. Random state control was employed to ensure reproducibility of the results. This computational experiment was executed within a Python programming environment on Google Colab, renowned for its extensive support for data analysis and machine learning. The meticulous selection and tuning of hyperparameters were pivotal in harnessing the full potential of the Gradient Boosting method. The resultant model not only encapsulates the intricate relationships within the data but also stands as a testament to the meticulous computational methodology, paving the way for future endeavors in process optimization.
Figure 13 shows a scatter plot that illustrates the relationship between the actual and predicted removal percentages derived from a Gradient Boosting regression model. The green dots represent predictions on the training dataset, while the red dots denote predictions on the testing dataset. A dashed blue line, described as the “fit line”, suggests a linear relationship indicative of the model’s predictions compared to the true values. Generally, the proximity of the dots to the fit line reflects the accuracy of the model: the closer the dots are to the line, the more accurate the predictions. In this case, both green and red dots are close to the fit line, indicating good model performance on both the training and testing sets. The continuation of the dots along the fit line without deviating significantly in either direction suggests that the model does not have a systemic bias toward over- or under-prediction. The blue fit line likely represents the ideal scenario where the predicted values exactly match the actual values, often referred to as the line of perfect fit. If the model’s predictions were perfect, all points would lie on this line, with predicted values equaling actual values. The training and testing data are clustered around a trend line defined by the equation y = 0.9583x + 1.5813. This equation embodies the model’s predictive behavior, where y represents the predicted removal percentage and x represents the actual removal percentage. The slope of the trend line, 0.9583, indicates that for every unit increase in the actual removal percentage, the predicted removal percentage increases by approximately 0.9583 units. This value, being close to 1, signifies that the model predictions are nearly proportional to the actual values, highlighting the model’s accuracy. The y-intercept of 1.5813 suggests that when the actual removal percentage is zero, the model predicts a baseline removal percentage of 1.5813, which might be interpreted as the model’s expected performance when no contributing factors are present. Overall, the proximity of the data points to the trend line, along with the slope nearing unity, reinforces the model’s precision and suggests minimal bias in prediction. However, the slight deviation from a perfect 1:1 slope implies that the model may slightly underestimate the removal percentage since the slope is less than 1. The y-intercept also suggests a small systematic offset in the model’s predictions. This detailed interpretation of the trend line equation allows for a nuanced understanding of the model’s characteristics and provides an avenue for further refining its predictions for use in optimization processes.
The Gradient Boosting (GB) model’s performance metrics for both training and testing phases offer a comprehensive evaluation of its predictive accuracy, as displayed in Table 11. The mean squared error (MSE) for the training set stands at approximately 4.05, suggesting a satisfactory model fit, whereas a notably higher MSE of around 33.17 for the testing set indicates potential overfitting, with the model not generalizing as well to unseen data. The sum of squares error (SSE) corroborates this, with values of 340.39 and 696.54 for the training and testing sets, respectively, indicating a greater cumulative deviation in the testing phase. The mean absolute percentage error (MAPE) values, 5.95% for training and 23.09% for testing, further reflect the model’s accuracy deterioration from training to testing. Root mean squared error (RMSE) follows the same pattern, with a lower training error of 2.01, contrasted with a higher testing error of 5.76, mirroring the MSE’s suggestion of overfitting. Similarly, the mean absolute error (MAE) is relatively low for training data at 1.49 but escalates to 3.70 for the testing data, highlighting larger average prediction errors when the model encounters new data. The mean percentage error (MPE) indicates an underprediction bias in both datasets, with more pronounced underestimation in the testing data at −14.27%, compared to −0.92% for training data. The coefficient of determination (COD), or R², which assesses the variance in the dependent variable explained by the independent variables, presents high values of 0.992 for training and 0.914 for testing. Despite the reduction in the testing phase, these values imply the model’s strong explanatory power. Overall, while the GB model exhibits robust predictive capabilities for the training data, evidenced by high R² and low error metrics, the escalation of error measures in the testing data suggests an imperative for model refinement to enhance its generalization and mitigate overfitting, ensuring more reliable predictions across new datasets.
Figure 14 presents a parallel coordinates plot, a multidimensional visualization that allows for the inspection of individual data points across several quantitative variables simultaneously. The axes represent the variables under consideration, namely Initial Concentration, Time, Temperature, pH, Dosis, and Removal Percentage. Each line in the plot corresponds to a single observation from the dataset, with the position on each axis reflecting its value for the respective variable. Color gradation from turquoise to deep red indicates the spectrum of Removal Percentage values, serving as the classification criterion in this context. Observations with higher removal efficiencies are denoted in warmer hues, which suggests a tendency towards lower Initial Concentrations and possibly lower Dosis levels. Notably, there is no discernible trend relating Temperature or Time to removal efficiency, as lines are largely intermixed along these axes. Meanwhile, pH values predominantly cluster around neutral, with only a few observations indicating more extreme acidic or basic conditions. This visualization technique reveals complex relationships within the dataset that may not be apparent through traditional scatter plots or correlation matrices. It can be instrumental in identifying patterns, outliers, and potential areas of interest for further statistical analysis or experimental investigation. For instance, one could hypothesize from the visualization that for this particular removal process, optimal efficiency may be more sensitive to variations in initial concentrations and less affected by temperature or dosis variations within the examined range. Such insights are valuable for optimizing the parameters of the removal process to achieve higher efficiencies.
The problem is then solved for the optimization problem [68,69]. For the maximization of the removal percentage, the developed code implements a process known as Bayesian optimization using the scikit-optimize library to fine-tune the hyperparameters of a Gradient Boosting Regressor. The aim is to find the set of hyperparameters that maximizes the removal percentage predicted by the model.
Hyperparameters: These are the adjustable parameters that control the model training process. For Gradient Boosting, typical hyperparameters can include the number of trees (n_estimators), the learning rate, the depth of each tree, etc. The specific hyperparameters being optimized in this snippet are not directly shown, but we assume that they are encapsulated within the GradientBoostingRegressor with n_estimators = 100 and random_state = 42. The n_estimators parameter determines the number of sequential trees to be modeled, and the random_state ensures reproducibility of results.
Search Space: The search space is defined by a set of ranges for the input variables the model will use to make predictions. This space is represented as a list of Real or Integer objects, which are types provided by scikit-optimize to specify continuous and discrete parameters, respectively.
  • ‘Initial concentration’ varies between 10 and 1000 (continuous).
  • ‘Time’ varies between 15 and 2880 (discrete).
  • ‘Temperature’ varies between 300.15 and 330.15 (continuous).
  • ‘pH’ varies between 4 and 10 (continuous).
  • ‘Dosis’ varies between 1 and 10 (continuous).
These parameters are part of the input feature set for the Gradient Boosting model, and the ranges provide boundaries for the Bayesian optimization algorithm to search within.
Optimization Algorithm: The following specific parameters are used for the Bayesian optimization algorithm.
  • gp_minimize is the function from scikit-optimize that performs the Bayesian optimization. It utilizes Gaussian Processes to model the probability distribution of the objective function and makes educated guesses where the function might achieve optimal values.
  • The objective function is defined to return the negative predicted removal percentage by the model given a set of parameters. The negative sign is used because gp_minimize by default searches for the minimum value of the function, but since we want to maximize the removal percentage, we need to minimize its negative.
  • n_calls = 50 specifies the number of evaluations of the objective function, or how many times the algorithm will try different sets of parameters.
  • random_state = 42 ensures that the results are reproducible; the algorithm will start from the same random seed.
The convergence plot, which is shown in Figure 15, visualizes the optimization algorithm’s progress over successive iterations, depicting the minimum objective function value attained after each function call. In the context of Bayesian optimization, each function call equates to an evaluation of the objective function with a particular set of hyperparameters. The y-axis, representing the minimum of the objective function, shows an initially steep decline, indicating a rapid improvement in the model’s performance with the first few evaluations. As the optimization proceeds, the curve flattens out, suggesting that subsequent iterations provide incremental improvements and eventually reach a plateau. This behavior is characteristic of optimization algorithms converging towards an optimal set of parameters. Around the 20th call, the curve levels significantly, implying that the algorithm has likely identified a region close to the optimum. Subsequent iterations refine the solution, but the objective function’s value changes only marginally, signifying that further searches yield little benefit, and the algorithm is approaching the maximum predicted removal percentage. This plateau indicates the optimal solution space within the predefined hyperparameter boundaries has been thoroughly explored, and additional evaluations do not substantially enhance the predictive capacity of the model. The figure encapsulates the efficiency and effectiveness of the Bayesian optimization process, demonstrating a pronounced convergence within 50 evaluations, which is critical in the context of resource-intensive computational tasks typical in scientific research.
Table 12 presents the set of input variables that yield the optimal removal percentage as determined by the optimization algorithm. The table outlines that an initial concentration of 10.0 units, a lengthy process time of 2880 units, a mid-range temperature of 313.15 units, a relatively acidic pH of 4.0, and the lowest dosis level of 1.0 unit culminate in a removal percentage of approximately 90.4733%. This optimal scenario likely reflects a specific operational sweet spot for the system under study, balancing the chemical and physical inputs to maximize efficiency. The removal percentage achieved underscores the model’s capability to identify conditions under which the process efficiency is near its peak within the constraints and ranges specified in the study.

3. Materials and Methods

3.1. Adsorbate

Congo red dye (C.I. 22120, CAS number: 573-58-0, chemical formula: C32H22N6Na2O6S2) was reagent grade (99.92%) and procured from Fagalab. Figure 16 shows the CR molecule (drawn with ChemDraw professional software 17 version 17.0.0.206), and Table 13 presents the physical properties of CR dye.

3.2. Biochar Preparation

The bark of Haematoxylum campechianum was collected from the Faculty of Engineering, Autonomous University of Carmen in Campeche, Mexico. First, a pretreatment was carried out. The bark was chopped, milled, and sieved between 0.2 and 0.5 mm. Then, it was washed with abundant distilled water at 50 °C to eliminate any residue from the surface and oven-dried at 70 °C for 12 h.
For the chemical activation, 250 mL of H3PO4 was mixed with 50 g of Haematoxylum campechianum (pretreated, previously) for 3 h at 50 °C, and then, the mixture was filtered and the solids dried at 70 °C for 12 h. For the thermal treatment, 50 g of phosphoric acid-treated Haematoxylum campechianum was introduced into a muffle at 500 °C for 60 min, at 10 °C/min. Then, when the sample was cold, in order to remove H3PO4 residue, the activated biochar was washed with a 5% NaHCO3 solution and then with abundant distilled water until the filtrate reached a pH value between 6 and 7. Finally, activated biochar from Haematoxylum campechianum (ABHC) was dried at 110 °C for 12 h and then stored in a closed glass bottle and placed in a desiccator.

3.3. Characterization Techniques

Scanning electron microscope (SEM) coupled with EDS system (HITACHI S-3400N), Brunauer Emmett Teller method (BET) (BELSORP MAX), X-ray diffractometer (XRD) (APD 2000 PRO), and Fourier-transform infrared (FTIR) spectrophotometer (Nicolet Nexus 670) techniques were used to characterize the surface morphology, shape, elemental composition, mean pore diameter, surface area, crystalline structure, and chemical composition of the activated biochar from Haematoxylum campechianum (ABHC). Additionally, the pH point of zero charge pHpzc was determined using 0.01 M NaCl as the background electrolyte solution; the pH of the working solutions was fixed between 1 and 12 (2, 4, 5, 6, 8, 10, and 12) using 0.1 M NaOH and/or HCl. First, the initial pH of each solution was measured with a Hanna instruments (model HI2020-01) pH meter. Subsequently, solutions were put in contact with the adsorbent at 120 rpm at 300.15 K for 24 h. Then, the solutions were decanted, and the final pH values were measured. The point of intersection between the curve was obtained by plotting pHfinal vs. pHinitial, and the diagonal curve corresponds to the pHpzc of the adsorbent [70].

3.4. Sorption

Sorption of Congo red is affected by several parameters, including the pH of the aqueous medium, temperature, adsorbent dose, adsorbate concentration, and contact time.
Different pH values (4.0, 5.4, 7.0, 8.4, and 10.2) of CR solutions were adjusted by adding 0.1 M HCl and 0.1 M NaOH, using a pH sensor (model HI2020-01, Hanna instruments) at different initial concentrations (10, 20, 40, 50, 60, 80, and 100 mg/L) and put in contact with a dose of activated biochar =1 g/L in a rotating incubator at 150 rpm for 24 h at 300.15K.
The effect of the temperature was evaluated by putting in contact with different initial concentration CR solutions (20, 50, 100, 250, 500, 750, and 1000 mg/L) with a dose of ABHC =1g/L, pH solution = 7.0 at three different temperatures: 300.15K, 313.15K, and 330.15K for 6 h.
In order to evaluate the effect of the dose of adsorbent material on the dye adsorption, four doses of ABHC (1, 2, 5, and 10 g/L) were put in contact for 24 h with different initial concentration CR solutions (20, 50, 100, 250, 500, 750, and 1000 mg/L) at pH solution = 7.0 and temperature = 300.15 K.
After reaching equilibrium, all samples were centrifugated for 5 min at 3500 rpm and then analyzed by a UV-vis spectrophotometer (model EVOLUTION 220, Thermo Scientific) using calibration curves (with an interval from 2 to 20 mg/L) and R2 > 0.995 at the wavelength of 497 nm.
The adsorption capacity (qe) was calculated using the following Equation:
q e = C i C e · V m
where Ci (mg/L) is the initial concentration, Ce (mg/L) is the equilibrium concentration, qe (mg/g) is the adsorption capacity, V (L) is the volume, and m (g) is the mass of the adsorbent.
Experimental equilibrium data were fitted with the following nonlinear isotherms models: Langmuir (Equation (3)), Freundlich (Equation (4)), and Redlich–Peterson (Equation (5)), using the solver add-in from Microsoft’s spreadsheet tool of Microsoft Excel for Office 365 version 2403.
q e = Q m a x K L C e 1 + K L C e
where KL (L/mg) is the Langmuir equilibrium constant and Qmax (mg/g) is the maximum adsorption capacity.
q e = K F C e 1 / n
where n (dimensionless) and KF (mg1−1/n·L1/n·g−1) are the exponent and the Freundlich parameter, respectively.
q e = K R P C e 1 + α R P C e β
where KRP (L/g) and αRP (L/mg) are the constant of the Redlich–Peterson and β is a dimensionless parameter between 0 and 1.
Kinetic studies were conducted from two different initial concentrations: 25 mg/L and 50 mg/L at pH = 7.0, an adsorbent dose = 1 g/L, temperature = 300.15 K, and the following contact times: 15, 60, 180, 360, 720, 1440, and 2880 min. The obtained experimental data were fitted with the pseudo-first-order (PFO), pseudo-second-order (PSO), and Elovich kinetic nonlinear models (Equation (6), Equation (7), and Equation (8), respectively) using the solver add-in from Microsoft’s spreadsheet tool of Microsoft Excel.
q t = q e , c a l 1 e x p ( k 1 t )
q t = k 2 q e , c a l 2 t 1 + k 2 q e , c a l t
where qe,cal is the theoretical adsorption capacity (mg/g) and k1 (1/min) and k2 [g/(mg·min)] are the pseudo-first and pseudo-second-order rate constants, respectively.
q t = 1 β l n ( α β t )
where β is the constant related to the extent of surface coverage (g/mg) and α is the theoretical adsorption capacity [mg/(g·min)].
In addition to the regression coefficient R2, the following error functions were calculated with the aim of confirming accurate measurement results and a good fit with the proposed models: absolute error (EABS), sum squared error (SSE), the average percentage error (ARE), nonlinear chi-square test (χ2), residual root mean square error (RMSE), normalized standard deviation (Δq (%)) (Equation (9), Equation (10), Equation (11), Equation (12), Equation (13) and Equation (14), respectively).
S S E = i = 1 N ( q e , e x p q e , c a l c ) 2
E A B S = 1 n i = 1 N ( q e , e x p q e , c a l c
A R E % = 100 N i = 1 N q e , e x p q e , c a l c q e , e x p
χ 2 = i = 1 N ( q e , e x p q e , c a l c ) 2 q e , c a l c
R M S E = 1 N 2 i = 1 N ( q e , e x p q e , c a l c ) 2
Δ q ( % ) = 100 1 N 1 i = 1 N ( q e , e x p q e , c a l c q e , e x p ) 2
where qe,exp is the adsorption capacity obtained from the batch experiment (mg/g), qe,calc is the adsorption capacity obtained with a mathematical model corresponding (mg/g), and N is the corresponding number of observations in the experiment.

3.5. Desorption and Regeneration Study

In order to determine the desorption behavior and reusability of ABHC, 0.01 M and 0.2 M NaOH solutions were used as desorbing agents. A total of 50 mL of CR at the initial concentration of 100 mg/L was put in contact with 50 mg of ABHC for 24 h. Then, CR-loaded ABHC was regenerated using 50 mL of NaOH solution for CR desorption. Regeneration was carried out at 300.15 for 24 h. After a time, the sample was centrifuged to obtain the adsorbent. The residual concentration was measured by UV-vis spectrophotometric analysis. The regenerated ABHC was rinsed with distilled water and dried at 105 °C for 2 h before the next cycle of adsorption took place. Two cycles of consecutive adsorption–desorption studies were carried out using two different concentrations of NaOH (0.01 M and 0.2 M).
Qdes (mg/g) is the quantity of dye desorbed that was calculated using Equation (15) [71].
Q d e s = C d e s · V m
where Cdes (mg/L) is the concentration of dye left after desorption, V (L) is the volume of the dye solution, and m (g) is the ABHC mass.
Qads (mg/g) is the quantity of dye adsorbed after sorption that was calculated using Equation (2).

4. Conclusions

According to the characterization, isotherm, and kinetic results, the interactions involved in the CR removal process onto biochar prepared from Haematoxylum campechianum waste (ABHC) were electrostatic attraction, hydrophobic attraction, hydrogen bonding, and n–π and π–π interactions. The results of the evaluated effects (pH, temperature, and dose) give a maximum adsorption capacity (Qmax) equal to 114.8 mg/g at pH = 5.4, 300.15K, and a dose of ABHC of 1 g/L.
Advanced machine learning techniques including Gradient Boosting Model combined with an optimization algorithm were applied to optimize a chemical removal process. The resulting model predicted a removal percentage of approximately 90.47% under optimal conditions: an initial concentration of 10 units, time of 2880 units, temperature at 313.15 units, pH of 4.0, and dosis of 1.0.
Finally, this study underscores the vital role of computational thinking in addressing environmental challenges through the implementation of computational packages such as data science, machine learning, and optimization. This work not only contributes to the field of environmental science but also provides valuable insights for integrating computational thinking skills into educational curricula. We can foster innovation and sustainability, empowering individuals to tackle real-world environmental problems with confidence and efficacy.

Author Contributions

D.M.P.G. and M.A. conceived and planned the experiments. D.M.P.G., M.A., and R.T. wrote the manuscript. R.T. also conducted all the data analytics, machine learning, and optimization of the process. J.V. contributed to the characterization of the samples. E.L., F.A.F., C.A.U. and M.A.R.E. contributed to the interpretation of the results. All authors provided critical feedback and helped shape the research, analysis, and manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank the financial support from Tecnológico de Monterrey through the “Challenge-Based Research Funding Program 2022”. Project ID # I003-IFE001-C2-T3-T.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We acknowledge the supports of the Faculty of Engineering and Chemistry, Universidad Autónoma del Carmen, Campeche, Mexico.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kumari, D.; Mazumder, P.; Kumar, M.; Deka, J.P.; Shim, J. Simultaneous Removal of Cong Red and Cr (VI) in Aqueous Solution by Using Mn Powder Extracted from Battery Waste Solution. Groundw. Sustain. Dev. 2018, 7, 459–464. [Google Scholar] [CrossRef]
  2. Lakshmi, U.R.; Srivastava, V.C.; Mall, I.D.; Lataye, D.H. Rice Husk Ash as an Effective Adsorbent: Evaluation of Adsorptive Characteristics for Indigo Carmine Dye. J. Environ. Manag. 2009, 90, 710–720. [Google Scholar] [CrossRef] [PubMed]
  3. Gregory, P. Classification of Dyes by Chemical Structure. In The Chemistry and Application of Dyes. Topics in Applied Chemistry; Waring, D.R., Hallas, G., Eds.; Springer: Boston, MA, USA, 1990. [Google Scholar]
  4. Binupriya, A.R.; Sathishkumar, M.; Swaminathan, K.; Kuz, C.S.; Yun, S.E. Comparative Studies on Removal of Congo Red by Native and Modified Mycelial Pellets of Trametes Versicolor in Various Reactor Modes. Bioresour. Technol. 2008, 99, 1080–1088. [Google Scholar] [CrossRef] [PubMed]
  5. Wekoye, J.N.; Wanyonyi, W.C.; Wangila, P.T.; Tonui, M.K. Kinetic and Equilibrium Studies of Congo Red Dye Adsorption on Cabbage Waste Powder. Environ. Chem. Ecotoxicol. 2020, 2, 24–31. [Google Scholar] [CrossRef]
  6. Csillag, K.; Emri, T.; Rangel, D.E.N.; Pócsi, I. PH-Dependent Effect of Congo Red on the Growth of Aspergillus Nidulans and Aspergillus Niger. Fungal Biol. 2023, 127, 1180–1186. [Google Scholar] [CrossRef] [PubMed]
  7. Tarek, M.; Saleem, A. Methods for Staining Amyloid in Tissues. A Review. Stain Technol. 1988, 63, 201–212. [Google Scholar]
  8. Xue-Yun, Y.; Ji-Gang, C.; Yong-Ning, H. Studies on the Relation between Bladder Cancer and Benzidine or Its Derived Dyes in Shanghai. Br. J. Ind. Med. 1990, 47, 544–552. [Google Scholar]
  9. Leudjo Taka, A.; Fosso-Kankeu, E.; Pillay, K.; Yangkou Mbianda, X. Metal Nanoparticles Decorated Phosphorylated Carbon Nanotube/Cyclodextrin Nanosponge for Trichloroethylene and Congo Red Dye Adsorption from Wastewater. J. Environ. Chem. Eng. 2020, 8, 103602. [Google Scholar] [CrossRef]
  10. Olaru, S.; Andrei, L.; Hancu, E.; Costica, N.; Ivanescu, L.; Zamfirache, M. Toxicity Assessment of Congo Red Azo Dye towards Lemna minor L. 2016. Available online: https://www.researchgate.net/profile/Stefan-Olaru/publication/303626230_Toxicity_assessment_of_Congo_Red_azo_dye_towards_Lemna_minor_L/links/574a8e9008ae5c51e29e8f07/Toxicity-assessment-of-Congo-Red-azo-dye-towards-Lemna-minor-L.pdf (accessed on 8 April 2024).
  11. Hernández-Zamora, M.; Martínez-Jerónimo, F. Congo Red Dye Diversely Affects Organisms of Different Trophic Levels: A Comparative Study with Microalgae, Cladocerans, and Zebrafish Embryos. Environ. Sci. Pollut. Res. 2019, 26, 11743–11755. [Google Scholar] [CrossRef]
  12. Chatterjee, S.; Dey, S.; Sarma, M.; Chaudhuri, P.; Das, S. Biodegradation of Congo Red by Manglicolous Filamentous Fungus Aspergillus Flavus JKSC-7 Isolated from Indian Sundabaran Mangrove Ecosystem. Appl. Biochem. Microbiol. 2020, 56, 708–717. [Google Scholar] [CrossRef]
  13. Hamad, M.T.M.H.; Soliman, M.S.S. Application of Immobilized Aspergillus Niger in Alginate for Decolourization of Congo Red Dye by Using Kinetics Studies. J. Polym. Environ. 2020, 28, 3164–3180. [Google Scholar] [CrossRef]
  14. Ghany, T.; Abboud, M.; Alawlaqi, M.; Shater, A. Dead Biomass of Thermophilic Aspergillus Fumigatus for Congo Red Biosorption. Egypt. J. Exp. Biol. 2019, 15, 1–6. [Google Scholar] [CrossRef]
  15. Rafiq, A.; Ikram, M.; Ali, S.; Niaz, F.; Khan, M.; Khan, Q.; Maqbool, M. Photocatalytic Degradation of Dyes Using Semiconductor Photocatalysts to Clean Industrial Water Pollution. J. Ind. Eng. Chem. 2021, 97, 111–128. [Google Scholar] [CrossRef]
  16. Yashni, G.; Al-Gheethi, A.; Radin Mohamed, R.M.S.; Dai-Viet, N.V.; Al-Kahtani, A.A.; Al-Sahari, M.; Nor Hazhar, N.J.; Noman, E.; Alkhadher, S. Bio-Inspired ZnO NPs Synthesized from Citrus Sinensis Peels Extract for Congo Red Removal from Textile Wastewater via Photocatalysis: Optimization, Mechanisms, Techno-Economic Analysis. Chemosphere 2021, 281, 130661. [Google Scholar] [CrossRef]
  17. Pascariu, P.; Cojocaru, C.; Olaru, N.; Samoila, P.; Airinei, A.; Ignat, M.; Sacarescu, L.; Timpu, D. Novel Rare Earth (RE-La, Er, Sm) Metal Doped ZnO Photocatalysts for Degradation of Congo-Red Dye: Synthesis, Characterization and Kinetic Studies. J. Environ. Manag. 2019, 239, 225–234. [Google Scholar] [CrossRef] [PubMed]
  18. Devi, V.S.; Sudhakar, B.; Prasad, K.; Jeremiah Sunadh, P.; Krishna, M. Adsorption of Congo Red from Aqueous Solution onto Antigonon Leptopus Leaf Powder: Equilibrium and Kinetic Modeling. Mater. Today Proc. 2019, 26, 3197–3206. [Google Scholar] [CrossRef]
  19. Adelaja, O.A.; Bankole, A.C.; Oladipo, M.E.; Lene, D.B. Biosorption of Hg(II) Ions, Congo Red and Their Binary Mixture Using Raw and Chemically Activated Mango Leaves. Int. J. Energy Water Resour. 2019, 3, 1–12. [Google Scholar] [CrossRef]
  20. Jabar, J.M.; Odusote, Y.A.; Alabi, K.A.; Ahmed, I.B. Kinetics and Mechanisms of Congo-Red Dye Removal from Aqueous Solution Using Activated Moringa Oleifera Seed Coat as Adsorbent. Appl. Water Sci. 2020, 10, 136. [Google Scholar] [CrossRef]
  21. Lim, L.B.L.; Priyantha, N.; Latip, S.A.A.; Lu, Y.C.; Mahadi, A.H. Converting Hylocereus Undatus (White Dragon Fruit) Peel Waste into a Useful Potential Adsorbent for the Removal of Toxic Congo Red Dye. Desalination Water Treat. 2020, 185, 307–317. [Google Scholar] [CrossRef]
  22. Jabar, J.M.; Odusote, Y.A. Removal of Cibacron Blue 3G-A (CB) Dye from Aqueous Solution Using Chemo-Physically Activated Biochar from Oil Palm Empty Fruit Bunch Fiber. Arab. J. Chem. 2020, 13, 5417–5429. [Google Scholar] [CrossRef]
  23. Gupta, V.K.; Agarwal, S.; Ahmad, R.; Mirza, A.; Mittal, J. Sequestration of Toxic Congo Red Dye from Aqueous Solution Using Ecofriendly Guar Gum/ Activated Carbon Nanocomposite. Int. J. Biol. Macromol. 2020, 158, 1310–1318. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Z.; Hanafy, H.; Zhang, L.; Sellaoui, L.; Schadeck Netto, M.; Oliveira, M.L.S.; Seliem, M.K.; Luiz Dotto, G.; Bonilla-Petriciolet, A.; Li, Q. Adsorption of Congo Red and Methylene Blue Dyes on an Ashitaba Waste and a Walnut Shell-Based Activated Carbon from Aqueous Solutions: Experiments, Characterization and Physical Interpretations. Chem. Eng. J. 2020, 388, 124263. [Google Scholar] [CrossRef]
  25. Patawat, C.; Silakate, K.; Chuan-Udom, S.; Supanchaiyamat, N.; Hunt, A.J.; Ngernyen, Y. Preparation of Activated Carbon FromDipterocarpus Alatusfruit and Its Application for Methylene Blue Adsorption. RSC Adv. 2020, 10, 21082–21091. [Google Scholar] [CrossRef] [PubMed]
  26. Njewa, J.B.; Vunain, E.; Biswick, T. Synthesis and Characterization of Activated Carbons Prepared from Agro-Wastes by Chemical Activation. J. Chem. 2022, 2022, 9975444. [Google Scholar] [CrossRef]
  27. Ausavasukhi, A.; Kampoosaen, C.; Kengnok, O. Adsorption Characteristics of Congo Red on Carbonized Leonardite. J. Clean Prod. 2016, 134, 506–514. [Google Scholar] [CrossRef]
  28. Jawad, A.H.; Ismail, K.; Ishak, M.A.M.; Wilson, L.D. Conversion of Malaysian Low-Rank Coal to Mesoporous Activated Carbon: Structure Characterization and Adsorption Properties. Chin. J. Chem. Eng. 2019, 27, 1716–1727. [Google Scholar] [CrossRef]
  29. Osman, A.I.; Blewitt, J.; Abu-Dahrieh, J.K.; Farrell, C.; Al-Muhtaseb, A.H.; Harrison, J.; Rooney, D.W. Production and Characterisation of Activated Carbon and Carbon Nanotubes from Potato Peel Waste and Their Application in Heavy Metal Removal. Environ. Sci. Pollut. Res. 2019, 26, 37228–37241. [Google Scholar] [CrossRef] [PubMed]
  30. Homagai, P.L.; Poudel, R.; Poudel, S.; Bhattarai, A. Adsorption and Removal of Crystal Violet Dye from Aqueous Solution by Modified Rice Husk. Heliyon 2022, 8, e09261. [Google Scholar] [CrossRef] [PubMed]
  31. Reck, I.M.; Paixão, R.M.; Bergamasco, R.; Vieira, M.F.; Vieira, A.M.S. Removal of Tartrazine from Aqueous Solutions Using Adsorbents Based on Activated Carbon and Moringa Oleifera Seeds. J. Clean Prod. 2018, 171, 85–97. [Google Scholar] [CrossRef]
  32. Bansode, R.R.; Losso, J.N.; Marshall, W.E.; Rao, R.M.; Portier, R.J. Pecan Shell-Based Granular Activated Carbon for Treatment of Chemical Oxygen Demand (COD) in Municipal Wastewater. Bioresour. Technol. 2004, 94, 129–135. [Google Scholar] [CrossRef]
  33. Mandal, S.; Calderon, J.; Marpu, S.B.; Omary, M.A.; Shi, S.Q. Mesoporous Activated Carbon as a Green Adsorbent for the Removal of Heavy Metals and Congo Red: Characterization, Adsorption Kinetics, and Isotherm Studies. J. Contam. Hydrol. 2021, 243, 103869. [Google Scholar] [CrossRef] [PubMed]
  34. Lorenc-Grabowska, E.; Gryglewicz, G. Adsorption Characteristics of Congo Red on Coal-Based Mesoporous Activated Carbon. Dye. Pigment. 2007, 74, 34–40. [Google Scholar] [CrossRef]
  35. Lafi, R.; Montasser, I.; Hafiane, A. Adsorption of Congo Red Dye from Aqueous Solutions by Prepared Activated Carbon with Oxygen-Containing Functional Groups and Its Regeneration. Adsorpt. Sci. Technol. 2019, 37, 160–181. [Google Scholar] [CrossRef]
  36. Igwegbe, C.A.; Ighalo, J.O.; Onyechi, K.K.; Onukwuli, O.D. Adsorption of Congo Red and Malachite Green Using H3PO4 and NaCl-Modified Activated Carbon from Rubber (Hevea Brasiliensis) Seed Shells. Sustain. Water Resour. Manag. 2021, 7, 63. [Google Scholar] [CrossRef]
  37. Bessaha, F.; Mahrez, N.; Marouf-Khelifa, K.; Çoruh, A.; Khelifa, A. Removal of Congo Red by Thermally and Chemically Modified Halloysite: Equilibrium, FTIR Spectroscopy, and Mechanism Studies. Int. J. Environ. Sci. Technol. 2019, 16, 4253–4260. [Google Scholar] [CrossRef]
  38. Tran, H.N.; Tomul, F.; Thi Hoang Ha, N.; Nguyen, D.T.; Lima, E.C.; Le, G.T.; Chang, C.T.; Masindi, V.; Woo, S.H. Innovative Spherical Biochar for Pharmaceutical Removal from Water: Insight into Adsorption Mechanism. J. Hazard Mater. 2020, 394, 122255. [Google Scholar] [CrossRef]
  39. Wang, J.; Guo, X. Adsorption Kinetic Models: Physical Meanings, Applications, and Solving Methods. J. Hazard. Mater. 2020, 390, 122156. [Google Scholar] [CrossRef]
  40. Sabarinathan, C.; Karuppasamy, P.; Vijayakumar, C.T.; Arumuganathan, T. Development of Methylene Blue Removal Methodology by Adsorption Using Molecular Polyoxometalate: Kinetics, Thermodynamics and Mechanistic Study. Microchem. J. 2019, 146, 315–326. [Google Scholar] [CrossRef]
  41. Ho, Y.S.; Mckay, G. Sorption of Dye from Aqueous Solution by Peat. Chem. Eng. J. 1998, 70, 11–16. [Google Scholar] [CrossRef]
  42. Zhou, Q.; Xie, C.; Gong, W.; Xu, N.; Zhou, W. Comments on the Method of Using Maximum Absorption Wavelength to Calculate Congo Red Solution Concentration Published in J. Hazard. Mater. J. Hazard. Mater. 2011, 198, 381–382. [Google Scholar] [CrossRef]
  43. Nguyen, V.D.; Nguyen, H.T.H.; Vranova, V.; Nguyen, L.T.N.; Bui, Q.M.; Khieu, T.T. Artificial Neural Network Modeling for Congo Red Adsorption on Microwave-Synthesized Akaganeite Nanoparticles: Optimization, Kinetics, Mechanism, and Thermodynamics. Environ. Sci. Pollut. Res. 2021, 28, 9133–9145. [Google Scholar] [CrossRef] [PubMed]
  44. Batool, M.; Khurshid, S.; Qureshi, Z.; Daoush, W.M. Adsorption, Antimicrobial and Wound Healing Activities of Biosynthesised Zinc Oxide Nanoparticles. Chem. Pap. 2021, 75, 893–907. [Google Scholar] [CrossRef]
  45. Fawzy, M.A.; Gomaa, M. Low-Cost Biosorption of Methylene Blue and Congo Red from Single and Binary Systems Using Sargassum Latifolium Biorefinery Waste/Wastepaper Xerogel: An Optimization and Modeling Study. J. Appl. Phycol. 2021, 33, 675–691. [Google Scholar] [CrossRef]
  46. Yaneva, Z.; Georgieva, N.V.; Yaneva, Z.L. Insights into Congo Red Adsorption on Agro-Industrial Materials-Spectral, Equilibrium, Kinetic, Thermodynamic, Dynamic and Desorption Studies. A Review. Int. Rev. Chem. Eng. IRECHE 2012, 4, 127–146. [Google Scholar]
  47. Al-Ghouti, M.A.; Da’ana, D.A. Guidelines for the Use and Interpretation of Adsorption Isotherm Models: A Review. J. Hazard. Mater. 2020, 393, 122383. [Google Scholar] [CrossRef] [PubMed]
  48. Roy, T.K.; Mondal, N.K. Potentiality of Eichhornia Shoots Ash towards Removal of Congo Red from Aqueous Solution: Isotherms, Kinetics, Thermodynamics and Optimization Studies. Groundw. Sustain. Dev. 2019, 9, 100269. [Google Scholar] [CrossRef]
  49. Alghamdi, A.A.; Al-Odayni, A.B.; Saeed, W.S.; Almutairi, M.S.; Alharthi, F.A.; Aouak, T.; Al-Kahtani, A. Adsorption of Azo Dye Methyl Orange from Aqueous Solutions Using Alkali-Activated Polypyrrole-Based Graphene Oxide. Molecules 2019, 24, 3685. [Google Scholar] [CrossRef] [PubMed]
  50. Lawal, I.A.; Chetty, D.; Akpotu, S.O.; Moodley, B. Sorption of Congo Red and Reactive Blue on Biomass and Activated Carbon Derived from Biomass Modified by Ionic Liquid. Environ. Nanotechnol. Monit. Manag. 2017, 8, 83–91. [Google Scholar] [CrossRef]
  51. Amran, F.; Zaini, M.A.A. Sodium Hydroxide-Activated Casuarina Empty Fruit: Isotherm, Kinetics and Thermodynamics of Methylene Blue and Congo Red Adsorption. Environ. Technol. Innov. 2021, 23, 101727. [Google Scholar] [CrossRef]
  52. Ojedokun, A.T.; Bello, O.S. Kinetic Modeling of Liquid-Phase Adsorption of Congo Red Dye Using Guava Leaf-Based Activated Carbon. Appl. Water Sci. 2017, 7, 1965–1977. [Google Scholar] [CrossRef]
  53. Omidi Khaniabadi, Y.; Mohammadi, M.J.; Shegerd, M.; Sadeghi, S.; Basiri, H. Removal of Congo Red Dye from Aqueous Solutions by a Low-Cost Adsorbent: Activated Carbon Prepared from Aloe Vera Leaves Shell. Environ. Health Eng. Manag. 2016, 4, 29–35. [Google Scholar] [CrossRef]
  54. Sharma, A.; Siddiqui, Z.M.; Dhar, S.; Mehta, P.; Pathania, D. Adsorptive Removal of Congo Red Dye (CR) from Aqueous Solution by Cornulaca Monacantha Stem and Biomass-Based Activated Carbon: Isotherm, Kinetics and Thermodynamics. Sep. Sci. Technol. 2019, 54, 916–929. [Google Scholar] [CrossRef]
  55. Latinwo, G.K.; Alade, A.O.; Agarry, S.E.; Dada, E.O. Process Optimization and Modeling the Adsorption of Polycyclic Aromatic-Congo Red Dye onto Delonix Regia Pod-Derived Activated Carbon. Polycycl. Aromat. Compd. 2021, 41, 400–418. [Google Scholar] [CrossRef]
  56. Ojedokun, A.T.; Bello, O.S. Liquid Phase Adsorption of Congo Red Dye on Functionalized Corn Cobs. J. Dispers. Sci. Technol. 2017, 38, 1285–1294. [Google Scholar] [CrossRef]
  57. El Messaoudi, N.; El Khomri, M.; Chlif, N.; Chegini, Z.G.; Dbik, A.; Bentahar, S.; Lacherai, A. Desorption of Congo Red from Dye-Loaded Phoenix Dactylifera Date Stones and Ziziphus Lotus Jujube Shells. Groundw. Sustain. Dev. 2021, 12, 100552. [Google Scholar] [CrossRef]
  58. Shahzad, M.W.; Nguyen, V.H.; Xu, B.B.; Tariq, R.; Imran, M.; Ashraf, W.M.; Ng, K.C.; Jamil, M.A.; Ijaz, A.; Sheikh, N.A. Machine Learning Assisted Prediction of Solar to Liquid Fuel Production: A Case Study. Process Saf. Environ. Prot. 2024, 184, 1119–1130. [Google Scholar] [CrossRef]
  59. Yin, G.; Jameel Ibrahim Alazzawi, F.; Mironov, S.; Reegu, F.; El-Shafay, A.S.; Lutfor Rahman, M.; Su, C.H.; Lu, Y.Z.; Chinh Nguyen, H. Machine Learning Method for Simulation of Adsorption Separation: Comparisons of Model’s Performance in Predicting Equilibrium Concentrations. Arab. J. Chem. 2022, 15, 103612. [Google Scholar] [CrossRef]
  60. Zhao, C.; Zhang, W.; Zhang, Y.; Yang, Y.; Guo, D.; Liu, W.; Liu, L. Influence of Multivalent Background Ions Competition Adsorption on the Adsorption Behavior of Azo Dye Molecules and Removal Mechanism: Based on Machine Learning, DFT and Experiments. Sep. Purif. Technol. 2024, 341, 126810. [Google Scholar] [CrossRef]
  61. Zhao, S.; Chen, K.; Xiong, B.; Guo, C.; Dang, Z. Prediction of Adsorption of Metal Cations by Clay Minerals Using Machine Learning. Sci. Total Environ. 2024, 924, 171733. [Google Scholar] [CrossRef]
  62. Cai, Q.Y.; Qiao, L.Z.; Yao, S.J.; Lin, D.Q. Machine Learning Assisted QSAR Analysis to Predict Protein Adsorption Capacities on Mixed-Mode Resins. Sep. Purif. Technol. 2024, 340, 126762. [Google Scholar] [CrossRef]
  63. Lee, H.; Choi, Y. Predicting Apparent Adsorption Capacity of Sediment-Amended Activated Carbon for Hydrophobic Organic Contaminants Using Machine Learning. Chemosphere 2024, 350, 141003. [Google Scholar] [CrossRef]
  64. Guo, L.; Xu, X.; Niu, C.; Wang, Q.; Park, J.; Zhou, L.; Lei, H.; Wang, X.; Yuan, X. Machine Learning-Based Prediction and Experimental Validation of Heavy Metal Adsorption Capacity of Bentonite. Sci. Total Environ. 2024, 926, 171986. [Google Scholar] [CrossRef] [PubMed]
  65. Chen, K.; Guo, C.; Wang, C.; Zhao, S.; Lu, G.; Dang, Z. Using Machine Learning to Explore Oxyanion Adsorption Ability of Goethite with Different Specific Surface Area. Environ. Pollut. 2024, 343, 123162. [Google Scholar] [CrossRef]
  66. Selvaraj, R.; Jogi, S.; Murugesan, G.; Srinivasan, N.R.; Goveas, L.C.; Varadavenkatesan, T.; Samanth, A.; Vinayagam, R.; Ali Alshehri, M.; Pugazhendhi, A. Machine Learning and Statistical Physics Modeling of Tetracycline Adsorption Using Activated Carbon Derived from Cynometra Ramiflora Fruit Biomass. Environ. Res. 2024, 252, 118816. [Google Scholar] [CrossRef] [PubMed]
  67. Shboul, B.; Zayed, M.E.; Tariq, R.; Ashraf, W.M.; Odat, A.S.; Rehman, S.; Abdelrazik, A.S.; Krzywanski, J. New Hybrid Photovoltaic-Fuel Cell System for Green Hydrogen and Power Production: Performance Optimization Assisted with Gaussian Process Regression Method. Int. J. Hydrogen Energy 2024, 59, 1214–1229. [Google Scholar] [CrossRef]
  68. Ashraf, W.M.; Uddin, G.M.; Tariq, R.; Ahmed, A.; Farhan, M.; Nazeer, M.A.; Hassan, R.U.; Naeem, A.; Jamil, H.; Krzywanski, J.; et al. Artificial Intelligence Modeling-Based Optimization of an Industrial-Scale Steam Turbine for Moving toward Net-Zero in the Energy Sector. ACS Omega 2023, 8, 21709–21725. [Google Scholar] [CrossRef]
  69. Muhammad Ashraf, W.; Moeen Uddin, G.; Afroze Ahmad, H.; Ahmad Jamil, M.; Tariq, R.; Wakil Shahzad, M.; Dua, V. Artificial Intelligence Enabled Efficient Power Generation and Emissions Reduction Underpinning Net-Zero Goal from the Coal-Based Power Plants. Energy Convers. Manag. 2022, 268, 116025. [Google Scholar] [CrossRef]
  70. Ji, Y.; Zhang, X.; Chen, Z.; Xiao, Y.; Li, S.; Gu, J.; Hu, H.; Cheng, G. Silk Sericin Enrichment through Electrodeposition and Carbonous Materials for the Removal of Methylene Blue from Aqueous Solution. Int. J. Mol. Sci. 2022, 23, 1668. [Google Scholar] [CrossRef]
  71. Dada, A.O.; Adekola, F.A.; Odebunmi, E.O.; Ogunlaja, A.S.; Bello, O.S. Two–Three Parameters Isotherm Modeling, Kinetics with Statistical Validity, Desorption and Thermodynamic Studies of Adsorption of Cu(II) Ions onto Zerovalent Iron Nanoparticles. Sci. Rep. 2021, 11, 16454. [Google Scholar] [CrossRef]
Figure 1. pHpzc of ABHC.
Figure 1. pHpzc of ABHC.
Ijms 25 04771 g001
Figure 2. XRD pattern of ABHC and ABHC-CR.
Figure 2. XRD pattern of ABHC and ABHC-CR.
Ijms 25 04771 g002
Figure 3. SEM images of (a) ABHC and (b) ABHC-CR. EDS results of (c) ABHC and (d) ABHC-CR.
Figure 3. SEM images of (a) ABHC and (b) ABHC-CR. EDS results of (c) ABHC and (d) ABHC-CR.
Ijms 25 04771 g003
Figure 4. Adsorption–desorption isotherm of N2 on the ABHC.
Figure 4. Adsorption–desorption isotherm of N2 on the ABHC.
Ijms 25 04771 g004
Figure 5. FT-IR of ABHC before and after CR adsorption.
Figure 5. FT-IR of ABHC before and after CR adsorption.
Ijms 25 04771 g005
Figure 6. Kinetic modeling of Congo red adsorption (a) Ci = 25 mg/L and (b) Ci = 50 mg/L pseudo-first-order, pseudo-second-order, and Elovich nonlinear models by ABHC (pH = 7; D ABHC = 1 g/L; T = 300.15 K; contact time = 15–2880 min).
Figure 6. Kinetic modeling of Congo red adsorption (a) Ci = 25 mg/L and (b) Ci = 50 mg/L pseudo-first-order, pseudo-second-order, and Elovich nonlinear models by ABHC (pH = 7; D ABHC = 1 g/L; T = 300.15 K; contact time = 15–2880 min).
Ijms 25 04771 g006
Figure 7. Adsorption isotherms modeling of Congo red at different pH values: (a) pH = 4.0, (b) pH = 5.4, (c) pH = 7.0, (d) pH = 8.4, and (e) pH = 10.2 (T = 300.15 K; CiCR = 10 mg/L–100 mg/L; dose = 1 g/L; contact time = 1440 min).
Figure 7. Adsorption isotherms modeling of Congo red at different pH values: (a) pH = 4.0, (b) pH = 5.4, (c) pH = 7.0, (d) pH = 8.4, and (e) pH = 10.2 (T = 300.15 K; CiCR = 10 mg/L–100 mg/L; dose = 1 g/L; contact time = 1440 min).
Ijms 25 04771 g007aIjms 25 04771 g007b
Figure 8. Adsorption isotherms modeling of Congo red at different ABHC doses: (a) dose = 1 g/L, (b) dose = 2 g/L, (c) dose = 5 g/L, and (d) dose = 10 g/L (pH = 7; CiCR = 20–1000 mg/L; T = 300.15 K; contact time = 1440 min).
Figure 8. Adsorption isotherms modeling of Congo red at different ABHC doses: (a) dose = 1 g/L, (b) dose = 2 g/L, (c) dose = 5 g/L, and (d) dose = 10 g/L (pH = 7; CiCR = 20–1000 mg/L; T = 300.15 K; contact time = 1440 min).
Ijms 25 04771 g008
Figure 9. Adsorption isotherms modeling of Congo red at different temperatures: (a) T = 300.15, (b) T = 313.15, (c) T = 330.15 (pH = 7; CiCR = 20–1000 mg/L; dose = 1 g/L; contact time = 360 min).
Figure 9. Adsorption isotherms modeling of Congo red at different temperatures: (a) T = 300.15, (b) T = 313.15, (c) T = 330.15 (pH = 7; CiCR = 20–1000 mg/L; dose = 1 g/L; contact time = 360 min).
Ijms 25 04771 g009
Figure 10. Effect of different effluent concentrations on biosorption/desorption cycles of CR onto ABHC (a) 0.01M NaOH and (b) 0.2M NaOH.
Figure 10. Effect of different effluent concentrations on biosorption/desorption cycles of CR onto ABHC (a) 0.01M NaOH and (b) 0.2M NaOH.
Ijms 25 04771 g010
Figure 11. Scatter matrix showing the relationship between the input variables and the removal percentage.
Figure 11. Scatter matrix showing the relationship between the input variables and the removal percentage.
Ijms 25 04771 g011
Figure 12. Heatmap based upon coefficient of correlation.
Figure 12. Heatmap based upon coefficient of correlation.
Ijms 25 04771 g012
Figure 13. Regression fit model between the actual removal percentage and the predicted removal percentage obtained through the application of Gradient Boosting method.
Figure 13. Regression fit model between the actual removal percentage and the predicted removal percentage obtained through the application of Gradient Boosting method.
Ijms 25 04771 g013
Figure 14. Parallel coordinate plot showing the multi-dimensional visualization of the data variables.
Figure 14. Parallel coordinate plot showing the multi-dimensional visualization of the data variables.
Ijms 25 04771 g014
Figure 15. Convergence plot for the optimization process using Bayesian optimization algorithm.
Figure 15. Convergence plot for the optimization process using Bayesian optimization algorithm.
Ijms 25 04771 g015
Figure 16. Chemical structure of Congo red (CR) (ChemDraw professional software).
Figure 16. Chemical structure of Congo red (CR) (ChemDraw professional software).
Ijms 25 04771 g016
Table 1. FT-IR bands before and after the adsorption process.
Table 1. FT-IR bands before and after the adsorption process.
Functional GroupWavenumber (cm−1)
ABHCABHC-CR
O-H and N-H33583368
C=O17051705
C=C15931600
C-C (in aromatic ring)-1421
C-O-1360
C-O12241196
S=O-1045
N-H779769
Table 2. Kinetic parameters of CR adsorption onto (ABHC) at different Ci values.
Table 2. Kinetic parameters of CR adsorption onto (ABHC) at different Ci values.
Ci
(mg/L)
Pseudo-First-OrderPseudo-Second-OrderElovich
qe,exp
(mg/g)
qe
(mg/g)
k1/10−2
(1/min)
R2qe
(mg/g)
k2/10−3
[g/(mg·min)]
R2α
[mg/(g·min)]
β
(g/mg)
R2
2515.2113.871.1870.95715.091.0560.9880.81000.41830.928
5023.7421.950.56930.87723.810.37240.9390.78800.25900.968
Table 3. Error functions data of the kinetic models.
Table 3. Error functions data of the kinetic models.
Ci
(mg/L)
ModelError Functions
ARESSE∆q (%)χ2EABSRMSE
25PFO7.5935.9188.9960.5120.8191.088
PSO7.6691.59915.0850.3730.3890.565
Elovich20.5189.89141.8922.6171.0331.406
50PFO19.16037.76332.0244.4761.9742.748
PSO12.08118.76223.5872.3311.1981.937
Elovich7.81810.00110.3370.6931.0401.414
Table 4. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto activated biochar from Haematoxylum campechianum (ABHC) at different pH values.
Table 4. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto activated biochar from Haematoxylum campechianum (ABHC) at different pH values.
Isotherm ModelsParameterspH
45.478.410.2
LangmuirQmax (mg/g)53.93114.868.3330.8310.51
KL/10−2 (L/mg)7.2521.1083.6790.8329.695
R20.9820.9940.9720.9880.962
n2.4791.3351.8301.3623.4315
FreundlichKF (mg/g)(L/mg)1/n8.9452.1895.3930.4982.684
R20.9400.9900.9920.9830.941
Redlich–PetersonαRP (L/mg)0.0850.0011.7260.0150.214
KRP (L/g)4.0881.14112.7610.2701.395
β0.9731.4530.5160.8870.895
R20.9820.9950.9920.9880.969
Table 5. Error functions data of the pH effect.
Table 5. Error functions data of the pH effect.
pHModelError Functions
ARESSE∆q (%)χ2EABSRMSE
4.0Langmuir4.44011.515.9220.3941.1611.696
Freundlich11.4937.9919.432.4392.2143.082
Redlich–Peterson4.73111.416.4430.4211.1861.689
5.4Langmuir6.0506.1199.4950.4950.8161.106
Freundlich5.24610.647.9100.5571.0921.459
Redlich–Peterson5.9975.69210.970.5650.7101.067
7.0Langmuir10.3232.0915.521.9331.8112.533
Freundlich6.0608.5349.4910.6020.9471.306
Redlich–Peterson5.3178.6158.0120.5200.8991.313
8.4Langmuir5.6121.2467.8020.1530.3540.499
Freundlich9.3241.67216.820.3820.4070.578
Redlich–Peterson5.8181.2378.7310.1660.3500.497
10.2Langmuir4.5871.0495.8990.1430.3250.458
Freundlich6.9171.61711.290.3340.3860.569
Redlich–Peterson4.6110.8366.5310.1390.2870.409
Table 6. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto different doses of activated biochar from Haematoxylum campechianum (ABHC).
Table 6. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto different doses of activated biochar from Haematoxylum campechianum (ABHC).
Isotherm ModelsParametersDose (g/L)
12510
LangmuirQmax (mg/g)92.8659.5347.6212.22
KL/10−2 (L/mg)1.8870.3690.1250.685
R20.9400.9880.9700.966
n3.0792.0511.5792.551
FreundlichKF (mg/g)(L/mg)1/n10.851.7470.3460.787
R20.9760.9930.9850.982
Redlich–PetersonαRP (L/mg)0.4010.091252.83.338
KRP (L/g)6.7880.50687.642.811
β0.7420.6690.3670.618
R20.9840.9970.9850.982
Table 7. Error functions data of the dose effect.
Table 7. Error functions data of the dose effect.
Dose (g/L)ModelError Functions
ARESSE∆q (%)χ2EABSRMSE
1Langmuir14.52552.422.1710.835.0047.835
Freundlich16.67224.932.1410.603.9504.999
Redlich–Peterson7.262151.710.662.8132.8864.105
2Langmuir12.2219.5520.211.8811.3551.977
Freundlich9.94512.4619.841.2151.0831.579
Redlich–Peterson4.8304.9126.9410.3320.6980.991
5Langmuir26.8914.4143.573.4371.2151.698
Freundlich13.677.35226.161.2640.7771.213
Redlich–Peterson13.687.35626.181.2660.7781.213
10Langmuir20.932.99636.481.3460.5730.774
Freundlich9.0201.54913.550.3150.3630.557
Redlich–Peterson9.9771.54314.890.3420.3750.556
Table 8. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto activated biochar from Haematoxylum campechianum (ABHC) at different temperatures.
Table 8. Parameters of the Langmuir, Freundlich, and Redlich–Peterson nonlinear models for the adsorption of CR onto activated biochar from Haematoxylum campechianum (ABHC) at different temperatures.
Isotherm ModelsParametersTemperature (K)
300.15313.15330.15
LangmuirQmax (mg/g)92.8627.4429.20
KL/10−2 (L/mg)1.88726.781.917
R20.9400.7090.980
n3.0797.4113.878
FreundlichKF (mg/g)(L/mg)1/n10.8512.385.106
R20.9760.9990.951
Redlich–PetersonαRP (L/mg)0.40119.660.057
KRP (L/g)6.788247.20.852
β0.7420.8670.899
R20.9840.9990.988
Table 9. Error functions data of temperature effect.
Table 9. Error functions data of temperature effect.
Temperature (K)ModelError Functions
ARESSE∆q (%)χ2EABSRMSE
300.15Langmuir14.52552.422.1710.835.0047.835
Freundlich16.67224.932.1410.603.9504.999
Redlich–Peterson7.262151.710.662.8132.8864.105
313.15Langmuir13.6471.5716.613.3162.9783.783
Freundlich0.8130.3211.2110.0170.1620.254
Redlich–Peterson0.6730.2891.0540.0140.1450.240
330.15Langmuir6.9218.18314.160.9580.7461.279
Freundlich9.95020.1313.091.2451.5802.006
Redlich–Peterson5.7204.8078.7740.4330.7520.980
Table 10. Comparison of the maximum adsorption capacity of Congo red onto different activated carbons.
Table 10. Comparison of the maximum adsorption capacity of Congo red onto different activated carbons.
SubstrateQmax (mg/g)T (°C)pHCi (mg/L)References
Coffee waste90.90253.020–120[35]
Kenaf fiber (Hibiscus cannabinus)14.202775–25[33]
Guava leaves47.6230310–50[52]
Rubber (Hevea brasiliensis)55.87302100–500[36]
Aloe vera leaves91.00252100[53]
Cornulaca monacantha78.19552.020–160[54]
Peanut shell153.4--20–200[50]
Casuarinas waste232.025-5–1000[51]
Delonix regia17.1230-200–1200[55]
Corn cobs41.6750310–50[56]
Haematoxylum campechianum114.827510–100This study
Table 11. Performance metrices of Gradient Boosting method during the training and testing phase.
Table 11. Performance metrices of Gradient Boosting method during the training and testing phase.
Type of DataMSESSEMAPERMSEMPECOD
Training data4.052313340.39435.9470512.013036−0.924480.992258
Testing data33.1685696.538623.088195.75921−14.27490.91353
Table 12. Optimal removal percentage and the input variables.
Table 12. Optimal removal percentage and the input variables.
Initial Concentration
(mg/L)
Time
(min)
Temperature
(K)
pHRemoval Percentage
10.02880313.154.090.4733
Table 13. Physical properties of Congo red (CR).
Table 13. Physical properties of Congo red (CR).
ParameterValue
Molecular weight696.66 g/mol
Density0.995 g/mL at 25 °C
SolubilityH2O: 25 g/L
Water solubility
pKa
Soluble
4.5 [44]
pH6.7 (10 g/L, H2O and at 20 °C)
Color and pH range3 (blue)–5.2 (red)
λmax567 nm at pH 2.18–3.16 and 497 nm at pH ≥ 3.86
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gamboa, D.M.P.; Abatal, M.; Lima, E.; Franseschi, F.A.; Ucán, C.A.; Tariq, R.; Elías, M.A.R.; Vargas, J. Sorption Behavior of Azo Dye Congo Red onto Activated Biochar from Haematoxylum campechianum Waste: Gradient Boosting Machine Learning-Assisted Bayesian Optimization for Improved Adsorption Process. Int. J. Mol. Sci. 2024, 25, 4771. https://doi.org/10.3390/ijms25094771

AMA Style

Gamboa DMP, Abatal M, Lima E, Franseschi FA, Ucán CA, Tariq R, Elías MAR, Vargas J. Sorption Behavior of Azo Dye Congo Red onto Activated Biochar from Haematoxylum campechianum Waste: Gradient Boosting Machine Learning-Assisted Bayesian Optimization for Improved Adsorption Process. International Journal of Molecular Sciences. 2024; 25(9):4771. https://doi.org/10.3390/ijms25094771

Chicago/Turabian Style

Gamboa, Diego Melchor Polanco, Mohamed Abatal, Eder Lima, Francisco Anguebes Franseschi, Claudia Aguilar Ucán, Rasikh Tariq, Miguel Angel Ramírez Elías, and Joel Vargas. 2024. "Sorption Behavior of Azo Dye Congo Red onto Activated Biochar from Haematoxylum campechianum Waste: Gradient Boosting Machine Learning-Assisted Bayesian Optimization for Improved Adsorption Process" International Journal of Molecular Sciences 25, no. 9: 4771. https://doi.org/10.3390/ijms25094771

APA Style

Gamboa, D. M. P., Abatal, M., Lima, E., Franseschi, F. A., Ucán, C. A., Tariq, R., Elías, M. A. R., & Vargas, J. (2024). Sorption Behavior of Azo Dye Congo Red onto Activated Biochar from Haematoxylum campechianum Waste: Gradient Boosting Machine Learning-Assisted Bayesian Optimization for Improved Adsorption Process. International Journal of Molecular Sciences, 25(9), 4771. https://doi.org/10.3390/ijms25094771

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop