Removal of Methyl Violet from Aqueous Solution by Adsorption onto Halloysite Nanoclay: Experiment and Theory

Methyl Violet (MV) was removed from aqueous solutions by adsorption onto halloysite nanoclay (HNC) employing equilibrium, kinetics, thermodynamic data, molecular modellingR (MD), and Monte Carlo (MC) simulations. The chosen experimental variables were pH, temperature, starting MV concentration, contact time, and adsorbent dosage. The adsorption rate was determined to increase with increasing contact time, initial dye concentration, pH, and temperature. The Langmuir, Freundlich, Temkin, and Dubinin–Radushkevich (D-R) isotherms were utilized to determine the adsorption capacity of HNC. The Langmuir equation matched equilibrium data better than the other models, whereas the pseudo-second-order model better described kinetic data, and thermodynamic analyses revealed that the adsorption process was spontaneous, endothermic, and physisorption-based. This study focused on two distinct molecular mechanics-based theoretical approaches (MC and MD). These techniques enabled a molecular comprehension of the interaction between the MV molecule and the halloysite surface. Theoretical results were consistent with experimental findings. The outcomes revealed that HNC is an excellent dye adsorbent for industrial effluents.


Introduction
Water is a crucial component for life, human evolution, and biodiversity conservation. Unfortunately, human activities have caused a serious degeneration in the quality of water resources. Even though a number of regulations are designed to ensure the safe disposal of industrial effluents, over 80 percent of the world's wastewater is still disposed of without sufficient treatment [1][2][3]. The release of dye-containing effluents into the natural water bodies can cause harmful effects on the living systems because of their toxic, allergenic, mutagenic, and carcinogenic nature. Dyes obstruct light penetration, impede photosynthetic activity, and as a result, hinder the growth of biota and also cause micro-toxicity to fish and other organisms because of their predisposition to chelate metal ions [4,5].
In several industries such as food, plastic, textile, leather, paper, cosmetics, etc., organic dye effluent is one of the main water pollution problems [6]. Methyl violet (MV) is vital because of its extensive use in textiles, paints, and print inks, and the dying of cotton, silk, paper, and leather, among other applications [7]. MV is also used in medicine as an active ingredient in Gram's biological stain for bacteria categorization [8][9][10][11]. Sometimes it can be used as a second-rate disinfectant but has been found to be toxic to most animals. Inhaling MV may irritate the respiratory tract, while ingesting it typically irritates the digestive tract [12,13].
As adsorbate, Methyl Violet dye purchased from Fluka was utilized and was used without further purification. The stock solution (1 g L −1 ) was prepared by dissolving 1 g of the dye into distilled water. The experimental dye solution was made by diluting the stock solution with the necessary amount of distilled water.

Batch Experimental Procedure
Batch experiments were performed by stirring 25 mL of known concentration of MV solution with 0.2 g of HNC using magnetic stirrer.
The pH was adjusted using 0.1 M HCl and 0.1 M NaOH. After adsorption, the adsorbent and the supernatants were separated by centrifugation at 5000 rpm for 10 min and samples were investigated for residual dye concentration using a UV-Visible Spectrophotometer (type T70+), at wavelength 579 nm, using a 1 cm quartz cell.
All experiments were repeated three times and mean values were presented with a maximum deviation of 5%.
The amount of MV adsorbed per unit adsorbent (mg/g) was calculated using Equation (1): where C i is the initial dye concentration (mg/L), C e is the equilibrium dye concentration (mg/L), V is the volume of MV solution used (L), and m(g) is the mass of HNC used. The MV percent removal was calculated using Equation (2).
To construct and evaluate an adsorption process, it is necessary to fit equilibrium adsorption data with various adsorption isotherm and kinetic models. Thus, several theoretical models (Table 1) are used for experimental data in order to identify a model that predicts kinetic and isotherm data adequately. The validity of the models was determined by calculating the regression coefficient (R 2 ) and the root mean square errors (RMSE): where q exp is the experimental value while q cal is the calculated value and N is the number of observations in the experiment. The smaller the RMSE values, the better the curve fitting [50,51]. Exel 2016 was used to create all of the graphs. Table 1. Isotherm, kinetic, and other equations used in this study.

Monte Carlo and Molecular Dynamic Simulation Details
Six cycles of simulated annealing, each with 15,000 steps, were used to do the MC computations. The annealing temperature was automatically selected between 10 5 and 10 2 K for each cycle. As the temperature steadily decreased, potential adsorption structures were uncovered [52,53]. The MD is performed under an NVT ensemble at 25 °C with a 1 fs time step and 1 ns total simulation period. The temperature control is realized using the Berendsen thermostat MD [54]. MC and MD simulations use COMPASS II's force field. Radial Distribution Function (RDF) is computed using 1 ns of MD trajectory [55,56].

Effect of pH
Solution pH affects adsorption. It affects the adsorbent's surface charge, the solution's ionization, and the dissociation of functional groups on the active sites, as well as solution dye chemistry [57,58]. Figure 2 shows the effect of pH on MV removal effectiveness from 4.0 to 10.5. As illustrated in Figure 2, the pH rise from 4 to 10.5, has an insignificant effect on the adsorption capacity. This is in line with other studies such as the

Monte Carlo and Molecular Dynamic Simulation Details
Six cycles of simulated annealing, each with 15,000 steps, were used to do the MC computations. The annealing temperature was automatically selected between 10 5 and 10 2 K for each cycle. As the temperature steadily decreased, potential adsorption structures were uncovered [52,53]. The MD is performed under an NVT ensemble at 25 • C with a 1 fs time step and 1 ns total simulation period. The temperature control is realized using the Berendsen thermostat MD [54]. MC and MD simulations use COMPASS II's force field. Radial Distribution Function (RDF) is computed using 1 ns of MD trajectory [55,56].

Effect of pH
Solution pH affects adsorption. It affects the adsorbent's surface charge, the solution's ionization, and the dissociation of functional groups on the active sites, as well as solution dye chemistry [57,58]. Figure 2 shows the effect of pH on MV removal effectiveness from 4.0 to 10.5. As illustrated in Figure 2, the pH rise from 4 to 10.5, has an insignificant effect on the adsorption capacity. This is in line with other studies such as the adsorption of methyl violet on halloysite nanotube [49] and the adsorption of MG on treated ginger waste [59].

Effect of Contact Time and Initial Concentration
In Figure 3  adsorption of methyl violet on halloysite nanotube [49] and the adsorption of MG on treated ginger waste [59].

Effect of Contact Time and Initial Concentration
In Figure   MV adsorption was faster in the first 10 min, but slowed when the halloysite became saturated, reaching equilibrium after 240 min. High beginning adsorption rates may be owing to a large number of binding sites, whereas lower rates at the end are due to saturation and equilibrium. Similar conclusions have been proposed [60] for the adsorption of MG by neem sawdust. treated ginger waste [59].

Effect of Contact Time and Initial Concentration
In Figure   MV adsorption was faster in the first 10 min, but slowed when the halloysite became saturated, reaching equilibrium after 240 min. High beginning adsorption rates may be owing to a large number of binding sites, whereas lower rates at the end are due to saturation and equilibrium. Similar conclusions have been proposed [60] for the adsorption of MG by neem sawdust. MV adsorption was faster in the first 10 min, but slowed when the halloysite became saturated, reaching equilibrium after 240 min. High beginning adsorption rates may be owing to a large number of binding sites, whereas lower rates at the end are due to saturation and equilibrium. Similar conclusions have been proposed [60] for the adsorption of MG by neem sawdust.
Also, adsorption increases with increasing initial dye concentration. According to the results obtained, adsorption capacity increased from 12.46 mg/g to 24.66 mg/g as the initial concentration of MV dye increased from 100 mg/L to 200 mg/L. As initial concentration increases, mass transfer driving force overcomes resistances to dye molecule mass transfer from solution to the solid phase. Also, increasing concentration promotes dye adsorption by increasing dye-halloysite interaction [61,62].

Effect of Adsorbent Dosage
It was noticed (Figure 4), that the percent of MV adsorption grew significantly from 91-99.9%, with the increase of adsorbent dosage from 0.025 to 0.2 g. Thereafter, by increasing the adsorbent dosage up to 0.4 g, an insignificant increase in the removal efficiency was observed [63].
Also, adsorption increases with increasing initial dye concentration. According to the results obtained, adsorption capacity increased from 12.46 mg/g to 24.66 mg/g as the initial concentration of MV dye increased from 100 mg/L to 200 mg/L. As initial concentration increases, mass transfer driving force overcomes resistances to dye molecule mass transfer from solution to the solid phase. Also, increasing concentration promotes dye adsorption by increasing dye-halloysite interaction [61,62].

Effect of Adsorbent Dosage
It was noticed (Figure 4), that the percent of MV adsorption grew significantly from 91-99.9%, with the increase of adsorbent dosage from 0.025 to 0.2 g. Thereafter, by increasing the adsorbent dosage up to 0.4 g, an insignificant increase in the removal efficiency was observed [63]. It was also found that the enhancement of adsorbent dosage, resulted in a decrease in the amount of adsorbed dye per unit mass of adsorbent, from 91.3 to 6.3 mg/g.
Due to a divide in the concentration gradient between dye concentration in the solution and dye concentration on the HNC surface, qe (mg/g) decreases with increasing adsorbent mass [64].
A similar way of behaving was reported for mercury (II) removal on EB [65] and MB adsorption on guava leaf [66], on gulmohar plant leaf [67], and on cashew nutshell activated carbon [68].

Effect of Temperature
The effect of temperature on MV dye adsorption onto HNC is introduced in Figure  5. It was also found that the enhancement of adsorbent dosage, resulted in a decrease in the amount of adsorbed dye per unit mass of adsorbent, from 91.3 to 6.3 mg/g.
Due to a divide in the concentration gradient between dye concentration in the solution and dye concentration on the HNC surface, q e (mg/g) decreases with increasing adsorbent mass [64].
A similar way of behaving was reported for mercury (II) removal on EB [65] and MB adsorption on guava leaf [66], on gulmohar plant leaf [67], and on cashew nutshell activated carbon [68].

Effect of Temperature
The effect of temperature on MV dye adsorption onto HNC is introduced in Figure 5. The equilibrium adsorption capacity of MV onto HNC was found to increase slightly with increasing temperature, from 18.63 mg/g at 284.15 K to 18.71 mg/g at 313.15 K. With increasing temperature, the viscosity of the solution decreases, and consequently the rate of diffusion of dye molecules increases. This leads to an increase in mobility of molecules from the bulk solution to the surface of the adsorbent, and as a result, an increase in the amount of MV adsorbed [69,70].
A similar temperature effect was reported for adsorption of MV adsorption by clinoptilolite [71] and adsorption of the dyestuff astrazon red violet 3 rn (basic violet 16) on montmorillonite clay [72]. The equilibrium adsorption capacity of MV onto HNC was found to increase slightly with increasing temperature, from 18.63 mg/g at 284.15 K to 18.71 mg/g at 313.15 K. With increasing temperature, the viscosity of the solution decreases, and consequently the rate of diffusion of dye molecules increases. This leads to an increase in mobility of molecules from the bulk solution to the surface of the adsorbent, and as a result, an increase in the amount of MV adsorbed [69,70].
A similar temperature effect was reported for adsorption of MV adsorption by clinoptilolite [71] and adsorption of the dyestuff astrazon red violet 3 rn (basic violet 16) on montmorillonite clay [72].

Adsorption Isotherms
Four isotherm equations specifically Langmuir, Freundlich, Temkin and Dubinin-Radushkevich (D-R) ( Table 1) were employed to interpret the experimental data for MV adsorption on HNC at different temperatures.
The results are shown in Table 2.

Adsorption Isotherms
Four isotherm equations specifically Langmuir, Freundlich, Temkin and Dubinin-Radushkevich (D-R) ( Table 1) were employed to interpret the experimental data for MV adsorption on HNC at different temperatures.
The results are shown in Table 2. Temkin Based on the R 2 , RMSE values, q e,calc values, and modelled isotherms, the Langmuir isotherm model fitted best amongst all the isotherm models investigated. Langmuir isotherm model assumes similar, energetically equivalent adsorption sites and monolayer adsorption [73].
The Langmuir isotherm constants, K L and q m , were calculated from the slope and intercept respectively, of linear plot of C e /q e against C e ( Figure 5 and Table 2).
The separation factor (R L ), an essential parameter of the Langmuir isotherm, is generally used to show if the adsorption process is unfavorable (R L > 1), linear (R L = 1), favorable (0 < R L > 1), or irreversible (R L = 0). It can be calculated by Equation (4).
The values of R L in Table 2 indicate that the adsorption of MV on HNC is favorable. The Freundlich isotherm which assumes multilayer sorption on the heterogeneous adsorbent surface was also utilized to evaluate experimental data. K F and 1/n were calculated from the intercept and slope of the linear plot of log q e versus log C e . The values of n greater than unity in Freundlich adsorption isotherm, for all temperatures, indicated that the HNC is suitable for the adsorption of MV. The Temkin isotherm model includes a factor that takes into account the interactions between the adsorbent and adsorbate. The heat of molecular adsorption in the layer would reduce due to these interactions and the adsorption is identified by a dissipation of binding energies. The determined Temkin parameters (b T and K T ) exhibited that the interactions between the adsorbent surface and the MV dyes are weak and may be a physical adsorption process. The increased value of b T with temperature also supports increased adsorption efficiency with increasing temperature [74,75].
Dubinin-Radushkevich (D-R) isotherm model determines whether adsorption is physical or chemical. The values of constants q m and β are given in Table 2. The parameter β could be used to calculate the mean free energy (E = 1/ √ 2β), which could distinguish the type of the adsorption process. If the value of E is less than 8 kJ mol −1 , the adsorption process is physical, and when E is between 8 kJ mol −1 and 16 kJ mol −1 , the process is chemical adsorption [76].
The value of apparent energy E (less than 8.0 kJ mol −1 ) indicated that the adsorption of MV onto the HNC is a physical process and is consonant with the ∆G parameter results [77].
The nonlinear relationship of the isotherm models used, is shown in Figure 6.

Adsorption Kinetics
In order to investigate the adsorption process and potential rate controlling step of MV adsorption onto HNC first order, pseudo-first-order, second order, pseudo-secondorder, and intra-particle diffusion kinetic models were studied (Table 1). Figure 7 shows the linear relationship of kinetic models for the MV adsorption process and the value of kinetics parameters determined using these models are given in Ta-

Adsorption Kinetics
In order to investigate the adsorption process and potential rate controlling step of MV adsorption onto HNC first order, pseudo-first-order, second order, pseudo-second-order, and intra-particle diffusion kinetic models were studied (Table 1). Figure 7 shows the linear relationship of kinetic models for the MV adsorption process and the value of kinetics parameters determined using these models are given in Table 3.  On the basis of the high values of the regression coefficient, (R = 1) and very close values of q e,calc with q e,exp , the pseudo-second-order rate model fitted the experimental data better than any other model studied, providing evidence that the adsorption of MV on HNC followed the pseudo-second-order kinetic model. Considerable reduction of dye amount during batch adsorption experiments and longer time needed for adsorbing species to diffuse to remote locations deep within a network of fine pores are the two contributing explanations for these results [79,80].
Similar findings are also reported for adsorption of methylene blue by coconut husks/polylactide blended films [81], adsorption of CV dye on zeolites from coal fly and bottom ashes [78], and adsorption of CV dye on coffee husks [82].
Adsorbate molecules travel from the aqueous phase to the adsorbent surface, then diffuse into pores. External mass transfer (boundary layer diffusion) or intraparticle diffusion determines the adsorption rate. The experimental adsorption kinetics data were fitted using the intraparticle diffusion model and the liquid film diffusion model to determine the rate-determining step.
If the intraparticle diffusion is the sole rate-limiting step, the plot of q t vs. t 0.5 gives a straight line with zero intercept. As the intraparticle diffusion plot did not pass through the origin some other mechanism might be included. For that reason, the kinetic data were examined by the liquid film diffusion model. Linear graphs of ln (1 − F) vs. t with zero intercept suggest that adsorption kinetics are driven by diffusion through the liquid film surrounding solid adsorbents. The low values of R 2 (Table 4) and the non-zero intercept plots of the liquid film diffusion model indicate that the liquid film diffusion is also not the sole rate-determining step. So the adsorption mechanism is a combination of the two processes mentioned above [51,77].

Adsorption Thermodynamics
The free energy change (∆G 0 ), enthalpy change (∆H 0 ), and entropy change (∆S 0 ) were evaluated using the following equations: where Kc is the equilibrium constant (dimensionless), [83,84] R is the gas constant (J/Kmol), and T is the temperature (K).
The ∆H 0 and ∆S 0 were determined from the slope and intercept of the plot of lnK c versus 1/T (Figure 8), and the results are shown in Table 5. where Kc is the equilibrium constant (dimensionless), [83,84] R is the gas constant (J/Kmol), and T is the temperature (K).
The ΔH 0 and ΔS 0 were determined from the slope and intercept of the plot of lnKc versus 1/T (Figure 8), and the results are shown in Table 5. The positive value of ΔH 0 (+37.84 kJ/mol) indicated that the adsorption process was endothermic and since this value is less than 80 kJ/mol (80 kJ/mol is the upper limit of the change of enthalpy for physisorption), it shows that the adsorption follows a physisorption.
The positive value of ΔS 0 indicated increased disorder and randomness at the solidliquid interface; however, the negative value of ΔG 0 suggested spontaneous and thermodynamically favorable adsorption of MV onto HNC. Meantime, the value of ΔG 0 became more negative with increasing temperature, which exhibits that higher temperature is contributory to adsorption [19].
The activation energy Ea was obtained using the Arrhenius equation: where k2 is the pseudo-second-order rate constant, A is the Arrhenius constant, Ea refers to the energy of activation (J mol −1 ), R is the ideal gas constant (8.314 J mol −1 K −1 ), and T is the temperature (K). The slope of the plot of ln k2 versus 1/T (Figure 9) is used to evaluate Ea, which was found to be 44 kJ/mol.  The positive value of ∆H 0 (+37.84 kJ/mol) indicated that the adsorption process was endothermic and since this value is less than 80 kJ/mol (80 kJ/mol is the upper limit of the change of enthalpy for physisorption), it shows that the adsorption follows a physisorption.
The positive value of ∆S 0 indicated increased disorder and randomness at the solidliquid interface; however, the negative value of ∆G 0 suggested spontaneous and thermodynamically favorable adsorption of MV onto HNC.
Meantime, the value of ∆G 0 became more negative with increasing temperature, which exhibits that higher temperature is contributory to adsorption [19].
The activation energy Ea was obtained using the Arrhenius equation: where k 2 is the pseudo-second-order rate constant, A is the Arrhenius constant, E a refers to the energy of activation (J mol −1 ), R is the ideal gas constant (8.314 J mol −1 K −1 ), and T is the temperature (K).
The slope of the plot of ln k 2 versus 1/T (Figure 9) is used to evaluate E a , which was found to be 44 kJ/mol.  Figure 10 depicts the lowest energy configurations for the MV on the HNC surface under the simulated adsorption circumstances that were employed (as chosen above). As evidenced by Mulliken charges (Figure 1), the adsorption geometries of the MV demonstrate that nitrogen atoms are mostly responsible for this behavior. The following equation may be used to calculate the adsorption energies of the adsorbate molecule on the halloysite surface.

Monte Carlo and Molecular Dynamic Simulations
where: Etotal is the total energy of the system as a result of adsorbent-adsorbate; Esurface+water and EMV+water is system energy in the absence and presence of MV. This value is consistent with the values of the activation energy (43.0 kJ/mol) for the adsorption of reactive red 189 on cross-linked chitosan beads [85], and (5.6-49.1 kJ/mol) for the adsorption of polychlorinated biphenyls on fly ash [86]. Figure 10 depicts the lowest energy configurations for the MV on the HNC surface under the simulated adsorption circumstances that were employed (as chosen above). As evidenced by Mulliken charges (Figure 1), the adsorption geometries of the MV demonstrate that nitrogen atoms are mostly responsible for this behavior.  Figure 10 depicts the lowest energy configurations for the MV on the under the simulated adsorption circumstances that were employed (as cho evidenced by Mulliken charges (Figure 1), the adsorption geometries of th strate that nitrogen atoms are mostly responsible for this behavior. The following equation may be used to calculate the adsorption energies of the adsorbate molecule on the halloysite surface. E ads = E total − E sur f ace+water + E MV+sur f ace + E water (9) where: E total is the total energy of the system as a result of adsorbent-adsorbate; E surface+water and E MV+water is system energy in the absence and presence of MV. Figure 11 depicts the distribution of the adsorption energies for the large number of adsorptive configurations developed and computed by the Monte Carlo approach for the MV adsorbate, which was based on the results of the Monte Carlo method.  The Adsorption Energy (Eads) values for the adsorption of MV onto the HNC surface are presented in Table 6. As indicated by the high Eads, the MV has a strong contact with the halloysite surface, resulting in good adsorptive properties. Adsorption of MV onto the halloysite surface was determined by both the MC and MD calculations, respectively. In the case of the adsorption process, Monte Carlo simulations (see Figures 10 and 12) have indicated that the ensuing negative values of adsorption are indicative of the spontaneity with which the process occurs (Figure 11). It is important to employ MD simulations because they give a very simple means of tracking and recording the kinetics of inhibitor adsorption on the metal surface under study. Figure 12 shows the final arrangements of MV on the halloysite surface as depicted in the previous figure. The Adsorption Energy (E ads ) values for the adsorption of MV onto the HNC surface are presented in Table 6. As indicated by the high E ads , the MV has a strong contact with the halloysite surface, resulting in good adsorptive properties. Adsorption of MV onto the halloysite surface was determined by both the MC and MD calculations, respectively. In the case of the adsorption process, Monte Carlo simulations (see Figures 10 and 12) have indicated that the ensuing negative values of adsorption are indicative of the spontaneity with which the process occurs ( Figure 11). It is important to employ MD simulations because they give a very simple means of tracking and recording the kinetics of inhibitor adsorption on the metal surface under study. Figure 12 shows the final arrangements of MV on the halloysite surface as depicted in the previous figure. The presence of an RDF peak between 1 and 3.5 Å from the surface plane terial and the adsorbate atom confirms adsorption, but for physisorption, RD greater distances confirm physisorption, and vice versa.

Monte Carlo and Molecular Dynamic Simulations
The RDF of the nitrogen atoms in the MV (Figure 13 is close to the sur suggesting that these components have an increased interaction with the adso face. The findings of the MD and RDF analyses confirm that these inhibitors h proclivity to adsorb onto the surface, owing to their peculiar attraction to share electrons with the surface in question. [55].  The presence of an RDF peak between 1 and 3.5 Å from the surface plane of the material and the adsorbate atom confirms adsorption, but for physisorption, RDF peaks at greater distances confirm physisorption, and vice versa.
The RDF of the nitrogen atoms in the MV (Figure 13 is close to the surface plane, suggesting that these components have an increased interaction with the adsorbent surface. The findings of the MD and RDF analyses confirm that these inhibitors have a high proclivity to adsorb onto the surface, owing to their peculiar attraction to share and accept electrons with the surface in question. [55]. The presence of an RDF peak between 1 and 3.5 Å from the surface plane of the material and the adsorbate atom confirms adsorption, but for physisorption, RDF peaks at greater distances confirm physisorption, and vice versa.
The RDF of the nitrogen atoms in the MV (Figure 13 is close to the surface plane, suggesting that these components have an increased interaction with the adsorbent surface. The findings of the MD and RDF analyses confirm that these inhibitors have a high proclivity to adsorb onto the surface, owing to their peculiar attraction to share and accept electrons with the surface in question. [55].

FTIR Spectroscopy
FTIR-8400S was used for Total Attenuated Reflection (ATR) measurements with the following parameters: resolution 2 cm −1 , 100 scans, 500-4000 cm −1 . Figure 14 exhibits the peaks that appeared at certain wave numbers of MV FTIR spectra before and after adsorption onto HNC.

FTIR Spectroscopy
FTIR-8400S was used for Total Attenuated Reflection (ATR) measurements with the following parameters: resolution 2 cm −1 , 100 scans, 500-4000 cm −1 . Figure 14 exhibits the peaks that appeared at certain wave numbers of MV FTIR spectra before and after adsorption onto HNC.
The following peaks: 3082 and 1496 cm −1 show -CH-aromatic and 2927 cm −1 -CH3 methyl vibrating vibrations, respectively 1688 cm −1 , aromatic ring vibrations -C=C-vibrations from the skeleton of aromatic ring structures of MV, while we have new displays of vibrations of the group -OH at 3345 cm −1 stretching, 1127.2 cm −1 bending by the group connections -Si-O-Si-bending and 990.3 cm −1 from the -Al-O-OH-bond arising from the Halloysite and MV groups with distinctive features (as shown in Figure 14). The interaction of MV molecules with Halloysite functional groups was confirmed by the emergence and disappearance or depletion of various peaks.

Conclusions
This study established the viability of halloysite nanoclay as a dye-removal adsorbent in wastewater. The findings of the batch trials revealed that adsorption was affected by adsorbent dose, contact time, initial dye concentration, and temperature.
Adsorption of MV is very little affected by changing the temperature and pH. The optimal adsorption capacity (27.7 mg/g) was reached at ambient temperature (298 K) and pH (4.26).
Among the investigated isotherm models, the Langmuir isotherm model demonstrated the best fit.
The experimental data were best suited by the pseudo-second-order rate model, indicating that the adsorption of MV on HNC followed the pseudo-second-order kinetic model. The kinetic analysis indicated that the adsorption mechanism is a combination of intra-particle diffusion and film diffusion.
Thermodynamic investigations demonstrated that the removal of MV by HNC is possible, exothermic, and spontaneous.  Figure 14).
The interaction of MV molecules with Halloysite functional groups was confirmed by the emergence and disappearance or depletion of various peaks.

Conclusions
This study established the viability of halloysite nanoclay as a dye-removal adsorbent in wastewater. The findings of the batch trials revealed that adsorption was affected by adsorbent dose, contact time, initial dye concentration, and temperature.
Adsorption of MV is very little affected by changing the temperature and pH. The optimal adsorption capacity (27.7 mg/g) was reached at ambient temperature (298 K) and pH (4.26).
Among the investigated isotherm models, the Langmuir isotherm model demonstrated the best fit.
The experimental data were best suited by the pseudo-second-order rate model, indicating that the adsorption of MV on HNC followed the pseudo-second-order kinetic model. The kinetic analysis indicated that the adsorption mechanism is a combination of intra-particle diffusion and film diffusion.
Thermodynamic investigations demonstrated that the removal of MV by HNC is possible, exothermic, and spontaneous.
This study focuses on two distinct molecular mechanics-based theoretical approaches (MC and MD). These techniques enabled a molecular comprehension of the interaction between the MV molecule and the halloysite surface. Simulations using Monte Carlo and molecular dynamics revealed that this molecule spontaneously adsorbs onto surfaces. The negative value of the adsorption energies supports a strong interaction between the MV and HNC surface, which is consistent with experimental findings.