Cerium Oxide Nanorods Synthesized by Dalbergia sissoo Extract for Antioxidant, Cytotoxicity, and Photocatalytic Applications

In this study, cerium oxide nanorods (CeO2-NRs) were synthesized by using the phytochemicals present in the Dalbergia sissoo extract. The physiochemical characteristics of the as-prepared CeO2-NRs were investigated by using ultraviolet-visible spectroscopy (UV-VIS), scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR), and X-ray diffraction analysis (XRD). The SEM and UV-VIS analyses revealed that the acquired nanomaterials possessed a rod-like morphology while the XRD results further confirmed that the synthesized NRs exhibited a cubic crystal lattice system. The antioxidant capacity of the synthesized CeO2-NRs was investigated by using several in vitro biochemical assays. It was observed that the synthesized NRs exhibited better antioxidant potential in comparison to the industrial antioxidant of the butylated hydroxyanisole (BHA) in 1,1-diphenyl-2-picrylhydrazyl (DPPH) assay. The biochemical assays, including lipid peroxidation (LPO), total antioxidant capacity (TAC), and catalase activity (CAT), were also performed in the human lymphocytes incubated with the CeO2-NRs to investigate the impact of the NRs on these oxidative biomarkers. Enhanced reductive capabilities were observed in all the assays, revealing that the NRs possess excellent antioxidant properties. Moreover, the cytotoxic potential of the CeO2-NRs was also investigated with the MTT assay. The CeO2-NRs were found to effectively kill off the cancerous cells (MCF-7 human breast cancer cell line), further indicating that the synthesized NRs exhibit anticancer potential as well. One of the major applications studied for the prepared CeO2-NRs was performing the statistical optimization of the photocatalytic degradation reaction of the methyl orange (MO) dye. The reaction was optimized by using the technique of response surface methodology (RSM). This advanced approach facilitates the development of the predictive model on the basis of central composite design (CCD) for this degradation reaction. The maximum degradation of 99.31% was achieved at the experimental optimized conditions, which corresponded rather well with the predicted percentage degradation values of 99.58%. These results indicate that the developed predictive model can effectively explain the performed experimental reaction. To conclude, the CeO2-NRs exhibited excellent results for multiple applications.

The production of reactive oxidative groups as byproducts of (normal/pathological) human metabolic activities is regarded as a naturally occurring phenomenon. Recent research has linked the overproduction of these reactive species to a variety of serious health problems, including cardiovascular diseases, weakened immune systems, neurological diseases, the onset of cancer, and other degenerative diseases [16]. Every living organism is equipped with some type of defense mechanism to counteract these potential oxidative stresses. [17]. Similarly, synthetic antioxidants such as butylated hydroxyanisole (BHA) and butylated hydroxytoluene (BHT) are used as food additives to prevent oxidative decomposition. However, nutritionists are somewhat concerned about the safety issues surrounding the use of these compounds, and the focus has shifted toward identifying natural anti-oxidants and nutraceuticals to improve the nutritional quality of consumers [16]. Consequently, the development of new antioxidant materials that can be used in place of these harmful synthetic antioxidants is of the utmost importance.
Another problem that has attracted the interest of the scientific community in recent years is the presence of noxious pollutants in aqueous environments, as these pollutants ultimately end up directly/indirectly impacting human beings. Numerous removal methodologies (including adsorption [16], catalytic reduction, Fenton oxidation, flocculation, coagulation, etc.) have been utilized to remove these pollutants from water reservoirs. However, these methodologies suffer from the drawbacks of incomplete removal, economically unviable, requirement of additional setups, and extensive recovery procedures [17]. Here, we have utilized the process of photocatalysis for the removal of the noxious pollutant of methyl orange (MO) by using the synthesized CeO 2 -NRs. The complete removal/mineralization of the pollutant into less harmful products of water (H 2 O) and carbon dioxide (CO 2 ), with economical time consumption, and high removal rate are some of the major advantages of photocatalysis over other conventional techniques [18]. The NRs were prepared by using the phytochemicals present in the D. sissoo extract as the stabilization and fabrication medium. The synthesized CeO 2 -NRs were characterized by using different analytical techniques. Furthermore, the antioxidant, antibacterial, and photocatalytic potential of the synthesized CeO 2 was also investigated in this study.

Preparation of D. sissoo-Stabilized CeO 2 -NRs
For the preparation of the D. sissoo bioextract, fresh leaves of D. sissoo were collected from the botanical garden of the University of Gujrat. The collected leaves were properly cleaned, washed, dried, and ground into fine powder for further experimentation. A total of 10 g of the powdered D. sissoo leaves was added to 150 mL of the double-distillated water and was continuously stirred for 24 h. The acquired extract was filtered and the filtrate was stored in air-tight sample bottles at 4 • C. For the preparation of the CeO 2 -NRs, 10 mL of 3 mM of Ce(NO 3 ) 3 .6H 2 O was added to 100 mL of D. sissoo extract. The reaction mixture was continuously stirred for 6 h at the temperature of 75 ± 3 • C until white precipitates were obtained. The formed precipitates were separated by using centrifugation and washed with water and ethanol several times to remove any impurity present in the sample. The precipitates were further placed in an oven for 2 h at a temperature of 150 • C. The acquired CeO 2 -NRs were stored in air-tight sample bottles for further experimentation [19].

Characterization
The ultraviolet-visible (UV-VIS) spectroscopy of the synthesized CeO2-NRs was carried out by using double-beam PerkinElmer, Incorporation®(LAMBDA-365+; Waltham Massachusetts, USA) spectrophotometer. The Bruker Corporation (Aubrey, Texas, USA) diffractometer was utilized for performing the X-ray diffraction (XRD) analysis. The LYRA-3 (TESCAN; Edgerton, Missouri, USA) microscope was used for identifying the morphological characteristics of the synthesized CeO2-NRs via scanning electron microscopy (SEM). The Fourier transform infrared spectroscopy (FTIR, Carry-630 Agilent; Santa Clara, California, USA), in the range of 4000-400 cm−1, was also performed for identifying the functional group of the phytochemicals involved in the stabilization of the CeO 2 -NRs.

DPPH Radical Scavenging Assay
The DPPH radical scavenging assay was used for investigating the antioxidant potential of CeO 2 -NRs. In brief, the solution containing the 23 mg/mL solution of DPPH in ethanol was prepared and subsequent absorbance of the DPPH radical ion was measured at 517 nm. The BHA was used as the positive control for the test. All the tests were replicated thrice to ensure the precision of the acquired results. The prevention capacity of the active radical was measured by using the standard of ascorbic acid [20].

Biochemical Assays for Human Lymphocytes Incubated with CeO 2 -NRs
The human blood sample (2 mL), from a healthy donor, was collected and centrifuged at 3000 rpm for 15 min. The acquired white-ish layer (that was formed just below the upper plasma layer) was removed and collected in a separate sample bottle. The collected lymphocytes were washed thrice with ammonium chloride to remove any impurity (particularly red blood cells) and the material was suspended in PBS to be stored at 4 • C [21]. The acquired lymphocytes were incubated with the synthesized CeO 2 -NRs in varying concentrations of 15, 30, 60, and 120 µmol/mL. The treated lymphocytes were tested for numerous oxidative stress biomarkers.
The total antioxidant assay (TAC) for the CeO 2 -NRs was investigated by using the ferric reducing ability of plasma (FRAP) methodology. This methodology determines the capability of the plasma (incubated with CeO 2 -NRs) to reduce Fe 3+ ions into Fe 2+ ions. The reduced Fe 2+ ions exhibited a blue color complex with the TPTZ whose amount can be spectroscopically measured owing to the maximum absorbance of this complex observed to be at 593 nm [22].
The lipid peroxidation (LPO) methodology was also utilized to determine the antioxidant capacity of the lymphocytes incubated with CeO 2 -NRs. The extent of MDA produced during the acid heating reaction serves as the main species utilized for spectroscopically measuring the antioxidant activity. In a typical assay, 10 µL of lymphocyte samples containing CeO 2 -NRs in varying concentrations or a standard MAD solution and 40 µL of PBS were added to an Eppendorf tube placed in an ice bath. In each sample tube, TBA reagent (containing 100 µL of TBA, 30 µL of phosphotungstic acid, 50 µL of sodium dodecyl sulphate, and 200 µL of HCl) was added and these air-tight tubes were boiled for 25 min at 100 • C in the water bath. Next, 400 µL of n-butanol was added in the sample and the tubes were centrifuged at 3000 rpm for 10 min. The separated supernatants were added in a 96-well plate. The excitation/emission fluorescence wavelengths of 515/555 nm were measured using the microplate reader [23].
The catalase activity (CAT) was also assessed by observing the decrease in the value of the absorption intensity at 240 nm associated with H 2 O 2 owing to the presence of the lymphocyte medium containing the CeO 2 -NRs. The specific CAT activity was measured in units/mL with one unit of enzyme equal to the enzyme required for consuming1 mol of H 2 O 2 per minute [24].

Cytotoxicity of CeO 2 -NRs
The cytotoxicity potential of the CeO 2 -NRs was investigated with the MTT assay. This quantitative colorimetric evaluation was used to observe the activity of the enzymes. The MCF-7 (human breast cancer cell line) was seeded in a 96-well plate (with the density of 10,000 cells/well) and these plates were kept at 37 • C for 24 h. The CeO 2 -NRs (in varying concentrations of 0, 25, 50, 75, 100, and 250 µg/mL) were inoculated into the grown cell with 100 µL of the medium. After every 24 h of incubation, 5 mg/mL (20 µL) MTT dissolved in the buffer was added into every well. After incubation, the formazan crystals were secured while the remaining media were discarded. The crystals (shaped via MTT metabolism) were liquefied and dissolved into 100 µL of DMSO. The optical absorbance of the acquired plates was measured at 590 nm [25].

Photocatalytic Efficacy of CeO 2 -NRs
The photocatalytic potential of the CeO 2 -NRs was calculated by carrying out the model degradation reaction of MO. All the photocatalytic reactions were carried out under a 10 W ultraviolet (UV) lamp. The response surface methodology (RSM) was utilized as a means to develop the design of the experiment for this photocatalytic reaction. The central composite design (CCD) was developed by using the three variables of MO dose, photocatalyst dose, and pH while the parameter of the percentage degradation (D%) was used as the main response of the CCD. Twenty experiments were performed and the details of these experiments are provided in the Supplementary File. In a typical photocatalytic reaction, the specific amount of MO and photocatalysts (as specified by the design of the experiment CCD model) was added in a photocatalyst reactor. The reaction was left to be homogenized for 30 min to ensure the presence of adsorption/de-sorption equilibrium. The reaction mixture was then placed under a UV lamp for the progression of the reaction. A total of 3 mL of the reaction mixture was collected after regular intervals for measuring the absorbance (at a wavelength of 463 nm) for calculating the main response of D% as given by Equation 1.
C o and C res represent the concentration of MO at the start of the photocatalytic reaction and at any time (t), respectively.

Synthesis
The phytochemicals present in the D. sissoo extract were utilized as fabricating as well as stabilizing agents for the production of CeO 2 -NRs. The physiochemical profiling detailing the identification of the phytochemicals involved has already been reported by our research group [14]. The extract was found to essentially contain numerous natural oxidants and reductants including tannis, reducing sugars, terpenoids, glycosides, proteins, and saponins [14]. In another work, we identified the following specific biomolecules including sissoic acid, biochanin A, quercetin-3-O-rutinoside, kaempferol-3-O-rutinoside, and dalsissoside which were present in abundance in the D. sissoo extract. The most significant biomolecule was found to be quercetin, which was observed in the UV-VIS spectrum of the D. sissoo extract [12,15]. The cumulative interaction of the phytochemicals associated with the precursor salt of the cerium was found to be responsible for the preparation of CeO 2 -NRs as indicated in Figure 1.
nificant biomolecule was found to be quercetin, which was observed in the UV-VIS spectrum of the D. sissoo extract [12,15]. The cumulative interaction of the phytochemicals associated with the precursor salt of the cerium was found to be responsible for the preparation of CeO2-NRs as indicated in Figure 1.

Characterization
The successful fabrication of CeO 2 -NRs was primarily validated by using the XRD technique, as shown in Figure 2. The XRD spectrum indicated that the acquired sample possessed the cubic crystal lattice system with the observed lattice parameters of a = b = c = 5.4113Å. The indexing of the XRD spectrum revealed it to have diffraction peaks at the 2θ values of 28.55 • , 33.08 • , 47.47. • , 56.33 • , 59.08 • , 69.40 • , 76.70 • , and 79.07 • associated with the diffraction planes of (111), (200), (220), (311), (222), (400), (331), and (420), respectively. The XRD pattern correlated well with the standard CeO 2 -NRs JPCDS Card no. 00-034-0394. It was also found that the CeO 2 -NRs sample was extremely pure as no unindexed peak was observed in the case of the CeO 2 -NRs' XRD spectrum [26]. The UV-VIS spectrum and the Tauc plot of the CeO2-NRs is presented in Figure 3. As indicated in Figure 3A, the CeO2-NRs do not exhibit the optical absorbance band in the visible region. The characteristics band at 285 nm of the CeO2-NRs was observed in the The UV-VIS spectrum and the Tauc plot of the CeO 2 -NRs is presented in Figure 3. As indicated in Figure 3A, the CeO 2 -NRs do not exhibit the optical absorbance band in the visible region. The characteristics band at 285 nm of the CeO 2 -NRs was observed in the ultraviolet (UV) region, indicating that the UV lamp should be preferred for the activation of the CeO 2 -NRs for photocatalytic applications [25]. The Tauc equation (as presented in Equation 2) was also utilized to determine the band gap (E g ) values of the synthesized NRs.
where hv represents the energy of the incident photon, h exhibits the Planck's constant, α is the optical coefficient, v represents the frequency of the photon, and A represents the transition constant. The direct E g value of the CeO 2 -NRs was found to be 3.4541 eV, which is again indicative of the fact that the engineered NMs are capable of utilizing the UV region effectively [27]. The UV-VIS spectrum and the Tauc plot of the CeO2-NRs is presented in Figure 3. As indicated in Figure 3A, the CeO2-NRs do not exhibit the optical absorbance band in the visible region. The characteristics band at 285 nm of the CeO2-NRs was observed in the ultraviolet (UV) region, indicating that the UV lamp should be preferred for the activation of the CeO2-NRs for photocatalytic applications [25]. The Tauc equation (as presented in Equation 2) was also utilized to determine the band gap (Eg) values of the synthesized NRs.
where hv represents the energy of the incident photon, h exhibits the Planck's constant, is the optical coefficient, v represents the frequency of the photon, and A represents the transition constant. The direct Eg value of the CeO2-NRs was found to be 3.4541 eV, which is again indicative of the fact that the engineered NMs are capable of utilizing the UV region effectively [27].  The morphological characteristics of the CeO 2 -NRs were studied by using SEM as presented in Figure 4. The SEM micrograph confirms that the engineered material exhibited rod-like morphology. However, the size of the synthesized CeO 2 -NRs was not uniform, which is in accordance with the preparation methodology as the green phytochemicalbased synthetic methodology is a crude technique for the preparation of NMs. The rough surface observed in the case of the CeO 2 -NRs further suggests that the synthesized NRs will be effective for surface-based applications. The morphological characteristics of the CeO2-NRs were studied by using SEM as presented in Figure 4. The SEM micrograph confirms that the engineered material exhibited rod-like morphology. However, the size of the synthesized CeO2-NRs was not uniform, which is in accordance with the preparation methodology as the green phytochemical-based synthetic methodology is a crude technique for the preparation of NMs. The rough surface observed in the case of the CeO2-NRs further suggests that the synthesized NRs will be effective for surface-based applications. The FTIR spectrum of the D. sissoo extract and the CeO2-NRs is elaborated in Figure  5. The presence of all the major characteristic peaks of the D. sissoo extract in the FTIR spectrum of the CeO2-NRs reveals that the phytochemicals of the D. sissoo extract act as a stabilization medium for the fabricated NRs. The presence of a strong band at 3350 cm −1 represents the vibrational frequency of the hydroxyl group, indicating that the D. sissoo extract is predominantly composed of polyphenols. This band was reduced in intensity in The FTIR spectrum of the D. sissoo extract and the CeO 2 -NRs is elaborated in Figure 5. The presence of all the major characteristic peaks of the D. sissoo extract in the FTIR spectrum of the CeO 2 -NRs reveals that the phytochemicals of the D. sissoo extract act as a stabilization medium for the fabricated NRs. The presence of a strong band at 3350 cm −1 represents the vibrational frequency of the hydroxyl group, indicating that the D. sissoo extract is predominantly composed of polyphenols. This band was reduced in intensity in the case of the CeO 2 -NRs, which is indicative of the fact that these groups strongly interacted with the CeO 2 -NRs, resulting in the displacement of this peak from 3350 cm −1 to 3368 cm −1 [28]. The decrease in the intensity can be further associated with the process of oven drying, as high temperatures also removed any attached water from the CeO 2 -NRs and the remaining band only indicated that the polyphenols of the phytochemicals adhered to the surface of the NRs as the stabilization medium [12]. The decrease in the intensity can be further associated with the process of oven drying as high temperatures also removed any attached water from the CeO 2 -NRs. The subsequent bands at 2920 cm −1 (CH group; stretching), 1724 cm −1 (C=O group), 1649 cm −1 (C=C group), and 1382 cm −1 (CH group; bending) were present in the D. sissoo extract. These bands were slightly displaced in the case of the CeO 2 -NRs, indicating that the phytochemicals indeed interacted with the CeO 2 -NRs [29]. The presence of the characteristics additional band at 573 cm −1 associated with the vibrational frequency of the Ce-O bond further validates the formation of CeO 2 -NRs [29].

DPPH Assay for CeO2-NRs
The antioxidant capacity of the D. sissoo extract and the CeO2-NRs was tested by using the DPPH assay as indicated in Figure 6. It was observed that the CeO2-NRs exhibited even better antioxidant results in comparison to the industrial antioxidant of BHA at all the concentrations. Moreover, the observed antioxidant activity of the CeO2-NRs was found to be dose-dependent, as with the increase in concentration from the series 1 (15 g/mL) to series 6 (250 g/mL), the percentage inhibition of the DPPH free radicals by CeO2-NRs also increases. This is to be expected as CeO2 (apart from having the antioxidant potential of the phytochemicals present on the surface) also possesses lower ratios of Ce 3+ /Ce 4+ on the surface that contribute to enhance its antioxidant potential in comparison to the extract. Dutta et al. [30] documented similar results for the CeO2-NRs. Furthermore, it should also be observed that the D. sissoo extract itself acted as a rather good antioxidant material and surpassed the antioxidant efficacy of the BHA at the higher concentrations. At lower concentrations of the extract, the BHA exhibited better results than the D. sissoo extract. This could be further interpreted as an interesting insight regarding the effectiveness of the D. sissoo extract, revealing that a concentrated extract is more effective for achieving better pharmacological applications of this traditional ethno-medicinal plant

DPPH Assay for CeO 2 -NRs
The antioxidant capacity of the D. sissoo extract and the CeO 2 -NRs was tested by using the DPPH assay as indicated in Figure 6. It was observed that the CeO 2 -NRs exhibited even better antioxidant results in comparison to the industrial antioxidant of BHA at all the concentrations. Moreover, the observed antioxidant activity of the CeO 2 -NRs was found to be dose-dependent, as with the increase in concentration from the series 1 (15 µg/mL) to series 6 (250 µg/mL), the percentage inhibition of the DPPH free radicals by CeO 2 -NRs also increases. This is to be expected as CeO 2 (apart from having the antioxidant potential of the phytochemicals present on the surface) also possesses lower ratios of Ce 3+ /Ce 4+ on the surface that contribute to enhance its antioxidant potential in comparison to the extract. Dutta et al. [30] documented similar results for the CeO 2 -NRs. Furthermore, it should also be observed that the D. sissoo extract itself acted as a rather good antioxidant material and surpassed the antioxidant efficacy of the BHA at the higher concentrations. At lower concentrations of the extract, the BHA exhibited better results than the D. sissoo extract.
This could be further interpreted as an interesting insight regarding the effectiveness of the D. sissoo extract, revealing that a concentrated extract is more effective for achieving better pharmacological applications of this traditional ethno-medicinal plant [31].  Table 1 provides the results of the biochemical assays documented for the human lymphocytes incubated with the varying concentrations of CeO2-NRs. The TAC assay revealed that the increasing concentration of CeO2-NRs in the lymphocytes enhanced the FRAP values. This is indicative of the fact that the CeO2-NRs possessed antioxidant potential and effectively facilitated the Fe 3+ ions reduction into Fe 2+ ions. The reduced Fe 2+ ions then combined with the TPTZ to provide a complex which was spectroscopically measured. The enhanced concentration of the Fe 2+ ions indicated that from the concentration of 15-60 mol/mL, an increase in antioxidant potential with the concentration was observed, while at the concentration of 120 mol/mL, the trend was reversed. This is crucial as it reveals that the lower concentrations of CeO2-NRs were found to be more effective in terms of antioxidant potential [32]. The significant decrease in the concentration of peroxide radical was observed in the case of the human lymphocytes incubated with the CeO2-NRs. The optimum concentration for this assay was also found to be 60 mol/mL, which is again indicative of the fact that the antioxidant potential of the CeO2-NRs was effective at lower concentration of ROS [24]. Increasing the dose of the CeO2-NRs in the system also resulted in an increase in the peroxide radicals scavenging ability of the NRs as well [33].  Table 1 provides the results of the biochemical assays documented for the human lymphocytes incubated with the varying concentrations of CeO 2 -NRs. The TAC assay revealed that the increasing concentration of CeO 2 -NRs in the lymphocytes enhanced the FRAP values. This is indicative of the fact that the CeO 2 -NRs possessed antioxidant potential and effectively facilitated the Fe 3+ ions reduction into Fe 2+ ions. The reduced Fe 2+ ions then combined with the TPTZ to provide a complex which was spectroscopically measured. The enhanced concentration of the Fe 2+ ions indicated that from the concentration of 15-60 µmol/mL, an increase in antioxidant potential with the concentration was observed, while at the concentration of 120 µmol/mL, the trend was reversed. This is crucial as it reveals that the lower concentrations of CeO 2 -NRs were found to be more effective in terms of antioxidant potential [32]. The significant decrease in the concentration of peroxide radical was observed in the case of the human lymphocytes incubated with the CeO 2 -NRs. The optimum concentration for this assay was also found to be 60 µmol/mL, which is again indicative of the fact that the antioxidant potential of the CeO 2 -NRs was effective at lower concentration of ROS [24]. Increasing the dose of the CeO 2 -NRs in the system also resulted in an increase in the peroxide radicals scavenging ability of the NRs as well [33].

Cytotoxicity Effects of CeO 2 -NRs on the MCF-7 Cell Line
The cytotoxic effect of CeO 2 -NRs was investigated by using the human breast cancer MCF-7 cell line as indicated in Figure 7. It was observed that the addition of the CeO 2 -NRs significantly reduced the percentage viability of the MCF-7 cell line. This cytotoxic effect is attributed to the extremely small size of the NRs as the dangling bonds available at the surface of the NRs make them chemically reactive particularly for biochemical catalytic reactions [34]. These NMs can penetrate into the cells owing to their high chemical reactivity and cause the lysis of the cell. The acquired results correspond rather well with the findings of Navada et al. [35]. The decrease in the viability of the cancerous cells was also found to be dose-as well as time-dependent. The higher doses (i.e., 100 µg/mL) and the treatment time duration of 72 h were found to be most effective with a maximum decrease of 35.124% in the viability of the MCF-7 cells.

Cytotoxicity Effects of CeO2-NRs on the MCF-7 Cell Line
The cytotoxic effect of CeO2-NRs was investigated by using the human breast cancer MCF-7 cell line as indicated in Figure 7. It was observed that the addition of the CeO2-NRs significantly reduced the percentage viability of the MCF-7 cell line. This cytotoxic effect is attributed to the extremely small size of the NRs as the dangling bonds available at the surface of the NRs make them chemically reactive particularly for biochemical catalytic reactions [34]. These NMs can penetrate into the cells owing to their high chemical reactivity and cause the lysis of the cell. The acquired results correspond rather well with the findings of Navada et al. [35]. The decrease in the viability of the cancerous cells was also found to be dose-as well as time-dependent. The higher doses (i.e., 100 g/mL) and the treatment time duration of 72 h were found to be most effective with a maximum decrease of 35.124% in the viability of the MCF-7 cells.

Assortment of Predictive Model
The three-factor CCD matrix and the acquired photocatalytic results acquired for the photocatalytic degradation reaction of MO are presented in Figure 8. The detailed design of the experiment as well as the selection tests performed to identify the right model are presented in the supplementary information. The response factor (percentage degradation of MO) was found to be a second-order polynomial (quadratic model), as presented in Equation 3. The selection of the model was achieved by using Fischer's exact test (p-value), the F-statistics test (F-value), regression square (R 2 ), adjusted R 2 , and predicted R 2 as detailed in the Supplementary File. The analysis of variance (ANOVA) test was also performed for further experimentation. The three-factor CCD matrix and the acquired photocatalytic results acquired for the photocatalytic degradation reaction of MO are presented in Figure 8. The detailed design of the experiment as well as the selection tests performed to identify the right model are presented in the supplementary information. The response factor (percentage degradation of MO) was found to be a second-order polynomial (quadratic model), as presented in Equation 3. The selection of the model was achieved by using Fischer's exact test (p-value), the F-statistics test (F-value), regression square (R 2 ), adjusted R 2 , and predicted R 2 as detailed in the Supplementary File. The analysis of variance (ANOVA) test was also performed for further experimentation.
where the factors A, B, and C represent the coded individual factors of MO dose, photocatalyst dose, and pH units, respectively. The factors of AB, AC, and BC represent the interaction terms while the remaining A 2 , B 2 , and C 2 represent the squared terms.

ANOVA
The ANOVA test was applied on the quadratic model to verify the validity of the selected model as explained in Table 2. The statistically significant p-value (i.e., <0.0001) indicated that the predictive model selected for explaining the photocatalytic degradation of MO correlated rather well with the observed experimental results. Furthermore, the individual terms of A, B, C, and the interaction term of AB were found to be statistically significant. The other interaction terms of AC and AD were not found to be statistically significant, indicating that the correlation among these factors did not impact the %D of the reaction as much. In the case of the ANOVA test, it was further observed that our predictive model provided the F-value of 149.5, which was found to be rather impactful. Conventionally, the model is considered significant when the tabulated F-value (15.11) is almost four or five times smaller in comparison to the calculated F-value [36]. Moreover, the p-value for the confidence interval of 99% and 95% (i.e., under the significance levels of 1% and 5%, respectively) is also an important way to assess the validity of the model. A p-value of less than 0.05 (for a confidence interval of 95%), or even better with a p-value of less than 0.01 (for a confidence interval of 99%) is typically considered statistically significant. The p-value for the elected model was found to be <0.0001, which indicated that

ANOVA
The ANOVA test was applied on the quadratic model to verify the validity of the selected model as explained in Table 2. The statistically significant p-value (i.e., <0.0001) indicated that the predictive model selected for explaining the photocatalytic degradation of MO correlated rather well with the observed experimental results. Furthermore, the individual terms of A, B, C, and the interaction term of AB were found to be statistically significant. The other interaction terms of AC and AD were not found to be statistically significant, indicating that the correlation among these factors did not impact the %D of the reaction as much. In the case of the ANOVA test, it was further observed that our predictive model provided the F-value of 149.5, which was found to be rather impactful. Conventionally, the model is considered significant when the tabulated F-value (15.11) is almost four or five times smaller in comparison to the calculated F-value [36]. Moreover, the p-value for the confidence interval of 99% and 95% (i.e., under the significance levels of 1% and 5%, respectively) is also an important way to assess the validity of the model. A p-value of less than 0.05 (for a confidence interval of 95%), or even better with a p-value of less than 0.01 (for a confidence interval of 99%) is typically considered statistically significant. The p-value for the elected model was found to be <0.0001, which indicated that the model is rather good and lies within the significance level of 1% [37]. Combining both the acquired model and ANOVA results revealed that the factors of A, B, C, and AB should be given special emphasis for modulating the percentage degradation of MO. The factor A (MO dose) had a negative sign in the model, indicating that any increment in the factor A will negatively impact the D%. The same results were observed in the case of factor AB. However, the value of the interaction coefficients of the individual terms A, B, and C was found to be rather higher in comparison to the interaction term of factor AB. This further indicated that the optimization of this reaction is heavily dependent on the individual factors rather than on the interaction terms [38].

Diagnostics Plots for the Selected Model
The validity of the selected quadratic model was further investigated by numerous validity plots, as indicated in Figure 9. A normal plot associated with the residuals of the photocatalytic degradation of MO is presented in Figure 8A. The slight deviation of the acquired experimental values from the straight line is expected as the errors can never be completely eradicated from the experiments [39]. The acquired response in terms of the D% values for the photocatalytic reaction of MO in the presence of the CeO 2 -NRs was also found to be rather comparable with the predicted values ( Figure 8B), indicating that the selected model was significantly valid for explaining the behavior of this particular reaction [37].

RSM-Based Three-Dimensional and Contour Plot
The interaction terms can be better explained on the basis of their contour or threedimensional surface plots. Figure 10A,B explain the contour and three-dimensional surface plot observed in the case of the interaction factor of AB. Other interaction factors (AC and BC) were found to be statistically insignificant, indicating that the variation in these factor does not affect the D% values of the reaction and will therefore not be discussed here. It is a typical convention to only discus the statistically significant interaction factors via contour or surface plots [12]. As indicated by both these graphs, the maximum degradation indicated by the red colored region was observed at the high CeO2-NRs of about 6 mg/L. This is justified as the increment in the concentration of photocatalyst ultimately increased the number of adsorption sites available for the dye molecules to get attached. The adsorbed dye molecule is then able to interact with the free radicals generated in the medium owing to the activation of the photocatalyst (generation of electron-hole pair) by UV light irradiation. These interactions will result in the degradation of the MO dye into less harmful products. As increasing the photocatalyst concentration positively impacts the D% values, higher values of catalyst should be preferred for achieving maximum degradation. Moreover, the lowest MO dose was found to be effective for achieving the maximum degradation. This is again expected as the lower the concentration of MO, the easier it will be to degrade it. There will be no diffusion barrier and the dye molecules will easily interact with the NMs and ultimately be degraded [39].

RSM-Based Three-Dimensional and Contour Plot
The interaction terms can be better explained on the basis of their contour or threedimensional surface plots. Figure 10A,B explain the contour and three-dimensional surface plot observed in the case of the interaction factor of AB. Other interaction factors (AC and BC) were found to be statistically insignificant, indicating that the variation in these factor does not affect the D% values of the reaction and will therefore not be discussed here. It is a typical convention to only discus the statistically significant interaction factors via contour or surface plots [12]. As indicated by both these graphs, the maximum degradation indicated by the red colored region was observed at the high CeO 2 -NRs of about 6 mg/L. This is justified as the increment in the concentration of photocatalyst ultimately increased the number of adsorption sites available for the dye molecules to get attached. The adsorbed dye molecule is then able to interact with the free radicals generated in the medium owing to the activation of the photocatalyst (generation of electron-hole pair) by UV light irradiation. These interactions will result in the degradation of the MO dye into less harmful products. As increasing the photocatalyst concentration positively impacts the D% values, higher values of catalyst should be preferred for achieving maximum degradation. Moreover, the lowest MO dose was found to be effective for achieving the maximum degradation. This is again expected as the lower the concentration of MO, the easier it will be to degrade it. There will be no diffusion barrier and the dye molecules will easily interact with the NMs and ultimately be degraded [39].

Optimization
The RSM-based optimization approach revealed that the maximum degradation value of 99.58% should be achieved with the optimized reaction parameter values of MO dose = 20.1371 ppm; photocatalyst dose = 5.88582 mg/L, and pH units of 8.9729. The results documented by Kaur et al. [40] (optimum pH of 9.00 units), Sithole et al. [41] (optimum MO dose of 20 ppm), and Khalilian et al. [42] (optimum photocatalyst dose of 7 mg/L) relate strongly with the acquired results in our experiment. In order to further validate the optimized parameters provided by the RSM, the typical photocatalytic reaction operating under the specified optimization parameters was carried out. A D% value of 99.317% was observed by using this point prediction analysis, indicating that the model was successful in explaining this model. The ramp graph showing the predictive results at the optimized conditions is exhibited in Figure 11.

Optimization
The RSM-based optimization approach revealed that the maximum degradation value of 99.58% should be achieved with the optimized reaction parameter values of MO dose = 20.1371 ppm; photocatalyst dose = 5.88582 mg/L, and pH units of 8.9729. The results documented by Kaur et al. [40] (optimum pH of 9.00 units), Sithole et al. [41] (optimum MO dose of 20 ppm), and Khalilian et al. [42] (optimum photocatalyst dose of 7 mg/L) relate strongly with the acquired results in our experiment. In order to further validate the optimized parameters provided by the RSM, the typical photocatalytic reaction operating under the specified optimization parameters was carried out. A D% value of 99.317% was observed by using this point prediction analysis, indicating that the model was successful in explaining this model. The ramp graph showing the predictive results at the optimized conditions is exhibited in Figure 11.

Optimization
The RSM-based optimization approach revealed that the maximum degrad value of 99.58% should be achieved with the optimized reaction parameter values o dose = 20.1371 ppm; photocatalyst dose = 5.88582 mg/L, and pH units of 8.9729. The r documented by Kaur et al. [40] (optimum pH of 9.00 units), Sithole et al. [41] (opti MO dose of 20 ppm), and Khalilian et al. [42] (optimum photocatalyst dose of 7 m relate strongly with the acquired results in our experiment. In order to further va the optimized parameters provided by the RSM, the typical photocatalytic reaction ating under the specified optimization parameters was carried out. A D% value of 99. was observed by using this point prediction analysis, indicating that the model wa cessful in explaining this model. The ramp graph showing the predictive results a optimized conditions is exhibited in Figure 11.

Comparison with the Literature
The acquired results regarding the photocatalytic degradation of MO performed in the presence of the CeO 2 -NRs are presented in Table 3.

Conclusions
In this study, the phytochemicals present in the D. sissoo extract were utilized to engineer CeO 2 -NRs. Not only was this technique found to be effective for the preparation of morphology-specific NMs (i.e., rod-like CeO 2 -NRs), but the NMs were also prepared at rather lower temperatures as opposed to other methodologies. The acquired CeO 2 -NRs were characterized by using XRD, UV-VIS, FTIR, and SEM techniques. All these techniques validated the formation of the NRs by using this phytochemistry-based approach. The CeO 2 -NRs exhibited excellent antioxidant, cytotoxic, and photocatalytic results, indicating that the engineered NRs can be effectively utilized to perform a multitude of practical applications. The best results were found in the case of the antioxidant activity, as the CeO 2 -NRs expressed better antioxidant activity in comparison to industrial antioxidant of BHA. The RSM-based photocatalytic degradation of MO carried out by using the CeO 2 -NRs exhibited the percentage degradation value of 99.317%, which is rather higher in comparison to the degradation percentages documented in the academic literature for the CeO 2 -NRs. Our study signifies that the synthesized CeO 2 -NRs exhibit copacetic potential for multiple applications.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules27238188/s1, Table S1. Reaction variables utilized for developing 3-factor based Central Composite Design (CCD) for the photocatalytic degradation of methyl orange (MO) in the presence of CeO2-NRs. Table S2. CCD based design of experiment for the photocatalytic degradation reaction of MO. Table S3. Statistical analysis for investigating the best predictive model for the photocatalytic degradation of MO performed in the presence of CeO2-NRs. Data Availability Statement: The data will be made available upon request.