Abstract

The adsorption of cyclohexene with two sp2 and four sp3 carbon atoms in graphitic slit pores was studied by performing grand canonical Monte Carlo simulation. The molecular arrangement of the cyclohexene on the graphitic carbon wall depends on the pore width. The distribution peak of the sp2 carbon is closer to the pore wall than that of the sp3 carbon except for the pore width of 0.7 nm, even though the Lennard-Jones size of the sp2 carbon is larger than that of the sp3 carbon. Thus, the difference in the interactions of the sp2 and sp3 carbon atoms of cyclohexene with the carbon pore walls is clearly observed in this study. The preferential interaction of sp2 carbon gives rise to a slight tilting of the cyclohexene molecule against the graphitic wall. This is suggestive of a π-π interaction between the sp2 carbon in the cyclohexene molecule and graphitic carbon.

1. Introduction

There has been an increase in the use of activated carbons in the construction of environmentally friendly technologies with sufficient safety. Activated carbon has basically slit-shaped micropores, which provide much better accessibility for most molecules than the cylindrical pores of zeolites [1, 2]. Recent researches have succeeded in controlling the pore width in the range from subnanometers to several nanometers; this has resulted in improvements in a wide range of applications, such as air separation (or air purification), solvent recovery, and manufacture of automobile canisters [37]. Furthermore, activated carbons have high electronic conductivities, allowing them to be used in electrochemical applications such as supercapacitors [810]. Another striking advantage of activated carbons is the fact that their production from natural products that fix atmospheric CO2 can contribute to CO2 sequestration [5, 11, 12]. However, there are still many research issues that need to be resolved to improve the performance of activated carbons. With the exception of the surface functional groups involved in the adsorption of polar molecules, the specific interaction involved in the gas adsorbability of activated carbons has not been studied sufficiently. Smith et al. [13] conducted the first study on the comparison between benzene and cyclohexane adsorptions on graphite; in this study, stronger cyclohexane adsorption was observed than benzene adsorption. On the other hand, the strong π-electron interaction of benzene with nanopore surfaces was observed by Dosseh et al. [14]. Moreover, comparison studies on adsorbed molecules of sp2 and sp3 carbons have been conducted for benzene and cyclohexane and for ethylene and ethane [1517]. The adsorption isotherms of ethane and ethylene on activated carbons showed larger amounts of adsorbed ethylene than ethane [18], indicating stronger interaction of sp2 carbon with pore walls. Radovic and Bockrath used molecular orbital calculations to emphasise the essential importance of the conjugated π-electron nature of the basal plane of the graphite in the adsorption of organic molecules [19]. The characteristic interaction, the π-π interaction of the sp2 carbons in an adsorbed molecule with graphite, is caused by the electron cloud delocalization between the sp2 carbons in the molecule and graphite. We still need to understand the difference between the adsorbed structures of the sp2 and sp3 carbon atoms of the adsorbate molecules upon adsorption on graphitic pores along with the dependence on pore width in order to obtain a better carbon adsorbent of high specificity for use in future technologies for the selective reaction and separation of unsaturated molecules. Vernov and Steele showed the orientation of adsorbed benzene on graphite using a 12-site model [2022]. Do et al. studied the orientation structures of organic molecules on graphite and porous carbons using grand canonical Monte Carlo (GCMC) simulation, taking into account the molecular structures [23, 24]. Such molecular simulation studies can reveal the interactions between sp2 carbon and a graphitic carbon surface along with the molecular structures.

Cyclic hydrocarbons are important intermediates in many chemical reactions on transition metal surfaces [25, 26]. Cyclohexene is one of the critical intermediates in catalytic reforming. Many studies on the reforming process have been reported [2527]. The cyclohexene molecule has two sp2 carbon atoms in addition to four sp3 carbon atoms, making it an appropriate probe molecule for research on the comparison of the adsorptions of sp2 and sp3 carbon atoms. To understand the difference between the adsorbed structures of the sp2 and sp3 carbon atoms, we adopted GCMC simulation for studying cyclohexene adsorption on activated carbon. In addition, this study could serve as a guide in the development of better adsorption media for separating cyclohexene from commercial natural gas [4, 6, 15].

2. Simulation Procedures

Figure 1 shows a space-filled model of cyclohexene; the structural parameters are given in the paper by Faller et al. [28]. The cyclohexene molecule has a distorted plane structure that consists of five single bonds (C–C) and one double bond (C=C). The sp3 and sp2 carbon atoms must be distinguished by using interaction potential calculation and GCMC simulation. Faller et al. determine the parameters of cyclohexene for the flexible model [28]. As rigid and flexible models give no significant difference in structure [29, 30], a 16-centre model for the rigid cyclohexene structure in this study was adopted for the calculation ( K for H and 35.6 K for sp3 and sp2 carbons;  nm for H, 0.311 nm for sp3 carbon, and 0.321 nm for sp2 carbon) [28, 31]. The availability is verified by the comparison between the rigid and flexible models before the calculation; the assembly structures of cyclohexene in the rigid and flexible models are almost agreeing with each other, although the structure distribution in the flexible model is slightly dispersed due to the flexibility. The intermolecular interactions were described by the Lennard-Jones potentials:

Here, is the distance between the th atom in molecule and the th atom in molecule . The Lorentz-Berthelot mixing rules were applied to obtain the interaction energy and size parameters for the heteroatomic interaction:

Here, each component atom in the molecule was assumed to be electrically neutral; the stable molecular geometry of a single molecule in the flexible model was used for the calculation.

The slit-shaped pore was modelled using the interface between two semi-infinite graphite slabs; the molecule-carbon wall interaction was approximated by the 10-4-3 Steele potential [32]:

Here, is 76.38 π, and is the vertical distance of the th atom in a molecule from the centre of a carbon atom on one side of a carbon surface; and were derived from the Lorentz-Berthelot mixing rules. and are the Lennard-Jones parameters for the carbon atom in graphite ( nm and K), and and are the Lennard-Jones parameters for the th atom in the molecule [33]. Strictly speaking, the molecule-graphite slab interaction was summed and smoothed by all of the interactions between the component atoms in the molecule and the graphite slab using the Steele potential. In the case of the graphite slit pore, the interaction of the cyclohexene molecule with the pore was expressed by the sum of the interaction potentials of the cyclohexene molecule with both graphite walls, as given by (4):

Here, is the physical slit pore width, which is the internuclear distance between opposite graphite walls. The physical slit pore width was associated with the pore width, , which can be experimentally measured as follows [34]:

Here, is the closest contact distance between the sp2 carbon and the graphite pore wall.

The electrostatic interaction between a polar molecule and a graphitic slab must be taken into account for a detailed analysis of the adsorption of polar molecules, as reported in previous papers [20, 35]. The electrically neutral cyclohexene atoms assumed here nominally give no electrostatic contribution to the interaction between the cyclohexene and the pore wall. However, the electrostatic interaction is roughly included in the Lennard-Jones parameters obtained from the ab initio calculation by Faller et al. [28]. Grand canonical Monte Carlo simulations of the cyclohexene were performed using 3 × 106 steps at 298 K, with three equivalent trials for the creation, deletion, and movement of the molecules. The pressure was calculated from the bulk molecular density in the unit cell of 6 × 6 × 6 nm3 by using the van der Waals equation after more than 1 × 108 calculation steps and the accumulation of  9 × 107 steps. A unit cell size of 6 × 6 × H nm3 and a 2-dimensional periodic boundary condition were used in this calculation.

3. Results and Discussion

The molecule-graphite pore interaction energy was calculated for three arrangements, as shown in Figure 2. The molecular plane is perpendicular to the graphite wall in arrangements (a) and (b). The double bond is head-on toward a plus-side wall for (a), whereas there is side-on adsorption on the wall for (b). Hence, (a) and (b) are called head-on and side-on arrangements, respectively. On the other hand, because the molecular plane is parallel to the wall in (c), arrangement (c) is called an in-plane configuration. The molecule-pore interaction potentials were calculated as a function of the vertical distance of the molecular centre of gravity from the pore centre for the three arrangements.

Figure 3 shows the molecule-pore interaction potentials for the three arrangements shown in Figure 2. For the smallest pore ( nm), the in-plane arrangement gets a large stabilization energy of −8000 K, suggesting the presence of adsorption from an extremely low pressure. The side-on arrangement provides the considerably large stabilization energy of −4000 K, whereas the head-on arrangement leads only to repulsive interaction energy. The potential profiles for  nm are remarkably different from those for  nm. The most stable arrangement (−7000 K) is found to be the side-on type with  nm. The in-plane arrangement gives double potential minima, and the head-on and side-on arrangements give a single minimum. All the potential minima are close to each other, within 500 K, for  nm. The head-on arrangement leads to the single deepest potential minimum for  nm, and the side-on and in-plane arrangements have double minima. In the case of  nm, the in-plane arrangement gave the deepest potential double minima, which are situated near the pore wall. Then, cyclohexene molecules on the pore walls will be adsorbed on the pore walls with  nm in the in-plane structure. Similar double potential minima were obtained for the side-on arrangement, and the depth was slightly shallower than that for the in-plane arrangement. It is worth noting that the head-on arrangement gave an asymmetrical potential profile for  nm; sp2 carbon leads to a stronger interaction than sp3 carbon. Even the single potential minima for and 0.8 nm shift slightly in the positive direction due to the contribution by the sp2 carbons (see Figures 3(b) and 3(c)). This interaction potential difference in the sp2 and sp3 carbon atoms is clearly observed in the adsorbed structure obtained from the GCMC simulation. The potential profiles for  nm have features that are almost similar to those for  nm (Figure 3(d)). Figure 4 shows the changes in the interaction potential minimum with the pore width for different adsorption structures. The interaction energy depends sensitively on the pore width and the molecular arrangement in the pore; the most favourable arrangement predicted from the interaction potential varies with the pore width. The deepest potential energy is obtained in the in-plane arrangement for  nm. The side-on arrangement provides the deepest potential energy in the range of 0.64 to 0.75 nm. The head-on arrangement gives the deepest interaction potential energy for  nm, indicating a contribution by the sp2 carbon atoms in the pores when is approximately 0.8 nm, although the energy difference from the side-on arrangement is not marked compared with the difference for smaller pores.

Figure 5 shows the adsorption isotherms of cyclohexene in pores with  nm at 298 K. The ordinate represents the number of cyclohexene molecules adsorbed in the pore, with wider pores resulting in larger numbers of molecules above 10−5 MPa. As the pore with  nm is too narrow to accommodate cyclohexene molecules easily, adsorption begins above 10−8 MPa. On the other hand, the pore with  nm, which has the deepest potential minimum (−8000 K) for the in-plane arrangement, exhibits a quite intensive adsorption for cyclohexene; this adsorption begins even below 10−10 MPa. As the pore with  nm has double minima, which are in the range from −4000 to −4500 K, adsorption begins at 5 × 10−9 MPa. The adsorption for  nm begins above 5 × 10−9 MPa. The cyclohexene molecule has a distorted ring structure, making dense packing in a smaller slit pore space very difficult. This can be understood from the following snapshot analysis.

Figure 6 shows the snapshots at 10−3 MPa for , and 1.0 nm. Cyclohexene molecules form a single adsorbed layer that is inherent to the pore width for  nm, as shown in Figure 4. In the pore with  nm, the cyclohexene molecules basically exhibit the in-plane arrangement. Strictly speaking, the cyclohexene molecules are tilted against the pore wall due to the stronger sp2 carbon-pore wall interaction. The tilt angle will be shown later. Single adsorbed layers having the side-on and head-on arrangements are observed in the pores with and 0.8 nm, respectively, which is expected from the results shown in Figure 4. Two adsorbed layers of cyclohexene molecules are formed in the pore with  nm. Three layers and four layers are observed in the pores with and 2.0 nm, respectively, which are not shown in this paper. In the case of the pores with and 2.0 nm, molecules at the monolayer position orientate along the pore walls; molecules in the central space of the pore have almost a perpendicular configuration against the pore wall with less molecular orientation.

The tilting of the cyclohexene molecules toward the pore walls clearly shows the molecular orientation to a pore wall. Figure 7 shows the angular distribution between the normal vector of a cyclohexene molecule and the pore wall surface. Here, the normal vector is defined as a perpendicular vector to the plane of two sp2 carbons and two sp3 carbons next to the sp2 carbon; that is, the in-plane arrangement is at 0° and the head-on or side-on arrangement is at 90°. The probable angles for and 0.7 nm are 20–30° and 50–90°, respectively. Thus, the in-plane arrangement is expected for  nm, while the head-on or side-on arrangement should be formed in the pore with  nm. For and 1.0 nm, angles of 20–30° are observed in the contacts with the pore walls. Angle distributions above 50° are also observed at some pore-centred positions, suggesting the head-on or side-on arrangement. The head-on or side-on arrangement is mainly formed for  nm. On the other hand, the in-plane arrangement is the most common for  nm.

As the exact arrangement of the molecules facing the pore wall can provide useful information on the role of the sp2 carbon atoms in the molecule-pore wall interaction, distribution profiles were determined for the sp2 and sp3 carbons in the adsorbed cyclohexene against the perpendicular axis of the pore walls, as shown in Figure 8. The snapshots in Figure 6 and angular distribution in Figure 7 cannot provide a clear insight into the role of the sp2 carbon atoms in the cyclohexene-pore wall interaction. Cruz and Mota found a slight difference between the distances of the sp2 and sp3 carbons from pore walls [16]. In contrast, Figure 8 distinctly shows the difference in the positions of the sp2 and sp3 carbon atoms and the hydrogen atoms bonding to the sp2 and sp3 carbons. The sp2 carbon atoms are actually in contact with a pore wall for and 1.0 nm, although the remarkable difference between the sp2 and sp3 carbon atoms near the walls is not observed well for and 0.8 nm. The sp2 carbon has a considerably uniform distribution for  nm, and a few sp2 carbons are closer to the wall than the sp3 carbon, whereas the sp3 carbon prefers the monolayer position. The distributions for and 2.0 nm had tendencies similar to that for  nm; these distributions for nm have sharp peaks near the pore walls, and no clear peak exists at around the centre of the pore. The distribution peak of the sp2 carbon is closer to the pore walls than that of the sp3 carbon atoms for , and 2.0 nm, despite the fact that the Lennard-Jones diameter of the sp2 carbon is larger than that of the sp3 carbon by 0.01 nm. Therefore, the sp2 carbon atoms can interact preferentially with the graphite walls in comparison with the interaction of the sp3 carbon atoms for the specific pore widths. The hydrogen atoms bound to the sp2 carbons are also distributed on the pore walls, whereas those bound to the sp3 carbons are orderly distributed in the pores. In comparison to the carbon atoms, there exist nearer hydrogen atoms to the pore walls.

4. Conclusion

The preferential interaction of sp2 carbon atoms with a graphite wall is explicitly observed in the adsorption of cyclohexene in graphite slit pores, with the exception of pores with and 0.8 nm, as cyclohexene has both sp2 and sp3 carbon atoms. The preferential interactions of the sp2 carbon atoms and the hydrogen atoms bound to the sp2 carbon atoms in cyclohexene are enhanced in a restricted nanoscale pore space; this results in cyclohexene having a slightly tilted conformation. The application of the preferential interaction nature of the sp2 carbon atoms to the design of a better adsorbent for aromatic compounds is expected to be quite useful, because the orientation of cyclohexene to the graphite walls depends on the pore width.

Acknowledgments

The financial support provided by the Kurata Memorial Hitachi Science and Technology Foundation, Kurita Water and Environment Foundation, Nippon Sheet Glass Foundation, and Global COE program, MEXT, Japan, is gratefully acknowledged.