Performance assessment of activated carbon thermally modified with iron in the desulfurization of biogas in a static batch system supported by headspace gas chromatography

A static batch arrangement composed of anti-leak vials coupled to gas chromatography is proposed as a complementary system for performance assessment of biogas desulfurization by adsorption. For testing, a modified commercial activated carbon produced by controlled thermal treatment in the presence of iron(III) species improved biogas desulfurization. The adsorbents showed a superior hydrogen sulfide removal compared to ordinary one. Pseudo-first-order, pseudo-second-order, and Bangham’s kinetic models were used to fit experimental data. All studied samples followed pseudo-first-order model, indicating the predominance of physisorption, and Bangham’s model, confirming that the micropores structure played an important role for gases diffusion and adsorbent capacity. Additionally, the materials were characterized by N 2 adsorption–desorption, X-ray diffraction, infrared spectroscopy, scanning electron microscopy and energy-dispersive spectroscopy. The thermal treatment associated with iron impregnation caused significant modifications in the surface of the materials, and the iron species showed two main benefits: an expressive increase in the specific area and the formation of specific adsorption sites for hydrogen sulfide removal. The results reinforce the advantages of iron-modified adsorbents in relation to their non-modified counterparts. The analytical methodology based on the confinement of multiple gases contributes to improving the understanding of the hydrogen sulfide adsorption process using pressure swing adsorption technology.


Introduction
Hydrogen sulfide (H 2 S) is one of the most harmful contaminants of biogas, and its removal is mandatory to produce safe and sustainable biofuels based on renewable hydrocarbons.The main reasons for this are its high potential to cause severe corrosion in pipelines and metallic equipment, and the significant amount of sulfur oxides (SO x ) that can be released into the atmosphere during the combustion of H 2 S-containing gaseous mixtures (Becker et al. 2022).
The biogas composition is quite variable and depends on the type of biomass used in the anaerobic digestion process (Calbry-Muzyka et al. 2022); however, in any case, H 2 S is the primary source of sulfur in biogas, and its level usually ranges from 700 to 2,000 ppmv but can reach up to 20,000 ppmv (Dumont 2015;Egger et al. 2023).
The European Committee of Standards specifies that biomethane, which is a gas with methane content usually higher than 90%, should have an H 2 S concentration < 20 ppmv or ~ 0.0015% (vol/vol) (Okoro and Sun 2019), as well as the Brazilian Agency for Oil, Natural Gas and Biofuels (ANP) determines that the maximum limit for H 2 S in biomethane is 10 ppmv (ANP, Resolution n° 685/2017).In fact, some technologies compatible with (bio)methane require very low concentration of impurities, such as solid oxide fuel cells (SOFCs), in which the previous removal of volatile organic compounds (VOCs), from sulfur species to aromatic, chlorine and siloxane compounds, is mandatory in order to maximize the reliability of the equipment and increase its long-term service life (Illathukandy et al. 2022;Papurello et al. 2018;Zhu et al. 2020).
Therefore, raw biogas is typically submitted to a refining process to remove impurities and meet quality standards.This process is usually divided into two steps: a cleaning step, in which H 2 S, H 2 O, VOCs, siloxanes and other minor impurities are removed from the biogas (Calbry-Muzyka et al. 2022;Lecharlier et al. 2022;Tiwari et al. 2022), and an upgrading step, in which carbon dioxide, the primary contaminant, is removed to increase the methane content and raise the heat value of the biofuel (Pinto et al. 2022).Although siloxanes and other minor compounds are also potentially corrosive to metal components and can cause lubricant leaks (Boldrin et al. 2016;Illathukandy et al. 2022;Papurello et al. 2019), the concern with H 2 S is relatively greater mainly due to its high concentration.
Cleaning and upgrading of biogas are usually supported by four technologies: pressure swing adsorption (PSA), water/chemical scrubbing, membrane separation, and cryogenic distillation (Barbera et al. 2019).Additionally, membrane biofilm or microorganism-based filters can be used instead of physicochemical approaches, or even in association with them, to maximize the efficiency of the desulfurization step, despite its limitations due to the fragility of the microorganisms (Dumont 2015).
The choice of technology depends on several factors, such as the type of biomass available for anaerobic digestion, specific biomethane end-use segment, and even the producer's investment capacity (Ardolino et al. 2021).Nonetheless, many biogas producers, on both small and large scales, use pressure swing adsorption technology or fixed bed columns stuffed with an adsorbent material and operating at room pressure for biogas cleaning and upgrading.Adsorption processes rule these two cases, and the surface of the filling material must have high capacity and selectivity for a given target contaminant.
Depending on the interaction strength between the target and the sorbent surface, the adsorption processes can occur by physisorption or chemisorption.While in physisorption, the sorbate and sorbent interact by weak forces and allow a more straightforward regeneration, in chemisorption, this interaction involves more energy and is generally irreversible (Bezerra et al. 2011).
Activated carbon (AC) is employed as an adsorbent for molecules and ions, both in gaseous and liquid systems.AC is obtained from chemical or physical activation of the biochar, which is the solid portion resultant from the pyrolysis of lignocellulosic biomass (Indayaningsih et al. 2017;Qu et al. 2018).Usually, AC shows a high surface area, good commercial availability, and relatively low cost, which are significant advantages for biogas refining worldwide (Sinha et al. 2020).
Its surface properties, such as surface charge distribution, functional group availability, pore size, and specific surface area, are directly correlated with H 2 S removal from biogas (Zulkefli et al. 2023).An accepted mechanism for H 2 S removal by AC is based on the direct reaction between oxygen-containing groups localized on the AC surface and H 2 S molecules to form oxidized elemental sulfur due to the presence of pyrone groups (Bian et al. 2019).Therefore, researchers have sought to increase the reactive sites in AC to overcome the H 2 S fluctuations generated in the anaerobic digestion that demands a longer-term fixed-bed column in order to minimize stoppages for adsorbent replacement (Ou et al. 2020).
An adsorbent with a large capacity would lead to fewer operational stoppages for changing the column packing, and the producer's profit could be positively impacted if biogas processing had fewer interruptions.At least for the small biogas producers worldwide, an adsorbent such as AC, with enhanced capacity for H 2 S removal, would be essential.
A systematic strategy to increase the adsorption capacity of materials has been to focus on modifying the AC surface by impregnation with chemical species of high affinity for H 2 S, such as iron, zinc, copper, and other metal-based compounds (Barroso-Bogeat et al. 2019;Choudhury and Lansing 2021).
Like other d-block metals, iron presents an outstanding thermodynamic energy to adsorb H 2 S. Therefore, an intuitive strategy has been to promote the contact between metal-containing materials and biogas, aiming at its purification (Cheremisina et al. 2022).For instance, Choudhury and Lansing (2021) synthesized two carbon-based biosorbents under 500 °C and promoted their impregnation with iron(III) chloride from an aqueous solution.The H 2 S adsorption capacity in the fixed-bed column increased by 2.5-3.9 compared to the unmodified one.Cimino et al. (2020) developed a CuO-ZnO mixture dispersed onto commercial activated carbon, oxidizing H 2 S to form sulfates or elemental sulfur.According to the authors, the presence of O 2 or H 2 O favors oxidative reaction and the formation of elemental sulfur on copper sites but not in zinc ones.On the other hand, sulfate species are predominant in the absence of oxygen or moisture, and they are formed utilizing lattice oxygen from metal oxide clusters or superficial oxygen species.
Another aspect that has been extensively addressed is related to functional groups on the AC surface and their role in H 2 S adsorption.Shen et al. (2019) performed a density functional theory (DFT) study and found that groups such as carbonyl, ester, and, preferably, pyrone present an important effect on H 2 S adsorption.Pyrone groups are considered basic and are produced by exposing heat-treated activated carbon to air.Authors inferred that the combination of carbonyl and epoxy oxygen atoms from the pyrone group creates an appropriate chemical environment for dissociating H 2 S and its adsorption on pyrone groups predominantly ruled by chemisorption with a highly exothermic adsorption enthalpy.
These examples depict that modifying carbon-based materials by d-block metals and thermal treatment are important strategies to enhance adsorption-based technologies, i.e., PSA and fix-bed columns, to meet biogas production standards.In any adsorption case, the assessment of the mass capacity of the adsorbent material is crucial for process planning and viability.
In this work, a static batch system composed of antileak vials were used for gas confinement in order to evaluate the synthetic biogas desulfurization by headspace gas chromatography using a thermal conductivity detector.The studied adsorbent was a commercial activated carbon impregnated with iron species and submitted to thermal treatment to boost its performance.

Preparation of the adsorbents
Adsorbents were prepared by impregnation with FeCl 3 •6H 2 O and thermal treatment.The activated carbon (AC) was impregnated with iron according to the method described by Choudhury and Lansing (2021) with slight modifications.In brief, approximately 2.0 g of AC was added to an Erlenmeyer flask containing 60 mL of a 1000 mg L −1 of Fe(III) dissolved in 1.0% (vol/vol) HNO 3 aqueous solution.The magnetic stirring (~ 500 rpm) was controlled from 5 to 300 min.Afterward, the impregnated AC was filtered using a 0.45 µm nylon syringe filter without washing to maximize iron ions availability during the thermal treatment step.The supernatant was analyzed by UV-Vis spectroscopy to monitor the iron(III) ions concentration through the formation of the complex red pentaaqua(thiocyanate)iron(III) [Fe(NCS)(OH 2 ) 5 ] 2+ , detected at λ max = 460 nm.
For the calculation of the iron(III) concentration remaining in the supernatant, a calibration curve with the dynamic range of 1.0-10.0mgFe 3+ L −1 and Eq. 1 was considered: in which q e is the maximum adsorption capacity (mgFe g −1 ); C 0 is the initial concentration of Fe(III) (mg L −1 ); C 1 is the concentration of Fe(III) at equilibrium (mg L −1 ); V is the volume of Fe(III) aqueous solution (L); and m is the mass of AC (g).
Then, a portion of each material was dried at 100 °C for 4 h and treated in a muffle furnace at different temperatures (300, 500, and 700 °C) for 60 min with a thermal increment of 10 °C min −1 .Upon controlled cooling in a desiccator, the samples were stored and labeled AC3Fe, AC5Fe, and AC7Fe, referring to the thermal treatment at 300, 500, and 700 °C, respectively.A sample labeled AC (commercial activated carbon submitted only to drying at 100 °C for 60 min for moisture removal) was considered a reference for the desulfurization tests.For language standardization purposes, in this work, a treated material means chemically or thermally treated; therefore, it does not comprise a simple drying to remove moisture.An additional material labeled AC3 was obtained by thermal treatment of the AC at 300 °C for 60 min without previous impregnation with iron.

Equipment for analysis and characterization
Powder X-ray diffraction analysis (PXRD) was conducted in a Bruker D2 Phaser diffractometer operating at 30 kV and 40 mA with copper radiation (Cu-Kα = 0.15418 nm) in the range 10-75° with dwell a time of 2°min −1 .
Scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS) were recorded in an FE-SEM TESCAN model VEGA 4, corresponding to a tungsten thermionic emission system.The measurements were conducted at 15-30 kV acceleration potential under a high vacuum.The surfaces of the samples were coated with gold (sputtering-Denton Vacuum Desk V, using a current of 30 mA for 2 min).The absorption peaks due to the gold layer were attenuated in the data processing software Tescan Essence.
Attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectra were registered on a Perki-nElmer spectrometer model Spectrum Two in the range ( 1) 400-4000 cm −1 .The spectra were collected with a resolution of 4 cm −1 and an accumulation of 16 scans.Adsorption/desorption of N 2 analysis was performed to measure the specific surface area, the average pore diameter, and the pore volume of the ACs.For this purpose, the BET (Brunauer et al. 1938) and the BJH (Barret et al. 1951) methods were applied to the corresponding nitrogen desorption isotherms.The method de Boer (De Boer et al. 1966) was used for micropore and external area calculations.Data were collected on an ASAP 2020 N Automatic Physisorption Analyzer from Micrometrics at 77 K.The samples were previously degassed by treatment at 150 °C until the system pressure had achieved 10 μmHg.The nitrogen adsorption data were obtained using approximately 0.1 g of the sample.
UV-Vis analysis was carried out on a spectrophotometer Biochrom model S60 Libra at λ max = 460 nm with a resolution of 2 nm.
Gas chromatography (GC) was used to determine H 2 S using experimental parameters according to the column supplier (Agilent 2023).An equipment Agilent model 7820A, equipped with a thermal conductivity detector (TCD) and an HP-Plot Q column, 0.53 mm × 30 m, was employed in the experiments.The main conditions of analysis were: carrier gas He, flow rate 8.6 mL min −1 ; oven temperature 60 °C; initial temperature (T i ) 60 °C (by 2 min), temperature rate 30 °C min −1 , final temperature (T f ) 240 °C (by 1 min); the temperature of the detector 250 °C; splitless injection mode; sample volume 250 μL.

Desulfurization experiments
H 2 S removal from the gaseous mixture (synthetic biogas) was evaluated in static batch experiments, and the analytical steps are depicted in Fig. 1.Although the continuous flow approach in a fixed bed column can be a better way to represent the most used operational processes, the capacity of H 2 S removal of the samples could be determined in batch experiments, which represents, in a simplified way, the mass transfer equilibrium expected in a PSA system.In addition, this approach can allow a significant cost reduction owing to the lower consumption of synthetic gaseous mixture.
For the experiments, all samples were grounded in a ceramic mortar, sieved in the range, 0.044-0.063mm followed by drying at 100 °C for 60 min and cooling down to room temperature in a vacuum desiccator.Each sample was diluted before weighing by a factor of a hundred (10.0 mg of sample to 1.0 g of diluent) by addition of powdered dry sodium chloride to enable weighing in a reliable range of the analytical balance with resolution of 0.1 mg (Ohaus PR224BR/E).Both sample and diluent were placed together in a flask with a lid and vigorously shaken until a homogeneous greyish mixture was formed.Afterward, 50.0 mg of the diluted material was weighed in a 12-mL glass tube purchased from Labco Exetainer sealed with plastic screw-caps containing chlorobutyl septum (anti-leak vials).The air pressure inside the tube was lowered to ~ 3.0 kPa for 1 min with a vacuum pump (Tube A).In parallel, another sealed tube (Tube B) had its internal content purged by a gentle stream of synthetic biogas.After a few seconds, 12.0 mL of synthetic biogas were carefully sampled during the purge to ensure room pressure stability and immediately transferred to Tube A. The contact time between diluted materials and synthetic biogas started to be registered, and the adsorption process was monitored with periodic headspace sampling in intervals of 15 min using a 250-μL gastight syringe by gas chromatography until the H 2 S concentration remained constant for at least three consecutive chromatographic runs.A blank test was performed following the same protocols described above, however, using an excess of 50.0 mg of the diluent (NaCl) in a control tube to prove the diluent's inertia.For this reason, some aliquots were taken from the control tube for GC analysis after 250 min of contact with the synthetic biogas.
The calculations of maximum adsorption capacity ( q H 2 S ) were carried out considering the ideal gas behavior of the gaseous mixture with initial pressure ~ 101.3 kPa and constant temperature at 25 °C (298.15K), which was possible in an air-conditioned room.According to Clapeyron's equation for ideal gases, the pressure inside the tube must decrease with decreasing number of moles of the gaseous mixture, which is expected to have happened after each new aliquot withdrawn from the tube.From the difference between the initial and remaining ratio between area of chromatographic peaks of H 2 S and total gases, it was possible to calculate the mass of H 2 S adsorbed on the evaluated materials according to Eq. ( 2). in which q H 2 S is the adsorption capacity of H 2 S at equilib- rium (mg g −1 ); n H 2 S initial = moles of H 2 S contained in 12.0 mL of standard gaseous mixture at 101.3 kPa and 298.15 K; = ratio between area of H 2 S and total area of a given reaction; = ratio between area of H 2 S and total area of the standard mixture; MM = molar mass of H 2 S (34,100 mg mol −1 ); m ads = mass of adsorbent (g).
Analysis of variance test (ANOVA) with 95% confidence level was carried out on the results of performance.The null hypothesis was that the average of the H 2 S removal from synthetic biogas in the presence of the modified adsorbents would be the same.
in which q t and q e are the amount of H 2 S adsorbed (mg g −1 ) at time t (min) and after equilibrium of adsorption, respectively; k 1 is the pseudo-first-order rate constant (min −1 ); t is the elapsed time for the adsorption reaction (min). in which q t and q e are the amount of H 2 S adsorbed (mg g −1 ) at time t (min) and after equilibrium of adsorption, respectively; k 2 is the pseudo-second-order rate constant (g mg −1 min −1 ); t is the elapsed time for the adsorption reaction (min). (2) (5) q t = q e 1 − e −k b t n in which q t and q e are the amount of H 2 S adsorbed (mg g −1 ) at time t (min) and after equilibrium of adsorption, respectively; k b (min −n ) and n are Bangham's constants; t is the elapsed time for the adsorption reaction (min).

Preparation of the adsorbents
The iron concentration of the ACs was indirectly determined by UV-Vis spectroscopy analysis of the iron(III) ions remaining in the supernatant as [Fe(NCS)(OH 2 ) 5 ] 2+ complex.This method was chosen instead of the direct determination of iron on the solid to avoid an underestimated result due to a possible difficulty to total solubilization of the iron by acidic digestion.Considering that the impregnated materials were not submitted to washing during their preparation, the difference between the concentration of iron(III) ions in the supernatant and the initial iron(III) concentration corresponds to the amount of iron adsorbed on the ACs.The AC showed a high affinity for iron(III) ions, and only 5 min was required to reach a high uptake from the acid solution.All the reaction times led to a similar iron content (average adsorption of 28.8 ± 0.20 mg of Fe by g of AC), representing 97.6 ± 0.25% of the total iron(III) ions available in the acid solution.Although there was no significant difference in iron load in ACs concerning reaction time, extended contact usually helps to create better control of experimental parameters in the impregnation system.In addition, the iron load was comparable to those obtained by Chang et al. (2010) and Kalaruban et al. (2019), who found values ranging from 19 to 26 mgFe g −1 using granular activated carbon.Nevertheless, Martins et al. (2022) obtained a higher value, 178 mgFe g −1 , using commercial activated carbon under strong acid conditions.These disparities indicate that iron impregnation has a high dependence on the experimental procedures, as well as on characteristics of the material used, such as specific area, pore size/volume, the intensity of the thermal treatment before the impregnation step, weight ratio metal/ AC, among others.
An important consideration must be made regarding the iron load in each material after thermal treatment.Because mass loss was detected as the temperature was raised due to the thermal decomposition of activated carbon under an O 2 -rich atmosphere, an increase in the iron load in materials can be expected (Table 1).The mass loss of 15.1% depicted by AC3Fe was responsible for an estimated iron load increasing from 28.8 to 33.9 mg g −1 .In addition, the sample AC3 (without iron) presented 9.89% mass loss after thermal treatment, confirming the role of the iron species in the modification of the carbonaceous materials.On the other hand, owing to the treatment intensity performed over AC5Fe and AC7Fe, these materials depicted a more substantial mass loss, reaching 72.5 and 79.0%, respectively, indicating intensive mineralization.Considering that the residence time in the muffle was relatively short (60 min), the thermal treatment may not have been able to decompose all the organic matter, leading to an increasing loading of iron in the presence of the remaining carbonaceous matrix, which reached 104.7 mgFe g −1 for AC5Fe and 137.1 mgFe g −1 for AC7Fe.

Characterization of the modified ACs
Often, two crucial characteristics to consider for activated carbon applied to adsorption processes are a high specific area (high pore volume) and a well-formed pore structure.Table 2 shows that AC exhibited a specific area of 325.3 m 2 g −1 , which is a little lower than the typical range reported by the literature for activated carbons 500-2000 m 2 g −1 (Fahmi et al. 2019), as well as compared to the exceptional areas that this type of material can reach (~ 3500 m 2 g −1 ) if the synthesis conditions are optimized (Ma et al. 2020).Additionally, AC showed 172.8 m 2 g −1 (53% of the total area), referring to micropores (S micro ) and 153.0 m 2 g −1 (47.1% of the total area), referring to the external area (S ext ), which is interpreted as the available area outside the pores.The average pore size was 1.53 nm, confirming the predominance of micropores, and the average pore volume was 0.25 cm 3 g −1 .Table 2 shows that the increase in the temperature of the thermal treatment up to 300 °C under air led to a remarkable increase in the specific area, with AC3 reaching 479.9 m 2 g −1 , accompanied by a slight decrease of its average pore size (1.36 nm), and an increase of its pore volume (0.32 cm 3 g −1 ).Considering the different pore size distributions and morphologies possible in activated carbon, including conical-shaped pores, these results suggest that the thermal treatment caused an abrasion on the surface of the activated carbon particles, exposing regions with smaller opening micropores.Moreover, new micropores may also have appeared.The gain in the specific area was even more evident for AC3Fe material (533.4 m 2 g −1 , 64 wt% higher than AC), which was treated in the presence of iron species.Its average pore size (1.37 nm) was too close to those presented by AC3; however, its pore volume was slightly higher than for AC3, indicating that despite the 15.1% mass loss (Table 1), the iron species contributed to the formation of new pores and deepening of the existing ones.The thermal treatment at 500 °C for AC5Fe did not cause a significant decrease in its S BET (317.2 m 2 g −1 ) and its pore volume (0.22 cm 3 g −1 ) in comparison with AC; it seems that the high mass loss (72.5 wt%) had been offset by the gain in porosity.Conversely, a significant drop at higher temperatures was observed for AC7Fe, with the total specific area (S BET ) decreasing to 236.6 m 2 g −1 and pore volume to 0.17 cm 3 g −1 .Moreover, the decrease in S ext /S BET ratio presented by modified ACs (from 47.1 to ~ 30%), and consequently, the increase in S micro /S BET ratio (from 53.1 to ~ 70%), occurred due to the appearance of micrometric pores, depicted in Fig. 3 (see SEM analysis), of an exceptionally large size of about ~ 5-15 μm.The walls of these "big pores" are expected to be full of interconnected micropores (< 2 nm), originally blind and closed pores, which were connected due to the intense abrasion caused by the combination of high temperature and iron species.
In general, the average pore volume was directly proportional to the specific area of the materials.All nitrogen adsorption-desorption isotherms registered at 77 K showed a similar profile (Fig. 2).All isotherms were classified as Type I, referring to a micropore-containing material (Pires et al. 2021); however, the hysteresis loop type H4 seen ranging from P/ P 0 0.45 to 0.92, that usually arises from capillary condensation, indicates the existence of additional small mesopores forming a heterogeneous porous matrix (Bläker et al. 2019).The insets in Fig. 2 show the average pore size in alignment with the results presented in Table 2, which means a predominance of micropores and small mesopores.
SEM analysis aids in comprehending the surface characteristics of a given material proposed as an adsorbent and its adsorption dynamic.Therefore, information about roughness and agglomeration, mainly the existence of pores, their shape, size, and distribution over the surface of the particles, are fundamental to increase the adsorbent performance.
Figure 3a-j shows a set of images that depict the profile of ACs samples based on typical pyrolyzed biomass composed of graphite and probably a few contributions of unchanged cellulose fragments, in which lignin might remain stable even at about 900 °C (Waters et al. 2017).The first row of images (Fig. 3a-e) shows a comparison for at the same magnification (3kx), while the second row (Fig. 3f-j) shows some characteristic information about each material at different magnifications.The presence of cellulosic fragments is common to all materials, however in those thermal treated materials the mineralization process caused changes in the surface texture.Mineralization can be evidenced by the appearance of fracture and accumulation of small aggregates (Fig. 3b), as well as formation of macropores owing to the breakup of thinner pore walls already existing in the starting material (Fig. 3i).The images in Fig. 3c and h give a notion of the cross section of the cellulosic fragments in different directions.In both cases, the views allow to notice the massive presence of pores with ~ 10 μm in size; they are well organized and oriented parallel to each other along the analyzed fragment.The thermal treated material AC3 (without iron chloride) presented an accumulation of particles apparently inorganic in nature, with a degree of organization that resembles small quartz crystallites.Interestingly, the presence of these crystallites as well as the decrease in the availability of cellulosic fragments was more intense for the AC3 material than for the AC3Fe material, suggesting that the iron chloride modifier has a relevant role during the thermal degradation of activated carbon.
Figure 4a-j shows the sample's EDS spectrum, element mapping, and weight fraction.As expected, the micrometric particles are predominantly formed by carbon.However, oxygen and silicon were found in all materials, suggesting the presence of a significant content of silicon oxide.For AC3Fe and AC5Fe, it was possible to detect iron and chlorine, which is in accordance with the starting material FeCl 3 .The EDS mapping also shows that iron species are highly dispersed over the surface of the materials, and their content is lower than those obtained by UV-Vis (Table 1), which is typical for surface probing techniques such as EDS analysis.Another expected information is about the increase of the silicon content with the increase of the temperature used in the thermal treatment.Although EDS analysis results are highly dependent on the specific point of the analyzed particle and may not represent the bulk material, the image in Fig. 4g represents the accumulation of silicon oxide species.The particles chosen for scanning are the same as shown in Fig. 3g (SEM), showing the significant amount of silicon oxide species embedded in the carbonaceous matrix.Unlike AC3Fe, AC3 material did not present iron or chlorine in the EDS mapping; indeed, it showed a large predominance of carbon and a significant amount of silicon and oxygen as a result of the thermal decomposition of the carbonaceous matrix (Fig. 4i and j).The oxygen content increased from 6.04 wt% for AC (only dried at 100 °C) to 31.8 wt% for AC7Fe (700 °C), confirming that the ACs surface was strongly modified during the thermal treatment.Additionally, trace levels of certain elements such as aluminum (Al), potassium (K), and magnesium (Mg) were detected, which can be attributed to the typical composition of biomass.A persistent broad peak of relatively low intensity at ~ 2.1 kV referring to gold spectral line from samples metalization can be observe in all spectra after attenuation of unwanted signals.This peak was not accounted for centesimal composition.
In order to further characterize the materials, especially the presence of inorganic-thermally stable compounds and to enhance the comprehension on the role of the iron(III) chloride during the thermal treatment, the PXRD analysis was carried out, and the results are shown in Fig. 5.All samples presented diffraction peaks that stood out compared to the carbonaceous phase.The predominant signals were related to silica-based crystals, which agrees with the EDS results, indicating the presence of silicon and oxygen.Despite the high load of iron species (Table 1), it was impossible to identify diffraction peaks related to iron oxides or other iron crystalline Fig. 2 Nitrogen adsorption-desorption loops obtained for the evaluated materials at 77 K.The inset shows the pore size distribution remaining between micropores and small mesopores phases, confirming that iron compounds may be highly dispersed onto the AC surface, avoiding crystallization.
The samples AC and AC3 presented some narrow diffraction peaks overlaying a discrete amorphous halo ranging from 20 to 30° (2θ).As a reference for the other materials, the diffraction peaks for AC were ticked by Q and assigned to quartz as follows: 2θ = 21.45°(Q1), 27.09° (Q2), and 36.81°(Q3), referring to the planes (100), (101), and ( 200), respectively, according to ICSD Collection code 079634.Quartz can be considered a common inorganic impurity to biomass (Bhattacharjya and Yu 2014).Silicon from the soil converts into amorphous silica within the plant (Luyckx et al. 2017), forming crystalline phases under heating.Moreover, Yang et al. (2022) showed that silica-phase transformation could occur at lower temperatures than those widely reported, for instance, if metals or metal compounds are present.So, the peak at 2θ = 22.26° was assigned to tridymite, a polymorph of silica stably formed above 870 °C according to ICSD Collection code 408,281 (Li et al. 2021), which is a temperature usually reached in pyrolytic processes for activated carbon production.
Also, as expected for a carbon-based material, the presence of graphite is expected to be detected.As graphite usually shows an intense diffraction peak at 2θ = 26-27° according to ICSD Collection code 031829 (Ganash et al. 2019), the peak at 2θ = 27.17° in Fig. 5 was considered as overlapping between graphite and quartz signals, in which the former is partially hidden by the last.Despite that, a discrete splitting of this diffraction peak can be observed in sample AC5Fe due to long-range atomic order changes, shifting the graphite peak to 2θ = 26.6°.In addition, Fig. 5 shows that the thermal treatment allowed slight changes in crystallinity and the number of detectable peaks of quartz due to the intense decomposition of graphite and the gradual disappearance of the amorphous halo.
It is notable that, due to the thermal treatment at 300 °C, the AC3 material presented relatively intense quartz-based peaks clearly higher than the AC, but nevertheless, also higher than the AC3Fe, whose thermal treatment was performed in the presence of iron(III) chloride.In addition, by inspection of Fig. 5, the profile of the AC3 is similar to that obtained for AC7Fe, which is known to have undergone a high degree of thermal decomposition and shows intense quartz peaks.As aforementioned, the iron(III) chloride played an important role in the production of modified activated carbons, resulting in more interesting textural properties, such as higher specific area, higher pore volume, increase of the microporosity, and deepening of pores.At the same time, iron(III) chloride seems to have mitigated the abrupt thermal decomposition of the ACs and promoted a controlled increase in porosity.
ATR-FTIR spectra of the materials are shown in Fig. 6.Although thermal treatment applied on the samples was enough to cause a significant mass loss, in general, the spectra show slight changes when the non-modified AC is compared to the other materials.Despite that, a remarkable change occurred in the region below 1,000 cm −1 , in which the ACs treated at higher temperatures (500 and 700 °C) presented two intense absorption bands centered at 790 and 460 cm −1 .These bands can be related to stretches involving heavier atoms, such as Fig. 3 SEM images depicting the roughness of the surface of the studied materials AC (a, f), AC3Fe (b, g), AC5Fe (c, h), AC7Fe (d, i), and AC3 (e, j).The upper row shows comparative images for the materials at the same magnification scale (3kx).The down row shows some relevant details for each material under relevant magnification Fig. 4 EDS results showing spectrum, element mapping and the weight fraction for each sample.Inset depicts the increase of the oxygen, silicon and iron contents as the thermal treatment temperature was raised those caused by Fe-O due to the decrease of the carbonaceous content during the thermal treatment.The other absorption bands are in accordance with the expected for activated carbon produced from vegetal biomass, including the broad absorption band centered at 3450 cm −1 , which was attributed to O-H stretching modes from H 2 O and hydroxyl groups.The counterpart bending mode absorption band from H 2 O was detected at around 1630 cm −1 .Alkyl structures of aliphatic groups, probably from a remaining portion of lignin, were evidenced by the set of absorption bands ranging from 2969 to 2924 cm −1 and were assigned to C(sp 3 )-H stretching modes (Waters et al. 2017).Moreover, the inorganic portion formed by silicon compounds was identified by the broad absorption band centered at 1090 cm −1 , assigned to Si-O vibration modes from siloxane (Si-O-Si), probably overlapped on a less intense band referring to silanol (Si-OH) groups.The significant content of silica depicted in ATR-FTIR spectra is in accordance with results obtained by PXRD and EDS.

Evaluation of the materials in desulfurization
The performance of the modified and non-modified ACs was comparatively assessed by adsorption batch experiments using a synthetic gaseous mixture composed of CH 4 , CO 2 , and H 2 S. The major gases (methane and carbon dioxide) had their concentrations practically constant throughout the experiments.Figure 7a shows five typical chromatograms obtained for the standard gaseous mixture sampled by the headspace technique and injected to monitor the GC system response, presenting an apparent profile of overlapping peaks and analytical repeatability.The average deviation calculated considering the area ratio of the peaks for these five injections of the standard mixture was 11.8%, which was considered a reasonable deviation for the experimental arrangement used.In the open-window image, it is possible to notice only the two peaks corresponding to major components, CH 4 and CO 2 , centered at 1.12 and 1.18 min, respectively, and a third peak attributed to H 2 O centered at 3.21 min (Agilent 2023), which is expected due to small contamination of the system with atmospheric moisture during the sampling.The peak corresponding to H 2 S centered at 2.44 min can only be observed in the inset of Fig. 7a owing to its relatively low concentration in the mixture.Although all peaks were detected before 5 min, each run was performed up to 9 min to ensure total elution of the components.These four peaks were considered for area integration.Figure 7b shows the decrease in the area ratio between H 2 S and the total area along the batch adsorption experiments, with the curves presenting similar profiles but reaching a plateau at different levels.All materials were able to adsorb H 2 S, resulting in a decrease in peak area ratio, including AC, AC3, and AC7Fe, with similar results, corresponding to 53.1, 55.6 and 49.3%, respectively.However, those ACs modified with iron and submitted to thermal treatment at 300 and 500 °C presented better performances in removing H 2 S, with AC5Fe and AC3Fe reaching 83.9 and 97.8% of the decrease in their area ratios.These results reinforce the important role of the thermal treatment carried out in the presence of iron species.Additionally, the performance of the AC7Fe was affected by the excessive temperature applied to it, and the decrease in its specific area (Table 2) was responsible for a more discreet result.Nevertheless, the high iron content on the AC7Fe surface Fig. 5 PXRD patterns of the evaluated materials evidencing the increase in quartz peak intensity with increase of the temperature Fig. 6 ATR-FTIR spectra of the produced materials evidencing the appearance of bands due to the presence of Fe-O stretches contributed to the feasibility of the H 2 S adsorption.However, under more drastic thermal treatment conditions, this material has a small contribution from the structure of the original activated carbon.
The adsorption capacities of the materials were calculated based on H 2 S removal as a function of time (Fig. 8).The materials AC, AC3, and AC7Fe presented similar H 2 S adsorption capacities, reaching 38.1, 40.0, and 37.9 mgH 2 S g −1 of adsorbent, respectively.
A preliminary analysis considering only the average deviation indicates that these results are equivalent and confirm that the modification caused by the thermal treatment in the absence of iron, that is, gain in micro/ mesopore area and micro/mesopore volume, did not favor a significant increase in the adsorption capacity for AC3.
Indeed, the results obtained from textural analysis and EDS analysis provided relevant information on the performance of the studied materials.For instance, the AC7Fe seems to have offset its specific area and micro/ mesopore volume losses with the high concentration of iron species on its surface because its adsorption capacity was similar to those materials with higher specific area and micro/mesopore volume (AC and AC3).The iron-impregnated materials AC5Fe and AC3Fe showed a greater slope of the kinetic curve and, consequently, a higher desulfurization capacity than the first group (AC, AC3, and AC7Fe), ranging from 52.1 to 66.3 mgH 2 S g −1 of adsorbent.These results are approximately 37 and 74% higher than the one obtained for AC, respectively.Although a deeper study to identify what iron species were formed on the AC surface upon thermal treatment was not carried out, these results confirm our expectation regarding boosting the activity of the AC porous surface.Moreover, this suggests that the higher surface area and micro/mesopore volume of iron-containing AC3Fe were not the only reasons for a superior H 2 S removal performance because AC3, which also showed a gain in the specific area and micro/mesopore volume, had more discreet performance in desulfurization, reaching only 40.0 mgH 2 S g −1 .In addition, the material AC5Fe, which presented a relatively lower specific  2) area (317.2 m 2 g −1 ) and high iron content, showed an H 2 S removal only ~ 27% lower than those obtained for AC3Fe, suggesting once again the relevant role of the iron species for H 2 S removal.Indeed, the adsorbent's carbonaceous matrix received an activation boost due to a simultaneous combined effect between improving its textural properties (increase in specific area and micro/ mesopore volume, Table 2), and formation of sites containing iron species (Fig. 4).As a first insight, iron species might have had their oxidizing power increased because the heat would have favored its attack on the activated carbon wall, causing a partial collapse of the carbonaceous matrix (abrasion, appearance of new pores, and interconnection of blind and closed pores to open ones) and formation of active sites as iron species found suitable places to settle down.Comparatively, a negative aspect concerning AC5Fe and AC7Fe is that, although they had high performances, the great weight loss depicted during thermal treatment (~ 75-79%, Table 1) makes their use unfeasible, that is, the only material that shows technical viability is AC3Fe.According to the literature, the adsorption capacity of activated carbons can vary intensely but usually ranges from 20 to 40 mg g −1 for ordinary commercial activated carbon and can reach 210 mg g −1 or more for modified ones (Yan et al. 2002;Sitthikhankaew et al. 2011).However, literature usually reports adsorption capacities obtained in ous flow fixed bed columns.Therefore, due to the special confinement conditions achieved in the static system, the adsorption capacity of the materials studied in this work may be lowered if evaluated in a continuous flow system.Despite that, the analytical method worked well for comparative purposes between modified activated carbon and non-modified one submitted to the same experimental conditions.
Table 3 shows the results obtained from the singlefactor ANOVA analyses considering the performances of the modified and non-modified activated carbons.The first analysis was performed with a complete dataset containing the five materials.The obtained p-value 0.005155 (< < 0.05) allows to reject the null hypothesis and confirm the distinction between the adsorption capacities.In order to evaluate the behavior of the material AC5Fe (pink line in Fig. 8), a second dataset was considered without the AC3Fe material and the p-value increased to 0.055737 (> 0.05), indicating that, statistically, the adsorption capacity of the AC5Fe material cannot be clearly distinct from others.In fact, the p-value is slightly higher than the established limit of 0.05.A third analysis was performed without samples AC3Fe and AC5Fe (not included in Table 3) and the p-value found was 0.941795 (> > 0.05), confirming that the materials AC, AC3 and AC7Fe presented practically the same performance for H 2 S removal.
In order to get an insight into the adsorption process and mainly the different performances displayed by the studied materials, three kinetic models were used to assess the fit of the experimental data.Figure 9 shows the result of fitting the nonlinear models to the studied materials.All materials had their data fitted to three kinetic models and showed good linear correlation (R 2 ).However, the pseudo-first order (PFO) model was the one that showed the closest proximity between experimentally determined adsorption capacity (q exp ) and calculated adsorption capacity (q calc ).Bangham's model, a derivation of the PFO model, was included to verify if the microporosity of the materials could have interfered with the adsorption process.Additionally, the pseudo-second order model (PSO) was used to evaluate if the presence of iron adsorption sites could have favored a second-order mechanism.Table 4 presents the kinetic parameters obtained for all the materials; for PFO and Bangham's models, all q exp were very close to q calc .Conversely, the PSO model showed a significant deviation between q exp and q calc , presenting a trend toward a steeper slope of the curve and a supposedly higher adsorption capacity, which is inconsistent with the observed results.The depicted behavior indicates that, despite the intense thermal treatment performed on the ACs, especially in the presence of iron species, the adsorption mechanism did not change significantly for each modified material, and the PFO model satisfactorily represents the adsorption dynamics, that is, a mechanism predominantly based on the formation of multilayers of H 2 S on the adsorbent's surface.Notwithstanding, the good fitting presented by Bangham's model may indicates the important role of the micropores comprising each material, which represents 68-71 wt% of the total specific area (Table 2).
Although the PFO model is closer to a physisorption process, literature suggests that ordinary activated carbon surfaces contain oxygen-based groups that allow the formation of highly exothermic bonds regarding complexes between H 2 S and C-O, indicating that this process can be ruled by chemisorption (Kazmierczak-Razna 2015;Li et al. 2016).The results presented in Table 4 suggest that the non-modified AC material follows the same trend as the modified ones, possibly as an effect of coexistence of chemi-and physisorption processes due to the intrinsic characteristics of this material, as aforementioned in XRD discussion, due to the presence of oxygen-based species.In the case of AC3Fe, AC5Fe and AC7Fe, it is possible to infer that the thermal  4 Kinetic parameters obtained for the ACs studied using three different models q exp = experimental adsorption capacity (mg g −1 ) q calc = calculated adsorption capacity (mg g −1 ) k 1 and k b = rate constant for PFO and Bangham's models, respectively, (min −1 ) k 2 = rate constant for PSO model (mg g −1 min −1 ) PFO pseudo-first-order model, PSO pseudo-second-order model

Sample ID PFO model PSO model
Bangham's model q exp q calc k 1 (× 10 -2 ) R 2 q calc k 2 (× 10 -4 ) R 2 q calc k b (× 10 treatment in the presence of iron species and an oxidizing atmosphere (air), which caused a more drastic oxidation, led to the increase of stronger and non-dissociative adsorption sites for H 2 S chemisorption (Shen et al. 2019), supporting a scenario of coexistence of physisorption and chemisorption.Table 4 also reveals that the rate constant value for the adsorption involving the AC material is one of the lowest (k 1 = 1.65 min −1 ), almost as low as that obtained for the AC7Fe (1.58 min −1 ), which is a highly mineralized material, and practically half of that obtained for the AC3Fe material.Actually, this is opposite what would be expected if the reaction rate were dependent only on the type of adsorption site (i.e., weak attachment-physisorption; strong attachment-chemisorption); this seemingly contradictory information is explained by the fact that the adsorption kinetic on the studied ACs is ruled by both type of adsorption sites (presence of active species spread on its surface) and textural characteristics (mainly, specific area and pore volume).The gain in specific area and volume of pores containing iron active sites for H 2 S adsorption are responsible for the higher adsorption rate of the modified materials.

Conclusion
Although activated carbon is one of the most known materials for adsorption purposes, surface modification is yet an alternative to improve its performance in separation processes.The modified materials showed higher specific area, pore volume, and iron species dispersed over its surface, which led to a significant increase in its adsorption capacity compared to the ordinary AC.The thermal treatment associated to iron(III) chloride allows mitigate the abrupt mineralization of the activated carbons and favored the formation of controlled microporosity.Both physical and chemical sorption may coexist until the saturation of the adsorbent's surface.
The static batch system used to evaluate H 2 S adsorption showed potential to be used as a technique for determining the adsorption capacity of an adsorbent material.However, the experimental arrangement requires adjustments to increase the level of reliability of the data collected.A clear advantage of the proposed method is to allow the use of a standard synthetic biogas mixture, which excludes fluctuations in the concentrations of gases that form biogas. Also, the use of standard gas mixture in a continuous flow column could present an impeditive cost for multiple long-term analysis.On the other hand, the proposed system is sensitive to the dilution step with the inert solid, to the transfer of gas aliquots between flasks, and moisture from air due to inefficient evacuation.Nevertheless, this analytical method may be an alternative to measure total adsorption capacity of materials with potential for desulfurization.

Fig. 1
Fig. 1 Schematic representation of the steps considered in the materials assessment during adsorption of H 2 S from synthetic biogas

Fig. 7 Fig. 8
Fig. 7 (a) GC chromatogram of the three-component gaseous standard mixture plus water from room moisture and their retention times.The inset shows the hydrogen sulfide peak.(b) Decrease in the area ratio evidencing the high capacity of H 2 S adsorption on studied ACs

Fig. 9
Fig. 9 Fittings of nonlinear models pseudo-first-order, pseudo-second-order and Bangham's model performed in each experimental kinetic curve

Table 1
Iron loading estimate in the ACs after mass loss observed during thermal treatment a Obtained by gravimetric analysis b Estimated according to mass loss during thermal treatmentThe initial iron content was 28.8 ± 0.20 mg Fe g −1 by UV-Vis spectroscopy

Table 2
Textural parameters obtained for the materials, including micropore, external and BET specific areas, average pore size and average pore volume Determined by a BET method, b de Boer method and c BJH method

Table 3
Single-factor ANOVA analysis for the materials performance results