Searching for DFT-based methods that include dispersion interactions to calculate the physisorption of H2 on benzene and graphene

Studying the adsorption of emerging organic contaminants in zeolites with dispersion-corrected density functional theory calculations: From numbers to recommendations

Adsorption energies obtained from dispersion-corrected density functional theory (DFT) calculations show a considerable dependence on the choice of exchange-correlation functional and dispersion correction. A number of investigations have employed different approaches to compute adsorption energies of small molecules in zeolites, using reference values from high-level calculations and/or experiments. Such comparative studies are lacking for larger functional organic molecules such as pharmaceuticals or personal care products, despite their potential relevance for applications, e. g., in contaminant removal or drug delivery. The present study aims to fill this gap by comparing adsorption energies and, for selected cases, equilibrium structures of emerging organic contaminants adsorbed in MOR- and FAU-type all-silica zeolites. A total of 13 dispersion-corrected DFT approaches are compared, including methods using a pairwise dispersion correction as well as non-local van der Waals density functionals. While absolute values of adsorption energies vary widely, qualitative trends across the set of zeolite-guest combinations are not strongly dependent on the choice of functional. For selected cluster models, DFT adsorption energies are compared to reference values from coupled cluster (DLPNO-CCSD(T)) calculations. Although all DFT approaches deliver systematically more negative adsorption energies than the coupled cluster reference, this tendency is least pronounced for the rev-vdW-DF2 functional.

Density Functional Theory-Based Approaches to Improving Hydrogen Storage in Graphene-Based Materials

Various technologies have been developed for the safe and efficient storage of hydrogen. Hydrogen storage in its solid form is an attractive option to overcome challenges such as storage and cost. Specifically, hydrogen storage in carbon-based structures is a good solution. To date, numerous theoretical studies have explored hydrogen storage in different carbon structures. Consequently, in this review, density functional theory (DFT) studies on hydrogen storage in graphene-based structures are examined in detail. Different modifications of graphene structures to improve their hydrogen storage properties are comprehensively reviewed. To date, various modified graphene structures, such as decorated graphene, doped graphene, graphene with vacancies, graphene with vacancies-doping, as well as decorated-doped graphene, have been explored to modify the reactivity of pristine graphene. Most of these modified graphene structures are good candidates for hydrogen storage. The DFT-based theoretical studies analyzed in this review should motivate experimental groups to experimentally validate the theoretical predictions as many modified graphene systems are shown to be good candidates for hydrogen storage.

Adsorption of Carbamazepine in All-Silica Zeolites Studied with Density Functional Theory Calculations

The anticonvulsant drug carbamazepine (-) is an emerging contaminant of considerable concern due to its hazard potential and environmental persistence. Previous experimental studies proposed hydrophobic zeolites as promising adsorbents for the removal of carbamazepine from water, but only a few framework types were considered in those investigations. In the present work, electronic structure calculations based on dispersion-corrected density functional theory (DFT) were used to study the adsorption of CBZ in eleven all-silica zeolites having different pore sizes and connectivities of the pore system (AFI, ATS, BEA, CFI, DON, FAU, IFR, ISV, MOR, SFH, SSF framework types). It was found that some zeolites with one-dimensional channels formed by twelve-membered rings (IFR, AFI) exhibit the highest affinity towards CBZ. A "good fit" of CBZ into the zeolite pores that maximizes dispersion interactions was identified as the dominant factor determining the interaction strength. Further calculations addressed the role of temperature (for selected systems) and of guest-guest interactions between coadsorbed CBZ molecules. In addition to predicting zeolite frameworks of particular interest as materials for selective CBZ removal, the calculations presented here also contribute to the atomic-level understanding of the interaction of functional organic molecules with all-silica zeolites.

Assessment of density functional approximations for N2 and CO2 physisorption on benzene and graphene

Experimental isotherms of N₂ and CO₂ on carbon‐based porous materials and models of the physisorption of gases on surfaces are used to obtain the pore size distribution (PSD). An accurate modelization of the physisorption of N₂ and CO₂ on the surface of carbon‐based porous materials is important to obtain accurate N₂ and CO₂ storage capacities and reliable PSDs. Physisorption depends on the dispersion interactions. High precision ab initio methods, such as CCSD(T), consider accurately the dispersion interactions, but they are computationally expensive. Double hybrid, hybrid and DFT‐based methods are much less expensive. In the case of graphene, there are experimental data of the adsorption of N₂ and CO₂ on graphite that can be used to build the Steele interaction potential of these gases on graphene. The goal is to find out hybrid and/or DFT methods that are as accurate as the CCSD(T) on benzene and as accurate as the experimental results on graphene. Calculations of the interaction energy curves of N₂ and CO₂ on benzene and graphene have been carried out using the CCSD(T) method and several double hybrid, hybrid, and DFT methods that consider the dispersion interactions. The energy curves on benzene have been compared to the CCSD(T) and the energy curves on graphene have been compared with the Steele energy curves. The comparisons indicate that double hybrids with dispersion corrections and ωB97 based DFT methods are accurate enough for benzene. For graphene, only the PBE‐XDM functional has a good agreement with the Steele energy curves.

Physi-Sorption of H2 on Pure and Boron–Doped Graphene Monolayers: A Dispersion–Corrected DFT Study

High-surface-area carbons are of interest as potential candidates to store H2 for fuel–cell power applications. Earlier work has been ambiguous and inconclusive on the effect of boron doping on H2 binding energy. Here, we describe a systematic dispersion–corrected density functional theory study to evaluate the effect of boron doping. We observe some enhancement in H2 binding, due to the presence of a defect, such as terminal hydrogen or distortion from planarity, introduced by the inclusion of boron into a graphene ring, which creates hydrogen adsorption sites with slightly increased binding energy. The increase is from −5 kJ/mol H2 for the pure carbon matrix to −7 kJ/mol H2 for the boron–doped system with the boron content of ~7%. The H2 binding sites have little direct interaction with boron. However, the largest enhancement in physi-sorption energy is seen for systems, where H2 is confined between layers at a distance of about 7 Å, where the H2 binding nearly doubles to −11 kJ/mol H2. These findings suggest that interplanar nanoconfinement might be more effective in enhancing H2 binding. Smaller coronene model is shown to be beneficial for understanding the dependence of interaction energy on the structural configurations and preferential H2 binding sites.

Low-temperature hydrogen-graphite system revisited: Experimental study and Monte Carlo simulation

Hydrogen adsorption by microporous carbon materials attracts much attention for the past few decades, which has been stimulated by growing interest in hydrogen storage. Numerous studies of this topic based on molecular simulation technique have been reported. However, in many cases, the reliability of the results obtained with numerical methods is insufficient, which is a consequence of poor reference data used for fitting parameters of the fluid-fluid and fluid-solid potentials. This study is devoted to a detailed experimental investigation of the hydrogen-graphite system and its modeling with a kinetic Monte Carlo method at temperatures from 20 to 77 K and the bulk pressure from 0.1 Pa to 100 kPa. We found that the best fit of the bulk hydrogen equation of state corresponds to the 10-6 Lennard-Jones potential with the temperature dependent parameters to account for the quantum effects. The experimental hydrogen adsorption isotherms on the graphite surface were fitted with a high accuracy, which constitutes a firm basis of subsequent simulation of hydrogen adsorption in various nanoporous carbons and their pore size distribution analysis using a kernel generated with the developed approach.

Enhanced hydrogen storage performance of graphene nanoflakes doped with Cr atoms: a DFT study

The hydrogen storage performances of novel graphene nanoflakes doped with Cr atoms were systematically investigated using first-principles density functional theory. The calculated results showed that one Cr atom could be successfully doped into the graphene nanoflake with a binding energy of −4.402 eV. Different from the H₂ molecule moving away from the pristine graphene nanoflake surface, the built Cr-doped graphene nanoflake exhibited a high affinity to the H₂ molecule with a chemical adsorption energy of −0.574 eV. Moreover, the adsorptions of two to five H₂ molecules on the Cr-doped graphene nanoflake were studied as well. It was found that there were a maximum of three H₂ molecules stored on the graphene nanoflake doped with one Cr atom. Also, the further calculations showed that the numbers of the stored H₂ molecules were effectively improved to be six (or nine) when the graphene nanoflakes were doped with two (or three) Cr atoms. This research reveals that the graphene nanoflake doped with Cr atom could be a promising material to store H₂ molecules and its H₂ storage performance could be effectively enhanced through modifying the number of doped Cr atoms.

Our study reveals that the H₂ storage performance of a graphene nanoflake based material could be significantly enhanced through doping with Cr atoms.

Extent of Fock-exchange mixing for a hybrid van der Waals density functional?

The vdW-DF-cx0 exchange-correlation hybrid design [K. Berland et al., J. Chem. Phys. 146, 234106 (2017)] has a truly nonlocal correlation component and aims to facilitate concurrent descriptions of both covalent and non-covalent molecular interactions. The vdW-DF-cx0 design mixes a fixed ratio, a, of the Fock exchange into the consistent-exchange van der Waals density functional, vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)]. The mixing value a is sometimes taken as a semi-empirical parameter in hybrid formulations. Here, instead, we assert a plausible optimum average a value for the vdW-DF-cx0 design from a formal analysis; A new, independent determination of the mixing a is necessary since the Becke fit [A. D. Becke, J. Chem. Phys. 98, 5648 (1993)], yielding a' = 0.2, is restricted to semilocal correlation and does not reflect non-covalent interactions. To proceed, we adapt the so-called two-legged hybrid construction [K. Burke et al., Chem. Phys. Lett. 265, 115 (1997)] to a starting point in the vdW-DF-cx functional. For our approach, termed vdW-DF-tlh, we estimate the properties of the adiabatic-connection specification of the exact exchange-correlation functional, by combining calculations of the Fock exchange and of the coupling-constant variation in vdW-DF-cx. We find that such vdW-DF-tlh hybrid constructions yield accurate characterizations of molecular interactions (even if they lack self-consistency). The accuracy motivates trust in the vdW-DF-tlh determination of system-specific values of the Fock-exchange mixing. We find that an average value a' = 0.2 best characterizes the vdW-DF-tlh description of covalent and non-covalent interactions, although there exists some scatter. This finding suggests that the original Becke value, a' = 0.2, also represents an optimal average Fock-exchange mixing for the new, truly nonlocal-correlation hybrids. To enable self-consistent calculations, we furthermore define and test a zero-parameter hybrid functional vdW-DF-cx0p (having fixed mixing a' = 0.2) and document that this truly nonlocal correlation hybrid works for general molecular interactions (at reference and at relaxed geometries). It is encouraging that the vdW-DF-cx0p functional remains useful also for descriptions of some extended systems.