Methylparaben Adsorption onto Activated Carbon and Activated Olive Stones: Comparative Analysis of Efﬁciency, Equilibrium, Kinetics and Effect of Graphene-Based Nanomaterials Addition

Featured Application: The content of this paper allows a positive assessment of the possibilities of using activated olive stones as a substitute adsorbent for activated carbon in the elimination of organic pollutants


Introduction
During the last decades, there has been an important growth in the world population, which has generated a great industrial expansion and a considerable technological development that has led to an increase in the presence of substances in the environment that can cause harmful effects on living beings [1]. Among these substances are the so-called emerging pollutants or pollutants of emerging concern, substances that are not regulated or controlled by national or international environmental authorities, whose presence in the environment is not necessarily new but about which there is growing concern for their potential effects on human health and ecosystems [2,3]. They include compounds as varied in nature as pharmaceuticals, lifestyle products, drugs of abuse, personal care products, Methylparaben (methyl 4-hydroxybenzoate, 99%) was supplied by Alfa Aesar, Kandel (Germany). Activated carbon Darco AC 4-12 (surface area 600-675 m 2 /gm) was obtained from Sigma-Aldrich, Steinheim (Germany). Crushed raw olive stones (ranging in size from 1 to 4 mm) were supplied by the olive oil mill Valle de Ricote, Archena, Murcia (Spain), and 36% HCl was supplied by Panreac, Darmstadt (Germany). Graphene oxide (99%) and reduced graphene oxide (80% C) were obtained from Abalonyx (Oslo, Norway).

Adsorbents Preparation
Commercial activated carbon (AC) was repeatedly washed with distilled water to remove the carbon dust it contained until the washing water was completely transparent. Then, it was dried in an oven at 105 • C for 24 h and stored in a desiccator.
Crushed olive stones were washed with distilled water and dried at 60 • C in an oven. In order to increase their active surface area, remove soluble components (tannins, resins, reducing sugars and dyes) and generate oxygenated functional groups (carboxyl, phenol, ketone and lactone), the olive stones were subjected to an activation process consisting of calcination at 300 • C for 3 h in a closed muffle furnace and subsequent treatment with 1 M HCl at room temperature (25 ± 2 • C) for 8 h [27,28]. Finally, the activated olive stones (AOS) were filtered, washed thoroughly with distilled water, dried at 100 • C for 24 h and stored in a desiccator. The surface area of the prepared activated olive stones was 479 m 2 /g [27].
Graphene-based nanocomposite-modified adsorbents were prepared from suspensions of 0.1% GO and rGO in water, obtained by applying ultrasound for 30 min at 10 min intervals separated by 5 min in an ice bath to avoid heating. A total of 50 mL of these suspensions were added to 150 mg of the adsorbents and sonicated again in the same way. Once the process was finished, the adsorbents were separated by settling, washed several times with distilled water, dried in an oven for 24 h and stored in a desiccator.
Batch experiments were performed in 50 mL Erlenmeyer flasks containing 40 mL of MPB solution of selected concentrations and pH and the appropriate amount of sorbent. The flasks were shaken at the selected speed on a thermostated shaker (20 • C) for 24 h. At that time, 1 mL of sample was taken and 2 mL of 1 M NaOH was added (basic pH provides a higher intensity of the MPB absorption band), determining MPB concentration by UV-Vis spectrophotometry at 289 nm using an Agilent 8453 spectrophotometer (Agilent Technologies, Santa Clara, CA, USA) from the calibration curve (concentration range: 10-60 mg/L; R 2 = 0.9997). The results obtained showed a maximum deviation of 4%.

Equilibrium and Kinetic Experiments
In the equilibrium and kinetic experiments, the adsorption studies were carried out by shaking 0.02 g of AC or 0.06 g of AOS with 40 mL of MPB solutions of different concentrations (10-50 mg/L) for 24 h (equilibrium time) at pH 7 in a rotary shaker (50 mL Erlenmeyer flasks) at 240 rpm and at 20 • C, taking samples at different predetermined times and determining MPB concentration as described in the previous section.
The amount of MPB loaded onto the adsorbents at equilibrium, q e (mg/g), and at any time t, q t (mg/g), was determined by the following equations: where C 0 , C t and C e were the initial, time t and equilibrium MPB concentrations in the solution (mg/L), respectively; V was the volume of the MPB solution (L); and m was the mass of the adsorbent (g). The data obtained from these experiments were used to analyze their fit to the different isotherm, kinetic and adsorption mechanism models by determining the corresponding correlation coefficients.

Mathematical Models Equilibrium Studies
An adsorption isotherm represents an equilibrium relationship at a given temperature between the amount of adsorbate in the liquid phase and that on the adsorbent surface. In this paper, we analyze the fit of the equilibrium experimental data of MPB adsorption onto AC and AOS to six adsorption isotherm models (Langmuir, Freundlich, Elovich, Temkin, Jovanovic and Dubinin-Radushkevich) to determine which model best describes those processes (Table 1).

Kinetic Studies
Adsorption kinetics describes the rate of retention of an adsorbate in a fluid environment on the surface of the adsorbent, determining the time required to reach equilibrium.
Four kinetic models (Lagergren pseudo-first order, Ho pseudo-second order, Elovich and Avrami) were used to determine the one that best describes methylparaben adsorption onto activated carbon and activated olive stones (Table 2). Table 2. Adsorption kinetics mathematical models.

Adsorption Mechanism
The above kinetic models are not able to identify the diffusion mechanisms and the rate-controlling step that affect the adsorption process. Any solid-liquid sorption process is characterized by the following two mass transfer processes: 1. The adsorbate molecules diffuse from the bulk solution to the external surface of the adsorbent (boundary layer diffusion); 2. The adsorbate molecules diffuse through the interior of the adsorbent particles (intraparticle diffusion).
Any one of them, or a combination of both [51], can control the rate of the adsorption process.
The Weber-Morris intraparticle diffusion model and the Boyd model are used to identify the rate-controlling step of the adsorption process (Table 3). Table 3. Adsorption mechanism mathematical models.

Model Linearized Equation Plot Model Parameters
Weber-Morris [52] q t = k intp · t 1/2 + C i q t against t 1/2 k intp : intraparticle diffusion rate constant (mg/g·h 1/2 ) C i : contribution of the boundary layer diffusion B·t against t When pure intraparticle diffusion exists, the Weber-Morris representation shows a straight line passing through the origin. The presence of multilinearity in that representation means that a combination of both boundary layer diffusion and intraparticle diffusion controls the adsorption process.
If the latter occurs, the Boyd kinetic model is used to identify the rate-controlling step of the adsorption process. If the plot of B·t versus time shows a straight line passing through the origin, the adsorption process is mainly controlled by intraparticle diffusion; otherwise, it is mainly controlled by boundary layer diffusion.

Effect of the Initial Methylparaben Concentration
Increasing the initial concentration of methylparaben from 10 to 50 mg/L results in a significant increase in the adsorption capacity values of both adsorbents, from 18.87 mg/g to 77.18 mg/g in the case of activated carbon and from 4.33 mg/g to 15.95 mg/g in the case of activated olive stones ( Figure 1a). This increase must be a result of the increased driving force to overcome the resistance to methylparaben mass transfer between the aqueous and the solid phases, which enhances the interaction between the methylparaben and the adsorbents [54].

Effect of Adsorbent Dose
The amount of adsorbent used in the adsorption process gives the cost of the adsorbent per unit volume of the solution to be treated.
Increasing the adsorbent dose produces a decrease in the adsorption capacity values of both adsorbents (Figure 1b), from 95.58 to 39.73 (58.43% decrease) in the case of activated carbon (increase in dose from 10 mg to 50 mg) and from 27.23 to 12.25 (55.01% decrease) in the case of activated olive stones (increase in dose from 20 mg to 100 mg).
This decrease indicates that, in the adsorbent mass ranges used, although the mass of methylparaben adsorbed increases with the amount of adsorbent, as the total surface area available for adsorption (and, consequently, the total number of available adsorption sites) increases, the ratio of mass of methylparaben adsorbed per gram of adsorbent decreases as the initial concentration of methylparaben used remains unchanged [55].
For the same adsorbent dose (20 mg), the adsorption capacity of methylparaben on activated carbon (77.18 mg/g) is 4.05 times higher than its adsorption capacity on activated olive stones (19.07 mg/g) due to the greater surface area and the greater proportion of functional groups of the former with respect to the latter. Appl

Effect of Adsorbent Dose
The amount of adsorbent used in the adsorption process gives the cost of the adsorbent per unit volume of the solution to be treated.
Increasing the adsorbent dose produces a decrease in the adsorption capacity values of both adsorbents (Figure 1b), from 95.58 to 39.73 (58.43% decrease) in the case of activated carbon (increase in dose from 10 mg to 50 mg) and from 27.23 to 12.25 (55.01% decrease) in the case of activated olive stones (increase in dose from 20 mg to 100 mg).
This decrease indicates that, in the adsorbent mass ranges used, although the mass of methylparaben adsorbed increases with the amount of adsorbent, as the total surface area available for adsorption (and, consequently, the total number of available adsorption sites) increases, the ratio of mass of methylparaben adsorbed per gram of adsorbent decreases as the initial concentration of methylparaben used remains unchanged [55].
For the same adsorbent dose (20 mg), the adsorption capacity of methylparaben on activated carbon (77.18 mg/g) is 4.05 times higher than its adsorption capacity on activated olive stones (19.07 mg/g) due to the greater surface area and the greater proportion of functional groups of the former with respect to the latter.

Effect of Stirring Speed
The effect of stirring speed on the removal of methylparaben by adsorption on activated carbon and activated olive stones is shown in Figure 1c. As can be seen, by increasing the stirring speed from 80 rpm to 240 rpm, the adsorption capacity of methylparaben increases from 17.26 mg/g to 77.18 mg/g (4.47 times) in the case of activated carbon and from 2.75 to 15.95 (5.80 times) in the case of activated olive stones. This increase is a consequence of the increased diffusion of methylparaben from inside the solution to the surface and inside the pores of the adsorbent.

Effect of the Initial pH of Methylparaben Solution.
The pH of the adsorbate solution is a parameter that has a very important influence on the adsorption process by determining the ionization state of the functional groups and, consequently, the charge of both the adsorbate and the adsorbent. The effect of the

Effect of Stirring Speed
The effect of stirring speed on the removal of methylparaben by adsorption on activated carbon and activated olive stones is shown in Figure 1c. As can be seen, by increasing the stirring speed from 80 rpm to 240 rpm, the adsorption capacity of methylparaben increases from 17.26 mg/g to 77.18 mg/g (4.47 times) in the case of activated carbon and from 2.75 to 15.95 (5.80 times) in the case of activated olive stones. This increase is a consequence of the increased diffusion of methylparaben from inside the solution to the surface and inside the pores of the adsorbent.

Effect of the Initial pH of Methylparaben Solution
The pH of the adsorbate solution is a parameter that has a very important influence on the adsorption process by determining the ionization state of the functional groups and, consequently, the charge of both the adsorbate and the adsorbent. The effect of the initial pH of the methylparaben solution on its adsorption on activated carbon and activated olive stones is shown in Figure 1d. It can be observed that the adsorption capacity of methylparaben on both adsorbents decreases with the increase in pH, with this decrease being more pronounced at basic pH values, especially in the case of activated carbon.
Methylparaben has a pKa value of 8.17 [7], which means that at values below this pH, it has no charge, and at higher values, it has a negative charge as a consequence of the ionization of the phenolic hydroxyl group.
The commercial activated carbon used has a zero charge point of 4.7 [56], which means that it will have a positive charge at pH values below 4.7 and a negative charge at higher pH values. Activated olive stones have a zero charge point of 3.5 [28], being positively charged at pH values below 3.5 and negatively charged at pH values above 3.5. With regard to the functional groups present on the surface of these adsorbents, the presence of carboxyl, phenol, anhydride, lactone, quinone and ether groups has been described in both activated carbon [57,58] and activated olive stones [28,59].
At pH 3, as methylparaben is in a neutral state and the two adsorbents are positively charged, some electrostatic interactions may occur between the positive zones of the adsorbents and the oxygens of the hydroxyl and ester groups of methylparaben (oxygens with a certain density of negative charge due to their inductive-I effect). This results in adsorption capacity values of 89.66 mg/g for activated carbon and 18.26 mg/g for activated olive stones. At pH 7, methylparaben is in a neutral state and the two adsorbents are negatively charged, with the possibility of establishing hydrogen bridge interactions between the hydrogen of the hydroxyl group of the methylparaben and the negatively charged groups of the adsorbents, but with no possibility of the electrostatic interactions described in the previous paragraph. All this causes the adsorption capacity to decrease to values of 77.18 mg/g for AC (13.92% decrease) and 15.95 mg/g for AOS (12.65% decrease).
Finally, at pH 11, both methylparaben and the adsorbents are negatively charged, so there will be a significant electrostatic repulsion between the adsorbate and adsorbents, resulting in a significant decrease in the adsorption capacity of methylparaben, which decreases to values of 39.32 mg/g in AC (40.05% decrease in the value at pH 7 and 56.15% decrease in the value at pH 3) and of 10.03 mg/g in AOS (37.12% decrease in the value at pH 7 and 45.07% decrease in the value at pH 3).

Selected Standard Experimental Conditions
The selected standard experimental conditions were as follows: MPB initial concentration: 50 mg/L; adsorbent dose: 20 mg/L of AC and 60 mg/L of AOS; stirring speed: 240 rpm; initial methylparaben solution pH: 7.

Adsorption Equilibrium
In order to find out which isotherm model best describes the adsorption of methylparaben on activated carbon and activated olive stones, the experimental equilibrium data were fitted to the linearized mathematical equations of each of the six equilibrium models described above (Langmuir, Freundlich, Elovich, Temkin, Jovanovic and Dubinin-Radushkevich) by analyzing the correlation coefficients obtained in the corresponding linear regression. Figures 2 and 3 show the linear representation of these six models for activated carbon and activated olive stones, respectively, while the values of their model constants and correlation coefficients (R 2 ) are included in Tables 4 and 5, respectively.      The Langmuir isotherm best describes the adsorption process of methylparaben on both activated carbon and activated olive stones, showing the highest correlation coefficient values in all cases. The fit of the experimental data to the Langmuir isotherm model suggests that in the two adsorbents, the adsorption sites are homogeneously distributed and energetically equivalent, each adsorption site adsorbs only one methylparaben molecule and there is no interaction between the adsorbed methylparaben molecules. Moreover, the values of the R L parameter are less than 1 in both adsorbents, indicating the favorable character of methylparaben adsorption on both activated carbon and activated olive stones [31].
The analysis of the Langmuir model constants for both adsorbents allows some significant conclusions to be drawn. The maximum monolayer adsorption capacity (q m ) of methylparaben on activated carbon is 3.4 times higher than its value on activated olive stones, and the increase in temperature in both adsorbents indicates that the adsorption of methylparaben on both activated carbon and activated olive stones is endothermic.
The values of the Langmuir constant (K L ), related to the adsorbate-adsorbent affinity, allow us to affirm that the MPB/AC affinity is 6.85 times higher than the MPB/AOS affinity.
Although the Dubinin-Radushkevich model does not fit as well as the Langmuir model to the experimental data, the mean free energy of adsorption calculated from that model is less than eight in both adsorbents, suggesting the physical character of the adsorption process [42] and confirming the nature of the interactions described above in explaining the variations of adsorption with pH.

Adsorption Kinetics
As in the case of the isotherms, in order to determine the kinetic model that best describes the adsorption of methylparaben on activated carbon and activated olive stones, the experimental data on the variation of adsorption capacity with time were fitted to the linearized mathematical equations of the four kinetic models studied (Lagergren pseudofirst order, Ho pseudo-second order, Elovich and Avrami), and the correlation coefficients obtained were comparatively analyzed.
Linear representations of these four models for activated carbon and activated olive stones are shown in Figure 4, while the kinetic parameters and R 2 values calculated from these kinetic models are included in Tables 6 and 7.
These results allow us to state that while in the case of the activated olive stones, the experimental data clearly fit better to a pseudo-first-order kinetic model as the correlation coefficients obtained in this model are clearly higher than in the other three models, in the case of activated carbon, the values of the correlation coefficients of the pseudo-first and pseudo-second-order models are very similar, so it is not possible to establish a definitive conclusion. In order to discern which of these two models best defines the adsorption process of methylparaben on activated carbon, the greater/lesser coincidence of the experimental values of q e with the values of this parameter obtained from pseudo-first and pseudo-second-order models was analyzed. Figure 5 shows this coincidence for both activated carbon and activated olive stones. It is observed that the theoretical q e values obtained by the pseudo-first-order model practically coincide with the experimental values of q e ; therefore, it can be concluded that the pseudo-first-order model is the one that best describes the adsorption of methylparaben on both activated carbon and activated olive stones.
As in the case of the isotherms, in order to determine the kinetic model that best describes the adsorption of methylparaben on activated carbon and activated olive stones, the experimental data on the variation of adsorption capacity with time were fitted to the linearized mathematical equations of the four kinetic models studied (Lagergren pseudofirst order, Ho pseudo-second order, Elovich and Avrami), and the correlation coefficients obtained were comparatively analyzed.
Linear representations of these four models for activated carbon and activated olive stones are shown in Figure 4, while the kinetic parameters and R 2 values calculated from these kinetic models are included in Tables 6 and 7. These results allow us to state that while in the case of the activated olive stones, the experimental data clearly fit better to a pseudo-first-order kinetic model as the correlation coefficients obtained in this model are clearly higher than in the other three models, in the case of activated carbon, the values of the correlation coefficients of the pseudo-first and pseudo-second-order models are very similar, so it is not possible to establish a definitive conclusion. In order to discern which of these two models best defines the adsorption process of methylparaben on activated carbon, the greater/lesser coincidence of the experimental values of qe with the values of this parameter obtained from pseudo-first and pseudo-second-order models was analyzed.     The analysis of the model constants for both adsorbents shows that, in the adsorption of methylparaben on activated carbon, the pseudo-first-order rate constant and the equilibrium adsorption capacity calculated by the model are 1.19 and 4.68 times, respectively, higher than those obtained in the adsorption of methylparaben on activated olive stones, presumably because activated carbon has a larger surface area and a greater number of active sites (functional groups) than activated olive stones. Figure 5 shows this coincidence for both activated carbon and activated olive stones. It is observed that the theoretical qe values obtained by the pseudo-first-order model practically coincide with the experimental values of qe; therefore, it can be concluded that the pseudo-first-order model is the one that best describes the adsorption of methylparaben on both activated carbon and activated olive stones. The analysis of the model constants for both adsorbents shows that, in the adsorption of methylparaben on activated carbon, the pseudo-first-order rate constant and the equilibrium adsorption capacity calculated by the model are 1.19 and 4.68 times, respectively, higher than those obtained in the adsorption of methylparaben on activated olive stones, presumably because activated carbon has a larger surface area and a greater number of active sites (functional groups) than activated olive stones.

Adsorption Mechanism
In order to identify the mechanism of methylparaben adsorption onto activated carbon and activated olive stones, the kinetic experimental data were further analyzed by the Weber-Morris intraparticle diffusion and Boyd models.
As can be seen in Figure 6, the intraparticle model representations are not linear over the whole time range, with three distinct linear zones clearly distinguishable, which means that intraparticle diffusion is not the only controlling step of the adsorption rate, with more than one controlling step involved. That is, the trial nature of the intraparticle model plots confirms that both boundary layer diffusion and intraparticle diffusion control the adsorption process in both adsorbents [60][61][62].
To distinguish the main rate-controlling step of the adsorption process, the experimental kinetic data were fitted to Boyd's kinetic model. Figure 6 shows that for both adsorbents, the plots are linear but do not pass through the origin, suggesting that the adsorption processes of methylparaben on both activated carbon and activated olive stones are mainly controlled by boundary layer diffusion [60][61][62].

Adsorption Mechanism
In order to identify the mechanism of methylparaben adsorption onto activated carbon and activated olive stones, the kinetic experimental data were further analyzed by the Weber-Morris intraparticle diffusion and Boyd models.
As can be seen in Figure 6, the intraparticle model representations are not linear over the whole time range, with three distinct linear zones clearly distinguishable, which means that intraparticle diffusion is not the only controlling step of the adsorption rate, with more than one controlling step involved. That is, the trial nature of the intraparticle model plots confirms that both boundary layer diffusion and intraparticle diffusion control the adsorption process in both adsorbents [60][61][62].

Effect of the Presence of Graphene-Based Nanomaterials in Adsorbents on Methylparaben Adsorption
Modification of activated carbon and activated olive stones by the addition of graphene oxide or reduced graphene oxide leads to an increase in the adsorption capacity of methylparaben on both adsorbents (Figure 7). To distinguish the main rate-controlling step of the adsorption process, the experimental kinetic data were fitted to Boyd's kinetic model. Figure 6 shows that for both adsorbents, the plots are linear but do not pass through the origin, suggesting that the adsorption processes of methylparaben on both activated carbon and activated olive stones are mainly controlled by boundary layer diffusion [60][61][62].

Effect of the Presence of Graphene-Based Nanomaterials in Adsorbents on Methylparaben Adsorption
Modification of activated carbon and activated olive stones by the addition of graphene oxide or reduced graphene oxide leads to an increase in the adsorption capacity of methylparaben on both adsorbents (Figure 7).  Figure 6. Adsorption mechanism plots for the adsorption of MPB onto AC and AOS.

Effect of the Presence of Graphene-Based Nanomaterials in Adsorbents on Methylparaben Adsorption
Modification of activated carbon and activated olive stones by the addition of graphene oxide or reduced graphene oxide leads to an increase in the adsorption capacity of methylparaben on both adsorbents (Figure 7). In the case of activated carbon, the addition of graphene oxide or reduced graphene oxide leads to increases in the methylparaben adsorption capacity of 16.6% and 5.1%, respectively, while in the case of activated olive stones, the increases are 11.8% and 4.0%, respectively. These increases must be the result of the modifications that the addition of In the case of activated carbon, the addition of graphene oxide or reduced graphene oxide leads to increases in the methylparaben adsorption capacity of 16.6% and 5.1%, respectively, while in the case of activated olive stones, the increases are 11.8% and 4.0%, respectively. These increases must be the result of the modifications that the addition of GO and rGO produces in the surface characteristics of activated charcoal and activated olive stones [63].
On the one hand, graphene-based nanomaterials are characterized by their high porosity and high surface area [63], so their presence on the surface of the adsorbents under study should increase their surface area and adsorption capacity.
On the other hand, GO is the oxidized derivative of graphene that maintains the polyaromatic p-electron system of graphene and contains oxygenated polar functional groups, including hydroxyl, epoxide, carbonyl and carboxyl groups, with the carbonyl and carboxyl groups located at the edges and the epoxide and hydroxyl groups located in the basal plane of GO. Reduced graphene oxide (rGO), as a reduced form of GO, contains the same basic structure as GO but with fewer oxygenated polar functional groups [64][65][66]. The additional contribution of oxygenated functional groups and polyaromatic p-electron systems that GO and rGO make to the surface of AC and AOS increases both the hydrogen bridge interactions (oxygenated functional groups of adsorbents/hydroxyl group of MPB) and the hydrophobic interactions (polyaromatic p-electron system of adsorbent/aromatic ring of MPB) between the adsorbents and MPB.
All this leads to the above-mentioned increase in the adsorption capacity of methylparaben by the modified adsorbents, which, for the reasons indicated, is higher in the modification with GO than with rGO. The greater increase in adsorption capacity observed in the CA with respect to the AOS must be a consequence of the greater amount of nanomaterial retained by the former due to its greater surface area.

Conclusions
A comparative study of the adsorption of methylparaben onto commercial activated carbon and olive stones activated by calcination at 300 • C and treatment with 1 M HCl was carried out in this paper. For both adsorbents, the methylparaben adsorption capacity increases with the increase in methylparaben concentration and stirring rate and decreases with the increase in adsorbent dose and pH.
At the same experimental conditions, the adsorption capacity of methylparaben onto activated carbon is about four times higher than that onto activated olive stones due to both higher surface area and higher MPB-AC adsorption interactions.
For both adsorbents, the methylparaben adsorption process is best described by the Langmuir isotherm model and the pseudo-first-order kinetic model, and it is mainly controlled by the boundary layer diffusion step. Langmuir model constants show that the maximum monolayer adsorption capacity of MPB onto AC is 3.4 times higher than that onto AOS and that the MPB/AC affinity is 6.85 times higher than the MPB/AOS affinity.
Pseudo-first-order kinetic model constants show that MPB adsorption onto AC has an adsorption rate constant and a theoretical adsorption capacity of 1.19 and 4.68, respectively, times higher than MPB adsorption onto AOS.
The addition of graphene oxide and reduced graphene oxide to both adsorbents increases their methylparaben adsorption capacity due to modifications in the surface characteristics of the adsorbents related to the increase in both adsorption surface area as a consequence of the high porosity and high surface area of graphene nanomaterials and the number of functional groups that can interact with the methylparaben through hydrogen bridge interactions (oxygenated functional groups) and hydrophobic interactions (polyaromatic p-electron functional groups). These increases are more important in the case of graphene oxide, due to its higher oxygenated functional group presence, and in the case of activated carbon, due to its higher surface area and higher oxygenated functional group presence.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to they are part of a much larger study that is still underway.

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