1. Introduction
Since it was first identified (Scheele, 1777), gas adsorption process had been investigated for more than two hundred years. Among the large number of phenomena nowadays recognized involved in adsorption, the attention of the scientific community was focused on two main issues. First, from the applied point of view, characterization of solid surfaces through the observed behavior of the adsorbed phase was pursued. Second, and from the more basic science perspective, elucidation of the nature and magnitude of the interaction forces for different experimental conditions and systems.
Taking into account that the behavior of adsorbed molecules depends on the properties of both the solid surface and the adsorbate itself, and for somehow exploring basic aspect of the problem, we focused our attention only in a set of basic carbonaceous substrates and a small groups of simple adsorbates.
We have organized this review in several sections, starting in Section 2 with a brief discussion of some carbonaceous surfaces associated to graphitic structures when they are considered at different scales. This section also deals with different carbon materials and their properties as well as how description of several basic structural units and orientation degrees emerge (Bandosz et al., 2003). In Section 3.1 the adsorption of different gases such as nitrogen, carbon dioxide, argon and methane are described, with special emphasis in the last one due to its environmental relevance (Beaver & Sircar, 2010; Wagner, 1996; Wuebbles, 2000). Three approach to study gas adsorption on graphite surfaces that tackle the problem from different points of view where discussed ( Albesa & Vicente, 2008 ): i) semi empirical quantum mechanical (SQM) model for gas – graphite interaction; ii) Grand Canonical Monte Carlo (GCMC) simulations (D. Do & H. Do, 2005; Sabzyan & Babajani, 2005); and Mean–Field Approximation (MFA) of the lattice model. To emulate by GCMC simulations the interaction gas – surface several models are compared ( Albesa et al., 2008 ), and then a unified description is utilized for the adsorption mechanism that therefore was not only merely based on the three above three mentioned theoretical tools, but also on experimental data. The results reveal changes occurred the adsorbed species, and also can be used to evaluate interaction energy between the gas and the graphite surface as the coverage increases; these allow elucidating somehow the role played by the substrate in adsorption process at different pressures and temperatures (Llanos et al., 2003; Steele, 1974). Section 3.2 is devoted to study the adsorption of simple gases using graphitic curved surfaces, for this end we used abs initio density functional theory (DFT) calculations for graphene sheets whose surfaces had different curvature, and discuss the relation between these model systems and experimental data for adsorption on single wall carbon nanotubes (Albesa et al., 2009; Albesa et al. 2010). In section 4, Gas Separation, we present simulated adsorption isotherms obtained for two characteristic carbonaceous structures, graphite and single wall carbon nanotubes.
2. Relation between different carbon and graphitic surfaces
The great interest in graphitic surfaces stems from their high surface activity, associated with strong adsorbate–surface forces and large surface area. From the commercial point of view, such structures are generally cheap to produce, and can be prepared with different special characteristics. As a result of the strong covalent sp2 bonding in graphite, the distance C – C in graphite (
The need to improve this activity, in order to have greater adsorption capacity and selectivity in mixture separation, has motivated the development of new carbon structures, related direct or indirectly to the graphitic ones. In few years the studies that began changing simple graphitic surfaces, rapidly grow beyond the broad word of fullerenes (Harris, 2003) and reach new and more exotic structures like nanotubes and nanohorns (Dresselhaus et al., 1996).
The two most common allotropes of carbon that occur naturally are graphite and diamond. Graphite is assembled from parallel sheets of hexagonal carbon atoms arrays, each atom being linked to three other sp2 hybridized carbon atoms by sigma bonds or length
The Can fullerenes are shell–like molecules consisting of n three–coordinated carbon atoms arranged as 12 pentagons and ½ (n – 20) hexagons where n is even, and greater than 20 except 22. In spite of great a number of these molecules have been detected experimentally, only a few high–purity fullerenes are produced in macroscopic quantities. Fullerene molecules are essentially polyhedral rather than spherical and, increasing the number of carbon atoms, each fullerene can exist as one of several different isomers.
Just like in order to obtain a fullerene one has to lose the original sheet constituted by hexagonal carbon rings (incorporating pentagonal rings to the atoms arrange), the single walled nanotube (SWNT) is other kind of structure where this not happen and graphene configuration is preserved. In other words, the unit cell of a SWNT is constructed from a portion of graphene sheet defined by the chiral vector
Like in the case of fullerenes, the graphenic membranes can take a conical shape only adopting arranges that are not purely formed by hexagonal carbon rings. The apex has to be defined by one or more ring atom removed in order to reach a conical hat that lead to the form of fivefold (or smaller) rings (Yudasaka et al., 2008). One class of such conical structures with a particularly sharp apical angle are the single–wall carbon nanohorns (SWNH). Conical graphenic structures with wider opening angles (corresponding to fewer pentagons at the apex) sometimes form multilayer structures). The diameter of an individual SWNH ranges from 2 to 4 nm, and the length is 40 to 50 nm. About 2000 of them assemble to form a spherical aggregate with a diameter of about 100 nm. Three types of spherical aggregates are known and, are called dahlia, bud and seed because of their appearances. The dahlia aggregate has long cone – shaped tips sticking through its surface, while the bud aggregate does not. The seed aggregates have lower graphitization that the dahlia and bud aggregates, and their tubules are corrugates (Azami et al., 2008).
In so far as increase the complexity of the carbonaceous structure the long range order of the graphene sheet is lost, because of the curvature, the presence of apexes, etc. Considering, according to the vast amount of experimental data, the hierarchical structure of the activated carbons (Bandosz et al., 2003), they would constitute perhaps a greater example of the previous mentioned loss. The first level of the hierarchy are the so-called basic structural units (BSU) that consist of a few roughly aligned polyaromatic–like molecules or “layers” that generally not preserve the graphite order between them, and due to the presence of functional groups, interlayer spacing are generally greater than that of graphite. The BSUs are assembled to form regions of local molecular orientation (the second level of the hierarchy) which are in turn assembled in space to yield the complex structure of activated carbons.
3. Adsorption on carbon surfaces
3.1. Simple gas adsorption on open graphite surfaces
The literature offers numerous experimental and theoretical studies of the gas adsorption on open graphite surfaces, considering the case of spherical molecules, e.g., rare gases (Cheng & Steele, 1990; Nicholson & Parsonage, 1982), simple molecules, e.g., nitrogen, oxygen, carbon dioxide and, methane (Bottani & Bakaev, 1994; Kowalczyk et al., 2005), and more complexes like ethylene or phenol (Bertoncini et al., 2000; Bottani, 1999). However, from a technological point of view, most of these studies restrict their descriptions to special conditions, i.e., low pressures or temperatures, soft surfaces, weak interactions, etc. In spite the powerful theoretical tools developed in the last decades, from semi empirical and DFT calculations to canonical (CMC) and grand canonical (GCMC) Monte Carlo simulations, nowadays two significant constrains limit computer performances to deal with the adsorption on graphite. One is the number of molecules considered in the calculations, and this frontier moves as fast as new computer facilities are developed. The other problem comes from the approximations adopted in each model. The structure of different adsorbates on the basal plane of graphite has been extensively studied using Monte Carlo and molecular dynamics, however, the validity of these results depend on the interaction potentials adopted. The interaction potentials used in almost all computer simulations are composed of two parts: fluid–fluid and fluid–surface interactions. Some years ago the proposal (Steele, 1974) of an approximate interaction potential between a spherical molecule and a graphite surface (averaged over carbon atoms) reduced and became accessible the expensive computation time. On the other hand many systems and/or conditions require the knowledge of the real potential beyond how well the adsorbate–adsorbate interactions are described. As in example, the variation of the entropy of adsorption with the coverage in the basal plane of graphite determines the shape of the adsorption isotherm (variation of the chemical potential), because in the submonolayer region of this homogeneous surface the adsorption energy varies slightly due to the adsorbate–adsorbate interaction and the main contribution to the change on the chemical potential comes from the entropy of adsorption. Different corrections were proposed to improve the approximate gas-surface interaction potential, considering two dimensional Fourier series expansions (Kim & Steele, 1992) or simply introducing adjusting prefactors.32 In spite it was argue, according to different reasons (Bottani & Bakaev, 1994; Kowalczyk et al., 2005), that neglecting the energetic inhomogeneity of the surface along the graphite plane is not expected to affect the results significantly. However, considering all kind of results that are possible to be obtained from homogeneous surface potentials, one can hardly agree with this affirmation, at least at low temperatures. It’s because energetic homogeneity assumption means all sites on the surface are equivalent then, during the simulations, the probability that a molecule remains in one site or another is the same. On the contrary, considering energetic inhomogeneity, a molecule adsorbed on a site of the surface can reach other more favorable site later, during the simulation steps, in order to reach equilibrium that means inhomogeneity requires more trials to obtain each simulated point, even at low coverage. Of course, while the differences can be considered negligible at very low coverage, as the pressure increase they become significant, i.e. all sites are not equivalent on the surface and there will be fewer molecules adsorbed in the monolayer, and so on.
Among the simple gases that can adsorb on graphite, now we concentrate our study on the case of methane. Different reasons have leading the interest in the adsorption of methane on carbonaceous surfaces in general and, on graphite in particular. Practical reasons come from the fact that natural gas, which is composed mostly of methane, provides an alternative to traditional liquid petroleum fuels because their environmental advantages and natural abundance. From the theoretical point of view, next to the numerous and complete descriptions about rare gases adsorption process, the interaction between methane molecules and graphite surfaces remain with many unknown subjects to study.
According to the previous arguments, all description of the adsorption process, even on a simple graphite surface, requires a previous microscopic knowledge of the interactions between gas molecules and the substrate and between adsorbed molecules, as well as the different configurations of these molecules on the surface. But describing the interaction potential between methane molecules and graphene surface by using classical quantum mechanical methods, for real systems (i.e. many molecules involved), becomes a task beyond regular computational costs. On the other hand we are interested more than absolute values of the adsorption energies on each site of the surface, about the relative differences between one site and another and, most of all how the presence of neighboring adsorbed molecules modified the adsorbed system geometrically and energetically. To tackle these problems we performed PM3 semi–empirical and density functional theory (DFT) calculations of increasing numbers of adsorbed methane molecules on a graphite surface. To mimic the graphene plane we have chosen a planar polyaromatic molecule which has chemical characteristic similar to graphite sheet in many aspects, as is shown in Fig. 1. The interaction energy,
where
The interaction energy between methane and the graphite surface, at low coverage, is almost the same from one adsorption site to another. Although it seems to indicate that graphene structure does not play an important role in methane adsorption, one has to take into account that this result tell us only what happened at very low coverage, but says nothing about the consequences of neighboring adsorbed molecules, in other words, when the coverage increases how near the first adsorbed molecule next molecules can be placed?.
The interaction energy, corresponding to the most favorable configuration, between one, two, and three methane molecules adsorbed on a graphene surface, is almost the same (
By increasing the coverage (i.e. four methane molecules), the structure is modified, the distance between methane molecules changes to 3.95 Å and between methane and graphene 3.80 Å, the interaction energy becomes
Computer Monte Carlo simulations (MCS) with the grand canonical ensemble (
Albesa et al., 2008
) is one of the more useful techniques that gives microscopic descriptions without expensive computational cost and allowing to over comes the quantum mechanical limitation of have to consider only a few adsorbed molecules. MCS offers not only the possibility of studying the adsorption phenomena at molecular scales but also making direct comparisons with experimental results. From the theoretical point of view MCS allow to compare how homogenous or non homogenous surface potential can be affect the results. We performed this comparison for the basal plane of graphite by taking the fluid–surface interaction (
where
being
At low pressures there is a good agreement between experimental and simulated isotherms for both models. In this pressure ranges, when the temperature changes, the curves exhibit different shapes. At temperatures lower than 113 K isotherms shows two clear horizontal steps located at relative pressures of 13/40 and 29/40, respectively (see Fig. 2a), due to the completion of the first and second layers, that confirms an adsorption mechanism layer by layer (Hamilton & Goodstein, 1983). As the temperature increases, this kind of ordered adsorption is less significant, the steps are absent (see Fig. 2b) and all layers are available to be filled.
In Figure 3 the density profile
Among the thermodynamic aspects of gas adsorption on carbonaceous surfaces, enthalpy of adsorption study gives interesting and useful information about the process being accessible to experimental and theoretical explorations. From the experimental side both: direct methods like adsorption (Gravelle, 1978), flow (Groszek, 1998), differential scanning (Baudu et al., 1993), and immersion calorimetry (Zettlemoyer & Narayan, 1966); and indirect methods as inverse gas chromatography and adsorption isotherm have been widely employed last years. When the adsorbed phase is a rare gas, nitrogen, carbon dioxide or, methane on simple surfaces (graphene, single walled nanotubes, or nanocones), the use of adsorption isotherm measurements to determine the variation of adsorbed amount with temperature gives enough results to obtain accurate and reliable values of the heat of adsorption
From two isotherms determined at similar but different temperatures,
where
From the theoretical side the Computer Monte Carlo simulations (MCS) with the grand canonical ensemble not only gives the average number of molecules in the simulation box but also yields information about the isosteric heat, discriminating the two contributions, i.e. the solid–fluid interaction, and the fluid–fluid interaction. The isosteric heat of adsorption,
The square–bracketed term in eq. 5 is the contribution of the fluid–fluid interaction to the isosteric heat of adsorption, whereas the last term is the contribution from the fluid–solid interaction. The isosteric heat of adsorption
To describe the adsorbed system, aside quantum mechanical calculations, that gives answers at atomic level, where only a few molecules are analyzed; and Monte Carlo simulations (MCS) with the grand canonical ensemble, which explore at molecular level, where a more representative number of molecules is taken into account; the mean–field approximation (MFA) of the lattice model (De Oliveira & Griffiths, 1978) is other theoretical tool that that tell, at mesoscopic level, the behavior of adsorption isotherms of several layer of rare or simple gases adsorbed on graphite. According to the differences in the molecular density of the first layer obtain with the two mentioned models employed in MCS, the MFA could elucidate about this controversy. The mesoscopic level, where MFA operates, gives a global description of both the surface and the gas, but ignores what happens on the molecular level, and also assumes, a priori, a layer by layer adsorption process, being each layer constitutes by molecules adsorbed on a fixed lattice. These three mentioned techniques, in their respective ranges, give complementary information of the same phenomena.
Following the description of methane – graphite system, in the MFA one assumes, as in the lattice gas model (Burley, 1972), that the region above the substrate accessible to methane molecules can be divided into cells whose centers forma a regular lattice. No more that one molecule is permitted in a cell, and each par of gas molecules in adjacent cells contributes an amount
The occupation number,
For each configuration {
β =
Calling
is satisfied for
To compute adsorption isotherms, eq. (9) is truncated at
The equation can be solved numerically for different values of
Figures 5 a and 5b show the results obtained taking
3.2. Simple gas adsorption on graphitic curved surfaces
Carbon nanotubes and in general microporous carbonaceous materials offer interesting properties from the technological point of view, as a mean to store hydrogen and methane (Bhatia & Myers, 2006), that are energy sources environmentally favorable, or as sieves to gas separations like oxygen and nitrogen (Arora & Sandler). Adsorptive properties of nanotubes and fullerenes differ from other graphitic carbons because of their carbon surface curvature (C–C–C bonding angles). Carbon atoms have to adopt quasi sp2 hybridization due to their highly curved structure. The level of curvature of valence orbital depends on the material radius: great radius lead to hybridization near sp2 pure, while small radius leads to sp3 (Niyogi, 2002).
To describe phenomena that occur at molecular level one employs simulations, Monte Carlo grand canonical ( Albesa et al., 2008 ) to study adsorption at equilibrium, or Molecular Dynamics (Cheng et al., 2005) when one is interested transport properties far from the equilibrium. However, the accuracy of these descriptions depend on the potential models adopted, choosing generally the 12–6 Lennard Jones potential, or the Crowell – Brown potential that don’t take into account the surface curvature (Wang & Johnson, 1999). A correct potential has to identify the equilibrium position and to give a good approximation in its neighborhood, at least up to the curve inflexion point. To reach this it is necessary to have with a well description of the molecular bonds, because the approximations would neglect the effect of the hybridization differences of the carbon atomic orbitals, induced by the curvature of the surface. To solve this problem force fields were derived (Kostov et al., 2002) where the curvature effect was taking into account, however, the parameters calculated for these fields overestimate the interaction of the nanotubes of smaller radius, because they are considered as a model molecule with free radicals.
In order to describe the potential energy that suffer a molecule near the surface of different nanotubes, theoretical calculations at DFT level were performed (Albesa et al., 2009) for hydrogen, nitrogen and methane, were the surface was mimicked by deforming a coronene molecule (see Fig 6). The nanotube diameters
The exploration involves diameters according to experimental values, that means,
The experimental isotherms for the adsorption of methane, measured at 77.3 K, unlike the case of the adsorption on graphite that has only one step in the monolayer regime (0.2 kPa), have one sub step below this pressure and other one above it (Albesa et al., 2010). As in the case of graphite when the temperature increases, this stepwise adsorption behavior is less pronounced, and above 103 K the substeps disappear.
To explain this behavior, beyond the theoretical interest to explore the adsorption of simple molecules on the curved surface, that mimic the external environment of a nanotube, these kind of carbonaceous structures don’t appear alone but forming different bundles of numerous tubes. These kind of arrange generate new potential adsorption sites, besides the mentioned sites on curved graphene, corresponding to sites on the external phase of the bundle external tubes. For single walled nanotubes (SWNT) four different adsorption sites has been identified (see Fig. 8), i.e. the previous mentioned convex external walls or outer surface sites (S), the interstitial channels (ICs), the grooves between two adjacent outer nanotubes of the bundle (G), and, in case of SWNTs open at the ends, the inner sites (T).
To describe the adsorption process, in closed SWNT, it has been proposed (Bienfait et al., 2004; Kuznetsova et al., 2000) that it starts as linear chains at the strongest binding energy sites: grooves on the outside surface of the bundles and some larger, accessible, defect – induced interstitial channels. After these sites are filled, adsorption proceeds on the external surface of the bundles. The 2D adsorbate structure on the external bundle surface initially builds up adjacent to the occupied grooves until the entire external surface is covered by a single monolayer. The binding energies for this latter stage of adsorption are comparable but somewhat smaller than for adsorption on the basal plane of graphite, whereas the binding energy on the preferential adsorption sites (G, IC) populated during the initial stage of adsorption is considerably larger that tat on the planar graphite. In order to elucidate between the more attractive sites like groove (G) and interstitial channels (ICs), which one begins to fill up first, it has been performed Ab Initio and Molecular Mechanical Calculations, for the case of methane. In order to mimic different sites (S, IC, G, and T) of the substrate, a triangular array consisting of three identical tubes was used, this is the minimum bundle configuration that exhibits the four (S, IC, G) characteristic sites (see Fig. 8). The distance between the tubes or van der Waals gap always is taken equal to
A complementary picture at molecular level, can be performed by Monte Carlo simulations (MCS) with the grand canonical ensemble, in order to take into account the behavior of a representative number of molecules in the adsorbed phase. Results showed in Figure 9, for tubes with a diameter of 2.034 nm, confirm the sequence of filling of different sites in the bundles (is spite of not being realistic, this diameter allows adsorption in IC sites, and the other sites follow the same steps, independent of the kind of nanotubes employed).
The isotherms of nanotubes of 1.628 nm of diameter at different temperatures (83, 93, and 113 K) are compared between closed and opened end in Figures 10 and 11. Unlike the closed end (Fig. 10), the opened end tube (Fig. 11) shows a knee when a monolayer starts forming on the inside surface of the nanotubes, and also when this interior layer is completed (Tasca et al., 2002), and a new high energy site appears in the axial phase (see Fig. 12). It can also mentioned the isotherms corresponding to smaller tubes (0.949 nm of diameter) exhibits only the first knee because there is no possibility of an inner axial phase.
This dependence, of the adsorption process, on the nanotube diameter, suggests that it would be useful to classify the bundles according to the diameter of the tubes involved, analogous to the method used for porous substrates classification employing porous sizes.
Comparing simulated isotherms with the adsorption on graphite we note that the monolayer on graphite is completed at a pressure below that corresponding to the nanotubes, and that the isosteric heat of adsorption is greater for graphite. This is due to the greater gas – solid interaction for the graphite case, because the nanotubes are single – walled, while graphite is formed by several graphene layers. The gas – gas interaction on the nanotubes is weaker than on graphite, because the curvature of the nanotubes reduces the number of neighboring adsorbate molecules in this case, relative to those present in the graphite case.
4. Gas separation of nitrogen-oxygen and ethane-ethylene mixtures
Most commercial processes used in air separation employed synthetic zeolites. These materials are selective with respect to nitrogen, that is, they are able to adsorb more nitrogen than oxygen in a rate 4:1, due to the interaction between nitrogen molecule and oxygen cation to compensate charges. However find some kind of adsorbent that prefer oxygen is also interesting, because to reach the same separation than zeolites will only require a quarter of work (Jayaraman & Yang, 2005). Among different substrates that have greater affinity with oxygen, carbonaceous materials are promising options (Arora & Sandler, 2007), being their principal use as selective sieves. Beside the isotherms, and the isosteric heat of adsorption, selectivity is other quantity of interest, selectivity of specie 1 respect to specie 2 is defined as:
where θi, and pi, (i = 1, 2) are the relative coverage and partial pressure respectively, and, when there is no interaction between adsorbed species, it only depends on temperature, T, otherwise it is function of θ1, θ2 and T. In previous sections it was mention that adsorption explorations at molecular levels, i.e. by MCS, imply to select a model, and in case of gas mixtures, choose one approach enough simply to reduce computing time, but preserving the essential characteristics, becomes critical. One possibility is consider the molecules as pseudo spheres, that is only one Lennard – Jones interaction site per molecule (1CLJ), and a more realistic one would be consider molecules as composed by two interaction sites (2CLJ). We make a comparison of these two models for the adsorption of oxygen and nitrogen mix similar to the air composition, i.e. 80% and 20 % respectively. A remarkable difference appears in selectivity (Fig.13) because 1CLJ approximation gives fewer adsorption of nitrogen than oxygen, although the difference between adsorbed quantities is not high, reflects an important increase of the 1CLJ model selectivity. In Figure 14 we also note that the adsorption for 2CLJ model begins a lesser pressures than in the case of 1CLJ model. The density profiles (Figs. 15 and 16) show that the increase of pressure in 1CLJ model moves some molecules from the first layer to the second layer, so for nitrogen as for oxygen. In 2CLJ approximation this moving occurs basically with oxygen, while the density of the nitrogen first layer increase, this behavior, opposite to what one would hope from the entropic point of view according to the molecules size, is due to temperature considered, because greater temperatures favor nitrogen adsorption
Given that carbon nanotubes have internal sites that, after activation, are available for adsorption, and these have a narrow pore size distribution, these materials are presented as a promising alternative as adsorbents capable of achieving a good separation between N2 and O2.
By comparing the results obtained by Monte Carlo simulations in nanotubes (Fig. 17 a-c) with those obtained in graphite is observed that both materials show similar trends: there is a maximum degree of separation every time a monolayer is completed. As the nanotubes can complete an inner and outer monolayer, there are two peaks in the separation. Nanotubes with greater separation power are the
The degree of separation on the outer surface is negligible. From this arises the importance of obtaining good methods of activation, since the presence of interstitial defects (which are the size of the order of the nanotubes
Olefin – paraffin separation represents one of the most important separations in chemical and petrochemical industries. The production of plastics, rubbers, films, and other chemicals from olefins such as ethylene, requires high purity olefins (> 99%). Ethylene is generally produced by cracking or thermal decomposition of ethane, and during these processes the gas is an ethane and ethylene mix. The conventional cryogenic distillation is as efficient and reliable, and continues being the leading technology for ethane/ethylene separation; however, it is also energetic expensive due to the similar volatility of ethane and ethylene. In a typical ethylene plant, the cracking apparatus represent approximately 25% of the unit const, while the other 75% is due to compression, calefaction, dehydration, recovering and refrigeration systems (Anson et al., 2008).
An effective separation method that produces raw materials with highly ethylene enrichment without using cryodistillation would reduce the energetic and equipment costs associated with paraffin/olefin separation. Therefore exist interest to develop new separation techniques (Fuentes & Mendez, 2002). Many alternative separations have been investigated, including chemical, extractive distillation, physical adsorption, and separation of the base membrane (Eldridge, 1993). Adsorption is one alternative to separation, which is till under develop. Adsorbent materials are generally incorporated joined to transition metals (copper or silver), and as a result the preferential adsorption of olefin is obtained. This selective adsorption is due to the strong interaction between non saturated olefins bonds and the metal ion on the surface, forming a π complex. Because of the preferential way in which olefin is adsorbed (Anson et al., 2010) by these substance it is possible to obtain paraffin pure. However, when paraffin selective adsorbed, olefin pure can be obtained, which becomes adsorption as efficient alternative. Unfortunately there are few materials able to reach this separation (Herdes et al., 2007). This preferential adsorption is due fundamentally to the interaction between methyl sites of the adsorbent and the adsorbat (Kroon & Vega, 2009).
It is well known that for only one component, the adsorption on carbonaceous materials is greater for ethane that ethylene. That is why these materials are studying when one is looking for a highly pure ethylene.
The simulation results show that the selectivity (Fig. 18), as a function of the number of adsorbed molecules, reaches a maximum when the monolayer is completed, as it happen with nitrogen – oxygen mixtures.
The profile is the same for all the compositions studied. The decrease in the degree of separation occurs because, once filled the monolayer, the adsorption of ethylene occurs more rapidly than ethane (Fig. 19 a-c). At high concentrations of ethylene, the number of molecules of ethane decreases, i.e., ethylene displaces ethane (Fig. 19 a). Here we see two phenomena, the principle of selective adsorption of ethane is due to energetic effects. As pressure increases the entropic effects are greater and that is why, because of their size, increases the adsorption of ethylene
Isosteric heat of adsorption of the mix is quite similar to the corresponding values of pure substances. It is also observed that there isn’t a great variation of isosteric heat of adsorption with composition.(Fig 20 a-b)
A similar situation is observed in the adsorption of mixtures of nanotubes
In this substrate the displacement of ethane molecules by ethylene which is observed in graphite is not present. The effect can be attributed to the curvature of this kind of nanotubes; this is, the small diameter provokes a less compact arrangement of the molecules, and therefore a less important effect of the entropic factors.
If we now consider nanotubes with larger diameters, so that adsorption occurs in the interior of them, we note that, as in the case of mixtures of nitrogen and oxygen, the
separation on the outside is very small, being slightly greater than one. Also note that the degree of separation in nanotubes
Of particular interest is that the results obtained in systems nitrogen / oxygen and ethane / ethylene for the case of graphite, show similar trends to those obtained from more complex substrates, such as carbon nanotubes. This is why the above mentioned studies provide rich information on new systems, because, thanks to its structural simplicity, experiments can be carried out with great precision.
5. Conclusions
Experimental results of methane adsorption at low temperature, confirm a layer–by–layer filling mechanism, where the first two layers are clearly defined, but above 113 K this mechanism disappears and all layers are available to be filled. MCS confirm the necessity to employ atomist models in order to describe the adsorption a higher pressures, beyond that at very low pressures simple analytic models give good results with cheaper computational cost. The production of changes on the arrangement of adsorbed molecules as surface coverage increases is also corroborated, at microscopic level, by quantum mechanical descriptions. Although the heat of adsorption constituted an experimental and theoretical useful tool to explore thermodynamic aspects of the adsorption process, and from can be easily determined from the experimental side, theoretical descriptions have to be considered with some details. For instance if the models deal with isotherms associated to ideal monolayer, with molecules laying in a plane at a fix distance of the surface, heat of adsorptions values obtained are well defined, but by using other models, perhaps more realistic, when molecules are not arranged on a plane but around a mean distance of the surface, the number of molecules needed to have values with low dispersion is higher than the usual employed in MCS. The comparison, between experimental results and those obtained from MCS with different model approximations, confirms previous asseveration. In addition, Mean Field Approximation confirms that lower density values, similar to those derived from MCS by employing atomistic model, fit better experimental isotherms than the higher densities obtained from other models.
By considering curved graphene surfaces one can explore the adsorption on other types of graphitic structures (carbon nanotubes, nanocones, and others kinds of fullerenes). In spite DFT calculations of the interaction energy between the curved graphene substrates and simple molecules, reveals that surface potential energy is almost independent on the curvature, at least when these curvatures range values corresponding to common nanotubes; these kind of carbonaceous structures have many others adsorption sites that have to be explored. As an example, single walled nanotubes (SWNT) form bundles of various tubes, and then, besides the previous mentioned convex external wall sites (S), where the adsorption is a two dimensional process, other kinds of sites can appear. One are sites lying between two adjacent outer nanotubes of the bundle or groove sites (G), associated to one dimensional adsorption, other sites are similar to those founded in cylindrical pores, and emerge when tubes diameter are so higher that they left interstitial channels (IC) with enough space in order to accommodate adsorbed molecules in them; and finally, and eventually they have inner sites (T) when there are SWNTs with their ends open. Beyond DFT calculations, where the energy of each sites confirm the postulated adsorption sequence: starts at the linear chains of G sites, and then in IC, if there are any one with enough space to accommodate molecules inside them, and finally on S sites. In case of open sites T sites have to be filled before G sites. MCS gives information about the shape of the isotherms, for instance, in case of open tubes, on the curve emerge a knee when a monolayer starts forming on the inside surface of the tubes, and also when this interior layer is completed.
By comparing these results with adsorption on graphite, one observes two facts: the monolayer of graphite is completed at a pressure below that corresponding to the nanotubes, and the isosteric heat of adsorption is greater for graphite, the reason is because the nanotubes considered are single walled (one graphene layer, rolled up) while graphite is formed by several graphene layers.
Results obtained from the study of gas mixtures adsorption (nitrogen / oxygen and ethane / ethylene) on different carbonaceous structures, show that graphite and more complex surfaces (such as carbon nanotubes) present similar trends. In this way, the behavior above described will provide valuable information on new systems.
Acknowledgments
The authors gratefully acknowledge financial support from the UNLP (Universidad Nacional de La Plata), CICPBA (Comisión de Investigaciones Científicas de la Provincia de Buenos Aires) and CONICET (Consejo de Investigaciones Científicas y Tecnológicas)
References
- 1.
Albesa A. G. Llanos J. L. Vicente J. L. 2008 Comparative study of methane adsorption on graphite. ,24 3836 3840 - 2.
Albesa A. G. Vicente J. L. 2008 Theoretical study of methane adsorption on graphite. ,95 1-2 ,10 22 - 3.
Albesa A. G. Filippin A. Vicente J. L. 2009 Effects of the Surface Curvature on the Adsorption of Simple Gases over Carbonaceous Materials using the Density Functional Theory (DFT) ,29 31 37 - 4.
Albesa A. G. Fertitta E. A. Vicente J. L. 2010 Comparative Study of Methane Adsorption on Single-Walled Carbon Nanotubes.,26 786 795 - 5.
Anson A. Wang Y. Lin C. C. H. Kuznicki T. M. Kuznicki S. M. 2008 Adsorption of ethane and ethylene on modified ETS-10. ,63 4171 4175 - 6.
Anson A. Lin C. C. H. Kuznicki T. M. Kuznicki S. M. 2010 Separation of ethylene/ethane mixtures by adsorption on small-pored titanosilicate molecular sieves.65 807 811 - 7.
Arora G. Sandler S. 2007 Nanoporous carbon membranes for separation of nitrogen and oxygen: Insight from molecular simulations. ,259 3 8 - 8.
Azami T. Kasuya D. Yuge R. Yudasaka M. Iijima S. Yoshitake T. Kubo Y. 2008 Large-Scale Production of Single-Wall Carbon Nanohorns with High Purity.112 1330 1334 - 9.
Bandosz T. J. Biggs M. J. Gubbins K. E. Hattori Y. Iiyima T. Kaneko K. Pikunic J. Thomson K. T. 2003 Molecular Models of Porous Carbons, In:28 Marcel Dekker, Inc., L.R. Radovic, (Ed),41 228 ,0-8247-0987-X New York, USA - 10.
Baudu M. Le Cloirec P. Martin G. 1993 First approach of desorption energies of water and organic molecules onto activated carbon by differential scanning calorimetry studies.,27 69 76 - 11.
Beaver M. Sircar S. 2010 Methane Reforming Using a Novel Sorption Enhanced Reaction Concept, In: , InTech, K. Eguchi (Ed.), Viena, Austria - 12.
Bertoncini C. Odetti H. Bottani E. J. 2000 Computer Simulation of Phenol Physisorption on Graphite.16 7457 7463 - 13.
Bhatia S. K. Myers A. L. 2006 Optimal Conditions for Adsorptive Storage,22 1688 1700 - 14.
Bienfait M. Zeppenfield P. Dupont-Pavlovsky N. Muris M. Johnson M. R. De Pie M. Vilches O. E. 2004 Thermodynamics and structure of hydrogen, methane, argon, oxygen, and carbon dioxide adsorbed on single-wall carbon nanotube bundles. ,70 035410 - 15.
Bottani E. J. Bakaev V. A. 1994 The Grand Canonical Ensemble Monte Carlo Simulation of Nitrogen on Graphite.,10 1550 1555 - 16.
Burley D. M. 1972 Col.2,329 In; Academic Press, C. Domb & M.S. Green (Eds.), London, Great Britain - 17.
Cheng A. Steele W. A. 1990 Computer Simulation Study of the Chemical Potential of Argon Adsorbed on Graphite. ,4 349 359 - 18.
Cheng H. Cooper A. C. Pez G. P. Kostov M. K. Piotrowski P. Stuart S. J. 2005 Molecular Dynamics Simulations on the Effects of Diameter and Chirality on Hydrogen Adsorption in Single Walled Carbon Nanotubes,109 3780 3786 - 19.
De Oliveira M. J. Griffiths R. S. 1978 Lattice-gas model of multiple layer adsorption. .,71 687 694 - 20.
Do D. D. Do H. D. 2005 Evaluation of 1-Site and 5-Site Models of Methane on Its Adsorption on Graphite and in Graphitic Slit Pores. ,109 19288 19295 - 21.
Dresselhaus M. S. Dresselhaus G. D. Eklund P. C. 1996 Academic Press (Ed.), New York, USA - 22.
Eldridge A. 1993 Olefin/paraffin separation technology: a review.32 2208 - 23.
Fuentes A. B. Mendez I. 2002 Separation and Purification Technology,28 29 - 24.
Gravelle P. C. J. 1978 Methods for the determination of heats of adsorption..,14 53 77 - 25.
Groszek A. J. 1998 Flow adsorption microcalorimetry ,312 133 143 - 26.
Hamilton J. J. Goodstein D. L. (1983) 28 28 3838 3848 - 27.
Harris P. J. F. 2003 Impact of the Discovery of Fullerenes on Carbon Science, In:28 Marcel Dekker, Inc., L.R. Radovic, (Ed),1 40 ,0-8247-0987-X New York, USA - 28.
Herdes C. Valente A. Lin Z. Rocha J. Coutinho J. A. P. Medina R. Vega L. F. 2007 Selective Adsorption of Volatile Organic Compounds in Micropore Aluminum Methylphosphonate-α: A Combined Molecular Simulation−Experimental Approach,23 7299 - 29.
Jayaraman A. Yang R. T. 2005 Stable oxygen-selective sorbents for air separation. ,60 625 634 - 30.
Kim H. Y. Steele W. A. 1992 Computer-simulation study of the phase diagram of the CH4 monolayer on graphite: Corrugation effects. ,45 6226 6233 - 31.
Kostov M. K. Cheng H. Cooper A. C. Pez G. P. 2002 Influence of Carbon Curvature on Molecular Adsorptions in Carbon-Based Materials: A Force Field Approach.,89 146105 146101 -4) - 32.
Kroon M. C. Vega L. F. 2009 Selective Paraffin Removal from Ethane/Ethylene Mixtures by Adsorption into Aluminum Methylphosphonate-α: A Molecular Simulation Study. ,25 2148 2152 - 33.
Krungleviciute V. Calvi M. M. Wagner J. A. Migone A. D. Yudasaka M. Iijima S. 2008 Probing the Structure of Carbon Nanohorn Aggregates by Adsorbing Gases of Different Sizes.,112 5742 5746 - 34.
Krungleviciute V. Migone A. D. Pepka M. 2009 Characterization of single-walled carbon nanohorns using neon adsorption isotherms.,47 769 774 - 35.
Kuznetsova A. Yates J. Liu J. Smalley R. E. 2000 112, 9590−9598 - 36.
Llanos J. L. Fertita A. E. Flores E. S. Bottani E. J. 2003 SO2 Physisorption on Exfoliated Graphite,107 8448 8453 - 37.
Myers A. L. Prausnitz J. M. 1965 Thermodynamics of mixed-gas adsorption.AIChE, J.,11 121 127 - 38.
Nicholson D. Parosonage N. G. 1982 , Academic Press (Ed.), London, Great Britain - 39.
Niyogi S. Hamon M. A. . Hu H. . Zhao B. Bownik S. 2002 Chemistry of Single-Walled Carbon Nanotubes. .,35 1105 1113 - 40.
Okamoto Y. Miyamoto Y. 2001 Ab Initio Investigation of Physisorption of Molecular Hydrogen on Planar and Curved Graphenes,105 17),3470 3474 - 41.
Parr R. G. Perdew J. P. 1989 Oxford University Press (Ed.), Oxford, Great Britain - 42.
Pascual P. Ungerer P. Taviatan B. Pernot P. Boutin A. 2003 Development of a transferable guest-host force field for adsorption of hydrocarbons in zeolitesI. Reinvestigation of alkane adsorption in silicalite by grand canonical Monte Carlo simulation. .,5 3684 3693 - 43.
Piper J. Morrison J. A. 1984 Heats of adsorption of methane multilayers on graphite,30 3486 3489 - 44.
Sabzyan H. Babajani M. 2005 A semiempirical quantum mechanical study of methane-graphite interaction.726 155 160 - 45.
Scheele C. W. 1777 Akad.d.Wiessenschaft Mitgliedes, Leipzig, Germany - 46.
Steele W. A. 1973 The physical interaction of gases with crystalline solids: I. Gas-solid energies and properties of isolated adsorbed atoms..,36 317 352 - 47.
Steele W. A. 1974 , Pergamon Press (Ed.),10-0080-177247 Great Britain - 48.
Talapatra S. Migone A. D. 2002 Adsorption of methane on bundles of closed-ended single-wall carbon nanotubes65 045416 1-6) - 49.
Tasca R. A. Calbi M. M. Cole M. W. 2002 Lattice model of gas condensation within nanopores.,65 061607 1-9) - 50.
Vilaplana A. F. 2005 Ab initio computational investigation of physisorption of molecular hydrogen on achiral single- walled carbon nanotubes.: 122, 214724-1-7 - 51.
Wagner K. 1996 Methane, In:13 IUPAC chemical data series, Blacwell Sciece (Ed.), Boston, USA - 52.
Wang Y. Perdew J. P. 1991 Correlation hole of the spin-polarized electron gas, with exact small-wave-vector y high-density scaling, : 44(24), 13298-13307 - 53.
Wang Y. Johnson J. 1999 Molecular simulation of hydrogen adsorption in single-walled carbon nanotubes y idealized carbon slit pores,110 577 - 54.
Wuebbles D. J. Hyhoe K. 2000 Atmospheric Methane: Trends and Impacts, In: Springer Verlag, J. van Ham (Ed.), Amsterdam, Holland - 55.
Yudasaka M. Iijima S. Crespi V. H. 2008 Single Wall Carbon Nanohorns and Nanocones, In: 111, 605, Springer Verlag, A. Jorio; G.D. Dresselhaus & M.S. Dresselhaus (Eds.), Berlin, Germany - 56.
Zettlemoyer A. C. Narayan K. S. 1966 ,2 197