Enhanced Adsorptive Removal of β -Estradiol from Aqueous and Wastewater Samples by Magnetic Nano-Akaganeite: Adsorption Isotherms, Kinetics, and Mechanism

: A surfactant-free method was used to synthesize iron oxyhydroxide (akaganeite, β -FeOOH) nanorods and characterized using Fourier transform infrared spectroscopy (FTIR), X-ray di ﬀ raction (XRD), scanning electron microscopy combined with energy-dispersive X-ray spectroscopy (SEM-EDS), and transmission electron microscopy (TEM). The synthesized nanoadsorbent was applied for the adsorptive removal of β -estradiol from aqueous solutions. The parameters a ﬀ ecting the adsorption were optimized using a multivariate approach based on the Box–Behnken design with the desirability function. Under the optimum conditions, the equilibrium data were investigated using two and three parameter isotherms, such as the Langmuir, Freundlich, Dubinin–Radushkevich, Redlich–Peterson, and Sips models. The adsorption data were described as Langmuir and Sips isotherm models and the maximum adsorption capacities in Langmuir and Sips of the β -FeOOH nanorods were 97.0 and 103 mg g − 1 , respectively. The adjusted non-linear adsorption capacities were 102 and 104 mg g − 1 for Langmuir and Sips, respectively. The kinetics data were analyzed by ﬁve di ﬀ erent kinetic models, such as the pseudo-ﬁrst order, pseudo-second order, intraparticle, as well as Boyd and Elovich models. The method was applied for the removal β -estradiol in spiked recoveries of wastewater, river, and tap water samples, and the removal e ﬃ ciency ranged from 93–100%. The adsorbent could be reused up to six times after regeneration with acetonitrile without an obvious loss in the removal e ﬃ ciency (%RE = 95.4 ± 1.9%). Based on the results obtained, it was concluded that the β -FeOOH nanorods proved to be suitable for the e ﬃ cient removal of β -estradiol from environmental matrices.


Introduction
Water is the most important natural resource for the survival of all organisms [1]. Due to the increase in urbanization, climate change, industrial production, and population growth, the quality of water can be negatively affected. Pollutants entering soil or water pose a variety of threats to the

Instrumentation
All pH measurements were done using an OHAUS ST series pen pH meter (Nanikon, Switzerland). The adsorption studies were carried out using the Branson 5800 Ultrasonic Cleaner (Danbury, CT, USA). The FTIR measurements were taken using the KBr pellet technique and recorded in the region 400-4000 cm −1 using a Perkin Elmer Spectrum 100 spectrometer. Scanning electron microscopy (SEM, TESCAN VEGA 3 XMU, LMH instrument, Czech Republic) coupled with energy-dispersive X-ray spectroscopy (EDS) was used to study the morphology and elemental composition of the adsorbent at an accelerating voltage of 20 kV. The transmission electron microscopic image was captured using TEM (JEM-2100, JEOL, Tokyo, Japan). An Agilent high-performance liquid chromatography (HPLC) 1200 Infinity series, equipped with a photodiode array detector (Agilent Technologies, Waldbronn, Germany), was used for all analysis. The separation was done using an Agilent Zorbax Eclipse Plus C18 column (3.5 µm × 150 mm × 4.6 mm) (Agilent, Newport, CA, USA) operated at an oven temperature of 25 • C. The chromatograms were recorded using a 1.00 mL min −1 flow rate, a solvent mixture of 55% mobile phase A (water) and 45% mobile phase B (acetonitrile), and an adsorption wavelength of 260 nm using an isocratic elution system.

Sampling and Storage
Influent and effluent wastewater samples were collected from a local wastewater treatment plant (WWTP, Gauteng, South Africa) in winter 2019. River water samples were collected from Apies River, which is upstream from the wastewater treatment plant. Tap water samples were collected from the University of Johannesburg Doornforntein campus. All samples were collected in 1 L glass bottles and refrigerated until use.

Preparation of Akaganeite Nanorods
Akaganeite nanoparticles were prepared according to the method described by [23]. Briefly, a 5 mol L −1 solution of NaOH was prepared in 30 mL deionized water and stirred at room temperature in a beaker equipped with a magnetic stirrer. A 25 mL 2 mol L −1 FeCl 3 was rapidly added to the NaOH solution using a syringe; the stirring was continued for an hour. Ethanol was then added, and the resulting precipitate was washed several times with a 1:1 ethanol-water solution. The precipitate was then dried in an oven at 50 • C for 5 h.

Batch Adsorption of β-Estradiol
The adsorption was done using the batch incubation method described in [24]. The method was optimized using the Box-Behnken design. The levels of independent factors, such as sample pH, contact time (CT), and mass adsorbent (MA), are presented in Table S1. The appropriate amount of adsorbent was then weighed and added to a glass bottle. Aliquots of a 10 mL sample solution containing 6 mg L −1 β-estradiol were added into the sample bottle. The samples were then agitated using an ultrasound water bath; thereafter, the adsorbent and supernatant were separated using an external magnet. Thereafter, 1 mL of sample was filtered via a 0.22 µm PVDF filter before analysis with HPLC. All the experiments were carried out in triplicates, after which the adsorption removal efficiency was calculated using Equation (1): where C 0 and C e are the initial and equilibrium concentrations, respectively. An optimization process is a critical stage in analytical chemistry. This is because of regulations due to the high cost of reagents and quantity of variables that affect the development and validation of analytical methods [25]. Multivariate optimization is an important technique because it takes less time, resources, and effort than the univariate approach. The design of experiments is important in analytical applications [25]. Multivariate optimization strategies allow to efficiently minimize the parameter space-saving time and resources when determining the optimum conditions, here for the deployment of an adsorbent material [25]. For the adsorption of β-estradiol onto akaganeite, the (1) pH of the sample, (2) mass of adsorbent (MA), and (3) contact time (CT) were optimized using a multivariate approach, as presented in Table S1.
Under optimum conditions, the equilibrium isotherms and kinetics of adsorption were studied. The β-estradiol solutions had concentrations ranging from 1 to 9 mg L −1 . The equilibrium isotherm studies were done in triplicates and the β-estradiol concentration adsorbed onto akaganeite nanorods was calculated using Equation (2): where C 0 and C e are the initial and equilibrium concentrations (mg. L −1 ) of β-estradiol, respectively, and m (g) and V (L) represent the mass of the adsorbent and volume of the sample used. The adsorption kinetics was investigated by the addition of 10 mL aliquots (6 mg. L −1 ) into 7 glass bottles. The bottles were then sonicated for 1-60 min for each sample before magnetic separation and HPLC analysis.

Method Validation and Real Sample Analysis
The collected wastewater (influent and effluent), river, and tap water samples were analyzed before adsorption and β-estradiol could not be detected in the sample. As a result, the method was validated using spiked sample recoveries; this was due to the absence of certified reference materials for emerging pollutants. To validate the method, wastewater effluent samples were spiked at two levels (4 and 8 mg. L −1 ) with β-estradiol by appropriate dilution with the effluent water. The validation experiments were carried out by adding 10 mL of the sample (adjusted to pH 5.25) onto 42.5 mg of akaganeite in a capped glass bottle. The mixture was then agitated in an ultrasonic bath for 60 min, after which the supernatant and adsorbent were separated using an external magnet before filtering and HPLC analysis of the supernatant. Figure 1 shows the Fourier transform infrared spectrum of the as-synthesized akaganeite. The peaks at 3437 and 3135 cm −1 were ascribed for OH stretching and the bending at 1620 cm −1 and 1009 cm −1 were assigned to the absorption bands of the iron oxyhydroxides arising from the Fe-OH and Fe-O vibration of the akaganeite [19,26]. The bending peak at 650 cm −1 is reported to be the characteristic vibrations of Fe-O in β-FeOOH [26]. Lastly, the OH bending band around the 1431 cm −1 region corresponded to the β-FeOOH FTIR spectrum reported in the literature [19,26].

Characterisation of the Akaganeite Nanoparticles
Akaganeite nanoparticles are known to be low-crystallinity materials, which during XRD characterization are indicated by rather low-intensity peaks [27]. In Figure 2, an exemplary XRD data set for the synthesized akaganeite nanoparticles are presented. From the XRD pattern, the characteristic low-intensity peaks were observed, confirming low crystallinity. According to [26], decreasing the pH of the synthesis media results in more crystalline akageneite, while increasing the pH of the synthesis media results in less crystalline akaganeite nanoparticles, which was a similar case with the present study (i.e., sodium hydroxide was used). According to [27], the diffraction peaks at around 2θ = 20 • and 40 • are characteristic nano-akaganeite peaks and are evident in Figure 2. Akaganeite nanoparticles are known to be low-crystallinity materials, which during XRD characterization are indicated by rather low-intensity peaks [27]. In Figure 2, an exemplary XRD data set for the synthesized akaganeite nanoparticles are presented. From the XRD pattern, the characteristic low-intensity peaks were observed, confirming low crystallinity. According to [26], decreasing the pH of the synthesis media results in more crystalline akageneite, while increasing the pH of the synthesis media results in less crystalline akaganeite nanoparticles, which was a similar case with the present study (i.e., sodium hydroxide was used). According to [27], the diffraction peaks at around 2θ = 20° and 40° are characteristic nano-akaganeite peaks and are evident in Figure  2.   Akaganeite nanoparticles are known to be low-crystallinity materials, which during XRD characterization are indicated by rather low-intensity peaks [27]. In Figure 2, an exemplary XRD data set for the synthesized akaganeite nanoparticles are presented. From the XRD pattern, the characteristic low-intensity peaks were observed, confirming low crystallinity. According to [26], decreasing the pH of the synthesis media results in more crystalline akageneite, while increasing the pH of the synthesis media results in less crystalline akaganeite nanoparticles, which was a similar case with the present study (i.e., sodium hydroxide was used). According to [27], the diffraction peaks at around 2θ = 20° and 40° are characteristic nano-akaganeite peaks and are evident in Figure  2.  Scanning electron microscopy (SEM) was used to investigate the surface morphology of the akaganeite nanomaterial combined with energy-dispersive X-ray (EDX) for elemental analysis. The major components of the akaganeite nanomaterial were determined by EDX and given in the spectrum represented in Figure 3B. SEM images ( Figure 3A,C) revealed aggregated akaganeite nanoparticles, which resulted from the absence of a surfactant during the synthesis, as was reported in [19]. The yellow frame in Figure 3C indicates the area used for evaluation by EDX.
Scanning electron microscopy (SEM) was used to investigate the surface morphology of the akaganeite nanomaterial combined with energy-dispersive X-ray (EDX) for elemental analysis. The major components of the akaganeite nanomaterial were determined by EDX and given in the spectrum represented in Figure 3B. SEM images ( Figure 3A,C) revealed aggregated akaganeite nanoparticles, which resulted from the absence of a surfactant during the synthesis, as was reported in [19]. The yellow frame in Figure 3C indicates the area used for evaluation by EDX. Transmission electron microscopy (TEM) studies, as shown in Figure 4, indicate the characteristic rod-shaped akaganeite structures next to more spherical/octagonal particles, as described by [19]. While no surfactant was used herein, to better control the shape of the obtained akageneite material, a surfactant is needed, ensuring more ordered mesostructures. However, one aim of the present study was to ensure a straightforward and possibly simple synthesis route. Hence, a more uniform particle shape/size was traded against the simplicity of the synthesis strategy, which in absence of surfactants does not require further processing steps, including autoclaved heating and time-consuming Soxhlet extraction [19].  Transmission electron microscopy (TEM) studies, as shown in Figure 4, indicate the characteristic rod-shaped akaganeite structures next to more spherical/octagonal particles, as described by [19]. While no surfactant was used herein, to better control the shape of the obtained akageneite material, a surfactant is needed, ensuring more ordered mesostructures. However, one aim of the present study was to ensure a straightforward and possibly simple synthesis route. Hence, a more uniform particle shape/size was traded against the simplicity of the synthesis strategy, which in absence of surfactants does not require further processing steps, including autoclaved heating and time-consuming Soxhlet extraction [19].
Scanning electron microscopy (SEM) was used to investigate the surface morphology of the akaganeite nanomaterial combined with energy-dispersive X-ray (EDX) for elemental analysis. The major components of the akaganeite nanomaterial were determined by EDX and given in the spectrum represented in Figure 3B. SEM images ( Figure 3A,C) revealed aggregated akaganeite nanoparticles, which resulted from the absence of a surfactant during the synthesis, as was reported in [19]. The yellow frame in Figure 3C indicates the area used for evaluation by EDX. Transmission electron microscopy (TEM) studies, as shown in Figure 4, indicate the characteristic rod-shaped akaganeite structures next to more spherical/octagonal particles, as described by [19]. While no surfactant was used herein, to better control the shape of the obtained akageneite material, a surfactant is needed, ensuring more ordered mesostructures. However, one aim of the present study was to ensure a straightforward and possibly simple synthesis route. Hence, a more uniform particle shape/size was traded against the simplicity of the synthesis strategy, which in absence of surfactants does not require further processing steps, including autoclaved heating and time-consuming Soxhlet extraction [19].

Optimisation of Removal Procedure
A Pareto chart (see Figure S1) was used to identify critical factors during the adsorption process by evaluating the variance within the experimental design matrix [28]. The Pareto chart is based on analysis of variance (ANOVA), which is a linear model. As a result, during the analysis of the Pareto chart, only the linear factors (i.e., factors with the index L) were included, while quadratic factors (with the index Q) were not considered [29]. In the analysis of a Pareto chart, for a factor to be considered statistically significant, its bar should pass the redline that indicates a p-value of 0.05 (95% confidence level) [30,31]. The ultimately dominating factors were determined to include the (1) pH of the sample, (2) the mass of adsorbent (MA), and (3) the contact time (CT). Of these three factors, the only one beyond the 95% confidence line (red line in Figure S1) was the pH.
The response surface methodology (RSM) was used to investigate (a) the interaction and (b) the quadratic effects of the main parameters pH, MA, and CT, respectively, using the data generated via a Box-Behnken experimental design, yielding 3D surface plots ( Figure S2) for these three parameters. As it can be seen from Figure S2A, increasing the mass of the adsorbent increases the percentage removal (% RE) with any pH below 7, and at a pH of around 7 the maximum % RE could be observed. The same could be observed from the interaction of pH and contact time. Figure S2C shows the interaction between CT and MA; masses between 25 and 10 together with a CT of 20 min gives a minimum 5 RE while a CT and MA above 30 min and 30 mg results in higher percentage removal (% RE). The plots were used in conjunction with the profiles of desirability to determine the optimum conditions.
Using the desirability function (DF, Figure S3) allows the simultaneous estimation of the optimal values for all the investigated factors. Desirability always takes values within the 0-1 range where 0 indicates the least desired result (33.4% removal), 0.5 being the central point 69.4%, and 1 being the most desirable value assigned a % removal of 105 [32,33]. According to the desirability profiles, the optimal conditions for the adsorption process were a pH of 5.25, MA of 42.5 mg, and CT of 60 min. The DF predicted values were then used to confirm the optimum conditions of the adsorption method.
To understand the effect of pH on the adsorption, the point of zero charge of akaganeite was considered. At a pH lower than the point of zero charge, the surface of the material is positively charged, while at a pH higher than the point of zero charge, the surface of the adsorbent is negatively charged [34][35][36]. The point of zero charge for the akaganeite nanorods as described in the literature is around 7-7.5 [21]. The optimum pH for the adsorption was found to be 5.25, which is lower than the point of zero charge of akageneite. This means that the surface of the material was positively charged. β-estradiol also has a pKa of 10.71, meaning that at the optimum pH, the analyte is in its molecular state. This translates to the presence of the OH groups on the analyte to be available for ionic interaction between the β-FeOOH and β-estradiol.

Adsorption Isotherms
The ratio between the concentration of adsorbed analyte and the amount of analyte remaining in the supernatant solution at equilibrium conditions provides the absolute removal efficiency [24]. The interaction of the analyte species with the adsorbent material reaches a dynamic equilibrium after a certain incubation time, which may be described via so-called isotherm models [37]. In the present study, the isotherms for β-estradiol adsorbing at the surface of β-FeOOH nanorods were studied at a pH of 5.25 with a β-estradiol concentration of 1-10 mgL −1 , an adsorbent mass (MA) of 42.5 mg, and contact time (CT) of 60 min at 25 • C. It was observed that the adsorption capacity increased with an increase in initial β-estradiol concentration ( Figure 5A). In this study, five isotherm models were used to describe the equilibrium. These include the Langmuir, Freundlich, Dubinin-Radushkevich, Redlich-Peterson, and Sips isotherm models, and the resulting plots are shown in Figure 5B-F. The R 2 values for the different models were used to determine which model best fits the adsorption process. The Langmuir and Freundlich models were used to determine the primary adsorption mechanism. As seen in Table 1 the adsorption process fitted the Langmuir model with an R 2 of 0.9996, which was higher than the 0.8811 for the Freundlich model.
The β-estradiol adsorption data were fitted with Langmuir, Freundlich, Dubinin-Radushkevich, and Redlich-Peterson isotherms models, as summarized in Table 1. Comparing the obtained goodness-of-fit values in Table 1, the experimental data followed a Langmuir isotherm (R 2 = 0.9996), which was higher than for any other isotherm. This implies that adsorption occurs homogeneously at the surface sites of the β-FeOOH nanomaterial and that the adsorption of β-estradiol may be assumed occurring in a monolayer fashion. Thus, a Langmuir constant (K L ) of 16 L g −1 and the maximum adsorption capacity of 97.0 mg g −1 was obtained. The calculated separation factor (R L ) of 0.10-0.25 for the β-estradiol adsorption on the β-FeOOH nanomaterial, according to the literature, demonstrates a favorable adsorption [38,39]. The β-estradiol adsorption data were fitted with Langmuir, Freundlich, Dubinin-Radushkevich, and Redlich-Peterson isotherms models, as summarized in Table 1. Comparing the obtained goodness-of-fit values in Table 1, the experimental data followed a Langmuir isotherm (R 2 = 0.9996), which was higher than for any other isotherm. This implies that adsorption occurs homogeneously at the surface sites of the β-FeOOH nanomaterial and that the adsorption of β-estradiol may be assumed occurring in a monolayer fashion. Thus, a Langmuir constant (KL) of 16 L g −1 and the maximum adsorption capacity of 97.0 mg g −1 was obtained. The calculated separation factor (RL) of 0.10-0.25 for the β-estradiol adsorption on the β-FeOOH nanomaterial, according to the literature, demonstrates a favorable adsorption [39,40].
The Dubinin-Radushkevich isotherm is conventionally used to differentiate between chemical and physical adsorption by evaluating the mean energy (E). It is said that an E value that is lower than 8 kJ mol indicates physisorption dominance while an E value above 8 indicates chemisorption dominance [41,42]. The calculated E value was 8.5 kJ mol −1 , which indicates a chemisorption process. The Redlich-Peterson model combines both the Freundlich and Langmuir isotherms, whereby B is the exponent from the linear plot ranging between 0 and 1. If B equals 1, the model reduces to the Langmuir equation, and if B = 0, it reduces to the Freundlich equation [43]. From the obtained results, B has a value of 0.88, which again favors the Langmuir model to describe the interaction between β-estradiol and the β-FeOOH nanorods. To further investigate monolayer adsorption, the Sips model was used. When the Sips model exponent nS = 1, the adsorption process follows a Langmuir The Dubinin-Radushkevich isotherm is conventionally used to differentiate between chemical and physical adsorption by evaluating the mean energy (E). It is said that an E value that is lower than 8 kJ mol indicates physisorption dominance while an E value above 8 indicates chemisorption dominance [40,41]. The calculated E value was 8.5 kJ mol −1 , which indicates a chemisorption process. The Redlich-Peterson model combines both the Freundlich and Langmuir isotherms, whereby B is the exponent from the linear plot ranging between 0 and 1. If B equals 1, the model reduces to the Langmuir equation, and if B = 0, it reduces to the Freundlich equation [42]. From the obtained results, B has a value of 0.88, which again favors the Langmuir model to describe the interaction between β-estradiol and the β-FeOOH nanorods. To further investigate monolayer adsorption, the Sips model was used. When the Sips model exponent n S = 1, the adsorption process follows a Langmuir isotherm, and when n S >1, the data is considered to fit the Freundlich model [43][44][45]. In this work, the Sips model exponent was equal to 1, further confirming monolayer adsorption. Table 1. Adsorption models and parameters investigated in the present study.

Model Parameter R 2
Langmuir C e q e = 1 q max ·K L · C e + 1 q max ,R L = The adsorption data were also fitted to the nonlinear isotherm models. This because the major problem encountered when using linearized models is the violation of theories behind each model using the modification of the original equation [46][47][48]. In such cases, the linearized model may give best-fitting parameters as opposed to the original nonlinear model [46][47][48]. Therefore, to compare the linearised and nonlinearized regression models, the data were fitted to the original nonlinear isotherm equation (Table 2). The better-fitting kinetic of the nonlinear equations was investigated using residual standard error (RSE) and R 2 . According to the RSE and R 2 , the best model is Redlich-Peterson, followed by the Sips, Langmuir, Dubinin-Radushkevich, and Freundlich isotherms. In this study, the isotherm parameters obtained for both linear and nonlinear models were comparable, indicating that the two regression models were not significantly different at the 95% confidence level. Besides, these findings proved that the linearized model did not violate the theory behind the original models.

Adsorption Kinetics
Batch experiments were carried out to investigate the adsorption kinetics of β-estradiol ( Figure 6A) using 42.5 mg of β-FeOOH nanorods, a concentration of 5 mg L −1 β-estradiol, and a pH of 5.25 at room temperature (25 • C). As seen, the rapid uptake of β-estradiol by the adsorbent happened from 5-20 min and the equilibrium was attained between 20 and 60 min. The kinetics data were fitted into various kinetics equations, as seen in Table 3. The plots for the pseudo-first order (ln (q e -q t ) vs. t), pseudo second order (t/q t vs. t), and Elovich (q t vs. ln t) kinetic models are shown in Figure 6B-D. The derived parameters together with the correlation coefficient are summarized in Table 3. As seen in Figure 6B and Table 3, the correlation coefficients for the pseudo-first-order model was 0.6555. The experimental adsorption capacity (97.2 mg g −1 ) was almost four times higher than the calculated (26.8 mg g −1 ) adsorption capacity. These finding suggested that the adsorption of β-estradiol onto β-FeOOH nanorods was not an ideal pseudo-first-order reaction. The R 2 value for the pseudo-second-order model is higher than that of the first-order model. This suggests that the kinetic equilibrium data were best described by a pseudo-second-order model. The calculated q e value of the pseudo-second-order model is also 100 mg g −1 , which is in agreement with the experimentally obtained value of q e at 97.2 mg g −1 . This agreement suggested that the adsorption process was driven by chemisorption involving the electrostatic interaction between the positively charged adsorbent and lone pairs of electrons of β-estradiol. The Elovich kinetic model was used to describe the chemisorption process on the surface of the β-FeOOH nanorods. The kinetics data were fitted to the Elovich equation and gave a relatively good R 2 value (0.9595), which shows that the uptake of β-estradiol onto the β-FeOOH nanorods involved chemisorption [49]. Table 2. Non-linear adjusted adsorption isotherms.

Model
Parameter Adjusted R 2 R 2 Langmuir q e =  Table 3. Kinetic models and parameters investigated in the present study.

Parameter R 2
Pseudo-First order ln(q e − q t ) = ln q e − K 2 t q e (mg/g) = 26.8 k 1 (mg/g min 1/2 ) = 0.063 0.6555 Pseudo-Second order 1 q t = 1 K 2 q 2 e + 1 q e t q e = 100 k t (mg/g min 1/2 ) = 0.004 0.9989 α (mg/g min) = 104 β (g/mg) = 0.058 0.9595    To understand the rate-determining step, various intraparticle diffusion models were applied, as summarized in Table 3. From the plots of the intraparticle diffusion ( Figure 6E), it is evident a two-step adsorption process occurred during the uptake of β-estradiol by the akageneite nanorods. The first step is based on the so-called film diffusion or surface adsorption, i.e., diffusion of β-estradiol from the solution to the external surface of the adsorbent [24]. The plots also do not go through the origin, which indicates that the intraparticle diffusion is not the rate-determining step. The second step suggested that the adsorption process proceeds by intraparticle diffusion and the intraparticle diffusion happened at a higher concentration. Kid 1 derived from the film diffusion is exceeding the Kid 2 of the intraparticle diffusion, which indicates that the intraparticle diffusion is a slow process. The value of the intercept (C) is also larger for the intraparticle diffusion, thus providing a more relevant contribution to the adsorption process [50]. Because the intraparticle diffusion plot shows a non-linear nature, the Boyd model ( Figure 6F) was used to investigate the definite rate-controlling step. According to Hu et al. [49], when the Boyd plot passes through the origin it means that particle diffusion is the dominant mechanism that governs the adsorption process. However, as seen in Figure 6F, the plot did not pass through the origin, suggesting that the adsorption process might be governed by a boundary layer diffusion mechanism [51]. However, the high correlation coefficient of the Elovich model meant that between the two steps of the rate-determining steps, intraparticle diffusion was more prominent when compared to film diffusion [52,53]. This validates the results obtained during the optimization, which essentially results in the fact that electrostatic interactions are driving the β-estradiol adsorption onto akageneite. This is also confirmed in the literature for the adsorption of β-estradiol onto carbon nanotubes [16].
Nonlinear regression of the kinetics data for each model was investigated using the adjusted R 2 , R 2 , and residual standard error. The calculated kinetics parameter and the correlation coefficients are illustrated in Table 4. Unlike in the linearized regression model where the trend was pseudo second order > Elovich > pseudo-first order, the data was best fitted by the pseudo-second order, followed by the pseudo-first-order and Elovich models. When comparing the kinetic parameters between the non-linear and linear models, it could be seen that the calculated q e for the linearised pseudo-first order was the only model that violated the theory of the original non-linear model. Table 4. Adjusted non-linear adsorption kinetics.

Model
Parameter Adjusted R 2 R 2 Pseudo-First order q t = q e 1 − e −k 1 t q e (mg/g) = 94.

Adsorption Thermodynamics Studies
Thermodynamic studies were investigated to ascertain the dominant adsorption mechanisms (that is, physisorption or chemisorption). The adsorption process was carried at different temperatures (298, 303, 308, and 313 K) and the thermodynamic parameters, such as enthalpy (∆H • ), entropy (∆S • ), and Gibbs energy (∆G • ), were calculated according to the expressions reported by References [46,48,50,52]. The graph of lnK C against 1/T(K −1 ) is presented in Figure 7. The values of ∆H • and ∆S • were estimated from the slope and intercept of the plot. The ∆H • and ∆S • were found to be 45.5 kJ mol −1 and 186 J mol− 1 K −1 while the Gibbs energies were ranged from −9.86 to −12.6 kJ mol −1 . The positive ∆H • reveals that the adsorption process was endothermic in nature and the positive ∆S • suggested an increase in the randomness at the boundary of the solid/liquid phases during the adsorption process [46,48,50,52]. The negative ∆G • indicate the spontaneity of the adsorption process. Moreover, the ∆H • value was higher than 20.9 kJ/mol, confirming that the chemisorption was dominant [52].

Adsorption Mechanism
In order to gain further insight into the adsorption process of β-estradiol onto akaganeite, the adsorption was characterized by performing FTIR analysis (Figure 8a-c) on akaganeite before and after adsorption, as well as on the analyte. The adsorption of β-estradiol resulted in peaks in the 1500-1000 cm −1 region, which were ascribed to various C-O, C-H, and C-OH bonds. As evident in Figure 8c, the emergence of the prominent peaks at 2800-2900 cm −1 , corresponding to the stretching vibration of the C-H bonds of the β-estradiol, and the broadening of the OH vibration peak at 3135 cm −1 confirmed the presence β-estradiol on the surface of the adsorbent [54]. Figure 8a revealed that akageneite has distinctive vibration peaks of hydroxyl groups and according to the literature the pKa of -OH is around 9.5-13 [55,56]. This implies that at pH values below the pKa values, the hydroxyl groups of the adsorbent are protonated [55,56]. In addition, the pH pzc of the adsorbent was reported to be between 7 and 7.5, suggesting that adsorbent is positively charged below the pH pzc . Moreover, the pKa value of β-estradiol is reported to be 10.71, suggesting the analyte existed as a neutral species. There, the adsorption mechanism was dominated by the electrostatic attraction between the positively charged adsorbent and the lone pair of electrons of the oxygen atoms present in the analyte. Furthermore, the presence of the electron-rich aromatic ring in the structure of the analyte resulted in cation-π interactions between the positively charged adsorbent surface and the π-system of the β-estradiol molecule. These finding are in line with the observation in the optimization data where pH influenced the interactions between the positively charged adsorbent and the electron-rich β-estradiol molecule [6]. Moreover, the enhanced peak intensities, shifting of adsorption bands at 1700-1500 cm −1 , and disappearance of some peaks at 650-600 cm −1 , as seen in Figure 8b, confirmed that the adsorption mechanism was dominated by the chemisorption process, which is in line with the adsorption isotherms, kinetics data, and thermodynamics studies.

Adsorption Thermodynamics Studies
Thermodynamic studies were investigated to ascertain the dominant adsorption mechanisms (that is, physisorption or chemisorption). The adsorption process was carried at different temperatures (298, 303, 308, and 313 K) and the thermodynamic parameters, such as enthalpy (ΔH°), entropy (ΔS°), and Gibbs energy (ΔG°), were calculated according to the expressions reported by References [47,49,51,53]. The graph of lnKC against 1/T(K −1 ) is presented in Figure 7. The values of ΔH° and ΔS° were estimated from the slope and intercept of the plot. The ΔH° and ΔS° were found to be 45.5 kJ mol −1 and 186 J mol− 1 K −1 while the Gibbs energies were ranged from −9.86 to −12.6 kJ mol −1 . The positive ΔH° reveals that the adsorption process was endothermic in nature and the positive ΔS° suggested an increase in the randomness at the boundary of the solid/liquid phases during the adsorption process [47,49,51,53]. The negative ΔG° indicate the spontaneity of the adsorption process. Moreover, the ΔH° value was higher than 20.9 kJ/mol, confirming that the chemisorption was dominant [53].

Adsorption Mechanism
In order to gain further insight into the adsorption process of β-estradiol onto akaganeite, the adsorption was characterized by performing FTIR analysis (Figure 8a-c) on akaganeite before and after adsorption, as well as on the analyte. The adsorption of β-estradiol resulted in peaks in the 1500-1000 cm −1 region, which were ascribed to various C-O, C-H, and C-OH bonds. As evident in Figure 8c, the emergence of the prominent peaks at 2800-2900 cm −1 , corresponding to the stretching vibration of the C-H bonds of the β-estradiol, and the broadening of the OH vibration peak at 3135 cm −1 confirmed the presence β-estradiol on the surface of the adsorbent [55]. Figure 8a revealed that akageneite has distinctive vibration peaks of hydroxyl groups and according to the literature the pKa of -OH is around 9.5-13 [56,57]. This implies that at pH values below the pKa values, the hydroxyl groups of the adsorbent are protonated [56,57]. In addition, the pHpzc of the adsorbent was reported to be between 7 and 7.5, suggesting that adsorbent is positively charged below the pHpzc. Moreover, the pKa value of β-estradiol is reported to be 10.71, suggesting the analyte existed as a neutral species. There, the adsorption mechanism was dominated by the electrostatic attraction between the positively charged adsorbent and the lone pair of electrons of the oxygen atoms present in the analyte. Furthermore, the presence of the electron-rich aromatic ring in the structure of the analyte resulted in cation-π interactions between the positively charged adsorbent surface and the π-system of the β-estradiol molecule. These finding are in line with the observation in the optimization data where pH influenced the interactions between the positively charged adsorbent and the electron-rich β-estradiol molecule [6]. Moreover, the enhanced peak intensities, shifting of

Analysis of Real Samples
During the experimental studies, the actual removal of β-estradiol from the spiked aqueous samples was investigated. The two-level spiked effluent wastewater samples showed that the adsorption method could remove between 93 and 100%. This proved that no significant interferences on the adsorption properties were experienced due to the presence of other wastewater components. The performance parameters of the method developed in the present study were in fact comparable with other reported β-estradiol removal methods based on adsorption, except for References [18,57], as summarized in Table 5. Table 5. Comparison of materials used in previous reports for the removal of β-estradiol from water samples.

Regeneration Studies
The reusability and regeneration of the prepared akaganeite nanorods were evaluated by performing a series of adsorption-desorption cycles ( Figure S4). Desorption of the analyte was achieved using 100% acetonitrile. Thereafter, the adsorbent was washed with an ethanol-water solution and dried in an oven at 50 • C. The removal efficiency of the nanoadsorbent toward β-estradiol slightly decreased after the 6th usage/regeneration cycle; however, it remained at >90%. These findings suggested that the akaganeite nanoadsorbent could be regenerated and subjected to multiple usages with suitable robustness and reusability properties.

Conclusions
Akaganeite nanoparticles generated via a facile and simple synthesis route using a precipitation technique are highly useful for the adsorptive removal of β-estradiol. A variety of characterization techniques, including XRD, SEM, EDX, and TEM, confirmed the formation of akaganeite nanoparticles that are predominantly rod-shaped. The method was demonstrated to remove β-estradiol at concentrations up to 8 mg L −1 from a variety of spiked aqueous solutions with removal efficiencies of 93-100%. The adsorption process was characterized using both linear and adjusted non-linear models. The adsorption process followed pseudo-second order kinetics and a Langmuir isotherm model with a maximum adsorption capacity of 97.0 mg g −1 for the linearized isotherm model. The adjusted non-linear Langmuir isotherm resulted in a maximum adsorption capacity of 102 mg g −1 , which was comparable to the linear model. The akaganeite nanoparticles also could be reused during up to six adsorption/desorption cycles maintaining an approx. 90% removal efficiency. Hence, a low-cost and environmentally friendly adsorbent was demonstrated, which is a suitable candidate for water treatment via adsorption for hormone-like substances such as β-estradiol.