Adsorption on Edge-Functionalized Bilayer Graphene Nanoribbons: Assessing the Role of Functional Groups in Methane Uptake

A combination of ab initio and classical Monte Carlo simulations is used to investigate the effects of functional groups on methane binding. Using Moller-Plesset (MP2) calculations, we obtain the binding energies for benzene functionalized with NH2, OH, CH3, COOH, and H2PO3 and identify the methane binding sites. In all cases, the preferred binding sites are located above the benzene plane in the vicinity of the benzene carbon atom attached to the functional group. Functional groups enhance methane binding relative to benzene (-6.39 kJ/mol), with the largest enhancement observed for H2PO3 (-8.37 kJ/mol) followed by COOH and CH3 (-7.77 kJ/mol). Adsorption isotherms are obtained for edge-functionalized bilayer graphene nanoribbons using grand canonical Monte Carlo simulations with a five-site methane model. Adsorbed excess and heats of adsorption for pressures up to 40 bar and 298 K are obtained with functional group concentrations ranging from 3.125 to 6.25 mol 96 for graphene edges functionalized with OH, NH2, and COOH. The functional groups are found to act as preferred adsorption sites, and in the case of COOH the local methane density in the vicinity of the functional group is found to exceed that of bare graphene. The largest enhancement of 44.5% in the methane excess adsorbed is observed for COOH-functionalized nanoribbons when compared to H terminated ribbons. The corresponding enhancements for OH- and NH2-functionalized ribbons are 10.5% and 3.7%, respectively. The excess adsorption across functional groups reflects the trends observed in the binding energies from MP2 calculations. Our study reveals that specific site functionalization can have a significant effect on the local adsorption characteristics and can be used as a design strategy to tailor materials with enhanced methane storage capacity.

[1]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[2]  Pavel Hobza,et al.  Stabilization and structure calculations for noncovalent interactions in extended molecular systems based on wave function and density functional theories. , 2010, Chemical reviews.

[3]  Ashley L. Ringer,et al.  Aliphatic C-H/pi interactions: Methane-benzene, methane-phenol, and methane-indole complexes. , 2006, The journal of physical chemistry. A.

[4]  Christian J. Doonan,et al.  Multiple Functional Groups of Varying Ratios in Metal-Organic Frameworks , 2010, Science.

[5]  Sridhar Komarneni,et al.  Porous Adsorbents for Vehicular Natural Gas Storage: A Review , 1998 .

[6]  B. Saha,et al.  Surface modification and characterisation of a coal-based activated carbon , 2005 .

[7]  Taner Yildirim,et al.  Graphene oxide derived carbons (GODCs): synthesis and gas adsorption properties , 2012 .

[8]  J. Ilja Siepmann,et al.  Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-Alkanes , 1998 .

[9]  D. Do,et al.  Evaluation of 1-site and 5-site models of methane on its adsorption on graphite and in graphitic slit pores. , 2005, The journal of physical chemistry. B.

[10]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[11]  S. Tsuzuki,et al.  Magnitude of the CH/pi interaction in the gas phase: experimental and theoretical determination of the accurate interaction energy in benzene-methane. , 2006, The journal of physical chemistry. A.

[12]  S. Narasimhan,et al.  Methane and carbon dioxide adsorption on edge-functionalized graphene: a comparative DFT study. , 2012, The Journal of chemical physics.

[13]  Kazumasa Honda,et al.  The Magnitude of the CH/π Interaction between Benzene and Some Model Hydrocarbons , 2000 .

[14]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[15]  Jacopo Tomasi,et al.  Counterpoise corrections to the interaction energy components in bimolecular complexes , 1985 .

[16]  S. Komarneni,et al.  Probing the pore structure and dehumidifying characteristics of chemically modified zeolite Y , 1995 .

[17]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[18]  Dolores Lozano-Castelló,et al.  Advances in the study of methane storage in porous carbonaceous materials , 2002 .

[19]  C. Wilmer,et al.  Large-scale screening of hypothetical metal-organic frameworks. , 2012, Nature chemistry.

[20]  Edward G Hohenstein,et al.  An assessment of theoretical methods for nonbonded interactions: comparison to complete basis set limit coupled-cluster potential energy curves for the benzene dimer, the methane dimer, benzene-methane, and benzene-H2S. , 2009, The journal of physical chemistry. A.

[21]  Duncan H. Gregory,et al.  Hydrogen storage materials: present scenarios and future directions , 2009 .

[22]  I. Martín-Gullón,et al.  Chapter 1 Types of carbon adsorbents and their production , 2006 .

[23]  Peter A. Kollman,et al.  Simulation of the solvation free energies for methane, ethane, and propane and corresponding amino acid dipeptides : a critical test of the bond-PMF correction, a new set of hydrocarbon parameters, and the gas phase-water hydrophobicity scale , 1992 .