From Rice Husk Ash to Silica-Supported Carbon Nanomaterials: Characterization and Analytical Application for Pre-Concentration of Steroid Hormones from Environmental Waters

Rice husk (RH) in the rice industry is often air-burnt to obtain energy in the form of heat and RH ash (RHA) residue. In this work, RHA was applied as a starting material to obtain silica-supported carbon nanomaterials, resulting in a new reuse of a globally produced industrial waste product, in a circular economy approach. The preparation involves ultrasound-assisted one-pot oxidation with a sulfonitric mixture followed by wet oven treatment in a closed vessel. A study of oxidation times and RHA amount/acid volume ratio led to a solid material (nC-RHA@SiO2) and a solution containing silica-supported carbon quantum dots (CQD-RHA@SiO2). TEM analyses evidenced that nC-RHA@SiO2 consists of nanoparticle aggregates, while CQD-RHA@SiO2 are carbon-coated spherical silica nanoparticles. The presence of oxygenated carbon functional groups, highlighted by XPS analyses, makes these materials suitable for a wide range of analytical applications. As the main product, nC-RHA@SiO2 was tested for its affinity towards steroid hormones. Solid-phase extractions were carried out on environmental waters for the determination of target analytes at different concentrations (10, 50, and 200 ng L−1), achieving quantitative adsorption and recoveries (RSD < 20%, n = 3). The method was successfully employed for monitoring lake, river, and wastewater treatment plant water samples collected in Northern Italy.


Introduction
Rice is the second most produced cereal in the world, and its cultivation is the main activity and source of income for several countries in the world. At the beginning of the 1960s, its annual production was around 215 million tons; in 2020, its annual production was around 750 million tons [1]. After harvesting the rice, the first step in the production chain is milling. During this process, approximately 80% by weight is retained as rice, broken rice, and bran, while the remaining 20% is retained as rice husk (RH) [2]. RH is non-food biomass composed of 15-28% silica and 72-85% other organic components such as lignin, cellulose, and hemicellulose [3]. Its high silica content makes it difficult to compost or use as manure or animal feed due to its low nutritional value [2]. As a residue, RH in the rice industry is often air-burnt, obtaining, in line with the principles of circular economy, energy in the form of heat and RH ash (RHA), made of ca. 10% unburned carbon and around 90% amorphous silica [2]. Some other ways to reduce its accumulation are uses as a supplementary cementitious material, fuel in kilns to produce bricks and clay products or in rice mills to generate steam for the parboiling process, additive in fertilizers, flooring, and the production of biochar [2,4]. On the other hand, RH has great potential as a low-cost precursor to produce valuable materials, and its advantages are sustainability and environmental friendliness, despite the high diversity of the agricultural RH around the world leading to a certain variation in its properties and composition [5]. RH proves to be a good precursor for other sorbent materials as it has a granular structure, insolubility in water, chemical stability, and high mechanical strength. Adsorption is favored both by its irregular surface, which makes adsorption possible in different parts of the material, and by the presence of carboxylic functional groups and silanols that promote bonding and give considerable cation exchange capacity. It also has various cavities of different sizes, which are evidence of the porous structure and the large specific internal area that further facilitates analyte adsorption [6,7].
In the last few years, environmental pollution and the need to conserve energy and resources have led to the search for new, inexpensive, and available multipurpose sorbent materials that can be found in agro-food waste and by-products as primary sources [6]. Urbanization, industrialization, and population growth lead to the release of various types of pollutants from different sources into the environment. Natural activity contributes in small part to water pollution, but basically, the main sources are anthropogenic sources such as untreated municipal wastes and effluents, agricultural releases, and industrial chemical wastes [8]. Besides the conventional water pollutants, emerging contaminants encompass a very heterogeneous group of natural and synthetic compounds, such as pesticides, antimicrobials, steroid hormones (SHs), personal care products, and a variety of pharmacologically active drug residuals. They are defined as natural and synthetic compounds, generally detected at concentrations of nanograms per liter to micrograms per liter in water, and their presence in aqueous bodies is of particular concern because their behavior and toxicity are still largely unknown and unregulated [9].
The aim of this work is to develop a carbon nanomaterial endowed with sorption affinity for organic compounds that can be present in environmental waters as emerging contaminants. For the preparation of the sorbent, RHA provided by a local company has been chosen as a raw precursor to directly obtain silica-supported carbon nanomaterials, thus also allowing a different reuse of an industrial by-product. Silica-supported C-based materials prepared via grafting onto SiO 2 substrates, including polymers [10] and different carbon nanomaterials [11][12][13], are currently used for the extraction of contaminants from environmental, food, and biological matrices.
Besides nC-RHA@SiO 2 , the procedure leads also to a solution containing silicasupported carbon quantum dots (CQD-RHA@SiO 2 ), which were characterized (morphology, particle size, chemical composition) for future application.

Evaluation of Operational Factors for RH and RHA Conversion into Carbon Nanomaterials
The first part of this study was aimed at developing a preparation procedure for the synthesis of silica-supported carbon nanomaterials starting from RH, as described by Wang et al. [3].
In detail, 1 g of RH was pyrolyzed in a tubular furnace under nitrogen flow by following these conditions: heating ramp of 10 • C min −1 to 700 • C, isothermal condition for two hours, and cooling ramp of 10 • C min −1 to room temperature.
After verifying that pyrolyzed RH did not show any adsorption capability towards the target SHs, the char (1 g) was subjected to oxidation [3] by ultrasound-assisted one-pot treatment in 30 mL of a sulfonitric mixture (H 2 SO 4 :HNO 3 1:2, v/v) [14], and different oxidation times were investigated, specifically 2, 4, and 8 h. Then, 450 mg of the oxidized char were re-dispersed in 65 mL of tridistilled water and submitted to a thermal treatment in an oven at 200 • C for 10 h inside a 100 mL Teflon-lined autoclave (see details in Section 3.2). The solid materials obtained by filtration (0.22 µm) were investigated in terms of morphology by TEM analysis, highlighting the presence of micrometric aggregates of spherical nanoparticles with a diameter between 20 and 40 nm for all materials ( Figure S1, Supplementary Materials). To evaluate their adsorption capacity and extraction efficiency towards SHs, preconcentration tests were carried out in tap water enriched with the target analytes. In detail, 25 mg of each material packed in a 1 mL SPE column was conditioned with 5 mL MeOH, 3 mL tridistilled water, and 2 mL tap water, and then 50 mL of tap water spiked with 2 µg L −1 of each steroid was passed through. Elution was carried out with 2 × 1 mL MeOH (5% v/v FA). Quantitative analysis was performed using HPLC-MS/MS in multiple reaction monitoring (MRM) mode with quantification by standard additions on each SPE eluate. These experiments showed that most analytes were quantitatively eluted from the materials prepared at any oxidation times, but the percentage of the adsorbed analytes increased as oxidation times increased (Table 1); thus, the 8 h oxidative treatment was selected. Table 1. Adsorbed percentage (% Ads) and mean recovery of analytes in the sample (% R, n = 2) observed testing the oxidized material at different times (2 h, 4 h, and 8 h) on 50 mL of tap water spiked with 2 µg L −1 of each analyte. Moreover, SPE control tests were also carried out on RH char after oxidation without the thermal treatment in a closed vessel; these tests resulted in non-quantitative adsorption that proved the effectiveness of the thermal treatment.

Oxidation
These encouraging results led to the possibility of directly employing RHA-the bottom ash from RH combustion that is conveniently practiced by the rice industry to reduce waste and to obtain energy-as starting material for the preparation of the sorbent phase. Synthesis was carried out as described above, resulting in a solid material (nC-RHA@SiO 2 ) as the major product and a filtered solution containing CQDs grafted to silica nanoparticles (CQD-RHA@SiO 2 , fully characterized in Section 2.2).
The TEM analyses on nC-RHA@SiO 2 showed dimensions similar to those of the material obtained from pyrolyzed RH, and the pre-concentration tests (25 mg of sorbent material, 50 mL tap water enriched with 2 µg L −1 ) highlighted quantitative adsorption and recoveries (65-110%, n = 3). Given these satisfactory results, it was decided to further study RHA, as a major goal of this work. The amount of ash in the oxidation process, i.e., the liquid-to-solid phase ratio, was evaluated to reduce the acid volume per gram RHA. In detail 1, 3, and 5 g of char were oxidized with 30 mL of sulfonitric mixture for 8 h. Full adsorption and quantitative recoveries (60-110%, n = 3) were observed for all three materials obtained, and hence the liquid-to-solid phase ratio 30 mL:5 g RHA was chosen. This is a remarkable point that makes the original preparation method by Wang et al. [3] quicker and more sustainable and at the same time allows a scale-up of the synthesis, with an RHA-to-product conversion yield of 70%.
The batch-to-batch reproducibility was assessed by recovery tests on eight batches independently synthesized according to the final procedure (see Section 2.2); recovery was quantitative for all drugs, and the RSDs did not exceed 15%, thus attesting a satisfying reproducibility of sorbent preparation.
Driven by these appealing results, nC-RHA@SiO 2 was deeply characterized (Section 2.2) and applied as an SPE sorbent in environmental waters (Sections 2.3 and 2.5).

Characterization of nC-RHA@SiO 2 and CQD-RHA@SiO 2
The BET analysis of nC-RHA@SiO 2 showed a surface area of 16 m 2 g −1 . The TGA analysis ( Figure S2, Supplementary Materials) shows that from 100 • C to 1000 • C, the material is stable with a weight loss of about 2% around 500 • C, in accordance with the small amount of carbon in nC-RHA@SiO 2 evidenced by XPS analysis, reported in the following.
The zero-charge point (PZC) of nC-RHA@SiO 2 was determined using the immersion technique [15], and the PZC value was observed at~pH = 7.
nC-RHA@SiO 2 was investigated by TEM-based techniques, showing few nanometersize nanoparticles, often forming aggregates (Figure 1a,b). Compositional analysis (Figure 1c) shows the co-localization of Si, O, and C.
A similar elemental distribution is observed for CDQ-RHA@SiO 2 , where the relative amount of C is increased compared to the sample nC-RHA@SiO 2 ( Figure S3 and Table S1, Supplementary Materials), and a closer look shows that each of the particles consists of a core of oxidized silicon (≈10 nm diameter) coated by a thin (≈3 nm thick) carbonaceous layer ( Figure 2). High-resolution TEM analysis (not reported here) shows that the particles are amorphous.
XPS analysis was carried out to shed light on the chemical state of the carbon fraction of nC-RHA@SiO 2 . Analysis of the survey scan ( Figure S4, Supplementary Materials) suggests that carbon represents approximately 10 at.% of the surface composition of the nC-RHA@SiO 2 . The remaining 90 at.% is attributed to Si (approximately 27 at.%) and O (approximately 59 at.%), as well as traces of other metals (mainly Mg, Ca, Fe, and Mn). To infer the hybridization state of carbon, the X-ray-excited C KLL Auger spectrum was acquired and analyzed. As reported in the literature, the analysis of this spectrum is indeed more diagnostic of the carbon sp 2 and sp 3 configuration than the XPS C 1s spectrum, in which very similar position values are expected for the different hybridization of C [16]. The so-called D-parameter, defined as the distance in energy between the most positive maximum and most negative minimum of the first derivative of the C KLL spectrum, was therefore estimated in the case of nC-RHA@SiO 2 , resulting in a value of approx. 14.5 eV, indicative of an sp 2 content in the 12-15% range [16]. Moreover, by analyzing the XPS C 1s spectrum, we were able to identify the presence of oxygenated carbon functionalities. As shown in Figure 3, besides the main component centred at 284.8 ± 0.2 eV and ascribable to C-C and C-H bonds, three other C species are present. These are ascribable to C-O (peak at 286.1 ± 0.2 eV), C=O/O-C-O (peak at 287.5 ± 0.2 eV), and -CO 3 (peak at 287.5 ± 0.2 eV) groups [17]. A similar elemental distribution is observed for CDQ-RHA@SiO2, where the relative amount of C is increased compared to the sample nC-RHA@SiO2 ( Figure S3 and Table S1, Supplementary Material), and a closer look shows that each of the particles consists of a core of oxidized silicon (≈10 nm diameter) coated by a thin (≈3 nm thick) carbonaceous layer ( Figure 2). High-resolution TEM analysis (not reported here) shows that the particles are amorphous.   Similar XPS analyses were also conducted on CQD-RHA@SiO2. In this case, the sample is characterized by a higher carbon content (around 18 at.% of the surface composition), but also by a relatively high K content. The presence of K effectively precludes the possibility of analyzing the C KLL Auger spectrum, because of the overlapping with the K LMM Auger peak. Analysis of the XPS C 1s spectrum (not reported) revealed the presence of the same oxygenated carbon functionalities as seen in the nC-RHA@SiO2 case, even if in different relative amounts (C-C/C-H: 53.5%; C-O: 11.5%; C=O/O-C-O: 13.6%; -CO3: 21.4%).
The photoluminescence spectra (Figure 4) of the CQD-RHA@SiO2 solution acquired in the range λex 225-275 nm show a maximum emission at 425 nm with a single symmetrical band that is evidence of homogeneity of the suspended nanoparticles. The optical features in terms of incident radiation absorption and emission in the visible range are in agreement with those of CQDs prepared by valuable precursors [18,19]. Similar XPS analyses were also conducted on CQD-RHA@SiO 2 . In this case, the sample is characterized by a higher carbon content (around 18 at.% of the surface composition), but also by a relatively high K content. The presence of K effectively precludes the possibility of analyzing the C KLL Auger spectrum, because of the overlapping with the K LMM Auger peak. Analysis of the XPS C 1s spectrum (not reported) revealed the presence of the same oxygenated carbon functionalities as seen in the nC-RHA@SiO 2 case, even if in different relative amounts (C-C/C-H: 53.5%; C-O: 11.5%; C=O/O-C-O: 13.6%; -CO 3 : 21.4%).
The photoluminescence spectra (Figure 4) of the CQD-RHA@SiO 2 solution acquired in the range λ ex 225-275 nm show a maximum emission at 425 nm with a single symmetrical band that is evidence of homogeneity of the suspended nanoparticles. The optical features in terms of incident radiation absorption and emission in the visible range are in agreement with those of CQDs prepared by valuable precursors [18,19].

Analytical Application of nC-RHA@SiO 2 as SPE Sorbent in Environmental Waters
As described in Section 2.1, the explorative SPE experiments highlighted that adsorption and recoveries were quantitative on nC-RHA@SiO 2 , well supporting the use of an agro-food by-product as RHA (instead of RH char from pyrolysis) as the sorbent precursor.
After evaluating the key role of the carbon aggregates in terms of adsorption capability by control test on pure silica (Ads < 40%), SPE tests on nC-RHA@SiO 2 were undertaken to evaluate the role of different parameters affecting the SPE process. Firstly, the sample volume was increased up to 100 mL, but semi-quantitative adsorption was observed using 25 mg sorbent. Working with 50 mg nC-RHA@SiO 2 , quantitative adsorption and recovery (R > 75%, RSD < 15%, n = 3) were achieved.

Analytical Application of nC-RHA@SiO2 as SPE Sorbent in Environmental Waters
As described in Section 2.1, the explorative SPE experiments highlighted that adsorption and recoveries were quantitative on nC-RHA@SiO2, well supporting the use of an agro-food by-product as RHA (instead of RH char from pyrolysis) as the sorbent precursor.
After evaluating the key role of the carbon aggregates in terms of adsorption capability by control test on pure silica (Ads < 40%), SPE tests on nC-RHA@SiO2 were undertaken to evaluate the role of different parameters affecting the SPE process. Firstly, the sample volume was increased up to 100 mL, but semi-quantitative adsorption was observed using 25 mg sorbent. Working with 50 mg nC-RHA@SiO2, quantitative adsorption and recovery (R > 75%, RSD < 15%, n = 3) were achieved.
Elution, as another SPE key factor, was then investigated by evaluating the percent acidity of MeOH. For this, MeOH, MeOH (1% v/v FA), and MeOH (5% v/v FA) were tested (2 × 1 mL). Lower recoveries were obtained with pure MeOH compared to acidified MeOH, and no significant variations were found between the two concentrations of FA; thus, 1% FA was selected. The higher elution efficiency under acidic conditions compared to neat MeOH can be justified considering the sorbent PZC (~pH = 7), which resulted in a positively charged surface during acidic elution promoting the desorption of the SHs. Since these are neutral compounds, as it is well known, their affinity for the material is favoured when the sorbent surface is globally neutral (not charged), whereas it is reasonably weakened when the sorbent possesses a neat charge.
This outcome is also in agreement with the good adsorption observed in tap water at the native pH, which is around 7.5 with no need for pH adjustment.
Under neutral conditions, the sorption of SHs is driven by interaction with the carbon fraction of the sorbent, as experimentally evidenced by SPE control tests on bare silica (<40%). Besides apolar (hydrophobic) interactions, other interplays with the oxygenated Elution, as another SPE key factor, was then investigated by evaluating the percent acidity of MeOH. For this, MeOH, MeOH (1% v/v FA), and MeOH (5% v/v FA) were tested (2 × 1 mL). Lower recoveries were obtained with pure MeOH compared to acidified MeOH, and no significant variations were found between the two concentrations of FA; thus, 1% FA was selected. The higher elution efficiency under acidic conditions compared to neat MeOH can be justified considering the sorbent PZC (~pH = 7), which resulted in a positively charged surface during acidic elution promoting the desorption of the SHs. Since these are neutral compounds, as it is well known, their affinity for the material is favoured when the sorbent surface is globally neutral (not charged), whereas it is reasonably weakened when the sorbent possesses a neat charge.
This outcome is also in agreement with the good adsorption observed in tap water at the native pH, which is around 7.5 with no need for pH adjustment.
Under neutral conditions, the sorption of SHs is driven by interaction with the carbon fraction of the sorbent, as experimentally evidenced by SPE control tests on bare silica (<40%). Besides apolar (hydrophobic) interactions, other interplays with the oxygenated functionalities of the carbon fraction, confirmed by XPS, and the analyte molecular structures cannot be excluded ( Figure S4, Supplementary Materials).
Guided by these excellent results, the final procedure (50 mg sorbent, 100 mL sample, elution with 2 × 1 mL 1% v/v FA MeOH) was extended to different environmental waters such as lake, river, and UWWTP effluent waters. Recovery tests were carried out in triplicate at three concentration levels (200 ng L −1 , 50 ng L −1 , and 10 ng L −1 ), and unspiked samples were simultaneously analyzed, showing no peaks at the retention time of the target analytes. For the 10 ng L −1 and 50 ng L −1 concentrations, the eluates were evaporated (see Section 3.4). As shown in Table 2, satisfactory recoveries were obtained for all the different matrices also at the lowest concentration (10 ng L −1 ), with an overall enrichment factor up to 200, thus strengthening the applicability of the material for extraction of trace pollutants in real-contaminated waters. Good results are achieved by a small amount of sorbent, with operation in the micro-SPE mode [20] further reducing solvent use. A representative MRM chromatogram from the analysis of the SPE lake eluates is presented in Figure 5.
Driven by these satisfactory results and the inherent characteristics of the developed sorbent, its applicability will be extended to the pre-concentration of other classes of emerging organic pollutants, e.g., fluoroquinolone antimicrobials. Additional work is in progress to enhance the sample throughput and to increase the sustainability of the procedure, further minimizing waste generation.

Analytical Evaluation of the Method
The entire analytical procedure was within-laboratory evaluated, according to the key features described in detail in Section 3.7.
Selectivity was guaranteed by MRM detection, which affords identification/quantification of the target compounds by using the most intense transition precursor/product ions of each compound (Table S2, Supplementary Materials). Furthermore, analysis of the eluates of the samples without spikes was carried out, showing that the target analytes of this method were not present.
The matrix-matched calibration, for quantitation of the concentrations expected, was performed in the range 1-50 µg L −1 by five independent calibration curves in the eluate from blank environmental waters and provided good linearity (R 2 > 0.9994 for lake sample, R 2 > 0.9992 for river and UWWTP effluent samples). Table 3 reports the method detection limits (MDLs) and method quantification limits (MQLs) obtained from each matrix-matched calibration curve for all the analytes, taking into consideration the whole sample treatment (see Section 3.7). Experimental MQL (10 ng L −1 ) was verified for all analytes in environmental waters.   The presence of matrix effect (ME) was assessed by comparing the slopes of matrixmatched calibration curves with those of the calibration lines obtained in the pure solvent, achieving the results reported in Table 3, and it was considered significative when not included within the method precision (RSD < 20%) [21].

Determination of SHs in Environmental Waters
The final procedure was used to investigate the presence of target SHs in environmental water samples, viz. from Garda Lake, Ticino River, and the wastewater effluent downstream the Vigevano urban wastewater treatment plant. Most analytes were detected (3-10 ng L −1 ) or quantified at concentrations in the range 10-23 ng L −1 (RSD < 20%, n = 3) (see Table 4). Table 4. Nanograms per liter of SHs quantified in Garda Lake, Ticino River, and Vigevano UWWTP effluent water samples.

Materials Preparation
Five grams of rice husk ash (RHA) was sieved (100 mesh) and subjected to an oxidative treatment in an ultrasonic bath (Ultrasonic cleaner, 45 kHz, 80 W USC 200-2600, VWR International, Milan, Italy) with 30 mL of a sulfonitric mixture (H 2 SO 4 :HNO 3 , 1:2, v/v) for 8 h. The obtained solid char was vacuum filtered using a nylon filter (0.22 µm) and washed with plenty of tridistilled water (product yield 80% wt%). Four hundred fifty milligrams of the obtained material dispersed in 65 mL of tridistilled water was subjected to a thermal treatment in a high-pressure Teflon vessel (Acid Digestion Bomb with Teflon liner, bomb volume 23 mL, stainless steel body, Parr Instruments, Moline, IL, USA) at 200 • C (oven) for 10 h. The obtained suspension was then filtered through a 0.22 µm nylon filter and CQD-RHA@SiO 2 were collected in the filtrate, while the final material used as sorbent, nC-RHA@SiO 2 , was dried at 60 • C for 2 h (yield 85% wt%). The overall yield of RHA conversion into nC-RHA@SiO 2 was 70 wt%.

Silica-Supported Carbon Nanomaterial Characterization
Specific surface area measurements on the samples were performed using a Sorptomatic 1990 instrument produced by Thermo Fisher Scientific. About 350 mg of powder was charged in the glass sample holder and degassed at 25 • C for 24 h. Subsequently, samples were cooled down at −196 • C and 2 adsorption-desorption cycles followed by a last adsorption run were performed (BET method, analyzing gas N 2 , 50 points for run; blank performed in He before the first adsorption run).
TGA measurements were performed with a TGA 1 Star System from Mettler Toledo (Novate Milanese, Italy). The heating rate was 10 • C min −1 , and the range of temperature was from 25 • C to 1000 • C. Reactive gas was air (Sapio, Monza, Italy) at a flow rate of 100 mL min −1 . For each experiment, about 10 mg of material was weighed and put into an alumina crucible. The TGA furnace was conditioned for 2 h to provide a suitable atmosphere for combustion.
Overview BF-TEM images were recorded on a Jeol JEM-1200 EX II instrument; a small amount of material (ca. 5 mg) was dispersed in water (ca. 1 mL) and sonicated for a few minutes, and a certain amount (ca. 5 µL) of the suspension was taken and deposited on the grid (400 mesh) through a micropipette.
Energy-filtered TEM (EF-TEM) and high-angle annular dark-field scanning TEM (HAADF-STEM) images and energy-dispersive X-ray spectroscopy data were acquired on an image-Cs-corrected Jeol JEM-2200FS S/TEM equipped with a Bruker X-Flash 5060 silicon drift detector-based EDS system.
The high silica content of CQD-RHA@SiO 2 obtained from RH, evidenced by TEM analysis, was reduced by a NaOH treatment [3] and further characterized as follows: the CQD-RHA@SiO 2 suspensions were drop-cast onto Cu grids coated with a holey carbon film. EDS mapping reported here was obtained by integration of intensities on the Kα peaks of C, O, and Si, and quantification was performed using the Cliff-Lorimer method and the same peaks. EF-TEM maps were obtained using the three-window method at the K core-loss edges of C and O, and the L 23 edge of Si.
X-ray photoelectron spectroscopy (XPS) analyses were conducted with a Kratos Axis Ultra DLD spectrometer, using a monochromatic Al Kα source operated at 15 kV and 20 mA. Survey scans were acquired with an analysis area of 300 × 700 microns and a pass energy of 160 eV. High-resolution analyses were carried out with the same analysis area and a pass energy of 20 eV. Spectra were charge-corrected to the main line of the carbon 1s spectrum set to 284.8 eV. Spectra were analyzed using CasaXPS software (version 2.3.24).

Environmental Waters
Environmental waters were collected directly in amber glass bottles in the spring and summer of 2020 at 30-50 cm depth from Garda Lake, Ticino River, and UWWTP effluent in Vigevano (Pavia, Italy). Samples were stored in the dark (4 • C) and the physical-chemical parameters were measured for each sample, as reported in Table S3, Supplementary Materials.

HPLC-ESI-MS/MS
The HPLC apparatus Agilent 1260 Infinity coupled with an Agilent 6460C MS spectrometer ESI-MS/MS system (Cernusco sul Naviglio, Italy) was used for the analysis of the target analytes, according to previous papers [11,22]. Chromatographic separation was achieved with ZORBAX Eclipse Plus C18 Rapid Resolution (100 mm × 4.6 mm, 3.5 µm) preceded by a pre-column Supelco Supelguard Ascentis C18 (2 cm × 2.1 mm, 5 µm), both maintained at 25 • C (±0.8 • C). The sample injection volume was 10 µL, robotically injected by the HPLC autosampler. The flow rate was 0.5 mL min −1 , with an analysis time equal to 20 min and a re-equilibration time of 7 min. Elution was performed with (A) aqueous 1 mM NH 4 F in H 2 O and (B) ACN, according to the following elution program: 30% B for 4 min, linear gradient to 85% B in 8 min, maintained for 5 min, and to 98% B in 0.5 min, maintained for 2 min. The MS operating parameters optimized by Agilent Mass Hunter Source Optimizer Software (Agilent, Santa Clara, CA, USA) were as follows: drying gas (N 2 ) temperature 350 • C; drying gas flow 12 L min −1 ; nebulizer 50 psi; sheath gas temperature 400 • C; sheath gas flow 12 L min −1 ; capillary voltage 4000 V positive and 3000 V negative; nozzle voltage 0 V positive and 1500 V negative; electron multiplier voltage (EMV) 200 V positive and 0 V negative; cell accelerated voltage (CAV) 4 V positive and 1 V negative. Quantitative analysis was performed in multiple-reaction monitoring (MRM) mode, using the two most intense transitions of each compound (see Table S2, Supplementary Materials). All analytes were detected in the same chromatographic run by setting a polarity-switching tool. Typical MRM chromatograms of methanol standard solutions (10 µg L −1 of each analyte) are presented in Figure S6, Supplementary Materials.

Analytical Evaluation of the Method
Selectivity was evaluated by the analysis of the blank environmental water with and without analyte spikes to evaluate the role of co-extracted components from the matrix.
To assess linearity in the matrix elutes, calibration curves were prepared for each analyte in pure solvent (MeOH) and in the elutes from each blank environmental water sample (matrix-matched calibration curves) [23,24] in the concentration range of 1-50 µg L −1 . In ESI, the ionization efficiency of the analytes may be strongly altered by the co-eluting compounds, and different matrices can lead to different matrix effects; thus, ME was investigated in each extract from lake, river, and wastewater effluent samples, and it was calculated as follows: where "b m " and "b s " are the slopes of the matrix-matched calibration curve and of the calibration line obtained in pure solvent, respectively, obtained by ordinary linear lowest squares regression (OLLSR). MDLs and MQLs (µg L −1 ) were calculated from the matrix-matched calibration curves as 3.3 and 10 times the ratio between the baseline noise away from the peak tail and the regression line slope, respectively, taking into consideration the whole sample treatment [11,22]. Independent recovery tests (n = 3) were performed on environmental water samples enriched with each analyte at three different concentration levels, namely low-quality control (LQC, 10 ng L −1 ), medium-quality control (MQC, 50 ng L −1 ), and high-quality control (HQC, 200 ng L −1 ), to evaluate trueness and precision. Trueness was calculated as the ratio of the concentration obtained after the entire procedure on the spiked sample and the concentration expected. The inter-day precision was expressed as RSD% associated with the mean recovery.
Robustness was evaluated by small and deliberate variations in operating conditions, whereas instrumental carry-over was checked by injections of SPE eluting solution (MeOH 1% v/v FA) and SPE eluates from each blank water sample as control blanks after consecutively analyzed samples [22].

Conclusions
RHA, obtained in the rice industry as a bottom ash of the RH combustion, has been successfully transformed into a working fixed-bed SPE sorbent (nC-RHA@SiO 2 ) for the preconcentration of emerging and hazardous pollutants, such as SHs. Indeed, the presence of carbon embedding oxygenated groups, shown by XPS and HAADF-STEM analyses, exerts a key role in the quantitative adsorption/desorption of the target SHs in natural waters.
A simple and quick analytical method, based on a SPE procedure followed by HPLC-MS/MS quantification, was developed for the multiclass determination of SHs at environmentally relevant concentrations.
Another important goal achieved in this work is the one-pot treatment involved in the RHA oxidation, with a remarkable reduction in acid consumption. Noteworthily, this procedure also led to added-value products such as CQD-RHA@SiO 2 which will be tested in future work.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules28020745/s1, Figure S1: TEM images at different magnifications (25,000, 200,000, and 250,000) of the material after oven treatment from pyrolyzed RH oxidized at different times: (a) 2 h, (b) 4 h, and (c) 8 h; Figure S2: TGA profiles recorded on nC-RHA@SiO 2 ; Figure S3: (a) HAADF-STEM image of a region in the sample CQD-RHA@SiO 2 , and corresponding STEM-EDS maps for Si (b), C (c) and O (d). For these regions, partly suspended on holes in the support carbon film, there is a clear co-localization of the three elements in the nanostructured material; Figure S4: XPS survey spectrum acquired on nC-RHA@SiO 2 ; Figure S5: Molecular structures and LogP values of the target analytes; Figure S6: Typical MRM chromatogram of a SH standard solution (10 µg L −1 of each analyte) prepared in methanol.; Table S1: STEM-EDS compositional analysis of (a) nC-RHA@SiO 2 and (b) CQD-RHA@SiO 2 suspended on holes in the carbon support film; Table S2: MRM conditions for the HPLC-ESI-MS/MS analysis of SHs; Table S3: Physicochemical parameters of the water samples.  Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.