Potential models for the simulation of methane adsorption on graphene: development and CCSD(T) benchmarks

Theoretical insights into the methane catalytic decomposition on graphene nanoribbons edges

Catalytic methane decomposition (CMD) is receiving much attention as a promising application for hydrogen production. Due to the high energy required for breaking the C-H bonds of methane, the choice of catalyst is crucial to the viability of this process. However, atomistic insights for the CMD mechanism on carbon-based materials are still limited. Here, we investigate the viability of CMD under reaction conditions on the zigzag (12-ZGNR) and armchair (AGRN) edges of graphene nanoribbons employing dispersion-corrected density functional theory (DFT). First, we investigated the desorption of H and H₂ at 1200 K on the passivated 12-ZGNR and 12-AGNR edges. The diffusion of hydrogen atom on the passivated edges is the rate determinant step for the most favourable H₂ desorption pathway, with a activation free energy of 4.17 eV and 3.45 eV on 12-ZGNR and 12-AGNR, respectively. The most favourable H₂ desorption occurs on the 12-AGNR edges with a free energy barrier of 1.56 eV, reflecting the availability of bare carbon active sites on the catalytic application. The direct dissociative chemisorption of CH₄ is the preferred pathway on the non-passivated 12-ZGNR edges, with an activation free energy of 0.56 eV. We also present the reaction steps for the complete catalytic dehydrogenation of methane on 12-ZGNR and 12-AGNR edges, proposing a mechanism in which the solid carbon formed on the edges act as new active sites. The active sites on the 12-AGNR edges show more propensity to be regenerated due lower free energy barrier of 2.71 eV for the H₂ desorption from the newly grown active site. Comparison is made between the results obtained here and experimental and computational data available in the literature. We provide fundamental insights for the engineering of carbon-based catalysts for the CMD, showing that the bare carbon edges of graphene nanoribbons have performance comparable to commonly used metallic and bi-metallic catalysts for methane decomposition.

Development of accurate potentials for the physisorption of water on graphene

From coupled-cluster singles and doubles model including connected triples corrections [CCSD(T)] calculations on the water dimer and B97D/CC on the water-circumcoronene complex at a large number of randomly generated conformations, interaction potentials for the physisorption of water on graphene are built, accomplishing almost sub-chemical accuracy. The force fields were constructed by decomposing the interaction into electrostatic and van der Waals contributions, the latter represented through improved Lennard-Jones potentials. Besides, a Chemistry at Harvard Macromolecular Mechanics (CHARMM)-like term was included in the water-water potential to improve the description of hydrogen bonds, and an induction term was added to model the polarization effects in the interaction between water and polyaromatic hydrocarbons (PAHs) or graphene. Two schemes with three and six point charges were considered for the interactions water-water and water-PAH, as Coulomb contributions are zero in the water-graphene system. The proposed fitted potentials reproduce the ab initio data used to build them in the whole range of distances and conformations and provide results for selected points very close to CCSD(T) benchmarks. When applied to the water-graphene system, the obtained results are in excellent agreement with p-CCSD(T), revised symmetry-adapted perturbation theory based on density functional theory monomer properties (DFT-SAPT), and diffusion Monte Carlo reference values. Furthermore, the stability of the various conformers water-PAH and water-graphene, as well as the different trends observed between these systems are rationalized in terms of the modifications of the electrostatic contribution.

Mechanical and gas adsorption properties of graphene and graphynes under biaxial strain

The exceptional properties of two-dimensional (2D) solids have motivated extensive research, which revealed the possibility of controlling many characteristics of these materials through strain. For instance, previous investigations demonstrated that compressive deformation could be used to direct the chemisorption of atomic hydrogen and oxygen. Still, to our knowledge, there is no work detailing how strain affects the adsorption isotherms of 2D materials and the adsorption properties of materials such as the graphynes, which are monolayers composed of sp and sp[Formula: see text] carbon atoms. In the present work, we analyze how biaxial tensile deformation changes the adsorption properties of four 2D materials (graphene, [Formula: see text]-graphyne, [Formula: see text]-graphyne, and [Formula: see text]-graphyne). To achieve this, we perform Monte Carlo Grand Canonical calculations to obtain the adsorption isotherms of H[Formula: see text], CO[Formula: see text], and CH[Formula: see text] on the monolayers with and without strain. And, to apply the deformation, we carry out Molecular Dynamics simulations. We find a substantial reduction in the amount of gas adsorbed on the monolayers for nearly all gas-solid combinations. This is particularly true for graphene, where 14.5% strain reduces the quantity of H[Formula: see text]/CO[Formula: see text]/CH[Formula: see text] by 44.7/64.1/41.7% at P [Formula: see text] 1 atm. To understand the results, we calculate adsorption enthalpies and analyze the gas distribution above the monolayers. We also characterize the mechanical properties of the considered solids under biaxial deformation. Finally, a comparison of pore sizes with the kinetic diameters of various gases suggests applications for the graphynes, with and without strain, in gas separation.

Multilayer Graphtriyne Membranes for Separation and Storage of CO2: Molecular Dynamics Simulations of Post-Combustion Model Mixtures

The ability to remove carbon dioxide from gaseous mixtures is a necessary step toward the reduction of greenhouse gas emissions. As a contribution to this field of research, we performed a molecular dynamics study assessing the separation and adsorption properties of multi-layered graphtriyne membranes on gaseous mixtures of CO₂, N₂, and H₂O. These mixtures closely resemble post-combustion gaseous products and are, therefore, suitable prototypes with which to model possible technological applications in the field of CO₂ removal methodologies. The molecular dynamics simulations rely on a fairly accurate description of involved force fields, providing reliable predictions of selectivity and adsorption coefficients. The characterization of the interplay between molecules and membrane structure also permitted us to elucidate the adsorption and crossing processes at an atomistic level of detail. The work is intended as a continuation and a strong enhancement of the modeling research and characterization of such materials as molecular sieves for CO₂ storage and removal.

Extended Line Defect Graphene Modified by the Adsorption of Mn Atoms and Its Properties of Adsorbing CH4

Extended line defect (ELD) graphene is a two-dimensional (2D) topologically defective graphene with alternate octagonal and quadrilateral carbon rings as basic defective units. This paper reports on the CH4 adsorption properties of ELD graphene according to the first principles of density functional theory (DFT). The effects on the CH4 adsorption of ELD graphene when modified by a single Mn atom or two Mn atoms were investigated, respectively. An ELD-42C graphene configuration consisting of 42 C atoms was first constructed. Then, the ELD-42C graphene configuration was used as a substrate, and a Mn-ELD-42C graphene configuration was obtained by modifying it with a single Mn atom. The results showed that the most stable adsorption site for Mn atoms was above the quadrilateral carbon ring. This Mn-ELD-42C graphene configuration could only stably adsorb up to 30 CH4 molecules on each side, with an average adsorption energy of −0.867 eV/CH4 and an adsorption capacity of 46.25 wt%. Three 2Mn-ELD-42C graphene configurations were then obtained by modifying the ELD-42C graphene substrate with two Mn atoms. When the two Mn atoms were located on either side of a 2Mn-ELD-42C graphene configuration and above the two octagonal carbon rings adjacent to the same quadrilateral carbon ring, it was able to adsorb up to 40 CH4 molecules on each side, with an average adsorption energy of −0.862 eV/CH4 and a CH4 adsorption capacity of 51.09 wt%.

Grand Canonical Monte Carlo Simulations to Determine the Optimal Interlayer Distance of a Graphene Slit-Shaped Pore for Adsorption of Methane, Hydrogen and their Equimolar Mixture

The adsorption—for separation, storage and transportation—of methane, hydrogen and their mixture is important for a sustainable energy consumption in present-day society. Graphene derivatives have proven to be very promising for such an application, yet for a good design a better understanding of the optimal pore size is needed. In this work, grand canonical Monte Carlo simulations, employing Improved Lennard–Jones potentials, are performed to determine the ideal interlayer distance for a slit-shaped graphene pore in a large pressure range. A detailed study of the adsorption behavior of methane, hydrogen and their equimolar mixture in different sizes of graphene pores is obtained through calculation of absolute and excess adsorption isotherms, isosteric heats and the selectivity. Moreover, a molecular picture is provided through z-density profiles at low and high pressure. It is found that an interlayer distance of about twice the van der Waals distance of the adsorbate is recommended to enhance the adsorbing ability. Furthermore, the graphene structures with slit-shaped pores were found to be very capable of adsorbing methane and separating methane from hydrogen in a mixture at reasonable working conditions (300 K and well below 15 atm).

A DFT investigation of lithium adsorption on graphenes as a potential anode material in lithium-ion batteries

We present a detailed study of the Li⁺ ion adsorption on two different hydrogenated carbon nanostructures, namely as pristine graphene (PG) and topologic Stone-Wales defective graphene (SWG) using the density functional theory (DFT). The studies are focused to analyze the structure-stability relationship with the estimated electronic and electrical properties for lithium-ion batteries (LIB) formed with an anode based on the Li/Li⁺#PG and Li/Li⁺#SWG systems. In addition, the electronic effects induced due to Li⁺ adsorption and the presence of SW defect on the graphene models were analyzed by the frontier molecular orbitals, ChelpG charges, Raman and UV-Vis spectra. It was verified that Li⁺ is more stably adsorbed on the edges on both graphene structures through an electrostatic interaction between cation and more negatively charged edges of nanostructures. TD-DFT calculations showed that the metallic nature of isolated graphene is disturbed after the adsorption of Li⁺, and this was demonstrated from the calculated HOMO-LUMO gap. The same Li⁺-Graphene geometries were optimized by introducing neutral charge in order to enable the calculation of ionization potentials. I was also found that such systems potentially contributed to the modeling of graphene-based anodes with reasonable electrical voltage responses estimated for a LIB. The simulation of Raman and UV-Vis spectra revealed significant variations in intensity and shifts the typical bands of graphene due to the presence of the Li⁺ ion that can contribute to point out new experiments to the spectroscopic characterization of these systems. Our results suggest that these carbon nanostructures are potential candidates for efficient applications in electrochemical systems, mainly dealing with LIB.

Two-dimensional diamine-linked covalent organic frameworks for CO2/N2 capture and separation: theoretical modeling and simulations

Two-dimensional covalent organic frameworks (2D-COFs) with diamine-based linkers have been designed and investigated for CO2/N2 gaseous mixture adsorption and separation via a systematic theoretical study by combining density functional theory (DFT) calculations and force field-based molecular dynamics (MD) simulations. We explored the adsorption sites and adsorption energies of CO2/N2 on 2D-COFs. The gas uptake capacity, adsorption isotherms, permeability, and selectivity were simulated based on an improved formulation of force fields for mixture separation in post-combustion conditions. This theoretical approach provided atomistic understanding and quantitative description of intermolecular interactions governing the physisorption dynamics of the considered systems. The results suggest that 2D-COFs investigated in this study are competitive with other 2D materials for carbon capture and separation and can be considered as alternative molecular sieving materials offering efficient and rapid separation and adsorption of different molecules.

Molecular Dynamics of CH4/N2 Mixtures on a Flexible Graphene Layer: Adsorption and Selectivity Case Study

We theoretically investigate graphene layers, proposing them as membranes of subnanometer size suitable for CH₄/N₂ separation and gas uptake. The observed potential energy surfaces, representing the intermolecular interactions within the CH₄/N₂ gaseous mixtures and between these and the graphene layers, have been formulated by adopting the so-called Improved Lennard-Jones (ILJ) potential, which is far more accurate than the traditional Lennard-Jones potential. Previously derived ILJ force fields are used to perform extensive molecular dynamics simulations on graphene's ability to separate and adsorb the CH₄/N₂ mixture. Furthermore, the intramolecular interactions within graphene were explicitly considered since they are responsible for its flexibility and the consequent out-of-plane movements of the constituting carbon atoms. The effects on the adsorption capacity of graphene caused by introducing its flexibility in the simulations are assessed via comparison of different intramolecular force fields giving account of flexibility against a simplified less realistic model that considers graphene to be rigid. The accuracy of the potentials guarantees a quantitative description of the interactions and trustable results for the dynamics, as long as the appropriate set of intramolecular and intermolecular force fields is chosen. In particular it is shown that only if the flexibility of graphene is explicitly taken into account, a simple united-atom interaction potential can provide correct predictions. Conversely, when using an atomistic model, neglecting in the simulations the intrinsic flexibility of the graphene sheet has a minor effect. From a practical point of view, the global analysis of the whole set of results proves that these nanostructures are versatile materials competitive with other carbon-based adsorbing membranes suitable to cope with CH₄ and N₂ adsorption.