Search

Zirconium-based metal–organic frameworks (Zr-MOFs) have been increasingly studied over the past two decades as heterogeneous catalysts due to their synthetic tunability, well-defined nature, and ch...

Through quantum-chemical calculations, we investigate a family of metal-organic frameworks (MOFs) containing triazolate linkers, M2 X2 (BBTA) (M=metal, X=bridging anion, H2 BBTA=1H,5H-benzo(1,2-d:4,5-d')bistriazole), for their ability to form terminal metal-oxo sites and subsequently activate the C-H bond of methane. By varying the metal and bridging anion in the framework, we show how to significantly tune the reactivity of this series of MOFs. The electronic structure of the metal-oxo active site is analyzed for each combination of metal and bridging ligand, and we find that spin density localized on the oxo ligand is not an inherent requirement for low C-H activation barriers. For the Mn- and Fe-containing frameworks, a transition from ferromagnetic to antiferromagnetic coupling between the metal binding site and terminal oxo ligand during the C-H activation process can greatly reduce the kinetic barrier, a unique case of two-state reactivity without a change in the net spin multiplicity.

Metal-organic frameworks (MOFs) with coordinatively unsaturated metal sites are appealing as adsorbent materials due to their tunable functionality and ability to selectively bind small molecules. Through the use of computational screening methods based on periodic density functional theory, we investigate O2 and N2 adsorption at the coordinatively unsaturated metal sites of several MOF families. A variety of design handles are identified that can be used to modify the redox activity of the metal centers, including changing the functionalization of the linkers (replacing oxido donors with sulfido donors), anion exchange of bridging ligands (considering μ-Br-, μ-Cl-, μ-F-, μ-SH-, or μ-OH- groups), and altering the formal oxidation state of the metal. As a result, we show that it is possible to tune the O2 affinity at the open metal sites of MOFs for applications involving the strong and/or selective binding of O2. In contrast with O2 adsorption, N2 adsorption at open metal sites is predicted to be relatively weak across the MOF dataset, with the exception of MOFs containing synthetically elusive V2+ open metal sites. As one example from the screening study, we predicted that exchanging the μ-Cl- ligands of M2Cl2(BBTA) (H2BBTA = 1H,5H-benzo(1,2-d:4,5-d')bistriazole) with μ-OH- groups would significantly enhance the strength of O2 adsorption at the open metal sites without a corresponding increase in the N2 affinity. Experimental investigation of Co2Cl2(BBTA) and Co2(OH)2(BBTA) confirms that the former exhibits weak physisorption of both N2 and O2, whereas the latter is capable of chemisorbing O2 at room temperature in a highly selective manner. The O2 chemisorption behavior is attributed to the greater electron-donating character of the μ-OH- ligands and the presence of H-bonding interactions between the μ-OH- bridging ligands and the reduced O2 adsorbate.

The development of adsorbents with molecular precision offers a promising strategy to enhance storage of hydrogen and methane─considered the fuel of the future and a transitional fuel, respectively─and to realize a carbon-neutral energy cycle. Herein we employ a postsynthetic modification strategy on a robust metal-organic framework (MOF), MFU-4l, to boost its storage capacity toward these clean energy gases. MFU-4l-Li displays one of the best volumetric deliverable hydrogen capacities of 50.2 g L-1 under combined temperature and pressure swing conditions (77 K/100 bar → 160 K/5 bar) while maintaining a moderately high gravimetric capacity of 9.4 wt %. Moreover, MFU-4l-Li demonstrates impressive methane storage performance with a 5-100 bar usable capacity of 251 cm3 (STP) cm-3 (0.38 g g-1) and 220 cm3 (STP) cm-3 (0.30 g g-1) at 270 and 296 K, respectively. Notably, these hydrogen and methane storage capacities are significantly improved compared to those of its isoreticular analogue, MFU-4l, and place MFU-4l-Li among the best MOF-based materials for this application.

In this work, we leverage advances in computational screening based on periodic density functional theory (DFT) to study a diverse set of experimentally derived metal–organic frameworks (MOFs) with...

Metal−organic frameworks (MOFs) with coordinatively unsaturated metal sites are appealing as adsorbent materials due to their tunable functionality and ability to selectively bind small molecules. Through the use of computational screening methods based on periodic density functional theory, we investigate O2 and N2 adsorption at the coordinatively unsaturated metal sites of several MOF families. A variety of design handles are identified that can be used to modify the redox activity of the metal centers, including changing the functionalization of the linkers (replacing oxido donors with sulfido donors), anion exchange of bridging ligands (considering μ-Br-, μ-Cl-, μ-F-, μ-SH-, or μ-OH- groups), and altering the formal oxidation state of the metal. As a result, we show that it is possible to tune the O2 affinity at the open metal sites of MOFs for applications involving the strong and/or selective binding of O2. In contrast with O2 adsorption, N2 adsorption at open metal sites is predicted to be relatively weak across the MOF dataset, with the exception of MOFs containing synthetically elusive V2+ open metal sites. As one example from the screening study, we predict that exchanging the μ-Cl- ligands of M2Cl2(BBTA) (H2BBTA = 1H,5H-benzo(1,2-d:4,5-d′)bistriazole) with μ-OH- groups would significantly enhance the strength of O2 adsorption at the open metal sites without a corresponding increase in the N2 affinity. Experimental investigation of Co2Cl2(BBTA) and Co2(OH)2(BBTA) confirms that the former exhibits only weak physisorption, whereas the latter is capable of chemisorbing O2 at room temperature. The chemisorption behavior is attributed to the greater electron-donating character of the μ-OH- ligands and the presence of H-bonding interactions between the μ-OH- bridging ligands and the O2 adsorbate.

Polymorphism in metal−organic frameworks (MOFs) means that the same chemical building blocks (nodes and linkers) can be used to construct isomeric MOFs with different topological networks. The choice of topology can substantially impact the pore network of the MOF, changing the sizes and shapes of the pores, which has implications for adsorption and separation applications. In this work, we look at the influence of topology in 38 polymorphic MOFs on the separation of linear and branched C4–C6 alkane isomers, a separation of great importance to the petrochemical industry. We find that the MOF Cu2(1,4-benzenedicarboxylate) in nbo topology (nbo-Cu2BDC) has particularly high affinity for linear alkanes due to its small pore size, which excludes the branched isomers. Upon studying this MOF in further detail, we find that it can take either of two conformations: a cubic conformation, which is typical of nbo MOFs, and a unique star conformation that contains 1D triangular and hexagonal channels. The determination of which conformation the MOF will adopt depends on steric effects between the nodes and linkers.

Metal-organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O2 and N2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal-adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure-property relationships.

Copper oxide catalysts supported on KIT-6 silica were evaluated for cyclohexane oxidative dehydrogenation (ODH) to determine the effects of copper oxide domain size on ODH activity and selectivity. The catalysts were prepared by incipient wetness impregnation of KIT-6 at copper surface densities spanning 0.01–0.7 Cu/nm2 with carefully controlled drying and calcination conditions to systematically vary the average local copper oxide domain size. A distinct copper oxide active site exhibiting an order of magnitude higher activity than large copper oxide domains was identified by model cyclohexane ODH studies coupled with in situ X-ray absorption and UV–visible spectroscopies during reduction in H2. The structure of this site is experimentally identified by a combination of extended X-ray absorption fine structure analysis, resonant Raman studies, and modeling by density functional theory. All constraints imposed by these techniques indicate the active site is a mono(μ-oxo)dicopper(II) structure with copper ...

A new method is presented to calculate binding energies and eigenfunctions for molecules, using the Dirac–Slater Hamiltonian. A numerical basis set of four component wavefunctions is obtained from atom‐like Dirac–Slater wavefunctions. A discrete variational method (DVM) has been applied to generate the binding energies and eigenfunctions for the molecule. Results are given for a series of molecules, including dihydrides H2X (X=O, S, Se, Te), diatomic indium halides InX (X=F, Cl, Br, I), and metal chlorides XCl (X=B, Al, Ga, In, Tl). Comparison is made with results from nonrelativistic calculations using the DVM with numerical Hartree–Fock–Slater‐type wavefunctions and with other types of nonrelativistic calculations. In particular, relativistic level shifts and spin–orbit splitting have been analyzed. The theoretical ionization energies are compared with experimental results. Generally a very good agreement is obtained between the experimental and theoretical binding energies for the valence levels, calcu...

A numerical-variational method for performing self-consistent molecular calculations in the Hartree-Fock-Slater (HFS) model is presented. Molecular wavefunctions are expanded in terms of basis sets constructed from numerical HFS solutions of selected one-center atomlike problems. Binding energies and wavefunctions for the molecules are generated using a discrete variational method for a given molecular potential. In the self-consistent-charge (SCC) approximation to the complete self-consistent-field (SCF) method, results of a Mulliken population analysis of the molecular eigenfunctions are used in each iteration to produce 'atomic' occupation numbers. The simplest SCC potential is then obtained from overlapping spherical atomlike charge distributions. Molecular ionization energies are calculated using the transition-state procedure; results are given for CO, H2O, H2S, AlCl, InCl, and the Ni5O surface complex. Agreement between experimental and theoretical ionization energies for the free-molecule valence levels is generally within 1 eV. The simple SCC procedure gives a reasonably good approximation to the molecular potential, as shown by comparison with experiment, and with complete SCF calculations for CO, H2O, and H2S.