Cutting-edge molecular modelling to unveil new microscopic insights into the guest-controlled flexibility of metal-organic frameworks

Metal-organic frameworks are a class of porous solids that exhibit intriguing flexibility under stimuli, leading often to reversible giant structural changes upon guest adsorption. DUT-49(Cu) and MIL-53(Cr) are fascinating flexible MOFs owing to their guest-induced breathing and negative gas adsorption behaviors respectively. Molecular simulation is one of the most relevant tools to examine these phenomena at the atomistic scale and gain a unique understanding of the physics behind them. Although molecular dynamics and Monte Carlo simulations are widely used in the field of porous materials, these methods hardly consider the structural deformation of a soft material upon guest adsorption. In this work, a cutting-edge osmotic molecular dynamics approach is developed to consider simultaneously the fluid adsorption process and material flexibility. We demonstrate that this newly developed computational strategy offers a unique opportunity to gain unprecedented molecular insights into the flexibility of this class of materials.

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH₄, C₂H₄, C₂H₆, C₃H₆, C₃H₈, and CO₂ across a temperature range of 250-400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure-temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption-desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300-550 K. A temperature threshold of ∼400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

Molecular Characterization of Membrane Gas Separation under Very High Temperatures and Pressure: Single- and Mixed-Gas CO2/CH4 and CO2/N2 Permselectivities in Hybrid Networks

This work illustrates the potential of using atomistic molecular dynamics (MD) and grand-canonical Monte Carlo (GCMC) simulations prior to experiments in order to pre-screen candidate membrane structures for gas separation, under harsh conditions of temperature and pressure. It compares at 300 °C and 400 °C the CO₂/CH₄ and CO₂/N₂ sieving properties of a series of hybrid networks based on inorganic silsesquioxanes hyper-cross-linked with small organic PMDA or 6FDA imides. The inorganic precursors are the octa(aminopropyl)silsesquioxane (POSS), which degrades above 300 °C, and the octa(aminophenyl)silsesquioxane (OAPS), which has three possible meta, para or ortho isomers and is expected to resist well above 400 °C. As such, the polyPOSS-imide networks were tested at 300 °C only, while the polyOAPS-imide networks were tested at both 300 °C and 400 °C. The feed gas pressure was set to 60 bar in all the simulations. The morphologies and densities of the pure model networks at 300 °C and 400 °C are strongly dependent on their precursors, with the amount of significant free volume ranging from ~2% to ~20%. Since measurements at high temperatures and pressures are difficult to carry out in a laboratory, six isomer-specific polyOAPS-imides and two polyPOSS-imides were simulated in order to assess their N₂, CH₄ and CO₂ permselectivities under such harsh conditions. The models were first analyzed under single-gas conditions, but to be closer to the real processes, the networks that maintained CO₂/CH₄ and CO₂/N₂ ideal permselectivities above 2 were also tested with binary-gas 90%/10% CH₄/CO₂ and N₂/CO₂ feeds. At very high temperatures, the single-gas solubility coefficients vary in the same order as their critical temperatures, but the differences between the penetrants are attenuated and the plasticizing effect of CO₂ is strongly reduced. The single-gas diffusion coefficients correlate well with the amount of available free volume in the matrices. Some OAPS-based networks exhibit a nanoporous behavior, while the others are less permeable and show higher ideal permselectivities. Four of the networks were further tested under mixed-gas conditions. The solubility coefficient improved for CO₂, while the diffusion selectivity remained similar for the CO₂/CH₄ pair and disappeared for the CO₂/N₂ pair. The real separation factor is, thus, mostly governed by the solubility. Two polyOAPS-imide networks, i.e., the polyorthoOAPS-PMDA and the polymetaOAPS-6FDA, seem to be able to maintain their CO₂/CH₄ and CO₂/N₂ sieving abilities above 2 at 400 °C. These are outstanding performances for polymer-based membranes, and consequently, it is important to be able to produce isomer-specific polyOAPS-imides for use as gas separation membranes under harsh conditions.