Biogenic engineered zinc oxide nanoparticle for sulfur black dye removal from contaminated wastewater: comparative optimization, simulation modeling, and isotherms

ABSTRACT This research work aimed to isolate and culture the bacterium Bacillus paramycoides for biogenic fabrication of zinc oxide nanoparticles, specifically ZnO and ZnO-ME nanoparticles (nanoparticles fabricated from bacterial extracts only – ZnO, and from bacterial cell mass including extract – ZnO-ME). SEM investigation revealed the spherical-shaped NPs with 22.33 and 39 nm in size for ZnO and ZnO-ME, respectively. The Brunauer, Emmett, and Teller (BET) studies revealed mesoporous structure with pore diameters of 13.839 and 13.88 nm and surface area of 7.617 and 33.635 m2/gm for ZnO and ZnO-ME, respectively. Various parameters for the adsorption of sulfur black dye onto both ZnO and ZnO-ME were screened and optimized using Plackett–Burman Design (PBD), Full Factorial Design (FFD) and Central Composite Design (CCD). The results of the optimization modeling study revealed that FFD yielded the most predictable and best-fitting results among all the models studied, with R2 values of 0.998 for ZnO and 0.993 for ZnO-ME. Notably, ZnO-ME exhibited a greater dye removal efficiency 80% than ZnO i.e., 71%, it may be due to the presence of amorphous carbon on the surface of ZnO-ME. Among the various isothermal models, the Freundlich model displayed the strongest correlation with the dye removal data, confirming the multilayer adsorption of dye on both nanoparticles and supporting physisorption. Therefore, ZnO and ZnO-ME nanoparticles have been proven as potential tools for mitigating environmental impacts associated with dye-containing wastewater.


Introduction
Every government must ensure that its citizens have the access to clean, safe drinking water, yet as economies develop, many nations experience a depletion of their natural resources, especially freshwater [1,2].Nowadays, water resources are continually contaminated by a diverse array of pollutants, for instance, dyes, micropollutants, pesticides, pharmaceuticals and heavy metals, creating a significant global environmental problem [3].This contamination is particularly prevalent in aquatic environments due to the widespread use of synthetic dyes in diverse industrial applications like plastics, paper, antiseptics, cosmetics, food, textiles, and leather [4].Among these dyes, sulfur dyes stand out as one of the most commonly used in the textile industry because of their affordability and capacity to produce darker, more subdued hues of color, including green, dark blue, black, brown, and olive [5].
However, the exposure and transport of wastewater containing these dyes from textile production can pose significant risks to aquatic life and plant ecosystems.Acute exposure to this toxic dye can result in various health issues in humans, including higher heart rate, nausea, shock, cyanosis, jaundice, quadriplegia, and tissue necrosis.Therefore, before discharging industrial effluents into the environment, it is imperative to remove color from effluent.Numerous methods such as solvent extraction, adsorption, membrane filtration, photocatalytic degradation using metal oxide, combined electrochemical degradation process, ozonation, advanced oxidation process, and ultrasonication have been adopted to treat industrial wastewater [5,6].Adsorption is one of these methods that has attracted particular attention because of its technological viability, versatility, and simplicity of use.Ion-exchange materials, bentonite, zeolites, and activated carbon clays are commonly used as adsorbents [7].Their broad range of uses can be attributed to their ordered pore size, vast surface area, excellent stability, ease of production, simplicity, rapidity, specificity, sensitivity, and cost-effectiveness [8,9] Furthermore, recent research has focused toward the fabrication of nanomaterials for wastewater treatment [10].The implementation of metal oxide nano-adsorbents in the adsorption process has gained attraction within the scientific community, representing the next generation of water treatment nanotechnology [12].However, traditional approaches to producing these adsorbents, such as pyrolysis and chemical synthesis, have some limitations, including high preparation costs, hazardous chemicals, and high energy demands.Therefore, the biogenic synthesis of nanomaterials has been widely studied as it offers a more environmentally friendly and economical approach.Highly effective bioreducing, biocapping, and biostabilizing agents that are derived from biological sources are essential for producing high-quality yields of nanomaterials [11].For instance, the biosynthesis of nanoparticle using Gelidium amasii (marine red alga) has exhibited antibacterial agents especially against prominent microfouling bacteria [12].While a recent study used Purpureocillium lilacinum (filamentous fungus) for the biosynthesis of CuO NPs exhibited photocatalysis of navy blue and safranin dye of 57.5% and 63% removal, respectively [13].
This research work focuses on the bacterial synthesis, characterization, and application of ZnO (nanoparticles fabricated from bacterial extracts) and ZnO-ME (bacterial cell mass) nanoparticles for the adsorption of sulfur black dye from synthetic solutions.Bacillus paramycoides bacteria was isolated and cultivated to produce these nanoparticles, which were subsequently analyzed using techniques such as zeta potential, X-ray diffraction (XRD), Fourier-transform infrared spectroscopy (FTIR), UV-visible spectrophotometry (UV-VIS), field emission scanning electron microscopy (FE-SEM), energy-dispersive X-ray spectroscopy (EDX), and Brunauer-Emmett-Teller (BET) analysis.The optimization of various factors affecting the adsorption efficiency of the synthesized nanoparticles was performed using Plackett-Burman Design (PBD), Full Factorial Design (FFD), and Central Composite Design (CCD).Lastly, isotherm models were employed to gain insights into the adsorption process.

Chemicals
Sulfur black dye, raw wastewater samples before treatment, and soil and sludge samples were obtained from RSWM Limited (a unit of LNJ denim), the textile industry located in Banswara, Rajasthan, India, with permission obtained from the landowner and local authorities.The wastewater was collected into an airtight bottle and placed in an ice-packed container.It was transported into the laboratory and transferred into the refrigerator at 5°C.The chemicals employed in this study, zinc sulfate heptahydrate (99%, CAS No.: 7446-20-0), absolute ethanol (>99%, CAS No.: 64-17-5), sodium hydroxide (98%, CAS No.: 1310-73-2), and hydrochloric acid (37%, CAS No.: 7647-01-0) were analytical grades and used as such supplied by Merck (Sigma-Aldrich, Chemicals Pvt. Ltd., Bangalore, India).Analytical-grade nutrient broth (SKU: M002) and agar (SKU: PCT0901) were used to perform the experiments and were purchased from Himedia Laboratories Pvt. Ltd., Mumbai, India.

Bacterial isolation and screening
Differential media was used for the isolation of bacterial species from the samples of wastewater, soil, and sludge.Serial dilution was performed for all the collected samples, with dilutions from 10 −1 to 10 −10 .The prepared dilution and undiluted samples were then spread onto nutrient agar plates supplemented with one mM metal salt (zinc sulfate) and 100 mg/L sulfur black dye.The experiments were performed in triplicates.The isolation of bacterial strain was achieved through repeated streaking process, and the bacterial strain was subsequently stored for further investigation and analysis.

Molecular characterisation and phylogenetic analysis
Molecular characterization of isolated bacterial strain was performed by 16S rRNA technique.Using the primers fD1 (5'-AGAGTTTGA TCCTGGCTCAG-3") and rD1 (5"-AAGGAGGT GATCCAGCC-3'), the 16S rRNA gene was amplified [14].The Microbial Culture Collection at the National Center for Cell Science in Pune, Maharashtra, India, sequenced the amplified 16S rRNA gene product on a commercial scale.The received sequence was aligned using MEGA 11 with homologous sequences discovered through NCBI's BLASTn similarity searching.Using the neighbor-joining method (NJ method) of the MEGA 11 software, a phylogenetic tree was constructed to accurately determine the phylogenetic identity of the strain.

Biosynthesis and characterization of ZnO and ZnO-ME nanoparticles
The identified bacteria were cultivated at 28 ± 1°C and 100 rpm for 24 h in a nutrient broth medium supplemented with 1 mM zinc sulfate and 100 mg/l of sulfur black dye.For the synthesis of ZnO-ME (bacterial cell mass including extract), 5 mM zinc sulfate was added to the bacterial culture after 24 h of incubation.The mixture was then incubated under the same conditions for an additional 48 h.On the other hand, the bacterial culture supernatant was extracted for the synthesis of ZnO (nanoparticles fabricated from bacterial extracts) by centrifugation at 6000 rpm for 5 min.The supernatant was then mixed with 5 mM zinc sulfate, and the mixture was incubated for an additional 48 h under the same conditions.The observation of a change in the color of the solution served to confirm the biological synthesis of ZnO and ZnO-ME nanoparticles, as shown in Figure 1.Organic/biomolecules are used as ligands (reducing and capping) in the reduction of metal ions to create metal NPs [15].Following synthesis, the nanomaterials were separated from both reaction mixtures by centrifuging them at 6-000 rpm for 5 min, with and without bacterial cell mass.The obtained nanomaterials were subjected to an oven at 150°C for 72 h and then using a pestle and mortar, crushed into a fine powder.The surface charge, wettability, and geometry of nanochannels can be precisely regulated in response to external stimulus, such as pH, temperature [16].The final ZnO and ZnO-ME nanoparticles were stored in plastic bottles in a desiccator for use in subsequent dye removal studies.
The physio-chemical characteristics of the nanoparticles were thoroughly investigated through zeta-potential, XRD, FTIR, UV-VIS, FE-SEM, EDX, and BET analysis.In order to analyze the absorption spectra of ZnO and ZnO-ME, a UV-VIS spectrophotometer was utilized.The colloidal solution of the synthesized NPs was used to measure the zeta potential using the Zetasizer (Malvern Instruments) and dynamic light scattering methods.FE-SEM imaging was used to determine nanoparticle surface morphology and size, while ESEM and EDX examination involved casting gold-coated solutions onto glass slides.In order to identify surface functional groups, FTIR analysis was carried out in the range 400-4000 cm −1 and for the analysis, a dried powder sample was mixed with KBr.XRD analysis was used to determine the nanoparticles' crystalline structure, which involved scanning in 2 theta degrees ranging from 20° to 80°.In the XRD analysis, the dried, powdered material was employed.CuKα radiation was used to study the nanoparticles, using 40 kV voltage and 40 mA of current.The specific surface area values were calculated using Autosorb IQ 2 Quantachrome equipment and the BET (Brunauer-Emmett-Teller) mathematical method employing cycles of N 2 adsorption and desorption.

Adsorption experiments
At room temperature, experiments were conducted in 100 mL conical flasks in an incubator shaker (Scigenics Biotech, Orbitek) with a 160rpm agitation speed.Sulfur black dye due to presence of negative charge on its surface precipitates out at higher pH that's why in the current study solution's pH was held constant at 7.0.Lim et al. (2021) showed the discolouration and precipitation of methyl blue and methyl orange dye at higher pH value [17].Room temperature is used in this study to make the adsorption process more cost-effective and reliable.The contact period ranged from 60 to 180 min, while the initial sulfur black dye concentration ranged from 10 to 100 mg/L in this research.Similarly, the amount of adsorbent used varied from 5 to 25 mg for each 20 mL of synthetic sulfur black dye solution.To assess the efficiency of dye's adsorption onto the ZnO and ZnO-ME nanomaterials, the initial and final concentrations of sulfur black dye were measured using a UV-VIS spectrophotometer.

Experimental design
This study examined the impact of various independent variables on the removal of sulfur black dye from simulated water, including initial dye concentration (10-100 mg/L), contact time (60-80 min), and ZnO or ZnO-ME dosage (5-25 mg).Minitab software version 21 was used for PBD and FFD statistical analysis, and Research Design Expert version 13 was utilized for CCD statistical analysis.To determine how well these models fit the experimental data, the variance analysis (ANOVA) was used.The Fischer (F-test) test with a 95% confidence interval was used to assess the models' significance as well as the interactions between the variables and the related responses.

Central Composite Design (CCD)
The CCD analyzes the relationships between process variables and their corresponding responses.It is useful for estimating secondorder polynomial equations over a larger portion of the design space.In CCD design, the fivelevel fractional factorial consisted of two axial designs, two central designs, and a two-level factorial design is included.To improve the measurement of reproductive capacity and the lack of fit in the models, the central point is repeated.Additionally, CCD designs feature orthogonality and rotatability characteristics, which facilitate curvature representation in 3D charts [18].(Equation 1) of CCD model is as follows: In this model equation, x i and x j are coded values of independent variables, a model constant (β 0 ), the expected response (Y), the number of independent variables (k), and regression coefficients for linear, quadratic, and interaction terms (β i , β ii , β ij ).

Plackett-Burman Design (PBD)
PBD was employed to determine the significance of the aspects under investigation.The design had various combinations of the independent elements in each run (A-C).Three levels of these parameters were examined: low (−), high (+), and center (0).The aim of this multivariate regression analysis was to identify the crucial variables in influencing the experimental response.PBD was established using the linear model expressed as (Equation 2): where the regression coefficients for intercept and linear expression are βo and βi, respectively; X i represents the coded independent factors and Y represents the experimental response [19].

Full Factorial Design (FFD)
Multiple regression analysis was used to determine the relationship between the response variable and the three factors.The model equation is expressed as follows: where initial sulfur black dye concentration, contact time, and ZnO or ZnO-ME amount are denoted as A, B, and C, respectively, with a response of R, global mean X 0 , and remaining regression coefficients X i [20].The full factorial design with two or more levels is considered as most effective design that may be used to screen numerous parameters, identify the critical components, and determine their suitable levels.

Determination of removal efficiency and adsorption capacity
Using UV-VIS spectroscopy at a wavelength of 600 nm, the initial and residual amounts of sulfur black dye were assessed.The standard curve of absorbance against a stock solution of sulfur black dye with concentrations ranging from 10 to 100 mg/L was used to measure the sulfur black dye concentration (refer to Supplementary Information (SI)).The sulfur black dye's adsorption capacity (mg/g) and removal effectiveness (percent) were expressed as follows: where m (g) is the mass of the ZnO and ZnO-ME, V (L) is the working solution volume used in the experiments, C o (mg/L) and C e (mg/L) is the initial concentration, and the final concentration in the solution, respectively [21].

Adsorption isotherm
Adsorption isotherms more accurately represent the equilibrium behavior by graphing the amounts of adsorbent adsorbed per unit amount of adsorbent against the corresponding equilibrium concentration in the solution phase.The equations of the Langmuir (Equation 6), Freundlich (Equation 7) and Temkin (Equation 8) isotherm adsorption models are given below: The values of Q 0 and b were determined using the slope and intercept of the Langmuir plot of C e versus C e /q e , and a dimensionless constant called separation factors (R L ) was used to test the favorability of the process [22].
The K f and 1/n values were determined by plotting the graph between log q e and log C e .High Where q e is the quantity of adsorbate sorbed per unit amount of sorbent (mg/g); C e is the concentration of sorbate at equilibrium (mg/L); R is the gas constant (0.0083 kJ/mol K); T is the temperature (K); b is the Temkin constant related to heat of sorption (kJ/mol); a is the Temkin isotherm constant called equilibrium binding constant (L/ g) corresponding to the maximum binding energy [24].A plot of q e versus log C e enables the determination of the isotherm constants a and b.

Result and discussion
This study investigates the synthesis, characterization, and application of ZnO and ZnO-ME nanoparticles from Bacillus paramycoides bacteria for the adsorption of sulfur black dye from synthetic solutions.Techniques including zeta potential, XRD, FTIR, UV-VIS, FE-SEM, EDX, and BET analysis were used.The adsorption efficiency was optimized using Plackett-Burman Design, Full Factorial Design, and Central Composite Design, and isotherm models were used for further understanding.

Molecular characterisation and phylogenetic analysis
Bacterial isolates were extracted using differential media from the textile wastewater, sludge, and soil samples.Using the NCBI's BLAST program, the isolated strain's 16S rRNA gene sequence was compared to previously reported 16S rRNA sequences of bacteria and archaea and revealed similarities with numerous Bacillus species, as shown in Figure 2. The isolated bacterial strain from the wastewater sample was identified as Bacillus paramycoides.
Furthermore, the 16S rDNA gene sequence is linked to the gene bank (https://www.ncbi.nlm.nih.gov).The accession number for strains was obtained as SAMN341003344.
ZnO nanoparticles are moderately stable, whereas ZnO-ME nanoparticles are reasonably stable, according to zeta potential measurements.These zinc nanoparticles are coated with negatively charged biomolecules, as indicated by the negative value [26].

FTIR analysis
The FTIR spectra of the nanoparticles are depicted in Figure 3a.In FTIR spectra of the synthesized zinc oxide nanoparticles (ZnO), a strong peak was observed at 484.04 and 619.72 cm −1 , suggesting the biogenic synthesis of zinc oxide nanoparticles.Metal oxides, including zinc oxide, have vibrational peaks between 400 and 600 cm −1 [27].The presence of water molecules can be inferred from the broadband that occurred between 3416 and 3458 cm −1 since the O-H stretch is associated with this peak.The -C=C stretching vibration of alkynes is responsible for the two prominent peaks at 2336 cm −1 .The C-H and HO-C=O groups may bend symmetrically and asymmetrically to produce the spectral peak at 2926.26 cm −1 [28].The N-H stretching vibration of amine groups on the sample's surface is suggested by high absorption peaks at 1633.39 cm −1 .Absorption peak at 1246 cm −1 due to C-N stretching of anime group.Inplane C-H bending is responsible for the absorption peaks at 1042 and 1406.09cm −1 [29].Similar peaks were observed in ZnO-ME; however, they were less intense compared to ZnO.

UV-VIS analysis
UV-VIS spectroscopy was used to describe the structural characteristics of ZnO and ZnO-ME nanomaterials.In the UV region, absorption peaks were seen for the ZnO.Both nanomaterials exhibited an average absorption wavelength of 364 nm, as shown in  2019), the absorption bands of zinc oxide nanoparticles typically fall between 330 and 370 nm [28,33].According to Vaseem et al. (2012), the quantum confinement effect of ZnO nanoparticles led to a blue-shifted near band-edge UV absorption at 340 nm [34].
When compared to the bulk absorption edge, which at ambient temperature appears at 373 nm, a larger excitonic absorption pattern at about 345 nm is blue-shifted by roughly 30 nm [35].The absorption peak underwent a blue shift and was found at about 229 nm which accordance with the smaller nanoparticle size [36].The ZnO's energy gap (Eg) was calculated using (Equation 9): Where the wavelength of light is λ, the Planck constant is h, and the speed of light is c.The energy gap of the nanomaterials,' denoted as Eg, was estimated to be 3.4 eV.

FE-SEM, EDX, and BET analysis
Figure 4a,b shows the FE-SEM images that illustrate the surface characteristics of biogenic ZnO and ZnO-ME, in which spherical-shaped aggregates of zinc oxide were formed.Previous research reported by Singh et al.. (2010) used Escherichia coli for the biosynthesis of aggregated micrometer copper oxide nanoparticles [37].This spherical shape of ZnO has also been observed in other biological synthesis methods using Pseudomonas fluorescens, where ZnO was found to aggregate on the surface of bacterial cells [38].The spherical ZnO and ZnO-ME nanoparticles' measured diameters were discovered to be approximately 22.33 and 39 nm, respectively, using imageJ software.Figure 4c,d show the findings of analyzing the elemental composition of zinc oxide nanoparticles using EDX spectroscopy.The peaks at various keV levels are clearly visible in the EDX spectrum, which suggests that the adsorbent contains zinc, oxygen, and carbon.
Relative pressure (P/P0) and adsorbed volume (cc/ g at STP) are plotted in the N 2 adsorption-desorption isotherms, along with the BJH pore size distribution of the ZnO-ME and ZnO nanomaterials were analyzed by BET analyzer, as shown in Figure 4e,f.The isotherm curve closely resembles a normal type IV isotherm graph, demonstrating the H3 hysteresis loop and supporting the mesoporous feature of the nanomaterials.Mesoporous zinc oxide nanoparticles were also obtained in a prior investigation by Varadavenkatesan et al. (2019) [39].This kind of loop exhibits mono-multilayer adsorption, and it is well known that the surface area of the adsorbent has a significant impact on adsorption.The adsorption processes are facilitated by a wide surface area and porosity, which offer a large number of active sites [40].BET analysis for nanomaterials indicated a specific surface area of 7.617 and 33.635 m 2 /g for ZnO and ZnO-ME, respectively.Additionally, the total pore volumes for ZnO, and ZnO-ME were 0.045 and 0.218 cc/g, respectively, and average pore diameters discovered to be 13.839 and 13.882 nm, respectively.

Central composite design
CCD is utilized to assess the most important interdependence between variables as well as the individual effects of each variable on the dye removal efficiency.The significant correlation between each variable and their responses is expressed in (Equation 10) and (Equation 11) below: Removal efficiency (ZnO-ME)

BIOENGINEERED
Where Y stands for dye removal efficiency (%), A stands for dose; B stands for dye concentration, and C stands for contact time between nanomaterials and sulfur black dye.
The ANOVA and model statistics summary utilized to analyze the observed experimental data for sulfur black dye adsorption are shown in Tables 1 and 2, respectively.High F-values     5.It shows the degree to which the experimental results and the theoretical distribution are comparable.It demonstrates that there is good agreement between the model predictions and the actual experimental values of sulfur black dye adsorption.The interaction of two variables and their combined effect on sulfur black dye removal is expressed in 3D model graph as shown in Fig. S3.Plots with curvature verified the degree of interaction between the independent process variables [41].

Plackett-Burman Design (PBD)
Unbiased screening of the variables that significantly influence the sulfur black dye's adsorption onto ZnO and ZnO-ME was done using a two-level PBD factorial design with 13 runs.The results of a three-factor experiment served as its foundation.The significance of dye removal experiment results is stated as an ANOVA, and 3D plots as shown in Table 1 and Figure S4, respectively.The surface plots of the response functions facilitate the combined impact of the investigated factors and their interaction.The correlation between the studied variables and their responses for both the nanoparticles for the sulfur black dye adsorption % are expressed as in (Equation 12) and (Equation 13) below: The high F-value [315.97 for ZnO and 461.12 for ZnO-ME] and probability p-value <0.0001 of the model, as determined by the analysis of the experimental data for Sulfur Black dye adsorption, make it clear that the model is highly significant.A nonsignificant lack of fit test was validated with a low F-value of >0.0500.The fact that the 'adjusted R 2 ' is quite close to the 'predicted R 2 ' indicates that the experimental data fits the model very well.The removal % data are perfectly agreed with the model's anticipated value, as indicated by R 2 value.The linear correlation was depicted by the normal probability plot, showing that the actual experimental values of sulfur black dye adsorption are in good agreement with those predicted by the model, as shown in Figure 6.The results of the PB design can be seen graphically using Pareto charts of standardized effects and normal plot of standardized effects and main effects.Based on standardized effect Pareto plots, the variables that are above the red reference line are statistically significant [42].The length of the depicted bar represents the parameter's weighting as illustrated in Figure S5 (Supplementary data).The following factors were impacted on the sulfur black dye's adsorption rate: In ZnO-ME nanoparticles, the initial dye conc.(B) impact was found higher on the dye adsorption rate which was determined by the 95% confidence interval followed by contact time (C) and dose (A), respectively.For ZnO-ME nanoparticles, the sulfur black dye adsorption rate followed the following order: C (Contact time) > A (dosage) > B (initial dye conc.).At the 95% confidence level, the effect of contact time on the dye adsorption rate was shown to be most significant.
The direction of the effect is additionally shown by the normal plot of the standardized effects (Figure S6).The standardized effects of factors A, B, and C were positive for both nanoparticles.The response increases when these shifts take place from the component's low level to its high level.The parameters' impact between high and low settings is depicted in the main impact plot (Figure S7).The blue circles (Corner point type) on either side of the line represent the mean of the values of the high and low sets.These two points are connected by the blue line.A line with a high slope denotes a considerable influence of the parameter on the response.

Full factorial design (FFD)
Based on the results of a three-factor experiment, a two-level FFD factorial design with nine runs was used to unbiased screen the variables that significantly influence the sulfur black dye adsorption onto ZnO and ZnO-ME nanoparticles.To prevent systemic errors, the experimentation was done in a order.Factorial plots, including the Pareto chart, normal probability, main effect, and surface plots, factors that affect the amount of dye that is adsorbed onto nanoparticles were assessed.The significance of the influence on the adsorption rate was examined using ANOVA and p-value significant levels (refer to Table 1).Additionally, the correlation coefficient (R 2 ) of the model for ZnO-ME and ZnO nanoparticles was 0.994 and 0.999, respectively, confirming the well-fitted statistical model.Also, the predicted R 2 value is close to the adjusted R 2 (adj).Equation 14and Equation 15) are used to represent the dye adsorption by both nanoparticles: Figure S8 illustrates the main effects of each parameter on dye adsorption, while the regression analysis findings are visualized through main effect plots.Only the significant factors are revealed at the 95% confidence level.The ANOVA was also used to create the main effect charts.The relative importance of the primary impacts and their interactions might be seen on the Pareto chart (Figure S9).These values are shown in the Pareto graphic's horizontal columns for each effect.A Student's t-test was used to see whether the calculated effects were statistically different from zero.The significant values are seen above a reference line and inside the 95% confidence interval, as shown in Figure S9 [43].The three main components (A, B, and C) and their interactions (AB only for ZnO-ME) were significant at the level of 0.05.
It is uncertain whether these results are 'real' or 'chance.'The 'real' impacts are identified using a normal probability plot in Figure 7.Each effect has a single point on the plot.In accordance with the normal probability plots, the estimated factors that do not show a substantial impact on the response variables are represented by the points that lie along a line that passes through the center of the middle group of points.According to Geyikci and Büyükgüngör (2013), points that are far from the line most likely indicate the 'real' effects of factor [44].The normal probability plot of residuals for the adsorption rate demonstrated that the set of experimental results closely matched the expected distribution.A normal distribution is suggested by the fact that the majority of the experimental points are reasonably aligned.
Surface plots of the response functions show the major impacts of the variables as well as their interactions.The response surface plots depicting the adsorption rate are presented in (Figure S10).The standardized effects plot's normal plot further displays the effect's direction (Figure S11).Factors A, B, and C (AB only in the case of ZnO) have favorable standardized impacts for both nanoparticles.When these changes from the low level to the high level of the component take place, the reaction increases.

Statistical comparison, performance of optimising models and comparison with other studies
In addition to analyzing the fit of each model using ANOVA, the models' propensities for forecasting the adsorption of sulfur black dye from the aqueous matrix were statistically assessed using the coefficient of determination (R 2 ) and standard deviation (SD).Based on the R 2 and SD values, the accuracy of each model was compared and summarized in Table 2.The FFD model clearly outperformed the other model for data fitting based on the R 2 value for ZnO (0.998) and ZnO-ME (0.993) nanoparticles, which means FFD well predicted the outcomes.All other models also offered good-quality predictions (R 2 ) for the parameters within the design range.In conclusion, FFD also performed well compared to other models while simulating the adsorption process.Interestingly, all the models show the same optimized conditions and removal efficiency on these optimized conditions (i.e.71% for ZnO and 80% for ZnO-ME).These results are comparable with other studies using zinc oxide nanoparticles for dye adsorption, as given in Table 3.

Adsorption isotherm model
Three isotherm models -Langmuir, Freundlich, and Temkin -were used to identify how the adsorption mechanism.These isotherm models were applied to  the experimental datasets generated using FFD.Equation 16, Equation 17, and Equation 18 offer the linear forms of the models under study.
Where C e (mg/L) is the amount of remaining adsorbate in solution and q e (mg/g) denotes the number of dye molecules that interact with the adsorbent.The maximum monolayer adsorption capacity is denoted by the symbol q max (mg/g), and the Langmuir constant b (J/mol) is related to the affinity of binding sites.B = RT/b, where b is the equilibrium binding constant, R is the gas constant, and K T (L/mg) is the Temkin constant related to the heat of adsorption.T (K) is the absolute temperature and Kf (mg/g (L/ mg)1/ n ) is the Freundlich constant.
In order to determine the values of these parameters and correlation coefficients (R 2 ), linear equations were generated using the corresponding plots of C e /q e versus C e , log q e versus log C e , and q e versus lnC e (Figure S12 and Figure S13).With the help of the regression coefficient and parameter values, the models' appropriateness is demonstrated.After comparison, it was shown that the Freundlich model better explained the behavior of dye adsorption, demonstrating multilayer sulfur black dye adsorption onto the ZnO and ZnO-ME nanoparticles.Its coefficient of determination was greater than other examined models for ZnO and ZnO-ME nanoparticles, respectively, at 0.9625 and 0.9549 as can be seen in Table 4. Furthermore, multilayer adsorption is thought to occur on heterogeneous surfaces and on the surface sites of the adsorbent that have different binding energies, according to the Freundlich isotherm, an empirical isotherm model.Additionally, it states that when the concentration of dye in an aqueous solution rises, adsorption increases indefinitely.A heterogeneous adsorption process, which is favored by the Freundlich isotherm, is assumed to occur when there is an uneven distribution of heat on the adsorbent's surface.More variability in the adsorbent surface is indicated by a slope value that is closer to 0 [49].

Adsorption mechanism
Analysis of the dye adsorption mechanism on adsorbent surfaces is challenging because of the complex interactions between the adsorbent and adsorbate.The features of the adsorbent surface, functional groups, and molecular structure of the dye molecules all play a key part in the interaction and binding of the dye molecules to the surface of the adsorbent.The hypothesized method may be demonstrated by comparing the FTIR spectra of ZnO and ZnO-ME nanoparticles before and after dye adsorption, as shown in Figure 8a.The IR peak of both nanoparticles prior to adsorption has already been covered in Section 3.2.3.The hydrogen bond interaction between the dye mlecules and adsorbent surface may be the cause of the rise in band strength at 3440 cm −1 in the FTIR spectra of sulfur black dye adsorbed ZnO and ZnO-ME nanoparticles.Kataria et al., (2022) have demonstrated the significance of the O-H group in the adsorption of dye molecules on zinc oxide nanoparticles [3].This study illustrates its function in the adsorption of color molecules as ZnO-ME exhibits more dye removal than ZnO due to a sharp increase in peaks at 3440 cm −1 .Strong peaks at 2926.26 cm −1 and 2336 cm −1 are associated to the -C=C stretching vibration of alkynes, and a minor decline in these peaks was also noticed, demonstrating the function of -C=C stretching in adsorption.The C-H bond is responsible for the peak between 1050 and 900 cm −1 , and completely diminishing of these peaks highlights its function in the adsorption of the sulfur black dye due to hydrogen bonding.Peak at 1633.39 cm −1  Zeta potential values show negative charge on the surface of nanomaterials and dye molecules show positive charge (N + ) on their surface, showing electrostatic attraction and this helps in the adsorption of dye [52].From the freundlich adsorption isotherm, it can be seen the K f value for ZnO and ZnO-ME NPs 2.9 and 3.8, respectively, which shows the physisorption through weak Van der Waals forces.Further, the role of physisorption is shown by the small size of the NPs which is shown by the blue shift of UV-VIS and the small size of NPs can also be seen in the SEM images.BET surface area of 7.617 and 33.635 m 2 /g for ZnO and ZnO-ME, respectively, also supports physisorption due to high surface area.A potential mechanism for color adsorption for both the nanocomposites is depicted in Figure 8b.

Conclusion
ZnO and ZnO-ME nanoparticles having diameters of 22.33 and 39 nm, respectively, with mesoporous structure were synthesized using Bacillus paramycoides for SB dye removal.Out of used optimized models (CCD, PBD, and FFD), FFD showed better predictability, with R 2 values of 0.9985 and 0.9937 for ZnO and ZnO-ME, respectively, and removal efficiency on these optimized conditions (i.e.71% for ZnO and 80% for ZnO-ME).The Freundlich isotherm model accurately predicted the physisorption of sulfur black dye, with R 2 values of 0.96 for ZnO and 0.95 for ZnO-ME, and maximum adsorption capacities of 2.77 mg/g for ZnO and 4.12 mg/g for ZnO-ME.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Figure 1 .
Figure 1.Synthesis of zinc oxide nanoparticles via isolated bacterial strain.
K f and high 'n' values indicate high sorption throughout the concentration range, while low K f and high 'n' values indicate low sorption.Low 'n' values indicate high sorption at strong solute concentration.Beneficial sorption occurs at n values between 1 and 10 [23].

Fig SI- 2
Fig SI-2).The correlation of ZP data with colloid stability is a widely used application in nanoparticle characterization.According to established classifications, nanomaterial dispersions with ZP values between ± 0 and ±10 mV, ±10 and ±20 mV, ±20 and ±30 mV, and > ±30 mV are, respectively, categorized as highly unstable, reasonably stable, moderately stable, and highly stable[25].ZnO nanoparticles are moderately stable, whereas ZnO-ME nanoparticles are reasonably stable, according to zeta potential measurements.These zinc nanoparticles are coated with negatively charged biomolecules, as indicated by the negative value[26].

Figure 3 .
Figure 3. Characterization of ZnO and ZnO-ME.(a) FT-IR of ZnO and ZnO-ME (before dye adsorption), (b) XRD of ZnO and ZnO-ME and (c) UV-VIS spectrum of ZnO and ZnO-ME.

Figure 4 .
Figure 4. (a,b) surface morphology analysis, (c,d) EDX spectrum and (e,f) nitrogen adsorption-desorption curve and pore size distribution for ZnO-ME and ZnO.

FFD-for
Analysis of variance (ANOVA) for the adsorption of suflur black dye using ZnO and ZnO-ZnO and ZnO-ME (303.79 and 431.

Figure 5 .
Figure 5. Predicted against actual values plots for (a) for ZnO and (b) for ZnO-ME and normal percentage probability plots for (c) for ZnO and (d) for ZnO-ME for sulfur black dye adsorption.

Figure 6 .
Figure 6.Residual plots for PBD for ZnO and ZnO-ME nanoparticles.

Figure 7 .
Figure 7. Residual plots for FFD for ZnO and ZnO-ME nanoparticles.
represents amine group and alteration in intensity and slight shift in position after dye adsorption shows its role in adsorption process and formation of amide group[50].π-π interaction is presented between the benzene rings of the dye and the adsorbent.A sharp broad peak at 619-480 cm −1 corresponds to the Zn-O stretching, FTIR of dye-laden ZnO and ZnO-ME NPs showed slight change in the peak, indicating the significance of it in the process of dye adsorption.Adsorption of dye molecules onto ZnO and ZnO-ME NPs showed hydrogen bonding, electrostatic attraction, π-π interaction in the adsorption of sulfur black dye which suggest chemisorption[51].

Figure 8 .
Figure 8.(a) FT-IR spectra of nanomaterials before and after adsorption; (b) potential adsorption mechanism.

Table 1 .
Analysis of variance (ANOVA) for the adsorption of sulfur black dye using ZnO and ZnO-ME nanoparticles

Table 2 .
Model summary and statistical analysis.

Table 3 .
Comparison of adsorption percentage of zinc oxide nanoparticles for various dyes.

Table 4 .
Isotherm parameters for sulfur black dye adsorption on ZnO and ZnO-ME.