We use density functional theory calculations with van der Waals corrections to study the role of dispersive interactions on the structure and binding of CO 2 within two distinct metal−organic frameworks (MOFs): Mg-MOF74 and Ca-BTT. For both classes of MOFs, we report calculations with standard gradient-corrected (PBE) and five van der Waals density functionals (vdW-DFs), also comparing with semiempirical pairwise corrections. The vdW-DFs explored here yield a large spread in CO 2 −MOF binding energies, about 50% (around 20 kJ/ mol), depending on the choice of exchange functional, which is significantly larger than our computed zero-point energies and thermal contributions (around 5 kJ/mol). However, two specific vdW-DFs result in excellent agreement with experiments within a few kilojoules per mole, at a reduced computational cost compared to quantum chemistry or many-body approaches. For Mg-MOF74, PBE underestimates adsorption enthalpies by about 50%, but enthalpies computed with vdW-DF, PBE+D2, and vdW-DF2 (40.5, 38.5, and 37.4 kJ/mol, respectively) compare extremely well with the experimental value of 40 kJ/mol. vdW-DF and vdW-DF2 CO 2 −MOF bond lengths are in the best agreement with experiments, while vdW-C09 x results in the best agreement with lattice parameters. On the basis of the similar behavior of the reduced density gradients around CO 2 for the two MOFs studied, comparable results can be expected for CO 2 adsorption in BTT-type MOFs. Our work demonstrates for this broad class of molecular adsorbate-periodic MOF systems that parameter-free and computationally efficient vdW-DF and vdW-DF2 approaches can predict adsorption enthalpies with chemical accuracy.
■ INTRODUCTIONMetal−organic frameworks (MOFs) are a broad class of threedimensional nanoporous materials that have attracted much attention during the past decade for CO 2 capture from flue gas. 1−6 MOFs consist of metal centers joined by organic molecular "linkers"; there are many possibilities for cation/ linker combinations, and thousands of different MOFs have already been synthesized. 3,5 Understanding mechanisms for CO 2 binding within these frameworks is a fundamental step toward the design of new materials for carbon capture.In order to explore and predict new materials that efficiently capture CO 2 , an accurate quantitative description of both its adsorption energy and geometry is needed. Whereas quantum chemistry methods, such as HF/MP2 or CCSD(T), scale unfavorably with the number of basis functions (N 5 and N 7 ) and are primarily restricted to cluster calculations, density functional theory (DFT) is a very promising framework for studying the mechanism of interaction between CO 2 and extended MOFs. Although DFT is a many-particle framework that includes, in principle, fully nonlocal interactions, its common approximations, such as the local density approximation (LDA) and the generalized gradient approximation (GGA), neglect attractive long-range contributions to van der Waals interactions, so-called London dispersion intera...