Competitive adsorption of water is an important issue in the adsorption-based industrial processes of bio- and flue gases separation. The dehumidification of gases prior to separation would increase process complexity and lower its economic interest. In this work, large-scale computational screening was applied to identify Metal-Organic Frameworks (MOFs) structures which exhibit high CO2/CH4 selectivity and total loading higher than 0.5 mol/kg (in the presence of water). High-throughput Grand Canonical Monte Carlo (GCMC) screening of nearly 3000 existing MOF materials was carried out. Initial selection assumed fixed values of pore limiting diameter (PLD) and Henry's constant for water and allowed one to preselect 764 structures. After GCMC simulations carried for 50/50 CO2/CH4 mixture, at ambient conditions (p = 1 bar, T = 298 K), and variable gas humidity (0%, 5%, 30% and 40%) the final selection revealed 13 most promising MOFs structures. We focused on analysis of the correlations between the properties of the selected MOFs and the separation selectivity. We show that the selectivity is a complex function of the porous materials characteristics and finding selective sorbent, performing well in dry and wet conditions requires careful analysis of available MOFs.

Soft porous crystals undergo large structural transformations under a variety of physical stimuli. Breathing-like transformations, occurring with a large volume change, have been associated with an existence of bi-stable or multi-stable crystal structures. Here, we provide a detailed, multiscale qualitative and quantitative analysis of the adsorption-induced “breathing” transformations in two metal organic frameworks (MOFs): MIL-53(Al) which is a reference case of our approach, and recently synthetized JUK-8, which does not show any bistability without adsorbate. The proposed approach is based on atomistic simulations and does not require any empirical or adjustable parameters. It allows for a prediction of potential structural transformations in MOF and, in the case of adsorption-induced deformations, a calculation of the adsorption stress. We also show that the quantitative agreement between calculated and experimental results critically depends on the quality of the dispersion energy correction.

We discuss the mechanism of the structural transformations of gas adsorbed in microporous crystalline solids. We show that it is possible to induce structural transformations in a confined system by simply varying the number of molecules adsorbed in the pore. We found that the mechanism of these novel, adsorption-induced structural transformation in nano-pores differs from the capillary condensation. First, the structure of the confined gas is determined by a competition between adsorption sites attractive forces and intermolecular interaction. Second, at low temperature, the transformation is discontinuous because it is defined by limited number of accessible adsorption sites [1,2]. In the case of methane adsorbed in IRMOFs porous structures the character of transformation depends on temperature but also on the IRMOF structure: it changes from strongly discontinuous (especially at low temperatures), to continuous transition. The mechanism of the transformation is also modified by the size of the gas molecules and the strength of interaction. However, even discontinuous transition can produce continuous isotherms. We will show that the continuous isotherm can be the effect of statistical dynamical switching between two phase, characterized by different number of adsorbed molecules. We show the simulated microscopic mechanism and experimental observations support such statistical interpretation. This work was supported by the Polish National Science Centre (NCN, grant no. 2015/17/B/ST8/00099). The calculations have been partially performed at the WCSS computer center of The Wroclaw University of Science and Technology, grant no 33.References[1] Heterogeneous melting of methane confined in nano-pores. Dundar, E.; Boulet, P.; Wexler, C.; et al. J. Chem. Phys. 145, 144704 (2016)[2] Adsorption-Induced Structural Phase Transformation in Nanopores. Kuchta, Bogdan; Dundar, Ege; Formalik, Filip; et al. Angewandte Chimie In. Ed. 56, 16243–16246. (2017)

Introduction/PurposeStorage and separation of gases are one of the most challenging processes when they have to be implemented at industrial scale. Hundreds of potential materials have been tested over the past years, especially for gas storage (hydrogen and methane for fueling applications) and separation of biogas (a mixture of 55 – 75% of methane 25 – 45% of carbon dioxide). Today, MOFs (Metal-Organic Frameworks) [1-2] are considered as the most promising materials for such applications. MOFs are hybrid, nanoporous materials with a very large surface area, in which pore sizes, shapes and adsorption energies can be almost unlimitedly tailored to fit the required adsorption characteristics.MethodsIn this work, we present simulations of methane and carbon dioxide adsorption in a wide range of MOF structures from CoRE-MOF database [3], and at variable content of humidity. Over 1000 structures have been tested using Grand Canonical Monte Carlo method (as implemented in RASPA code) to find the optimal material for the CH4/CO2 separation in conditions of increasing humidity.ResultsMOFs with pore limiting diameter (PLD) over 3.3 Å were chosen for the simulations. Due to the humid conditions, we chose structures from the middle range of hydrophilicity, kH between 5·10-2 mol/kg·Pa (CuBTC) and 5·10-6 mol/kg·Pa (ZIF-8). The initial screening allowed us to select the most promising structures for more detailed investigation. We show that 20 materials exhibit very promising selectivity (defined as the adsorbed amount of CO2­ to the amount of CH4) in the range from 100 to 900 and total uptake above 0.5 mol/kg.ConclusionsThe results provide clear information about different mechanisms of adsorption for selected structures. Furthermore, the adsorption of methane tends to be reduced with increasing humidity. Materials with a growing level of selectivity for CO2 and CH4 with increasing hydration degree are the most promising from the future application point of view.

Low frequency lattice vibrations (phonons) can be used as indicators of a variety of structural transformations. In particular, conformational changes in flexible metal-organic frameworks (MOFs) can be successfully explained by analyzing materials phonon spectrum in frequency range below 200 cm-1 [1,2]. It has been shown that such lattice vibrations are responsible for pressure- and adsorption-induced gate-opening in zeolitic imidazolate frameworks (ZIFs), and pore breathing in materials from MIL-53 family. In this work we present a phonon-based methodology to analyze MOFs deformations that can be detected in adsorption experiments but cannot be explained at the microscopic level without theoretical support based on numerical simulations.As a case study we have chosen to analyze adsorption-induced transformations in ZIF-8. The step-like adsorption isotherm observed in this material was interpreted in the literature as a signature of gate-opening (enlargement of the pore free volume) induced by the increasing number of guest molecules entering the pore. To analyze this phenomenon at the microscopic level, we first calculated and visualized normal modes of the lattice vibrations. For that, density functional theory (DFT) with PBE functional and empirical dispersion correction D3(BJ) was used. Only the modes which fulfill two criteria: (i) their frequency is below 200 cm-1, and (ii) the related lattice deformation leads to increase of the pore free volume, have been further analyzed. These deformed structures, (including one related to gate-opening, pore breathing and one which combines both deformations) were selected to serve as the adsorbent model. Grand Canonical Monte Carlo (GCMC) simulations have been performed to calculate the isotherms of adsorption in selected structures. Surprisingly, in the structure related to gate-opening, the mode, located at ~40 cm-1, and considered in the literature as the one responsible for step-like form of the isotherm did not provided an uptake larger than the undistorted structure. On the other hand, in the structure resulting from simultaneous pore breathing and gate-opening deformation the adsorption is significantly increased. Therefore, we consider that this deformation, induced by pore breathing and gate-opening phonons explains better the experimental observations..1. F. Formalik, et al., Microporous and Mesoporous Materials, (2018),2. M. R. Ryder, et al., Physical Review Letters, 113 (2014), 215502-21550

Storage and separation of gases is one of the most challenging processes when they have to be implemented at industrial scale. Hundreds of potential materials have been tested over the past years, especially for gas storage (hydrogen and methane for fueling applications) and separation of biogas (a mixture of 55 – 75% of methane 25 – 45% of carbon dioxide). Today, MOFs (Metal-Organic Frameworks) [1-2] are considered as the most promising materials for such applications. MOFs are hybrid, nanoporous materials with a very large surface area, in which pore sizes, shapes and adsorption energies can be almost unlimitedly tailored to fit the required adsorption characteristics. A particular category of these materials are the so-called breathing MOFs, in which a reversible structural transformation (variation of pore volume, pore breathing) can be induced by temperature, pressure or interaction with adsorbate molecules. Molecular modeling is an effective way to rapidly obtain microscopic information on the mechanism of adsorption in these materials (in many cases, not accessible using traditional experimental methods) and to shorten the time necessary to implement research results in an industrial process.In the this work we will present simulations of methane and carbon dioxide adsorption in a wide range of MOF structures from CoRE-MOF database [3], and at variable content of humidity. MOFs with pore limiting diffusion (PLD) over 3.3 Å were chosen for the simulations. Due to the humid conditions we chose structures from the middle range of hydrophilicity, kH between 5·10-2 mol/kg·Pa (CuBTC) and 5·10-6 mol/kg·Pa (ZIF-8). Over 1000 structures have been tested using Grand Canonical Monte Carlo method (as implemented in RASPA code) to find the optimal material for the CH4/CO2 separation in conditions of increasing humidity. The initial screening allowed us to select the most promising structures for more detailed investigation. We show that 16 materials exhibit very promising selectivity (defined as adsorbed amount of CO2¬ to amount of CH4) in range from 50 to 440. Furthermore, adsorption of methane tends to be reduced with increasing humidity. Reference[1]S. Li, Y. G. Chung, and R. Q. Snurr, “High-Throughput Screening of Metal–Organic Frameworks for CO 2 Capture in the Presence of Water,” Langmuir, vol. 32, no. 40, pp. 10368–10376, 2016.[2]R. Krishna, “Screening metal–organic frameworks for mixture separations in fixed-bed adsorbers using a combined selectivity/capacity metric,” RSC Adv., vol. 7, no. 57, pp. 35724–35737, 2017.[3]Y. G. Chung et al., “Computation-Ready, Experimental Metal–Organic Frameworks: A Tool To Enable High-Throughput Screening of Nanoporous Crystals,” Chem. Mater., vol. 26, no. 21, pp. 6185–6192, Nov. 2014.

Phonons are quantum elastic excitations ofcrystalline solids. Classically, they correspond to thecollective vibrations of atoms in ordered periodic structures.They determine the thermodynamic properties ofsolids and their stability in the case of structural transformations.Here we review for the first time the existingexamples of the phonon analysis of adsorption-inducedtransformations occurring in microporous crystallinematerials. We discuss the role of phonons in determiningthe mechanism of the deformations. We point outthat phonon-based methodology may be used as a predictivetool in characterization of flexible microporousstructures; therefore, relevant numerical tools must bedeveloped.