Application of Sludge-Based Activated Carbons for the Effective Adsorption of Neonicotinoid Pesticides

: The amount of sludge produced in wastewater treatment plants (WWTPs) has increased over the years, and the methods used to reduce this waste, such as incineration, agricultural use, or disposal in landﬁlls, cause problems of secondary pollution. For this reason, it is necessary to ﬁnd sustainable and low-cost solutions to manage this waste. Additionally, emerging and priority pollutants are attracting attention from the scientiﬁc community as they can generate health problems due to inadequate removal in conventional WWTPs. In this work, a pharmaceutical industry sludge was used as a precursor in the synthesis of four activated carbons (ACs) using different activating agents (ZnCl 2 , FeCl 3 · 6H 2 O, Fe(NO 3 ) 3 · 9H 2 O, and Fe(SO 4 ) 3 · H 2 O), to be used for the removal by adsorption of three neonicotinoid pesticides included in latest EU Watch List (Decision 2018/840): acetamiprid (ACT), thiamethoxam (THM), and imidacloprid (IMD). The prepared ACs showed micro–mesoporous properties, obtaining relatively slow adsorption kinetics to reach equilibrium, but despite this, high values of adsorption capacity (q e ) were obtained. For example, for AC-ZnCl 2 (S BET = 558 m 2 /g), high adsorption capacities of q e = 128.9, 126.8, and 166.1 mg/g for ACT, THM, and IMD, respectively, were found. In most cases, the adsorption isotherms showed a multilayer proﬁle, indicating an important contribution of the mesoporosity of the activated carbons in the adsorption process.


Introduction
Sewage sludge is defined as the unavoidable waste generated in wastewater treatment processes and consists of a highly heterogeneous mixture of water and solids. The solid phase consists mainly of heavy metals, organic pollutants, and pathogenic microorganisms [1]. The composition of the sludge varies considerably depending on its source, and it can therefore be classified mainly into two types: if it comes from domestic wastewater, it is known as urban sludge, and if, on the other hand, it comes from industrial water, it is known as industrial sludge. This latter is notable for its high organic matter content and is generally the most suitable precursor to produce carbonaceous adsorbents [2,3].
In the EU, around 13.0 million tons of sewage sludge (in dry matter) are produced annually [4]. The amount of sludge produced has increased enormously with the development of industrialization and urbanization [5,6], and it is considered to be a highly alarming wastewater treatment problem. Different methods have been selected for the removal of sewage sludge, such as disposal in landfills, incineration, and application in agriculture as a fertilizer; but each of these options causes problems of secondary pollution [7,8]. Furthermore, due to increasingly strict legislation (Directive 2008/98/EC) [9,10], sustainable solutions need to be found to manage this waste at low cost. As an alternative, a valorization route has been proposed within the concept of a circular economy; this route consists of using the sludge as a low-cost raw material for its transformation into at room temperature, followed by drying in an oven at 105 • C for 24 h. Following the activation process, the material was pyrolyzed in a vertical quartz reactor at 800 • C for 2 h, using a N 2 flow rate of 100 mL/min and a heating rate of 10 • C/min. Afterward, the activated carbon was washed first with a 5 M HCl solution and then with deionized water. The solid was again dried in an oven at 105 • C for 24 h to continue with a second pyrolysis using the same conditions as described previously. Finally, the resulting activated carbons were ground and sieved (Ø = 50-250 µm).

Characterization of the Sludge
In this work, the part of interest of the sludge was the solid fraction, from which the activated carbons were synthesized. Thus, the chemical oxygen demand (COD) of the aqueous fraction was determined using a PF-11 Macherey-Nagel photometer (Düren, Germany). The total (TS), fixed (FS), and volatile (VS) solids of the sludge were measured according to the Standard Methods [28]. The elemental analysis of the solid was carried out in a LECO CHNS-932 analyzer (Mönchengladbach, Germany), while the chemical composition was determined by X-ray fluorescence spectroscopy (XRF) using PANanalytical Axios equipment (Malvern, UK).

Characterization of the Activated Carbons
The textural properties of the activated carbons were determined by N 2 adsorptiondesorption at 77 K in a Micromeritics ASAP 2020 system. In the same way, the surface chemistry properties were determined by Fourier transform infrared spectroscopy (FT-IR) using a Thermo Nicolet AVATAR 360 spectrophotometer in the wavenumber range of 400 to 4000 cm −1 . The thermogravimetric studies of the solids (TGA) were carried out in a PerkinElmer STAR 6000 device, using a temperature range of 25 to 900 • C, a heating rate of 10 • C/min, and a nitrogen volumetric flow rate of 50 mL/min. The elemental analyses were performed in a LECO CHNS-932 analyzer. Finally, the isoelectric point (pH IEP ) of the samples was measured by zeta potential measurements of a dispersion of 0.05 g of carbon (Ø = 10-20 µm) in 20 mL of deionized water using a Zetasizer Nano ZS (Malvern Instruments, Ltd., Malvern, UK).

Batch Adsorption Studies
Batch adsorption tests were carried out using a LabMate orbital shaker at constant temperature (25 • C) and a shaking rate of 250 rpm. For the kinetic studies, 7.5 mg of adsorbent and 25 mL of each pesticide solution (carbon dose of 0.3 g/L) were maintained under constant agitation. Samples were collected at the determined time intervals until the equilibrium time was reached. The initial adsorbate concentration used in the experiments was 50 mg/L. All the equilibrium adsorption studies were carried out under the same operating conditions, but varying the mass of adsorbent between 1.5 and 37.5 mg (carbon dose ranging from 0.06 to 1.5 g/L). When the equilibrium time was reached, the equilibrium adsorption capacity, q e (mg/g), was calculated using Equation (1): where q e (mg/g) is the equilibrium adsorption capacity, C 0 (mg/L) is the initial pesticide concentration, C e (mg/L) is the pesticide concentration at equilibrium, V (L) is the volume of solution, and W (g) the weight of adsorbent.

Analytical Procedure
The pesticide concentration in aqueous samples was analyzed by a high-performance liquid chromatography (HPLC) technique, using Varian ProStar equipment with a "diode array" detector. A Teknokroma column (25 mm × 0.46 mm; 5 µm) was used as the stationary phase, and as the mobile phase, a 70-30% (v/v) acetonitrile/75 mM acetic acid solution was used. The THM, IMD, and ACT concentrations were determined at a wavelength of 270 nm, using a flow rate of 0.8 mL/min.

Characterization of the Sludge
The macroscopic properties of the sludge are reported in Table 1. It is worth highlighting the high COD value of the aqueous fraction (55.36 g/L), which is to be expected since it is an industrial sludge. In addition, both the COD and volatile solids values indirectly indicate the organic matter content of the waste; approximately 80% of the total solids are volatile. These solids are the most relevant, as the greater the amount of organic matter present in the precursor, the higher the efficiency obtained in the synthesis of activated carbon [29]. Another noteworthy value is the Ca content (7.25%), as this is a value that stands out when compared to the rest of the metals present in the sludge. This may be due not only to the origin of the sludge but also to the products used in the wastewater treatments applied in the industry, involving the use of calcium hydroxide, commonly known as slaked lime, which is the most-used calcium-based product in WWTPs.

Characterization of the Activated Carbons
The textural properties of the synthesized activated carbons were studied through N 2 adsorption-desorption analysis. The isotherms and pore size distributions of these materials are shown in Figure 1a,b. The textural parameters, specific surface area (S BET ), external surface area (S ext ), total pore volume (V Total ), and micropore volume (V Micro ) are reported in Table 2. According to the International Union of Pure and Applied Chemistry (IUPAC) classification, all the samples exhibited type IV isotherms, characteristic of mainly mesoporous materials. At medium pressure, due to capillary condensation, the desorption branch did not match with the adsorption isotherm, causing a hysteresis loop, which demonstrates the presence of wider pores (meso-and macropores) in the carbon structure [30]. In addition, the hysteresis loops presented by the activated carbons showed two parallel and almost horizontal branches (except for the case of AC-Fe 2 (SO 4 ) 3 ), which can be classified as H4-type hysteresis loops, according to the IUPAC classification.  According to the International Union of Pure and Applied Chemistry (IUPAC) classification, all the samples exhibited type IV isotherms, characteristic of mainly mesoporous materials. At medium pressure, due to capillary condensation, the desorption branch did not match with the adsorption isotherm, causing a hysteresis loop, which demonstrates the presence of wider pores (meso-and macropores) in the carbon structure [30]. In addition, the hysteresis loops presented by the activated carbons showed two parallel and almost horizontal branches (except for the case of AC-Fe2(SO4)3), which can be classified as H4-type hysteresis loops, according to the IUPAC classification.
Sludge is a heterogeneous waste and, therefore, generates activated carbons with heterogeneous textural properties-that is, with meso-microporous character. This phenomenon was previously observed in activated carbons synthesized from biomass [31,32]. In the case of AC-ZnCl2, it can be observed that part of the pore filling occurs at low pressure, since it could be considered a mainly microporous activated carbon. In contrast, in the case of AC-Fe2(SO4)3 carbon, a significant contribution of macroporosity was observed, showing capillary condensation at high pressure values.
The pore size distributions of the materials were evaluated by the Barrett-Joyner-Halenda (BJH) method. It can be seen that all the samples showed a pore distribution mainly centered in the micro-mesoporosity range, except for the AC-Fe2(SO4)3 sample, which showed wider pores (~50 nm). Table 2 shows that the AC-ZnCl2 carbon showed the highest specific surface area. The order according to development of the specific surface area for the different synthesized carbons was as follows: AC-ZnCl2 > AC-FeCl3 > AC-Fe(NO3)3 > AC-Fe2(SO4)3. Regarding the pore volume, AC-Fe2(SO4)3 carbon showed a practically negligible micropore volume, whereas with the AC-ZnCl2 sample, the opposite was observed, i.e., the mi- Sludge is a heterogeneous waste and, therefore, generates activated carbons with heterogeneous textural properties-that is, with meso-microporous character. This phenomenon was previously observed in activated carbons synthesized from biomass [31,32]. In the case of AC-ZnCl 2 , it can be observed that part of the pore filling occurs at low pressure, since it could be considered a mainly microporous activated carbon. In contrast, in the case of AC-Fe 2 (SO 4 ) 3 carbon, a significant contribution of macroporosity was observed, showing capillary condensation at high pressure values.
The pore size distributions of the materials were evaluated by the Barrett-Joyner-Halenda (BJH) method. It can be seen that all the samples showed a pore distribution mainly centered in the micro-mesoporosity range, except for the AC-Fe 2 (SO 4 ) 3 sample, which showed wider pores (~50 nm). Table 2 shows that the AC-ZnCl 2 carbon showed the highest specific surface area. The order according to development of the specific surface area for the different synthesized carbons was as follows: AC-ZnCl 2 > AC-FeCl 3 > AC-Fe(NO 3 ) 3 > AC-Fe 2 (SO 4 ) 3 . Regarding the pore volume, AC-Fe 2 (SO 4 ) 3 carbon showed a practically negligible micropore volume, whereas with the AC-ZnCl 2 sample, the opposite was observed, i.e., the micropore volume corresponded to practically half of the total pore volume. This behavior was corroborated by the N 2 adsorption-desorption isotherms, as discussed before.
The FT-IR spectra of the activated carbons are depicted in Figure 2a. A broad band located at~3400 cm −1 can be attributed to O-H stretching vibration due to the presence of water in the carbon samples. The band close to 2900 cm −1 , very strong for the AC-ZnCl 2 and AC-FeCl 3 samples, is associated with stretching vibration of aliphatic C-H bonds, whereas the peak at 2331 cm −1 , observed only in the AC-Fe(NO 3 ) 3 and AC-Fe 2 (SO 4 ) 3 spectra, indicates the presence of ketone functional groups. The band close to 1560 cm −1 corresponds to C=C stretching vibration of the aromatic rings [33]. Thus, the bands associated with alcohols, phenols, esters, acids, and ethers were observed at a wavelength of 1100 cm −1 . Finally, the bands attributed to out-of-plane bending vibrations of O-H and C-H were observed between 500 and 750 cm −1 [34].  Figure 2b shows the thermogravimetric (TG) analysis results of the synthesized activated carbons, which evaluated the thermal stability of the solids. In all the samples, a weight loss was observed at temperatures around 100 • C, which is attributed to the evaporation of the adsorbed water. The weight loss observed between 200 and 400 • C may be due to the decomposition of volatile organic and inorganic compounds, and that in the range between 400 and 600 • C may be attributed to the decomposition of high-molecular-weight compounds. At temperatures around 600 • C, a more abrupt weight loss was observed for activated carbons synthesized with Fe(NO 3 ) 3 and Fe 2 (SO 4 ) 3 , whereas this drop was more stable for activated carbons synthesized from impregnating agents with chloride, i.e., ZnCl 2 and FeCl 3 . Therefore, the weight loss observed at 600 • C can be mainly attributed to the thermal decomposition of the used activating agent.
The elemental analysis results of the samples and the corresponding isoelectric point (pH PIE ) values are collected in Table 3. Generally, activated carbons are mainly composed of fixed carbon and inorganic compounds. As can be seen in Table 3, the fixed carbon contents (C, wt.%) were 55.60, 59.77, 33.51, and 31.71% for AC-ZnCl 2 , AC-FeCl 3 , AC-Fe(NO 3 ) 3 , and AC-Fe 2 (SO 4 ) 3 , respectively-values significantly lower than those found in activated carbons obtained from fossil fuels [25]. Even when the material used as a precursor has a high organic matter content, at a high carbonization temperature, the organic matter decomposes; this causes a decrease in the fixed carbon percentage, transforming it into ash. Therefore, the higher the carbonization temperature, the lower the percentage of fixed carbon found in the final material, indicating greater carbonization of organic matter. For this reason, the carbonization temperature of sludge in industry usually varies between 400 and 600 • C, since in this temperature range, the product obtained is often suitable for use as an alternative fuel due to its high organic content [35]. Also, it is worth highlighting the S content of 13.15% measured in AC-Fe 2 (SO 4 ) 3 , indicating that impregnation with Fe 2 (SO 4 ) 3 was efficiently carried out. On the other hand, according to the literature, washing with HCl promotes the production of acidic adsorbents [36], although it is necessary to correctly remove ash residues that clog the pores of the final materials. Furthermore, independently of the washing conditions, the use of activating agents such as FeCl 3 and ZnCl 2 [37,38] conditioned the preparation of predominantly acidic carbon materials, as can be observed from their isoelectric point (pH IEP ) values.

Kinetic Adsorption Studies
The adsorption kinetics of the pesticides onto the synthesized activated carbons are depicted in Figure 3a-d.
carbonization of organic matter. For this reason, the carbonization temperature of sludge in industry usually varies between 400 and 600 °C, since in this temperature range, the product obtained is often suitable for use as an alternative fuel due to its high organic content [35]. Also, it is worth highlighting the S content of 13.15% measured in AC-Fe2(SO4)3, indicating that impregnation with Fe2(SO4)3 was efficiently carried out. On the other hand, according to the literature, washing with HCl promotes the production of acidic adsorbents [36], although it is necessary to correctly remove ash residues that clog the pores of the final materials. Furthermore, independently of the washing conditions, the use of activating agents such as FeCl3 and ZnCl2 [37,38] conditioned the preparation of predominantly acidic carbon materials, as can be observed from their isoelectric point (pHIEP) values.

Kinetic Adsorption Studies
The adsorption kinetics of the pesticides onto the synthesized activated carbons are depicted in Figure 3a-d.  Slow adsorption rates were observed for the three pesticides, achieving equilibrium times ranging from 24 to 144 h. As can be seen in Figure 3, in the first stage, rapid adsorption occurred because adsorption sites on the surface are easily accessed and occupied by the micropollutant; in the next stage, the diffusional effect comes into play (conditioned by the steric hindrance caused by the adsorbate molecules), due to the relatively high molecular size of the pesticides used compared to the medium available pore size of the activated carbons, requiring in this stage more time to diffuse until reaching equilibrium. According to this, longer equilibrium times were found when the molecular volume of the pesticides used was larger, following the order ACT (266 Å 3 ) < IMD (271 Å 3 ) < THM (303 Å 3 ). For this reason, the slowest equilibrium time (144 h) was obtained for the larger pesticide, THM.
In the same way, very high adsorption removal values were obtained at equilibrium time. The best results were achieved with the AC-ZnCl 2 sample. Generally, it was observed that better adsorption capacities were obtained with those activated carbons that showed a higher specific surface area (Table 2); therefore, it is clear that the textural properties of the materials played a key role in their adsorption performance. An exception was observed for THM removal using the AC-Fe(NO 3 ) 3 and AC-Fe 2 (SO 4 ) 3 samples, where a simultaneous contribution of adsorption and reaction phenomena was observed (Figure 3c,d).
In both cases, a first stage related to the adsorption process and a second stage where the adsorption process occurred simultaneously with a possible degradation-mediated by-reaction were observed; this phenomenon may be due to the high iron content found in these carbonaceous materials.
Pseudo-first-order and pseudo-second-order kinetic models were applied to assess the experimental kinetic data. The pseudo-first-order (PPO) kinetic model, also known as the Lagergren equation, assumes that physisorption limits the adsorption rate on the adsorbent particles and can be expressed by the following equation: where k 1 (h −1 ) is the first-order reaction rate equilibrium constant, and q (mg/g) and q e (mg/g) denote the adsorption capacity at any time and at the equilibrium time, respectively. On the other hand, the pseudo-second-order (PSO) kinetic model assumes that chemisorption is the rate-limiting step of the adsorption process, and it is expressed as follows: where k 2 (g/mg·h) is the rate constant of the pseudo-second-order kinetic model. As shown in Table 4, considering the experimental and estimated adsorption capacities and R 2 values, both the pseudo-first-order and pseudo-second-order kinetic models well described the adsorption of pesticides onto the synthesized activated carbons. This suggests that physical and chemical interactions simultaneously contribute to and control the adsorption of the pesticides on the surfaces of the activated carbons employed. Generally, physical adsorption, or physisorption, involves van der Waals forces, while chemical adsorption, or chemisorption, occurs when electrons are shared or transferred through covalent bonds [39].

Isotherm Adsorption Studies
The adsorption isotherms of the three studied pesticides onto the synthesized adsorbents are shown in Figure 4a-d. The initial concentration of micropollutant was fixed at 50 mg/L, whereas the dose of adsorbent was varied within the range of 0.05-1.5 mg/mL. From the results, the best adsorption capacities for the three contaminants were found for AC-ZnCl 2 : q e = 128.9, 126.8, and 166.1 for ACT, THM, and IMD, respectively. As can be observed, all the adsorption isotherms showed a multilayer profile. Thus, they could be classified as type S-3 and S-4, according to the classification established by Giles et al. [40]. The S-type isotherms are characterized by a steep slope, and as the adsorbate concentration increases, the adsorption capacity approaches a first plateau. The first plateau is attributed to saturation caused by the termination of the first monolayer. Then, as the adsorbate concentration increases, a vertical rearrangement of the adsorbed molecules occurs, and more active sites become available for adsorption [41]. Therefore, the newly available adsorption sites lead to a further increase of the adsorption capacity, and a second plateau occurs, which is attributed to a second layer. These S-type isotherms may also be the result of a competitive effect between the solvent and the adsorbate towards the active sites available for adsorption [42].

. Isotherm Adsorption Studies
The adsorption isotherms of the three studied pesticides onto the synthesized adsorbents are shown in Figure 4a-d. The initial concentration of micropollutant was fixed at 50 mg/L, whereas the dose of adsorbent was varied within the range of 0.05-1.5 mg/mL. From the results, the best adsorption capacities for the three contaminants were found for AC-ZnCl2: qe = 128.9, 126.8, and 166.1 for ACT, THM, and IMD, respectively. As can be observed, all the adsorption isotherms showed a multilayer profile. Thus, they could be classified as type S-3 and S-4, according to the classification established by Giles et al. [40]. The S-type isotherms are characterized by a steep slope, and as the adsorbate concentration increases, the adsorption capacity approaches a first plateau. The first plateau is attributed to saturation caused by the termination of the first monolayer. Then, as the adsorbate concentration increases, a vertical rearrangement of the adsorbed molecules occurs, and more active sites become available for adsorption [41]. Therefore, the newly available adsorption sites lead to a further increase of the adsorption capacity, and a second plateau occurs, which is attributed to a second layer. These S-type isotherms may also be the result of a competitive effect between the solvent and the adsorbate towards the active sites available for adsorption [42].  Correlation of the experimental equilibrium data to empirical models is essential for the practical interpretation of scientific results. In this work, the Langmuir, Freundlich, Sips, and Guggenheim-Anderson-de Boer (GAB) isotherm models were used to fit the  Correlation of the experimental equilibrium data to empirical models is essential for the practical interpretation of scientific results. In this work, the Langmuir, Freundlich, Sips, and Guggenheim-Anderson-de Boer (GAB) isotherm models were used to fit the experimental adsorption data. The Langmuir isotherm model supposes that the adsorption process takes place in a monolayer, that the solid surface is energetically homogeneous, and that the molecules are adsorbed on well-defined sites, i.e., without migration of the adsorbate from one active site to another [43]. This model is expressed by the following equation: where q e is the equilibrium adsorption capacity (mg/g), C e is the equilibrium concentration (mg/L), q sat is the maximum adsorption capacity in the monolayer (mg/g), and b is a constant related to the affinity between the adsorbate and adsorbent (mg/L). The Freundlich isotherm is an empirical model used when the adsorption mechanism appears to be mixed and the material surface can be considered energetically heterogeneous. As a consequence, this model can be applied to multilayer adsorption isotherms [44]. The Freundlich equation is as follows: where q e is the equilibrium adsorption capacity (mg/g), k f (L/g) is the Freundlich constant related to the adsorption capacity of the adsorbent, and 1/n f indicates the intensity of the adsorption process or the surface heterogeneity; when the value of n f is >1, the adsorption conditions can be considered as favorable.
The Sips model was developed to give the Freundlich model a finite limit at high equilibrium concentrations. The difference from the Langmuir isotherm model is an additional parameter, n; when n = 1, the Sips model simplifies to the Langmuir equation. The parameter n can be defined as a measure of the surface heterogeneity: this heterogeneity may be due to the adsorbent surface, the adsorbate, or both. The value of n is usually higher than unity, and the greater the value, the more heterogeneous the system. The Sips model can be expressed by the following equation: q e = q sat · (b · C e ) 1/n 1 + (b · C e ) 1/n (6) where b (L/mg) and 1/n are temperature-dependent Sips parameters, and q sat (mg/g) can be taken as a function (or not) of temperature. Finally, the Guggenheim-Anderson-De Boer (GAB) isotherm model is an extension of the Langmuir isotherm model that considers multilayer adsorption by assuming that each adsorbed molecule provides a site for the second one. The expression of the GAB model is as follows, assuming that only two layers are formed in the adsorption process: where K 1 and K 2 (L/mg) are equilibrium constants attributable to the first and second layer, respectively, and q m (mg/g) is the maximum adsorption capacity in the first layer. The isotherm model parameters for adsorption onto the synthesized activated carbons are shown in Tables 5 and 6. The best fit to the experimental data of the adsorption isotherms was found for the Freundlich model, with higher values of the correlation coefficient (R 2 ); this behavior was expected since it is a model suitable for fitting heterogeneous adsorption systems with an exponential distribution of active sites [45]. Also, it is worth noting that for the adsorption onto AC-Fe 2 (SO 4 ) 3 , the Sips isotherm model fitted equally as well as Freundlich model. Regarding the parameter n f of the Freundlich equation, generally, values higher than unity are indicative that significant adsorption takes place at low concentration values (favorable adsorption isotherms), since the increase in the adsorption capacity with aqueous concentration becomes less significant at a higher concentration (n f values lower than unity). Most of the obtained n f values were greater than unity, thus verifying good adsorption onto a heterogeneous carbon surface [46,47].
Although the Freundlich isotherm model achieved the best results, it did not provide a complete reproduction of the experimental data, mainly due to the complexity of reproducing multilayer adsorption isotherms and the large plateau region exhibited between the layers. Table 7 shows the adsorption capacities obtained in previous studies for the adsorption of the pesticides studied in this work, i.e., acetamiprid and imidacloprid, with adsorbents synthesized from different biomass sources. The best results in this work were obtained with the AC-ZnCl 2 activated carbon (S BET = 558 m 2 /g), with adsorption capacities of 128.9, 126.8, and 166.1 mg/g for ACT, IMD, and THM, respectively, at 24 h. When comparing the results obtained for AC-ZnCl 2 with those shown in Table 7, it can be seen that a higher adsorption capacity was obtained with the pesticide IMD when using sugarcane bagasse as a precursor, while for the pesticide ACT, lower adsorption capacity values were obtained; hence, it can be concluded that the adsorption capacity varies considerably depending on the origin of the activated carbon. On the other hand, if S BET values are considered, excluding that for the activated carbon from peach stones, the values were equal to or even higher than that found for the AC-ZnCl 2 activated carbon; however, the adsorption capacity and equilibrium time values were lower, except those for the activated carbon from sugarcane bagasse. Therefore, if the referenced studies did not allow enough contact time between the adsorbent and adsorbate, it is possible that saturation of the activated carbon did not occur, which would explain the very fast equilibrium times. For the activated carbon from sugarcane bagasse, the adsorption capacity was the highest (q e = 313 mg/g), and it took 720 min to reach the equilibrium. Furthermore, it should be noted that no literature results were found using the pesticide THM, which shows the innovativeness of this work. Furthermore, AC-ZnCl 2 outperformed all previously published adsorbents tested on ACT and outperformed three of four published adsorbents for IMD. In conclusion, more research is needed to determine the performance of this kind of adsorbent in pesticide removal.

Conclusions
The use of sewage sludge as a low-cost precursor for the preparation of activated carbons to removal neonicotinoid pesticides, e.g., acetamiprid, thiamethoxam, and imidacloprid, is a sustainable way to solve the serious problem of sludge management while being in accordance with the concept known as a circular economy. The synthesized activated carbons showed heterogeneous textural properties (micro-mesoporous character)-a characteristic of activated carbons obtained from waste and perhaps also as a result of the heterogeneous nature of the precursor used. The maximum specific surface area (S BET = 558 m 2 /g) was obtained for the AC synthesized using ZnCl 2 as the impregnation agent. Relatively slow adsorption kinetics were obtained for all the materials as a result of the steric hindrance caused by the molecules adsorbed, that is, due to the relatively large molecular size of the pesticides in relation to the medium pore size of the activated carbons; thus, the equilibrium times were longer when THM was used, since this is the largest of the tested molecules. Despite this, very high equilibrium adsorption capacities were obtained: using AC-ZnCl 2 , q e values of 128.9, 126.8, and 166.1 mg/g for ACT, THM, and IMD, respectively, were obtained. Besides this, it is important to highlight the behavior observed in the removal of THM onto the AC-Fe(NO 3 ) 3 and AC-Fe 2 (SO 4 ) 3 activated carbons, since simultaneous contributions of adsorption and reaction phenomena were observed; this is possibly due to the high content in Fe of these materials. Regarding the adsorption isotherms, it is interesting that, in general, they showed a multilayer profile, which indicates an important contribution of the mesoporosity in the adsorption process.
For all these reasons, the use of sewage sludge to produce activated carbons has great potential as it is a sustainable solution for the management of this kind of waste, and the resulting activated carbons can be considered capable of competing with conventional adsorbent materials synthesized from fossil fuels.