Abstract
The surface heterogeneity of the multiwalled carbon nanotubes is investigated on the basis of adsorption isotherms from dilute aqueous phenol solutions at various temperatures. The Dubinin–Astakhov isotherm connected with the non-symmetrical distribution function of adsorption energies has been applied to describe adsorption equilibrium. Theoretical isosteric heats of adsorption connected with the Dubinin–Astakhov isotherm have been estimated as well. It is known that the decrease of the isosteric heat of adsorption with adsorbate loading is characteristic for the energetically heterogeneous surfaces. Theoretical description of kinetics is based on the Statistical Rate Theory of Interfacial Transport. The Statistical Rate Theory approach links the rate of transport between two phases with the difference in the chemical potentials of the molecules in these phases. Theoretical studies were verified successfully using the literature experimental adsorption data.
Introduction
Adsorption at solid/solution interfaces plays essential role in a variety of technological processes applied in the protection of the environment. Multiwalled carbon nanotubes (MWCNTs) can be used for effective removal of aromatic organic compounds like phenol, nitrophenols, chlorophenols, etc., from contaminated water (Arasteh et al., 2010; Li et al., 2013; Liao et al., 2008; Pan and Zhang, 2014; Podkościelny et al., 2014; Shen et al., 2009; Tóth et al., 2012; Wiśniewski et al., 2012; Wu et al., 2012; Yang and Xing, 2010; Yang et al., 2016).
Due to strong van der Waals interactions, MWCNTs adhere to each other and form bundles as can be seen by TEM images. The interstitial channel space between the tubes within the bundles can be regarded as pores. Generally adsorption can occur (Shen et al., 2009; Sheng et al., 2010) on the external surfaces of the nanotubes, in the interstitial channels between nanotubes, in the hollow space inside nanotubes, and on the grooves present on the periphery of a nanotube bundle. So, the overall adsorption of phenols contains the contributions from the adsorption on the graphene surface and because of pore filling. Surface defects also play a key role in the properties of the surface (Lehman et al., 2011). Chemical heterogeneity of the nanotubes is associated with the presence of oxygen functional groups, traces of metal catalyst, and amorphous carbon. Above-mentioned heterogeneity of the surface generates the exceptional adsorption properties of MWCNTs.
Sorptive properties of MWCNTs can be modified by chemical oxidation using HNO3, H2O2, KMnO4, etc. After oxidization, the carbon nanotubes became more hydrophilic and have more functional groups on the surface. The increase of pore volume and the BET surface area are observed. It is attributed to removal of amorphous carbon, traces of metal catalyst, and hemispherical caps on the nanotubes (Shen et al., 2009; Sheng et al., 2010).
The main objective of the article is to investigate the surface heterogeneity of MWCNTs when placed in contact with dilute aqueous phenols solutions. The Dubinin–Astakhov (DA) isotherm connected with the non-symmetrical distribution function of adsorption energies has been applied to describe adsorption equilibrium. Theoretical isosteric heats of adsorption connected with the (DA) isotherm have been estimated as well. During adsorption at the solid/solution interface, a solid and a solution are usually brought into contact for a limited period of time. Theoretical description of kinetics is based on the Statistical Rate Theory (SRT) of Interfacial Transport (Azizian and Bashiri, 2008; Podkościelny and Nieszporek, 2011; Rudziński and Płaziński, 2008a, 2008b, Azizian et al. 2008; Rudziński et al., 2005; Ward et al., 1982). The SRT approach links the rate of transport between two phases with the difference in the chemical potentials of the molecules in these phases.
Theoretical studies were tested using the literature experimental adsorption data (Sheng et al., 2010). They included the data of phenol adsorption from aqueous solution on as-prepared and oxidized MWCNTs.
Theory
Adsorption equilibrium
The fundamental expression of the integral equation approach describing single-solute adsorption from dilute aqueous solutions on heterogeneous solids has the following form (Jaroniec and Madey, 1988):
Use of so-called non-symmetrical function in the condensation approximation (CA) and Langmuir isotherm (as the local isotherm) results in well-known DA isotherm in the form (Crittenden et al., 1999; Dobruskin, 1998; Podkościelny and László, 2007; Yang and Xing, 2010; Yang et al., 2008, 2016):
For special values of the parameter r, the (DA) isotherm is reduced to the following isotherms: the Dubinin–Radushkevich (DR) isotherm for r = 2 and the Freundlich (F) for r = 1.
Adsorption model based on the DA equation has been successfully used to describe the adsorption of simple aromatic compounds from aqueous solutions on MWCNTs (Pan and Zhang, 2014; Wu et al., 2012; Yang and Xing, 2010; Yang et al., 2008).
Heat effects of adsorption
The proper analysis regarding the heat effects accompanying the adsorption phenomenon provides additional information on the analyzed adsorption systems (Nieszporek et al., 2009; Podkościelny and Nieszporek, 2011; Podkościelny et al., 2014). The elementary equation which enables to obtain expressions for isosteric heats of adsorption qst has the following form:
The method to receive the theoretical expressions for
Combination of equation (3) and isotherm (4) leads to the expression that defines isosteric heat of adsorption qst:
The parameter
You may notice that the equation specifying the isosteric heat of adsorption (equation (5)), and the adsorption isotherm (equation (2)) both use the same parameters, which is a great advantage of this calculation.
Kinetics of adsorption (SRT theory)
In the publications on the theoretical description of the experimental kinetic isotherms in most cases, the authors use known expressions, for example, Lagergren, pseudo-second-order equation, etc. A relatively new theory is SRT of Interfacial Transport (Ward et al., 1982). This theory was successfully used to describe kinetics of many different experimental adsorption systems in the liquid/solid phase (Azizian and Bashiri, 2008; Azizian et al., 2008; Podkościelny and Nieszporek, 2011; Rudziński and Płaziński, 2008a, 2008b).
The SRT is based on the assumption that the chemical potentials of the adsorbed molecules
The amount of adsorbate in the bulk phase strongly prevails over the adsorbed amount, so after the system is isolated and equilibrated, the adsorbate concentration in the solution, c, does not change much,
The chemical potential of bulk phase
Therefore, the expression specifying rate of adsorption can be obtained by using equations (7) to (9) (Rudziński et al., 2005):
The actual adsorbate concentration, c, can be calculated based on the following equation (Rudziński and Płaziński, 2008b):
The obtained SRT kinetic equation (10) by using equilibrium (DA) isotherm (2) ensures the consistency of theoretical studies.
Results and discussion
The basic characteristics of the MWCNTs surface (Sheng et al., 2010).
Equilibrium adsorption
The adsorption isotherms of phenol from aqueous solutions on as-prepared (CNT) and oxidized MWCNTs at 298 K, 318 K, and 338 K are presented in Figure 1(a) and (b), respectively. The symbols denote the experimental data (Sheng et al., 2010), whereas the lines are theoretical isotherms estimated from DA equation (2). It can be seen that the adjustments of DA to the experimental isotherms are very good. The results of calculations for equation (2) are summarized in Table 2. The table includes the values M of monolayer capacity, E/kT values, and the heterogeneity parameter r. The quality of the fit of equation (2) to the experimental data is described by the following Error function (residual sum of squares): The adsorption isotherms of phenol from aqueous solutions on as-prepared (a) and oxidized (b) MWCNTs at various temperatures. The symbols are the measured values of isotherms (Sheng et al., 2010) and the lines are the theoretical isotherms calculated from (DA) equation. Parameters characterizing adsorption of phenol from aqueous solutions on as-prepared (CNT) and oxidized (CNTox) MWCNTs (Sheng et al., 2010) obtained from the (DA) equation (2).
The adsorption mechanism of phenol is complex. The driving forces are “π-π” electron-donor-acceptor interactions between aromatic molecules and the polarizable graphene sheets of MWCNTs (Chen et al., 2008; Pan and Xing, 2008; Podkościelny et al., 2014; Shen et al., 2009; Sheng et al., 2010). Both electron-donator and electron-acceptor groups on benzene ring can increase the adsorption on MWCNTs. The –OH group as a strong electron-donating group makes the benzene ring(s) electron rich, thus allowing the compound to interact strongly with the (electron-depleted) surfaces of carbon nanotubes.
After oxidization, hydrophilic oxygen-containing groups, including carboxylic, lactonic, and hydroxyl groups were introduced into the surface of MWCNTs. The addition of these groups results in a more negatively charged MWCNT surface due to deprotonation of carboxylic groups at the equilibrium adsorption (near neutral pH = 6.5). At this pH, adsorption of water is more favorable relative to the adsorption of simple aromatic compounds. Therefore, water is preferentially sorbed by the oxygen-containing surface groups (Tóth et al., 2012). The adsorbed water molecules can act as nuclei for the formation of larger H-bonded water clusters. As a result, the associated water can prevent migration of the phenol molecules to the hydrophobic parts of the surface and effectively reduce their uptake. Carboxylic groups on the surfaces of oxidized MWCNTs can act as electron withdrawing groups localizing electron from π system of MWCNTs, and in consequence—weakening the π–π interaction between aromatic ring of phenol and carbon nanotubes. Many studies show that oxygen functional groups decrease the adsorption of organic compounds on carbon materials such as activated carbons, carbon fibers, and MWCNTs.
Calorimetry
The next stage of our studies is analyzing heat effects accompanying adsorption of phenol on MWCNTs (Nieszporek et al., 2009; Podkościelny and Nieszporek, 2011; Podkościelny et al., 2014). Unfortunately, the experimental heats of adsorption were inaccessible for the analyzed adsorption systems, so we could calculate the theoretical ones only. From the theoretical analysis of the experimental equilibrium adsorption isotherms by DA equation (2), we obtained the quantities E/kT and r. These parameters can just be used to calculate the theoretical isosteric heats of adsorption qst in terms of equation (5). Thus the occurrence of common parameters is the significant advantage of both equations.
Figure 2(a) and (b) presents the theoretical isosteric heats of adsorption of phenol from aqueous solutions estimated based on equation (5). The calculations were carried out by using the parameters included in Table 2.
Theoretical isosteric heats of phenol adsorption from aqueous solutions on as-prepared (a) and oxidized (b) MWCNTs, calculated based on equation (5). The calculations were performed by using the parameters included in Table 2.
In all analyzed systems, the isosteric heats of adsorption decrease with the increasing surface coverage. The sharp decrease of the isosteric heat at low surface coverage followed by a slower drop at higher coverage is observed. It is known, the decrease of the isosteric heat of adsorption with adsorbate loading is characteristic for the energetically heterogeneous surfaces. The correctness of adjustment of (DA) equation (2) to the experimental adsorption isotherms is confirmed by the regular temperature dependence of the theoretical isosteric heats of adsorption for a given adsorption system (Figure 2(a) and (b)).
Kinetics of adsorption
Important issue of our research is theoretical description of adsorption kinetics. Figure 3 presents the comparisons of the experimental kinetic isotherms for the systems studied (Sheng et al., 2010) with the theoretical isotherms (lines) calculated from the SRT equation (10) at 298 K.
Comparison of the experimental kinetic isotherms for phenol adsorption from aqueous solutions on as-prepared and oxidized MWCNTs (Sheng et al., 2010) with the theoretical isotherms (lines) calculated from the SRT equation (10). The values of the obtained best-fit parameters are included in Table 3.
Values of parameters used while fitting kinetic experimental data (Sheng et al., 2010) presented in Figure 3 by SRT equation (10).
The difference in adsorption kinetics between as-prepared and oxidized MWCNTs may be due to the changes in the surface chemistry of MWCNTs after HNO3 oxidization, which was discussed in “Equilibrium adsorption” section.
The agreement between the theoretical kinetic isotherms and the experimental ones is good for all analyzed systems, which emphasizes the utility of the SRT model for the description of kinetics of phenols adsorption on MWCNTs surfaces. The result of the lower adsorption for HNO3-oxidized MWCNTs indicates that surface chemical properties, rather than specific surface areas or pore volume, are important factors to determine the adsorption of MWCNTs.
Originally, the kinetics for systems studied was described in terms of pseudo-second-order (PSO) model. It should be noted that a theoretical interpretation for (PSO) model based on SRT has been derived (Rudziński and Płaziński, 2006), which indicates the universal character of this theory.
Conclusions
The concise theoretical description of the phenol adsorption from aqueous solutions on MWCNTs has been presented. The (DA) isotherm equation has been applied to investigate the effect of the surface heterogeneity. The perfect adjustments of (DA) isotherm to the experimental data have been obtained. Then, theoretical isosteric heats of adsorption connected with the (DA) model have been estimated. For the all analyzed systems, the isosteric heats decrease with the increasing surface coverage. Correctness of fitting of (DA) isotherm to the experimental data is confirmed by the regular temperature dependence of the theoretical isosteric heats of adsorption for a given adsorption system.
The final step in our research was theoretical description of adsorption kinetics based on the SRT of Interfacial Transport. The chemical potential of the adsorbed phase μs was determined in terms of the DA isotherm. The agreement between the theoretical kinetic isotherms and experimental ones was good for all analyzed systems. Therefore, the utility of the SRT model for the description of kinetics of phenols adsorption on the MWCNTs surfaces has been confirmed. We can conclude that greatest advantage of the model of calculations is the set of common parameters occurring in each type of equations concerning adsorption equilibrium, heat effects, and kinetics of adsorption.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
