Porous liquids (PLs), which are solvent-based systems that contain permanent porosity due to the incorporation of a solid porous host, are of significant interest for the capture of greenhouse gases, including CO2. Type 3 PLs formed by using metal-organic frameworks (MOFs) as the nanoporous host provide a high degree of chemical turnability for gas capture. However, pore aperture fluctuation, such as gate-opening in zeolitic imidazole framework (ZIF) MOFs, complicates the ability to keep the MOF pores available for gas adsorption. Therefore, an understanding of the solvent molecular size required to ensure exclusion from MOFs in ZIF-based Type 3 PLs is needed. Through a combined computational and experimental approach, the solvent-pore accessibility of exemplar MOF ZIF-8 was examined. Density functional theory (DFT) calculations identified that the lowest-energy solvent-ZIF interaction occurred at the pore aperture. Experimental density measurements of ZIF-8 dispersed in various-sized solvents showed that ZIF-8 adsorbed solvent molecules up to 2 Å larger than the crystallographic pore aperture. Density analysis of ZIF dispersions was further applied to a series of possible ZIF-based PLs, including ZIF-67, −69, −71(RHO), and −71(SOD), to examine the structure-property relationships governing solvent exclusion, which identified eight new ZIF-based Type 3 PL compositions. Solvent exclusion was driven by pore aperture expansion across all ZIFs, and the degree of expansion, as well as water exclusion, was influenced by ligand functionalization. Using these results, a design principle was formulated to guide the formation of future ZIF-based Type 3 PLs that ensures solvent-free pores and availability for gas adsorption.
Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO2 molecules were unable to diffuse into the ZIF-8 pore, with CO2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO2 adsorption should be based on the strength of gas binding at the solvated particle surface.
Porous liquids (PLs) based on the zeolitic imidazole framework ZIF-8 are attractive systems for carbon capture since the hydrophobic ZIF framework can be solvated in aqueous solvent systems without porous host degradation. However, solid ZIF-8 is known to degrade when exposed to CO2 in wet environments, and therefore the long-term stability of ZIF-8-based PLs is unknown. Through aging experiments, the long-term stability of a ZIF-8 PL formed using the water, ethylene glycol, and 2-methylimidazole solvent system was systematically examined, and the mechanisms of degradation were elucidated. The PL was found to be stable for several weeks, with no ZIF framework degradation observed after aging in N2 or air. However, for PLs aged in a CO2 atmosphere, formation of a secondary phase occurred within 1 day from the degradation of the ZIF-8 framework. From the computational and structural evaluation of the effects of CO2 on the PL solvent mixture, it was identified that the basic environment of the PL caused ethylene glycol to react with CO2 forming carbonate species. These carbonate species further react within the PL to degrade ZIF-8. The mechanisms governing this process involves a multistep pathway for PL degradation and lays out a long-term evaluation strategy of PLs for carbon capture. Additionally, it clearly demonstrates the need to examine the reactivity and aging properties of all components in these complex PL systems in order to fully assess their stabilities and lifetimes.
Rare-earth terephthalic acid (BDC)-based metal-organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H2O, NO2, and SO2) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC-gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with −OH, −NH2, and −SH were investigated. The calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.
Understanding the selectivity of metal–organic frameworks (MOFs) to complex acid gas streams will enable their use in industrial applications. Herein, ab initio molecular dynamic simulations (AIMD) were used to simulate ternary gas mixtures (H2O-NO2-SO2) in rare earth 2,5-dihydroxyterephthalic acid (RE-DOBDC) MOFs. Stronger H2O gas-metal binding arose from thermal vibrations in the MOF sterically hindering access of SO2 and NO2 molecules to the metal sites. Gas-gas and gas-linker interactions within the MOF framework resulted in the formation of multiple secondary gas species including HONO, HNO2, NOSO, and HNO3−. Four gas adsorption sites were identified along with a new de-protonation reaction mechanism not observable through experiment. This study not only provides valuable information on competitive gas binding energies in the MOF, it also provides important chemical insights into transient chemical reactions and mechanisms.
Direct air capture (DAC) of CO2 is one of the negative emission technologies under development to limit the impacts of climate change. The dilute concentration of CO2 in the atmosphere (~400 ppm) requires new materials for carbon capture with increased CO2 selectivity that is not met with current materials. Porous liquids (PLs) are an emerging material that consist of a combination of solvents and porous hosts creating a liquid with permanent porosity. PLs have demonstrated excellent CO2 selectivity, but the features that control how and why PLs selectively capture CO2 is unknown. To elucidate these mechanisms, density functional theory (DFT) simulations were used to investigate two different PLs. The first is a ZIF-8 porous host in a water/glycol/2-methylimidazole solvent. The second is the CC13 porous organic cage with multiple bulky solvents. DFT simulations identified that in both systems, CO2 preferentially bound in the pore window rather than in the internal pore space, identifying that the solvent-porous host interface controls the CO2 selectivity. Additionally, SNL synthesized ZIF-8 based PL compositions. Evaluation of the long-term stability of the PL identified no change in the ZIF-8 crystallinity after multiple agitation cycles, identifying its potential for use in carbon capture systems. Through this project, SNL has developed a fundamental understanding of solvent-host interactions, as well as how and where CO2 binds in PLs. Through these results, future efforts will focus not on how CO2 behaves inside the pore, but on the porous host-solvent interface as the driving force for PL stability and CO2 selectivity.
Rare-earth polynuclear metal-organic frameworks (RE-MOFs) have demonstrated high durability for caustic acid gas adsorption and separation based on gas adsorption to the metal clusters. The metal clusters in the RE-MOFs traditionally contain RE metals bound by μ3-OH groups connected via organic linkers. Recent studies have suggested that these hydroxyl groups could be replaced by fluorine atoms during synthesis that includes a fluorine-containing modulator. Here, a combined modeling and experimental study was undertaken to elucidate the role of metal cluster fluorination on the thermodynamic stability, structure, and gas adsorption properties of RE-MOFs. Through systematic density-functional theory calculations, fluorinated clusters were found to be thermodynamically more stable than hydroxylated clusters by up to 8-16 kJ/mol per atom for 100% fluorination. The extent of fluorination in the metal clusters was validated through a 19F NMR characterization of 2,5-dihydroxyterepthalic acid (Y-DOBDC) MOF synthesized with a fluorine-containing modulator. 19F magic-angle spinning NMR identified two primary peaks in the isotropic chemical shift (δiso) spectra located at -64.2 and -69.6 ppm, matching calculated 19F NMR δiso peaks at -63.0 and -70.0 ppm for fluorinated systems. Calculations also indicate that fluorination of the Y-DOBDC MOF had negligible effects on the acid gas (SO2, NO2, H2O) binding energies, which decreased by only ∼4 kJ/mol for the 100% fluorinated structure relative to the hydroxylated structure. Additionally, fluorination did not change the relative gas binding strengths (SO2 > H2O > NO2). Therefore, for the first time the presence of fluorine in the metal clusters was found to significantly stabilize RE-MOFs without changing their acid-gas adsorption properties.