Investigating the Adsorption Kinetics of Dimethoate, Malathion and Chlorpyrifos on Cellulose-Derived Activated Carbons: Understanding the Inﬂuence of Physicochemical Properties

: In light of the escalating environmental concerns regarding pesticide accumulation, it is imperative to devise efﬁcient strategies for their removal. Among the various options, activated carbons have emerged as promising candidates for adsorptive pesticide removal due to their many advantages, such as large surface area, well-developed porosity, and cost-effectiveness. However, the intricate relationship between the properties of these materials and their performance in pesticide adsorption remains largely unexplored. This study primarily focuses on examining the adsorption kinetics of three organophosphate pesticides: dimethoate, malathion (aliphatic), and chlorpyrifos (aromatic), using a range of cellulose-based activated carbon ﬁbers with diverse speciﬁc surface areas, pore size distributions, and elemental compositions. By employing sophisticated data analysis tools, principal component analysis, and semi-empirical quantum chemical calculations, this study uncovers the importance of these distinct properties in efﬁciently removing structurally diverse pesticides. The results of the adsorption experiments suggested that these processes can be described using a pseudo-second-order kinetic model, which is conﬁrmed via multiple linear regression. The obtained data suggest that the most effective carbon material for pesticide removal should have a pore diameter of approximately 4 nm, low oxygen content, a unimodal pore size distribution, and a high presence of sp 2 domains. The insights from this research have the potential to guide the development of improved adsorbents and facilitate the rational selection of adsorbents tailored to speciﬁc pollutants based on their physicochemical properties and the pollutants’ chemical structure. By shedding light on the vital connection between adsorbent properties and performance, our ﬁndings signiﬁcantly advance sustainable and effective pesticide removal, thereby fostering a cleaner and healthier environment.


Introduction
Organophosphates (OPs) represent a widely used class of pesticides, posing a significant global threat to both the environment and human health [1][2][3].Their toxicity towards animals stems from their inhibition of acetylcholinesterase (AChE), a critical enzyme in The adsorption process is also influenced by solution conditions such as pH, ionic strength, and temperature [12][13][14].Notably, the pH of the solution plays a critical role in determining the surface charge of carbon materials and the dissociation or protonation of electrolytes [13].When the pH of the solution surpasses the carbon material's zero charge pH, the surface acquires a negative charge and becomes capable of attracting cations from the solution.Conversely, when the pH falls below the carbon material's zero charge pH, the surface turns positive, attracting anions [13,15,20,21].A higher pH pzc value results in a higher adsorption rate in basic solutions.Similarly, in acidic solutions, a lower pH pzc results in a higher adsorption rate [15,21].These findings have been demonstrated for reactive dyes, including methylene blue [21].Additionally, the solution's pH controls the dissociation or ionization of electrolytes through their respective pKa values, where acidic electrolytes dissociate at a pH > pKa [13].
In the aqueous phase, a variety of interactions, including van der Waals, induceddipole, dipole-dipole, and donor-acceptor hydrogen bonding forces, play crucial roles in binding and accumulating chemical compounds on carbonaceous adsorbents [14].Among these interactions, hydrogen and π-π bonds, covalent and electrostatic interactions, as well as the hydrophobic effect, contribute significantly to the adsorption process [14,22,23].
Ionic strength is another pivotal factor influencing electrostatic interactions.By increasing the solution's ionic strength, attractive and repulsive interactions can be reduced due to the screening effect of the surface charge produced by the added salt.As a result, increasing the ionic strength leads to enhanced adsorption when electrostatic interactions between the surface and adsorbate are repulsive or when the surface concentration is sufficiently high.Conversely, increasing ionic strength will reduce adsorption when electrostatic interactions are attractive or when the surface concentration is relatively low.Concerning the impact of temperature on the adsorption process, a decrease in adsorption temperature is generally expected to lead to an increased intake of organic molecules since adsorption is a spontaneous process.However, some instances demonstrate an increase in the adsorbed amount with temperature [13], which is not typical for adsorption.Primarily, adsorption is an exothermic process, meaning that the adsorption capacity decreases with the increase in temperature, but in some cases, high temperature can be more beneficial for the adsorption of molecules on the surface of the adsorbent.This can be explained by the fact that when the temperature increases, an enhancement in the diffusion rate of the molecules adsorbed on the adsorbent surface through the external boundary layer and into the internal pores can occur, leading to a higher adsorption rate.Also, temperature can affect the adsorbent and its adsorption sites and activity, e.g., an increase in the available pores for adsorption might occur due to deformation of the adsorbent with increasing temperature.
While numerous materials have been investigated as adsorbents, a comprehensive understanding of the critical material properties that define adsorption performance remains elusive.This complexity arises from the interplay of multiple parameters influencing the adsorption process.In this paper, we present a systematic analysis using a well-characterized series of activated carbon fibers to explore the impact of their properties on the kinetics of chlorpyrifos, dimethoate, and malathion adsorption from aqueous solutions.Specifically, we correlate adsorbent properties such as pore volume integrated up to a given pore diameter, total pore volume, surface area, and elemental content with the kinetics of pesticide adsorption onto the investigated materials.Multiple linear regression analysis is employed to accomplish this task, aided by theoretical calculations.By combining empirical data with theoretical insights, we propose an optimal adsorbent with an ideal combination of pore diameter, sp 2 domains, and heteroatom concentration for efficient chlorpyrifos, dimethoate, and malathion adsorption.This innovative approach provides valuable information for the rational design of materials with immense potential for applications in environmental protection.Through this study, we aim to contribute significantly to the understanding of the critical factors influencing adsorption performance, paving the way for developing highly effective and tailored materials for tackling pesticide contamination in aqueous environments.

Activated Carbon Fiber Preparation
The preparation of activated carbon fibers (ACF) is described in [24] and reproduced here for completeness.Viscose fibers 1.7 dtex, 38 mm (Lenzing AG, Lenzing, Austria), were impregnated with diammonium hydrogen phosphate (DAHP) solutions of different concentrations (0.0-75.7 mmol dm −3 , corresponding to 0.00-10.00wt.% DAHP) for 15 min.After a 24 h drying period, carbonization was conducted at 850 • C for 20 min under a nitrogen atmosphere, with a heating rate of 1.0 • C min −1 .The carbonized fibers were then activated at 870 • C for 165 min with CO 2 flow (80 dm 3 h −1 ).No additional washing was needed due to pure precursor fibers and non-altering impregnation agents [25][26][27].The samples were labeled as DAHP-X, where X represents the concentration of DAHP used in the impregnation step.

Characterization of Materials
The ACF samples' morphology and elemental composition were investigated using a scanning electron microscope (PhenomProX, Thermo Fisher Scientific, Waltham, MA, USA) equipped with energy-dispersive X-ray spectroscopy (EDX).The textural properties and specific surface area of materials were studied via N 2 isothermal adsorption (−196.15• C).For this purpose, a gas sorption system was used (AutosorbiQ, Quantachrome Instruments, Graz, Austria).The specific surface area and derived pore size distribution (PSD) were calculated using the Brunauer-Emmett-Teller (BET) method and the non-local density functional theory (NLDFT), respectively [28][29][30].

Pesticide Adsorption Measurements
Batch adsorption experiments were performed by dispersing the materials in double distilled water and adding the required pesticide stock solution (Pestanal, Sigma Aldrich, Søborg, Denmark) to achieve the desired concentrations of adsorbent and pesticide.The mixtures were then placed on an Orbital Shaker-Incubator at 25 • C for specific periods of time.After centrifugation and filtration, ultra performance liquid chromatography (UPLC) was used to determine the concentration of organophosphates (C eq ).Control experiments were conducted without materials to confirm no pesticide degradation.The pH of the dispersions remained consistent, indicating that differences in pesticide removal kinetics were solely due to adsorption on the materials.UPLC measurements were performed using an ACQUITY UPLC system (Waters, Milford, MA, USA) with a tunable UV detector and an ACQUITY UPLC™ BEH C18 column (1.7 µm, 100 mm × 2.1 mm, Waters, Milford, MA, USA) under isocratic conditions with acetonitrile and water as the mobile phase (see ref. [29] for details).

Principal Component Analysis
The principal component analysis (PCA) was conducted using the built-in functions provided by Scikit-learn.To ensure effective statistical analysis, the input variables were scaled using the StandardScaler function, accommodating various levels and ranges of the considered data.

Semi-Empirical Quantum Chemical Calculations
Semi-empirical calculations were performed using the MOPAC2016 code [31] with the PM7 method [32].A model of (7,7) SWCNT was constructed with 112 atoms to obtain the deformation energies of chlorpyrifos.Then, the chlorpyrifos molecule was placed into the tube, and full structural relaxation was carried out.After the relaxation, the system was split into nanotube and chlorpyrifos parts, and their molecular energies were calculated in a self-consistent cycle.Deformation energies were obtained by comparing the energy of fully relaxed chlorpyrifos/nanotube with the energy of chlorpyrifos/nanotubes with the mathematical result corresponding to the chlorpyrifos inserted into the nanotube.

Materials Properties
The materials studied in this work have been partially characterized and investigated, with respect to the thermodynamic properties' influence on adsorption capacities, before [24,29].One of the main conclusions was that the specific surface area does not play a dominant role in determining the adsorption capacity but rather pore volume and carbon and oxygen content do.However, the studied materials are micro-and mesoporous with different pore size distributions (Figure 1).Thus, it is still unclear which portion of the total pore volume is responsible for the overall uptake of pollutants, and particularly the kinetics of this process.This work aims to correlate the kinetics of pesticide adsorption to the adsorbents' physicochemical properties.Detailed analysis of the pore size distribution and elemental content of the materials are provided in Table 1.From the data presented, it is clear that the chemical composition of the adsorbents was influenced by the loading of DAHP during the impregnation process.Utilizing EDX, we observed a consistent rise in the phosphorus (P) content in the materials, while the levels of carbon (C) and oxygen (O) exhibited fluctuations but no obvious trend.Notably, the initial phosphorus-containing sample in the series exhibited a significant decrease in specific surface area (S tot ) when compared to the material produced without DAHP impregnation.Furthermore, as the concentration of DAHP increased, both the S tot and total volume (V tot ) increased.Additionally, the increase in phosphorus content led to a gradual increase in pore sizes.
the deformation energies of chlorpyrifos.Then, the chlorpyrifos molecule was placed into the tube, and full structural relaxation was carried out.After the relaxation, the system was split into nanotube and chlorpyrifos parts, and their molecular energies were calculated in a self-consistent cycle.Deformation energies were obtained by comparing the energy of fully relaxed chlorpyrifos/nanotube with the energy of chlorpyrifos/nanotubes with the mathematical result corresponding to the chlorpyrifos inserted into the nanotube.

Materials Properties
The materials studied in this work have been partially characterized and investigated, with respect to the thermodynamic properties' influence on adsorption capacities, before [24,29].One of the main conclusions was that the specific surface area does not play a dominant role in determining the adsorption capacity but rather pore volume and carbon and oxygen content do.However, the studied materials are micro-and mesoporous with different pore size distributions (Figure 1).Thus, it is still unclear which portion of the total pore volume is responsible for the overall uptake of pollutants, and particularly the kinetics of this process.This work aims to correlate the kinetics of pesticide adsorption to the adsorbents' physicochemical properties.Detailed analysis of the pore size distribution and elemental content of the materials are provided in Table 1.From the data presented, it is clear that the chemical composition of the adsorbents was influenced by the loading of DAHP during the impregnation process.Utilizing EDX, we observed a consistent rise in the phosphorus (P) content in the materials, while the levels of carbon (C) and oxygen (O) exhibited fluctuations but no obvious trend.Notably, the initial phosphorus-containing sample in the series exhibited a significant decrease in specific surface area (Stot) when compared to the material produced without DAHP impregnation.Furthermore, as the concentration of DAHP increased, both the Stot and total volume (Vtot) increased.Additionally, the increase in phosphorus content led to a gradual increase in pore sizes.Integral (Figure 1a) and differential pore volume (Figure 1b) distributions reveal that pore volumes increase with the amount of DAHP used, while pore sizes also increase.Table 1 shows the same trend.For example, the volume of pores with a diameter below 1 Integral (Figure 1a) and differential pore volume (Figure 1b) distributions reveal that pore volumes increase with the amount of DAHP used, while pore sizes also increase.Table 1 shows the same trend.For example, the volume of pores with a diameter below 1 nm decreases along the series, while the total pore volume increases.The exception is the first sample in the series, DAHP-0, which has a higher V tot than the samples DAHP-0.25 and DAHP-0.50, and also a higher specific surface area.We note that in these experiments, the ultramicropores were not able to be measured, but, based on the results presented later on, it is not likely that this fraction of pores can significantly contribute to the adsorption process.
As explained previously [24,29], these materials are particularly suitable for testing different hypotheses about the adsorption process.All these materials have been derived from the same precursor and have identical morphology (Figure 2).During the carbonization and activation steps, all samples kept their initially identical morphology, with the surface texture of the precursor being retained and the fiber diameter showing a similar C 2023, 9, 103 6 of 15 shrinkage of approx.25% for all samples.For this reason, the effects of the morphology on the rate of the adsorption processes can be safely excluded.
Table 1.Adsorbents' properties-pore volume integrated up to a given pore diameter, total pore volume (integrated up to pore diameter of 35 nm), surface area, and elemental content, determined using energy dispersive X-ray spectroscopy.The notation of samples is DAHP-X, where X represents the concentration of DAHP in wt.% used for the impregnation step.As explained previously [24,29], these materials are particularly suitable for testing different hypotheses about the adsorption process.All these materials have been derived from the same precursor and have identical morphology (Figure 2).During the carbonization and activation steps, all samples kept their initially identical morphology, with the surface texture of the precursor being retained and the fiber diameter showing a similar shrinkage of approx.25% for all samples.For this reason, the effects of the morphology on the rate of the adsorption processes can be safely excluded.Moreover, all the samples have undergone the same milling procedure, so the effects particle sizes on the adsorption performance can also be excluded.Additionally, due to the same carbonization and activation temperatures, the Raman spectra of all the samples were identical [24].Thus, the effects of the disorder degree can also be excluded.That leaves pore size distribution and elemental content (Table 1), effectively tuned by the amount of impregnating agent, DAHP, as the decisive properties that affect the adsorption kinetics.Due to the careful synthesis design, the parametric space in which the adsorption process can be investigated is significantly reduced, allowing for the derivation Moreover, all the samples have undergone the same milling procedure, so the effects particle sizes on the adsorption performance can also be excluded.Additionally, due to the same carbonization and activation temperatures, the Raman spectra of all the samples were identical [24].Thus, the effects of the disorder degree can also be excluded.That leaves pore size distribution and elemental content (Table 1), effectively tuned by the amount of impregnating agent, DAHP, as the decisive properties that affect the adsorption kinetics.Due to the careful synthesis design, the parametric space in which the adsorption process can be investigated is significantly reduced, allowing for the derivation of more straightforward conclusions about the effects of adsorbent properties on the adsorption kinetics.

Performance of the Materials
In the forthcoming sections, the adsorption data for chlorpyrifos, dimethoate, and malathion are analyzed.First, the adsorption was analyzed for up to 60 min, and the concentration of the pesticide remaining in the solution at a given moment was evaluated using UPLC analysis (Section 2, see Figures S1-S3, Supplementary Materials).Then, the amount of adsorbed pesticide in a given moment (q t ) was fitted into the kinetic equations, assuming kinetics obeying pseudo-first-(Equation ( 1)) and pseudo-second-order (Equation ( 2)): ln q e − q t ) = ln q e − k 1 t (1) where q e is the equilibrium adsorption capacity, while k 1 and k 2 are the pseudo-first and pseudo-second rate order constants.The obtained rate order constants for chlorpyrifos, dimethoate, and malathion are presented in Tables 2-4.From the value of k 1 , the adsorption half-time (t 1/2 ) can be calculated as C 2023, 9, 103 8 of 15 which gives the time at which half of the pesticide adsorption capacity is taken up by the adsorbent.

Chlorpyrifos Removal
The rate constants for chlorpyrifos removal are given in Table 2.It can be seen that the rate constants for pseudo-first-order reaction kinetics are determined with large uncertainties.In contrast, rate constants are determined more reliably for the pseudosecond-order kinetics.This is partially contributed to by the fact that a linearized form of the equation for the second-order kinetics was used.While linearized forms are easier to handle, they artificially increase the correlation coefficient [33].However, for this investigation, the overall trends and their connection with the materials' properties are of interest rather than determining the adsorption rate constants with high accuracy.Apparently, there is no association between the chlorpyrifos removal constants and the overall trends in the materials' properties (Table 1)-an overall increase in total pore volume and an increase in O and P atomic content.Although the reliability of the pseudo-first-order constants is not high, it is possible to conclude that the half-times for chlorpyrifos adsorption are of the order of 10 min.This result aligns well with previous adsorption kinetics studies of different pesticides on carbon materials [34][35][36][37].Minute ranges of chlorpyrifos adsorption halftimes were reported for gamma radiation-modified activated carbon [38], while graphitic carbon nitride (g-C 3 N 4 )-incorporated chitosan showed chlorpyrifos adsorption kinetics very similar to the presented ACFs [39].Nevertheless, the half-times obtained here were much higher compared to, for example, the one found for chlorpyrifos removal by an NU-1000 metal-organic framework, which was only 0.43 min [40].

Dimethoate and Malathion Removal
In contrast to chlorpyrifos, which possesses an aromatic ring in its structure, dimethoate and malathion are aliphatic molecules.Nevertheless, the same general observations that are made for chlorpyrifos hold for dimethoate and malathion adsorption kinetics, as well.There are no apparent trends, and k 2 is determined with higher reliability than k 1 .The half-times for dimethoate adsorption are in the range of 3.4-17.8min.For malathion, the range is 1-12.6 min.Similarly to chlorpyrifos, the estimated half-times for dimethoate adsorption on studied ACFs are in the same range or lower, compared to those found in the literature or calculated based on reported k 1 , as with the case of KOH-modified Thevetia peruviana shell-activated carbon [41].The studied ACFs are very efficient adsorbents for malathion removal and also possess better adsorption kinetics than some other carbons, like waste-derived activated carbons [26].

Discussion
Multiple linear regression analysis was used to investigate the impact of material properties on the adsorbents' performances.It is a relatively simple tool but can reveal some hidden links between dependent and independent variables, which can only be easily understood with the help of statistical tools.In short, it is assumed that the rate constants k 1 or k 2 can be presented as a linear combination of the properties (x i ) given in Table 1 with corresponding coefficients (A i ) as Counter i is 1-4 for pore volumes up to 1, 2, 3, and 4 nm, respectively, 5 for total pore volume, and 6-8 for C, O, and P atomic content.The results of the linear regression analysis are given in Figure 3. Before the analysis, k 1 and k 2 were multiplied by 100.At the same time, the C and O contents were divided by 100 and 10, respectively, to scale all the quantities used in the analysis to similar orders of magnitude.To analyze the correlation between the various input features given in Table 1 and the pseudo-second-order adsorption rate constants of three studied pesticides, Pearson correlation coefficients are calculated to depict the degree of correlation between the variables (Figure 4).To analyze the impact of different materials' properties, one should also consider the properties of the adsorbates.Malathion, dimethoate, and chlorpyrifos are of similar sizes, around 1 nm along the longest axis of a molecule, but with some differences in branching.Also, chlorpyrifos has an aromatic ring [29].Thus, the accommodation of all three molecules into the pore structure of studied adsorbents should start from pores with a diameter above 1 nm.However, it could be even a little larger as both adsorbent and adsorbate are solvated, which effectively reduces the pore diameter while increasing the dimensions of the adsorbate.Now, if the results of the linear regression analysis are considered (Figure 3 and Table A1), it is evident that generally good predictions for adsorption rate constants are obtained.This is particularly true for pseudo-second-order kinetics in the cases of dimethoate and malathion.In both cases, three independent variables stand out.First, the model is susceptible to the cumulative pore volumes for pores with diameters up to 2 nm and 4 nm, which positively correlate to the rate constants.Second, the carbon content (Table 1) also seems to be a determinant of fast adsorption kinetics.As discussed before [24,29], there are indications that for the studied materials, physisorption is operative.Thus, a positive correlation of carbon content with adsorption rate constants can be understood through the reduced number of highly solvated domains.However, O atomic content also positively correlates with the adsorption rate constants (much more so for malathion than dimethoate).It can be due to dipole-dipole [42] and electrostatic interactions, which can positively affect the adsorption rates in accordance with the charges of the adsorbate molecules and adsorbent.If their charges are opposite, the adsorption will be promoted and the adsorption rates will be enhanced [43,44].If some highly oxidized domains are present at the pore openings, they might block the entrance of pesticides into the pores, reducing the adsorption capacity and impeding the adsorption process rate.
To analyze the correlation between the various input features given in Table 1 and the pseudo-second-order adsorption rate constants of three studied pesticides, Pearson correlation coefficients are calculated to depict the degree of correlation between the variables (Figure 4).   1) and pseudo-second-order rate constants for pesticide adsorption (Tables 2-4).The intensity of red indicates the strength of a positive correlation, while the intensity of blue indicates the strength of a negative correlation.Abbreviations used: %DAHPpercentage (X) of DAHP in the impregnation step, VYnm-pore volume up to the width of Y nm (Y = 1, …, 4), elemental content in at.%, and k2-pseudo-second-order adsorption rate constant of malathion (MLT), dimethoate (DMT), and chlorpyrifos (CPF).
Regarding the correlation of the studied material-related properties, we find that all of them, except the V2nm, are strongly correlated and significantly correlated with the pseudo-second-order adsorption rate constant for chlorpyrifos (k2,CHP).However, we do not find any significant correlation to the other two studied pesticides with these features.Importantly, %DAHP is negatively correlated to V1nm, which might be explained by the increasing activation rate at higher DAHP content sizes of the molecule.The pore development in correlation with the impregnation ratio and activation yield has been discussed in detail previously [45].A moderate positive correlation between the used DAHP percentage and k2,CHP allows for the training of linear regression models with this feature as the only input.It depicts the important aspects of this pretreatment to the rate of chlorpyrifos removal.Based on the Pearson correlation coefficients ≥0.6, we can relate the success of DAHP as a modifier in this case to the decrease in V1nm and changes in the sample's elemental composition.The lack of significant correlation of any input feature to k2,DMT and k2,MLT might also be due to the low capacity of pesticide removal in these cases.
It is interesting to observe the role of the cumulative pore volumes based on the mul-  1) and pseudo-second-order rate constants for pesticide adsorption (Tables 2-4).The intensity of red indicates the strength of a positive correlation, while the intensity of blue indicates the strength of a negative correlation.Abbreviations used: %DAHP-percentage (X) of DAHP in the impregnation step, V Ynm -pore volume up to the width of Y nm (Y = 1, . .., 4), elemental content in at.%, and k 2 -pseudo-second-order adsorption rate constant of malathion (MLT), dimethoate (DMT), and chlorpyrifos (CPF).
Regarding the correlation of the studied material-related properties, we find that all of them, except the V 2nm , are strongly correlated and significantly correlated with the pseudo-second-order adsorption rate constant for chlorpyrifos (k 2,CHP ).However, we do not find any significant correlation to the other two studied pesticides with these features.Importantly, %DAHP is negatively correlated to V 1nm , which might be explained by the increasing activation rate at higher DAHP content sizes of the molecule.The pore development in correlation with the impregnation ratio and activation yield has been discussed in detail previously [45].A moderate positive correlation between the used DAHP percentage and k 2,CHP allows for the training of linear regression models with this C 2023, 9, 103 11 of 15 feature as the only input.It depicts the important aspects of this pretreatment to the rate of chlorpyrifos removal.Based on the Pearson correlation coefficients ≥0.6, we can relate the success of DAHP as a modifier in this case to the decrease in V 1nm and changes in the sample's elemental composition.The lack of significant correlation of any input feature to k 2,DMT and k 2,MLT might also be due to the low capacity of pesticide removal in these cases.
It is interesting to observe the role of the cumulative pore volumes based on the multiple linear regression models.Namely, there are indications that the first limit corresponds to the entrance of individual molecules (including their solvation shells) into the pores in such a way that one molecule interacts with two opposite sides of the pore interior along the pore size diameter (Figure 5a).The second situation (4 nm) corresponds to the case where the adsorbate molecule interacts with one side of the pore interior, but the pore diameter is large enough to accommodate two adsorbate molecules along the pore size diameter (Figure 5b).The diffusion of organophosphate molecules into the pore system of adsorbents is likely to be spontaneous and without the significant deformation of pesticide molecules, which would produce a large energy penalty during this process.To check for this assumption, semi-empirical quantum chemical calculations on the chlorpyrifos inserted into the interior of a (7,7) single-walled carbon nanotube (SWCNT) were performed, imitating the interior of a pore of the studied ACFs.Such a SWCNT has a diameter of 9.5 Å (if the length is infinite).For the finite-sized SWCNT studied in this case, the aliphatic part of the chlorpyrifos molecule is "expelled" from the tube (Figure 6a).At the same time, both chlorpyrifos and the nanotube undergo significant deformation (Figure 6b).The deformation energy of the chlorpyrifos molecule is 33 kJ mol −1 , while for SWCNT, it amounts to 710 kJ mol −1 (please note that this energy relates to 1 mol of the finite-size SWCNTs used in this study).Even without a solvent shell, the entrance of chlorpyrifos into pores smaller than 1 nm is unlikely.Moreover, one should acknowledge that the studied pesticides do not have molecules which are perfect spheres.For example, in the case of chlorpyrifos, presented here, the longer axis of the molecule is approx.1.08 nm, while the molecule dimension perpendicular to this axis is approx.0.85 nm.Moreover, polar parts of the molecules become more heavily solvated compared to non-polar ones, which additionally complicates the situation.A particular orientation of the molecules is absolutely necessary for them to enter the pores, but at this point it is unclear whether molecular reorientation and proper placement at the pore openings can be the rate-determining step for adsorption in the pores.An additional point relates to the surface functional groups at the pore openings, which can affect molecular reorientation and/or partial desolvation during adsorption in the pores.Finally, it should be noted that the rate constants obtained here generally do not correlate with the adsorption capacities [24,27].Hence, one should consider the properties that determine the materials' adsorption capacities and fast adsorption kinetics under the desired conditions to design proper adsorbents for different pesticides.The diffusion of organophosphate molecules into the pore system of adsorbents is likely to be spontaneous and without the significant deformation of pesticide molecules, which would produce a large energy penalty during this process.To check for this assumption, semi-empirical quantum chemical calculations on the chlorpyrifos inserted into the interior of a (7,7) single-walled carbon nanotube (SWCNT) were performed, imitating the interior of a pore of the studied ACFs.Such a SWCNT has a diameter of 9.5 Å (if the length is infinite).For the finite-sized SWCNT studied in this case, the aliphatic part of the chlorpyrifos molecule is "expelled" from the tube (Figure 6a).At the same time, both chlorpyrifos and the nanotube undergo significant deformation (Figure 6b).The deformation energy of the chlorpyrifos molecule is 33 kJ mol −1 , while for SWCNT, it amounts to 710 kJ mol −1 (please note that this energy relates to 1 mol of the finite-size SWCNTs used in this study).Even without a solvent shell, the entrance of chlorpyrifos into pores smaller than 1 nm is unlikely.Moreover, one should acknowledge that the studied pesticides do not have molecules which are perfect spheres.For example, in the case of chlorpyrifos, presented here, the longer axis of the molecule is approx.1.08 nm, while the molecule dimension perpendicular to this axis is approx.0.85 nm.Moreover, polar parts of the molecules become more heavily solvated compared to non-polar ones, which additionally complicates the situation.A particular orientation of the molecules is absolutely necessary for them to enter the pores, but at this point it is unclear whether molecular reorientation and proper placement at the pore openings can be the rate-determining step for adsorption in the pores.An additional point relates to the surface functional groups at the pore openings, which can affect molecular reorientation and/or partial desolvation during adsorption in the pores.Finally, it should be noted that the rate constants obtained here generally do not correlate with the adsorption capacities [24,27].Hence, one should consider the properties that determine the materials' adsorption capacities and fast adsorption kinetics under the desired conditions to design proper adsorbents for different pesticides.

Conclusions
In this study, the kinetics of chlorpyrifos, dimethoate, and malathion adsorption on a series of activated carbon fibers with diverse pore structures and chemical compositions were analyzed.The adsorption data were processed to obtain pseudo-first-and pseudosecond-order rate constants.No particular trend was found in connection with the studied materials' properties alone (integral pore volumes, the content of C, O, P).However, employing multiple linear regression allowed for the establishment of connections between the materials' properties and the adsorption kinetics.Specifically, for dimethoate and malathion, the regression analysis provided reliable predictions for the pseudo-second-order kinetics rate constants, while for chlorpyrifos, the performance of the regression model was lower.Among the considered material properties, cumulative pore volumes up to 2 nm and 4 nm and the carbon content showed strong positive correlations with the rate constants.This observation suggests an interplay between the sizes of pesticide molecules and pore diameters, enabling the accommodation of one or two pesticide molecules within the pore diameter.Theoretical calculations further indicated that, in the absence of a solvent, the penetration of pesticide molecules into pores with a diameter below 1 nm is unlikely without significant deformation of the molecule (and the pore).Such a process would require a considerable energy input, making it thermodynamically unfavorable.To optimize the design or selection of adsorbents for particular applications, it is essential to strike a balance between the properties governing the adsorption capacities and kinetics, as these two aspects are not directly correlated.The ideal adsorbent for the pesticides under study would be a mesoporous carbon material featuring a pore diameter of around 4 nm, along with a unimodal pore size distribution.Additionally, it should possess a relatively low concentration of oxygen (and potentially other heteroatoms) while maintaining a significant presence of sp 2 domains.The findings in this study provide valuable insights into the critical factors influencing the adsorption kinetics for these pesticides, which can contribute to the development and optimization of efficient adsorbents in various environmental and industrial applications.

Supplementary Materials:
The following supporting information can be downloaded at www.mdpi.com/xxx/s1, Figure S1: PDA signal of malathion (5 × 10 −4 mol dm −3 ) with chromatogram and the UV-Vis spectrum at the selected retention time.

Conclusions
In this study, the kinetics of chlorpyrifos, dimethoate, and malathion adsorption on a series of activated carbon fibers with diverse pore structures and chemical compositions were analyzed.The adsorption data were processed to obtain pseudo-first-and pseudosecond-order rate constants.No particular trend was found in connection with the studied materials' properties alone (integral pore volumes, the content of C, O, P).However, employing multiple linear regression allowed for the establishment of connections between the materials' properties and the adsorption kinetics.Specifically, for dimethoate and malathion, the regression analysis provided reliable predictions for the pseudo-secondorder kinetics rate constants, while for chlorpyrifos, the performance of the regression model was lower.Among the considered material properties, cumulative pore volumes up to 2 nm and 4 nm and the carbon content showed strong positive correlations with the rate constants.This observation suggests an interplay between the sizes of pesticide molecules and pore diameters, enabling the accommodation of one or two pesticide molecules within the pore diameter.Theoretical calculations further indicated that, in the absence of a solvent, the penetration of pesticide molecules into pores with a diameter below 1 nm is unlikely without significant deformation of the molecule (and the pore).Such a process would require a considerable energy input, making it thermodynamically unfavorable.To optimize the design or selection of adsorbents for particular applications, it is essential to strike a balance between the properties governing the adsorption capacities and kinetics, as these two aspects are not directly correlated.The ideal adsorbent for the pesticides under study would be a mesoporous carbon material featuring a pore diameter of around 4 nm, along with a unimodal pore size distribution.Additionally, it should possess a relatively low concentration of oxygen (and potentially other heteroatoms) while maintaining a significant presence of sp 2 domains.The findings in this study provide valuable insights into the critical factors influencing the adsorption kinetics for these pesticides, which can contribute to the development and optimization of efficient adsorbents in various environmental and industrial applications.

Figure 1 .
Figure 1.(a) Integral pore volumes up to a given pore diameter function; (b) common pore size distribution curves for studied adsorbents.The notation of samples is DAHP-X, where X is given in the figure legend.

Figure 1 .
Figure 1.(a) Integral pore volumes up to a given pore diameter function; (b) common pore size distribution curves for studied adsorbents.The notation of samples is DAHP-X, where X is given in the figure legend.

C
2023, 9, x FOR PEER REVIEW 11 of 16

Figure 4 .
Figure 4. Heat map of the Pearson correlation matrix coefficients between various features, including material-related properties (Table1) and pseudo-second-order rate constants for pesticide adsorption (Tables2-4).The intensity of red indicates the strength of a positive correlation, while the intensity of blue indicates the strength of a negative correlation.Abbreviations used: %DAHPpercentage (X) of DAHP in the impregnation step, VYnm-pore volume up to the width of Y nm (Y = 1, …, 4), elemental content in at.%, and k2-pseudo-second-order adsorption rate constant of malathion (MLT), dimethoate (DMT), and chlorpyrifos (CPF).

Figure 4 .
Figure 4. Heat map of the Pearson correlation matrix coefficients between various features, including material-related properties (Table1) and pseudo-second-order rate constants for pesticide adsorption (Tables2-4).The intensity of red indicates the strength of a positive correlation, while the intensity of blue indicates the strength of a negative correlation.Abbreviations used: %DAHP-percentage (X) of DAHP in the impregnation step, V Ynm -pore volume up to the width of Y nm (Y = 1, . .., 4), elemental content in at.%, and k 2 -pseudo-second-order adsorption rate constant of malathion (MLT), dimethoate (DMT), and chlorpyrifos (CPF).

C 16 Figure 5 .
Figure 5. Schematic representation of the accommodation of chlorpyrifos molecules into pores of different diameters: (a) when the pore diameter allows for the entrance of only one molecule, (b) when the pore diameter is sufficient to place two chlorpyrifos molecules in the pore.Any pores with diameters between the hypothetical cases of (a,b) can only accommodate one molecule along the pore diameter.

Figure 5 .
Figure 5. Schematic representation of the accommodation of chlorpyrifos molecules into pores of different diameters: (a) when the pore diameter allows for the entrance of only one molecule, (b) when the pore diameter is sufficient to place two chlorpyrifos molecules in the pore.Any pores with diameters between the hypothetical cases of (a,b) can only accommodate one molecule along the pore diameter.

C 16 Figure 6 .
Figure 6.Entrance of a chlorpyrifos molecule into the SWCNT: (a) side view, (b) a view along the SWCNT (pore) axis; the deformation is visible as the SWCNT deviates from the circular shape.

Figure 6 .
Figure 6.Entrance of a chlorpyrifos molecule into the SWCNT: (a) side view, (b) a view along the SWCNT (pore) axis; the deformation is visible as the SWCNT deviates from the circular shape.

Table 1 .
Adsorbents' properties-pore volume integrated up to a given pore diameter, total pore volume (integrated up to pore diameter of 35 nm), surface area, and elemental content, determined using energy dispersive X-ray spectroscopy.The notation of samples is DAHP-X, where X represents the concentration of DAHP in wt.% used for the impregnation step.

Table 2 .
Rate constants for chlorpyrifos removal.The notation of samples is DAHP-X, where X represents the concentration of DAHP in wt.% used for the impregnation step.The chlorpyrifos concentration was 1 × 10 −4 mol dm −3 , while the adsorbents' concentration was 0.1 mg mL −1 in each case.

Table 3 .
Rate constants for dimethoate removal.The notation of samples is DAHP-X, where X represents the concentration of DAHP in wt.% used for the impregnation step.The dimethoate concentration was 1 × 10 −4 mol dm −3 , while the adsorbents' concentration was 0.1 mg mL −1 in each case.

Table 4 .
Rate constants for malathion removal.The notation of samples is DAHP-X, where X represents the concentration of DAHP in wt.% used for the impregnation step.For the adsorption measurements, the concentration of malathion was 1 × 10 −4 mol dm −3 , while the adsorbents' concentration was 0.1 mg mL −1 in each case.