Ab initio molecular dynamics determination of competitive O₂ vs. N₂ adsorption at open metal sites of M₂(dobdc)

Ab Initio Molecular Dynamics Investigation of Water and Butanone Adsorption on UiO-66 with Defects

Volatile organic compounds (VOCs) are harmful chemicals that are found in minute quantities in the atmosphere and are emitted from a variety of industrial and biological processes. They can be harmful to breathe or serve as biomarkers for disease detection. Therefore, capture and detection of VOCs is important. Here, we have examined if the Zr-based UiO-66 metal-organic framework (MOF) can be used to capture butanone─a well-known VOC. Toward that end, we have performed Born-Oppenheimer ab initio molecular dynamics (AIMD) at 300 and 500 K to probe the energetics and molecular interactions between butanone [CH3C(O)CH2CH3] and open-cage Zr-UiO-66. Such interactions were systematically interrogated using three MOF structures: defective MOF with a missing 1,4-benzene-dicarboxylate linker and two H2O; pristine MOF with two H2O; and pristine dry MOF. These structures were loaded with one and four molecules of butanone to examine the effect of concentration as well. One-molecule loading interacted favorably with the defective structure at 300 K, only. In comparison, interactions with four-molecule loading were energetically favorable for all conditions. Persistent hydrogen bonds between the O atom of butanone, H2O, and μ3-OH hydroxyl attachments at Zr nodes substantially contributed to the intermolecular interactions. At higher loadings, butanone also showed a pronounced tendency to diffuse into the adjoining cages of Zr-UiO-66. The effect of such movement on interaction energies was rationalized using simple statistical mechanics-based models of interacting and noninteracting gases. Broadly, we learn that the presence of prior moisture within the interstitial cages of Zr-UiO-66 significantly impacts the adsorption behavior of butanone.

Porous Liquids: Computational Design for Targeted Gas Adsorption

In this Perspective, we present the unique gas adsorption capabilities of porous liquids (PLs) and the value of complex computational methods in the design of PL compositions. Traditionally, liquids only contain transient pore space between molecules that limit long-term gas capture. However, PLs are stable fluids that that contain permanent porosity due to the combination of a rigid porous host structure and a solvent. PLs exhibit remarkable adsorption and separation properties, including increased solubility and selectivity. The unique gas adsorption properties of PLs are based on their structure, which exhibits multiple gas binding sites in the pore and on the cage surface, varying binding mechanisms including hydrogen-bonding and π-π interactions, and selective diffusion in the solvent. Tunable PL compositions will require fundamental investigations of competitive gas binding mechanisms, thermal effects on binding site stability, and the role of nanoconfinement on gas and solvent diffusion that can be accelerated through molecular modeling. With these new insights PLs promise to be an exceptional material class with tunable properties for targeted gas adsorption.

Design Elements for Enhanced Hydrogen Isotope Separations in Barely Porous Organic Cages

Barely porous organic cages (POCs) successfully separate hydrogen isotopes (H₂/D₂) at temperatures below 100 K. Identifying the mechanisms that control the separation process is key to the design of next-generation hydrogen separation materials. Here, ab initio molecular dynamics (AIMD) simulations are used to elucidate the mechanisms that control D₂ and H₂ separation in barely POCs with varying functionalization. The temperature and pore size dependence were identified, including the selective capture of D₂ in three different CC3 structures (RCC3, CC3-S, and 6ET-RCC3). The temperature versus capture trend was reversed for the 6ET-RCC3 structure, identifying that the D₂ and H₂ escape mechanisms are unique in highly functionalized systems. Analysis of calculated isotope velocities identified effective pore sizes that extend beyond the pore opening distances, resulting in increased capture in minimally functionalized CC3-S and RCC3. In a highly functionalized POC, 6ET-RCC3, higher velocities of the H isotopes were calculated moving through the restricted pore compared to the rest of the system, identifying a unique molecular behavior in the barely nanoporous pore openings. By using AIMD, mechanisms of H₂ and D₂ separation were identified, allowing for the targeted design of future novel materials for hydrogen isotope separation.

Kinetically Controlled Linker Binding in Rare Earth-2,5-Dihydroxyterepthalic Acid Metal-Organic Frameworks and Its Predicted Effects on Acid Gas Adsorption

In the pursuit of highly stable and selective metal-organic frameworks (MOFs) for the adsorption of caustic acid gas species, an entire series of rare earth MOFs have been explored. Each of the MOFs in this series (RE-DOBDC; RE = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu; DOBDC = 2,5-dihydroxyterepthalic acid) was synthesized in the tetragonal space group I4/m. Crystallized MOF samples, specifically Eu-DOBDC, were seen to have a combination of monodentate and bidentate binding when synthesized under typical reaction conditions, resulting in a contortion of the structure. However, extended crystallization times determined that this binding is kinetically controlled and that the monodentate binding option was crystallographically eliminated by extended reaction times at higher temperatures. Furthermore, this series allows for the direct study of the effect of the metal center on the structure of the of the MOF; herein, the lanthanide metal ionic radii contraction across the periodic table results in a reduction of the MOF pore size and lattice parameters. Scanning electron microscopy-energy-dispersive spectroscopy was used to investigate the stages of crystal growth for these RE-DOBDC MOFs. All MOFs, except Er-DOBDC had a minimum of two stages of growth. These analogues were demonstrated by analysis of neutron diffraction (PND) to exhibit a cooperative rotational distortion of the secondary building unit, resulting in two crystallographically distinct linker sublattices. Computational modeling efforts were used to show distinct differences on acid gas (NO₂ and SO₂) binding energies for RE-DOBDC MOFs when comparing the monodentate/bidentate combined linker with the bidentate-only linker crystal structures.

Prediction of Reactive Nitrous Acid Formation in Rare-Earth MOFs via ab initio Molecular Dynamics

Reactive gas formation in pores of metal-organic frameworks (MOFs) is a known mechanism of framework destruction; understanding those mechanisms for future durability design is key to next generation adsorbents. Herein, an extensive set of ab initio molecular dynamics (AIMD) simulations are used for the first time to predict competitive adsorption of mixed acid gases (NO2 and H2 O) and the in-pore reaction mechanisms for a series of rare earth (RE)-DOBDC MOFs. Spontaneous formation of nitrous acid (HONO) is identified as a result of deprotonation of the MOF organic linker, DOBDC. The unique DOBDC coordination to the metal clusters allows for proton transfer from the linker to the NO2 without the presence of H2 O and may be a factor in DOBDC MOF durability. This is a previously unreported mechanisms of HONO formation in MOFs. With the presented methodology, prediction of future gas interactions in new nanoporous materials can be achieved.

Tuned Hydrogen Bonding in Rare-Earth Metal-Organic Frameworks for Design of Optical and Electronic Properties: An Exemplar Study of Y-2,5-Dihydroxyterephthalic Acid

Organic linkers in metal-organic framework (MOF) materials exhibit differences in hydrogen bonding (H-bonding), which can alter the geometric, electronic, and optical properties of the MOF. Density functional theory (DFT) simulations were performed on a photoluminescent Y-2,5-dihydroxyterephthalic acid (DOBDC) MOF with H-bonding concentrations between 0 and 100%; the H-bonds were located on both bidentate- and monodentate-bound DOBDC linkers. At 0% H-bond concentration in the framework, the lattice parameters contracted, the density increased, and simulated X-ray diffraction patterns shifted. Comparison with published experimental data identified that Y-DOBDC MOF structures must have a degree of H-bond concentration. The concentration of H-bonds in the system shifted the calculated band gap energy from 2.25 eV at 100% to 3.00 eV at 0%. The band gap energies also indicate a distinction of H-bonds formed on bidentate-coordinated linkers compared to those on monodentate linkers. Additionally, when the calculated optical spectra are compared with experimental data, the ligand-to-ligand charge-transfer luminescence in Y-DOBDC MOFs is expected to result from an average of 20-40% H-bonding with at least 50% of the bidentate linkers containing H-bonding. Therefore, the type of H-bonding within the DOBDC linker determines the electronic structure and the optical absorption of the MOF framework structure. Tuning of the H-bonding in rare-earth MOFs provides an opportunity to control the specific optical and adsorption properties of the MOF framework on the basis of reactions between the linker and the environment.

NOx Adsorption and Optical Detection in Rare Earth Metal-Organic Frameworks

Acid gases (e.g., NO_x and SO_x), commonly found in complex chemical and petrochemical streams, require material development for their selective adsorption and removal. Here, we report the NO_x adsorption properties in a family of rare earth (RE) metal-organic frameworks (MOFs) materials. Fundamental understanding of the structure-property relationship of NO_x adsorption in the RE-DOBDC materials platform was sought via a combined experimental and molecular modeling study. No structural change was noted following humid NO_x exposure. Density functional theory (DFT) simulations indicated that H₂O has a stronger affinity to bind with the metal center than NO₂, while NO₂ preferentially binds with the DOBDC ligands. Further modeling results indicate no change in binding energy across the RE elements investigated. Also, stabilization of the NO₂ and H₂O molecules following adsorption was noted, predicted to be due to hydrogen bonding between the framework ligands and the molecules and nanoconfinement within the MOF structure. This interaction also caused distinct changes in emission spectra, identified experimentally. Calculations indicated that this is due to the adsorption of NO₂ molecules onto the DOBDC ligand altering the electronic transitions and the resulting photoluminescent properties, a feature that has potential applications in future sensing technologies.

Structure and electronic properties of rare earth DOBDC metal-organic-frameworks

Here, we apply density functional theory (DFT) to investigate rare-earth metal organic frameworks (RE-MOFs), RE12(μ3-OH)16(C8O6H4)8(C8O6H5)4 (RE = Y, Eu, Tb, Yb), and characterize the level of theory needed to accurately predict structural and electronic properties in MOF materials with 4f-electrons. A two-step calculation approach of geometry optimization with spin-restricted DFT and large core potential (LCPs), and detailed electronic structures with spin-unrestricted DFT with a full valence potential + Hubbard U correction is investigated. Spin-restricted DFT with LCPs resulted in good agreement between experimental lattice parameters and optimized geometries, while a full valence potential is necessary for accurate representation of the electronic structure. The electronic structure of Eu-DOBDC MOF indicated a strong dependence on the treatment of highly localized 4f-electrons and spin polarization, as well as variation within a range of Hubbard corrections (U = 1-9 eV). For Hubbard corrected spin-unrestricted calculations, a U value of 1-4 eV maintains the non-metallic character of the band gap with slight deviations in f-orbital energetics. When compared with experimentally reported results, the importance of the full valence calculation and the Hubbard correction in correctly predicting the electronic structure is highlighted.

Biomimetic O2-Evolving metal-organic framework nanoplatform for highly efficient photodynamic therapy against hypoxic tumor

Improving the supply of O₂ and the circulation lifetime of photosensitizers for photodynamic therapy (PDT) in vivo would be a promising approach to eliminate hypoxic tumors. Herein, by taking advantage of the significant gas-adsorption capability of metal-organic frameworks (MOFs), a biomimetic O₂-evolving photodynamic therapy (PDT) nanoplatform with long circulating properties was fabricated. Zirconium (IV)-based MOF (UiO-66) was used as a vehicle for O₂ storing, then conjugated with indocyanine green (ICG) by coordination reaction, and further coated with red blood cell (RBC) membranes. Upon 808 nm laser irradiation, the initial singlet oxygen (¹O₂) generated by ICG would decompose RBC membranes. At the same time, The photothermal property of ICG could facilitate the burst release of O₂ from UiO-66. Subsequently, the generated O₂ could significantly improve the PDT effects on hypoxic tumor. Owing to the advantages of long circulation and O₂ self-sufficient, the designed nanotherapeutic agent can improve the efficiency of treatment against hypoxia tumor via PDT. Hence, this study presents a new paradigm for co-delivery of O₂ and photosensitizers, and provides a new avenue to eliminate hypoxic tumors.

Mathematical modeling of gas desorption from a metal–organic supercontainer cavity filled with stored N2gas at critical limits

Gas escape rates from within the cavity of a MOSC were predicted by molecular dynamics and analytical equations.

On the elusive nature of oxygen binding at coordinatively unsaturated 3d transition metal centers in metal–organic frameworks

Low and ambient temperature binding of oxygen, O₂, in MOF-74, CPO-27-M (M = Mn, Co, Ni, Cu, Zn) framework materials remains in the physisorption regime, with energetics very similar to that of nitrogen, N₂, sorption.