An Effective Potential for Adsorption of Polar Molecules on Graphite

Bottom-Up Approach to the Coarse-Grained Surface Model: Effective Solid-Fluid Potentials for Adsorption on Heterogeneous Surfaces

Coarse-grained surface models with a low-dimension positional dependence have great advantages in simplifying the theoretical adsorption model and speeding up molecular simulations. In this work, we present a bottom-up strategy, developing a new two-dimensional (2D) coarse-grained surface model from the "bottom-level" atomistic model, for adsorption on highly heterogeneous surfaces with various types of defects. The corresponding effective solid-fluid potential consists of a 2D hard wall potential representing the structure of the surface and a one-dimensional (1D) effective area-weighted free-energy-averaged (AW-FEA) potential representing the energetic strength of the substrate-adsorbate interaction. Within the conventional free-energy-averaged (FEA) framework, an accessible-area-related parameter is introduced into the equation of the 1D effective solid-fluid potential, which allows us not only to obtain the energy information from the fully atomistic system but also to get the structural dependence of the potential on any geometric defect on the surface. Grand canonical Monte Carlo simulations are carried out for argon adsorption at 87.3 K to test the validity of the new 2D surface model against the fully atomistic system. We test four graphitic substrates with different levels of geometric roughness for the top layer, including the widely used reference solid substrate Cabot BP-280. The simulation results show that adding one more dimension to the traditional 1D surface model is essential for adsorption on the geometrically heterogeneous surfaces. In particular, the 2D surface model with the AW-FEA solid-fluid potential significantly improves the adsorption isotherm and density profile over the 1D surface model with the FEA solid-fluid potential over a wide range of pressure. The method to construct an effective solid-fluid potential for an energetically heterogeneous surface is also discussed.

Assessment of Density Functional Theory in Predicting Interaction Energies between Water and Polycyclic Aromatic Hydrocarbons: from Water on Benzene to Water on Graphene

The interactions of water with polycyclic aromatic hydrocarbons, from benzene to graphene, are investigated using various exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. The accuracy of the different functionals is assessed through comparisons with random phase approximation (RPA) and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] calculations. Diffusion Monte Carlo (DMC) data reported in the literature are also used for comparison. Relatively large variations are found in interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules. The meta-GGA B97M-rV and range-separated hybrid, meta-GGA ωB97M-V functionals provide nearly quantitative agreement with CCSD(T) values for the water-benzene, water-coronene, and water-circumcoronene dimers, while RPA and DMC predict interaction energies that differ by up to ∼1 kcal/mol and ∼0.4 kcal/mol from the corresponding CCSD(T) values, respectively. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for larger polycyclic aromatic hydrocarbons. By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA), it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet.

Preparation of Graphite Oxide Containing Different Oxygen-Containing Functional Groups and the Study of Ammonia Gas Sensitivity

A series of graphite oxide samples were prepared using the modified Hummers method. Flake graphite was used as the raw material and the reaction temperature of the aqueous solution was changed (0 °C, 30 °C, 50 °C, 60 °C, 70 °C, 80 °C, and 100 °C). X-ray diffraction, Fourier-transform infrared spectroscopy, Raman spectral analysis, X-ray photoelectron spectroscopy, and contact angle tests were performed to characterize the structure, chemical bonding, type, and content of oxygen-containing functional groups of the graphite oxide samples. The results showed that the type and content of each oxygen-containing functional group could be controlled by changing the reaction temperature with the addition of water. As the temperature of the system increased, the degree of oxidation of the graphite oxide samples first increased and then decreased. Too high a temperature (100 °C) of the system led to the formation of epoxy groups by the decomposition of some hydroxyl groups in the samples, causing the reduction of oxygen-containing functional groups between the graphite layers, poor hydrophilic properties, and low moisture content. When the system temperature was 50 °C, the interlayer spacing of the graphite oxide samples was at its highest, the graphite was completely oxidized (C/O = 1.85), and the oxygen-containing functional groups were mainly composed of hydroxyl groups (accounting for approximately 28.88% of the total oxygen-containing functional groups). The high content of hydroxyl and carboxyl groups had good hydrophilic ability and showed the highest moisture content. The sample at 50 °C had better sensitivity to ammonia because of its high hydroxyl group and carboxyl group content, with the sample showing an excellent profile when the ammonia concentration was 20⁻60 ppm.

Adsorption and phase behavior of water-like fluid models with square-well attraction and site-site association in slit-like pores: Density functional approach

The adsorption and phase behavior of two model fluids, both with square well inter-particle attraction and site-site associative interaction, in slit-like pores have been studied in the framework of a density functional theory. The mean field approach and the first-order mean spherical approximation have been applied to account for the attractive interactions. The chemical association effects are taken into account by using the first-order thermodynamic perturbation theory of Wertheim. A set of parameters for each fluid model has been chosen according to the work of [Clark et al., Mol. Phys. 104, 3561 (2006)], to describe successfully the vapor-liquid coexistence of water in the bulk phase. The influence of the slit-like pore width and of the strength of gas-solid interaction energy on the vapor-liquid coexistence envelope under confinement has been explored in detail. The theory and the results of the present work are valuable for further exploration of a wide set of models of associating fluids and of fluids with complex molecular architecture in different adsorbents, and to deal with activated carbon surfaces.

Computational chemistry for graphene-based energy applications: progress and challenges

Research in graphene-based energy materials is a rapidly growing area. Many graphene-based energy applications involve interfacial processes. To enable advances in the design of these energy materials, such that their operation, economy, efficiency and durability is at least comparable with fossil-fuel based alternatives, connections between the molecular-scale structure and function of these interfaces are needed. While it is experimentally challenging to resolve this interfacial structure, molecular simulation and computational chemistry can help bridge these gaps. In this Review, we summarise recent progress in the application of computational chemistry to graphene-based materials for fuel cells, batteries, photovoltaics and supercapacitors. We also outline both the bright prospects and emerging challenges these techniques face for application to graphene-based energy materials in future.

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces

Molecular dynamics simulations of the aqueous biomolecule–graphene interface have predicted the free energy of adsorption of amino acids and the structure of peptides.

Efficient simulations of the aqueous bio-interface of graphitic nanostructures with a polarisable model

To fully harness the enormous potential offered by interfaces between graphitic nanostructures and biomolecules, detailed connections between adsorbed conformations and adsorption behaviour are needed. To elucidate these links, a key approach, in partnership with experimental techniques, is molecular simulation. For this, a force-field (FF) that can appropriately capture the relevant physics and chemistry of these complex bio-interfaces, while allowing extensive conformational sampling, and also supporting inter-operability with known biological FFs, is a pivotal requirement. Here, we present and apply such a force-field, GRAPPA, designed to work with the CHARMM FF. GRAPPA is an efficiently implemented polarisable force-field, informed by extensive plane-wave DFT calculations using the revPBE-vdW-DF functional. GRAPPA adequately recovers the spatial and orientational structuring of the aqueous interface of graphene and carbon nanotubes, compared with more sophisticated approaches. We apply GRAPPA to determine the free energy of adsorption for a range of amino acids, identifying Trp, Tyr and Arg to have the strongest binding affinity and Asp to be a weak binder. The GRAPPA FF can be readily incorporated into mainstream simulation packages, and will enable large-scale polarisable biointerfacial simulations at graphitic interfaces, that will aid the development of biomolecule-mediated, solution-based graphene processing and self-assembly strategies.

Effective coarse-grained solid–fluid potentials and their application to model adsorption of fluids on heterogeneous surfaces

In the present contribution we emphasise the necessity of using an adequate averaging procedure to obtain effective fluid–surface potentials. A procedure to develop free-energy-averaged fluid–surface potentials retaining the important temperature dependence of the coarse-grained particle-surface interaction is described.

Wetting transition in water

Optical images were used to study the wetting behavior of water on graphite, sapphire, and quartz along the liquid vapor coexistence curve from room temperature to 300 °C. Wetting transitions were identified by the temperature at which the contact angle decreased to zero and also by the disappearance of dropwise condensation. These two methods yielded consistent values for the wetting temperatures, which were 185 °C, 234 °C, and 271 °C for water on quartz, sapphire, and graphite, respectively. We compare our results with the theoretical predictions based on a simplified model of the water-substrate potential and sharp interfaces.

Methane storage in molecular nanostructures

We survey various molecular structures which have been proposed as possible nanocontainers for methane storage. These are molecular structures that have been investigated through either experiments, molecular dynamics simulations or mathematical modelling. Computational simulation and mathematical modelling play an important role in predicting and verifying experimental outcomes, but both have their limitations. Even though recent advances have greatly improved computations, due to the large number of atoms and force field calculations involved, computational simulations can still be time consuming as compared to an instantaneous mathematical modelling approach. On the other hand, underlying an ideal mathematical model, there are many assumptions and approximations, but such modelling often reveals the key physical parameters and optimal configurations. Here, we review methane adsorption for three conventional nanostructures, namely graphite, single and multi-walled carbon nanotubes, and nanotube bundles (including interstitial and groove sites), and we survey methane adsorption in other molecular structures including metal organic frameworks. We also include an examination of minimum binding energies, equilibrium distances, gravimetric and volumetric uptakes, volume available for adsorption, as well as the effects of temperature and pressure on the adsorption of methane onto these molecular structures.