Optimization of ZnO Nanoparticles’ Synthesis via Precipitation Method Applying Taguchi Robust Design

: Zinc oxide (ZnO) possesses exceptional potential to be utilized in water and wastewater treatment applications, either as a photocatalyst or in membrane incorporation. In the present study, ZnO nanoparticles were synthesized using the precipitation method. The Taguchi approach with the L 32b orthogonal array was utilized in order to optimize the experimental conditions for the synthesis of the nanoparticles and to ensure that relatively smaller-sized particles were obtained. The design was characterized by ten factors, where nine of them possessed four levels, while one had two levels. This study’s design factors were the type of Zn precursor, the concentration of the Zn precursor, the type of precipitating agent, the precipitation agent’s concentration, the type of utilized solvent, the pH value of the solvent, the temperature used during the synthetic procedure, the calcination temperature, the time of stirring during synthesis, as well as the stirring speed. The influences of those factors on the selected response parameters (the average crystallite size, degree of crystallinity, energy band gap (E g ), and photodegradation constant (k)) were then evaluated. XRD analysis and the calculated E g values indicated that the hexagonal wurtzite structure was the only crystalline phase present in the produced samples. The photocatalytic efficiency of all ZnO nano-particles was examined in the degradation of rhodamine B under UV light irradiation. The optimal conditions were achieved using zinc acetate dihydrate as the Zn precursor at a concentration equal to 0.3 M, sodium hydroxide as the precipitating agent (1.5 M), methanol as the solvent (the pH value of the solvent was equal to 13), a temperature during the synthetic procedure of 70 °C, 600 °C as calcination temperature, a 90 min stirring time, and 700 rpm as the stirring speed. The optimized ZnO sample was synthesized based on the aforementioned conditions and thoroughly characterized. The acquired results confirmed the prediction of the Taguchi approach, and the most enhanced k-value was observed.


Introduction
According to the findings presented in the World Water Development Report (2020) [1], the safeguarding of water resources in the context of global climate change has become a pivotal concern in today's prominent worldwide challenges.Currently, natural aquatic ecosystems are vulnerable to an array of hazardous chemical and biological substances stemming from various origins such as civil society, the public sector, and notably, the industrial sphere.
Multiple sectors, including but not limited to textiles, pharmaceuticals, paper production, printing, cosmetics, as well as food manufacturing, produce a significant volume of wastewater that contains a variety of organic pollutants, including dyes, resulting in contamination.In relation to the textile industry, it is approximated that approximately 200,000 tons of dye are released into aquatic bodies annually during dyeing and finishing processes, primarily due to inefficiencies in the dye application procedure [2].Even at concentrations lower than 1 mg/L, these dyes contribute color to industrial wastewater, disrupt the penetration of sunlight into water bodies, and disrupt the balance of eutrophication processes [3].
Rhodamine B, which is being widely employed, poses a significant threat to aquatic ecosystems and human health.This particular dye is recognized for its mutagenic, poisonous, chemically inert, and non-biodegradable characteristics, rendering it notably dangerous.It poses a significant potential for causing cancer [4].The existence of organic dyes, even in minimal quantities, has the potential to result in a range of health issues, such as gastrointestinal problems, irritations, anemia, and other conditions [5].Hence, it is essential to remove these contaminating dyes from wastewater.
At present, a variety of typical treatment techniques are in use, including methods like chemical precipitation, separation, adsorption, and coagulation [6].However, these methods often lead to incomplete decomposition of dyes and simply shift the pollutants from one phase to another, resulting in the creation of secondary pollutants that demand further management [7,8].In the present day, increased attention is directed toward advanced oxidation processes (AOPs), which utilize semiconducting materials as a viable substitute for traditional approaches [9].The aforementioned procedures offer several benefits, such as reduced equipment demands, indiscriminate oxidation, straightforward operational control, cost efficiency, and complete conversion of organic dyes into harmless byproducts [10].A unique feature of AOPs is the generation of reactive agents like OH, which facilitate the swift and non-selective oxidation of organic contaminants.Notably intriguing is heterogeneous photocatalysis, employing oxide-based nanomaterials, which proves effective in eliminating water-soluble organic pollutants from both water and wastewater under light exposure [11].
Broadly, the process of photocatalytic degradation involves the utilization of suspended particles in a water solution, acting as photocatalysts, which are exposed to illumination.The role of the photocatalyst is pivotal in this procedure.These materials are semiconductors characterized by an electronic band structure, demonstrating a distinct division between the highest occupied energy band, referred to as the valence band (VB), and the lowest unoccupied band, designated as the conduction band (CB), as determined by an energy band gap (Eg).When photons possessing energies equivalent to or surpassing the semiconductor's band gap are absorbed, electron-hole (e − -h + ) pairs are engendered within the semiconductor particles.As these generated carriers move, there is a separation of charges, leading to the generation of various compounds (e.g., H2O, OH, and 1 O2).It is imperative to acknowledge that electrons and holes are also able to recombine without engaging in chemical reactions.The oxidative agents initiate the organic pollutants' decomposition located on or close to the catalytic surface, ultimately culminating in their comprehensive transformation into harmless substances [12].
Metal oxide semiconductors (TiO2, ZnO, etc.) possess remarkable adsorption capacities and serve as effective catalysts due to their elevated reactivity, heightened photosensitivity, extensive specific surface area, cost-efficiency, lack of toxicity, and enhanced catalytic effectiveness in dye photodegradation [13].
In comparison with titanium dioxide (TiO2), zinc oxide (ZnO) possesses a notable advantage as it absorbs a greater portion of the UV spectrum owing to its efficient recombination of electron-hole pairs [13].ZnO exists in three distinct forms: wurtzite, zinc blend, and rock salt [14].Among them, the wurtzite structure is the thermodynamically stable phase under normal conditions.Structurally, wurtzite comprises a hexagonal closepacked arrangement of O 2− anions, with half of the tetrahedral sites occupied by Zn 2+ cations.The remaining half of the tetrahedral sites and all octahedral sites remain unoccupied.Undoped ZnO exhibits n-type conductivity, which is commonly attributed to inherent defects such as excess zinc at interstitial positions and oxygen deficiencies [15].
Numerous investigations have been conducted on the photodegradation of rhodamine B using ZnO [5,[16][17][18].However, directly comparing their photocatalytic performances is challenging, if not impossible, due to the significant variation in parameters employed during the synthesis of ZnO particles, which greatly impacts their effectiveness in dye degradation.For instance, Shidpour et al. [19] conducted a study where they assessed the impacts of various ZnO morphologies, including nanowires, nanorods, spherical particles, and welded nanoparticles, on the degradation of methylene blue dye.Their findings revealed that nanowire-shaped structures displayed the most enhanced degradation efficiency.Similarly, Kumaresan and co-researchers [20] noticed a subtle alteration in ZnO's optical band gap due to particle size variations, leading to a 94% degradation rate for rhodamine B dye and demonstrating reliable stability for reuse.Moreover, Xie and his team [21] reported that ZnO structures resembling sheets exhibited superior photocatalytic activity compared with ZnO's spherical structure.Wang et al. [22] synthesized ZnO with morphologies resembling flowers and rods, concluding that ZnO nanoflowers exhibited higher photocatalytic activity compared with ZnO nanorods.In a different study, Sun et al. [23] fabricated ZnO with 3D hierarchical structures consisting of nanosheets and nanoneedles, with the flower-shaped sample displaying an enhanced photocatalytic performance in comparison with the nanosheet and nanoneedle structures.Efforts have been made to enhance the photocatalytic performance of ZnO through adjustments in preparation condition parameters, such as the utilized Zn precursors, the concentration of Zn precursors, as well as the precipitating agents, the pH and temperature during synthesis, the calcination temperature, and the duration of the synthetic reaction.According to existing studies [24], different Zn precursors present slight or no effects on ZnO nanoparticles' crystallite size, yet they substantially alter their structure.In addition, the divergence in ZnO nanoparticles' crystallite sizes depends on the utilized synthetic approaches [25].The concentration of Zn precursors and the precipitating agents used are factors that possess a vital role in the synthetic route of ZnO nanoparticles, affecting their morphology and crystallite size [26].More specifically, as the concentration of the zinc salt precursors increases, the crystallite size of the ZnO nanoparticles produced also increases.Similarly, the crystallite size of the zinc oxide nanoparticles also increases as the concentration of the precipitating agents increases [24].In another study [27], it was observed that increasing the concentration of the Zn precursor salt beyond a certain limit (> 0.3 M) led to a decrease in the formation of zinc oxide nanoparticles.Similar observations for the optimization of nanoparticles with different metal ion concentrations have also been reported elsewhere [28].Furthermore, the pH value of the reaction solution has a fundamental role in ZnO nanoparticles' synthesis [24,27].Based on conducted studies, ZnO nanoparticles synthesized at lower solution pH values (in acidic conditions), irrespective of the synthetic approach and reaction conditions, are characterized by smaller crystallite sizes compared with those produced in neutral or alkaline conditions.This observation suggests that the acidic environment influences the crystal structure of ZnO, resulting in the formation of smaller crystallites.Conversely, at higher solution pH values, various authors have reported different trends due to the formation of an intermediate compound (zinc hydroxide), which facilitates the development of larger crystallites [24].In general, varying morphologies of ZnO nanoparticles have been observed at different solution pH levels, with the crystallite size being dependent on the pH of the solution [24].Temperature proves to be another crucial factor during the synthesis of ZnO nanoparticles [24,27].In most reported cases, the crystallite size of ZnO nanoparticles tends to grow larger as the synthesis temperature increases.However, it is worth mentioning that there is an exception in a study in which a green synthetic approach for ZnO nanoparticle synthesis was utilized [29].This exception highlights that variations in the crystallite sizes of ZnO nanoparticles are influenced not only by the reaction temperature but also by the specific synthetic approach employed.Additionally, it has been proven that a direct correlation between the calcination temperature and the size of ZnO nanoparticles persists [24].It has been observed that regardless of the chosen synthetic approach, solvent, or other conditions utilized during synthesis, increasing the calcination temperature leads to a proportional increase in the crystallite size of the produced ZnO nanoparticles.The calcination temperature also exerts an influence on the morphology of ZnO nanoparticles, with spherical shapes being predominant in the majority of studies' outcomes.The size of ZnO crystallites varies based on the nucleation time, as reported by various researchers [24].In particular, extended nucleation periods result in the creation of larger crystallite sizes [27].This may be attributed to particles having an ample amount of time to merge, leading to the generation of larger ZnO particles through subsequent fusion [24].
Previous studies have highlighted the influential roles of reducing particle size and increasing crystallinity in the photocatalytic activity of ZnO.However, these investigations have primarily utilized direct methods, limiting the evaluation of a wide range of parameters.In contrast, statistical methods offer the advantage of optimizing procedures without the need for exhaustive exploration of numerous parameters [30].
Various methods are applied to optimize an experimental procedure.The Taguchi statistical method, as well as data mining and artificial intelligence, are typically chosen in cases of material synthetic processes, saving time and money, and considering the research study from an ethical and ecological point of view [31].
The Taguchi statistical method is widely recognized as a powerful engineering approach used to investigate the chemistry of materials, with a focus on minimizing variation and sensitivity to noise.Taguchi designs maintain a balanced nature, ensuring that each factor is independently analyzed without any factor being given undue weight.In general, the Taguchi method may require a greater number of experiments than other approaches; however, it is more accurate and statistically reliable for the optimization of process parameters.In addition, the use of Minitab software, as in the present study, drastically reduces potentially complicated analysis [32].Key concepts in Taguchi experimental designs include factors (the variables under study), levels (the assigned values to the factors), the orthogonal array (OA) (a set of arrays used for the experimental design), such as the L32b array (one two-level factor and nine four-level factors) for a specific number of factors and levels, and the optimum conditions falling into three categories: the maximum value as the preferred outcome, the minimum value as the preferred outcome, and a typical value as the optimal choice.Upon completion of the experiments, analysis is conducted based on the OA to achieve various objectives: identifying the optimal conditions to achieve the highest value, assessing the contribution of each factor to the response, and predicting the response under the optimized conditions [33,34].
This study focused on investigating the synthesis of ZnO nanoparticles using the precipitation method.The Taguchi method was employed to optimize the preparation conditions, aiming to achieve a smaller particle size and higher crystallinity, which are known to enhance photocatalytic activity for the degradation of rhodamine B under UV irradiation.The optimized conditions were determined and applied to prepare a sample, which was subsequently examined to assess its characteristics.Moreover, the photocatalytic mechanism involved in the degradation of rhodamine B was investigated using the sample that exhibited the optimum characteristics.The present research is considered a turning point for all researchers who aim to produce pure zinc oxide nanoparticles using the precipitation method, characterized by an exceptionally small crystallite size (3.29 nm) and enhanced photocatalytic properties in the degradation of organic dyes under UV irradiation (100% RhB degradation within 90 min).In addition, the unique feature of this research is the optimization of zinc oxide nanoparticle synthesis, examining quite a large number of parameters involved in the utilized synthetic approach via a facile statistical method.
Moreover, the crystallinity index (CI %) was calculated based on Equation (2) [40]: Total area of all crystalline peaks Total area of all crystalline and amorphous peaks No significant shifts in the peaks' positions were observed, as supported by the acquired XRD diffractograms of the prepared ZnO samples.Nevertheless, alterations in the peaks' intensities in the different samples can be observed, affecting the FWHM and thus resulting in different crystallite sizes.
Among all the indexed planes, that of (101) indicated the most enhanced relative intensity in the entire XRD pattern, implying anisotropic growth and the preferred orientation of the crystallites, given the fact that epitaxial growth along the C-axis in the (001) direction constitutes a distinctive phenomenon in materials possessing a wurtzite structure [41].
Based on the obtained results, the sample annotated as S3 owned the smallest average crystallite size (3.29 nm) and the highest crystallinity (79.83 %) among the produced ZnO samples.

Diffuse Reflectance UV-Vis Spectroscopy analysis (DRS)
Band gap energy (Eg) comprises a critical parameter that needs to be taken into consideration in photocatalytic studies.All examined ZnO samples presented absorbance bands, characterized by an absorption edge of <385 nm, which is habitual for the wurtzite crystal phase of ZnO.These bands correspond to the O2p→Zn3d electron transition from the VB to the CB [42].The DRS measurements of the produced ZnO samples are illustrated in Figure 2a-d To assess the reflectance of the prepared ZnO samples, the Kubelka-Munk (K-M) method was employed, as shown in Figure 2a-d, using Equation (3) [43]: where R is the reflectance.Figure 3a-d illustrates the direct energy band gaps of the synthesized ZnO samples using the K-M model vs. energy via the extrapolation of the (F(R)hv) 1/2 vs. hv spectra's linear region.The energy band gap (Eg) was estimated using Tauc's equation (Equation ( 4)): where Eg represents the energy band gap, h stands for Planck's constant, v corresponds to the frequency, α is the absorption's coefficient, and n is equal to ½ [44].The band gaps of the prepared ZnO samples (S1-S32) were calculated to be between 3.37 eV and 3.41 eV (Table 4), respectively.As can be observed, the estimated values present negligible alterations.However, blue shifts from 3.37 eV (S3) of the ZnO nanospheres to 3.40 eV (S12) of the ZnO nanoflakes (see Figure 6) could be possibly attributed to the quantum confinement of the ZnO nanostructures due to the fact that the flakes' mean thickness is comparable to the Bohr radius of ZnO [21].The S3 sample achieved the lowest Eg value (3.37 eV) among all studied samples.

Photocatalytic Activity of the Prepared ZnO Samples
The photocatalytic efficiency of the prepared ZnO samples was evaluated by measuring their efficiency in decomposing rhodamine B in an aqueous solution upon UV light irradiation.The photocatalytic trials were performed at 25 °C and a pH equal to 6.64 ± 0.01.
Based on the results obtained via the RhB degradation measurements, the acquired data of the kinetic model investigations under UV light exposure are presented in Figure 4a-d, indicating the relationship between the variation in −ln(C/C0) versus time for the examined ZnO powders.The rate of photocatalytic RhB absorption on the surface of the powders was estimated according to the pseudo-first-order kinetics model, which can be described using the following equation (Equation ( 5)) [42]: where C0 and C correspond to the RhB's initial and reaction time concentrations, respectively, k represents the apparent rate constant of the photocatalytic oxidation, and t is the time of irradiation.The slope of the linear fitted plot stands for the apparent rate constants of the examined ZnO samples (S1-S32) (Table 5).

Results Analysis Based on the Taguchi Approach
Using the provided data in Tables 1 and 2 and following the principles of the Taguchi method, Figures 5 and 7-9 are provided to demonstrate the impacts of the specified variables on various outcomes: the average size of the crystallites, the crystallinity, the energy band gap (Eg), and the apparent rate constants (k) in the degradation of RhB.Additionally, the Minitab software 21.1.0(Triola Stats, Wappingers Falls, NY, USA) was employed to calculate the percentage of influence and rank the factors with respect to the analyzed parameters, as presented in Table 6.The size of crystallites signifies a group of atoms within a crystal that share the same orientation.It represents the dimension of a coherent diffracting domain in the crystalline structure.The formation of point defects and dislocations in the synthesized crystals plays a role in determining their size [45].Hence, the average crystallite size is a significant variable that impacts both the physical and photocatalytic properties [46].More specifically, based on existing studies, a small average crystallite size leads to enhanced photocatalytic performance [40].
According to Figure 5, the type and concentration of the Zn precursor, the concentration of the precipitation agent, the pH of the solvent, as well as the calcination temperature, exerted the greatest influence on the average crystallite size.Notably, zinc nitrate and zinc sulfate were responsible for the largest crystallite size, in accordance with our previous research [42].To validate this observation, FESEM images of four specific samples, namely, sample S3 (containing zinc acetate), sample S5 (containing zinc chloride), sample S12 (containing zinc nitrate), and sample S32 (containing zinc sulfate), were captured and are presented in Figure 6a-d.These images corroborate the findings regarding crystallite size.Samples S3 and S5 exhibit semi-spherical particles with a particle size distribution of 50-100 nm, while samples S12 and S32 exhibit significant changes, with the particles experiencing extensive growth and forming micron-sized hexagonal flakes.The pronounced growth of the particles originating from the zinc nitrate and zinc sulfate precursors contributed to the larger crystallite size observed in this series of samples.Pourrahimi et al. [47] conducted a study on the impact of zinc salts on the grain size of ZnO particles.Based on the obtained results, during ZnO crystallization, the oriented self-assembly and condensation of nanoprisms led to the formation of petals, which were further shaped into octahedrons that interconnected to form polygons.Subsequently, these polygons underwent further growth into larger particles in the case of nitrate and sulfate through the migration and addition of crystallites via dissolution and reprecipitation.Other zinc salts such as chloride and acetate exhibit a stronger ability to prolong the stabilization of the formed crystallites.Additionally, the average crystallite size was found to increase with the increasing concentration of the Zn precursor from 0.1 M to 0.4 M, as proposed in the relevant literature [48], as well as with the increase in the calcination temperature from 400 °C to 700 °C, due to the provision of more energy for nucleation and particle growth [48].
The pH value constitutes another parameter that significantly impacts the average size of the crystallites.Elevating the pH value from 11 to 13 resulted in a decrease in the average crystallite size.At a pH value equal to 11, the solution contained various zinc hydroxide complexes, with Zn(OH)2 being the predominant one, leading to the formation of larger particles.When the pH was increased to 13, the condensation process was delayed, resulting in the deceleration of the particles' growth [49].
Finally, it was also observed that the increasing concentration of the examined precipitation agents led to a significant decrease in the average crystallite size, which is in accordance with existing studies' results [50].

Crystallinity
Based on the results depicted in Figure 7 and the corresponding analysis, the factors exerting the greatest impact on the degree of crystallinity were (a) the type and concentration of the Zn precursor, (b) the concentration of the utilized precipitating agent, (c) the pH value of the solvent, and (d) the defined temperature during the calcination of the prepared ZnO samples (Table 3).As for the type of the utilized Zn precursor, the samples that derived from zinc acetate indicated the highest mean crystallinity (74.34%), followed by those prepared with zinc nitrate (73.88%), zinc chloride (73.58%), and zinc sulfate (72.09%), results that are in accordance with a previously conducted study [42].Additionally, the mean degree of crystallinity presented a gradual enhancement during the elevation of the Zn precursor's concentration, as proposed in the relevant literature [51].
Moreover, the optimal level of crystallinity was attained by adjusting the pH to 14. Increasing the pH value from 11 to 14 resulted in increased crystallinity.This behavior can be attributed to the alteration in the solution component.At higher pH values, ZnO crystallites can form and remain stable [52], which explains the higher crystallinity observed at a pH of 14.Under lower pH conditions, the synthesized ZnO crystallites become unstable, leading to the synthesis of particles with lower crystallinity [53].The highest pH value facilitates the formation of particles with the highest crystallinity, as it promotes the presence of more Zn(OH)x ions, which serve as the building blocks of the ZnO structure [52].
The observed increase in the calcination temperature promoted the decrease in the crystallinity degree [54].

Energy Band Gap (Eg)
Despite the minor variation in Eg values, the utilization of zinc acetate in a solution with a pH value equal to 13 and subjected to thermal treatment at 70 °C for 90 min yielded the lowest value.Among the factors considered, the type of Zn precursor exhibited the most significant influence on this particular parameter (Figure 8).

Apparent Rate Constant (k)
The goal of the present study was to produce ZnO nanoparticles via the precipitation method, characterized by enhanced photocatalytic activity under UV light irradiation.Thus, in order to assess the photocatalytic performance of the synthesized ZnO samples, kinetic model studies were conducted, and the photocatalytic capabilities of the samples were represented as k (Figure 9).The sample that exhibited the highest k value was obtained by utilizing zinc acetate (0.3 M) as the precursor in a methanol solution, with a pH equal to 13 that was adjusted using sodium hydroxide (1.5 M), followed by stirring at 70 °C and 700 rpm for 90 min, and, finally, calcination at 600 °C (see Table 2).The choice of zinc acetate as the Zn precursor significantly influenced the size of the crystallites, Eg values, and ultimately the value of k.Optimal results, including the smallest crystallite size, highest crystallinity, and lowest Eg value, were obtained when this particular precursor was used.These factors collectively contributed to an enhanced photocatalytic behavior, while the highest photoactivity could be achieved when all of them were simultaneously provided.The remaining examined parameters (see Table 1) indicated similar impact rates (see Table 5).

Synthesis of the Optimized ZnO Sample
Initially, a 0.3 M zinc acetate solution was prepared via the addition of zinc acetate dihydrate to methanol under constant magnetic stirring (700 rpm).The starting pH was equal to 6.43 ± 0.01.Then, sodium hydroxide (1.5 M) was poured dropwise into the solution until a pH value equal to 13 was obtained.After 90 min at 70 °C, the obtained milky white solution was centrifuged to acquire the white precipitate, which was subsequently triturated and purified via rinsing and centrifugation in order to eliminate potential impurities.The acquired precipitate remained at 100 °C for 1 h to dry and was then calcinated at 600 °C for 3 h, finally leading to the production of a white powder (Figure 10).The reaction that occurred during the synthetic procedure is described using the following chemical reaction (Equation ( 6)) [42]:

Characterization of the Optimized ZnO Sample XRD Analysis
The annotated peaks in the derived spectrum (Figure 11) are in complete accordance with that of bulk ZnO ((JCPDS) card no.36-1451), as mentioned previously.The absence of impurities was confirmed, as no other peaks were noticed in the spectrum.The high intensity as well as narrowness of the depicted diffraction peaks confirm the highly crystalline structure of the produced ZnO particles [55].The indexed peaks at 2θ = 31.76°,34.43°, 36.27°,47.56°, 56.64°, 62.89°, 66.53°, 67.99°, 69.11°, 72.63°, and 76.98° are attributed to the (100), ( 002), ( 101), ( 102), ( 110), ( 103), ( 200), ( 112), ( 201), (202), and (004) crystal planes, respectively.Calculations for the determination of the crystal lattice indices and interplanar d-spacing were additionally conducted, besides the ascertainment of the crystallinity and average crystallite size that were previously studied.The obtained data are presented in detail in Tables 7 and 8.The average crystallite size as well as the degree of crystallinity of the optimized ZnO sample were determined utilizing Equations ( 1) and (2), while the interplanar d-spacing was estimated using Bragg's law equation (Equation ( 7)), as described below: Furthermore, the crystal lattice indices were estimated using Equation ( 8) [42]:

Micro-Raman Analysis
The Raman spectrum acquired for the optimized ZnO sample is illustrated in Figure 12.The Raman features in the studied ZnO sample correspond to the Raman active modes of the ZnO wurtzite crystal [56].The Raman peak presenting the highest intensity in the obtained spectrum is spotted at approximately 440 cm −1 (E2H mode), imputed to oxygen vibration [56].The pronounced lack of symmetry is connected to the disorder in the crystal structure and the interactions between non-harmonic phonons (P-P) [36].Unlike the E2 phonons, the polar phonons A1 and E1 split into transverse optical (TO) and longitudinal optical (LO) phonons.Consequently, this distinct peak represents a Raman active feature specific to the wurtzite hexagonal phase of ZnO [57].
The second-order mode observed in the low-wavenumber region around 335 cm −1 is attributed to the difference between E2H and E2L [36].In the middle range of the spectrum for non-activated ZnO, the second-order Raman modes possessing A1 symmetry (484, 701, 725, and 742 cm −1 ) are relatively weak and not easily distinguishable, except for the mode observed at 541 cm −1 .Additionally, acoustical and optical modes integration occur at ≈ 660 cm −1 .

FT-IR Analysis
In the FT-IR spectrum (Figure 13), the presence of ZnO is evident based on the peak detected at 554 cm −1 .The peak at 3478 cm −1 corresponds to inter-hydrogen bonding, indicating the presence of water molecules and hydroxyl groups.It was likely influenced by atmospheric humidity during the measurement of the sample [58].
The peak observed at 1113 cm −1 is associated with the stretching of C-O bonds in primary alcohols.Two additional distinct peaks are particularly noticeable, at 1520 cm −1 , which corresponds to the vibration of a C=O bond and could be indicative of the presence of organic residues that may have persisted after the reaction with acetate, as well as the peak at 1427 cm −1 , possibly corresponding to the bending of carboxylic acid O-H bonds.Furthermore, the identified peak at 2351 cm −1 is attributed to the stretching of carbon dioxide O=C=O bonds [58].

N2-sorption Analysis
The N2-sorption isotherm for sample S3 is presented in Figure 14.It is evident that the synthesized optimized ZnO sample exhibits a type-IV isotherm, with a narrow hysteresis loop and no saturation plateau, which indicates the presence of mesopores and macropores.The inset in Figure 14 depicts the corresponding pore size distribution, which is determined from the desorption curve using the BJH method.The pore size distribution is quite broad and extends in the mesopore (2-50 nm) as well as in the macropore area (>50 nm), in agreement with the N2-sorption isotherm [59][60][61].
The data of the physical parameters, including the BET surface area, micropore surface area, cumulative volume, and average pore diameter, are compiled in Table 9.The optimized ZnO sample demonstrates an enhanced BET surface area that aligns with the findings of the small crystallite size, as determined from the XRD results (see Table 7).Table 9. Results acquired using the BET method.(a) Specific surface area calculated using Brunauer-Emmett-Teller theory, (b) micropore surface area via t-plot analysis, according to the Harkins and Jura model, (c) cumulative volume of pores between 1.7 and 300 nm from N2-sorption data and the BJH desorption method, and (d) average pore diameter, calculated using the 4 V/σ method.V was set equal to the maximum volume of N2 adsorbed along the isotherm as P/Po → 1.0.The DLS measurements were conducted at a pH value equal to 6.43 ± 0.01. Figure 15 illustrates the hydrodynamic radius distribution as a function of the scattering light intensity for the synthesized optimized ZnO sample.The sample exhibited a favorable distribution of particle sizes, showing a nearly monodisperse nature (PdI ≈ 0.1) (Table 10).The size distribution of the particles ranges between 10 and 100 nm, with a peak value of approximately 26 nm.High-purity ZnO nanoparticles, prepared using acetate salt as the precursor and characterized by a mean size of approximately 25 nm, have been previously reported [62], aligning with the experimental findings of the present study.Table 10 presents the obtained data of the DLS measurements.Diffuse Reflectance UV-Vis Spectroscopy Analysis (DRS)

Sample
The DRS results for the synthesized sample are depicted in Figure 16a,b.The reflectance of the optimized ZnO sample was evaluated using the Kubelka-Munk (K-M) method, as shown in Figure 16a.Equation ( 3) was employed for this purpose.
As presented in Figure 16a, the absorption edge is below 385 nm, which is typical for the wurtzite crystal phase of ZnO, as previously mentioned.Furthermore, Figure 16b demonstrates the direct energy band gap of the produced ZnO sample, utilizing the K-M model vs. energy via the extrapolation of the linear region of the spectrum (F(R)hv) 1/2 vs. hv.Eg was evaluated via Tauc's equation (see Equation ( 4)).

FESEM Analysis
The main morphological aspects of the prepared optimized ZnO sample were evaluated via FESEM observation and are depicted in Figure 17a-c.According to the acquired observation results, the obtained nanoparticles exhibit a combination of spherical and hexagonal shapes.It is widely recognized that zinc oxide (ZnO) can exist in various structures, such as zero-dimensional (0D), one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) forms.Among these, hexagonal and diamond-shaped structures are the most commonly observed, which aligns with the findings reported in the relevant literature [63].

Study of the Photocatalytic Efficiency in Degradation of Rhodamine B
The effectiveness of the optimized ZnO sample in photocatalysis was evaluated by measuring the degradation of rhodamine B (RhB) dye in an aqueous solution under UV light exposure.All photocatalytic experiments were conducted at ambient temperature (25 °C) and a pH value of 6.43 ± 0.01. Figure 18 illustrates the photocatalytic performance of the investigated sample.Control experiments were also performed to ascertain photolysis (RhB photolysis) and adsorption-desorption equilibrium (RhB dark) without light irradiation (in dark conditions) and under continuous stirring, maintaining the same duration as the photocatalytic process.The observed degradation of rhodamine B under UV light irradiation was found to be <1%, indicating that RhB's degradation rate in the absence of the studied photocatalyst is remarkably low.Furthermore, consistent outcomes were acquired from the experiments that took place under dark conditions, confirming the stability of the dye [42].It can be observed that during the photocatalytic trials, the studied optimized sample indicated a high efficiency, attaining total degradation (100%) of the RhB dye within 90 min under UV light irradiation.
Furthermore, Figure 19 displays the UV-visible spectra recorded during the photocatalytic experiments.These spectra were utilized to monitor the degradation rate of RhB over time, investigate its degradation mechanism, as well as assess the photocatalytic efficiency of the sample.Generally, two pathways for RhB degradation are known: (a) Nde-ethylation and (b) disruption of its conjugated structure.When the first degradation mechanism occurs, a blue shift in the maximum absorption is observed, while in the second mechanism, there is no significant blue shift but rather a gradual decrease in absorption over time [64].In the second proposed degradation pathway, the maximum absorption shifts from 554 nm (which is a characteristic peak of RhB) to 498 nm (blue shift).However, the relevant literature suggests that both pathways contribute simultaneously to the degradation of RhB [65].The real-time UV-visible spectra obtained for the optimized ZnO sample during the photocatalytic trials under UV light irradiation (Figure 19) evidently demonstrate the presence of both pathways, highlighting their significant roles in the degradation of RhB.To validate the findings of RhB's photocatalytic degradation investigations, an additional analysis using TOC was performed in order to assess the percentage of RhB mineralization achieved during the photocatalytic experiments.The mineralization percentage of RhB was calculated using Equation ( 9): where TOCinitial comprises the total organic carbon concentration in the medium before the photocatalytic study, and TOCfinal stands for the total organic carbon concentration in the medium after the end of photocatalysis [41].

Kinetic Model Study
Figure 20a depicts the results obtained from the pseudo-first-order kinetic model study under UV light irradiation, as a variation of −ln(C/C0) vs. time for the mentioned sample that is described using Equation (5).However, the kinetics of photocatalytic trials may also be described using the pseudo-second-order equation (Equation ( 10)).Compared with the pseudo-first-order kinetics (Figure 20a), the R 2 values of the pseudo-second-order kinetic model study (Figure 20b) are significantly poorer.Table 11 illustrates the kinetic parameters of the studied optimized ZnO sample.
where qt and qe represent the pollutant amount adsorbed at time t and equilibrium, respectively (mg/g), and k2 comprises the rate constant (g•mg −1 •min −1 ).The photocatalytic procedure involving RhB dye can be outlined using the following Equations ( 11)-( 16) [42]: In the photocatalytic oxidation procedure, various key oxidative species play significant roles, which include superoxide radicals (O 2− ), electrons (e − ), holes (h + ), and hydroxyl radicals (OH).To achieve a deeper understanding of the fundamental mechanism underlying photocatalysis, extensive investigations were conducted to pinpoint the active species.This was accomplished through a series of experiments designed to capture and isolate specific species.Notably, p-benzonquinone (p-BQ), silver nitrate (AgNO3), disodiumethylene diaminetetraacetate (EDTA-2Na), and t-butanol (t-BuOH) were introduced into the RhB dye solution.These substances were specifically employed to capture and identify superoxide radicals (O 2− ), electrons (e − ), holes (h + ), and hydroxyl radicals (OH), respectively [66].
According to the findings presented in Figure 21, the RhB degradation efficiency of the optimized ZnO sample underwent a significant decrease to 34.08 ± 1.25% upon the addition of t-BuOH into the reaction solution.These data provided confirmation that the hydroxyl radicals (OH) exerted a key influence on RhB's photocatalytic degradation.Additionally, the inclusion of p-BQ, AgNO3, or EDTA-2Na in the photocatalytic reaction solution led to a decrease in RhB's degradation effectiveness from 100% to 96.34 ± 1.01%, 98.06 ± 1.13%, and 95.47 ± 1.08% respectively, indicating that the superoxide radicals (O2 − ), photogenerated electrons (e − ), and holes (h + ) were not the primary reactive species involved in the process.Based on the scavenging experiment's results, a plausible mechanism is proposed (Figure 22).When exposed to UV light irradiation, ZnO becomes excited, leading to the simultaneous production of e -and h + .The holes reside in the VB, while the electrons are located in the CB.The electrons in the CB readily take part in reactions with dissolved O2 in H2O, leading to the generation of O2 − , H2O2, or OH radicals.At the same time, the photoproduced h + in the VB is able to directly interact with OH -ions or H2O molecules bound to the surface, leading to the formation of OH radicals.These highly reactive species, such as OH radicals, relentlessly attack the molecules of rhodamine B until thorough decomposition is achieved, promoting the effective separation of e − -h + pairs and, as a result, enhancing the photocatalytic activity.In the process of photodegrading organic pollutants, •OH radicals participate in a redox reaction, transforming long-chain molecules into H2O and CO2 [67].Nair and colleagues conducted experiments involving both undoped and cobalt-doped ZnO nanoparticles for photocatalysis against a methylene blue (MB) solution.Their findings indicated that pure ZnO demonstrated significantly greater photocatalytic efficiency compared with cobalt-doped ZnO.This enhanced photocatalytic activity can be attributed to the presence of zinc vacancies and interstitial oxygen atoms, creating acceptor states, as well as oxygen vacancies and interstitial zinc atoms, forming donor states [68].In another study, Kaur and Singhal synthesized various metal-doped ZnO nanoparticles and assessed their photocatalytic performance against methyl orange (MO).The results revealed that the degradation of MO followed first-order reaction kinetics, with undoped ZnO nanoparticles outperforming their metal-doped counterparts [69].Furthermore, nanoscale zinc oxide photocatalyst powders were synthesized using the precipitation method, and their photocatalytic effectiveness was evaluated by observing the degradation of MB dye as a test case.The maximum degradation of MB occurred under UV radiation, reaching 81% within 3 h when using photocatalysts prepared via the precipitation method.Assuming a pseudo-first-order reaction kinetics model, this equated to a rate constant of 8.4 × 10 −3 min −1 [4].In a study by Sultana and colleagues [70], the generation of hydroxyl radicals using zinc oxide (ZnO) synthesized from zinc nitrate (Zn(NO3)2) and sawdust was investigated to degrade MO via photocatalysis.Their experimental setup involved a 5 × 10 −3 M tannic acid (TA) solution, and they exposed a suspension containing 1 mg of photocatalyst/mL to UVB light (290-320 nm).The results of their research revealed that pure ZnO, without the inclusion of sawdust, exhibited negligible photoactivity and a considerably lower formation of hydroxyl radicals, even when subjected to high-energy UVB irradiation.Conversely, in another investigation [71], ZnO nanowires synthesized via co-precipitation were examined for the degradation of Acid Red 57 dye.This evaluation was conducted using UVA, UVB, and UVC irradiation while varying the calcination temperature.The researchers determined that the optimal degradation performance occurred when calcining at 400 °C and using UVC light (280-100 nm).Moreover, they identified that the primary mechanism responsible for the photodegradation was the generation of hydroxyl radicals.

Reusability Study
As depicted in Figure 23, the capability of the optimized ZnO sample to be reused was demonstrated by subjecting it to UV light exposure for five consecutive photocatalytic cycles (catalyst loading = 5 mg; pH = 6.43 ± 0.01; C0(RhB) = 10 mg/L).After each cycle of degradation, the sample was treated with centrifugation and multiple rinses with distilled water.Following this, it was dried in a vacuum oven at 70 °C for 24 h in preparation for the subsequent trial, without any additional treatment [72].Importantly, the ZnO sample synthesized in this manner exhibited remarkable resistance to degradation under UV light, showing a ≈ 7% reduction (6.78 ± 0.91%) in its photocatalytic efficiency after undergoing five consecutive cycles.These findings provide proof of the robust stability of the optimized ZnO photocatalyst throughout the repeated cyclic photocatalysis process under UV light irradiation.

Experimental Procedure Design Based on the Taguchi Approach
The L32b array of the Taguchi method was utilized to design the experimental plan of the present study.For this reason, ten variable parameters (the precipitating agent type (PAT), precursor of Zn-Znprec, the concentration of Zn precursor (CZn), concentration of precipitating agent (CPA), solvent type (Stype), solvent pH (SpH), temperature during synthesis (Tsyn), temperature during calcination (Tcalc), stirring time (Stt), and stirring speed (Sts)) were selected.Nine were adjusted in four levels and one was adjusted in two levels.These variable parameters and their corresponding values are presented in Table 1.
Given the L32b array of the Taguchi method, the synthesis conditions were specified and are presented in Table 2.The average crystallite size, crystallinity, and band gap energy (Eg) were used as intermediate responses, while the photocatalysis degradation rate constant (k) was estimated as the response parameter.The optimization of the mentioned parameters was conducted via the comparison of the mean value.The optimized conditions were utilized to synthesize and fully characterize a ZnO sample.Figure 24 describes the main steps of the optimization process using a flowchart.

Figure 24.
Flowchart that describes the main steps of the optimization process.

Preparation of the ZnO Samples
Thirty-two (32) samples were prepared according to the conditions described in Table 11.Firstly, the Zn precursor and the proper solvent (double-distilled water (DDI-H2O), distilled water (DI-H2O), methanol (CH3OH), and ethanol (CH3CH2OH)) were mixed and magnetically stirred at the proper speed, according to Table 2, to gain a uniform solution with a concentration equal to 0.15 M. According to the procedure [42], the appropriate precipitating agent (NaOH or KOH) was added dropwise to the point where the solution reached the relevant pH value.After stirring at the proper time, speed, and temperature, the solution was centrifuged for 10 min at 9000 rpm in order to recover the white precipitate.The precipitate was washed with distilled water and centrifuged thrice in order to remove any organic remains.Each precipitate was left at 100 °C for 1 h to dry, and then, calcination was realized at the appropriate temperature (see Table 2) for 3 h to obtain powder samples.The final products, for convenience purposes, will be called hereinafter S1-S32 (see Table 2).

Characterization of the Prepared ZnO Samples
The X-ray diffractometer that was utilized to study the average crystallite size and crystallinity of the synthesized ZnO samples was a Brucker D8 Advance (Brucker, Berlin, Germany) with CuΚα radiation (λ = 1.5406Å) at 40 kV and 40 mA.The experimental conditions used were as follows: (a) a 2-theta (2θ) angle ranging from 20 and 80° and (b) a counting diffraction intensity step equal to 0.01°/1.0s.Additionally, measurements under the same experimental conditions were conducted using a NIST standard LaB6 powder (Malvern Panalytical Ltd., Malvern, UK) as the reference material in order to calculate the instrumental broadening considered for the Scherrer equation calculations.
The energy band gap (Eg) values for all prepared ZnO samples were estimated using an ultraviolet-visible (UV-Vis) spectrometer (Jasco UV/Vis/NIR Model name V-770, Interlab, Athens, Greece) equipped with an integrating sphere (wavelength accuracy = ±0.3nm), permitting the recording of diffuse reflectance measurements (DRS).For the measurements, the samples were used as received without any further pre-treatment.
Then, the optimized ZnO sample was thoroughly characterized using various techniques.For the XRD and DRS analyses, the same conditions and equipment were used as mentioned above.The morphology of the optimized ZnO sample was analyzed using field-emission SEM (FESEM, JSM-7401F, JEOL, Tokyo, Japan).Observations via fieldemission SEM were also conducted to complement the obtained statistical results based on the Taguchi analysis.For the measurements, the samples were used as received without any additional preparation.A beam voltage equal to 1.5 kV was set, while the resolution was equal to 1 nm.
A micro-Raman apparatus (inVia, Renishaw, Wotton-under-Edge, Gloucester, UK) provided with a green laser line (λ = 532 nm) as the excitation source with an average output power of approximately 50 mW was utilized in the present study.Raman spectra were recorded at ambient temperature, and an internal Si reference was used as the calibration of the frequency shift value.For each measured powder, 3 to 10 spots were acquired, and for recording, the following conditions were set: exposure time, 30 s; accumulations equal to 3; power, 1%; spectral range, 100-1000 cm −1 ; and the spectral resolution for all measurements was 2 cm −1 .
In order to complement the micro-Raman results, FTIR spectra were also acquired in the 400-4000 cm −1 range, with a resolution equal to 4 cm −1 at 25 °C with an FTIR JASCO4200 (Oklahoma City, OK, USA) instrument possessing a Ge crystal.
The N2 adsorption of the synthesized powders was conducted with a ChemBET 3000 Instrument (Yumpu, Diepoldsau, Switzerland) to ascertain the specific surface area BET.Prior to each measurement, the examined samples underwent degassing at 80 °C for 24 h.Dynamic light scattering (DLS) was used to estimate both the hydrodynamic diameter as well as its distribution for the studied powders' particles that were immersed in distilled water solutions.The incident light was a 633 nm laser, and a 173° scattering angle was arranged for the scattering intensity recording (Malvern Zetasizer Nano ZS, Malvern Panalytical Ltd., Malvern, UK).
RhB's mineralization was estimated by conducting total organic carbon (TOC) analysis (TOC-LCSH/CSN, Shimadzu Scientific Instruments, Columbia, MD, USA) for the exact time points as the ones used during the photocatalytic experiments.

Photocatalytic Activity Study of the Prepared ZnO Samples
The UV light photocatalytic efficiency of the prepared ZnO samples was studied for rhodamine B decomposition, utilizing 0.005 g of each examined sample in a 10 mg•L −1 rhodamine B (C28H31CIN2O3, Penta-Chemicals Unlimited, Prague, Czech Republic) aqueous solution (250 mL) at ambient temperature and a pH equal to 6.64 ± 0.01.Before each photocatalytic experiment (cycles), saturation of the RhB solution for 1 h was carried out through extra pure O2 (99.999%) flow.
The utilized photoreactor possessed a system of four parallel lamps at a distance of 10 cm above each powder's surface [73].Blacklight lamps of 368 nm, 830 lumens, and an incident light flux equal to 0.184 μmol quanta/s made by Sylvania (Sylvania, Wilmington, NC, USA) were used as the source of UV light, and all the experiments were carried out under ambient temperature conditions.
The degradation of RhB via photocatalysis was evaluated in the presence of UV light.The absorbance of the ZnO samples prepared in this manner was measured at 554 nm using a Thermo Fisher Scientific Evolution 200 Spectrometer (Thermo Fisher Scientific, Waltham, MA, USA).The C/C0 ratio, where C represents the concentration of RhB after a specific period of photocatalysis and C0 is the initial RhB concentration, was determined indirectly by comparing the measured A (absorption at each time) with the initial A (Ainitial) [73].

Conclusions
The primary objective of this research was to examine the production of ZnO nanoparticles using the precipitation method.To attain a smaller particle size and improved crystallinity, which are recognized as enhancing the photocatalytic effectiveness in degrading rhodamine B under UV irradiation, the Taguchi method was employed to optimize the preparation conditions, more specifically, the Taguchi L32b orthogonal array using input variables (Zn precursor, Zn precursor's concentration, precipitating agent type, precipitation agent's concentration, solvent type, solvent pH, synthesis temperature, calcination temperature, stirring time, and stirring speed) and response variables (average crystallite size, crystallinity, energy band gap, and photodegradation constant).Subsequently, the determined optimal conditions were implemented to create a ZnO sample, which was thoroughly characterized utilizing XRD, micro-Raman analysis, FTIR, BET, DLS, DRS, and FESEM to evaluate its properties.
Furthermore, this study delved into investigating the photocatalytic mechanism associated with the degradation of rhodamine B, utilizing the sample that demonstrated the most favorable characteristics.The optimized ZnO sample possessing a hexagonal wurtzite crystalline structure achieved 100% degradation of the RhB dye within 90 min under UV light irradiation, while the pseudo-first-order reaction kinetics provided a realistic description of the photocatalytic degradation of rhodamine B given the calculated R 2 values (0.989).

Figure 5 .Figure 6 .
Figure 5. Effects of the design parameters on the average crystallite sizes of the prepared ZnO samples.

Figure 7 .
Figure 7. Effects of the design parameters on the crystallinity of the prepared ZnO samples.

Figure 8 .
Figure 8. Effects of the design parameters on the energy band gap (Eg) values of the prepared ZnO samples.

Figure 9 .
Figure 9. Effects of the design parameters on the apparent rate constant (k) values of the prepared ZnO samples.

Figure 10 .
Figure 10.Schematic representation of the synthesis of the optimized ZnO sample.

Figure 11 .
Figure 11.XRD diffractogram of the optimized ZnO sample (S3) obtained using the Taguchi statistical approach.

Figure 12 .
Figure 12.Raman spectrum of the optimized ZnO sample.

Figure 13 .
Figure 13.FT-IR spectrum of the optimized ZnO sample.

Figure 14 .
Figure 14.N2-sorption diagram of the optimized ZnO sample, as indicated (sorption: black line; desorption: red line).The inset shows the pore size distribution using the BJH method.

Figure 16 .
Figure 16.(a) F(R) reflectance as a function of the wavelength for the optimized ZnO sample and (b) energy band gap (Eg) of the same sample.

Figure 18 .
Figure 18.Degradation curve of RhB for the optimized ZnO sample vs. time under UV light irradiation.In addition, photolysis and degradation of RhB in the dark are encompassed.

Figure 19 .
Figure 19.Real-time UV-visible spectra obtained under UV light photocatalytic degradation of RhB in the presence of the optimized ZnO sample.

Figure 20 .Table 11 .
Figure 20.Photocatalytic kinetic model studies for the optimized ZnO sample, following (a) the pseudo-first-order model and (b) the pseudo-second-order model, under UV light photocatalytic trials.Table 11.Kinetic parameters of the studied sample under UV light photocatalysis.

Figure 21 .
Figure 21.Scavenging experiments for the degradation of RhB in the optimized ZnO sample under UV light irradiation.

Figure 22 .
Figure 22.The proposed photocatalytic mechanism of the optimized ZnO sample toward degradation of RhB under UV light irradiation.

Figure 23 .
Figure 23.Reusability efficiency of the optimized ZnO sample after 5 experimental photocatalytic cycles under UV light irradiation.

Table 1 .
Parameters and their relative values for the experimental procedure design.

Table 2 .
Synthesized ZnO samples based on the Taguchi L32b orthogonal array.

Table 3 .
Average crystallite size, FWHM, and crystallinity of the synthesized ZnO samples.

Table 5 .
Apparent rate constant (k) values of the examined ZnO samples upon UV light photocatalytic trials based on the pseudo-first-order kinetic model study.

Table 6 .
Influence percentages and ranks of the factors for the studied parameters.

Table 7 .
Crystal lattice indices, average crystallite size, FWHM, and crystallinity of the optimized ZnO sample.
* The crystallite size was estimated based on the (101) plane's peak.

Table 8 .
D-spacing calculations for the optimized ZnO sample.