Spent Coffee Grounds as an Adsorbent for Malathion and Chlorpyrifos—Kinetics, Thermodynamics, and Eco-Neurotoxicity

Coffee is one of the most popular beverages, with around 10.5 million tons manufactured annually. The same amount of spent coffee grounds (SCGs) might harm the environment if disposed of carelessly. On the other hand, pesticide contamination in food and biowaste is a rising problem. Because pesticides are hazardous and can cause serious health consequences, it is critical to understand how they interact with food biowaste materials. However, it is also a question if biowaste can be used to remediate rising pesticide residues in the environment. This study investigated the interactions of SCGs with the organophosphate pesticides malathion (MLT) and chlorpyrifos (CHP) and addressed the possibility of using SCGs as adsorbents for the removal of these pesticides from water and fruit extracts. The kinetics of MLT and CHP adsorption on SCGs fits well with the pseudo-first-order kinetic model. The Langmuir isotherm model best describes the adsorption process, giving the maximal adsorption capacity for MLT as 7.16 mg g−1 and 7.00 mg g−1 for CHP. Based on the thermodynamic analysis, it can be deduced that MLT adsorption on SCGs is exothermic, while CHP adsorption is an endothermic process. The adsorption efficiency of MLT and CHP using SCGs in a complicated matrix of fruit extracts remained constant. The neurotoxicity results showed that no more toxic products were formed during adsorption, indicating that SCGs are a safe-to-use adsorbent for pesticide removal in water and fruit extracts.


Introduction
Food biowaste represents an organic waste material derived from food processing, preparation, and serving activities, including agricultural by-products and kitchen scraps [1]. It is becoming a significant global issue, as the components of biowaste release CO 2 and other greenhouse gases in the process of breaking down, in addition to the economic and moral factors of throwing food away. Up to 30% of food biowaste globally is leftovers from restaurants, grocery stores, and households [2].
Coffee is one of the most popular beverages, with around 10.5 million tons produced yearly [3]. However, the same amount of spent coffee grounds (SCGs) could negatively affect the environment if disposed of improperly. When exposed to different environmental influences, coffee grounds release nitrogen and carbon dioxide into the atmosphere [4,5]. Nitrogen reacts with oxygen to form smog and ozone, which are both detrimental to the environment, while carbon dioxide is a greenhouse gas [6]. If disposed of in water, spent coffee grounds could increase acidity, thus decreasing oxygenation. These factors could lead to an increase in algae and bacteria growth [7]. To reduce the negative environmental The adsorption of pesticides on food biowaste materials has many important implications, leading to the possible accumulation of pesticides or indicating the potential applicability of biowaste in pesticide remediation.
This study investigated the potential of using SCGs as an adsorbent for malathion and chlorpyrifos removal from water and fruit extracts. Our goal was to use SCGs as received, with minimal additional treatment, and without an intense use of energy and/or chemicals and release of greenhouse gases. In this way, we aimed to investigate the possible accumulation of OPs in waste SCGs and assess the direct reuse of SCGs in compliance with contemporary environmental protection strategies [27]. First, the results of the physicochemical characterization of SCGs are presented. Then, CHP and MLT adsorption on SCGs in water is analyzed. The kinetics and thermodynamics of MLT and CHP adsorption on SCGs are discussed. The adsorption of OP from realistic samples (plant extract) on SCGs is also addressed to determine the possibility of using SCGs as an adsorbent and detect if the complex matrixes influence the process. In addition, the eco-neurotoxicity of the spiked fruit extracts was monitored during the remediation process to determine if more toxic products were formed, such as pesticides' oxo-forms. Finally, this method's feasibility and potential benefits in food processing are discussed. At the time of submitting this article, to the best of our knowledge, no research articles were dedicated to the investigation of remediation of MLT and CHP using SCGs as a sustainable material. Moreover, in this study, the eco-neurotoxicity testing of the samples before and after the adsorption process was performed for the first time.

Adsorbent Preparation
Coffee (purchased from the local market, 80% Arabica and 20% Robusta) was brewed (treated with boiling water) and left for 2 h at room temperature until it became to cold. Next, the coffee grounds were separated with filtration and left to dry at room temperature for 24 h. Then, in order to release the leftover moisture, the obtained spent coffee grounds were dried in the oven at 80 • C for 1 h and ground for 15 min, using an agate stone mortar. Finally, 100 mg of the material was rinsed with 50 mL of HCl, NaOH (Centrochem, Stara Pazova, Srbija), and H 2 O and eventually dispersed in 50 mL 50% EtOH (J.T. Baker, Phillipsburg, NJ, USA), obtaining the stock dispersion concentration of 2 mg mL −1 .

Adsorbent Characterization
The fractional sieving method was used to determine the particle size distribution in the SCG sample (sieve sizes 200, 100, and 63 µm; ROTH, Karlsruhe, Germany). For the investigation of samples' morphology and elemental composition, Scanning Electron Microscopy (SEM) and Energy Dispersive X-Ray Analysis (EDX) were performed using a scanning electron microscope PhenomProX (Thermo Fisher Scientific, Waltham, MA, USA). Fourier-transform infrared (FTIR) spectra were recorded using a Nicolet iS20 FT-IR spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The applied wavenumber range was from 4000 to 500 cm −1 , with 64 scans and 4 cm −1 . A TA Instruments SDT 2960 thermoanalytical device (TA Instruments, Inc. New Castle, DE, USA) was used for thermogravimetric analysis (TGA). The analysis was performed with a heating rate of 10 • C min −1 up to the temperature of 900 • C and under purging helium gas (Messer, Belgrade, Serbia). For the temperature-programmed desorption (TPD) analysis, the SCG sample was heated in a vacuum (starting pressure 1 × 10 −7 mbar), with a constant heating rate of 10 • C min −1 , from room temperature to 1000 • C (on the heater). Desorbed gaseous products were detected using a quadrupole mass spectrometer EXTORR XT300 (Extorr Inc., New Kensington, PA, USA).

Adsorption Experiments
A mixture of 1 mL SCG stock dispersion and the designated amounts of chlorpyrifos and malathion (Pestanal, Sigma Aldrich, Søborg, Denmark) stock solutions (made in  50 vol.% ethanol in water) was made to deliver the targeted concentration of adsorbent and OP. After that, the mixtures were put in a laboratory shaker and left for the specified period. Subsequently, they underwent centrifugation at 14,500 rpm, with their supernatant filtered through a nylon membrane. Ultra-Performance Liquid Chromatography (UPLC) analysis was then conducted to determine the concentrations of CHP and MLT with a Waters ACQUITY UPLC system and a Photodiode array (PDA) detector managed by Empower 3 software. An ACQUITY UPLC™ BEH C18 column (1.7 µm, 100 mm × 2.1 mm) was used under isocratic conditions of 10% acetonitrile (J.T. Baker, Phillipsburg, NJ, USA) in water (v/v) for mobile phase A and pure acetonitrile for mobile phase B. The eluent flow rate was 0.2 mL min −1 in all cases, with an injection volume of 5 µL. The mobile phase used for chlorpyrifos contained 20% A and 80% B, and 40% A and 60% B for malathion. The retention time for malathion was 3.2 min, and for chlorpyrifos, it was 2.7 min. Both OPs were detected at 200 nm. Control experiments were performed identically but without an adsorbent.

Plant Extracts Preparation
Lemon juice and mint extract were used as food samples. Lemon juice was made by squeezing one lemon (75 g) and diluting it with 500 mL of tap water (pH = 4.5) before adding MLT to the desired concentration. Next, the sample was filtered. The mint extract was prepared by mixing 7 g of Mentha spicata leaves with 45 mL of 50% ethanol and leaving it for 72 h at room temperature; this extract was then filtered and diluted with 200 mL of 50% ethanol (pH = 6.0). CHP was then added in the required amount, and the resulting solution was further filtered through a nylon filter. The prepared extracts were used for sample analysis.

Eco-Neurotoxicity Assessment
The physiological effects of the treated solutions were analyzed using AChE inhibition measurements. By employing modified Ellman's procedure [28,29], 2.5 IU commercially purified AChE (Sigma Aldrich, Taufkirchen, Germany) from an electric eel was exposed to the OP solutions in 50 mM phosphate buffer, pH = 8.0, at 37 • C in a 0.650 mL final volume. The combination of acetylcholine-iodide (AChI, Sigma Aldrich, Taufkirchen, Germany) and DTNB (Sigma Aldrich, Taufkirchen, Germany) as a chromogenic reagent triggered the enzymatic reaction. The reaction was allowed to proceed for 8 min before being stopped with 10% sodium dodecyl sulfate (SDS). Thiocholine, the reaction product, reacts with DTNB, forming 5-thio-2-nitrobenzoate, whose optical adsorption was then read at 412 nm. The enzyme concentration was kept constant and set to produce an optimal spectrophotometric signal. The physiological effects were quantified as the AChE inhibition, which is given as follows: where A 0 and A stand for the AChE activity in the absence of OP and the one measured after exposure to a given OP.

Results and Discussion
3.1. Physicochemical Characterization of SCGs 3.1.1. Particle Sizes, Morphology, and Chemical Composition The particle size distribution, determined by the fractional sieving method (Figure 1a), shows that the largest fraction of SCG particles is in the 100-200 µm range. Particle sizes in the order of 10 2 µm generally agree with previous reports of SCG properties, but the SCG particles reported here are somewhat smaller than in earlier reports [30,31], likely due to the grounding step we performed. We note that particle sizes of the same order of magnitude are frequently reported for biobased carbon materials used as adsorbents (see, for example, [32]). However, carbon materials typically have large specific surface areas (SSA) due to the developed internal pore structure. SSA is usually determined by gas 1. The SEM micrograph showed that SCGs have a sponge-like, porous, heterogeneous structure ( Figure 1(b1-b3)). Similar morphologies and large pores with diameters of approx. 30 μm were reported previously [34,35]. It can be seen that the cavities are randomly distributed on the materials' surface, with different size openings. SEM images were taken without deposition of the conductive layer on the SCGs sample, resulting in charge accumulation at higher magnifications. However, smooth pore walls can be observed without bright spots, indicating heavy elements in the sample. The surface 3D SEM reconstruction confirms these findings (Figure 1(c1,c2)). At a lower magnification (×2000), a rough surface profile is seen, with macropores entering the interior of the particles and contributing to the overall internal pore structure. However, higher magnification at the pore wall (×5000) indicates a rather smooth surface.  To investigate the morphology of SCGs, SEM was used. In addition, the chemical composition of the SCGs was investigated using EDX. The results are presented in Figure 1. The SEM micrograph showed that SCGs have a sponge-like, porous, heterogeneous structure ( Figure 1(b1-b3)). Similar morphologies and large pores with diameters of approx. 30 µm were reported previously [34,35]. It can be seen that the cavities are randomly distributed on the materials' surface, with different size openings. SEM images were taken without deposition of the conductive layer on the SCGs sample, resulting in charge accumulation at higher magnifications. However, smooth pore walls can be observed without bright spots, indicating heavy elements in the sample. The surface 3D SEM reconstruction confirms these findings (Figure 1(c1,c2)). At a lower magnification (×2000), a rough surface profile is seen, with macropores entering the interior of the particles and contributing to the overall internal pore structure. However, higher magnification at the pore wall (×5000) indicates a rather smooth surface.
The EDX analysis showed that C, O, and N are present in high percentages in SCGs, namely 62.6 at. %, 28.1 at. %, and 9.0 at. %, respectively. Besides those mentioned, Mg, P, Ca, K, and S are present in traces ( Figure 2, top row), constituting less than 0.3 at. % of the SCG sample. These elements are uniformly distributed in the sample, as indicated by the EDX maps ( Figure 2, bottom row). The EDX analysis agrees with the literature and elemental composition of SCGs, as the other research reported a similar elemental composition with a different order of trace elements from source to source [36,37]. Here, we found that the trace elements' composition (in at. %) follows the order K (0.09) > S = Mg (0.05) > P = Ca (0.02). We note that a previously low N content (<2%) was reported for SCGs [10], but using the elemental analysis, which also gives the H content, thus effectively reducing the contents of the other elements in the samples; meanwhile, here, H cannot be determined using EDX. The EDX analysis showed that C, O, and N are present in high percentages in SCGs, namely 62.6 at. %, 28.1 at. %, and 9.0 at. %, respectively. Besides those mentioned, Mg, P, Ca, K, and S are present in traces ( Figure 2, top row), constituting less than 0.3 at. % of the SCG sample. These elements are uniformly distributed in the sample, as indicated by the EDX maps ( Figure 2, bottom row). The EDX analysis agrees with the literature and elemental composition of SCGs, as the other research reported a similar elemental composition with a different order of trace elements from source to source [36,37]. Here, we found that the trace elements' composition (in at. %) follows the order K (0.09) > S = Mg (0.05) > P = Ca (0.02). We note that a previously low N content (<2%) was reported for SCGs [10], but using the elemental analysis, which also gives the H content, thus effectively reducing the contents of the other elements in the samples; meanwhile, here, H cannot be determined using EDX.

Surface Functional Groups of SCGs
The FTIR spectrum of the SCGs is shown in Figure 3a. The obtained FTIR spectrum correlates well with the literature and assumed compounds in SCGs [36][37][38]. The bands present in the FTIR spectrum of the SCGs are assigned to a specific vibration. The broad band on 3307 cm −1 was attributed to the stretching of the O-H group due to hydrogen bonding. Bands at 3010 cm −1 and 2924 cm −1 indicated the presence of C-H bonds where C is sp 2 and sp 3 hybridized, respectively. The sharp band on 1745 cm −1 was assigned to C=O stretch vibration, while the 1635 cm −1 band was attributed to C=N stretch vibration. The presence of a condensed aromatic system is assumed due to bands on 1521 cm −1 and 1441 cm −1 representing aromatic skeletal stretch and aromatic vibrations coupled with aromatic C-H in-plane vibrations, respectively. It is a clear indication of ligning presence, as found in Ref. [10]. An aliphatic C-H stretch vibration was assigned to the band on 1376 cm −1 . The

Surface Functional Groups of SCGs
The FTIR spectrum of the SCGs is shown in Figure 3a. The obtained FTIR spectrum correlates well with the literature and assumed compounds in SCGs [36][37][38]. The bands present in the FTIR spectrum of the SCGs are assigned to a specific vibration. The broad band on 3307 cm −1 was attributed to the stretching of the O-H group due to hydrogen bonding. Bands at 3010 cm −1 and 2924 cm −1 indicated the presence of C-H bonds where C is sp 2 and sp 3 hybridized, respectively. The sharp band on 1745 cm −1 was assigned to C=O stretch vibration, while the 1635 cm −1 band was attributed to C=N stretch vibration. The presence of a condensed aromatic system is assumed due to bands on 1521 cm −1 and 1441 cm −1 representing aromatic skeletal stretch and aromatic vibrations coupled with aromatic C-H in-plane vibrations, respectively. It is a clear indication of ligning presence, as found in Ref. [10]. An aliphatic C-H stretch vibration was assigned to the band on 1376 cm −1 . The deformation vibration of the C-N and C-O bands in secondary alcohols and C-N are assigned to bands on 1237 cm −1 and 1152 cm −1 . The band on 1028 cm −1 was  The TGA curve ( Figure 3c) shows that by increasing the temperature from 20 to 900 °C in the helium atmosphere, SCGs go through three significant weight losses. In the first stage, from 20 to 85 °C, SCGs lose about 15% of their weight. The second stage starts at 200 and ends at 500 °C. It is the stage that represents carbonization, and during this stage, SCGs lose 48% of their initial weight. The third stage, during which SCGs lose another 15% of their initial weight, starts at 670 °C and lasts up to 900 °C. At the end of the heating process, 12% of the initial SCGs' weight is left.
The TPD curves for the most common functional groups and molecules desorbed during the heating process of SCGs are shown in Figure 3b. The same trend that TGA showed can be seen on these curves. In the correlation of TGA and TPD, it can be concluded that, during the first stage (20-85 °C), SCGs release physisorbed water. During the phase of SCG carbonization (200-500 °C), the weight loss is primarily due to the release of chemisorbed water, CO, and CO2 molecules. Besides these, it can be noticed that all of the investigated functional groups and molecules are desorbed in this stage. During the third stage (670-900 °C), additional CO and CO2 molecules are released, leading to the final weight loss. The evolved functional groups agree well with the elemental content found by EDX and also with the FTIR spectra of SCGs. A high concentration of H is expected in the sample, in line with Ref. [10], originating from physisorbed and strongly bound water, OH and NH2 groups, and terminal C-H bonds. The TGA curve ( Figure 3c) shows that by increasing the temperature from 20 to 900 • C in the helium atmosphere, SCGs go through three significant weight losses. In the first stage, from 20 to 85 • C, SCGs lose about 15% of their weight. The second stage starts at 200 and ends at 500 • C. It is the stage that represents carbonization, and during this stage, SCGs lose 48% of their initial weight. The third stage, during which SCGs lose another 15% of their initial weight, starts at 670 • C and lasts up to 900 • C. At the end of the heating process, 12% of the initial SCGs' weight is left.
The TPD curves for the most common functional groups and molecules desorbed during the heating process of SCGs are shown in Figure 3b. The same trend that TGA showed can be seen on these curves. In the correlation of TGA and TPD, it can be concluded that, during the first stage (20-85 • C), SCGs release physisorbed water. During the phase of SCG carbonization (200-500 • C), the weight loss is primarily due to the release of chemisorbed water, CO, and CO 2 molecules. Besides these, it can be noticed that all of the investigated functional groups and molecules are desorbed in this stage. During the third stage (670-900 • C), additional CO and CO 2 molecules are released, leading to the final weight loss. The evolved functional groups agree well with the elemental content found by EDX and also with the FTIR spectra of SCGs. A high concentration of H is expected in the sample, in line with Ref. [10], originating from physisorbed and strongly bound water, OH and NH 2 groups, and terminal C-H bonds.

Kinetics of Adsorption
To investigate the kinetic parameters of the adsorption process, 1 mg mL −1 of SCG suspension was incubated with CHP and MLT in the concentration 5 × 10 −5 mol dm −3 for various time intervals (from 1 min to 24 h) at 25 • C. The concentration of adsorbed pesticides was calculated as a difference between their initial concentration (C 0 ) and the equilibrium concentration (C e ) measured with UPLC after removing the adsorbent. The kinetics of CHP and MLT adsorption onto SCGs were analyzed using the nonlinear pseudo-first, pseudo-second, Elovich, and intraparticle diffusion models. The corresponding equations are listed in Table 1. Table 1. Kinetic models' equations.

Kinetic Model Equation
Pseudo-first-order model In Table 1, parameter q t represents the amount of adsorbed OP at time t (mg g −1 ), while q e is the amount of adsorbed OP at equilibrium (mg g −1 ). Parameter k 1 stands for the adsorption rate of pseudo-first order (min −1 ), k 2 is the adsorption rate of pseudo-second order (g mg −1 min −1 ), α is the initial adsorption rate in the Elovich model (mg g −1 min −1 ), and β is desorption constant in the Elovich model (g mg −1 ). Parameter k id represents the adsorption rate constant of the intraparticle diffusion model (mg g −1 min −0.5 ), and C is connected to a boundary layer (mg g −1 ).
The nonlinear forms of the mentioned kinetic models and the intraparticle diffusion model fittings are shown in Figure 4. The kinetic parameters and corresponding R 2 , χ 2 , and Root Mean Square Error (RMSE) values are given in Table 2.

Kinetics of Adsorption
To investigate the kinetic parameters of the adsorption process, 1 mg mL −1 of SCG suspension was incubated with CHP and MLT in the concentration 5 × 10 −5 mol dm −3 for various time intervals (from 1 min to 24 h) at 25 °C. The concentration of adsorbed pesticides was calculated as a difference between their initial concentration (C0) and the equilibrium concentration (Ce) measured with UPLC after removing the adsorbent. The kinetics of CHP and MLT adsorption onto SCGs were analyzed using the nonlinear pseudofirst, pseudo-second, Elovich, and intraparticle diffusion models. The corresponding equations are listed in Table 1.

Kinetic Model
Equation Intraparticle diffusion model q t = k id t 0.5 + C In Table 1, parameter qt represents the amount of adsorbed OP at time t (mg g −1 ), while qe is the amount of adsorbed OP at equilibrium (mg g −1 ). Parameter k1 stands for the adsorption rate of pseudo-first order (min −1 ), k2 is the adsorption rate of pseudo-second order (g mg −1 min −1 ), α is the initial adsorption rate in the Elovich model (mg g −1 min −1 ), and β is desorption constant in the Elovich model (g mg −1 ). Parameter kid represents the adsorption rate constant of the intraparticle diffusion model (mg g −1 min −0.5 ), and C is connected to a boundary layer (mg g −1 ).
The nonlinear forms of the mentioned kinetic models and the intraparticle diffusion model fittings are shown in Figure 4. The kinetic parameters and corresponding R 2 , χ 2 , and Root Mean Square Error (RMSE) values are given in Table 2. By observing the results presented in Figure 4, it can be concluded that the adsorption equilibrium is reached after 10 min for MLT and 1440 min for CHP. In addition, it can be supposed that in the case of both MLT and CHP adsorption, the pseudo-first-order kinetic model fitted the experimental data better compared to the pseudo-second and Elovich   Figure 4, it can be concluded that the adsorption equilibrium is reached after 10 min for MLT and 1440 min for CHP. In addition, it can be supposed that in the case of both MLT and CHP adsorption, the pseudo-first-order kinetic model fitted the experimental data better compared to the pseudo-second and Elovich models, according to the R 2 and χ 2 values. The calculated adsorption capacities of SCGs under experimental conditions are 3.31 mg g −1 for MLT and 2.26 mg g −1 for CHP. A higher k 1 value for MLT adsorption indicates a higher adsorption rate than that of CHP.
The value α, representing the initial adsorption rate, and the value β, the desorption constant, obtained from the Elovich model constant also indicate that the adsorption of CHP is noticeably slower than the adsorption of MLT on SCGs, as the α value is lower than the β value in the case of CHP adsorption.
By plotting q t against t 0.5 , three linear stages for MLT and two linear stages for CHP could be clearly observed in the intraparticle kinetic model (Figure 3c). The breakpoints for MLT adsorption onto SCGs are 10 min and 60 min, and for CHP, 10 min. The first rapid stage is attributed to the diffusion of MLT and CHP molecules through the solution to the external surface of the SCGs [39]. The second stage represents intraparticle diffusion, in which MLT and CHP molecules diffuse into the pores of the SCGs. The third stage seen for MLT is attributed to the final equilibrium stage. The K id values decrease by one order of magnitude after the breakpoints, indicating slower adsorption in every subsequent stage. The C value, representing the boundary layer, increases, suggesting that the boundary layer has a high significance in the adsorption of MLT and CHP onto SCGs. It can be concluded that the three mentioned processes control the rate of molecule adsorption for MLT and two for CHP, but for each time range during the adsorption process, one process at a time determines the adsorption kinetics. We believe that differences in the molecular structure of MLT and CHP, the first being aliphatic, while the other one possessing aromatic moiety (Figure 4, insets), renders drastically different adsorption kinetics. As CHP is a less polar molecule, also reflected in much lower solubility compared to MLT, this property likely reduces its adsorption rate on the solvated surface of SCGs. It also slows down its diffusion into pores, as stage two in the intraparticle diffusion model for CHP extends over 1000 min without reaching the final equilibration stage.

Adsorption Isotherms
To further study the adsorption process, 1 mg mL −1 of SCGs was incubated with MLT and CHP in the concentration range from 5 × 10 −6 to 1 × 10 −4 mol dm −3 for 24 h at 25, 30, and 35 • C. The obtained data were fitted with several nonlinear isotherm models (Freundlich, Langmuir, Temkin, and Dubinin-Radushkevich (DR)). The corresponding equations for used isotherms are given in Table 3. Table 3. Adsorption isotherms' equations.

Model Nonlinear Form
Freundlich isotherm q t = K F C 1 n e Langmuir isotherm q t = q max K L C e 1+K L C e Temkin isotherm q t = RT b T lnK T C e Dubinin-Radushkevich isotherm q t = q DR e −K DR ε 2 In Table 3, the q t parameter represents the amount of adsorbed OP for a given C e (mg g −1 ), C e is the equilibrium concentration (mg dm −3 ), and K F is the Freundlich isotherm constant related to the adsorption capacity ((mg g −1 )(dm 3 g −1 ) 1/n ). Parameter n stands for the so-called adsorption intensity, q max represents the maximum amount of adsorbed OP (mg g −1 ), and K L is the Langmuir isotherm constant (dm 3 mg −1 ). R is the universal gas constant (8314 J K −1 mol −1 ), T is the temperature (K), and b T represents the Temkin isotherm constant (J g mol −1 mg −1 ), and K T is the Temkin isotherm equilibrium binding constant (dm 3 mg −1 ). At the same time, q DR is the theoretical isotherm saturation capacity (mg g −1 ), while K DR (mol 2 J −2 ) is the Dubinin-Radushkevich isotherm constant. The Polanyi potential, ε (J mol −1 ), is calculated from C e , while K DR is related to the adsorption mean free energy (E), as E = 1/(2K DR ) 1/2 .
Nonlinear forms fitting of all investigated isotherms are shown in Figure 5. The parameters obtained from the different models provide essential information on the adsorption mechanisms and the surface properties and affinities of SCGs towards MLT and CHP. Table 4 summarizes the parameters obtained from the nonlinear plots of each isotherm. The fitting of the models to the experimental data was evaluated by χ 2 and R 2 . The R 2 values are the highest (R 2 > 0.9), while the χ 2 and RMSE values are the lowest for the Langmuir model, suggesting that that model could best describe the experimental data. Still, most of the used isotherm models fit with most experimental data, so the calculated parameters are relevant for discussion. The n-value of the Freundlich isotherm model for MLT adsorption is higher than 1, indicating that the adsorption is a favorable process. The decrease of the n-value with the increase in temperature indicates that the affinity of the adsorbent for MLT declines with the rising temperature. As mentioned, the Langmuir model showed the best agreement with the experimental results, indicating that the monolayer of MLT is adsorbed on an energetically homogenous surface without the interaction between adsorbed molecules. It indicates that all active sites are energetically equivalent and that equilibrium is reached when a monolayer is formed. Even though the temperature increase had almost no impact on the maximum amount (qmax = 7.16 mg g −1 ) of MLT that could be adsorbed by SCGs, with the decrease of the KL value, it can be concluded that the interaction between MLT and SCGs was weakened. It is also confirmed by the decrease of the E value obtained from the Dubinin-Radushkevich isotherm model. As E << 8000 J mol −1 , the proposed mechanism of MLT adsorption on SCGs is likely to be physisorption [40]. The Temkin isotherm model showed that the bT values, which are related to adsorption heat, increase with the temperature. It could indicate that the molecules of MLT are more readily adsorbed onto the SCGs' surface, contrasting with the previous parameters obtained from the rest of the tested isotherms. However, since the χ 2 and RMSE values are higher for the Temkin isotherm than other investigated models, it is unsurprising that the findings provided different indications. The n-value of the Freundlich isotherm model for MLT adsorption is higher than 1, indicating that the adsorption is a favorable process. The decrease of the n-value with the increase in temperature indicates that the affinity of the adsorbent for MLT declines with the rising temperature. As mentioned, the Langmuir model showed the best agreement with the experimental results, indicating that the monolayer of MLT is adsorbed on an energetically homogenous surface without the interaction between adsorbed molecules. It indicates that all active sites are energetically equivalent and that equilibrium is reached when a monolayer is formed. Even though the temperature increase had almost no impact on the maximum amount (q max = 7.16 mg g −1 ) of MLT that could be adsorbed by SCGs, with the decrease of the K L value, it can be concluded that the interaction between MLT and SCGs was weakened. It is also confirmed by the decrease of the E value obtained from the Dubinin-Radushkevich isotherm model. As E << 8000 J mol −1 , the proposed mechanism of MLT adsorption on SCGs is likely to be physisorption [40]. The Temkin isotherm model showed that the b T values, which are related to adsorption heat, increase with the temperature. It could indicate that the molecules of MLT are more readily adsorbed onto the SCGs' surface, contrasting with the previous parameters obtained from the rest of the tested isotherms. However, since the χ 2 and RMSE values are higher for the Temkin isotherm than other investigated models, it is unsurprising that the findings provided different indications. In the case of CHP adsorption, the n-values obtained from Freundlich isotherm models are also above 1, again indicating the favorability of the adsorption process. In addition, the n-value decreases with the temperature increase. The affinity of the adsorbent for CHP is declining with the rising temperature. Still, the n-values should be taken cautiously regarding the not-high R 2 values for all temperatures (Table 4). Similar to the results for MLT adsorption, the Langmuir isotherm model provides the best fit for adsorption at all temperatures. In theory, that indicates that CHP is adsorbed on an energetically homogenous surface without the interaction between adsorbed molecules in the form of the monolayer. Adsorption of CHP onto SCGs at 25 • C provides the highest K L value (K L = 9.23 dm 3 mg −1 ). It is in agreement with the n-values from the Freundlich isotherm. The obtained qmax values rise with the increasing temperature, and the maximum adsorption capacity is calculated to be 7.00 mg g −1 . The energy of adsorption obtained from the DR isotherm model indicates that the adsorption mechanism is physisorption. Still, at the temperature of 25 • C, there is a strong indication that the chemisorption also occurs. The Temkin isotherm model showed that the b T values decreased with the temperature increase, meaning that the readiness of the CHP molecule for adsorption on SCGs was decreasing, which agrees with the other calculated parameters.

Thermodynamic Analysis
The thermodynamic parameters of adsorption provide information about the interactions between the adsorbate molecules and the surface. The most important thermodynamic adsorption parameters are the enthalpy (∆H), entropy (∆S), and Gibbs free energy (∆G) of the adsorption process. These parameters give an indication of the heat, randomness changes, and spontaneity of the process. However, we note that the exact values of the thermodynamic parameters strongly depend on the choice of the standard state, as discussed in Ref. [41]. Thus, we adopted the recommendations for the standard state by Chen et al. [41].
The mentioned thermodynamic parameters are linked through the Gibbs-Helmholtz equation: To determine the ∆H 0 and ∆S 0 values to calculate ∆G 0 , it is necessary to plot the Van't Hoff equation, assuming that ∆G 0 = −RTlnK 0 dist : where the standard distribution coefficient is defined by the following: The values of C 0 and q 0 define the standard state for OP in solution (1 mol dm −3 ) and in the adsorbed state (1 mol kg −1 ), rendering the standard distribution coefficient to dimensionless quantity. Figure 6a shows graphic representations of the amount of MLT and CHP dependence on temperature and Van't Hoff plot (Figure 6b). The extracted thermodynamic parameters for the adsorption of MLT and CHP onto SCGs are summarized in Table 5. (ΔG) of the adsorption process. These parameters give an indication of the heat, randomness changes, and spontaneity of the process. However, we note that the exact values of the thermodynamic parameters strongly depend on the choice of the standard state, as discussed in Ref. [41]. Thus, we adopted the recommendations for the standard state by Chen et al. [41]. The mentioned thermodynamic parameters are linked through the Gibbs-Helmholtz equation: To determine the ΔH 0 and ΔS 0 values to calculate ΔG 0 , it is necessary to plot the Van't Hoff equation, assuming that ΔG 0 = −RTlnK dist 0 : where the standard distribution coefficient is defined by the following: The values of C 0 and q 0 define the standard state for OP in solution (1 mol dm −3 ) and in the adsorbed state (1 mol kg −1 ), rendering the standard distribution coefficient to dimensionless quantity. Figure 6a shows graphic representations of the amount of MLT and CHP dependence on temperature and Van't Hoff plot (Figure 6b). The extracted thermodynamic parameters for the adsorption of MLT and CHP onto SCGs are summarized in Table 5.  Figure 6a shows that increasing the temperature decreases the distribution coefficient of MLT, while the distribution coefficient of CHP increases. Temperature dependence is linear in both cases, but the effect of temperature is more pronounced for CHP. The obtained thermodynamic parameters (Table 5) explain this behavior and indicate that the adsorption of MLT is an exothermic process and that the adsorption is followed by a small negative change of entropy. Thus, for MLT, the adsorption on SCGs is an enthalpy-driven   Figure 6a shows that increasing the temperature decreases the distribution coefficient of MLT, while the distribution coefficient of CHP increases. Temperature dependence is linear in both cases, but the effect of temperature is more pronounced for CHP. The obtained thermodynamic parameters (Table 5) explain this behavior and indicate that the adsorption of MLT is an exothermic process and that the adsorption is followed by a small negative change of entropy. Thus, for MLT, the adsorption on SCGs is an enthalpydriven process, compensating for the reduction in entropy, giving ∆G 0 negative, and increasing with temperature. On the other hand, the case of CHP adsorption onto SCGs is the opposite, as adsorption is an endothermic process, but the adsorption process is followed by a large positive entropy. Hence, CHP adsorption onto SCGs is an entropydriven process. Again, the Gibbs free energy change is negative (slightly more negative than for MLT) and decreases with temperature, suggesting spontaneous adsorption, in line with experimental findings. Considering that a large fraction of the SCGs is constituted of lignin [10], we believe that the differences in thermodynamics of MLT and CHP adsorption can be explained by different types of interactions between SCGs and the studied OPs, keeping in mind that the adsorption process was classified as physisorption. We suggest that, in the case of MLT, dipole and dispersion interactions (Figure 7a) dominate the adsorption, thus becoming less prominent as the temperature increases. Such interactions can be placed between MLT molecules and different functional groups observed by FTIR (Figure 3a) and TPD. In the case of CHP, relatively strong π-π stacking between the aromatic ring of CHP and phenolic moieties of lignin can also take place (Figure 7b). However, in this case, desolvation of the surface has to take place, resulting in a positive change of the enthalpy and a large increase in entropy, resulting in negative Gibbs free energy change. As mentioned in Section 3.2.1., CHP and MLT show drastically different kinetic behavior, which was hypothesized to originate from different structures of these two pesticides. The argument mentioned above and the results regarding thermodynamic functions of the adsorption process support previous assumptions. In the case of CHP, it is proposed that desolvation of the SCG surface slows the adsorption down when compared to the MLT adsorption. As the dipole interactions are electrostatic, thus decaying with the distance but generally reaching longer ranges then π-π stacking (operative at distances below 4 Å [42]), it is assumed that MLT adsorption can proceed without significant slowdown due to desolvation. At the same time, the molecular forces capture MLT molecules at more considerable distances compared to CHP. two CHP molecules via π-π stacking interactions. Structures were optimized by molecular mechanics in Avogadro [44] with an MMFF94 force field [45].

Application of SCGs in Food Processing-Plant Extracts Testing and Eco-Neurotoxicology Assessment
To assess the potential for SCGs' application in food processing, the material was tested as an adsorbent in plant extracts. As mentioned, MLT is used to treat lemon crops, and CHP is used in mint crops. Therefore, diluted lemon juice for MLT and mint extract for CHP were prepared as described in Section 2.4. The samples were spiked with respective OP in the concentration of 5 × 10 −5 mol dm −3 . For comparison, the adsorptions of MLT and CHP from deionized water and 50% ethanol, respectively, were investigated under identical conditions. The results are shown in Table 6.
Our primary goal was to reduce the toxicity of plant extract samples and confirm that no respective oxo-analogs are formed during adsorption. To estimate the eco-neurotoxicity of the samples before and after the adsorption of OPs on SCGs, the AChE inhibition was measured, as described in Section 2.5. AChE inhibition is an important benchmark for eco-neurotoxicity, as it indicates a general impact of contaminants on the species containing AChE in their nervous systems [46]. The results are also summarized in Table 6. Table 6. OP concentrations and AChE inhibition in the tested samples (MLT in deionized water and lemon juice and CHP in 50% ethanol and mint ethanol extract) before and after the adsorption. Results were obtained when 1 mg mL −1 of adsorbent (SCG) was used at 25 °C. From the data presented in Table 6, it is obvious that the matrix effects in plant extracts tested in this work have a minor impact on the OPs' uptake. The amount of the adsorbed pesticides did not significantly change in the plant-extract analysis compared to Structures were optimized by molecular mechanics in Avogadro [44] with an MMFF94 force field [45].

Application of SCGs in Food Processing-Plant Extracts Testing and Eco-Neurotoxicology Assessment
To assess the potential for SCGs' application in food processing, the material was tested as an adsorbent in plant extracts. As mentioned, MLT is used to treat lemon crops, and CHP is used in mint crops. Therefore, diluted lemon juice for MLT and mint extract for CHP were prepared as described in Section 2.4. The samples were spiked with respective OP in the concentration of 5 × 10 −5 mol dm −3 . For comparison, the adsorptions of MLT and CHP from deionized water and 50% ethanol, respectively, were investigated under identical conditions. The results are shown in Table 6. Table 6. OP concentrations and AChE inhibition in the tested samples (MLT in deionized water and lemon juice and CHP in 50% ethanol and mint ethanol extract) before and after the adsorption. Results were obtained when 1 mg mL −1 of adsorbent (SCG) was used at 25 • C. Our primary goal was to reduce the toxicity of plant extract samples and confirm that no respective oxo-analogs are formed during adsorption. To estimate the eco-neurotoxicity of the samples before and after the adsorption of OPs on SCGs, the AChE inhibition was measured, as described in Section 2.5. AChE inhibition is an important benchmark for eco-neurotoxicity, as it indicates a general impact of contaminants on the species containing AChE in their nervous systems [46]. The results are also summarized in Table 6.

Sample
From the data presented in Table 6, it is obvious that the matrix effects in plant extracts tested in this work have a minor impact on the OPs' uptake. The amount of the adsorbed pesticides did not significantly change in the plant-extract analysis compared to deionized water and ethanol, indicating that SCGs as the adsorbent maintain a good performance and can be successfully used in plant extracts' treatment. Moreover, the AChE inhibition due to the contact with samples was reduced after the adsorption in all cases, indicating that no more toxic products, such as OPs' oxo-forms, were formed during this process. While the reduction of toxicity is not as prominent as in the case of some much more potent adsorbents, such as carbon materials [47,48], we note that SCGs' performance falls well within the range of previously reported biowaste adsorbents and surpasses many of them in terms of MLT and CHP adsorption capacity, as is discussed further. Moreover, an important advantage is that such an adsorbent is truly green and sustainable, as it does not require a carbonization step during which large amounts of CO 2 are released. In fact, with adsorption capacities in the order of 10 0 mg g −1 , SCGs perform quite similarly to carbons derived from biomass if the capacity is calculated versus raw material (precursor). Namely, during the carbonization and activation steps, a large amount of the precursor is lost, giving a yield in the order of 10 0 -10 1 % [49]. Thus, carbon materials with typical capacities for OPs' removal around 100 mg g −1 under high and realistic adsorbent loadings [50] have adsorption capacities in the order of 10 0 mg g −1 when referring to the raw precursor mass. Another critical point is the actual affinity of adsorbents towards contaminants. Here, the presented SCGs have a very low SSA (the direct measurement was not possible), but the adsorption capacities are in the order of 10 0 mg g −1 . Carbon-based adsorbents, on the other hand, have SSAs in the order of hundreds of m 2 g −1 . Thus, the affinity of the SCG surface towards MLT and CHP is quite similar to that of carbon materials, resulting from specific chemical composition and surface chemistry ( Figure 7). Hence, as the adsorbent surface, not the adsorbent mass, is responsible for the adsorption process, the guidelines for improving SCGs' (mass normalize) adsorption capacities are straightforward, increasing SSA while maintaining surface chemistry unchanged. In the next section, in order to put the obtained results in a proper context, we compare the performance of the here-described SCGs with previously reported biowaste-based adsorbents for MLT and CHP.

Biowaste as Adsorbent for Different Pollutants-Evaluation of SCGs' Performance Compared to the Data Available in the Literature
SCGs have been studied for the adsorption of various water contaminants. The removal of Cd using SCGs was investigated by Patterer et al. [51] and Kim and Kim [52]. Patterer et al. concluded that the adsorption of Cd is an exothermic and spontaneous process and that the adsorption capacity is 4.484 mg g −1 [51]. In their research, Kim and Kim reported that the adsorption of Cd is irrelevant to the pH values of the solution in the range of 4-8, while the adsorption capacity of the SCGs was 19.32 mg g −1 [52]. Besides heavy metals, SCGs are shown to be excellent adsorbents of dyes. In their work, Block et al. found that non-treated SCGs are a promising adsorbent for methylene blue, removing 80% of the dye in 3 h of contact but adsorbing only 20% of methyl orange, making it a partly selective adsorbent [53].
Many sustainable materials can be effectively used as adsorbents for CHP. For example, Memon et al. successfully used walnut shell powder to remediate chlorpyrifos. The adsorption efficiency increased with the increase in temperature, similar to the herereported behavior of SCGs. The used material had an adsorption capacity of 0.99 mmol g −1 (347 mg g −1 ) [54]. Moreover, Rojas et al. investigated the adsorption of chlorpyrifos on sunflower seed shells, rice husk, composted sewage sludge, and soil and concluded that their adsorption capacities are 0.035 mg g −1 , 0.014 mg g −1 , 0.010 mg g −1 , and 0.02 mg g −1 respectively [55]; thus, two orders of magnitude below here reported capacities of SCGs towards CHP.
Apart from CHP, research has been conducted on MLT remediation, but to a much lesser extent. The biosorption of MLT from aqueous solutions was investigated using two types of herbal leaf powder by Yadamari et al. The adsorption capacity for Achyranthes aspera was 3.401 mg g −1 , and for Phyllanthus niruri, it was 2.664 mg g −1 [56]. Veličković et al. reported the adsorption capacity of malathion on pulverized river shellfish shells of 46.462 mg g −1 . The adsorption of malathion in the mentioned investigation was assessed as exothermic, similar to the results presented in this work [57]. A comparative analysis of the maximum adsorption capacity obtained in this work and the previously mentioned research available in the literature is given in Table 7. From the presented data, it can be seen that SCGs are a reasonably good adsorbent for MLT and CHP remediation. Moreover, the obtained adsorption capacities are comparable with the previously reported MLT and CHP adsorption capacities on different types of biowastes and raw materials. Considering the discussion provided in the previous section regarding the performance of carbon-based adsorbents, losses of precursor mass, and significant CO 2 evolution during the carbonization process, we believe that SCGs have the potential for practical application in remediation processes, as they effectively stand side by side with apparently much more efficient adsorbents. However, we acknowledge that a detailed (socio-)techno-economical assessment is needed to correctly judge the possible use of SCGs in environmental remediation processes.

Conclusions
In this work, MLT and CHP adsorption on SCGs from aqueous solutions was analyzed. The adsorption kinetics was best described by the pseudo-first-order kinetic model for both MLT and CHP adsorption, while MLT adsorption was found to be significantly faster compared to CHP. The Langmuir isotherm model best describes the adsorption process for both pesticides, giving the maximal adsorption capacity of 7.16 mg g −1 for MLT and 7.00 mg g −1 for CHP; while based on the Dubinin-Radushkevich isotherm model, the adsorption process can be safely classified as physisorption. MLT adsorption on SCGs was found to be exothermic, whereas CHP adsorption is endothermic and entropy driven. The adsorption efficiency of MLT and CHP using SCGs in a complicated matrix of fruit extracts remained constant. The eco-neurotoxicity results showed that no more toxic products were formed during adsorption, indicating that SCGs are a safe-to-use adsorbent for pesticide removal. Their benefits as an adsorbent should be considered in line with the fact that their adsorption performance is similar to carbon materials if referred to the precursor mass used for carbon production while completely removing complicated steps involving carbonization and activation, resulting in significant CO 2 emissions. Considering that SCGs are globally available, their use in the remediation process also alleviates the problem of waste management, being particularly beneficial for less-developed countries, which, as a rule, suffer from extensive pesticide usage.

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