Parabens Adsorption onto Activated Carbon: Relation with Chemical and Structural Properties

Parabens (alkyl-p-hidroxybenzoates) are antimicrobial preservatives used in personal care products, classified as an endocrine disruptor, so they are considered emerging contaminants. A raw version of activated carbons obtained from African palm shell (Elaeis guineensis) modified chemically by impregnation with salts of CaCl2 (GC2), MgCl2 (GM2) and Cu(NO3)2 (GCu2) at 2% wt/v and carbonized in CO2 atmosphere at 1173 K was prepared. The process of adsorption of methyl (MePB) and ethylparaben (EtPB) from aqueous solution on the activated carbons at 18 °C was studied and related to the interactions between the adsorbate and the adsorbent, which can be quantified through the determination of immersion enthalpies in aqueous solutions of corresponding paraben, showing the lowest-value carbon GM2, which has a surface area of 608 m2 × g−1, while the highest values correspond to the activated carbon GCu2, with a surface area of 896 m2 × g−1 and the highest content of surface acid sites (0.42 mmol × g−1), such as lactonic and phenolic compounds, which indicates that the adsorbate–adsorbent interactions are favored by the presence of these, with interaction enthalpies that vary between 5.72 and 51.95 J × g−1 for MePB adsorption and 1.24 and 52.38 J × g−1 for EtPB adsorption showing that the process is endothermic.


Introduction
Parabens (alkyl-p-hydroxybenzoates) are antimicrobial preservatives whose use is widely distributed, mainly in personal care products. They are referred to as emerging contaminants whose harmful effects are related to inhibition of normal hormonal system activity in animals and presumably in humans, so they are classified as an endocrine disruptor [1]. Methyl and ethylparaben are generally considered to have a much lower potential to cause endocrine disrupting effects compared to other larger parabens [2,3]. Błędzka et al. [4] in their study showed the limited accessible data on the presence of parabens in surface waters. In general, among all measured parabens, the highest frequency and concentration for methylparaben was recorded, while the presence of ethylparaben in water samples was recorded less frequently.
Currently, there are several processes that are used in the removal of endocrine disruptors; these include coagulation/flocculation [5,6], biological treatment [5][6][7], oxidation [5, 6,8], and adsorption [6,9,10]. Some of these techniques are more traditional and used than others. While new methods or improvements are continually proposed, each process has advantages, disadvantages, and limitations in the removal of contaminant traces according to local conditions, which does not ensure an adequate technical-economic-environmental outcome in all cases. The particular case of adsorption on activated carbon of different organic compounds has been extensively studied or has become a viable alternative for the removal of various contaminants in water [11][12][13]. An important

Activated Carbon Characterization
Nitrogen adsorption isotherms at 77 K are present in Figure 1A, where the modification of the textural characteristics produced by chemical impregnation treatments with metal salts is evident. Samples of GM2 and GCu2 display a typical type I adsorption isotherm, while the GC2 sample exhibits a behavior composed of isotherms type I and II according to the IUPAC classification, showing a more pronounced adsorption at low pressures relative to micropore filling, while at high pressures, capillary condensation is presented, accompanied in all isotherms by small hysteresis loops of desorption, indicating the poor development of mesoporosity in activated carbons.
The pore size distribution calculated by the density functional theory (DFT) method is shown in Figure 1B. The best fitting was obtained using the quenched solid density functional theory (QSDFT), where the interaction between the adsorbate and the adsorbent is explicitly taken into account and is considered a "rough" or "heterogeneous" surface [22,23]. For the samples studied, the best results were obtained considering a combined porous system slit/cylinder with average error percentages between 0.023% and 0.158%. In this figure, the microporous structure of activated carbons is confirmed where the average pore size is between 0.3 and 1 nm, having a greater contribution of micropores around The pore size distribution calculated by the density functional theory (DFT) method is shown in Figure 1B. The best fitting was obtained using the quenched solid density functional theory (QSDFT), where the interaction between the adsorbate and the adsorbent is explicitly taken into account and is considered a "rough" or "heterogeneous" surface [22,23]. For the samples studied, the best results were obtained considering a combined porous system slit/cylinder with average error percentages between 0.023% and 0.158%. In this figure, the microporous structure of activated carbons is confirmed where the average pore size is between 0.3 and 1 nm, having a greater contribution of micropores around 0.4-0.6 nm, and the almost null development of mesoporosity is confirmed. Therefore, the mesopore volume calculation usually determined by the BJH (Barrett, Joyner, and Halenda) method was omitted.
The pore structure parameters are present in Table 1, showing that impregnation with metallic salts moderately modifies the textural characteristics of activated carbons, where it is shown that, all carbons are microporous, because the microporosity corresponds practically to the total pore volume in the activated carbons. The relationship between the surface area of the carbons and the volume of micropores developed is also evident, considering that most of the total surface area is in micropores. These results can be explained by considering the smaller atomic radius that is in copper (0.128 nm), which favors the formation of microporosity, compared to those of magnesium (0.160 nm) and calcium (0.197 nm). This confirms that the reagent used in the impregnation is used as a template for the creation of porosity [24]. The results of functional group analysis of activated carbons determined by the technique proposed by Boehm are presented in Table 2. In this case, it was observed that in all cases, the treatments carried out changed the chemical characteristics of the activated carbons. In the case of total acidity, it was observed that impregnation with copper salt develops a higher content of acid icing groups on the periphery of the graphenic layers of GCu2 activated carbon produced by some oxygenated groups of carboxylic, lactonic, and phenolic, particularly the latter, because the African palm shell is constituted by lignin, cellulose, and hemicellulose [25], which are high in this group. As for total basicity, it is evident that impregnation with calcium salt favored the formation of basic The pore structure parameters are present in Table 1, showing that impregnation with metallic salts moderately modifies the textural characteristics of activated carbons, where it is shown that, all carbons are microporous, because the microporosity corresponds practically to the total pore volume in the activated carbons. The relationship between the surface area of the carbons and the volume of micropores developed is also evident, considering that most of the total surface area is in micropores. These results can be explained by considering the smaller atomic radius that is in copper (0.128 nm), which favors the formation of microporosity, compared to those of magnesium (0.160 nm) and calcium (0.197 nm). This confirms that the reagent used in the impregnation is used as a template for the creation of porosity [24]. The results of functional group analysis of activated carbons determined by the technique proposed by Boehm are presented in Table 2. In this case, it was observed that in all cases, the treatments carried out changed the chemical characteristics of the activated carbons. In the case of total acidity, it was observed that impregnation with copper salt develops a higher content of acid icing groups on the periphery of the graphenic layers of GCu2 activated carbon produced by some oxygenated groups of carboxylic, lactonic, and phenolic, particularly the latter, because the African palm shell is constituted by lignin, cellulose, and hemicellulose [25], which are high in this group. As for total basicity, it is evident that impregnation with calcium salt favored the formation of basic groups associated with pyrone and chromene structures, as well as the presence of electrons π on basal graphene layers [26] in carbon GC2.

Adsorption Studies
It is important that adsorption is a very complex process, depending on the characteristics of the adsorbent as well as adsorbates, as well as affinity that is presented in said system, affecting the adsorption capacity of a given adsorbent. Table 3 presents the physicochemical properties of the parabens under study; they have low water solubility and increasing the length of the carbonated chain leads to decreased the solubility, so sodium salts are frequently used also in formulations. The values of the acid constant of dissociation (pK a ) is around 8.2, so in aqueous solutions, they are in their protonate form. groups associated with pyrone and chromene structures, as well as the presence of electrons π on basal graphene layers [26] in carbon GC2.

Adsorption Studies
It is important that adsorption is a very complex process, depending on the characteristics of the adsorbent as well as adsorbates, as well as affinity that is presented in said system, affecting the adsorption capacity of a given adsorbent. Table 3 presents the physicochemical properties of the parabens under study; they have low water solubility and increasing the length of the carbonated chain leads to decreased the solubility, so sodium salts are frequently used also in formulations. The values of the acid constant of dissociation (pKa) is around 8.2, so in aqueous solutions, they are in their protonate form. Giles et al. [28] related the shape of the adsorption isotherms in the aqueous phase with the nature of the adsorption. According to their classification, the adsorption isotherms of methyl and ethylparaben presented in Figure 2 are identified as L-type or Langmuir isotherms, with a concave shape so that as the concentration in the liquid phase increases, the amount absorbed slowly increases.  groups associated with pyrone and chromene structures, as well as the presence of electrons π on basal graphene layers [26] in carbon GC2.

Adsorption Studies
It is important that adsorption is a very complex process, depending on the characteristics of the adsorbent as well as adsorbates, as well as affinity that is presented in said system, affecting the adsorption capacity of a given adsorbent. Table 3 presents the physicochemical properties of the parabens under study; they have low water solubility and increasing the length of the carbonated chain leads to decreased the solubility, so sodium salts are frequently used also in formulations. The values of the acid constant of dissociation (pKa) is around 8.2, so in aqueous solutions, they are in their protonate form. Giles et al. [28] related the shape of the adsorption isotherms in the aqueous phase with the nature of the adsorption. According to their classification, the adsorption isotherms of methyl and ethylparaben presented in Figure 2 are identified as L-type or Langmuir isotherms, with a concave shape so that as the concentration in the liquid phase increases, the amount absorbed slowly increases.  Giles et al. [28] related the shape of the adsorption isotherms in the aqueous phase with the nature of the adsorption. According to their classification, the adsorption isotherms of methyl and ethylparaben presented in Figure 2 are identified as L-type or Langmuir isotherms, with a concave shape so that as the concentration in the liquid phase increases, the amount absorbed slowly increases. groups associated with pyrone and chromene structures, as well as the presence of electrons π on basal graphene layers [26] in carbon GC2.

Adsorption Studies
It is important that adsorption is a very complex process, depending on the characteristics of the adsorbent as well as adsorbates, as well as affinity that is presented in said system, affecting the adsorption capacity of a given adsorbent. Table 3 presents the physicochemical properties of the parabens under study; they have low water solubility and increasing the length of the carbonated chain leads to decreased the solubility, so sodium salts are frequently used also in formulations. The values of the acid constant of dissociation (pKa) is around 8.2, so in aqueous solutions, they are in their protonate form. Giles et al. [28] related the shape of the adsorption isotherms in the aqueous phase with the nature of the adsorption. According to their classification, the adsorption isotherms of methyl and ethylparaben presented in Figure 2 are identified as L-type or Langmuir isotherms, with a concave shape so that as the concentration in the liquid phase increases, the amount absorbed slowly increases.  To carry out the study of the interactions between parabens and activated carbons from an aqueous solution, adsorption data were described using the Sips model (Equation (1)). This three-parameter model is an extension of the Freundlich model that is appropriate for describing the adsorption of heterogeneous systems, but in its mathematical form, it is like the Langmuir model, which is applicable to ideal surfaces. The main difference between the Langmuir model and the Sips model is the n S parameter that describes the characteristic heterogeneity of the system derived from the adsorbent, adsorbate or a combination of both [29]: Q e is the adsorbed amount in the equilibrium (mg × g −1 ), C e is the concentration of the adsorbate in the equilibrium (mg × L −1 ), Q mS is Sips maximum adsorption capacity of (mg × g −1 ), and K S is the Sips equilibrium constant (L × mg −1 ) for heterogeneous solids that relates to the characteristic energy of micropores.
The parameters obtained from the adjustment to the Sips model are presented in Table 4, showing that the methylparaben adsorption values are lower than those obtained for ethylparaben, as is the case for the adsorption of EtPB with the carbon GCu2 that, possessing the highest value of the constant K S , indicates that there is a greater energy interaction between EtPB with the surface of this activated carbon, which also has the highest concentration of oxygenated functional groups on its surface. With respect to the n S parameter, it is indicated that the systems formed by the EtPB with the GC2 and GM2 carbons are energetically homogeneous because their value is exactly or is approximately to unity [29], while the system formed by EtPB with the GCu2 carbon by moving away from this value suggests that the system is energetically heterogeneous. As is well known, two regions are distinguished on the surface: a nonpolar and hydrophobic one that has basic properties associated with regions rich in π electrons located in the graphene basal planes and another hydrophilic region containing various functional groups [30,31]. Thus, the carbons GC2 and GM2 possess the lowest acidity and highest basicity, as shown in Figure 3, which is usually associated with π-delocalized electrons in the graphenic layers of the carbonous structure, and have dispersive hydrophobic interactions [32] with the electrons in the aromatic ring of the ethylparaben [33,34]. Carbon GCu2, meanwhile, having a higher concentration of acid groups that decrease the electron density of the basal plans of the activated carbon, favors the specific interactions that are established between these groups and the phenol present in the structure EtPB. These results suggest that in general, the presence of functional groups promotes heterogeneous interactions [35]. However, this heterogeneity of the system indicates not only adsorbent-adsorbate interactions but possible adsorbate-adsorbate associations on the surface that increase adsorption.  Similarly, the systems formed by MePB with the carbons GC2 and GM2 have an energetically homogeneous character, their value being less than the unit and for the reasons mentioned above with the EtPB, given the similarity in their structures. However, the system formed by MePB-GCu2 is also homogeneous, suggesting that a type of interaction prevails. Hence, surface chemistry may be deciding that either specific or nonspecific interaction is the prevailing mechanism for adsorption of a given adsorbate [31]. Figure 3 also shows the influence of surface chemistry on paraben adsorption, showing an almost direct interaction between the adsorbed paraben molecules and the acidic groups present on the surface, indicating that these groups are responsible for the adsorption evidenced in the different activated carbons, through the formation of specific interactions between the acidic groups on the periphery of the graphenic layers of activated carbon and the phenol groups present in paraben structures. With regard to adsorption capacities, it is shown that this is higher for EtPB, which is confirmed by various materials studied [36,37]. This is due to the decrease in solubility, as the size of the carbon chain in the ester group increases, so interaction with water is favored by decreasing the polarity of the EtPB. Therefore, paraben-adsorbent interaction will be stronger than that of parabensolvent as the molecular size increases.
As shown in the figure above, adsorption depends on the presence of surface groups, particularly of acid character, and expanding this information in Figure 4, it is evidenced that it depends on the interactions that are generated with phenolic and lactonic groups, through the formation of hydrogen bonds with the phenol of parabens and, on the other hand, a "donoracceptor" complex mechanism similar to that mentioned by Mattson and co-workers [17] between the carbonyl group on the carbon surface and the aromatic ring of the solute, with the carbonyl oxygen of lactone acting as the electron donor and the aromatic ring of the paraben structure acting as the acceptor. This indicates that the paraben adsorption mechanism involves several kinds of electron-donors as carbonyl-oxygen-containing functional groups (carboxylic/lactone) on the activated carbon surface and its interaction with the paraben aromatic ring. Therefore, given the similarity in the content of lactonic groups of the GC2 and GM2 carbons, similarities are presented with the adsorption capacity between these carbons. Similarly, the systems formed by MePB with the carbons GC2 and GM2 have an energetically homogeneous character, their value being less than the unit and for the reasons mentioned above with the EtPB, given the similarity in their structures. However, the system formed by MePB-GCu2 is also homogeneous, suggesting that a type of interaction prevails. Hence, surface chemistry may be deciding that either specific or nonspecific interaction is the prevailing mechanism for adsorption of a given adsorbate [31]. Figure 3 also shows the influence of surface chemistry on paraben adsorption, showing an almost direct interaction between the adsorbed paraben molecules and the acidic groups present on the surface, indicating that these groups are responsible for the adsorption evidenced in the different activated carbons, through the formation of specific interactions between the acidic groups on the periphery of the graphenic layers of activated carbon and the phenol groups present in paraben structures. With regard to adsorption capacities, it is shown that this is higher for EtPB, which is confirmed by various materials studied [36,37]. This is due to the decrease in solubility, as the size of the carbon chain in the ester group increases, so interaction with water is favored by decreasing the polarity of the EtPB. Therefore, paraben-adsorbent interaction will be stronger than that of paraben-solvent as the molecular size increases.
As shown in the figure above, adsorption depends on the presence of surface groups, particularly of acid character, and expanding this information in Figure 4, it is evidenced that it depends on the interactions that are generated with phenolic and lactonic groups, through the formation of hydrogen bonds with the phenol of parabens and, on the other hand, a "donor-acceptor" complex mechanism similar to that mentioned by Mattson and co-workers [17] between the carbonyl group on the carbon surface and the aromatic ring of the solute, with the carbonyl oxygen of lactone acting as the electron donor and the aromatic ring of the paraben structure acting as the acceptor. This indicates that the paraben adsorption mechanism involves several kinds of electron-donors as carbonyl-oxygen-containing functional groups (carboxylic/lactone) on the activated carbon surface and its interaction with the paraben aromatic ring. Therefore, given the similarity in the content of lactonic groups of the GC2 and GM2 carbons, similarities are presented with the adsorption capacity between these carbons.

Calorimetric Study
The calorimetric study of activated carbons was carried out through the determination of immersion enthalpy in water and paraben solutions, in order to observe changes in enthalpy by concentration as well as impregnation with the different metal salts. Table 5 presents the results obtained from the immersion enthalpy, the values of which correspond to exothermic processes of the order of physisorption. Regarding water immersion enthalpy, its relationship is observed with the increase in hydrophilic functional groups at the edges of the graphenic layers, which increase the polarity of activated carbons in the order of GCu2 > GC2 > GM2. By associating these results with those of chemical characterization, it is evident that this order relates to the content of phenolic groups, as shown in Figure 5, so it suggests that adsorbent-solvent interaction is associated with the formation of a hydrogen bond between phenol and water.

Calorimetric Study
The calorimetric study of activated carbons was carried out through the determination of immersion enthalpy in water and paraben solutions, in order to observe changes in enthalpy by concentration as well as impregnation with the different metal salts. Table 5 presents the results obtained from the immersion enthalpy, the values of which correspond to exothermic processes of the order of physisorption. Regarding water immersion enthalpy, its relationship is observed with the increase in hydrophilic functional groups at the edges of the graphenic layers, which increase the polarity of activated carbons in the order of GCu2 > GC2 > GM2. By associating these results with those of chemical characterization, it is evident that this order relates to the content of phenolic groups, as shown in Figure 5, so it suggests that adsorbent-solvent interaction is associated with the formation of a hydrogen bond between phenol and water.
With regard to immersion enthalpy in paraben solutions, which correspond to the result of total energy produced in adsorbent-adsorbate, adsorbent-solvent, and adsorbent-solvent interactions, it is evident that it decreases with the increase in paraben concentration, suggesting that the energy transferred in the form of heat in each system decreases as the surface of the adsorbent becomes saturated [35], as shown in Figure 6A. With regard to immersion enthalpy in paraben solutions, which correspond to the result of total energy produced in adsorbent-adsorbate, adsorbent-solvent, and adsorbent-solvent interactions, it is evident that it decreases with the increase in paraben concentration, suggesting that the energy transferred in the form of heat in each system decreases as the surface of the adsorbent becomes saturated [35], as shown in Figure 6A. On the other hand, in order to evaluate the calorimetric effect of the specific interactions between parabens and activated carbons, the interaction enthalpy (∆Hint) to both concentrations was determined, finding values from 5.72 to 51.95 J × g −1 for MePB and 1.24 to 52.38 J × g −1 for EtPB. As shown in Figure 6B, changes in the interaction enthalpy have positive values in all cases, suggesting that the process is endothermic and requires energy to be carried out; hence, it is usually associated with the breaking of interactions and water displacement by paraben molecules, which are much more apolar, and as the size of paraben molecule increases, the polarity of this decreases from the carbonaceous surface and the de-solvatation of the solute before adsorption [33]. Therefore, the EtPBadsorbent interaction is stronger than that of MePB-adsorbent, so that π-π-type interactions are favored. With regard to immersion enthalpy in paraben solutions, which correspond to the result of total energy produced in adsorbent-adsorbate, adsorbent-solvent, and adsorbent-solvent interactions, it is evident that it decreases with the increase in paraben concentration, suggesting that the energy transferred in the form of heat in each system decreases as the surface of the adsorbent becomes saturated [35], as shown in Figure 6A. On the other hand, in order to evaluate the calorimetric effect of the specific interactions between parabens and activated carbons, the interaction enthalpy (∆Hint) to both concentrations was determined, finding values from 5.72 to 51.95 J × g −1 for MePB and 1.24 to 52.38 J × g −1 for EtPB. As shown in Figure 6B, changes in the interaction enthalpy have positive values in all cases, suggesting that the process is endothermic and requires energy to be carried out; hence, it is usually associated with the breaking of interactions and water displacement by paraben molecules, which are much more apolar, and as the size of paraben molecule increases, the polarity of this decreases from the carbonaceous surface and the de-solvatation of the solute before adsorption [33]. Therefore, the EtPBadsorbent interaction is stronger than that of MePB-adsorbent, so that π-π-type interactions are favored. On the other hand, in order to evaluate the calorimetric effect of the specific interactions between parabens and activated carbons, the interaction enthalpy (∆H int ) to both concentrations was determined, finding values from 5.72 to 51.95 J × g −1 for MePB and 1.24 to 52.38 J × g −1 for EtPB. As shown in Figure 6B, changes in the interaction enthalpy have positive values in all cases, suggesting that the process is endothermic and requires energy to be carried out; hence, it is usually associated with the breaking of interactions and water displacement by paraben molecules, which are much more apolar, and as the size of paraben molecule increases, the polarity of this decreases from the carbonaceous surface and the de-solvatation of the solute before adsorption [33]. Therefore, the EtPB-adsorbent interaction is stronger than that of MePB-adsorbent, so that π-π-type interactions are favored. Figure 7 shows the relation between the lactonic groups and the microporosity developed in activated carbons, where it is shown that the microporous structure favors interactions between parabens, because they allow the distance between the solute and the active sites on the surface of the activated carbon to be shortened, thus increasing the possibility of interactions between these, and therefore, the intensity of interactions increases. These active sites on the surface may be associated with the presence of lactonic groups, which are susceptible to hydrolysis in acidic medium, forming a carboxylic acid and an alcohol, which would favor the formation of a greater number of hydrogen bridges with the phenol of the adsorbate, while phenol formation favors that the solvent not compete for the adsorption sites created as specific solvent-phenol interactions would be formed [33]. However, although the relationship with microporosity is shown, it should be noted that given the pore size distribution ( Figure 1B) of activated carbons that had a concentration of micropores in the region between 0.3 and 0.9 nm, it is suggested that there can be some diffuse restrictions if paraben molecules do not have adequate guidance to be able to enter micropores, in which case adsorption will depend more on surface chemistry.
Molecules 2019, 24, x; doi: FOR PEER REVIEW www.mdpi.com/journal/molecules Figure 7 shows the relation between the lactonic groups and the microporosity developed in activated carbons, where it is shown that the microporous structure favors interactions between parabens, because they allow the distance between the solute and the active sites on the surface of the activated carbon to be shortened, thus increasing the possibility of interactions between these, and therefore, the intensity of interactions increases. These active sites on the surface may be associated with the presence of lactonic groups, which are susceptible to hydrolysis in acidic medium, forming a carboxylic acid and an alcohol, which would favor the formation of a greater number of hydrogen bridges with the phenol of the adsorbate, while phenol formation favors that the solvent not compete for the adsorption sites created as specific solvent-phenol interactions would be formed [33]. However, although the relationship with microporosity is shown, it should be noted that given the pore size distribution ( Figure 1B) of activated carbons that had a concentration of micropores in the region between 0.3 and 0.9 nm, it is suggested that there can be some diffuse restrictions if paraben molecules do not have adequate guidance to be able to enter micropores, in which case adsorption will depend more on surface chemistry.

Preparation of Activated Carbons
A raw version of activated carbons obtained from African palm shell (Elaeis guineensis) was modified chemically by impregnation with salts of CaCl2 (GC2), MgCl2 (GM2), and Cu(NO3)2 (GCu2) at 2% wt/v in a thermostat bath at 358 K for 48 h. The carbonization of the samples was carried out in a CARBOLITE ® horizontal tubular furnace (Carbolite-Gero, Hope, UK) in CO2 atmosphere (110 mL × min −1 ) with a 5 K × min −1 heating ramp for 6 h at 1173 K [38,39]. Subsequently, these were washed with HCl 5% and water to remove remain salts until the pH was constant. Finally, the activated carbons were dried in an oven for 24 h at 363 K and stored in plastic hermetically sealed containers in a nitrogen atmosphere.

Characterization of Activated Carbons
The textural characteristics of activated carbons were determined from N2 adsorption isotherms at 77 K, in an automatic equipment Autosorb 3B Quantachrome (Quantachrome Instruments, Boynton Beach, FL, USA).The apparent surface area (SBET) was evaluated via application of the BET

Preparation of Activated Carbons
A raw version of activated carbons obtained from African palm shell (Elaeis guineensis) was modified chemically by impregnation with salts of CaCl 2 (GC2), MgCl 2 (GM2), and Cu(NO 3 ) 2 (GCu2) at 2% wt/v in a thermostat bath at 358 K for 48 h. The carbonization of the samples was carried out in a CARBOLITE ® horizontal tubular furnace (Carbolite-Gero, Hope, UK) in CO 2 atmosphere (110 mL × min −1 ) with a 5 K × min −1 heating ramp for 6 h at 1173 K [38,39]. Subsequently, these were washed with HCl 5% and water to remove remain salts until the pH was constant. Finally, the activated carbons were dried in an oven for 24 h at 363 K and stored in plastic hermetically sealed containers in a nitrogen atmosphere.

Characterization of Activated Carbons
The textural characteristics of activated carbons were determined from N 2 adsorption isotherms at 77 K, in an automatic equipment Autosorb 3B Quantachrome (Quantachrome Instruments, Boynton Beach, FL, USA).The apparent surface area (S BET ) was evaluated via application of the BET model, the volume of micropore (V o ) was calculated using the Dubinin-Radushkevich equation (D-R), and the QSDFT model was applied to nitrogen isotherms for the determination of total volume (V T ) and the pore size distribution using ASiQWin software (Quantachrome Instruments, Boynton Beach, FL, USA, version 3.01), assuming cylinder/slit geometry over the adsorption branch.
Chemical characterization was realized using the Boehm method [30,40], with bases of different strength that allow neutralizing acids with different pKa and then determine their individual contributions. For this, 0.500 g of activated carbon was weighed and mixed with 25 mL of NaOH, HCl, Na 2 CO 3 , and NaHCO 3 0.05 M solutions (Merck Millipore, Bedford, MA, USA), respectively, in hermetic plastic containers. The solutions were left in agitation at 18 • C for 5 days and under 100 rpm constant agitation, with the occasional bubbling of N 2 in order to remove atmospheric CO 2 . After some time, an aliquot of each supernatant solution was certified with the respective acid or base solutions standardized in an automatic titrator (TitroLine Alpha-Plus, Schott Instruments GnbH, Mainz, Germany). The acidic functional groups were determined from the titration of NaOH, Na 2 CO 3 , and NaHCO 3 supernatant solutions with standardized HCl. For the quantification of the basis groups, the titration of the HCl solution with standardized NaOH solution was continued. The HCl and NaOH solutions were standardized with boric acid and potassium biphthalate, respectively. Titration experiments were performed in triplicate, and mean values were presented.

Adsorption Studies
Adsorption studies were carried out using 0.025 g of activated carbons prepared previously weighted on an analytical balance with an accuracy of 0.001 g (Ohaus Pioneer PA 114, Ohaus Corporation, Parsippany, NJ, USA), with 50 mL of methylparaben (MePB) and ethylparaben (EtPB) solutions, analyte grade > 99% purity (Alfa Aesar, Tewksbury, MA, USA), in a range of 20 to 200 mg × L −1 for 3 weeks (based on a previous kinetic study, where it was shown that the equilibrium of the system was reached after this time) at 18 • C, with sporadic agitation (100 rpm). The residual concentration of paraben was measured via UV-Vis Spectroscopy (GENESYS 10uv, spectrophotometer, Thermo Fisher Scientific, Madison, WI, USA) at 254 nm [33,37,41]: The paraben adsorbed amount (Q e ) on the activated carbons was calculated from the initial concentration (C o ), equilibrium concentration (C e ), paraben solution total volume (V), and mass of activated carbon (m), as presented in the balance of Equation (2).

Immersion Enthalpy Determination
Immersion enthalpies of activated carbons were determined, using as wet liquids water and aqueous solutions of MePB and EtPB (20 y 200 mg × L −1 , respectively) at 18 • C on a locally built Calvet-type heat conduction microcalorimeter [42], using a stainless-steel calorimetric cell. About 0.100 g of each activated carbon was weighed into a glass vial located in the calorimetric cell containing 10 mL of the wet liquid; the electrical potential was measured for approximately 40 min until a stable baseline was obtained. The immersion of the sample and the increase in potential was recorded according to the time and product of the wet of the solid until it returned back to the baseline. The determinations were performed in triplicate, the standard deviations were calculated, and the results are shown as means values. Finally, electrical calibration was performed by providing known electrical work to the calorimetric cell [43,44]. This calibration allows calculating such electrical work (W elect ), in addition to the calorimeter constant (K cal ), immersion energy or heat (E imm ), and finally, the immersion enthalpy (∆H imm ), as shown in Equations (3)- (6).
K cal = W elect Area under calibration peak (4) E imm = K cal * Area under immersion peak (5) − ∆H imm = E imm (J) activated carbon mass (g) On the other hand, the enthalpy of interaction (∆H int ) corresponds to energy produced by the contact between the adsorbate and the adsorbent (∆H ads−AC ), regardless of solvent-solvent interactions, and is determined by the Hess Law from the immersion enthalpy in paraben solutions (∆H imm ) and solvent-activated carbon interactions (∆H solv−AC ): Therefore, the interaction enthalpy depends not only on the adsorbent-adsorbate interactions but also on the breakdown of adsorbate-solvent and solvent-adsorbent interactions, as presented in Equation (7).

Conclusions
In this work were presented the results of the adsorption of two parabens, methyl and ethylparaben, in three activated carbons obtained from African palm shell through chemical activation with different metallic salts. The adsorption of the parabens under study can be related to the presence of acid groups, particularly lactonic groups by a complex "donor-acceptor" mechanism between the carbonyl oxygen of lactone and the aromatic ring of paraben, while the concentration of phenolic groups is more associated with the formation of hydrogen bonds with water, favoring adsorbent-solvent interactions. It was shown that both homogeneous and heterogeneous interactions may occur, due to the presence of adsorbate-adsorbent interactions, as well as adsorbent-solvent interactions. However, it was also shown that the adsorption is favored with the development of microporosity, because they allow the distance between the solute and the active sites on the surface of the activated carbon to be shortened, thus increasing the possibility of interactions, in the event that no diffuse restrictions are introduced due to the size of the solutes into the smaller pores, in which case, the surface chemistry takes on greater importance in adsorption.
The calorimetric study determined that the values of immersion enthalpy at both low and high paraben concentrations are lower than those found for the solvent (from −13.10 to −68.45 J × g −1 ), indicating the formation of solvent-adsorbent interactions and for the adsorption process to be carried out on the adsorbents, the displacement of the solvent by the adsorbate is needed, so it requires energy to carry out the rupture of the interactions with water; thus, positive interaction enthalpy values were calculated.