Adsorption of Malachite Green Dye onto Mesoporous Natural Inorganic Clays: Their Equilibrium Isotherm and Kinetics Studies

: Contamination of water with organic dyes is a major environmental concern as it causes serious life-threatening environmental problems. The present research was designed to evaluate the potential of three different natural inorganic clays (NICs) i.e., Pakistani bentonite clay (PB), bentonite purchased from Alfa Aesar (BT), and Turkish red mud (RM) for malachite green (MG) dye removal from an aqueous solution. Various analytical techniques, namely X-ray ﬂuorescence spectrometry (XRF), X-ray diffractometry (XRD), Fourier transform infrared spectroscopy (FTIR), ﬁeld emission scanning electron microscopy (FESEM), Brunauer–Emmett–Teller surface area measurement (BET), and thermogravimetric analysis (TGA), were used to investigate the physicochemical properties of the NICs samples. The effect of adsorption operational parameters such as contact time, aqueous phase pH, dye concentration, and amount of NICs on the adsorption behavior of MG onto NICs samples were investigated under the batch adsorption system. The equilibrium and kinetic inspection reﬂected the best description of MG adsorption behavior by the Langmuir isotherm model and pseudo-ﬁrst-order kinetic model, respectively. The results indicated that the adsorption was favorable at higher pH. The maximum adsorption capacities calculated by Langmuir isotherm for PB, BT, and RM were found to be 243.90 mg/g, 188.68 mg/g, and 172.41 mg/g, respectively. It can be concluded that natural inorganic clays with a higher surface area can be used as an effective adsorbent material to remove the MG dye from an aqueous solution. Bruker Tiger S8 XRF Spectrometer. The X-ray diffraction (XRD) patterns using Bruker D-2 Phaser, Cu K α radiation over a 2 θ interval of 4 ◦ –40 ◦ were performed to investigate the crystallinity and phase composition of the NICs adsorbents. A Fourier transform infrared spectroscopic (FTIR) analysis was performed by using Bruker-Tensor-27 between 400 to 4000 cm − 1 for the examination of surface functional groups of NICs. A ﬁeld emission scanning electron microscopic (FESEM) analysis was performed to determine the NICs adsorbent morphologies using FESEM-EDS (Zeiss Ultra Plus). The Brunauer–Emmett–Teller (BET) nitrogen physisorption measurements were performed using a BET Micromeritics ASAP 2020 instrument for the calculations of pore sizes and speciﬁc surface area of the NICs adsorbents. The thermogravimetric analysis (TGA) measurements were carried out using the TGA Q500 model instrument, a nitrogen atmosphere (60 mL/min of N 2 ), and the samples were heated at 10 ◦ C per minute rise over the temperature range of 25–800 ◦ C.


Introduction
In recent years, planners, environmental scientists, and decision-makers have been paying attention to sustainable resource development [1]. Within these resources, water is the most precious renewable natural resource, is the essential enabler, and a major source of survival of life [2,3]. Groundwater is the most important source of drinking water in the world [4]. Globally, groundwater accounts for about 43% of total irrigation use and provides potable water for about 1.5 billion people [5]. However, the rapid development of industrialization, manifold increase in the human population growth, and uncontrolled usage of ground and washed with boiled distilled water followed by a filtration process to remove the soluble impurities. After that, NICs samples were then dried in an oven for 24 h at 80 • C and ground. Sodium hydroxide, hydrochloric acid, and malachite green is a water-soluble cationic dye that belongs to the triphenylmethane category (chemical formula: C 23 H 25 ClN 2 , color index number = 42,000, molar mass: 364.91 g/mol, abbreviated as MG) were obtained from Sigma Aldrich, and used without further pre-treatment. The structure of the MG is shown in Figure 1.
Water 2021, 13, x FOR PEER REVIEW 3 of 17

Materials
The NICs adsorbents were obtained: bentonite from Pakistan (PB), bentonite purchased from Alfa Aesar (BT), and red mud from Turkey (RM). The NICs samples were ground and washed with boiled distilled water followed by a filtration process to remove the soluble impurities. After that, NICs samples were then dried in an oven for 24 h at 80 °C and ground. Sodium hydroxide, hydrochloric acid, and malachite green is a watersoluble cationic dye that belongs to the triphenylmethane category (chemical formula: C23H25ClN2, color index number = 42,000, molar mass: 364.91 g/mol, abbreviated as MG) were obtained from Sigma Aldrich, and used without further pre-treatment. The structure of the MG is shown in Figure 1.

Pretreatment of NICs Samples
All of the NICs samples were pretreated as follows: the NICs samples were suspended in distilled water for 4 h and cleaned many times with freshly prepared distilled water. The clean NICs were then oven-dried for 24 h at 105 °C. Finally, the dried NICs materials were stored in sealed jars for further use, without any further physical or chemical treatment.

Characterization of NICs
To evaluate the mechanism of dye adsorption, it is imperative to examine the characteristics of the adsorbent materials. Therefore, the physicochemical characterizations of the adsorbent materials were analyzed by using various analytical techniques. The chemical composition and elemental analysis of NICs were determined by X-ray fluorescence spectroscopy (XRF) using Bruker Tiger S8 XRF Spectrometer. The X-ray diffraction (XRD) patterns using Bruker D-2 Phaser, Cu Kα radiation over a 2θ interval of 4°-40° were performed to investigate the crystallinity and phase composition of the NICs adsorbents. A Fourier transform infrared spectroscopic (FTIR) analysis was performed by using Bruker-Tensor-27 between 400 to 4000 cm −1 for the examination of surface functional groups of NICs. A field emission scanning electron microscopic (FESEM) analysis was performed to determine the NICs adsorbent morphologies using FESEM-EDS (Zeiss Ultra Plus). The Brunauer-Emmett-Teller (BET) nitrogen physisorption measurements were performed using a BET Micromeritics ASAP 2020 instrument for the calculations of pore sizes and specific surface area of the NICs adsorbents. The thermogravimetric analysis (TGA) measurements

Pretreatment of NICs Samples
All of the NICs samples were pretreated as follows: the NICs samples were suspended in distilled water for 4 h and cleaned many times with freshly prepared distilled water. The clean NICs were then oven-dried for 24 h at 105 • C. Finally, the dried NICs materials were stored in sealed jars for further use, without any further physical or chemical treatment.

Characterization of NICs
To evaluate the mechanism of dye adsorption, it is imperative to examine the characteristics of the adsorbent materials. Therefore, the physicochemical characterizations of the adsorbent materials were analyzed by using various analytical techniques. The chemical composition and elemental analysis of NICs were determined by X-ray fluorescence spectroscopy (XRF) using Bruker Tiger S8 XRF Spectrometer. The X-ray diffraction (XRD) patterns using Bruker D-2 Phaser, Cu Kα radiation over a 2θ interval of 4 • -40 • were performed to investigate the crystallinity and phase composition of the NICs adsorbents. A Fourier transform infrared spectroscopic (FTIR) analysis was performed by using Bruker-Tensor-27 between 400 to 4000 cm −1 for the examination of surface functional groups of NICs. A field emission scanning electron microscopic (FESEM) analysis was performed to determine the NICs adsorbent morphologies using FESEM-EDS (Zeiss Ultra Plus). The Brunauer-Emmett-Teller (BET) nitrogen physisorption measurements were performed using a BET Micromeritics ASAP 2020 instrument for the calculations of pore sizes and specific surface area of the NICs adsorbents. The thermogravimetric analysis (TGA) measurements were carried out using the TGA Q500 model instrument, a nitrogen atmosphere (60 mL/min of N 2 ), and the samples were heated at 10 • C per minute rise over the temperature range of 25-800 • C.

Batch Adsorption Studies
To evaluate the adsorption abilities of NICs adsorbents, batch adsorption experiments were performed using 250 mL conical flasks containing 100 mL of working solutions. The influence of operating parameters such as time (10-120 min), NICs amount (0.1-0.5 g), dye concentration (100-350 mg/L), and pH (3-11) on MG dye remediation was investigated in batch studies. Typically, 0.1 g of NICs and MG solution (100 mL) with the desired pH and concentration were added to a 250 mL conical flask and stirred continuously. At pre-determined intervals, the NICs samples were separated. The remaining concentration of MG in the filtrate was spectrophotometrically measured by a UV-vis spectrophotometer (Shimadzu UV1700 Japan).
The adsorbed amount of MG (mg/g) was calculated using the following equations: where V (L) represents the volume of MG solution, C o (mg/L) is the initial MG concentration, C e (mg/L) is the concentration at equilibrium, and M (g) is the amount of NICs.

Characterization of NICs
The XRF analysis was carried out to investigate the chemical compositions of the NICs used adsorbents. Table 1 shows that PB and BT clays are mainly composed of alumina and silica, while the major components of the RM clay are iron oxide, alumina, and silica, with other oxides present in trace amount in all of the NICs samples. Figure 2 shows the major chemical constituents present in NICs samples. were carried out using the TGA Q500 model instrument, a nitrogen atmosphere (60 mL/min of N2), and the samples were heated at 10 °C per minute rise over the temperature range of 25-800 °C.

Batch Adsorption Studies
To evaluate the adsorption abilities of NICs adsorbents, batch adsorption experiments were performed using 250 mL conical flasks containing 100 mL of working solutions. The influence of operating parameters such as time (10-120 min), NICs amount (0.1-0.5 g), dye concentration (100-350 mg/L), and pH (3-11) on MG dye remediation was investigated in batch studies. Typically, 0.1 g of NICs and MG solution (100 mL) with the desired pH and concentration were added to a 250 mL conical flask and stirred continuously. At pre-determined intervals, the NICs samples were separated. The remaining concentration of MG in the filtrate was spectrophotometrically measured by a UV-vis spectrophotometer (Shimadzu UV1700 Japan).
The adsorbed amount of MG (mg/g) was calculated using the following equations: where V (L) represents the volume of MG solution, Co (mg/L) is the initial MG concentration, Ce (mg/L) is the concentration at equilibrium, and M (g) is the amount of NICs.

Characterization of NICs
The XRF analysis was carried out to investigate the chemical compositions of the NICs used adsorbents. Table 1 shows that PB and BT clays are mainly composed of alumina and silica, while the major components of the RM clay are iron oxide, alumina, and silica, with other oxides present in trace amount in all of the NICs samples. Figure 2 shows the major chemical constituents present in NICs samples.

Parameter
Chemical Composition (%) PB Clay BT Clay RM Clay XRD analysis was employed to examine the mineralogical composition and crystalline nature of the NICs samples, as shown in Figure 3. The results indicate that PB and BT clay has the following mineral phases: illite, kaolinite, quartz, and calcite. The dominant diffraction peaks for the PB and BT clays were found at Bragg's angle (2θ) =~9 • ,~12.5 • ,~27 • , 28 • , and~29 • which corresponds to illite, kaolinite, quartz, feldspars, and calcite, respectively [43,44]. According to the XRD data, RM is mainly composed of hematite, with other minerals such as gibbsite, diaspore, and calcite being present as minor constituents [45]. K2O 3. XRD analysis was employed to examine the mineralogical composition and crystalline nature of the NICs samples, as shown in Figure 3. The results indicate that PB and BT clay has the following mineral phases: illite, kaolinite, quartz, and calcite. The dominant diffraction peaks for the PB and BT clays were found at Bragg's angle (2θ) = ~9°, ~12.5°, ~27°, ~28°, and ~29° which corresponds to illite, kaolinite, quartz, feldspars, and calcite, respectively [43,44]. According to the XRD data, RM is mainly composed of hematite, with other minerals such as gibbsite, diaspore, and calcite being present as minor constituents [45].  [46]. The broadbands 3386-3422 cm −1 were due to the stretching vibration of H-O-H of hydrogen-bonded inter-layer water molecules [47], while the bands of OH deformation mode of coordinated water molecules appeared at 1632-1641 cm −1 [48]. The strong bands at 970-998 cm −1 were referred to as the Si-O-Si stretching vibrations [49].   [46]. The broadbands 3386-3422 cm −1 were due to the stretching vibration of H-O-H of hydrogen-bonded inter-layer water molecules [47], while the bands of OH deformation mode of coordinated water molecules appeared at 1632-1641 cm −1 [48]. The strong bands at 970-998 cm −1 were referred to as the Si-O-Si stretching vibrations [49]. The NICs adsorbent morphologies were analyzed by SEM analysis ( Figure 5). It was observed that PB and BT clays have quite a rough porous surface with blunt edges due to the agglomeration of small particle size, while the RM clay has a smooth surface compared to the other two clays and a fluffy appearance because of the closely packed flakes. The NICs adsorbent morphologies were analyzed by SEM analysis ( Figure 5). It was observed that PB and BT clays have quite a rough porous surface with blunt edges due to the agglomeration of small particle size, while the RM clay has a smooth surface compared to the other two clays and a fluffy appearance because of the closely packed flakes. The NICs adsorbent morphologies were analyzed by SEM analysis ( Figure 5). It was observed that PB and BT clays have quite a rough porous surface with blunt edges due to the agglomeration of small particle size, while the RM clay has a smooth surface compared to the other two clays and a fluffy appearance because of the closely packed flakes. The pore size distribution and surface area of the NICs materials were determined by BET N 2 sorption isotherm analysis. Figure 6 shows that the observed sorption isotherms are of type IV with an H3 type of hysteresis loops. The results confirm the mesoporous nature of the used NICs adsorbents [52]. Both the adsorption-desorption isotherms completely overlapped at low relative pressure P/P o < 0.4, but showed a distinct hysteresis loop at relatively high pressure (P/P o > 0.4), which is the typical characteristics of layered materials. The calculated pore sizes and surface area derived from BET N 2 sorption isotherm of the NICs are given in Table 2. The thermal behaviors of NICs were studied by TG analysis. It is clear from Figure 7 that a continuous weight loss was observed for all the clay minerals in the test temperature range. A two-step weight loss of 9.8%, 8.7%, and 7.5% was observed for PB, BT, and RM adsorbents, respectively. The first weight loss due to the dehydration of adsorbed water occurred over the 30 • C to 160 • C temperature range, while the second weight loss over the 200 • C to 650 • C temperature range occurred due to the loss of hydrated cations on the exchangeable sites and interlayer water dehydration of the NICs adsorbents. The pore size distribution and surface area of the NICs materials were determined by BET N2 sorption isotherm analysis. Figure 6 shows that the observed sorption isotherms are of type IV with an H3 type of hysteresis loops. The results confirm the mesoporous nature of the used NICs adsorbents [52]. Both the adsorption-desorption isotherms completely overlapped at low relative pressure P/Po < 0.4, but showed a distinct hysteresis loop at relatively high pressure (P/Po > 0.4), which is the typical characteristics of layered materials. The calculated pore sizes and surface area derived from BET N2 sorption isotherm of the NICs are given in Table 2. The thermal behaviors of NICs were studied by TG analysis.  It is clear from Figure 7 that a continuous weight loss was observed for all the clay minerals in the test temperature range. A two-step weight loss of 9.8%, 8.7%, and 7.5% was observed for PB, BT, and RM adsorbents, respectively. The first weight loss due to the dehydration of adsorbed water occurred over the 30 °C to 160 °C temperature range, while the second weight loss over the 200 °C to 650 °C temperature range occurred due to the loss of hydrated cations on the exchangeable sites and interlayer water dehydration of the NICs adsorbents.

Adsorption Studies of MG Dye
The adsorption behaviors of the NICs were examined to assess their potential applications in the decontamination of wastewater containing organic dyes. The influence of adsorption time, aqueous phase pH, MG concentration, and the amount of NICs on the decolorization capabilities of NICs were evaluated under batch conditions.
The adsorption uptake of MG by the NICs under varied periods of adsorption contact time (10 to 120 min) is shown in Figure 8. During the experiments, the adsorption capacities increased with the increase in the adsorption contact time. Initially, the rate of dye uptake was high, and afterward, a gradual decrease was observed until the equilibrium was attained [53]. This trend of variation in the MG uptake rates was related to the availability of adsorbent surface active sites. In the beginning, the faster adsorption rates were attributed to the abundant surface-active sites that saturated over time and opposed further adsorption [54,55].

Adsorption Studies of MG Dye
The adsorption behaviors of the NICs were examined to assess their potential applications in the decontamination of wastewater containing organic dyes. The influence of adsorption time, aqueous phase pH, MG concentration, and the amount of NICs on the decolorization capabilities of NICs were evaluated under batch conditions.
The adsorption uptake of MG by the NICs under varied periods of adsorption contact time (10 to 120 min) is shown in Figure 8. During the experiments, the adsorption capacities increased with the increase in the adsorption contact time. Initially, the rate of dye uptake was high, and afterward, a gradual decrease was observed until the equilib-rium was attained [53]. This trend of variation in the MG uptake rates was related to the availability of adsorbent surface active sites. In the beginning, the faster adsorption rates were attributed to the abundant surface-active sites that saturated over time and opposed further adsorption [54,55].
The adsorption behaviors of the NICs were examined to assess their potential applications in the decontamination of wastewater containing organic dyes. The influence of adsorption time, aqueous phase pH, MG concentration, and the amount of NICs on the decolorization capabilities of NICs were evaluated under batch conditions.
The adsorption uptake of MG by the NICs under varied periods of adsorption contact time (10 to 120 min) is shown in Figure 8. During the experiments, the adsorption capacities increased with the increase in the adsorption contact time. Initially, the rate of dye uptake was high, and afterward, a gradual decrease was observed until the equilibrium was attained [53]. This trend of variation in the MG uptake rates was related to the availability of adsorbent surface active sites. In the beginning, the faster adsorption rates were attributed to the abundant surface-active sites that saturated over time and opposed further adsorption [54,55].  The initial solution pH is a foremost significant controlling parameter in the process of dye adsorption, as it can directly influence the adsorbent surface charges, the solute's dissociation, and as the adsorption mechanism [56]. It is also associated with the structural changes and color intensities of dye molecules and is directly related to the competitive adsorption process between the MG molecules and hydrogen ions [57]. The influence of initial pH on the adsorption characteristics of MG by NICs adsorbents was examined in the pH range of 3-11. According to the results (Figure 9), the adsorption capacities of NICs adsorbents increased from 45.6 to 130.8 mg/g for RM, 62.4 to 142.8 mg/g for BT, and 82.8 to 170.6 mg/g for PB, with the increase of pH of MG solution from 3 to 11, respectively. In an alkaline environment, the surface functional groups (Si-OH and Al-OH) of the NICs were deprotonated by -OH, due to which the surface became negatively charged; as a result, the electrostatic interaction between anionic NICs surfaced and cationic MG molecules increased [1], whereas at low solution pH, the H+ ions concentration was very high, thus competing intensely with the cationic dye to occupy the adsorption site, therefore causing a reduction in the adsorption of MG dye [58].
It is essential to analyze the effect of the adsorbent amount, to optimize and select the best-required dose of an adsorbent for scaling-up and designing large-scale equipment [59]. The influence of adsorbent dosage on dye adsorption is presented in Figure 10. The obtained results show that increasing the NICs dose from 0.1 to 0.5 g resulted in decreases in the adsorption capacities from 76 to 23.4 mg/g for RM, 92 to 27.4 mg/g for BT, and 134 to 29.6 mg/g for PB. The resulted decrease in adsorption capacities is due to the increasing NICs dose (0.1 to 0.5 g), while the number of MG molecules remained fixed (Ce = 100 mg/L). Hence, at a higher NICs dose, some of the surface-active sites of the NICs remained empty [60].
To determine the influence of the initial MG concentration on the adsorption capacities, the experiments were performed by changing the concentration in the range of 100-350 mg/L, keeping other conditions constant. The results are shown in Figure 11. It was determined that the adsorption capacities of MG onto NICs adsorbents increased from 92 to 223 mg/g for PB, 69 to 149 mg/g for BT, and 57 to 125 mg/g for RM as the initial concentration increased from 100 to 350 mg/L. Moreover, at higher concentrations, the adsorption uptake rate was also higher due to the higher concentration gradient [61].
to 170.6 mg/g for PB, with the increase of pH of MG solution from 3 to 11, respectively. In an alkaline environment, the surface functional groups (Si-OH and Al-OH) of the NICs were deprotonated by -OH, due to which the surface became negatively charged; as a result, the electrostatic interaction between anionic NICs surfaced and cationic MG molecules increased [1], whereas at low solution pH, the H+ ions concentration was very high, thus competing intensely with the cationic dye to occupy the adsorption site, therefore causing a reduction in the adsorption of MG dye [58]. It is essential to analyze the effect of the adsorbent amount, to optimize and select the best-required dose of an adsorbent for scaling-up and designing large-scale equipment [59]. The influence of adsorbent dosage on dye adsorption is presented in Figure 10. The obtained results show that increasing the NICs dose from 0.1 to 0.5 g resulted in decreases in the adsorption capacities from 76 to 23.4 mg/g for RM, 92 to 27.4 mg/g for BT, and 134 to 29.6 mg/g for PB. The resulted decrease in adsorption capacities is due to the increasing NICs dose (0.1 to 0.5 g), while the number of MG molecules remained fixed (Ce = 100 mg/L). Hence, at a higher NICs dose, some of the surface-active sites of the NICs remained empty [60]. To determine the influence of the initial MG concentration on the adsorption capacities, the experiments were performed by changing the concentration in the range of 100-350 mg/L, keeping other conditions constant. The results are shown in Figure 11. It was determined that the adsorption capacities of MG onto NICs adsorbents increased from 92 to 223 mg/g for PB, 69 to 149 mg/g for BT, and 57 to 125 mg/g for RM as the initial concentration increased from 100 to 350 mg/L. Moreover, at higher concentrations, the adsorption uptake rate was also higher due to the higher concentration gradient [61]. To determine the influence of the initial MG concentration on the adsorption capacities, the experiments were performed by changing the concentration in the range of 100-350 mg/L, keeping other conditions constant. The results are shown in Figure 11. It was determined that the adsorption capacities of MG onto NICs adsorbents increased from 92 to 223 mg/g for PB, 69 to 149 mg/g for BT, and 57 to 125 mg/g for RM as the initial concentration increased from 100 to 350 mg/L. Moreover, at higher concentrations, the adsorption uptake rate was also higher due to the higher concentration gradient [61]. Figure 11. Effect of initial dye concentration on MG adsorption. By taking into account the above experimental data and analysis, a mechanism for the electrostatic interactions between the MG and NICs was proposed, as illustrated in Figure 12. By taking into account the above experimental data and analysis, a mechanism for the electrostatic interactions between the MG and NICs was proposed, as illustrated in Figure 12.
Water 2021, 13, x FOR PEER REVIEW 11 of 17 Figure 12. The adsorption mechanism of MG onto NICs.

Adsorption Kinetics
Adsorption kinetics analysis was carried out to describe the rate of adsorbate uptake and evaluate the adsorption mechanism. The adsorption process occurs mainly in three steps; (i) the external mass transfer of the adsorbate molecules from the bulk solution to the external adsorbent surface, (ii) transfer of adsorbed molecules to the adsorption sites, (iii) and finally, retention via sorption itself [62]. Therefore, the experimental data were inspected by the two different kinetic models, namely the pseudo-first-order and pseudosecond-order kinetics model.
The pseudo-first-order kinetic model can be represented as follows [55]: where k1 (min −1 ) is the pseudo-first-order rate constant, qe and qt are the adsorption capacities (mg/g) at equilibrium and at any time t (min), respectively. The empirical equation for the pseudo-second-order kinetic model is [63]: Integrating and rearranging Equation (3), we received = 1 + 1 (4) where k2 (g/mg. min) is the pseudo-second-order model rate constant. The linear fitting kinetic results are presented in Figure 13a,b and the calculated kinetic parameters are given in Table 3. When the regression coefficients (R 2 ) values of the used models were compared, it was found that they have almost the same values for R 2 . However, the experimental (qexp) values of the pseudo-first-order model were found close to the calculated (qe) as compared to the pseudo-second-order model, which confirms that the adsorption of MG was followed by the pseudo-first-order model.

Adsorption Kinetics
Adsorption kinetics analysis was carried out to describe the rate of adsorbate uptake and evaluate the adsorption mechanism. The adsorption process occurs mainly in three steps; (i) the external mass transfer of the adsorbate molecules from the bulk solution to the external adsorbent surface, (ii) transfer of adsorbed molecules to the adsorption sites, (iii) and finally, retention via sorption itself [62]. Therefore, the experimental data were inspected by the two different kinetic models, namely the pseudo-first-order and pseudo-second-order kinetics model.
The pseudo-first-order kinetic model can be represented as follows [55]: where k 1 (min −1 ) is the pseudo-first-order rate constant, q e and q t are the adsorption capacities (mg/g) at equilibrium and at any time t (min), respectively. The empirical equation for the pseudo-second-order kinetic model is [63]: Integrating and rearranging Equation (3), we received where k 2 (g/mg. min) is the pseudo-second-order model rate constant. The linear fitting kinetic results are presented in Figure 13a,b and the calculated kinetic parameters are given in Table 3. When the regression coefficients (R 2 ) values of the used models were compared, it was found that they have almost the same values for R 2 . However, the experimental (q exp ) values of the pseudo-first-order model were found close to the calculated (q e ) as compared to the pseudo-second-order model, which confirms that the adsorption of MG was followed by the pseudo-first-order model.  A direct comparison of the adsorption capacities obtained in this study and earlier reported in the literature is difficult due to the varying experimental conditions employed in those studies. However, from Table 5, it can be concluded that the NICs used in this study provide better results for MG adsorption than the others reported in the literature.

Adsorption Isotherm
The adsorption isotherm studies provide helpful information to understand the nature of the interaction between the adsorbed matter and adsorbent, and to evaluate the efficiency of the adsorbent material used for adsorption. The experimental data were evaluated by two commonly used adsorption isotherm models, namely the Freundlich and Langmuir models. According to the Langmuir model, adsorption of analyte takes place on the homogenous sites the adsorbent with the monolayer formation [63], and can be expressed as: where K L (L/mg) represents the adsorption equilibrium constant for the Langmuir model and is related to the adsorption energy, q e and q m are the equilibrium adsorption amount of adsorbate and monolayer adsorption capacity of adsorbents in the experiments (mg/g), respectively, and Ce (mg/L) is equilibrium adsorbate concentration.
The Freundlich isotherm, applied to non-ideal and reversible adsorption, is valid for multilayer adsorption of analyte and can be expressed as: l n q e = l n K F + 1 n F l n C e (6) where K F (mg/g) represents the Freundlich constant and 1/n is the heterogeneity factor. All of the correlation coefficients and constants obtained from the used adsorption isotherm models were summarized in Table 4. The results show that both of the isotherms have very high R 2 ; however, experimental maximum adsorption capacities are close to the q m calculated from the Langmuir isotherm model. The results show that adsorption of MG onto the surface of NICs adsorbents occurs via monolayer formation. A direct comparison of the adsorption capacities obtained in this study and earlier reported in the literature is difficult due to the varying experimental conditions employed in those studies. However, from Table 5, it can be concluded that the NICs used in this study provide better results for MG adsorption than the others reported in the literature.

Conclusions
In the current study, three different mesoporous natural clays were used as an environment-friendly, efficient, easily available, and low-cost material for the adsorption of cationic MG dye from an aqueous environment. The chemical composition, mineralogical composition, and texture of NICs were determined by using XRF, XRD, FESEM, FTIR, BET, and TGA analysis. It was observed that PB clay has a quite rough surface with blunt edges and higher surface area as compared to the BT clay and RM clay. The experimental results indicated that adsorption removal of MG was highly dependent on the adsorption contact time, aqueous pH, NICs dose, and MG concentration. The adsorption capacities of MG onto NICs increased with the increase in the adsorption contact time, pH, and MG concentration, while it decreased with the increase in NICs dose. Furthermore, the kinetics data fitted the pseudo-first-order model well. The isotherm data illustrated the suitability of employing the Langmuir isotherm model. The adsorption capacity of PB clay (223 mg/g) was found higher than the BT clay (149 mg/g) and RM clay (125 mg/g). The adsorption of MG dye was more correlated to the BET surface area and available binding sites of the