Zinc Oxide Nanoparticles and Their Application in Adsorption of Toxic Dye from Aqueous Solution

Dye waste is one of the most serious types of pollution in natural water bodies, since its presence can be easily detected by the naked eye, and it is not easily biodegradable. In this study, zinc oxide nanoparticles (ZnO-NPs) were generated using a chemical reduction approach involving the zinc nitrate procedure. Fourier transform infrared (FTIR), scanning electron microscopy (SEM), Brunauer-Emmett-Teller (BET), and UV-vis techniques were used to analyse the surface of ZnO-NPs. The results indicate the creation of ZnO-NPs with a surface area of 95.83 m2 g−1 and a pore volume of 0.058 cm3 g−1, as well as an average pore size of 1.22 nm. In addition, the ZnO-NPs were used as an adsorbent for the removal of Ismate violet 2R (IV2R) dye from aqueous solutions under various conditions (dye concentration, pH, contact time, temperature, and adsorbent dosage) using a batch adsorption technique. Furthermore, FTIR and SEM examinations performed before and after the adsorption process indicated that the surface functionalisation and shape of the ZnO-NP nanocomposites changed significantly. A batch adsorption analysis was used to examine the extent to which operating parameters, the equilibrium isotherm, adsorption kinetics, and thermodynamics affected the results. The results of the batch technique revealed that the best results were obtained in the treatment with 0.04 g of ZnO-NP nanoparticles at 30 °C and pH 2 with an initial dye concentration of 10 mg L−1, which removed 91.5% and 65.6% of dye from synthetic and textile industry effluents, respectively. Additionally, six adsorption isotherm models were investigated by mathematical modelling and were validated for the adsorption process, and error function equations were applied to the isotherm model results in order to find the best-fit isotherm model. Likewise, the pseudo-second-order kinetic model fit well. A thermodynamic study revealed that IV2R adsorption on ZnO-NPs is a spontaneous, endothermic, and feasible sorption process. Finally, the synthesised nanocomposites prove to be excellent candidates for IV2R removal from water and real wastewater systems.


Introduction
In terms of industrial pollution, textile industries are a major source of pollution since they use a lot of water and chemicals and discharge high levels of toxic and nonbiodegradable dye effluents [1,2]. Dyes are widely applied in a variety of textile-based industries due to their advantageous properties [3,4], such as bright colours, water resistance, and ease of application [5]. Several sectors, including textile dyeing (60%), paper (10%), and plastic materials (10%), utilise large amounts of synthetic dyes (10%). According to some reports, there are approximately 100,000 commercially accessible dyes with a production capacity of more than 7 × 10 5 metric tons per year, with the textile industry using about 10,000 compounds [6]. Dye industry effluents, in particular, necessitate not only the treatment of problematic wastewater with high chemical and biological oxygen demands, suspended particles, and hazardous chemicals but also the treatment of dyes that are perceived by human eyes at very low concentrations [5]. When dyes are released In addition, the experimental equilibrium was applied using several different adsorption isotherm models to assess the adsorption mechanism, as well as thermodynamic and kinetic analyses. Furthermore, error functions were applied to reveal the most suitable model.

Materials
Polyvinyl pyrrolidone (PVP), H2O2, zinc nitrate, and sodium hydrosulphite were provided by M/s. Himedia Laboratories Pvt. Ltd., Mumbai, India. Table 1 shows the chemical structure of IV2R dye (molecular formula, molecular weight, and λ max). All of the used compounds were of analytical grade.

Synthesis of ZnO Nanoparticles
An improved approach was used to create ZnO-NPs, as described previously [26]. Briefly, 0.1 M Zn (NO3)2 was hydrolyzed with 250 mL of 0.2 M NaOH, and 1 % PVP was added and stirred continuously for 2 h using 100 mL of deionised water (Millipore, Milli-Q, Buenos Aires, Argentina). The pellet was produced by centrifuging the suspension at 3000 rpm for 5 min at 4 °C. At 75 °C, 1 M H2O2 was added and agitated for 1 h. The sample was also dried for 3 h at 65 °C in an oven before being maintained at 350 °C for 6 h.

Preparation of Dye Solution
The IV2R dye stock solution was made by weighing 1.00 g of powdered dye. The dye was transferred quantitatively into a 1 L measuring flask, which was then filled with distilled water to achieve a dye concentration of 1000 mg L −1 in the solution. The stock dye solutions (1000 mg L −1 ) were prepared separately and stored at 4 °C in distilled water. Dilution of the stock solutions with distilled water generated the working solutions.

ZnO-NP Characterisation
A scanning electron microscope was used to examine the morphologies of the adsorbents (JEOL JSM 6360 LA, Austin, USA). Brunauer-Emmett-Teller (BET) desorption-adsorption studies were also performed using a sorptiometer (Quantachrome TouchWinTM v1.2, Downers Grove, IL, USA) and N2 adsorption isotherms. Using N2 desorption-adsorption measurements in a N2 solution with a saturated vapour pressure of 33.5 atm and an adsorption temperature of 77.3 K, the average particle radius, mean pore diameter, and pore volume of the ZnO-NPs were estimated. Furthermore, a Shimadzu FTIR-8400 S was employed to analyse the functional groups of the pre-prepared ZnO-NP materials using

Synthesis of ZnO Nanoparticles
An improved approach was used to create ZnO-NPs, as described previously [26]. Briefly, 0.1 M Zn (NO 3 ) 2 was hydrolyzed with 250 mL of 0.2 M NaOH, and 1 % PVP was added and stirred continuously for 2 h using 100 mL of deionised water (Millipore, Milli-Q, Buenos Aires, Argentina). The pellet was produced by centrifuging the suspension at 3000 rpm for 5 min at 4 • C. At 75 • C, 1 M H 2 O 2 was added and agitated for 1 h. The sample was also dried for 3 h at 65 • C in an oven before being maintained at 350 • C for 6 h.

Preparation of Dye Solution
The IV2R dye stock solution was made by weighing 1.00 g of powdered dye. The dye was transferred quantitatively into a 1 L measuring flask, which was then filled with distilled water to achieve a dye concentration of 1000 mg L −1 in the solution. The stock dye solutions (1000 mg L −1 ) were prepared separately and stored at 4 • C in distilled water. Dilution of the stock solutions with distilled water generated the working solutions.

ZnO-NP Characterisation
A scanning electron microscope was used to examine the morphologies of the adsorbents (JEOL JSM 6360 LA, Austin, TX, USA). Brunauer-Emmett-Teller (BET) desorptionadsorption studies were also performed using a sorptiometer (Quantachrome TouchWinTM v1.2, Downers Grove, IL, USA) and N 2 adsorption isotherms. Using N 2 desorptionadsorption measurements in a N 2 solution with a saturated vapour pressure of 33.5 atm and an adsorption temperature of 77.3 K, the average particle radius, mean pore diameter, and pore volume of the ZnO-NPs were estimated. Furthermore, a Shimadzu FTIR-8400 S was employed to analyse the functional groups of the pre-prepared ZnO-NP materials using FTIR analysis spectrophotometry (Japan), and a UV-vis spectrophotometer was used to take spectrophotometric readings (UV 4000, MRI, Germany).

Adsorption Experiments 2.5.1. Batch Adsorption Experiments
In a 100 mL airtight Erlenmeyer flask, batch studies for IV2R elimination were carried out with 50 mL of 10 mg L −1 dye solution. The temperature of the system was kept constant at 25 • C. During the experiment, a weighed amount of adsorbent was placed in the flask and mixed with various concentrations of IV2R at constant and moderate mixing rates. At predetermined contact times, samples were removed and filtered to separate the adsorbent from the dye solution. The concentration of the dye solution was assessed using spectrometric techniques, and all experiments were carried out in triplicate, with the mean result reported. Furthermore, the effects of contact time (15 to 180 min) and adsorbent dosage (0.005 to 0.08 g) were investigated at a constant initial dye concentration of 10 mg L −1 . Experiments were conducted by altering the initial dye concentrations from 10 to 80 mg L −1 . The temperature (25 to 55 • C) was evaluated at a constant dose of 0.02 g of ZnO-NP adsorbent to determine the effects of the initial dye concentration on dye uptake. Moreover, the effects of final pH on IV2R adsorption were investigated by altering the dye's initial pH from 2 to 10. NaOH and HCl solutions were used to adjust the dye solutions' initial pH. A UV-visible spectrophotometer was used to measure the quantity of IV2R dye in the clear supernatant at 550 nm at any time after shaking.

Analytical Techniques
The dye concentrations were measured using a Spectronic Genesy 2PC UV-vis spectrophotometer at the wavelength of its maximum absorbance, λ max. Using the Beer-Lambert equation, the final dye concentration was measured spectrometrically to correspond to the dye's maximum concentration: where ε is the molar absorptivity, CS is the sample concentration, and l is the thickness of the absorbing material (1 cm). A pH/ion meter (WTW Inolab pH/ION Level 2, Germany) was used to determine the pH of the dye solution. Equations (2) and (3) were used to calculate the adsorption capacity (q e ) and dye elimination percentages [4]: where C i and C e (mg L −1 ) indicate the initial and final concentrations of IV2R at a given time, respectively, while V represents the volume of the dye mixture (L), and W represents the weight of the dry adsorbent (g).

Experiment with Isotherms
Experiments were carried out to develop adsorption isotherms for the ZnO-NP adsorbent by adjusting the adsorbent dose to 0.02 g at an initial dye concentration of 10 to 80 mg L −1 and an ambient pH of 2 for 3 h at 150 rpm, which was mixed with 50 mL of dye solution [36,37]. The following isotherm models were used to determine the most suitable one and calculate the data.

Freundlich Model
The ability of the Freundlich model to fit the experimental data was used to calculate the slope of n and the intercept value of K f by displaying a curve of log q e with respect to log C e . By visualizing the Freundlich model in logarithmic form [38], it is simple to linearise it: Log q e = Log K f + 1/n Log C e (4) The isotherm was used to determine the Freundlich constants K f and n.

The Langmuir Model
The Langmuir model is represented by the following mathematical formula [39]: q e = q max bC e /(1 + bC e ) where q max (mg g −1 ) is the maximum sorption capacity corresponding to the saturation capacity, and b (L mg −1 ) is a coefficient relating to the affinity between ZnO-NP and IV2R dye ions. The linear relationship obtained by graphing the curve (1/q e ) vs.
(1/C e ) is given in Equation (6). 1/q e = 1/(bq max C e ) + 1/q max (6) The slope and intercept, respectively, are used to calculate b and q max .

Henderson and Halsey Models
These models perform well with heteroporous substances and multilayer sorption. Using the equation below, the Halsey model [40] was calculated, as given in Equation (7).
where n and K are Halsey constants. Meanwhile, the Henderson model was obtained from the following equation: where the Henderson constants are nh and Kh.

The Harkins-Jura Model
This model describes multilayer adsorption, as well as the presence of heterogeneous pore scattering in an adsorbent [41]. This model was obtained from Equation (9).
where the isotherm constants are A and B.

Error Function Test
Various error functions were investigated in order to find the best and most appropriate model for investigating the equilibrium data. The error function tests were employed as follows: 2.7.1. Hybrid Fractional Error (HYBRID) Because it solves for low concentrations by balancing absolute deviation against fractional error and is more dependable than other error functions, the hybrid fractional error function was used. Equation (10) gives the hybrid error: where N is the total number of data points, and q e,exp and q e,cal (mg g −1 ) are the respective experimental and calculated adsorption capacity values. In addition, P is the number of isotherm factors.

Error Percentage Average (APE)
The APE model exhibits the suitability or trend between the predicted and observed values of the sorption capacity used to create model curves (APE) and can be designed using Equation (11).
where N is the number of data points under investigation.
The nonlinear chi-square test is a statistical method for determining which treatment system is most suitable. The following Equation (12) defines the approach to determining the chi-square error:

Sum of Absolute Errors (EABS)
An increase in errors improves the fit, resulting in a bias toward high-concentration data. The following Equation (13) can be used to evaluate EABS tests:

Adsorption Kinetics
Adsorption kinetic experiments were carried out using 0.05 g of adsorbent mixed separately with 50 mL of IV2R solution containing 10 mg L −1 concentrations in 100 mL conical flasks at a solution pH of 2, and the mixture was agitated at room temperature for 15, 30, 60, 120, and 180 min. The clear solutions were examined for any remaining IV2R concentrations.

Pseudo-First-Order Kinetics
The following equation gives the linear version of the generalised pseudo-first-order kinetics of dyes adsorbed at time t (mg g −1 ).
The linear formula of the pseudo-first-order kinetics is expressed in Equation (15).
Log (q e − q t ) = log (q e ) − k 1 t/2.303 (15) where q e is the amount of dye adsorbed at equilibrium (mg g −1 ), q t denotes the amount of time t, and K 1 denotes the pseudo-first-order rate constant (min −1 ).

Pseudo-Second-Order Kinetics
The pseudo-second-order Equation (16) is as follows: T/qt = 1/K 2 q e 2 + t/q e (16) where K 2 denotes the second-order rate constant (g mg −1 min). Plotting (t/q t ) versus (t) yields a linear relationship, and the slope and intercept can be used to derive the q e and K 2 parameters, respectively.

The Intraparticle Diffusion Model
The following is the intraparticle diffusion equation: K dif t 1/2 + C = q t (17) where C is the intercept, and K dif (mg g −1 min 0.5 ) is the intraparticle diffusion rate constant, which is calculated from the regression line's slope.

Thermodynamics of Adsorption
The value of Gibbs free energy change (∆G • ) is a fundamental principle of nonspontaneity, and a thermodynamic analysis is required to determine whether the nature of the sorption process is spontaneous or not. The following nonlinear forms (18)- (21) can be used to compute the Gibbs free energy change (∆G • kJ mol −1 ), enthalpy change (∆H • kJ mol −1 ), and entropy change (∆S • , J mol −1 K −1 ) parameters for the sorption process at various temperatures (e.g., 25, 30, 45, and 55 • C): In addition to ∆G • , which can be obtained from Equation (20) or (21), the values of ∆H • and ∆S • were computed using the intercept and slope of the plotted curve of T vs. ∆G • from Equation (20) or (21). The produced ZnO-NPs were identified using FTIR measurements, as depicted in Figure 1. Figure 1A shows the existence of IR signals before adsorption in the range of 3600-3000 cm −1 in dye samples, which indicate the presence of -OH and -NH 2 groups. Similarly, CH 2 stretching vibration is indicated by the peak at 2935 cm −1 . The stretching vibration of the C-H bond in -CH 2 groups is attributed to the band at 2860 cm −1 . According to the data, the peak at 1546 cm −1 is O-H bending vibration. Correspondingly, the peaks around 1394 and 1508 cm −1 indicate stretching of the C-H bond, which is attributed to CH 3 vibration. Likewise, the absorption bands in the range of 1000-1300 cm −1 indicate the presence of the C-O group. In addition, the spectra of ZnO-NPs reveal a strong peak at 426 cm −1 , which corresponds to the zinc-oxygen stretching mode [42]. The peaks at 949, 834, and 700 cm −1 generally correlate to C-H bending vibration.

Results and Discussion
On the other hand, Figure 1B shows that the bands at 3737, 3294, 2924, 2860, 1546, 1508, 1394, 1044, and 949 cm −1 were shifted to 3343, 2929, 2637, 1545, 1408, 1037, and 948 cm −1 after IV2R adsorption. According to this hypothesis, OH, C-H, C=C, and C=O groups may play a role in the adsorption of IV2R onto ZnO-NPs, while the C=O bond may be seen in the peaks between 1400 and 1500 cm −1 [43]. Additionally, the new peak at 1633 cm −1 is compatible with the C=O stretching of proteins [44]. C-H bending, C-O, or C-C stretching vibrations are represented by absorption bands that are located between 1100 and 1000 cm −1 . The presence of an aromatic ring on the dye compound is indicated by the band seen at 885 cm −1 , which is caused by aromatic C-H out-of-plane vibration [45]. absorption band with a centre radius of 667 cm −1 . Finally, the stretching vibrations of Zn-O are responsible for the peak range of 400-600 cm −1 [46].
It is well known that the production of ZnO particles from the hydrolysis of Zn 2+ ions in aqueous medium is a complicated process. Many polyvalent cationic species can form when Zn 2+ ions interact with OH ions, and their formation is significantly influenced by the pH of the solution. However, based on the pH, temperature, and synthesis processes, the precipitation of ZnO particles is usually defined by a development unit that can be either Zn (OH) 2 or Zn(OH) 2 4 ions. The dissolution-reprecipitation mechanism has been proposed as a mechanism for the formation of ZnO from Zn(OH) 2 [47]. ZnO will be produced by the chemical reaction [48]: stretching vibrations are represented by absorption bands that are located between 1100 and 1000 cm −1 . The presence of an aromatic ring on the dye compound is indicated by the band seen at 885 cm −1 , which is caused by aromatic C-H out-of-plane vibration [45]. Free O-H and N-H stretching vibrations are responsible for the robust and wide absorption band located at 3371 cm −1 . While the observed band at 1408 cm −h is related to the C-N stretching bond of amino acids, the band observed at a wavelength of 1037 cm −h can be ascribed to C-O-C. Moreover, C-H bending vibration is responsible for the weak absorption band with a centre radius of 667 cm −1 . Finally, the stretching vibrations of Zn-O are responsible for the peak range of 400-600 cm −1 [46].
It is well known that the production of ZnO particles from the hydrolysis of Zn 2+ ions in aqueous medium is a complicated process. Many polyvalent cationic species can form when Zn 2+ ions interact with OH ions, and their formation is significantly influenced by the pH of the solution. However, based on the pH, temperature, and synthesis processes, the precipitation of ZnO particles is usually defined by a development unit that can be either Zn (OH)2 or Zn(OH) 2 4 ions. The dissolution-reprecipitation mechanism has been proposed as a mechanism for the formation of ZnO from Zn(OH)2 [47]. ZnO will be produced by the chemical reaction [48]: ZnO + H2O.

BET Surface Analysis
The pore diameter, pore size distribution, and specific surface area of the ZnO-NPs were measured using BET and the nitrogen adsorption-desorption isotherm technique, as shown in Figure 2. It was revealed that an isotherm existed when the BET analysis was conducted at high P/P0 levels. Table 2 shows the BET results. The total pore volume and specific surface area of the ZnO-NPs were 95.83 m 2 g −1 and 0.058 cm 3 g −1 , respectively. The specific surface area calculated using the Langmuir technique was 140.692 m 2 g −1 . This superior property could provide ZnO-NPs with a larger surface area, allowing more active sorption sites to occur. The micropore volume of ZnO-NPs was determined to be 0.015 cm 3 g −1 using a cumulative BJH adsorption experiment. Furthermore, the average particle radius was 1.42 nm.
The pore diameter, pore size distribution, and specific surface area of the ZnO-NPs were measured using BET and the nitrogen adsorption-desorption isotherm technique, as shown in Figure 2. It was revealed that an isotherm existed when the BET analysis was conducted at high P/P0 levels. Table 2 shows the BET results. The total pore volume and specific surface area of the ZnO-NPs were 95.83 m 2 g −1 and 0.058 cm 3 g −1 , respectively. The specific surface area calculated using the Langmuir technique was 140.692 m 2 g −1 . This superior property could provide ZnO-NPs with a larger surface area, allowing more active sorption sites to occur. The micropore volume of ZnO-NPs was determined to be 0.015 cm 3 g −1 using a cumulative BJH adsorption experiment. Furthermore, the average particle radius was 1.42 nm.   Figure 3, which shows the particle morphologies. The nanostructure, which resembles nanoflowers and has   Figure 3, which shows the particle morphologies. The nanostructure, which resembles nanoflowers and has agglomerated nanoparticles with an average pore size of approximately 1.22 nm, refers to ZnO-NP powders that result in formations that resemble flowers in uniformly sized nanoporous channels. The length and thickness of these formations are in nanometres. SEM micrographs of pure ZnO-NPs reveal that they have a porous character and a large surface area. These flowers and swollen structures are preferable for the absorption of dye contaminants.
agglomerated nanoparticles with an average pore size of approximately 1.22 nm, refers to ZnO-NP powders that result in formations that resemble flowers in uniformly sized nanoporous channels. The length and thickness of these formations are in nanometres. SEM micrographs of pure ZnO-NPs reveal that they have a porous character and a large surface area. These flowers and swollen structures are preferable for the absorption of dye contaminants.

UV-Visible Spectra
The wavelength range of metal nanoparticles was determined by UV, and the results are shown in Table 3. The maximum absorbance is shown by the highest peak. The production of zinc oxide nanoparticles is confirmed by UV spectrometer absorption peaks in the 300-400 nm (390 nm) range. These peaks also indicate that the particles are nanosized and have a narrow particle size distribution. The absorption peak at 243 nm, which is associated with π→π * transitions in the sesquiterpene system, can hardly be seen in the UV spectrum of ZnO-NPs, which could be owing to a change in the sesquiterpene structure and the absence of these π→π * transitions [49]

Influence of pH
The shape and chemistry of the target dye ions, as well as the binding sites on the adsorbent, can be affected by the pH of a solution. In addition, the speciation of IV2R dye affects leads to an alteration in the reaction kinetics in addition to the equilibrium features of the sorption process [50]. Experiments were investigated by using different initial pH values of 2, 4, 6, 8, and 10 for the adsorbent to explore the effect of pH on the adsorption of IV2R. The percentage of adsorption tended to increase as the pH increased from 2 to 6. Furthermore, when the pH of the original dye solution was 2 and 6, the percentage of dye removal was at its maximum (90.7% and 89.17%, respectively) for the ZnO-NP adsorbent ( Figure 4). The maximum adsorption capacity was observed at pH 2 and 6 with 4.31 and 4.24 mg g −1 , respectively. These results indicate that solutions that are acidic in nature are

UV-Visible Spectra
The wavelength range of metal nanoparticles was determined by UV, and the results are shown in Table 3. The maximum absorbance is shown by the highest peak. The production of zinc oxide nanoparticles is confirmed by UV spectrometer absorption peaks in the 300-400 nm (390 nm) range. These peaks also indicate that the particles are nanosized and have a narrow particle size distribution. The absorption peak at 243 nm, which is associated with π→π * transitions in the sesquiterpene system, can hardly be seen in the UV spectrum of ZnO-NPs, which could be owing to a change in the sesquiterpene structure and the absence of these π→π * transitions [49]. The shape and chemistry of the target dye ions, as well as the binding sites on the adsorbent, can be affected by the pH of a solution. In addition, the speciation of IV2R dye affects leads to an alteration in the reaction kinetics in addition to the equilibrium features of the sorption process [50]. Experiments were investigated by using different initial pH values of 2, 4, 6, 8, and 10 for the adsorbent to explore the effect of pH on the adsorption of IV2R. The percentage of adsorption tended to increase as the pH increased from 2 to 6. Furthermore, when the pH of the original dye solution was 2 and 6, the percentage of dye removal was at its maximum (90.7% and 89.17%, respectively) for the ZnO-NP adsorbent ( Figure 4). The maximum adsorption capacity was observed at pH 2 and 6 with 4.31 and 4.24 mg g −1 , respectively. These results indicate that solutions that are acidic in nature are efficient for the adsorption of IV2R dye when the adsorption is below pH 6. In addition, these results prove that IV2R removal was slightly decreased at pH 4 with a percentage removal of 86.28% and an adsorption capacity of 4.10 mg g −1 . This result implies that there is a considerable electrostatic attraction between the adsorbent surface and the dye; this is likely owing to an increase in positively charged sites on the adsorbent surface when pH decreases, as reported by Netpradist et al. [51]. efficient for the adsorption of IV2R dye when the adsorption is below pH 6. In addition, these results prove that IV2R removal was slightly decreased at pH 4 with a percentage removal of 86.28% and an adsorption capacity of 4.10 mg g −1 . This result implies that there is a considerable electrostatic attraction between the adsorbent surface and the dye; this is likely owing to an increase in positively charged sites on the adsorbent surface when pH decreases, as reported by Netpradist et al. [51]. This study also revealed that anionic dye IV2R sorption is best in an acidic solution (positively charged adsorbent surface). Increased solution pH caused deprotonation and, as a result, a negatively charged surface of the adsorbent; consequently, this event may cause IV2R adsorption to decrease. Moreover, the nature of the surface charge of the adsorbent is affected by the pH of the solution in any adsorption system. In an acidic solution, the adsorbent's oxide surface acquires a net positive charge. Therefore, anionic dyes have a stronger electrostatic attraction in acidic solutions than in basic media [52]. Other probable reasons for such observations include dye-adsorbent interactions owing to hydrogen bonding and hydrophobic-hydrophobic interaction mechanisms. Furthermore, the surface area and pore size, which remain unaffected by pH changes, play a significant role during the process [53,54].
Ghoneim et al. [55] found that at a higher pH, elimination is reduced in comparison to the maximum condition. This can be explained by the binding site's ability to activate under normal conditions.

Influence of ZnO-NP Dose
The adsorbent dose is one of the most important factors for examining the impacts of the adsorption process to achieve the maximum adsorption capacity of the adsorbent by measuring the amount of ZnO-NP adsorbent. Experiments were performed with different adsorbent doses from 0.005 g to 0.08 g at a constant initial dye concentration of 10 mg L −1 to investigate the impact of the ZnO-NP dose on the adsorption system. At their respective equilibrium contact times, the percent adsorption increased from roughly 83.5 to 91.4 % as the adsorbent dose increased from 0.005 g to 0.04 g ( Figure 5). It is self-evident that as the adsorbent dose increases, the dye uptake is enhanced. It is well established that as the adsorbent dosage increases, the % adsorption increases. The adsorbed amount per unit mass decreases. The amount adsorbed per unit mass was only 3.4 mg g −1 adsorbent when the dose was 0.08 g. The removal percentage was roughly 89%. The adsorption density decreases as the adsorbent dose increases due to unsaturated adsorption sites during the adsorption process [56]. Another cause could be intraparticle interactions, such as aggregation, as a result of a high adsorbent dose. The total surface area of the adsorbent is reduced because of this aggregation, while the diffusion path length rises [57].  This study also revealed that anionic dye IV2R sorption is best in an acidic solution (positively charged adsorbent surface). Increased solution pH caused deprotonation and, as a result, a negatively charged surface of the adsorbent; consequently, this event may cause IV2R adsorption to decrease. Moreover, the nature of the surface charge of the adsorbent is affected by the pH of the solution in any adsorption system. In an acidic solution, the adsorbent's oxide surface acquires a net positive charge. Therefore, anionic dyes have a stronger electrostatic attraction in acidic solutions than in basic media [52]. Other probable reasons for such observations include dye-adsorbent interactions owing to hydrogen bonding and hydrophobic-hydrophobic interaction mechanisms. Furthermore, the surface area and pore size, which remain unaffected by pH changes, play a significant role during the process [53,54].
Ghoneim et al. [55] found that at a higher pH, elimination is reduced in comparison to the maximum condition. This can be explained by the binding site's ability to activate under normal conditions.

Influence of ZnO-NP Dose
The adsorbent dose is one of the most important factors for examining the impacts of the adsorption process to achieve the maximum adsorption capacity of the adsorbent by measuring the amount of ZnO-NP adsorbent. Experiments were performed with different adsorbent doses from 0.005 g to 0.08 g at a constant initial dye concentration of 10 mg L −1 to investigate the impact of the ZnO-NP dose on the adsorption system. At their respective equilibrium contact times, the percent adsorption increased from roughly 83.5 to 91.4 % as the adsorbent dose increased from 0.005 g to 0.04 g ( Figure 5). It is self-evident that as the adsorbent dose increases, the dye uptake is enhanced. It is well established that as the adsorbent dosage increases, the % adsorption increases. The adsorbed amount per unit mass decreases. The amount adsorbed per unit mass was only 3.4 mg g −1 adsorbent when the dose was 0.08 g. The removal percentage was roughly 89%. The adsorption density decreases as the adsorbent dose increases due to unsaturated adsorption sites during the adsorption process [56]. Another cause could be intraparticle interactions, such as aggregation, as a result of a high adsorbent dose. The total surface area of the adsorbent is reduced because of this aggregation, while the diffusion path length rises [57].

Influence of Contact Time
The equilibrium time is essential when considering economical water and wastewater applications, and contact time is an important component in all transfer phenomena for the adsorption process. The adsorption process was studied to establish the best contact time between 15 and 180 min. Figure 6 shows the adsorption removal of IV2R dye by the ZnO-NP adsorbent. Equilibrium was reached, and the optimum contact time for IV2R dye was chosen to be 120 min for the sorbate-sorbent contact. The uptake of IV2R was observed to occur in two phases as a function of time. The first phase involved fast dye uptake during the first 10 min of sorbate-sorbent interactions, followed by a slow dye removal phase that lasted significantly longer (>120 min) until equilibrium was reached. The higher sorption value at the start of the process could be attributable to the abundance of active sites on the sorbent at this time. The sorption process becomes less efficient during the slower phase as these sites are gradually occupied [58]. The IV2R removal effectiveness by ZnO-NPs was found to be 88.9% (4.23 mg g −1 ) at 180 min. Two-stage sorption has been extensively documented in the literature, with the first being quick and quantitatively dominant and the second being slower and quantitatively insignificant [59]. Furthermore, Ananta et al. [60] reported that adsorption was rapid at the first contact time, which happens in the early stage of adsorption after a few minutes, since most of

Influence of Contact Time
The equilibrium time is essential when considering economical water and wastewater applications, and contact time is an important component in all transfer phenomena for the adsorption process. The adsorption process was studied to establish the best contact time between 15 and 180 min. Figure 6 shows the adsorption removal of IV2R dye by the ZnO-NP adsorbent. Equilibrium was reached, and the optimum contact time for IV2R dye was chosen to be 120 min for the sorbate-sorbent contact. The uptake of IV2R was observed to occur in two phases as a function of time. The first phase involved fast dye uptake during the first 10 min of sorbate-sorbent interactions, followed by a slow dye removal phase that lasted significantly longer (>120 min) until equilibrium was reached. The higher sorption value at the start of the process could be attributable to the abundance of active sites on the sorbent at this time. The sorption process becomes less efficient during the slower phase as these sites are gradually occupied [58]. The IV2R removal effectiveness by ZnO-NPs was found to be 88.9% (4.23 mg g −1 ) at 180 min. Two-stage sorption has been extensively documented in the literature, with the first being quick and quantitatively dominant and the second being slower and quantitatively insignificant [59].

Influence of Contact Time
The equilibrium time is essential when considering economical water and wastewater applications, and contact time is an important component in all transfer phenomena for the adsorption process. The adsorption process was studied to establish the best contact time between 15 and 180 min. Figure 6 shows the adsorption removal of IV2R dye by the ZnO-NP adsorbent. Equilibrium was reached, and the optimum contact time for IV2R dye was chosen to be 120 min for the sorbate-sorbent contact. The uptake of IV2R was observed to occur in two phases as a function of time. The first phase involved fast dye uptake during the first 10 min of sorbate-sorbent interactions, followed by a slow dye removal phase that lasted significantly longer (>120 min) until equilibrium was reached. The higher sorption value at the start of the process could be attributable to the abundance of active sites on the sorbent at this time. The sorption process becomes less efficient during the slower phase as these sites are gradually occupied [58]. The IV2R removal effectiveness by ZnO-NPs was found to be 88.9% (4.23 mg g −1 ) at 180 min. Two-stage sorption has been extensively documented in the literature, with the first being quick and quantitatively dominant and the second being slower and quantitatively insignificant [59]. Furthermore, Ananta et al. [60] reported that adsorption was rapid at the first contact time, which happens in the early stage of adsorption after a few minutes, since most of Furthermore, Ananta et al. [60] reported that adsorption was rapid at the first contact time, which happens in the early stage of adsorption after a few minutes, since most of the binding sites are free and the adsorbent sites are empty, allowing the dye ions to bind quickly to the adsorbent. David and Joseph [61] mentioned that adsorption occurs quickly and is often controlled by diffusion from the majority of the solution on the surface. Moreover, Kumar and Gayathri [62] stated that adsorption increases as the contact time increases, which is possible due to a larger surface area of ZnO-NPs being available at the start and the exhaustion of the conversion of external adsorption sites, where the adsorbate (dye particles) is transported from the external to the internal sites of ZnO-NP adsorbent molecules.

Influence of the Initial Dye Concentration
The amount of dye removed from an aqueous mixture is highly dependent on the dye concentration. At a constant temperature, the adsorption process for IV2R was examined at concentrations ranging from 10 to 80 mg L −1 . The effect of IV2R dye concentrations on adsorption is seen in Figure 7. As can be seen, raising the initial dye concentration from 10 to 80 mg L −1 improved the ZnO-NP adsorption capacity from 4.33 to 20.58 mg g −1 . This could be owing to the strong driving force that occurs when the initial concentration of the adsorbates is increased, overcoming the mass transfer resistance between the aqueous and solid phases [20]. Furthermore, the increased dye clearance at higher concentrations is most likely due to greater diffusion and decreased dye absorption resistance [63]. The availability of active sites on the adsorbent and the final occupancy of these sites are attributed to the rapid first stage of dye removal with this adsorbent; the sorption then becomes less efficient. In addition, the greatest uptake of IV2R dye was found at a concentration of 10 mg L −1 (91.14%). The equilibrium loading capacity and initial dye concentration have a strong linear relationship with correlation coefficients greater than 0.99 for the ZnO-NP adsorbent. A dye concentration of 10 mg L −1 was used for further experiments. In general, the decreasing percentage of dye removal with increasing dye concentration could be due to the increase in sorption sites on the adsorbent surface [64]. Furthermore, the high probability of dye ions colliding with the adsorbent surface and the high rate of dye ion diffusion onto the adsorbent surface could be linked to the large amount of dye adsorbed at a high dye concentration [65]. This finding could point to the potential of treating textile effluent with a higher dye concentration. the binding sites are free and the adsorbent sites are empty, allowing the dye ions to bind quickly to the adsorbent. David and Joseph [61] mentioned that adsorption occurs quickly and is often controlled by diffusion from the majority of the solution on the surface. Moreover, Kumar and Gayathri [62] stated that adsorption increases as the contact time increases, which is possible due to a larger surface area of ZnO-NPs being available at the start and the exhaustion of the conversion of external adsorption sites, where the adsorbate (dye particles) is transported from the external to the internal sites of ZnO-NP adsorbent molecules.

Influence of the Initial Dye Concentration
The amount of dye removed from an aqueous mixture is highly dependent on the dye concentration. At a constant temperature, the adsorption process for IV2R was examined at concentrations ranging from 10 to 80 mg L −1 . The effect of IV2R dye concentrations on adsorption is seen in Figure 7. As can be seen, raising the initial dye concentration from 10 to 80 mg L −1 improved the ZnO-NP adsorption capacity from 4.33 to 20.58 mg g −1 . This could be owing to the strong driving force that occurs when the initial concentration of the adsorbates is increased, overcoming the mass transfer resistance between the aqueous and solid phases [20]. Furthermore, the increased dye clearance at higher concentrations is most likely due to greater diffusion and decreased dye absorption resistance [63]. The availability of active sites on the adsorbent and the final occupancy of these sites are attributed to the rapid first stage of dye removal with this adsorbent; the sorption then becomes less efficient. In addition, the greatest uptake of IV2R dye was found at a concentration of 10 mg L −1 (91.14%). The equilibrium loading capacity and initial dye concentration have a strong linear relationship with correlation coefficients greater than 0.99 for the ZnO-NP adsorbent. A dye concentration of 10 mg L −1 was used for further experiments. In general, the decreasing percentage of dye removal with increasing dye concentration could be due to the increase in sorption sites on the adsorbent surface [64]. Furthermore, the high probability of dye ions colliding with the adsorbent surface and the high rate of dye ion diffusion onto the adsorbent surface could be linked to the large amount of dye adsorbed at a high dye concentration [65]. This finding could point to the potential of treating textile effluent with a higher dye concentration. 3.2.5. Influence of Temperature Figure 8 illustrates the percentage removal of IV2R dye ions by ZnO-NPs. The maximum percentage of dye removal was obtained at 45 • C with 94.9%, which rose when the temperature was increased from 25 to 45 • C for the IV2R dye on ZnO-NPs. Because of the decreased solution viscosity as the temperature rose from 25 to 45 • C, the diffusion rate of the adsorbate molecules within the pores changed, as did the equilibrium capacity of the ZnO-NPs for a specific adsorbate. Temperature increases (above 45 • C) resulted in a decrease in the percentage removed. This is attributed to a reduction in surface activity [66]. The decrease in adsorption efficiency can be attributed to a variety of factors, including deactivation of the adsorbent surface, a growing tendency for dyes to migrate from the solid to the bulk stage, and the destruction of specific active sites on the adsorbent surface due to bond ruptures. According to Sivaprakash et al. [67], this could be due to an increase in the mobility of the large dye ion as temperature rises. A growing number of molecules may be able to obtain enough energy to interact with active areas on the surface. Furthermore, rising temperatures may cause a swelling effect within the adsorbent's internal structure, allowing large dyes to penetrate further. Additionally, there are two possible explanations for this outcome. At high temperatures, the pore diameters of adsorbent particles will expand [68]. Second, due to the breakage of some internal bonds, such as hydrogen bonds between the dye ion and the hydroxyl groups on the surface of the adsorbent's active surface sites, the number of adsorption sites will increase [69,70]. temperature was increased from 25 to 45 °C for the IV2R dye on ZnO-NPs. Because of the decreased solution viscosity as the temperature rose from 25 to 45 °C, the diffusion rate of the adsorbate molecules within the pores changed, as did the equilibrium capacity of the ZnO-NPs for a specific adsorbate. Temperature increases (above 45 °C) resulted in a decrease in the percentage removed. This is attributed to a reduction in surface activity [66]. The decrease in adsorption efficiency can be attributed to a variety of factors, including deactivation of the adsorbent surface, a growing tendency for dyes to migrate from the solid to the bulk stage, and the destruction of specific active sites on the adsorbent surface due to bond ruptures. According to Sivaprakash et al. [67], this could be due to an increase in the mobility of the large dye ion as temperature rises. A growing number of molecules may be able to obtain enough energy to interact with active areas on the surface. Furthermore, rising temperatures may cause a swelling effect within the adsorbent's internal structure, allowing large dyes to penetrate further. Additionally, there are two possible explanations for this outcome. At high temperatures, the pore diameters of adsorbent particles will expand [68]. Second, due to the breakage of some internal bonds, such as hydrogen bonds between the dye ion and the hydroxyl groups on the surface of the adsorbent's active surface sites, the number of adsorption sites will increase [69,70].

Isothermal Analysis
The equilibrium isotherm is a plot of the quantity of sorbate extracted per unit sorbent (qe) as the sorbent's solid phase concentration against the sorbate's liquid phase concentration (Ce). For the design and optimisation of an adsorption system for the removal of a dye from an aqueous solution, equilibrium isotherm data are necessary [71]. As a consequence, the most appropriate correlation for the equilibrium curve must be determined. To test the validity of the experimental data, a number of isotherm models were applied [22]. Therefore, in the present study, the most commonly used models, namely, the Langmuir, Harkins-Jura, Freundlich, Halsey, Henderson, and Tempkin isotherms, were used to describe the adsorption equilibrium.

Freundlich Isotherm
On the basis of the assumption of energy surface heterogeneity, the Freundlich isotherm model is the earliest known relationship describing non-ideal and reversible adsorption, which can be extended to multilayer adsorption. The obtained results were fit to

Isothermal Analysis
The equilibrium isotherm is a plot of the quantity of sorbate extracted per unit sorbent (q e ) as the sorbent's solid phase concentration against the sorbate's liquid phase concentration (C e ). For the design and optimisation of an adsorption system for the removal of a dye from an aqueous solution, equilibrium isotherm data are necessary [71]. As a consequence, the most appropriate correlation for the equilibrium curve must be determined. To test the validity of the experimental data, a number of isotherm models were applied [22]. Therefore, in the present study, the most commonly used models, namely, the Langmuir, Harkins-Jura, Freundlich, Halsey, Henderson, and Tempkin isotherms, were used to describe the adsorption equilibrium.

Freundlich Isotherm
On the basis of the assumption of energy surface heterogeneity, the Freundlich isotherm model is the earliest known relationship describing non-ideal and reversible adsorption, which can be extended to multilayer adsorption. The obtained results were fit to the Freundlich isotherm model's experimental data, which was supported by a strong correlation coefficient of R 2 = 0.994 for the ZnO-NPs, showing that this model is beneficial for the adsorption process, as shown in Figure 9. The IV2R dye and the adsorbents formed a strong bond, as indicated by the value of 1/n, also known as the heterogeneity factor, which describes the divergence from sorption linearity as follows: When 1/n equals 1, the adsorption is linear, and the concentration of dye particles has no effect on the two stages. When 1/n is less than 1, chemical adsorption occurs; when 1/n is greater than 1, cooperative adsorption occurs, which is more physically advantageous and contains strong contacts among the adsorbate particles [72]. The value of the factor "1/n" was smaller than 1 in this study; the results indicate that a chemical sorption method on an external surface is preferable with this isotherm equation [22]. the Freundlich isotherm model's experimental data, which was supported by a strong correlation coefficient of R 2 = 0.994 for the ZnO-NPs, showing that this model is beneficial for the adsorption process, as shown in Figure 9. The IV2R dye and the adsorbents formed a strong bond, as indicated by the value of 1/n, also known as the heterogeneity factor, which describes the divergence from sorption linearity as follows: When 1/n equals 1, the adsorption is linear, and the concentration of dye particles has no effect on the two stages. When 1/n is less than 1, chemical adsorption occurs; when 1/n is greater than 1, cooperative adsorption occurs, which is more physically advantageous and contains strong contacts among the adsorbate particles [72]. The value of the factor "1/n" was smaller than 1 in this study; the results indicate that a chemical sorption method on an external surface is preferable with this isotherm equation [22].

Langmuir Isotherm
The most popular model for quantifying the amount of adsorbate on an adsorbent as a function of partial pressure or concentration at a particular temperature is the Langmuir adsorption model. The Langmuir isotherm is based on the assumption of monolayer adsorption on a structurally homogeneous adsorbent, in which all adsorption sites are similar and energetically equivalent, adsorption occurs at specific homogeneous sites on the adsorbent, and once a dye molecule occupies a site, no further adsorption can occur. Table  4 shows the estimated parameters. The results obtained using the Langmuir isotherm coincided with data obtained throughout the experiment, with a strong correlation coefficient of R 2 = 0.974 for the ZnO-NPs. In addition, the ZnO-NP dye has a high maximum absorption capacity (qmax) of 119.05 mg g −1 (Figure 10). This is in line with the creation of a full monolayer on the adsorbent surface. The Langmuir constant (b), which is related to the heat of adsorption, was found to be 0.119. The dimensionless separation factor (RL), which is described below, can also be used to predict the affinity between the sorbate and the sorbent using Langmuir parameters (22): According to the parameters in Table 3, the value of RL can be used to determine whether a sorption system is "favourable" or "unfavourable." The sorption of IV2R onto ZnO-NPs has an RL value of 0.597, which shows that adsorption of IV2R on ZnO-NPs was "favourable".

Langmuir Isotherm
The most popular model for quantifying the amount of adsorbate on an adsorbent as a function of partial pressure or concentration at a particular temperature is the Langmuir adsorption model. The Langmuir isotherm is based on the assumption of monolayer adsorption on a structurally homogeneous adsorbent, in which all adsorption sites are similar and energetically equivalent, adsorption occurs at specific homogeneous sites on the adsorbent, and once a dye molecule occupies a site, no further adsorption can occur. Table 4 shows the estimated parameters. The results obtained using the Langmuir isotherm coincided with data obtained throughout the experiment, with a strong correlation coefficient of R 2 = 0.974 for the ZnO-NPs. In addition, the ZnO-NP dye has a high maximum absorption capacity (q max ) of 119.05 mg g −1 (Figure 10). This is in line with the creation of a full monolayer on the adsorbent surface. The Langmuir constant (b), which is related to the heat of adsorption, was found to be 0.119. The dimensionless separation factor (R L ), which is described below, can also be used to predict the affinity between the sorbate and the sorbent using Langmuir parameters (22):

Harkins-Jura Isotherm
Multilayer adsorption is accounted for by the Harkins-Jura model, which can be explained by the existence of a heterogeneous pore distribution [41]. Equation (8) can be used to solve the Harkins-Jura adsorption isotherm, which can be seen in Figure 11   According to the parameters in Table 3, the value of R L can be used to determine whether a sorption system is "favourable" or "unfavourable." The sorption of IV2R onto ZnO-NPs has an R L value of 0.597, which shows that adsorption of IV2R on ZnO-NPs was "favourable".

Harkins-Jura Isotherm
Multilayer adsorption is accounted for by the Harkins-Jura model, which can be explained by the existence of a heterogeneous pore distribution [41]. Equation (8) can be used to solve the Harkins-Jura adsorption isotherm, which can be seen in Figure 11 as a plot of 1/q e vs. log C e . The presence of a heterogeneous pore distribution and multilayer adsorption can be explained by the Harkins-Jura model. The isotherm constants are presented in Table 4, and the correlation coefficient was found to be R 2 = 0.905. This could mean that the Harkins-Jura model is useful for adsorption data.

Isotherm Models of Halsey and Henderson
The Halsey and Henderson isotherm models are suitable for multilayer adsorption, and the fitting of the Halsey equation can be applied to heteroporous solids [73]. Figures 12  and 13 show plots of Ln q e versus Ln C e for Halsey and ln [ln(1C e )] versus Ln q e for Henderson isotherms, respectively. Table 4 summarises the isotherm constants and correlation coefficients. The correlation coefficient for Halsey was R 2 = 0.994, while Henderson's was R 2 = 0.995. The results of Halsey and Henderson suggest that both models can be used to predict IV2R adsorption on ZnO-NPs.

Isotherm Models of Halsey and Henderson
The Halsey and Henderson isotherm models are suitable for multilayer adsorption, and the fitting of the Halsey equation can be applied to heteroporous solids [73]. Figures  12 and 13 show plots of Ln qe versus Ln Ce for Halsey and ln [ln(1Ce)] versus Ln qe for Henderson isotherms, respectively. Table 4 summarises the isotherm constants and correlation coefficients. The correlation coefficient for Halsey was R 2 = 0.994, while Henderson's was R 2 = 0.995. The results of Halsey and Henderson suggest that both models can be used to predict IV2R adsorption on ZnO-NPs.

Tempkin Isotherm
The Tempkin isotherm equation includes a component that describes the interactions between the adsorbing species and adsorbate [74]. It is assumed that due to adsorbateadsorbate repulsions, the heat of adsorption of all molecules in the layer decreases linearly with coverage, and that adsorption is a uniform distribution of maximum binding energy [75]. The Tempkin adsorption isotherm model was used to evaluate the adsorption potentials of the ZnO-NPs for IV2R. The derivation of the Tempkin isotherm assumes that the fall in the heat of adsorption is linear rather than logarithmic, as implied in the Freundlich equation. The data indicate that the Tempkin isotherm model applies to the adsorption of IV2R dye onto ZnO-NPs, as shown by the high linear regression correlation coefficient (R 2 = 0.928), as presented in Figure 14, which indicates that the Tempkin isotherm is appropriate for the equilibrium data attained for the adsorption of IV2R on ZnO-NPs.

Tempkin Isotherm
The Tempkin isotherm equation includes a component that describes the interactions between the adsorbing species and adsorbate [74]. It is assumed that due to adsorbateadsorbate repulsions, the heat of adsorption of all molecules in the layer decreases linearly with coverage, and that adsorption is a uniform distribution of maximum binding energy [75]. The Tempkin adsorption isotherm model was used to evaluate the adsorption potentials of the ZnO-NPs for IV2R. The derivation of the Tempkin isotherm assumes that the fall in the heat of adsorption is linear rather than logarithmic, as implied in the Freundlich equation. The data indicate that the Tempkin isotherm model applies to the adsorption of IV2R dye onto ZnO-NPs, as shown by the high linear regression correlation coefficient (R 2 = 0.928), as presented in Figure 14, which indicates that the Tempkin isotherm is appropriate for the equilibrium data attained for the adsorption of IV2R on ZnO-NPs. [75]. The Tempkin adsorption isotherm model was used to evaluate the adsorption potentials of the ZnO-NPs for IV2R. The derivation of the Tempkin isotherm assumes that the fall in the heat of adsorption is linear rather than logarithmic, as implied in the Freundlich equation. The data indicate that the Tempkin isotherm model applies to the adsorption of IV2R dye onto ZnO-NPs, as shown by the high linear regression correlation coefficient (R 2 = 0.928), as presented in Figure 14, which indicates that the Tempkin isotherm is appropriate for the equilibrium data attained for the adsorption of IV2R on ZnO-NPs.

Examining Error Functions to Find the Most Appropriate Isotherm Model
Several error functions, such as the hybrid fractional error, average percentage error (APE) equation, chi-square error (X 2 ) equation, and the sum of absolute errors, were used to determine the best-fit model for the investigational data (EABS). Table 4 summarises the data gathered from the various error functions. For each isotherm model, the analysed error functions produced varying findings, and the comparison across isotherm models should focus on each error function independently [2]. Tempkin > Henderson >

Examining Error Functions to Find the Most Appropriate Isotherm Model
Several error functions, such as the hybrid fractional error, average percentage error (APE) equation, chi-square error (X 2 ) equation, and the sum of absolute errors, were used to determine the best-fit model for the investigational data (EABS). Table 4 summarises the data gathered from the various error functions. For each isotherm model, the analysed error functions produced varying findings, and the comparison across isotherm models should focus on each error function independently [2]. Tempkin > Henderson > Freundlich > Langmuir > Halsey > Harkins-Jura are the most appropriate isotherm models based on the observed data. Nonetheless, the error function test offered variable data for all models, and the model evaluation focused on each error function individually.

Model of Pseudo-First-Order Kinetics
The pseudo-first-order equation explains the adsorption rate based on the adsorption capacity. According to this model [76], the ratio of occupied to empty adsorption sites is proportional to the number of vacant sites. At a dye concentration of 10 mg L −1 , the linear figure of log (q e -q t ) versus t is shown in Figure 15. The intercept and the slope of the linear plots for the removal of the IV2R dye from ZnO-NPs were used to compute the q e and K 1 values. The values of K 1 , the experimental and calculated values of q e , and the correlation coefficients for the pseudo-first-order kinetic plots are provided in Table 5. The obtained experimental data do not agree with the theoretical values of q e . This implies that the kinetic data and the pseudo-first-order model do not match well. Moreover, because the correlation coefficients (R 2 ) for the current experimental results for IV2R dye on ZnO-NPs were small, the pseudo-first-order equation was ruled out.

Model of Pseudo-First-Order Kinetics
The pseudo-first-order equation explains the adsorption rate based on the adsorption capacity. According to this model [76], the ratio of occupied to empty adsorption sites is proportional to the number of vacant sites. At a dye concentration of 10 mg L −1 , the linear figure of log (qe-qt) versus t is shown in Figure 15. The intercept and the slope of the linear plots for the removal of the IV2R dye from ZnO-NPs were used to compute the qe and K1 values. The values of K1, the experimental and calculated values of qe, and the correlation coefficients for the pseudo-first-order kinetic plots are provided in Table 5. The obtained experimental data do not agree with the theoretical values of qe. This implies that the kinetic data and the pseudo-first-order model do not match well. Moreover, because the correlation coefficients (R 2 ) for the current experimental results for IV2R dye on ZnO-NPs were small, the pseudo-first-order equation was ruled out.  Time (min) Figure 15. Adsorption kinetics of the pseudo-first-order kinetics of IV2R adsorption onto the ZnO-NPs. The rate-limiting step is assumed to be due to chemical adsorption containing valence forces through the exchange and/or sharing of electrons between dye ions and the adsorbent in the pseudo-second-order equation (adsorbent) [76]. As shown in Figure 16 and Table 5, plotting (t/q t ) against (t) should yield a linear correlation from which the data for parameters q e and k 2 may be determined from the slope and intercept, respectively. The second-order kinetic model's correlation coefficient was more than 0.999, pointing to the fact that the pseudo-second-order kinetic model resulted in a good correlation for IV2R adsorption onto ZnO-NPs. Apart from that, it was clear that the k 2 parameter value was higher than the corresponding k 1 parameter value. This is because, according to David and Joseph [61], the adsorption rate is proportional to the square of the number of empty sites in the pseudo-second-order model. for parameters qe and k2 may be determined from the slope and intercept, respectively. The second-order kinetic model's correlation coefficient was more than 0.999, pointing to the fact that the pseudo-second-order kinetic model resulted in a good correlation for IV2R adsorption onto ZnO-NPs. Apart from that, it was clear that the k2 parameter value was higher than the corresponding k1 parameter value. This is because, according to David and Joseph [61], the adsorption rate is proportional to the square of the number of empty sites in the pseudo-second-order model.

The Intraparticle Diffusion Equation
The adsorption technique requires a number of steps, including the transport of solute particles from the aqueous part to the exterior of the solid molecules, followed by the diffusion of solute molecules into the cavities' interior decoration, which is likely to be a slow process and a rate-determining step [77]. The figures of qt inverse t 0.5 may indicate a multilinear correlation, indicating that the adsorption process occurs in two or more stages ( Figure 17 and Table 5). The slope directly estimates the rate constant Kdif, and the intercept is C, as shown in Table 4. Because the barrier to exterior mass transfer increases as the intercept increases, the value of the C factor provides information about the thickness of the border layer.

The Intraparticle Diffusion Equation
The adsorption technique requires a number of steps, including the transport of solute particles from the aqueous part to the exterior of the solid molecules, followed by the diffusion of solute molecules into the cavities' interior decoration, which is likely to be a slow process and a rate-determining step [77]. The figures of q t inverse t 0.5 may indicate a multilinear correlation, indicating that the adsorption process occurs in two or more stages ( Figure 17 and Table 5). The slope directly estimates the rate constant K dif , and the intercept is C, as shown in Table 4. Because the barrier to exterior mass transfer increases as the intercept increases, the value of the C factor provides information about the thickness of the border layer. diffusion of solute molecules into the cavities' interior decoration, which is likely to be a slow process and a rate-determining step [77]. The figures of qt inverse t 0.5 may indicate a multilinear correlation, indicating that the adsorption process occurs in two or more stages ( Figure 17 and Table 5). The slope directly estimates the rate constant Kdif, and the intercept is C, as shown in Table 4. Because the barrier to exterior mass transfer increases as the intercept increases, the value of the C factor provides information about the thickness of the border layer.  Table 6 shows the thermodynamic parameters of IV2R adsorption onto ZnO-NPs. The large negative value of ∆G o demonstrates that dye sorption was spontaneous and feasible. Nonetheless, the free energy values in Table 6 increase with increasing temperature, demonstrating that the adsorption method is endothermic [78,79] and suggesting that the adsorption method's spontaneity decreases at lower temperatures. For physisorption, the range of ∆G o values is between −20 and 0 kJ mol −1 , whereas, for chemisorption, the range is between −80 and −400 kJ mol −1 [3]. The activation energy measurements in this work indicate that IV2R sorption onto ZnO-NPs occurs via a physisorption process.

Application to Real-Life Wastewater
To test the validity of using ZnO-NPs as adsorbents, real wastewater was mixed with simulated dye samples to see if the adsorbent could remove IV2R under optimum conditions. The results revealed that changing the kind of water had a significant impact on dye removal, with deionised water having the least impact on the adsorption process, with 91.51% of dye removed at an acidic pH after 180 min. On the other hand, real wastewater contains very high quantities of interfering ions from a variety of contaminants, which had a major impact on the IV2R dye removal effectiveness, with 62.6% of dye removed after 180 min. Regardless, the results show that ZnO-NP adsorbents may be employed successfully to remove IV2R dye from aqueous mixes and wastewater at a reasonable cost.

Comparative Studies of ZnO-NP Sorption Capacity
To demonstrate the efficacy of ZnO-NPs, the results of this experiment were compared to those of other studies on the adsorption capabilities of various dyes. This adsorption