**1. Introduction**

Carbonaceous materials constitute an attractive group of adsorbents due to their properties, such as large surface area and thermal and chemical stability. In addition, activated carbons with different pore sizes can be obtained, thus adapting them to the functions they are designed for. On the other hand, functionalizing the surface of activated carbons is possible, giving them specific chemical properties adapted to a particular purpose. It is also of great importance that carbonaceous materials are considered safe, accessible, and affordable [1]. Various types of biomass residues, due to their low cost and abundance, have been evaluated to obtain functional carbonaceous materials. Many methods provide materials whose physicochemical properties can be tuned at will, thus obtaining materials that have shown promising results in different fields, such as adsorption or energy storage [2]. For this reason, they are considered suitable candidates for obtaining sustainable materials. These materials have been successfully applied in catalysis, water

**Citation:** Izquierdo, S.; Pacheco, N.; Durán-Valle, C.J.; López-Coca, I.M. From Waste to Resource: Utilizing Sweet Chestnut Waste to Produce Hydrothermal Carbon for Water Decontamination. *C* **2023**, *9*, 57. https://doi.org/10.3390/c9020057

Academic Editors: Indra Neel Pulidindi, Pankaj Sharma, Aharon Gedanken and Shu-Yuan Pan

Received: 30 March 2023 Revised: 8 May 2023 Accepted: 29 May 2023 Published: 1 June 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

purification, and soil remediation since they can remove heavy metals, agrochemicals, antibiotics, and other organic compounds. Excellent reviews on these aspects can be found in Soffian et al. [3], Adegoke et al. [4], G˛eca et al. [5], Borchardt et al. [6], and Ma et al. [2]. One such type of carbonaceous material is biochar, which can be produced from biomass waste through hydrothermal carbonization, chemical carbonization, gasification, torrefaction, and pyrolysis [3,7–9].

Sweet chestnut (*Castanea sativa*) is a deciduous tree of the Fagaceae family, native to humid temperate regions of the Northern Hemisphere and widely distributed throughout Europe. It can reach 35 m in height and 2 m (exceptionally more) in diameter. Sweet chestnut trees live to an age of 600 years or even more. By autumn, the female flowers develop into spiny cupules containing 3–7 brownish nuts, which are shed in October [10–12]. After the chestnuts are harvested, the spiny dome is disposed of as waste and used as fuel for fireplaces.

On the other hand, per- and poly-fluoroalkyl substances (PFASs) are a group of synthetic chemicals that have been used in a wide range of industrial and consumer products, such as firefighting foams, non-stick coatings, and water-repellent fabrics [13,14]. Due to their persistence, bioaccumulation, and toxicity, PFASs are considered Contaminants of Emerging Concern (CECs), with potentially adverse effects on human health and the environment [15–21]. These substances have been detected in various environmental media, including water, soil, air, and biota, and have attracted increasing attention from regulators, scientists, and the public worldwide [22]. Among the different methods for treating PFAS-contaminated water, adsorption has been recognized as a promising technology, given its effectiveness, versatility, and ease of implementation [23,24]. Carbonaceous materials are widely used as adsorbents for PFAS removal due to their high surface area, microporous structure, and affinity for hydrophobic compounds. One of the classic and best-known carbonaceous materials is activated carbon (AC) [25–27]. It is usually obtained at high temperatures in a double carbonization and activation process. This material is characterized by a highly developed porous structure, a generally alkaline character, and a predominantly hydrophobic surface [28,29]. It is well known that these properties can be modified by adjusting the conditions of the activation processes or employing suitable chemical or thermal treatment [30,31].

In recent years, numerous studies have investigated using these adsorbents for PFAS adsorption in water, aiming to optimize the adsorbent properties and conditions, understand the adsorption mechanisms, and evaluate the performance and feasibility of the technology. The effectiveness of AC as an adsorbent for PFAS removal depends on its physicochemical properties, which can be tailored by adjusting the activation process and conditions. The surface area, pore size, surface chemistry, and surface charge of AC determine its adsorption capacity, selectivity, and kinetics for PFAS removal. The larger the surface area and the more mesopores present, the higher the adsorption capacity; this is due to the increased availability of adsorption sites and the enhanced accessibility of PFAS molecules to the internal surface of the carbonaceous material. The surface chemistry can also affect PFAS adsorption, as hydrophobic interactions between PFASs and the nonpolar surface are favored over electrostatic interactions with the polar surface of water molecules. Therefore, adsorbents with a higher content of surface functional groups, such as oxygen-containing or nitrogen-containing groups, may have a lower adsorption capacity for PFASs than non-functionalized ones due to the reduced hydrophobicity and the increased competition for adsorption sites.

In the last decade or so, a new type of carbonaceous material, hydrothermal carbons (HC), has attracted considerable interest from researchers in this field [32,33]. It is obtained by treating carbon-rich materials with water at a moderate temperature in a closed vessel so that endogenous water vapor pressure acts on this material [34]. Generally, some biomass is used as raw material. This type of material is characterized by low pore development, acidity, and a hydrophilic surface, with many surface groups. As in the previous case, these properties can be modified by thermal or chemical treatments. This

has the advantage of being a more environmentally friendly manufacturing process, as it is carried out in hermetically sealed containers and at a lower temperature. This allows for easier management of possible waste and lower energy consumption. Additionally, despite its limited specific surface area, its high proportion of functional groups can give rise to adequate adsorption capacity. We have been working on the obtention and use of hydrothermal carbons from chestnut cupules for the last few years. Our preliminary results (unpublished) show that these carbons are suitable adsorbents in the decontamination of heavy metals, such as thorium, from water. Given these results, we thought it would be apt to check for the adsorption capacity of apolar organic molecules such as PFAS. Adsorption onto carbonaceous materials has been established as a promising technology for PFAS removal from water [20,35–42]. as it is carried out in hermetically sealed containers and at a lower temperature. This allows for easier management of possible waste and lower energy consumption. Additionally, despite its limited specific surface area, its high proportion of functional groups can give rise to adequate adsorption capacity. We have been working on the obtention and use of hydrothermal carbons from chestnut cupules for the last few years. Our preliminary results (unpublished) show that these carbons are suitable adsorbents in the decontamination of heavy metals, such as thorium, from water. Given these results, we thought it would be apt to check for the adsorption capacity of apolar organic molecules such as PFAS. Adsorption onto carbonaceous materials has been established as a promising technology for PFAS removal from water [20,35–42]. The objective of this work is to investigate the use of chestnut cupule waste as sustainable biomass raw material to obtain hydrothermal carbons—hence, adding value to

This has the advantage of being a more environmentally friendly manufacturing process,

*C* **2023**, *9*, x FOR PEER REVIEW 3 of 18

The objective of this work is to investigate the use of chestnut cupule waste as sustainable biomass raw material to obtain hydrothermal carbons—hence, adding value to it—and its use as adsorbent material for the removal from water of sodium perfluorooctanoate (PFONa), the anion of perfluorooctanoic acid, one of the most significant PFASs. it—and its use as adsorbent material for the removal from water of sodium perfluorooctanoate (PFONa), the anion of perfluorooctanoic acid, one of the most significant PFASs. **2. Materials and Methods** 

#### **2. Materials and Methods** *2.1. Materials*

#### *2.1. Materials* Sodium perfluorooctanoate 97%, CAS 335-95-5 (Figure 1), was purchased from

Sodium perfluorooctanoate 97%, CAS 335-95-5 (Figure 1), was purchased from Thermo Fisher. Milli-Q water was used as a solvent in all determinations. Sweet chestnut (*Castanea sativa*) cupules were collected from chestnut groves in Valle del Jerte, in Cáceres, Spain. Thermo Fisher. Milli-Q water was used as a solvent in all determinations. Sweet chestnut (*Castanea sativa*) cupules were collected from chestnut groves in Valle del Jerte, in Cáceres, Spain.

**Figure 1.** Structural formula of sodium perfluorooctanoate. **Figure 1.** Structural formula of sodium perfluorooctanoate.

#### *2.2. Methods 2.2. Methods*

2.2.1. Synthesis of the Hydrothermal Carbons 2.2.1. Synthesis of the Hydrothermal Carbons

The chestnut cupules were allowed to dry naturally. They were subsequently ground and sieved, collecting three fractions: large size (1–2 mm), medium size (0.5–1 mm), and small size (less than 0.5 mm). Then, they were subjected to hydrothermal treatment in a The chestnut cupules were allowed to dry naturally. They were subsequently ground and sieved, collecting three fractions: large size (1–2 mm), medium size (0.5–1 mm), and small size (less than 0.5 mm). Then, they were subjected to hydrothermal treatment in a duralumin-coated Teflon™ reactor.

duralumin-coated Teflon™ reactor. To obtain HC-A, in a hydrothermal reactor, 100 mL of distilled water and 15 g of small-size sieved cupule were added. The reactor was closed, placed in an oven previously programmed at 180 °C, and kept for 24 h. The reactor was then removed from the oven and, once cold, opened. The final product was filtered and washed with plenty of distilled water (ca. 200 mL) on a filter. Subsequently, it was dried in an oven at 110 °C for To obtain HC-A, in a hydrothermal reactor, 100 mL of distilled water and 15 g of small-size sieved cupule were added. The reactor was closed, placed in an oven previously programmed at 180 ◦C, and kept for 24 h. The reactor was then removed from the oven and, once cold, opened. The final product was filtered and washed with plenty of distilled water (ca. 200 mL) on a filter. Subsequently, it was dried in an oven at 110 ◦C for 8 h. The other carbons were prepared similarly, according to the data in Table 1.

8 h. The other carbons were prepared similarly, according to the data in Table 1.

HC-C 1–2 15 200

**Table 1.** Operating conditions for the preparation of hydrothermal carbons.


#### 2.2.2. Characterization of the Hydrothermal Carbons

2.2.2. Characterization of the Hydrothermal Carbons Textural characterization was performed by N2 adsorption isotherms at 77 K in a Quantachrome Quadrasorb Evo (Quantachrome Ltd., Hook, UK), and the specific surface Textural characterization was performed by N<sup>2</sup> adsorption isotherms at 77 K in a Quantachrome Quadrasorb Evo (Quantachrome Ltd., Hook, UK), and the specific surface area was calculated by applying the BET method [43]. Porosimetry was performed in a Quantachrome Poremaster 60 (Quantachrome, Hook, UK). Thermal analysis was carried out in an STA 449 F3 Jupiter thermobalance (Netzsch, Selb, Bavaria, Germany).

Fourier-transform infrared spectra were recorded in a Bruker Vertex70 spectrometer in the 400–4000 cm−<sup>1</sup> range, with a DLaTGS detector and a 0.4 cm−<sup>1</sup> spectral resolution. A KBr tablet was prepared with 1/250 dilution (Bruker, Billerica, Massachusetts, USA). Elemental analysis (C, H, N, S, O) was performed using a Leco CHMS-932 elemental analyzer (Leco, St. Joseph, MI, USA). C, H, N, and S were analyzed, and the difference was assigned to ash (measured from proximate analysis) and oxygen content. The point of zero charge (PZC) values were determined using the method proposed by Valente Nabais and Carrott [44]. To carry out this measurement, a 0.1 M solution of sodium nitrate was prepared, and 7% in weight of the carbon was added, followed by stirring for 48 h at 25 ◦C in a thermostatic bath. Subsequently, it was filtered, and the pH of the filtered solution was measured with a pH meter. This pH value corresponded to the point of zero charge. A field emission scanning electron microscope (SEM) Quanta 3D FEG (FEI Company, Hillsboro, OR, USA) was used to explore the surface morphological characteristics of all samples. The sample analysis was performed under high vacuum with a secondary electron. The surface of the samples was coated with gold due to their low conductivity. Energy dispersive X-ray analyses (EDX) were recorded with the same equipment.

#### 2.2.3. Adsorption Experiments

The conjugate base sodium salt (PFONa), instead of the perfluorooctanoic acid, was used to carry out the adsorption experiments, since the acid has a K<sup>a</sup> ca. -1 and dissociates completely in water.

1. Adsorption kinetics

To determine the PFONa adsorption over time, the following procedure was used:

A 200 mg/L PFONa water solution was prepared by dissolving 10 mg of PFONa in a 50 mL volumetric flask. It was then transferred to an Erlenmeyer flask, and 50 mg of hydrothermal carbon (HC) was added, so the HC to PFONa ratio was 5:1. All flasks were thermostated at 25 ◦C and magnetically stirred at 400 rpm. Assays were performed in duplicate; 2 mL samples were taken at 2, 4, 24, 48, 72, and 96 h, then filtered with a 0.22 µm pore filter, and, finally, aliquots were analyzed for PFONa concentration.

#### 2. Adsorption isotherm

Different PFONa solutions of 40, 70, 100, 130, 165, and 200 mg/L were prepared in 10 mL volumetric flasks. They were then transferred to separate Erlenmeyer flasks, and 10 mg of hydrothermal carbon was added to each one, so the HC to PFONa mass quotient was 25.0, 14.3, 10.0, 7.7, 6.1, and 5.0, respectively. All flasks were thermostated at 25 ◦C and magnetically stirred at 400 rpm. Assays were performed in duplicate; 2 mL samples were taken at 96 h, then filtered with a 0.22 µm pore filter, and, finally, aliquots were analyzed for PFONa concentration.

#### 2.2.4. Analytical Determinations

The samples were analyzed using liquid chromatography (Agilent 1290 Infinity II LC (Agilent Technologies, Santa Clara, California, USA)) coupled to a triple quadrupole mass spectrometer (Agilent 6460 triple quad LC/MS (Agilent Technologies)). Aliquots of 3 µL were injected, and elution was carried out by operating in gradient mode at 0.4 mL/min with a reverse phase HPLC Zorbax Extend C18 column (3 µm, 100 mm × 1.8 mm) at 30 ◦C. As the mobile phase, 0.1% formic acid in ultrapure water (A) and 0.1% formic acid in acetonitrile (B) were used, as shown in Table 2.

Ionization was carried out in negative mode. The equipment parameters were gas temperature 225 ◦C, gas flow 10 L/min, nebulizer 45 psi, sheath gas temperature 350 ◦C, sheath gas flow 11 L/min, capillary voltage 3.6 kV, nozzle voltage 1500 V, and delta EMV 400. The MassHunter B.07.00 (Agilent) software was used for data acquisition and equipment control. Multiple Reaction Monitoring (MRM) was used, with the transitions sought to be 412.9 → 368.9 and 412.9 → 169, the former used in the quantification of the compound. The fragmentation and collision energy voltages were 86 V and 5 V, respectively.


**Table 2.** Mobile phase composition.

#### **3. Results and Discussion**

#### *3.1. Characterization of the Hydrothermal Carbons*

The nitrogen adsorption isotherms (Figure 2) corresponded to type III, according to the IUPAC classification [45,46]; therefore, the interaction between the N<sup>2</sup> molecules and the surface must have been weak. The surface obtained by the BET method produced values that were higher than those obtained by the DFT method but in the same order of magnitude (Table 3) [43,47]. Nevertheless, since the isotherms were of Type III, the BET model was not to be considered highly reliable. All three carbons showed a minimal surface area, low microporosity, and a considerably greater volume of large pores, i.e., mesopores and macropores (Figures 3–5). The HC-A coal had the least-developed porosity. The DFT method was used to study the distribution of micropores and mesopores. The results show that the three carbons presented a similar pore distribution. A high volume of narrow mesopores was detected, whereas no micropores were present. The HC-B carbon had a larger volume of narrow mesopores, while the HC-C showed more volume when the pores widened. The HC-A carbon showed a more significant development of wide porosity than the other two adsorbents. It should be noted that HC-B and HC-C carbons showed a more developed narrow mesoporosity than HC-A. However, the mesopore and macropore structure was very similar in the three materials. *C* **2023**, *9*, x FOR PEER REVIEW 6 of 18

**Figure 2.** Isotherms of N2 adsorption at 77 K. **Figure 2.** Isotherms of N<sup>2</sup> adsorption at 77 K.

Functional Theory [47].

**Material BET 1 DR 2 Porosimetry DFT 3** 

porosimetry, is shown in Figure 5.

HC-B 53.2 0.007 0.061 39.13 0.122 HC-B 53.2 0.007 0.061 39.13 0.122

**Figure 3.** Distribution of porosity of the hydrothermal carbons.

1 Brunauer, Emmett, and Teller model [43]. 2 Dubinin and Radushkevich model [48,49]. 3 Density

The DFT method was used to study the distribution of porosity in micro and mesopores, and the results are shown in Figure 3. The cumulative volume plot obtained from porosimetry is shown in Figure 4, and the width pore distribution plot, also obtained from

**Table 3.** Porosity and specific surface area of carbons.

**Material BET <sup>1</sup> DR <sup>2</sup> Porosimetry DFT <sup>3</sup> S, m<sup>2</sup> g**−**<sup>1</sup> Vmicro, cm<sup>3</sup> g**−**<sup>1</sup> Vmeso, cm<sup>3</sup> g**−**<sup>1</sup> S, m<sup>2</sup> g**−**<sup>1</sup> Vtotal, cm<sup>3</sup> g**−**<sup>1</sup>** HC-A 42.3 0.007 0.078 27.68 0.085 HC-B 53.2 0.007 0.061 39.13 0.122 HC-B 53.2 0.007 0.061 39.13 0.122 1 Brunauer, Emmett, and Teller model [43]. 2 Dubinin and Radushkevich model [48,49]. 3 Density Functional Theory [47]. The DFT method was used to study the distribution of porosity in micro and mesopores, and the results are shown in Figure 3. The cumulative volume plot obtained from

*C* **2023**, *9*, x FOR PEER REVIEW 6 of 18

**Table 3.** Porosity and specific surface area of carbons. HC-B 53.2 0.007 0.061 39.13 0.122

**Figure 2.** Isotherms of N2 adsorption at 77 K.

**Material BET 1 DR 2 Porosimetry DFT 3** 

**Table 3.** Porosity and specific surface area of carbons.

HC-B 53.2 0.007 0.061 39.13 0.122

<sup>1</sup> Brunauer, Emmett, and Teller model [43]. <sup>2</sup> Dubinin and Radushkevich model [48,49]. <sup>3</sup> Density Functional Theory [47]. porosimetry is shown in Figure 4, and the width pore distribution plot, also obtained from porosimetry, is shown in Figure 5.

**S, m2 g−1 Vmicro, cm3 g−1 Vmeso, cm3 g−1 S, m2 g−1 Vtotal, cm3 g−<sup>1</sup>**

**Figure 3. Figure 3.**  Distribution of porosity of the hydrothermal carbons. Distribution of porosity of the hydrothermal carbons.

**Figure 4.** Cumulative pore volume of the hydrothermal carbons. **Figure 4.** Cumulative pore volume of the hydrothermal carbons.

**Figure 5.** Pore width distribution of the hydrothermal carbons.

tribution of basicity by the aromatic structures.

**Table 4.** Point of zero charge of the carbons.

The measurements of the PZC show that all three materials were acidic (Table 4), as

**Material PZC**  HC-A 4.0 HC-B 3.8 HC-C 4.8

The Fourier-transform infrared (FTIR) spectra obtained are shown in Figure 6. Nu-

merous observable bands in the FTIR spectra belonging to many functional groups indicate that the degree of graphitization of these materials was not high. The most intense band at 3400 cm−1 corresponded to either alcohol or carboxylic O-H bonds. The band at 2935 cm−1 was assigned to C-H bonds, but its low wavenumber suggested they were

is usual in hydrothermal carbons. It must be considered that basic oxygenated functional groups (carbonyl) are more labile than those that provide acidic properties (alcohol and carboxylic), so it is to be expected that the former will degrade in a more significant proportion than the latter. In addition, as the degree of graphitization was low, as observed from the elemental composition with low carbon content, there was not a substantial con-

**Figure 4.** Cumulative pore volume of the hydrothermal carbons.

The measurements of the PZC show that all three materials were acidic (Table 4), as is usual in hydrothermal carbons. It must be considered that basic oxygenated functional groups (carbonyl) are more labile than those that provide acidic properties (alcohol and The DFT method was used to study the distribution of porosity in micro and mesopores, and the results are shown in Figure 3. The cumulative volume plot obtained from porosimetry is shown in Figure 4, and the width pore distribution plot, also obtained from porosimetry, is shown in Figure 5.

carboxylic), so it is to be expected that the former will degrade in a more significant proportion than the latter. In addition, as the degree of graphitization was low, as observed from the elemental composition with low carbon content, there was not a substantial contribution of basicity by the aromatic structures. **Table 4.** Point of zero charge of the carbons. The measurements of the PZC show that all three materials were acidic (Table 4), as is usual in hydrothermal carbons. It must be considered that basic oxygenated functional groups (carbonyl) are more labile than those that provide acidic properties (alcohol and carboxylic), so it is to be expected that the former will degrade in a more significant proportion than the latter. In addition, as the degree of graphitization was low, as observed from the elemental composition with low carbon content, there was not a substantial contribution of basicity by the aromatic structures.

**Material PZC Table 4.** Point of zero charge of the carbons.


merous observable bands in the FTIR spectra belonging to many functional groups indicate that the degree of graphitization of these materials was not high. The most intense band at 3400 cm−1 corresponded to either alcohol or carboxylic O-H bonds. The band at 2935 cm−1 was assigned to C-H bonds, but its low wavenumber suggested they were The Fourier-transform infrared (FTIR) spectra obtained are shown in Figure 6. Numerous observable bands in the FTIR spectra belonging to many functional groups indicate that the degree of graphitization of these materials was not high. The most intense band at 3400 cm−<sup>1</sup> corresponded to either alcohol or carboxylic O-H bonds. The band at 2935 cm−<sup>1</sup> was assigned to C-H bonds, but its low wavenumber suggested they were mainly aliphatic bonds, since the aromatic ones give a signal at higher values. The 1710 cm−<sup>1</sup> band was due to C-O double bonds that can be carbonyl or carboxylic, whereas the neighboring band of 1610 cm−<sup>1</sup> was characteristic of C-C conjugated double bonds. Several bands between 1550 cm−<sup>1</sup> and 1300 cm−<sup>1</sup> can be attributed to the carbon backbone of these materials, for instance, those at 1515 cm−<sup>1</sup> and 1460 cm−<sup>1</sup> ; however, the presence of O-H and C-O bands in this region made assigning them somewhat problematic. Between 1300 cm−<sup>1</sup> and 100 cm−1, the bands were mainly due to oxygenated or nitrogenous functional groups (C-O or C-N bonds) [50].

groups (C-O or C-N bonds) [50].

**Figure 6.** Fourier-transform infrared spectra of the hydrothermal carbons. **Figure 6.** Fourier-transform infrared spectra of the hydrothermal carbons.

The scanning electron microscopy images of the samples in which hydrothermal carbonization had been carried out (Figures 7–9) showed a more defined structure than in the original material (Figure 10), in which the cells of the lignocellulosic material were more clearly observed. The treatment, therefore, removes matter from the original material, increasing the porosity. In some images, small spheres on the surface could be seen; it is known that obtaining hydrothermal coals from carbohydrates gives rise to spherical particles, so they could have come from dissolved substances of the chestnut cupule that underwent hydrothermal carbonization. Another option is that they were wax deposited on the surface of the biochar; given the similarity of composition in both cases, it is difficult to determine which of the processes gave rise to these structures. The scanning electron microscopy images of the samples in which hydrothermal carbonization had been carried out (Figures 7–9) showed a more defined structure than in the original material (Figure 10), in which the cells of the lignocellulosic material were more clearly observed. The treatment, therefore, removes matter from the original material, increasing the porosity. In some images, small spheres on the surface could be seen; it is known that obtaining hydrothermal coals from carbohydrates gives rise to spherical particles, so they could have come from dissolved substances of the chestnut cupule that underwent hydrothermal carbonization. Another option is that they were wax deposited on the surface of the biochar; given the similarity of composition in both cases, it is difficult to determine which of the processes gave rise to these structures. *C* **2023**, *9*, x FOR PEER REVIEW 9 of 18

mainly aliphatic bonds, since the aromatic ones give a signal at higher values. The 1710 cm−1 band was due to C-O double bonds that can be carbonyl or carboxylic, whereas the neighboring band of 1610 cm−1 was characteristic of C-C conjugated double bonds. Several bands between 1550 cm−1 and 1300 cm−1 can be attributed to the carbon backbone of these materials, for instance, those at 1515 cm−1 and 1460 cm−1; however, the presence of O-H and C-O bands in this region made assigning them somewhat problematic. Between 1300 cm−1 and 100 cm−1, the bands were mainly due to oxygenated or nitrogenous functional

**Figure 7.** SEM image of HC-A. **Figure 7.** SEM image of HC-A.

**Figure 8.** SEM image of HC-B.

**Figure 8. Figure 8.**  SEM image of HC-B. SEM image of HC-B.

**Figure 7.** SEM image of HC-A.

**Figure 9.** SEM image of HC-C. ered simply as carbonaceous materials with a low content of mineral matter. **Figure 9.** SEM image of HC-C.

**Figure 10.** SEM image of raw chestnut cupule.

The thermal analysis program can be found in the Supplementary Material (Table S1). Thermograms (Figure 11) show that the moisture of these materials was low, as well as the amount of residual ash, which was less than 1% in the three carbons. Regarding the volatile matter (from 50 to 80 min approximately), it followed the order HC-A > HC-B > HC-C, and the fixed carbon followed the opposite order. Therefore, they can be consid-

**Figure 10.** SEM image of raw chestnut cupule. **Figure 10.** SEM image of raw chestnut cupule.

**Figure 9.** SEM image of HC-C.

The thermal analysis program can be found in the Supplementary Material (Table S1). Thermograms (Figure 11) show that the moisture of these materials was low, as well as the amount of residual ash, which was less than 1% in the three carbons. Regarding the volatile matter (from 50 to 80 min approximately), it followed the order HC-A > HC-B > HC-C, and the fixed carbon followed the opposite order. Therefore, they can be consid-The thermal analysis program can be found in the Supplementary Material (Table S1). Thermograms (Figure 11) show that the moisture of these materials was low, as well as the amount of residual ash, which was less than 1% in the three carbons. Regarding the volatile matter (from 50 to 80 min approximately), it followed the order HC-A > HC-B > HC-C, and the fixed carbon followed the opposite order. Therefore, they can be considered simply as carbonaceous materials with a low content of mineral matter. *C* **2023**, *9*, x FOR PEER REVIEW 11 of 18

**Figure 11.** Thermograms of the hydrothermal carbons. **Figure 11.** Thermograms of the hydrothermal carbons.

The results of the organic elemental analysis (C, H, N, S, and O) and the ash content (Table 5) show that the samples prepared at a higher temperature (HC-B and HC-C) presented more carbon, nitrogen, and ash, and less hydrogen and oxygen. This indicates that the carbonization process was more intense when the temperature was raised. The size did not seem to be crucial, though, since the smaller the particle size, the more intense the The results of the organic elemental analysis (C, H, N, S, and O) and the ash content (Table 5) show that the samples prepared at a higher temperature (HC-B and HC-C) presented more carbon, nitrogen, and ash, and less hydrogen and oxygen. This indicates that the carbonization process was more intense when the temperature was raised. The size did not seem to be crucial, though, since the smaller the particle size, the more intense

effect of pressurized water should have been, and this was not observed in the results. The data in the table agree with the fixed carbon content, which tended to increase with

**Sample C, % H, % N, % S, % O, % Ash, %**  HC-A 57.50 6.56 0.53 0.05 35.31 0.05 HC-B 61.70 5.93 0.58 0.05 31.55 0.19 HC-C 63.70 5.89 1.06 0.05 28.46 0.84

On the other hand, EDX analyses were performed on all three hydrothermal carbons and the raw chestnut cupule (Table 6). This determination, which studied a more superficial area than the elemental analysis, confirmed the same trend. All carbons had been enriched in carbon and depleted in oxygen. It is worth noticing that the HC-C sample contained the most significant mineral matter. This effect is larger than the one found in the elemental analysis, which indicates that the hydrothermal treatment affected the surface

C 54.2 69.6 72.5 72.0 O 43.6 30.0 27.4 26.9 Mg 0.1 0.0 0.0 0.0 Al 0.1 0.1 0.0 0.0 Si 0.2 0.1 0.0 0.0 P 0.1 0.0 0.0 0.2 K 1.0 0.0 0.1 0.2 Ca 0.7 0.1 0.0 0.5 Fe 0.1 0.0 0.0 0.0 S 0.0 0.1 0.1 0.1

**Raw Material HC-A HC-B HC-C** 

**Table 5.** Elemental analyses of the hydrothermal carbons.

more than the interior of the material particles.

**Table 6.** EDX analyses of the hydrothermal carbons.

and hydrogen content.

the effect of pressurized water should have been, and this was not observed in the results. The data in the table agree with the fixed carbon content, which tended to increase with elemental carbon content, and the volatile matter, which tended to increase with oxygen and hydrogen content.


**Table 5.** Elemental analyses of the hydrothermal carbons.

On the other hand, EDX analyses were performed on all three hydrothermal carbons and the raw chestnut cupule (Table 6). This determination, which studied a more superficial area than the elemental analysis, confirmed the same trend. All carbons had been enriched in carbon and depleted in oxygen. It is worth noticing that the HC-C sample contained the most significant mineral matter. This effect is larger than the one found in the elemental analysis, which indicates that the hydrothermal treatment affected the surface more than the interior of the material particles.


**Table 6.** EDX analyses of the hydrothermal carbons.

### *3.2. Adsorption Experiments*

The kinetic study of adsorption processes makes it possible to relate the amount of adsorbate adsorbed per gram of adsorbent with the elapsed contact time between adsorbent and adsorbate. The variation in the excess PFONa concentration over time data is shown in Table 7. At 96 h, equilibrium can be considered established, having reached the maximum adsorption. Our preliminary studies show that HC-A and HC-B have adsorption of around 15%, while HC-C adsorbs 26%; therefore, further decontamination studies presented herein focused on the most promising adsorbent, HC-C.

**Table 7.** Adsorption kinetics results for HC-C.


<sup>1</sup> Average value of duplicate runs.

The Lagergren adsorption kinetic model assumes that the adsorption rate is proportional to the number of unoccupied sites available for adsorption. However, when we

applied it, we found it to be inappropriate for our data. On the other hand, the pseudosecond-order adsorption kinetic model assumes that the rate-limiting step is the chemical adsorption process, and it can be expressed mathematically, as shown in Equation (1):

$$\frac{t}{q\_t} = \frac{1}{k \cdot q\_\varepsilon^2} + \frac{t}{q\_\varepsilon} \tag{1}$$

where *t* is the time, *q<sup>t</sup>* is the amount of solute adsorbed at time *t*, *q<sup>e</sup>* is the amount of solute adsorbed at equilibrium, and k is the rate constant. The data are shown in Table 8, whereas the graph of *t/q<sup>t</sup>* vs. *t* is plotted in Figure 12.


**Table 8.** Data for the pseudo-second-order kinetic model.

**Figure 12.** Plot for pseudo-second-order kinetics. **Figure 12.** Plot for pseudo-second-order kinetics.

The data points follow a straight line, suggesting that the pseudo-second-order adsorption kinetic model applies to this system. The slope was 9.7137 × 10−4, and the intercept was 1.1612 × 10−3, so we can calculate the rate constant (*k* = 8.126·10−4 h−1·mg−1·g) and the amount of solute adsorbed at equilibrium (*qe* = 1029.47 mg·g−1). The data points follow a straight line, suggesting that the pseudo-second-order adsorption kinetic model applies to this system. The slope was 9.7137 <sup>×</sup> <sup>10</sup>−<sup>4</sup> , and the intercept was 1.1612 <sup>×</sup> <sup>10</sup>−<sup>3</sup> , so we can calculate the rate constant (*<sup>k</sup>* = 8.126·10−<sup>4</sup> <sup>h</sup> −1 ·mg−<sup>1</sup> ·g) and the amount of solute adsorbed at equilibrium (*q<sup>e</sup>* = 1029.47 mg·g −1 ).

Experiments were conducted with different HC-C to PFONa ratios to determine the carbon adsorption isotherm. The values obtained are shown in Table 9 and plotted in Figure 13. Experiments were conducted with different HC-C to PFONa ratios to determine the carbon adsorption isotherm. The values obtained are shown in Table 9 and plotted in Figure 13.


We found that the data obtained did not mathematically fit the Langmuir adsorption model, i.e., the adjustment of the model to the experimental data was inconsistent. On the other hand, the Freundlich model fit these results. Data for the Freundlich model are

**Table 9.** Adsorption isotherm results for HC-C 1. **Table 9.** Adsorption isotherm results for HC-C <sup>1</sup> .

1 Time = 96 h. <sup>1</sup> Time = 96 h.

shown in Table 10.

**Figure 13.** Adsorption isotherm for PFONa onto HC-C.

**Figure 13.** Adsorption isotherm for PFONa onto HC-C. **Figure 13.** Adsorption isotherm for PFONa onto HC-C.

**Figure 12.** Plot for pseudo-second-order kinetics.

**Table 9.** Adsorption isotherm results for HC-C 1.

ure 13.

1 Time = 96 h.

amount of solute adsorbed at equilibrium (*qe* = 1029.47 mg·g−1).

We found that the data obtained did not mathematically fit the Langmuir adsorption model, i.e., the adjustment of the model to the experimental data was inconsistent. On the other hand, the Freundlich model fit these results. Data for the Freundlich model are shown in Table 10. We found that the data obtained did not mathematically fit the Langmuir adsorption model, i.e., the adjustment of the model to the experimental data was inconsistent. On the other hand, the Freundlich model fit these results. Data for the Freundlich model are shown in Table 10.

The data points follow a straight line, suggesting that the pseudo-second-order adsorption kinetic model applies to this system. The slope was 9.7137 × 10−4, and the intercept was 1.1612 × 10−3, so we can calculate the rate constant (*k* = 8.126·10−4 h−1·mg−1·g) and the

Experiments were conducted with different HC-C to PFONa ratios to determine the carbon adsorption isotherm. The values obtained are shown in Table 9 and plotted in Fig-

> *c0* **(mg/L)** *ce* **(mg/L)** *qe* **(mg/g)**  40 36.6 3.40 70 61.1 8.95 100 81.5 18.55 130 103.2 26.80 165 126.6 38.40 200 148.9 51.15

**Table 10.** Data points for the Freundlich model.


Using Equation (2) and the adjusted values, the Freundlich model constants can be calculated:

$$
ln q\_{\varepsilon} = \frac{1}{n} \ln c\_{\varepsilon} + \ln \mathcal{K}\_{F} \tag{2}
$$

The calculated values of Freundlich isotherm constants were *K<sup>F</sup>* = 12.3099; and *1/n* = 1.9513 (r<sup>2</sup> = 0.9962).

The consistency with the Freundlich model indicates some degree of heterogeneity on the surface (not all active sites are identical) and the possibility that more than one layer can be adsorbed. Since the material used as adsorbent was amorphous, this relative heterogeneity can be expected. Moreover, the value of n indicates that the interaction between adsorbent and sorbate is mainly of the chemisorption type. Our results are promising, with a *q<sup>e</sup>* = 1029.47 mg·g −1 ; interestingly, they suggest that the maximum adsorption capacity had not been achieved at high concentrations. To put these results in perspective, we can refer to the excellent review by Liang et al.; in it, the authors compared the different adsorption methods for removing PFAS from water [37]. Analyzing the results obtained in water at pH 7, we found that, for granular activated carbon (GAC), a value of *q<sup>e</sup>* = 2.34–102.7 mg·g −1 is reported; and for powdered activated carbon (PAC), *q<sup>e</sup>* values of 16.9–500 mg·g <sup>−</sup><sup>1</sup> are found. On the other hand, for biochar derived from hardwood sawdust, the value obtained for *q<sup>e</sup>* at pH 7.2 is 2.4 mg·g −1 , and for pinewood at the same pH, *q<sup>e</sup>* is 2.1 mg·g −1 . Liang et al. found that ion exchange resins (*q<sup>e</sup>* = 525–1500 mg·g −1 ) and carbon nanotubes (*q<sup>e</sup>* = 140–1110 mg·g −1 ) show the highest adsorption capacities, yet they interestingly concluded that low-cost adsorbents with high adsorptive capacity are highly recommended in place of commercial AC due to their local availability, technical feasibility, and cost-effectiveness.

The mechanisms affecting PFAS adsorption have been reviewed by Liang et al. in an excellent review [37]. Diffusion, electrostatic interaction, hydrophobic interaction, ion exchange, and hydrogen bonds all can play a role in the adsorption process and are affected by the physicochemical characteristics of the adsorbent, such as pore size, pore distribution, and functional groups. In our case, HC-C is the carbon with the largest volume of mesopores; since PFONa is a large PFAS, it is essential for diffusion and adsorption that there are large numbers of mesopores. It is also the most graphitized hydrothermal carbon, i.e., the least polar; this also plays a vital role in the PFONa adsorption mechanism onto HC-C. Because of the hydrophobic nature of PFAS chains, the sorption of anionic PFAS can be dominated by hydrophobic interactions between the negatively-charged hydrophobic fluorinated tail and the surface of hydrophobic adsorbents, regardless of the electrostatic repulsion [51,52]. Particularly at high pH, hydrophobic interaction, which is not sensitive to changes in pH, should be the primary PFAS adsorption mechanism [53,54].

Given these results, and due to the advantages of the hydrothermal carbon herein reported, further work is underway to gain insight into the full potential of this adsorbent in water decontamination.

#### **4. Conclusions**

The preparation conditions of the hydrothermal carbons greatly influence the result; this allows for obtention of coals with tuned properties by optimizing the obtention process.

Carbonaceous material obtained by hydrothermal carbonization of sweet chestnut (*Castanea sativa*) cupule residues is a cost-effective, environmentally benign adsorbent, which effectively adds value to this biomass waste.

The carbons thus obtained show a limited surface area and a low microporosity while having a large volume of mesopores and macropores. All of them are acidic, with a low graphitization index, minimal moisture, and slight residual ash content; they are carbonaceous materials with a low mineral matter content.

It is possible to purify water contaminated with perfluorooctanoate by adsorption on this type of hydrothermal carbon. The adsorption process follows pseudo-second-order kinetics. Additionally, it fits the Freundlich model, showing that more than one layer can be adsorbed, which is expected in such an amorphous adsorbent.

**Supplementary Materials:** The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/c9020057/s1, Table S1. Thermal Analysis Program; Figure S1. SEM Images of HC-A; Figure S2. SEM Images of HC-B; Figure S3. SEM Images of HC-C; Figure S4. SEM Images of raw chestnut cupule.

**Author Contributions:** Conceptualization, I.M.L.-C. and S.I.; methodology, I.M.L.-C. and C.J.D.-V.; investigation, S.I., N.P., C.J.D.-V. and I.M.L.-C.; writing—original draft preparation, S.I., C.J.D.-V. and I.M.L.-C.; writing—review and editing, S.I., N.P., C.J.D.-V. and I.M.L.-C.; funding acquisition, C.J.D.-V. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Spanish Regional Government 'Junta de Extremadura' and European Regional Development Fund (ERDF), grants IB20026 and GR21107.

**Data Availability Statement:** Data may be obtained from the authors under reasonable request.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
