Performance of van der Waals Corrected Functionals for Guest Adsorption in the M2(dobdc) Metal-Organic Frameworks

Multilevel quantum mechanical calculations show the role of promoter molecules in the dehydration of methanol to dimethyl ether in H-ZSM-5

Methyl carboxylate esters promote the formation of dimethyl ether (DME) from the dehydration of methanol in H-ZSM-5 zeolite. We employ a multilevel quantum method to explore the possible associative and dissociative mechanisms in the presence, and absence, of six methyl ester promoters. This hybrid method combines density functional theory, with dispersion corrections (DFT-D3), for the full periodic system, with second-order Møller-Plesset perturbation theory (MP2) for small clusters representing the reaction site, and coupled cluster with single, double, and perturbative triple substitution (CCSD(T)) for the reacting molecules. The calculated adsorption enthalpy of methanol, and reaction enthalpies of the dehydration of methanol to DME within H-ZSM-5, agree with experiment to within chemical accuracy (∼4 kJ mol-1). For the promoters, a reaction pathway via an associative mechanism gives lower overall reaction enthalpies and barriers compared to the reaction with methanol only. Each stage of this mechanism is explored and related to experimental data. We provide evidence that suggests the promoter's adsorption to the Brønsted acid site is the most important factor dictating its efficiency.

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4'-dioxidobiphenyl-3,3'-dicarboxylate)-henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn-is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew-Burke-Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential's (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

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.

Understanding Carbon Monoxide Binding and Interactions in M-MOF-74 (M = Mg, Mn, Ni, Zn)

Small molecules may adsorb strongly in metal-organic frameworks (MOFs) through interactions with under-coordinated open metal sites (OMS) that often exist within these structures. Among adsorbates, CO is attractive to study both for its relevance in energy-related applications and for its ability to engage in both σ-donation and π-backbonding interactions with the OMS in MOFs. Concomitant with strong adsorption, structural changes arise due to modifications of the electronic structure of both the adsorbate and adsorbent. These structural changes affect the separation performance of materials, and accurately capturing these changes and the resulting energetics is critical for accurate predictive modeling of adsorption. Traditional approaches to modeling using classical force fields typically do not capture or account for changes at the electronic level. To characterize the structural and energetic effects of the local structural changes, we employed density functional theory (DFT) to study CO adsorption in M-MOF-74s. M-MOF-74s feature OMS at which CO is known to adsorb strongly and can be synthesized with a variety of divalent metal cations with distinct performance in adsorption. We considered M-MOF-74s with a range of metals of varied d-band occupations (Mg (3d0), Mn (3d5), Ni (3d8), and Zn (3d10)) with various structural constraints ranging from geometrically constrained adsorbent and adsorbate ions to fully optimized geometries to deconvolute the relative contributions of various structural effects to the adsorption energetics and binding distances observed. Our data indicate that the most significant structural changes during adsorption correlate with the greatest π-backbonding behaviors and commensurately result in a sizable binding energy change observed for CO adsorption. The insights built from this work are relevant to two longstanding research challenges within the MOF community: rational design of materials for separations and the design of force fields capable of accurately modeling adsorption.

Characterizing Grain Boundary Effects on Mg2+ Conduction in Metal-Organic Frameworks

Next-generation materials for fast ion conduction have the potential to revolutionize battery technology. Metal-organic frameworks (MOFs) are promising candidates for achieving this goal. Given their structural diversity, the design of efficient MOF-based conductors can be accelerated by a detailed understanding and accurate prediction of ion conductivity. However, the polycrystalline nature of solid-state materials requires consideration of grain boundary effects, which is complicated by challenges in characterizing grain boundary structures and simulating ensemble transport processes. To address this, we have developed an approach for modeling ion transport at grain boundaries and predicting their contribution to conductivity. Mg²⁺ conduction in the Mg-MOF-74 thin film was studied as a representative system. Using computational techniques and guided by experiments, we investigated the structural details of MOF grain boundary interfaces to determine accessible Mg²⁺ transport pathways. Computed transport kinetics were input into a simplified MOF nanocrystal model, which combines ion transport in the bulk structure and at grain boundaries. The model predicts Mg²⁺ conductivity in the MOF-74 film within chemical accuracy (<1 kcal/mol activation energy difference), validating our approach. Physically, Mg²⁺ conduction in MOF-74 is inhibited by strong Mg²⁺ binding at grain boundaries, such that only a small fraction of grain boundary alignments allow for fast Mg²⁺ transport. This results in a 2-3 order-of-magnitude reduction in conductivity, illustrating the critical impact of the grain boundary contribution. Overall, our work provides a computation-aided platform for molecular-level understanding of grain boundary effects and quantitative prediction of ion conductivity. Combined with experimental measurements, it can serve as a synergistic tool for characterizing the grain boundary composition of MOF-based conductors.

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.

Probing growth of metal-organic frameworks with X-ray scattering and vibrational spectroscopy

Nucleation and crystallization arising from liquid to solid phase are involved in a multitude of processes in fields ranging from materials science to biology. Controlling the thermodynamics and kinetics of growth is advantageous to help tune the formation of complex morphologies. Here, we harness wide-angle X-ray scattering and vibrational spectroscopy to elucidate the mechanism for crystallization and growth of the metal-organic framework Co-MOF-74 within microscopic volumes enclosed in a capillary and an attenuated total reflection microchip reactor. The experiments reveal molecular and structural details of the growth processes, while the results of plane wave density functional calculations allow identification of lattice and linker modes in the formed crystals. Synthesis of the metal-organic framework with microscopic volumes leads to monodisperse and micron-sized crystals, in contrast to those typically observed under bulk reaction conditions. Reduction in the volume of reagents within the microchip reactor was found to accelerate the reaction rate. The coupling of spectroscopy with scattering to probe reactions in microscopic volumes promises to be a useful tool in the synthetic chemist's kit to understand chemical bonding and has potential in designing complex materials.

Exploring mechanistic routes for light alkane oxidation with an iron-triazolate metal-organic framework

In this work, we computationally explore the formation and subsequent reactivity of various iron-oxo species in the iron-triazolate framework Fe₂(μ-OH)₂(bbta) (H₂bbta = 1H,5H-benzo(1,2-d:4,5-d')bistriazole) for the catalytic activation of strong C-H bonds. With the direct conversion of methane to methanol as the probe reaction of interest, we use density functional theory (DFT) calculations to evaluate multiple mechanistic pathways in the presence of either N₂O or H₂O₂ oxidants. These calculations reveal that a wide range of transition metal-oxo sites - both terminal and bridging - are plausible in this family of metal-organic frameworks, making it a unique platform for comparing the electronic structure and reactivity of different proposed active site motifs. Based on the DFT calculations, we predict that Fe₂(μ-OH)₂(bbta) would exhibit a relatively low barrier for N₂O activation and energetically favorable formation of an [Fe(O)]²⁺ species that is capable of oxidizing C-H bonds. In contrast, the use of H₂O₂ as the oxidant is predicted to yield an assortment of bridging iron-oxo sites that are less reactive. We also find that abstracting oxo ligands can exhibit a complex mixture of both positive and negative spin density, which may have broader implications for relating the degree of radical character to catalytic activity. In general, we consider the coordinatively unsaturated iron sites to be promising for oxidation catalysis, and we provide several recommendations on how to further tune the catalytic properties of this family of metal-triazolate frameworks.

RbSnX3 (X = Cl, Br, I): promising lead-free metal halide perovskites for photovoltaics and optoelectronics

Lead (Pb) free metal halide perovskites by atomistic design are of strong interest to photovoltaics and optoelectronics industries because of the pressing need to resolve Pb-related toxicity and instability challenges. In this study, structural, mechanical, electronic, and optical properties of Pb-free RbSnX₃ (X = Cl, Br, I) perovskites have been evaluated by using ab initio density functional theory (DFT) calculations. The computed elastic constants suggest that the Rb-based halide perovskites are mechanically stable and highly ductile, making them suitable as flexible thin films in optoelectronic devices. Besides, the investigated electronic band structures reveal that the RbSnX₃ compounds are direct bandgap semiconductors, suitable for photovoltaic and optoelectronic applications. Furthermore, several optical parameters such as dielectric functions, reflectivity, photon absorptions, refractive index, optical conductivity, and loss functions have been investigated and the results predict the excellent optoelectronic efficiency of RbSnX₃. Also, the computed mechanical and optical properties of RbSnX₃ (X = Cl, Br, I) have been compared with the previously studied CsBX₃ (B = Ge, Sn, Pb; X = Cl, Br, I) phases, revealing that the Rb-based perovskites are extremely ductile and possess excellent light absorption and optical conductivity compared to the Cs-based perovskites. Importantly, RbSnI₃ shows superior ductility, absorption coefficient, and optical conductivity compared to the CsBX₃ (B = Ge, Sn, Pb; X = Cl, Br, I) perovskites. Superior absorption at the ultraviolet region of RbSnI₃ holds great promise of this perovskite to be used in next-generation ultraviolet photodetectors.

This work summarizes that RbSnX₃ (X = Cl, Br, I) exhibits remarkable ductility and absorption in the ultraviolet (UV) region of the electromagnetic spectrum compared to those of CsBX₃ (B = Ge, Sn, Pb; X = Cl, Br, I) metal halide perovskites.

Catalytic properties of the ferryl ion in the solid state: a computational review

This review summarises the last findings in the emerging field of heterogeneous catalytic oxidation of light alkanes by ferryl species supported on solid-state systems such as the conversion of methane into methanol by FeO-MOF74.

Parametrization of Nonbonded Force Field Terms for Metal-Organic Frameworks Using Machine Learning Approach

The enormous structural and chemical diversity of metal-organic frameworks (MOFs) forces researchers to actively use simulation techniques as often as experiments. MOFs are widely known for their outstanding adsorption properties, so a precise description of the host-guest interactions is essential for high-throughput screening aimed at ranking the most promising candidates. However, highly accurate ab initio calculations cannot be routinely applied to model thousands of structures due to the demanding computational costs. Furthermore, methods based on force field (FF) parametrization suffer from low transferability. To resolve this accuracy-efficiency dilemma, we applied a machine learning (ML) approach: extreme gradient boosting. The trained models reproduced the atom-in-material quantities, including partial charges, polarizabilities, dispersion coefficients, quantum Drude oscillator, and electron cloud parameters, with accuracy similar to the reference data set. The aforementioned FF precursors make it possible to thoroughly describe noncovalent interactions typical for MOF-adsorbate systems: electrostatic, dispersion, polarization, and short-range repulsion. The presented approach can also readily facilitate hybrid atomistic simulation/ML workflows.

Reinvestigating oxygen adsorption on Ag(111) by using strongly constrained and appropriately normed semi-local density functional with the revised Vydrov van Voorhis van der Waals force correction

While density functional theory (DFT) at the generalized gradient approximation (GGA) level has made great success in catalysis, it fails in some important systems such as the adsorption of the oxygen molecule on the Ag(111) surface. Previous DFT studies at the GGA level revealed theoretical inconsistencies on the adsorption energies and dissociation barriers of O₂ on Ag(111) in comparison with the experimental conclusion. In this study, the strongly constrained and appropriately normed-revised Vydrov van Voorhis van der Waals correction functional (SCAN-rVV10) method at the meta-GGA level with the nonlocal van der Waals (vdW) force correction was used to reinvestigate the adsorption properties of O₂ on the Ag(111) surface. The SCAN-rVV10 results successfully confirm the experimental observation that both molecular and dissociative adsorptions can exist for oxygen on Ag(111). The calculated adsorption energy for the physisorption state and the relevant dissociation energy barrier are close to the experimental data. It demonstrates that SCAN-rVV10 can outperform functionals at the GGA level for O₂/Ag(111). Therefore, our findings suggest that SCAN-rVV10 can be the desired method for systems where the correct description of intermediate-ranged vdW forces is essential, such as the physisorption of small molecules on the solid surface.

Molecular-Level Insights into Selective Transport of Mg2+ in Metal-Organic Frameworks

Metal-organic frameworks (MOF) are promising media for achieving solid-state Mg2+ conduction and developing a magnesium-based battery. To this end, the chemical behavior and transport properties of an Mg(TFSI)2/DME electrolyte system inside Mg-MOF-74 were studied by density functional theory (DFT). We found that inside the MOF chemical environment, solvent and anion molecules occupy the coordinatively unsaturated open metal sites of Mg-MOF-74, while Mg2+ ions adsorb directly onto the carboxylate group of the MOF organic linker. These predicted binding geometries were further corroborated by IR spectroscopy. We computed the free energies of desolvation of Mg2+ ions inside MOF to investigate the capacity of Mg-MOF-74 thin film to act as a separator for selective Mg2+ transport. We showed that Mg-MOF-74 could facilitate partial, but not full, desolvation of Mg2+. We found that the dominant minimum-energy pathway (MEP) for Mg2+ conduction inside Mg-MOF-74 corresponds to a "solvent hopping" mechanism, with an energy barrier of 4.4 kcal/mol. The molar conductivity of Mg2+ associated with the idealized solvent hopping mechanism along the MOF one-dimensional channel was predicted to be 2.4 × 10-3 S cm-1 M-1, which is one to two orders of magnitude greater than the experimentally measured value of 1.2 × 10-4 S cm-1 M-1 (with an estimated Mg2+ concentration). We have discussed several possible factors contributing to this apparent discrepancy. The current work demonstrates the validity of the computational strategies applied and the structural models constructed for the understanding of fast and selective Mg2+ transport in Mg-MOF-74, which serves as a cornerstone for studying transport of multivalent ions in MOFs. Furthermore, it provides detailed molecular-level insights that are not yet accessible experimentally.

vdW-DF-ahcx: a range-separated van der Waals density functional hybrid

Hybrid density functionals replace a fraction of an underlying generalized-gradient approximation (GGA) exchange description with a Fock-exchange component. Range-separated hybrids (RSHs) also effectively screen the Fock-exchange component and thus open the door for characterizations of metals and adsorption at metal surfaces. The RSHs are traditionally based on a robust GGA, such as PBE (Perdew J Pet al1996Phys. Rev. Lett.773865), for example, as implemented in the HSE design (Heyd Jet al2003J. Chem. Phys.1188207). Here we define an analytical-hole (Henderson T Met al2008J. Chem. Phys.128194105) consistent-exchange RSH extension to the van der Waals density functional (vdW-DF) method (Berland Ket al2015Rep. Prog. Phys.78066501), launching vdW-DF-ahcx. We characterize the GGA-type exchange in the vdW-DF-cx version (Berland K and Hyldgaard P 2014Phys. Rev. B89035412), isolate the short-ranged exchange component, and define the new vdW-DF hybrid. We find that the performance vdW-DF-ahcx compares favorably to (dispersion-corrected) HSE for descriptions of bulk (broad molecular) properties. We also find that it provides accurate descriptions of noble-metal surface properties, including CO adsorption.

A Microporous Hydrogen-Bonded Organic Framework for the Efficient Capture and Purification of Propylene

Adsorptive separation of propylene/propane (C₃ H₆ /C₃ H₈ ) mixture is desired for its potential energy saving on replacing currently deployed and energy-intensive cryogenic distillation. Realizing efficient C₃ H₆ /C₃ H₈ separation in the emerging hydrogen-bonded organic frameworks (HOFs) is very challenging owing to the lack of functional sites for preferential gas binding. By virtue of crystal engineering, we herein report a functionalized HOF (HOF-16) with free -COOH sites for the efficient separation of C₃ H₆ /C₃ H₈ mixtures. Under ambient conditions, HOF-16 shows a significant C₃ H₆ /C₃ H₈ uptake difference (by 76 %) and selectivity (5.4) in contrast to other carboxylic acid-based HOFs. Modeling studies indicate that free -COOH groups together with the suitable pore confinement facilitate the recognition and high-density packing of gas molecules. The separation performance of HOF-16 was validated by breakthrough experiments. HOF-16 is stable towards strong acidity and water.

Defects as catalytic sites for the oxygen evolution reaction in Earth-abundant MOF-74 revealed by DFT

The oxygen evolution reaction (OER) is the bottleneck of hydrogen production via water splitting and understanding electrocatalysts at atomic level becomes paramount to enhance the efficiency of this process.

Theoretical understanding of the properties of stepped iron surfaces with van der Waals interaction corrections

The stepped surfaces in nanoscale zero-valent iron (nZVI) play an essential role for environmental application. However, there is still currently a deficiency in the atomic understanding of stepped surface properties due to the limitation of the computational methodology. In this study, stepped Fe(210) and (211) surfaces were theoretically investigated using density functional theory (DFT) computations in terms of the flat Fe(110) surface. Our results suggest that the consideration of van der Waals (vdW) interaction correction is beneficial for the DFT study on Fe-based systems. The DF-cx method is found to be the most promising vdW correction method. The DF-cx results reveal that the stepped Fe(210) and Fe(211) surfaces experience significant surface relaxation and abnormal trends in their work function. Their electronic properties and reactivities of the surface atoms are strongly affected by the Fe coordination numbers and the strong adsorption strengths of oxygen on the surfaces are dependent on both the coordination number of the adsorbed atoms and the geometry of the adsorption sites.

Impact of Small Adsorbates in the Vibrational Spectra of Mg- and Zn-MOF-74 Revealed by First-Principles Calculations

In this work, we analyze the influence of small adsorbates on the vibrational spectra of Mg- and Zn-metal-organic framework MOF-74 by means of first-principles calculations. In particular, we consider the adsorption of four representative species of different interaction strengths: Ar, CO₂, H₂O, and NH₃. Apart from a comprehensive characterization of the structural and energetic aspects of empty and loaded MOFs, we use a fully quantum ab initio approach to evaluate the Raman and IR activities of the normal modes, leading to the construction of the whole vibrational spectra. Under this approach, not only are we able to proceed with the complete assignment of the spectra in terms of the usual internal coordinates but also we can discern the most relevant vibrational fingerprints of the adsorbates and their impact on the whole MOF spectra. On the one hand, some of the typical vibrational modes of the small molecules are slightly shifted but still visible when adsorbed on the MOFs, especially those appearing at high wavenumbers where the empty MOFs lack IR/Raman signals. On the other hand, some bands arising from the organic ligands are affected by the presence of the absorbates, displaying non-negligible frequency shifts, in agreement with recent experiments. We find a strong correlation between all of these frequency shifts and the interaction strength of the adsorbate with the hosting framework. The findings presented in this work expand the capabilities of vibrational spectroscopy techniques to analyze porous materials and can be useful for the design of sensors and new devices based on MOF technology.

Intra- and Interchain Interactions in (Cu1/2Au1/2)CN, (Ag1/2Au1/2)CN, and (Cu1/3Ag1/3Au1/3)CN and Their Effect on One-, Two-, and Three-Dimensional Order

Mixed-metal cyanides (Cu1/2Au1/2)CN, (Ag1/2Au1/2)CN, and (Cu1/3Ag1/3Au1/3)CN adopt an AuCN-type structure in which metal-cyanide chains pack on a hexagonal lattice with metal atoms arranged in sheets. The interactions between and within the metal-cyanide chains are investigated using density functional theory (DFT) calculations, 13C solid-state NMR (SSNMR), and X-ray pair distribution function (PDF) measurements. Long-range metal and cyanide order is found within the chains: (-Cu-NC-Au-CN-)∞, (-Ag-NC-Au-CN-)∞, and (-Cu-NC-Ag-NC-Au-CN-)∞. Although Bragg diffraction studies establish that there is no long-range order between chains, X-ray PDF results show that there is local order between chains. In (Cu1/2Au1/2)CN and (Ag1/2Au1/2)CN, there is a preference for unlike metal atoms occurring as nearest neighbors within the metal sheets. A general mathematical proof shows that the maximum average number of heterometallic nearest-neighbor interactions on a hexagonal lattice with two types of metal atoms is four. Calculated energies of periodic structural models show that those with four unlike nearest neighbors are most favorable. Of these, models in space group Immm give the best fits to the X-ray PDF data out to 8 Å, providing good descriptions of the short- and medium-range structures. This result shows that interactions beyond those of nearest neighbors must be considered when determining the structures of these materials. Such interactions are also important in (Cu1/3Ag1/3Au1/3)CN, leading to the adoption of a structure in Pmm2 containing mixed Cu-Au and Ag-only sheets arranged to maximize the numbers of Cu···Au nearest- and next-nearest-neighbor interactions.

Coordinated Water as New Binding Sites for the Separation of Light Hydrocarbons in Metal-Organic Frameworks with Open Metal Sites

Metal-organic frameworks with open metal sites are promising materials for gas separations. Particularly, the M₂(dobdc) (dobdc^4- = 2,5-dioxidobenzenedicarboxylate, M²⁺ = Co²⁺, Mn²⁺, Fe²⁺, ...) framework has been the Drosophila of this research field and has delivered groundbreaking results in terms of sorption selectivity. However, many studies focus on perfect two-component mixtures and use theoretical models, e.g., the ideal adsorbed solution theory, to calculate selectivities. Within this work, we shed light on the comparability of these selectivities with values obtained from propane/propene multicomponent measurements on the prototypical Co₂(dobdc) framework, and we study the impact of impurities like water on the selectivity. Despite the expected capacity loss, the presence of water does not necessarily lead to a decreased selectivity. Density functional theory calculations of the binding energies prove that the water molecules adsorbed to the metal centers introduce new binding sites for the adsorbates.

Coordinatively unsaturated metal sites (open metal sites) in metal–organic frameworks: design and applications

The defined synthesis of OMS in MOFs is the basis for targeted functionalization through grafting, the coordination of weakly binding species and increased (supramolecular) interactions with guest molecules.

Disparate strain response of the thermal transport properties of bilayer penta-graphene as compared to that of monolayer penta-graphene

In this study, strain modulation of the lattice thermal conductivity of monolayer and bilayer penta-graphene (PG) at room temperature was investigated using first-principles calculations combined with the phonon Boltzmann transport equation.

Ab Initio Evaluation of Henry Coefficients Using Importance Sampling

We present a new algorithm that allows for an efficient evaluation of the Henry coefficient of a guest molecule inside a porous material, which permits to use ab initio energy calculations. The Widom insertion method, which is currently used to compute these Henry coefficients, typically requires millions of energy evaluations. Our new methodology reduces this number by more than 1 order of magnitude, enabling the use of an ab initio potential energy surface. The methodology we propose is reminiscent of the well-known importance sampling technique which is frequently used in Monte Carlo integrations. First, a conventional Widom insertion simulation is performed using a force field. In the second step, the Widom results are used to select a limited number of configurations and only for these configurations the ab initio evaluation of the energy is required. Finally, by appropriately reweighting the latter energies, an accurate estimation of the ab initio Henry coefficient is possible at a moderate computational cost. We apply our methodology to the adsorption of CO₂ in Mg-MOF-74, a prototypical case where interactions of a polar guest molecule with unsaturated metal sites dominate the adsorption mechanism. In this case generic force fields such as UFF or Dreiding are inappropriate and the use of ab initio methods is indispensable. In a second case study, we compute Henry coefficients of methane in UiO-66 using different levels of theory. We pay particular attention to the influence of the dispersion corrections and the role of many-body effects. For the final example, we qualitatively investigate adsorption features for a series of functionalized UiO-66 frameworks. Overall the cases we present show that accurate computations of Henry coefficients is extremely challenging, as different levels of theory provide strongly varying results. At the same time ab initio calculations have added value compared to force fields, as they provide a physically more sound description of the adsorption mechanism and in some cases clearly improve correspondence with experiment.

Enhancement of CO2 binding and mechanical properties upon diamine functionalization of M2(dobpdc) metal-organic frameworks

The family of diamine-appended metal-organic frameworks exemplified by compounds of the type mmen-M2(dobpdc) (mmen = N,N'-dimethylethylenediamine; M = Mg, Mn, Fe, Co, Zn; dobpdc4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) are adsorbents with significant potential for carbon capture, due to their high working capacities and strong selectivity for CO2 that stem from a cooperative adsorption mechanism. Herein, we use first-principles density functional theory (DFT) calculations to quantitatively investigate the role of mmen ligands in dictating the framework properties. Our van der Waals-corrected DFT calculations indicate that electrostatic interactions between ammonium carbamate units significantly enhance the CO2 binding strength relative to the unfunctionalized frameworks. Additionally, our computed energetics show that mmen-M2(dobpdc) materials can selectively adsorb CO2 under humid conditions, in agreement with experimental observations. The calculations further predict an increase of 112% and 124% in the orientationally-averaged Young's modulus E and shear modulus G, respectively, for mmen-Zn2(dobpdc) compared to Zn2(dobpdc), revealing a dramatic enhancement of mechanical properties associated with diamine functionalization. Taken together, our calculations demonstrate how functionalization with mmen ligands can enhance framework gas adsorption and mechanical properties.

CO2 adsorption on the copper surfaces: van der Waals density functional and TPD studies

We investigated the adsorption of CO₂ on the flat, stepped, and kinked copper surfaces from density functional theory calculations as well as the temperature programmed desorption and X-ray photoelectron spectroscopy. Several exchange-correlation functionals have been considered to characterize CO₂ adsorption on the copper surfaces. We used the van der Waals density functionals (vdW-DFs), i.e., the original vdW-DF (vdW-DF1), optB86b-vdW, and rev-vdW-DF2, as well as the Perdew-Burke-Ernzerhof (PBE) with dispersion correction (PBE-D2). We have found that vdW-DF1 and rev-vdW-DF2 functionals slightly underestimate the adsorption energy, while PBE-D2 and optB86b-vdW functionals give better agreement with the experimental estimation for CO₂ on Cu(111). The calculated CO₂ adsorption energies on the flat, stepped, and kinked Cu surfaces are 20-27 kJ/mol, which are compatible with the general notion of physisorbed species on solid surfaces. Our results provide a useful insight into appropriate vdW functionals for further investigation of related CO₂ activation on Cu surfaces such as methanol synthesis and higher alcohol production.

First principles Monte Carlo simulations of unary and binary adsorption: CO2, N2, and H2O in Mg-MOF-74

Dative bonding of adsorbate molecules onto coordinatively-unsaturated metal sites in metal–organic frameworks (MOFs) can lead to unique adsorption selectivities.