Physicochemical Characterization of Santa Barbara Amorphous-15 (SBA-15) and Its Functionalization with Polyaniline for Phenol Adsorption

The chemical surface of Santa Barbara Amorphous-15 (SBA-15) was functionalized with polyaniline (SBA-15/PA) and was prepared using aniline and ammonium persulfate as the modifying agent and oxidant. The samples were characterized by N2 physical adsorption at 77 K. SBA-15’s surface area was 654.5 m2 g−1, changing to 254.6 m2 g−1. The total acidity and basicity were determined using a modified Boehm’s method, and the point of zero charge pH (pHpzc) was determined. Immersion calorimetry in water, benzene, and NaOH, HCl, and phenol (100 mg L−1) solutions was performed, producing values ranging from −15.1 to −174.1 J g−1 related to the chemical surface, pHpzc, and acid and base group values. Finally, the phenol adsorption capacity was found to be higher for the modified polyaniline silica, at 138 mg g−1.


Introduction
According to the Environmental Protection Agency (EPA) [1], phenol is considered a priority pollutant, dangerous for the environment due to its high toxicity even at low concentrations. In Colombia, according to the Ministry of Environment and Sustainable Development in resolution 631 of 2015, the maximum permissible limit for phenol in wastewater discharges is 0.2 mg L −1 . The most common sources of phenol in water are effluents from industries that manufacture iron, rubber, resins, plastics, fibers, adhesives, and synthetic fuels. Phenol can be found in air and water because of the use and disposal of products that contain this substance; in the soil it tends to move into groundwater. Phenolcontaminated water sources have a bad taste and odor, so it is necessary to remove it in order to avoid repercussions for human health and the ecosystem [1,2]. According to the World Health Organization (WHO), the treatments to be carried out to produce safe water depend on the quality that is desired; a process to achieve basic levels of sanitation requires up to four barriers to eliminate traces of organic compounds: a typical water treatment scheme includes physical separation, passage through a membrane bioreactor or activated sludge, followed by treatment, such as advanced oxidation, activated carbon adsorption, or membrane filtration [2,3]. As mentioned in the first paragraph of this introduction, in Colombia, as in many Latin American countries, phenol is one of the compounds found in aquifer systems, generating major health problems for the entire population, especially the most unprotected, and mainly children. Several methods have been developed for the removal of phenol from water, such as photocatalytic degradation, electrochemical oxidation, solvent extraction, and biological degradation and adsorption, which is one of the most widely used methods for decontamination [4][5][6][7][8][9][10][11].
Adsorption on porous adsorbents is an alternative to processes like advanced oxidation, chlorination, and ozonation which involve reactive molecules, such as drugs, which can result in the production of more toxic compounds if agents such as hydrogen peroxide, ozone, and chlorine are present, giving rise to organochlorine compounds, furans, dioxins, and chloroform-molecules that have greater persistence and low degradability in the environment. Activated carbon is considered a non-reactive barrier that retains organic and inorganic compounds due to its physicochemical properties, such as high surface area, porosity, and the presence of surface chemical groups that allow it to interact with compounds that have several chemical characteristics [5][6][7].
According to scientific studies for Colombia, there is evidence for the presence of several kind of molecules, including those for pharmaceutical use, in water sources, generating a great problem that is concerning for state, environmental, political, and academic organizations in the context of evaluating solutions from different areas. One of the processes endorsed by the WHO in the guide for the production of safe water is adsorption on porous solids because it is considered an alternative to chemical processes that change the properties of pharmaceutical compounds, turning them into more toxic substances with greater difficulty of removal [7,10].
Colombia has a serious problem in the form of phenol contamination in water, which manifests itself by causing a noticeable unpleasant taste and odor change through the chlorination of water, despite the fact that its concentration is at ppm level and even less than that. Most phenols and their derivates are corrosive to skin and eye. In Colombia, as well as in other countries in the area, coal distillation and organic synthesis waste flows bearing coal entities contain large amounts of phenols and their derivatives, as well as phenolic compounds, which are present in wastewater from pulp and paper-bleaching facilities, and from resin, pesticide, insecticide, paint, and solvent industries. That is why it is very important to face this problem and generate solutions that can be scaled to the most vulnerable populations; it is a social duty for academia [8,9].
The adsorption process of phenolic compounds has been evaluated in several studies. The main conclusions shared by these works are that the adsorption process depends on factors such as the physical nature of the adsorbent (porous structure, surface chemistry), the physical nature of the adsorbate (pKa, polarity, size), and conditions of the solution (pH, ionic strength, adsorbate concentration) [2][3][4][5][6].
Regarding the adsorbent, factors such as geometric and chemical heterogeneity have been considered from the initial composition of the precursor, which affect the phenol adsorption process due to the presence of pores of different sizes that modify not only the kinetics of the process but also the adsorption capacity, since by increasing the volume of the micropores the amount of phenol adsorbed decreases [4]. On the other hand, the presence of narrow microporosity may not favor the process because phenol can form aggregates with water and increase its size, presenting diffusive restrictions in very small pores, as indicated by the kinetic data reported in the literature [5,7].
The analysis of the influence of surface chemistry on the adsorption capacity of phenol has also been evaluated. Studies conclude that activated carbon surfaces that have been subjected to oxidative treatments have lower adsorption capacities due to an increase in polarity in the surface which raises the affinity for the solvent, generating a blockage of the pores, so that the phenol interacts with the activated carbon through interactions with π electrons. However, an excess of oxidized groups can generate an inductive effect of electrons, also weakening the formation of interactions of this type [2,6].
Another of the postulated mechanisms for the adsorption of phenol on other adsorbents, such as activated carbons, is the electron donor-acceptor interaction, where the carbonyl groups of activated carbon are electron donors and the aromatic ring of phenol is the acceptor. Some authors have reported [7] the effect of heat treatment of activated carbons on phenol adsorption. According to the data, the process is favored in an activated carbon subjected to heat treatment at 1173 K, followed by treatment at 973 K, and finally at 1373 K. The author attributes this behavior to the increase in surface area caused by the increase in porous structures generated in the heat treatment; however, an excess of temperature leads to the collapse of the porous structure. When comparing the adsorption capacities of activated carbon without heat treatment and with treatment at 1173 K for different temperatures for the adsorption system (298 K, 313 K, 328 K), it was found that the adsorption capacity does not vary with temperature but the heat treatment for the sample favors the process (1.27 < 1.54 mmol g −1 ) [7][8][9][10][11].
On the other hand, other authors have reported [8] that the isosteric heat of adsorption decreased when the concentration increased from −78 kJ mol −1 to −5.5 kJ mol −1 ; likewise, the immersion enthalpies showed the same trend and varied between −2152 and 82.8 J g −1 .
At basic pH, the presence of oxygenated groups on activated carbon decreases phenol adsorption capacity due to a strong repulsive effect between phenolate and acidic chemical groups that are ionized with negative charge; however, oxygenated groups of a basic nature, such as pyrones, can contribute to form adsorbate-adsorbent interactions because they are not dissociated [2,9].
Another study [10] indicated that phenol adsorption involves three phases: adsorption at infinite dilution, filling of micropores, and filling of larger micropores. In the first phase, there is a correlation between the number of oxygenated groups and a decrease in adsorption capacity, which is associated with an increase in the affinity of the surface for the solvent, which is corroborated by an increase in the immersion enthalpy in water with respect to the amount of acid groups. This same author affirms that aromatic compounds, such as phenol, present three mechanisms by which adsorption on activated carbon is carried out from the aqueous phase: interactions by dispersion of π electrons, formation of hydrogen bonds, and formation of electron donor-acceptor complexes [11].
Most phenol-adsorbent interactions are due to physisorption; however, the oxidative coupling of phenol on activated carbon is associated with chemisorption.
From the thermodynamic point of view, the adsorption of phenol on activated carbon is considered exothermic, with enthalpy values between −6.03 and −30.4 kJ mol −1 , according to [5]. The process presents a low energy barrier to be carried out given that the activation energies are around 13 kJ mol −1 , which confirms that the process is a physisorption. The entropy change values are positive and depend on the surface chemistry of the adsorbent and the chemistry of the solution [10,11].
The adsorption processes have made it possible to address this type of research during the last few years. A great variety of materials, such as activated carbons [12][13][14][15][16], silica [12,17,18], polymeric resins [12,19,20], fly ash [21], clays including kaolinite [22], and zeolites [23][24][25][26][27], have been explored in detail for the removal of phenolic pollutants from wastewater. Among these materials, some authors [12] have reported that materials such as zeolite adsorbents have been regarded as an alternative to activated carbon due to their adsorption-regeneration properties. They say that although they have attracted special attention in sorption processes, their use in this field is limited by their pore size diameter with regards to the size of the organic pollutants [12]. As such, the use of mesoporous materials seems to be promising because of their larger pore volume and diameter, their high surface area, and their regular channel type structures [28,29]. The interest in the potential utilization of mesoporous materials as sorbents is increasing due to the need for economical and efficient adsorbents to remove organic contaminants from wastewater.
According to the literature, the chemical modification of mesoporous adsorbents by introducing functional groups into the matrixes enhances their adsorption capacity for phenol-containing materials [30]. For example, organically functionalized mesoporous silica materials, such as MCM41 and SBA-15, have already been reported by Brunel et al. and Stein et al. [12,31,32]. However, information on mesoporous materials interacting with toxic phenolic compounds is still limited.
In some research published by renowned researchers, Mangrulkar et al. [12,33] have used MCM-41 to remove phenol from an aqueous solution and observed that the removal of phenols from water is strongly dependent on the type of species (dissociated or molecular), as dissociated species are more polar than molecular compounds and interact preferentially with water rather than with MCM-41. Their works have allowed them to establish that the most important features involved in the investigation of adsorption phenomenon are: (i) interface characteristics, (ii) the adsorbate-adsorbent interactions, and (iii) adsorption isotherm [12]. It is worth mentioning that the physicochemical mechanisms in these materials are not very well described; no simple theory of adsorption could adequately describe experimental results. The adsorption phenomenon depends on interactions between the surface of the adsorbent and the adsorbed species. These interactions may be due to (i) chemical bonding, (ii) hydrogen bonding, (iii) hydrophobic bonding, or (iv) Van der Waals force. It was observed that the adsorption of substituted benzene depends on the electronic effect of the substituent. However, during the adsorption of dioxane on non-porous, mesoporous, and microporous silica, it was observed that adsorption on microporous silica is governed by pore volume and dimension.
SBA-15 material (Santa Barbara Amorphous No 15) is a mesoporous silica that has been increasing in importance due to its textural properties, regular mesoporous structure, high specific surface area, thermal and mechanical stability, distribution, and pore size. Its high uniformity allows for its potential use in the removal of inorganic and organic contaminants in aqueous solutions [34]. The functionalization of the silica generates an increase in the adsorption capacity because, upon including functional groups, the active sites increase and the conditions on the chemical surface improve [29].
This work was developed with the objective of preparing and functionalizing SBA-15 with polyaniline to increase the capacity to adsorb phenol from the liquid phase, enabling the study of its interactions through the enthalpy of the immersion of solids in benzene, water, and solutions of 0.1 M NaOH, 0.1 M HCl, and phenol, and the elucidation of the relationships that exist between the enthalpies of the immersion of silica and its modification, along with its physicochemical properties, such as surface area, total acidity and basicity, and pH at the point of charge zero. The phenol adsorption capacity of SBA-15 and SBA-15/PA was determined using the adsorption isotherms for aqueous solutions.

SBA-15 Synthesis
SBA-15 was synthesized according to the method reported by Zhao et al. [35]. Initially, 18 g of pluronic 123 (Sigma-Aldrich, St. Louis, MO, USA) and 135 g of water were weighed and dissolved, and the mixture was left while stirring for 24 h. This solution was then placed in a water bath at 37 • C, a solution of HCl (72 mL of concentrated hydrochloric acid in 450 mL of water) was added slowly, and then 42 mL of TEOS (tetraethylorthosilicate) (Aldrich) was added. Finally, the mixture was left while stirring for 20 h, washed with plenty of water, dried at 70 • C for 24 h, and calcined at 540 • C for 6 h.

Modification with Polyaniline
SBA-15 was modified with polyaniline, as reported by Weng et al. [36]. A solution of 42.9 mL of 0.55 M hydrochloric acid and 1.95 mL of aniline was prepared, followed by the addition of 2 g of SBA-15. This was stirred for 30 min in an ice-water bath at 0 • C, a solution of ammonium persulfate, APS (4.45 g of APS in 19.5 mL of water), was added dropwise, and the mixture was left to stand for 3 h. Subsequently, it was filtered and washed three times with water, acetone, and a 0.55 M hydrochloric acid solution. Finally, it was dried at 70 • C for 24 h. The solid obtained was called SBA-15/PA.

Textural Characterization
The surface area and pore volume were determined by weighing approximately 100 mg of the sample, which had been degassed at 250 • C for 24 h, followed by N 2 fission at 77 K on a Quantachrome Autosorb IQ2 sortometer. The surface area was determined by applying the Brunauer-Emmett-Teller (BET) model. The most common form of the BET equation is: where n is the specific amount of gas adsorbed at the relative pressure P/P o , n m is the adsorption capacity in the monolayer, P is the pressure, P o is the saturation pressure of a substance that is adsorbed at the adsorption temperature, and C is the BET constant that is exponentially related to the monolayer adsorption energy [37]. Micropore volume (V o ) was calculated using Dubinin-Radushkevich (D-R) equation; mesopore volume (V meso ), total volume, and pore diameter were calculated using the density functional theory (DFT) method.

Chemical Characterization
The total acidity and basicity were determined according to the method proposed by Boehm [38] with modifications. First, 0.1 g of the sample was mixed with 50 mL of 0.05 M NaOH and 50 mL 0.05 M HCl. The mixture was left under constant stirring for five days. Subsequently, an aliquot of 10 mL of the solution was titrated with previously standardized acid or base using a CG 840B Schott automatic titrator. The zero-charge point (pH pzc ) was determined using the mass titration method. Quantities between 0.1 and 0.6 g of the samples were weighed, and 10 mL of 0.1 M NaCl was added. They were left under constant stirring for 48 h; then, the pH of each solution was measured with a pH meter, a CG 840B Schott [39].

Determination of Immersion Enthalpy
The immersion enthalpies of SBA-15 and its modification in benzene, water, and solutions of 0.1 M NaOH, 0.1 M HCl, and phenol (100 mg L −1 ) were determined using a heat-conduction microcalorimeter [40][41][42][43][44][45]. For the determinations, 0.1 g of a sample was weighed and placed in a glass vial that was fitted inside the equipment cell. Then, 10 mL of the liquid was placed in a stainless steel cell, the temperature was stabilized, and thermal recording was initiated by an electrical potential signal from the sensors; then, the sample was immersed in the liquid. The resulting thermal changes were recorded, and, finally, electrical calibration was performed.

Adsorption Isotherms
In order to investigate the influence of adsorption conditions on the adsorption ability of mesoporous SBA-15 and SBA-15/PA composites, including the initial concentration of phenol, adsorption temperature and pH value, batch adsorption experiments were conducted previously (these results are not shown here; they were taken as a basis to establish the experimental conditions of this investigation). Firstly, a series of phenol aqueous solutions with concentrations ranging from 20 to 150 mg L −1 were prepared by dissolving different amounts of phenol in deionized water. Then, a series of phenol aqueous solutions with the same phenol concentration (150 mg L −1 ) and various pH values in the range of 2-10 were also prepared by introducing different amounts of HCl or NaOH solutions. Then, for each experiment, the adsorbent was separated from the mixture by centrifugation, and the residual concentration of phenol solution was measured using an ultraviolet-visible (UV-vis) spectrometer (Thermo Nicolet Nexus) at a characteristic adsorption wavelength of 271.0 nm. In addition, for studying the effect of adsorption temperatures on the uptake capacity of the mesoporous composite, 0.06 g of adsorbent was added to 100 mL of 150 mg L −1 phenol solution and stirred at 15, 25, 30, 35, and 40 • C, respectively. Then, the residual phenol concentrations were measured using the same method as mentioned above. Moreover, the adsorption capacity of pristine SBA-15 silica was also analyzed for studying the effect of polyaniline modification on the removal of phenol. Typically, 0.02 g of pristine SBA-15 and SBA-15/PA composite was, respectively, added into 100 mL of phenol solutions with initial concentrations of 20-150 mg L −1 under stirring at 25 • C. After 15 h, the adsorbent was separated by centrifugation, and the residual resorcinol concentrations were measured using the UV-vis spectrometer. Finally, they were left stirring for 58 hours until they reached equilibrium. The samples were filtered and the phenol concentrations were determined on a UV-Vis Thermo Nicolet Nexus spectrophotometer at 271 nm. The equilibrium concentration, C e (mg L −1 ) and the amounts adsorbed on the adsorbent at equilibrium, q e (mg g −1 ), were calculated. The initial pH of phenol solutions is 6.9 in SBA-15 and 3.0 in SBA-15/PA, and these values were maintained during each of the experiments. Figure 1 shows the nitrogen adsorption isotherms obtained for SBA-15 and SBA-15/PA; the samples presented type IV(a) isotherms according to the updated IUPAC classification [39][40][41]. At medium relative pressures, capillary condensation began and H1-type hysteresis loops were generated, which are associated with mesoporous solids with cylindrical pore channels. For SBA-15/PA, the hysteresis loop is H5 type, which is associated with pore structures containing open and partially blocked mesopores. The incorporation of polyaniline in the pores of the silica explains the decrease in the amount of nitrogen adsorbed; this behavior is due to the limitation of the diffusion of nitrogen molecules that generate the surface groups located in the openings of the pores [42]. The textural parameters of the BET surface area, micropore volume (V o ), mesopore volume (V meso ), total volume, and pore diameter were calculated from the experimental data of the N 2 adsorption-desorption isotherm and are shown in Table 1. The surface area value for SBA-15 was 654.5 m 2 g −1 and its pore diameter was 6.08 nm; after the surface modification, the surface area and pore diameter values were reduced to 215.3 m 2 g −1 and 5.28 nm, respectively. Likewise, the total pore volume changed from 0.84 to 0.33 cm 3 g −1 due to the partial blocking of the pore channels of SBA-15 by polyaniline, supporting the idea that polymerization took place primarily within the pores of SBA-15 [42][43][44]. The precision for the surface area values was ±1 m 2 g −1 , and that for the pore volume was ±0.01 cm 3 g −1 , according to the characteristics of the Quantachrome Autosorb IQ2 sortometer.

Chemical Characterization
Using the methodology proposed by Boehm, the acid-base character of the surface of the silica was determined depending on the amount of acidic or basic sites that it contained. Table 2 shows the results obtained for the acidity, basicity, and the point of zero-charge for SBA-15 and SBA-15/PA. SBA-15 has total acid groups of 0.026 meq g −1 and a pH at the point of zero charge of 6.5. According to these values, the acidity of the surface is related to the deprotonation of the silanols; after their modification, there is evidence of a decrease in the amount of acid sites to 0.011 meq g −1 and in the pH pzc to 3.0, these values being attributed to the decrease in silanol groups present on the surface of SBA-15. The pH at the point of zero charge of the two samples is acidic, such that if the solids were put in contact with solutions whose pHs were lower than the pH pzc , the charges on their surfaces would be positive, on average [45,46]. Table 2. Acidity, basicity, and pH at zero charge point of SBA-15 and its modification.

Immersion Enthalpy Measurements
Immersion enthalpy is a thermodynamic parameter that correlates the pore structure and the surface chemistry of the solid. Table 3 shows the values obtained in the calorimetric determinations; all the immersion enthalpies are negative, which indicates that these processes are exothermic. Figures 2-4 show the electric potential curves as a function of time for the immersion process of samples and their modification into benzene, water, NaOH (0.1 M,) HCl (0.1 M), and phenol solution (100 mgL −1 ) to evaluate the interaction between the silica, the solvents, and the solutions used. This information allows the magnitude of the interaction in the adsorption to be established.  Figure 2 shows the potentiograms from the immersion calorimetry in water and benzene of SBA-15 and its modification. The polarity of the surface of the silicas could be evaluated based on the fact that the proportion of hydrophobicity and hydrophilicity of a surface can be estimated using the relationship with the enthalpies of immersion in a polar and a nonpolar liquid, such as water and benzene. As the values of the enthalpy of immersion in water of SBA-15 and SBA-15/PA were greater than those of the enthalpy of immersion in benzene, the surfaces were considered hydrophilic. However, the formation of hydrogen bonds between the water molecules and the silanol groups of SBA-15 generated an enthalpy value of 31.2 J g −1 , and upon reducing the amount of silanol groups in SBA-15/PA by shielding the polyaniline molecules located in the polar sites of the surface of the silica, the enthalpy of immersion in water was reduced to 20.4 J g −1 [39][40][41][42][43][44][45][46][47][48][49].   The absolute values of the enthalpy of SBA-15 and SBA-15/PA were higher for immersion in the sodium hydroxide solution, as observable in Figure 3, confirming a greater interaction with the acidic sites of the silicas. As the enthalpy of immersion of the silicas in the solutions were related to the acidity and total basicity of the former, it is evident that, in the case of the immersion enthalpy being a function of the content of acid groups, there was an increase in the value of the enthalpy of immersion of SBA-15 when the acidity increased, with a value of 0.026 meq g −1 , and an immersion enthalpy in NaOH of 174.1 J g −1 , decreasing after functionalization with polyaniline to values of 0.011 meq g −1 and 167.7 J g −1 .
With respect to the basic groups, it should be noted that they were related to the values for the immersion enthalpy of the silicas in HCl solution, since the values for the basic groups were lower than those for the acid groups. The same was observed with the values of the immersion enthalpy in HCl being lower than the enthalpy in NaOH [41,[49][50][51].
Another characteristic of silicas is the pH at the point of zero charge, pH pcc , which establishes the charge on the surface of the materials and is related to surface chemistry, and which, in turn, affects the energetic interactions that occur between the solid and the immersion liquid. For this reason, it corresponds to the enthalpy of immersion in benzene and a linear behavior is obtained; as the value of the enthalpy of immersion increases, the pH at the point of zero charge also increases, with 28.1 J g −1 being the immersion enthalpy of SBA-15, with a pH pcc of 6.5; for SBA-15/PA, the enthalpy decreases to 15.1 J g −1 , as does the pH pcc to 3.0 [47][48][49].
There is a directly proportional relationship between the enthalpy of immersion and the surface area of the silicas, with the greater 28.1 J g −1 enthalpy of immersion for SBA-15, because when there is a wider space for the entrance of molecules, greater interaction of the adsorbate with the porous structure is generated. The immersion enthalpy value for SBA-15/PA decreased to -15.1 J g −1 . This shows that there was a reduction in porosity due to the fixation of the functional groups on the surface [45][46][47][48][49].
The degree and nature of the interactions at the solid-liquid interface in processes such as adsorption can be determined using the immersion enthalpy, because, for some molecules, it can vary according to the specific interactions between the immersion liquid and the surface of the solid. The immersion enthalpy of SBA-15 in the 100 mg L −1 phenol solution was higher than that of SBA-15/PA. Figure 4 shows that the thermal effect obtained for SBA-15 was high and corresponded to the sum of various interactions, such as the wetting of the solid surface by the solvent, the interaction of the solvents with the chemical groups on the surface, and the interaction of phenol with the surface and its groups. In SBA-15/PA, the access of the molecules to the surface decreased in such a way that there was a considerably low thermal effect; however, it continued to present an exothermic enthalpy, correlated with the processes of the interaction of the aniline monomers and the surface of the solid. The adsorbed amount and immersion enthalpy were affected by the contents of total acid and base groups on the surface of SBA-15 and its modification. It is known that, upon decreasing the content of acid groups, the adsorbed amount of the phenolic compound increased [45][46][47][48][49].

Adsorption Isotherms
In order to evaluate the effectiveness of the adsorbents, the equilibrium phenol adsorption (q e ) was studied as a function of the equilibrium concentration (C e ). Adsorption isotherms are modeled by different adsorption isotherms from aqueous solutions.
Adsorption isotherms from aqueous solutions There are several theoretical models for the description of adsorption isotherms with different considerations, such as the starting hypotheses and the number of characteristic parameters of each of them. The models used in adsorption studies in single-component solid-liquid systems are related to thermodynamic models obtained from theoretical adsorption models based on the thermodynamic properties of the adsorbent and the fluid phase. The others are semi-empirical models, which are the result of an adjustment of experimental data to certain mathematical functions, although with a thermodynamic basis. The models commonly used for the adjustment and interpretation of the experimental data obtained in adsorption studies for solid-liquid single-component systems are detailed below.
Two-parameter isotherm models.

Langmuir Model
This is an empirical model that assumes monolayer adsorption and uniform adsorption energy on the surface [50][51][52]. There is no transmigration of the adsorbate in the plane of the surface [31]. Adsorption only occurs in a finite and fixed number of defined and localized sites that are identically equivalent and in which there are no lateral interactions or steric hindrance between the same adsorbed molecules on the adjacent sites [50]. Adsorption is homogeneous and each molecule has constant enthalpy and activation energy of sorption (all sites have the same affinity for adsorbate) [50][51][52].
where q e = the amount of solute adsorbed per unit weight of adsorbent at equilibrium (mg g −1 ), C e = the equilibrium concentration of solute in solution volume (mg L −1 ), q m = the maximum adsorption capacity (mg g −1 ), and K L = the adsorption equilibrium constant of Langmuir (L mg −1 ).

Freundlich Model
The equation being exponential, the model assumes that the adsorbate concentration increases on the surface of the adsorbent and that it is a multilayer adsorption [53,54]. It describes reversible and non-ideal adsorption with non-uniform distribution of heat of adsorption and affinity on the heterogeneous surface [53]. The amount adsorbed is the sum of the adsorption of all the sites. q e = K F C 1/n e where q e = the amount of solute adsorbed per unit weight of adsorbent at equilibrium (mg g −1 ), C e = the equilibrium concentration of solute in solution volume (mg L −1 ), K F = the Freundlich constant related to the relative adsorption capacity of the adsorbent (mg 1−(1/n) L 1/n g −1 ), and n = the heterogeneity factor, related to the constant that indicates the intensity of the adsorption. If n = 1, it is a linear adsorption. If n > 1, it is a physical adsorption and it is favorable. If n < 1 it is a chemical adsorption. On the other hand, if the value of 1/n < 1 indicates a normal Langmuir-type isotherm and if 1/n > 1, cooperative adsorption is indicated. Two-parameter isotherm models.

Sips Model
The Sips model is a combined form of the Langmuir and Freundlich expressions deduced to predict systems with heterogeneous adsorption and overcome the limitations at high concentrations of the Freundlich equation [55][56][57][58]. At low concentrations, the adsorbate is reduced to the Freundlich isotherm, whereas at high concentrations a monolayer adsorption is predicted, characteristic of the Langmuir isotherm. As a general rule, the parameters of the equation are mainly governed by adsorption conditions, such as pH, temperature, and concentration [55,56].
where q e = the amount of solute adsorbed per unit weight of adsorbent at equilibrium (mg g −1 ), C e = the equilibrium concentration of solute in solution volume (mg L −1 ), q ms = the Sips maximum adsorption capacity (mg g −1 ), K s = the Sips equilibrium constant, and n s = the Sips model exponent.

Tóth Model
This is an empirical equation modified from the Langmuir equation that reduces the error between the experimental data and the predicted values of the equilibrium adsorption data [54][55][56][57][58]. The application of this equation is more appropriate for multilayer adsorption and for heterogeneous adsorption systems that satisfy the low and high concentration ranges [57,58]. Its correlation assumes an asymmetric quasi-Gaussian distribution, with most sites having an adsorption energy lower than the maximum peak or the mean value [57,58].
where q e (mmol g −1 ) and C e (mmol L −1 ) have the same meaning as in the above equations, q mT (mmol g −1 ) = the Tóth maximun adsorption capacity, K T = the Tóth equilibrium constant, and mT = the model exponent. When mT is equal to 1, Toth's equation is reduced to Langmuir's equation, with which it can be deduced that this parameter is related to the heterogeneity of the system.

Redlich-Peterson Model
Similar to the Sips model, this takes into account the characteristics of the Langmuir and Freundlich isotherms in a single equation, indicating that the adsorption mechanism is a combination of both models. Due to its versatility, it allows the representation of an adsorption system in equilibrium for a wide range of concentrations and can be applied for both homogeneous and heterogeneous systems. The model is represented by the following equation [57,58].
where K RP (L g −1 ), RP (L mg −1 ) -β and the exponent β (dimensionless) are Redlich-Peterson parameters. The β value is between 0 and 1 and it determines the degree of heterogeneity of the surface. This equation is transformed into a linear isotherm (Henry's equation) at low surface coverage (β = 0), it is reduced to the Langmuir isotherm when β = 1, and becomes the Freundlich isotherm when K RP y α RP >> 1 y β = 1.
The results shown in Figure 5 for the SBA-15 and SBA-15/PA samples were analyzed and evaluated according to the aforementioned models using the ANOVA analysis of variance optimization method, included in the Sigma-Plot ® v.10.0 software for each adjustment function according to each model. These results correspond to experimental data in the equilibrium data of phenol that were adjusted with the mathematical expressions which represent the corresponding isotherms in their linear form for each model, and their results are presented in Tables 4 and 5.
The associated error analysis shows that the Freundlich model is the most appropriate for representing the adsorption process, as it had excellent correlations of R 2 > 0.99 and low standard mean error percentages, ε < 1.2%. Likewise, the Langmuir model, with R 2 correlations greater than 0.97 and ε errors less than 2.9%, can represent the isotherms in an acceptable manner, as shown when representing the experimental data in Figure 5 by the dotted lines in the experimental adsorption data. The measure of the heterogeneity of the surface 1/n in the Freundlich model in SBA-15/PA is 0.3296; as this is less than 1, it is established that the adsorption process is favorable and indicates that the surface is heterogeneous; the relative adsorption capacity of the SBA-15 modified with polyaniline is 24.07 times higher than the relative adsorption capacity of SBA-15. The maximum adsorption capacity values calculated using the Langmuir model were 138 and 85.54 mg g −1 for SBA-15/PA and SBA-15, respectively, again indicating an increase in adsorption effectiveness upon the inclusion of the polyaniline groups in the structure.   The phenol isotherms in SBA-15 and SBA-15/PA are shown in Figure 5. It was observed that the adsorption order in terms of the amount adsorbed was greater in SBA-15/PA than in SBA-15. This can be explained by the fact that large amounts of polyaniline macromolecules were fixed on the walls of the pores of SBA-15, and since phenol is a weak acid and polyaniline a weak base, acid-base interactions occurred in the adsorption process. SBA-15/PA, having an acidic surface, interacts with the phenol solution, while the amino groups of the aniline interact with the hydroxyl group through hydrogen bonds [49][50][51][52][53][54][55][56][57][58]. When you have the value of n (0 < 1 n < 1) for the Freundlich model, this suggests that the adsorbent surface is heterogeneous. In this study, it is interesting that by modifying the SBA-15 with polyaniline, one went from a surface that was not homogeneous to a more homogeneous one. This had an impact on the adsorption capacity of phenol.
In summary, it can be said that for the two-parameter models, the sample with the best fit between these two models was the SBA-15 sample, and, under the experimental conditions of this research work, there is a greater fit to the Freundlich model compared to Langmuir's, taking from now on the regression coefficient R 2 as a criterion for the analysis of all models.
Regarding the three-parameter models, neither of the two shows a behavior that allows an affirmation that one fits better than one of the abovementioned models, despite having three parameters. If we slightly amplify the R 2 criterion in qualitative way, it can be said that the SBA-15/PA sample has a relative fit to the Sips model. This result is consistent if it is taken into account that the Sips model is derived from a mathematical treatment between the Langmuir and Freundlich models. For this reason, the adsorption of phenol on the SBA-15/PA samples and even that of SBA-15 can be explained by this model. This result turns out to be very interesting and should be further explored at different phenol concentrations and temperatures.
The experimental results with the adsorbents at a constant temperature in equilibrium conditions were tested with various isotherm models (Langmuir, Freundlich, Sips, Tòth, and Redlich-Peterson isotherms) to establish the adsorption mechanisms. The testing of Langmuir and Freundlich isotherms is essential, since it reveals whether the adsorbent surface has a specific number of active sites which come in contact with adsorbates to form a permanent attachment. The test of the Langmuir isotherm ascertains physical adsorption, whereas the test of the Freundlich isotherm determines chemical adsorption, and they also facilitate the determination of adsorption mechanisms. SBA-15 and SBA-15/PA had interesting adsorption capacities and comparable adsorption capacities with respect to other adsorbents; therefore, these adsorbents are suitable for practical use.

Conclusions
SBA-15 was modified using polyaniline on the pore wall according to textural and chemical characterization. Novel adsorbents were synthesized to effectively remove phenol pollutants from aqueous solutions. SBA-15 and SBA-15/PA mesoporous composites were synthesized by successfully modifying polyaniline on the pore wall of SBA-15 silica. The experimental results demonstrate that these composite materials are very promising candidates for the removal of phenol.
The experimental results demonstrated that SBA-15/PA is a promising candidate for the removal of phenol with an adsorption capacity of 138 mg g −1 compared with 85.54 mg g −1 for SBA-15. The surface area of SBA-15 was 654.5 m 2 g −1 , and when it was modified this decreased to 254.6 m 2 g −1 ; the micropore volume, mesopore volume, total volume, and pore diameter presented the same behavior due to the fact that large amounts of polyaniline macromolecules were fixed on the walls of the pores of SBA-15. Increasing behavior of the immersion enthalpy for SBA-15 and SBA-15/PA in benzene was observed with respect to the increase of BET surface area, explained by its same porous nature. According to the enthalpies in water and benzene, the surfaces of the silicas were considered hydrophilic. There was an increase in the enthalpy when the acidity increased. The adsorbed amount and immersion enthalpy were affected in such a way that, as the content of acid groups decreased, the adsorbed amount of the phenolic compound increased.