High pressure sorption of various hydrocarbons and carbon dioxide in Kimmeridge Blackstone and isolated kerogen

Revisiting methane absolute adsorption in organic nanopores from molecular simulation and Ono-Kondo lattice model

Tackling the challenges in the estimation of methane absolute adsorption in kerogen nanoporous media from molecular and analytical approaches

As the world moves towards carbon neutrality, polymer electrolyte fuel cells (PEFCs) have emerged as a vital technology in recent years. The cathode catalyst layer of PEFCs is characterized by its porous structure, which promotes the production of water molecules from oxygen and protons through chemical reactions. The ionomer thin films, predominantly composed of Nafion, are critical in regulating the transport and mobility of molecules and ions in the cathode catalyst layer. Sulfonic groups in the side chains of Nafion enable accumulation of water molecules, causing the ionomer thin films to swell.Previous studies have examined the swelling of Nafion under wet conditions[1][2]. However, precise swelling mechanisms remain unclear. While privious research has analyzed the relationship between relative humidity and the amount of water within bulk Nafion samples[3], the swelling behavior in films with thickness of several nanometers differs significantly, necessitating further investigation[4]. This study therefore investigates the swelling mechanisms and the properties of ionomer thin films used in the cathode catalyst layer of PEFCs. This study aims to establish guidelines for designing the molecular structures of ionomers by elucidating the mechanisms underlying ionomer swelling in the cathode catalyst layer. Grand canonical Monte Carlo (GCMC) simulations and molecular dynamics (MD) simulations were employed to reproduce swelling in Nafion thin films and examine the impact of Nafion's structure on water content and swelling. These molecular simulation techniques enable the elucidation of the influence exerted by polymer structure on water content and swelling at the molecular scale.In this study, we firstly created a simulation model of a Nafion system using LAMMPS[5], the DREIDING Force Field[6], and Mulliken charges[7] in which parameter settings were referred from literature values[8]. The Nafion system contains 10 Nafion molecules with a degree of polymerization of 10, water molecules, and 100 oxonium ions for electrical neutrality. The density of 1.9 g/cm3 achieved during simulations closely agrees with literature values[8].To reproduce Nafion swelling under humid conditions, GCMC simulations and MD simulations were employed. The rationale behind utilizing GCMC simulations lies in the fact that they enable the insertion of water molecules into the Nafion system under constant chemical potential conditions, thereby inducing Nafion swelling. Subsequent to the GCMC simulations, MD simulations were performed under the NPT ensemble, which facilitated the system to reach equilibrium. This alternating cycle of GCMC and MD simulations is repeated to realize the final swelling state. Preliminary results indicate lower water content values than those experimentally measured under similar conditions. Possible causes include inadequate force field parameters that lead to polymer morphology different from the actual Nafion structure.Future work will explore the causes of these discrepancies and examine the effects of Nafion's structure on water content and swelling under various conditions. This will involve testing alternative force field parameters and initial structures to better represent the molecular-level structure of ionomer thin films on carbon support particles in the cathode catalyst. Additionally, we will investigate how different relative humidity conditions, temperatures, and ionomer compositions influence swelling behavior, which could provide insights into optimizing PEFCs performance under various operating environments.[1] G. C. Abuin et al., J. Membr. Sci. 428, 507 (2013).[2] J. Catalano et al., Int. J. Hydrog. Energy 37, 6308 (2012).[3] T. A. Zawadzinski et al., J. Electrochem. Soc. 140, 1981 (1993).[4] M. A. Modestino et al., Macromolecules 46, 867 (2013).[5] A. P. Thompson et al., Comp. Phys. Comm. 271, 10817 (2022).[6] S. L. Mayo et al., J. Phys. Chem. 94, 8897 (1990).[7] S. S. Jang et al., J. Phys. Chem. B 108, 3149 (2004).[8] Y. Kurihara et al., J. Power Sources 414, 263 (2019).Fig. 1 An ionomer model composed of Nafion molecules, oxonium ions, and water molecules. The red and white spheres represent oxygen and hydrogen atoms making up water molecules. The size of the water molecules is emphasized for clarity. Figure 1

ReaxFF Grand Canonical Monte Carlo simulation of adsorption and dissociation of oxygen on platinum (111)

Tailoring Zirconium-based metal organic frameworks for enhancing Hydrophilic/Hydrophobic Characteristics: Simulation and experimental investigation

Determination of the absolute adsorption/desorption isotherms of CH4 and n-C4H10 on shale from a nano-scale perspective

Adsorption of carbon dioxide in slit-shaped carbon micropores at 273 K has been studied by means of the grand canonical Monte Carlo (GCMC) simulations and the nonlocal density functional theory (NLDFT). Three molecular models of CO2 have been used. Long-run GCMC simulations were performed with the three-center model of Harris and Yung (J. Phys. Chem. 1995, 99, 12021). For NLDFT calculations, we developed an effective Lennard-Jones (LJ) model. GCMC simulations of the effective LJ model of CO2 have been performed for comparison. For each model used, parameters of intermolecular potentials have been determined and validated against two-phase bulk equilibrium data and experimental adsorption isotherms on graphite at 273 and 195 K. In the range of pore widths from 3 to 15 Å, the NLDFT isotherms of CO2 adsorption are overall in a satisfactory agreement with the GCMC isotherms generated using the three-center model. Some deviations have been observed between 6.5 and 8.5 Å, where the adsorbate undergoes a transition from a single-layer to a two-layer structure. The models developed are recommended for studying carbon dioxide adsorption in microporous adsorbents and also for calculating pore size distributions in carbonaceous materials and soil particles. The NLDFT model has the advantage of being much less computationally demanding, whereas the three-center GCMC model serves as a benchmark for quantitative estimates and can be used for studying CO2 sorption at ambient conditions close to the critical temperature.

Evaluation of CH4 and CO2 adsorption on HKUST-1 and MIL-101(Cr) MOFs employing Monte Carlo simulation and comparison with experimental data

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindrical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with different numbers of adsorption and desorption steps are not common. To better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and found that these yielded comparable results. The transition probability provides a measure of the length of GCMC simulations to yield reliable results.

Grand canonical Monte Carlo (GCMC) and liquid-vapor molecular dynamics (LVMD) simulations are performed to investigate the squeezing and phase transition of a simple liquid argon film confined between two solid surfaces. Simulation results show that the LVMD simulation is capable of capturing the major thermodynamic equilibrium states of the confined film, as predicted by the GCMC simulations. Moreover, the LVMD simulations reveal the non-equilibrium squeeze out dynamics of the confined film. The study shows that the solvation force hysteresis, observed in many surface force experiments, is attributed to two major effects. The first is related to the unstable jumps during the laying transitions of the confined film, in which the gradient of force profile is larger than the driving spring constant. The second effect is related to the squeeze out dynamics of the confined film even though the first effect is absent. In general, these two dynamic processes are non-equilibrium in nature and involve significant energy dissipations, resulting in the force hysteresis.

Open AccessCCS ChemistryRESEARCH ARTICLE3 Oct 2022Fluorous Metal–Organic Frameworks with Unique Cage-in-Cage Structures Featuring Fluorophilic Pore Surfaces for Efficient C2H2/CO2 Separation Xing-Ping Fu†, Yu-Ling Wang†, Xue-Feng Zhang, Zhenjie Zhang, Chun-Ting He and Qing-Yan Liu Xing-Ping Fu† College of Chemistry and Chemical Engineering, Key Laboratory of Functional Small Molecules for Ministry of Education, Jiangxi Normal University, Nanchang, Jiangxi 330022 Department of Ecological and Resources Engineering, Fujian Key Laboratory of Eco-Industrial Green Technology, Wuyi University, Wuyishan, Fujian 354300 , Yu-Ling Wang† College of Chemistry and Chemical Engineering, Key Laboratory of Functional Small Molecules for Ministry of Education, Jiangxi Normal University, Nanchang, Jiangxi 330022 , Xue-Feng Zhang College of Chemistry and Chemical Engineering, Key Laboratory of Functional Small Molecules for Ministry of Education, Jiangxi Normal University, Nanchang, Jiangxi 330022 , Zhenjie Zhang College of Chemistry and State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071 , Chun-Ting He College of Chemistry and Chemical Engineering, Key Laboratory of Functional Small Molecules for Ministry of Education, Jiangxi Normal University, Nanchang, Jiangxi 330022 and Qing-Yan Liu *Corresponding author: E-mail Address: [email protected] College of Chemistry and Chemical Engineering, Key Laboratory of Functional Small Molecules for Ministry of Education, Jiangxi Normal University, Nanchang, Jiangxi 330022 https://doi.org/10.31635/ccschem.021.202101575 SectionsSupplemental MaterialAboutAbstractPDF ToolsAdd to favoritesDownload CitationsTrack Citations ShareFacebookTwitterLinked InEmail The similarities in molecular size and physical properties of acetylene (C2H2) and carbon dioxide (CO2) produce a formidable challenge for their separation. Herein, we report two isoreticular fluorinated metal–organic frameworks (MOFs), labeled as JXNU-11(Fe2M) (M = Ni and Co), featuring unique octahedral cages encapsulated by cuboctahedral cages. JXNU-11(Fe2Ni) shows a record high C2H2-capture amount of 4.8 mmol g−1 and a long C2H2/CO2 breakthrough interval time of 55 min g−1 in an actual breakthrough experiment based on the equimolar C2H2/CO2 mixture under ambient conditions, indicating a high-performance material for C2H2 capture and C2H2/CO2 separation. Computational simulations revealed that the nanosized double-shell cages decorated with abundant fluorine and carboxylate oxygen atoms afford the optimized pore spaces for preferentially trapping C2H2, which account for the remarkable C2H2-capture capacity and the highly efficient C2H2/CO2 separation for JXNU-11(Fe2M). This work not only develops a design strategy for building the cage-in-cage structures in MOFs, but also provides universal guidance on designing porous MOFs with highly polarized and fluorophilic pore cages to capture C2H2. Download figure Download PowerPoint Introduction Acetylene (C2H2), the simplest alkyne, is an essential chemical for many common chemical products, such as acetaldehyde, acetic acid, plastic, and synthetic rubber,1 and electronic materials in the electronics industry.2 Additionally, C2H2 is an important flammable gas widely utilized for welding and metal cutting. The production of C2H2 from the oxidative coupling of natural gas is a commonly used preparation method for C2H2 gas in the petrochemical industry.3 Such a C2H2 preparation process leads to a small amount of carbon dioxide (CO2) in the C2H2 product. The presence of CO2 in C2H2 starting material has a seriously negative effect on subsequent preparation processes for the down-stream chemical products. Thus separation of C2H2/CO2 mixtures is vital in the petrochemical industry. However, the similar molecular sizes, as well as the same kinetic diameter (3.3 Å) of the linear-shaped C2H2 and CO2 molecules4–9 ( Supporting Information Table S1), cause the C2H2/CO2 separation to be a formidable task. Organic solvent extraction and cryogenic distillation are the current technologies for purifying C2H2 gas but are energy-intensive and environmentally unfriendly. Therefore, it is important to develop a highly energy-efficient and environment-friendly technology for achieving C2H2/CO2 separation. Metal–organic frameworks (MOFs) with porous and crystalline structures are fascinating porous solid materials.5–14 MOFs show highly promising potential in gaseous mixture separation due to their designable frameworks' structures and tunable pore environments.15–19 Recently, many MOFs have been utilized for the separation of C2H2/CO2 mixtures,20–29 and a few of them exhibit high-performance C2H2/CO2 separation.30–34 In the column breakthrough separation of binary gases, the breakthrough time and the captured amount for the late eluted adsorbate in the column are the important parameters for evaluating the separation efficiency of an adsorbent. The long breakthrough time for the late eluted gas generally results in a large captured amount for the late eluted gas, leading to a high separation efficiency of the adsorbent. However, in practice breakthrough experiments based on the equimolar C2H2/CO2 mixture, C2H2 breakthrough time for MOFs is generally shorter than 50 min g−1 under ambient conditions.35–40 The top-performing MOFs are SNNU-45 and ZJU-74a with the C2H2 breakthrough times of 113 and 81 min g−1,41,42 respectively. In contrast, the C2H2-captured amounts for MOFs in the packed columns are commonly less than 3 mmol g−1. SNNU-45 and ZJU-74a have the highest C2H2-capture amounts of 3.5 and 3.64 mmol g−1, respectively. It has been well demonstrated that the fluorine atoms of the inorganic anions, including MF62− (M = Si4+ and Ti4+) and NbOF52,− in the hybrid MOFs are the preferential bind sites for the acidic hydrogen atoms of the alkyne.43–47 The multiple C–H···F interactions between the C2H2 molecule and inorganic MF62− or NbOF52− groups of MOFs have resulted in an efficient C2H2/CO2 separation.47–49 Thus the fluorinated MOFs are a highly appealing platform for C2H2/CO2 separation. Unfortunately, due to the synthetic inaccessibility of the suitable fluorinated organic ligands, the fluorinated MOFs have thus been largely unexplored.50–54 Herein, we present two isostructural fluorous MOFs [termed as JXNU-11(Fe2M), M represents the divalent metal ions of Ni2+ and Co2+] with a fluorinated organic linker, featuring unique cage-in-cage structures. Remarkably, JXNU-11(Fe2Ni) shows a record-high C2H2-capture amount and a long C2H2 breakthrough time in the practice breakthrough experiment, resulting in high-performance C2H2/CO2 separation. The double-shell cages possessing highly electronegative F and O atoms in JXNU-11(Fe2M) are suitable pore cages to match the size of C2H2. The C2H2 molecules are trapped in the highly fluorophilic and polarized double-shell cages through strong host–guest interactions, which are responsible for the highly efficient C2H2/CO2 separation by JXNU-11(Fe2M). Experimental Section Synthesis of JXNU-11(Fe2Ni) A mixture of FeCl3·6H2O (8.10 mg, 0.03 mmol), Ni(NO3)3·6H2O (4.36 mg, 0.015 mmol), 3,3′,5,5′-tetrakis(fluoro)biphenyl-4,4′-dicarboxylate acid (4.71 mg, 0.015 mmol), 1,3,5-tris(4-carboxyphenyl)benzene (6.58 mg, 0.015 mmol), N,N-dimethylformamide (DMF) (2 mL), and CH3COOH (0.19 mL) was capped in a 20 mL vial and heated at 100 °C for 2 days. After cooling to room temperature, brown triangle-shaped crystals were obtained. Elemental analysis for {[Fe2Ni(μ3-O)(TFBPDC)(BTB)4/3(H2O)3]·6DMF·2.5H2O}n (C68H77F4O24.5N6Fe2Ni: 1616.74). Calcd/found: H, 4.80/4.71; C, 50.51/50.41; N, 5.19/5.02. IR data (KBr, cm−1): 3427 (m), 1662 (s), 1615 (s), 1394 (s), 1253 (w), 1185 (w), 1147 (w), 1098 (w), 1035 (s), 858 (s), 807 (w), 781 (s), 683 (w), 663 (w), 589 (m), 498 (m), 443 (s). Synthesis of JXNU-11(Fe2Co) A mixture of FeCl3·6H2O (8.10 mg, 0.03 mmol), Co(NO3)3·6H2O (4.37 mg, 0.015 mmol), 3,3′,5,5′-tetrakis(fluoro)biphenyl-4,4′-dicarboxylate acid (4.71 mg, 0.015 mmol), 1,3,5-tris(4-carboxyphenyl)benzene (6.58 mg, 0.015 mmol), DMF (2 mL), and CH3COOH (0.17 mL) was capped in a 20 mL vial and heated at 100 °C for 2 days. After cooling to room temperature, brown triangle-shaped crystals were obtained. Elemental analysis for {[Fe2Co(μ3-O)(TFBPDC)(BTB)4/3(H2O)3]·4DMF·H2O}n (C62H60F4O21N4Fe2Co: 1443.77). Calcd/found: H, 4.18/4.21; C, 51.57/51.46; N, 3.88/3.92. IR data (KBr, cm−1): 3420 (m), 1655 (s), 1612 (s), 1394 (s), 1254 (w), 1185 (w), 1146 (w), 1098 (w), 1034 (s), 859 (s), 808 (w), 781 (s), 706 (w), 663 (w), 589 (m), 500 (m), 443 (s). Methods X-ray single-crystal diffraction experiments were carried out with a Rigaku Oxford SuperNova diffractometer, and powder X-ray diffraction patterns were recorded on a Rigaku DMax 2500 powder diffractometer. Gas sorption–desorption isotherms were measured on a Micromeritics ASAP 2020 HD88 adsorption analyzer. Breakthrough experiments for separation of C2H2/CO2 (v/v, 50/50) were carried out in a fixed bed with a gas chromatograph detection system. The detailed experimental methods are provided in Supporting Information. Results and Discussion Compounds JXNU-11(Fe2M) constructed from the oxygen-centered heterometallic trimeric [Fe2M(μ3-O)(COO)6] (M = Ni2+ and Co2+) clusters (Figure 1a) were prepared. Compound JXNU-11(Fe2M) formulated as [Fe2M(μ3-O)(TFBPDC)(BTB)4/3(H2O)3]n is based on the linear 3,3′,5,5′-tetrakis(fluoro)biphenyl-4,4′-dicarboxylate (TFBPDC2−) and the triangular 1,3,5-tris(4-carboxyphenyl)benzene (BTB3−) ligands (Figure 1b) and was characterized by single-crystal X-ray diffraction ( Supporting Information Table S2). The phase purity of the bulk samples was confirmed by powder X-ray diffraction ( Supporting Information Figure S2). The metal concentrations in both compounds were determined by inductively coupled plasma atomic emission spectroscopy ( Supporting Information Table S3). The two compounds are isostructural and crystallize in a trigonal R3m space group. The trimeric [Fe2M(μ3-O)] core is coordinated by two TFBPDC2−, four BTB3−, and three water ligands ( Supporting Information Figure S1). Each [Fe2M(μ3-O)(COO)6] cluster is linked by four BTB3− ligands and two TFBPDC2− ligands to generate a three-dimensional (3D) framework. Exploration of the structure indicates the 3D framework contains two kinds of cages: a small octahedral cage formed from six [Fe2M(μ3-O)(COO)6] clusters, three TFBPDC2−, and four BTB3− ligands with a diameter of ∼12 Å (Figure 1c), and a large cage composed of 12 [Fe2M(μ3-O)(COO)6] clusters, six TFBPDC2−, and ten BTB3− ligands with a diameter of ∼26 Å, which can be described as a cuboctahedron (Figure 1d). Thus the 3D framework has close packing of the octahedral and cuboctahedral cages (Figure 1e). Such large voids in the cages in a single 3D framework lead to the generation of a twofold interpenetrated network (Figure 1f). It is remarkable that each octahedral cage is encapsulated by a large cuboctahedral cage from the interpenetrating network (Figure 1g). Such a unique cage-in-cage structure in a MOF is distinctive. The double-shell nested cages are reminiscent of the interesting Russian-doll-like cage. In addition, because of the presence of the highly polar C–F and C–O bonds and dense F atoms on the pore surfaces for these cages, double-shell nested cages exhibit highly polarized and fluorophilic character. As shown in Figure 1g, part of the voids in the large cuboctahedral cage are occupied by the small octahedral cage. The remaining pores of the total 3D framework are filled by the disordered guest solvent molecules, which occupy 64.1% and 65.6% of the volumes of unit cell for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co), respectively. JXNU-11(Fe2M) materials exhibit excellent chemical stability after exposure to air for a long time or immersion in water with pH values ranging from 3 to 11 for 24 h ( Supporting Information Figure S3). The solvent-free JXNU-11(Fe2M) was obtained through solvent-exchange and heat under vacuum ( Supporting Information Figure S4). Figure 1 | Structures of JXNU-11(Fe2M). (a) [Fe2M(μ3-O)(COO)6] unit and (b) TFBPDC2− and BTB3− ligands. (c) 3D framework of JXNU-11(Fe2M). (d) Octahedral and (e) cuboctahedral cages. (f) Twofold interpenetrating frameworks of JXNU-11(Fe2M). (g) Cage-in-cage structure in JXNU-11(Fe2M). The fluorine atoms in (b–d) are represented as green balls. Download figure Download PowerPoint The permanent porosity of JXNU-11(Fe2M) was confirmed by N2 adsorption isotherms at 77 K. The N2 sorption isotherms for both compounds have typical type I sorption behavior with the saturation sorption amounts of 561 and 556 cm3 g−1 (Figure 2a). The pore volumes for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) from the experimental N2 data are 0.87 and 0.86 cm3 g−1, respectively, which are close to the corresponding crystal structure calculated values of 0.87 and 0.89 cm3 g−1. Additionally, the pore sizes are mainly around 6.8 and 11 Å (Figure 2a, inset), in good agreement with the obtained values from the crystal structures. The Brunauer–Emmett–Teller surface areas for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) were calculated to be 2208 m2 g−1 (Langmuir surface area of 2418 m2 g−1) and 2122 m2 g−1 (Langmuir surface area of 2392 m2 g−1) ( Supporting Information Figure S5), respectively. Figure 2 | Adsorption data of JXNU-11(Fe2M). (a) N2 adsorption–desorption isotherms at 77 K and pore size distribution for JXNU-11(Fe2M). C2H2 and CO2 adsorption–desorption isotherms of JXNU-11(Fe2M) at 273 K (b) and 298 K (c). Download figure Download PowerPoint The C2H2 and CO2 sorption isotherms of JXNU-11(Fe2M) (Figures 2b and 2c) were collected to evaluate the gas separation potentials for JXNU-11(Fe2M). JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) exhibit high C2H2 adsorption capacities of 191 and 180 cm3 g−1 at 273 K and 1 bar, respectively, which are higher than the best-performance MOFs, including FJU-6-TATB (160 cm3 g−1),55 ATU-Cu (134 cm3 g−1),33 ZJU-74a (107 cm3 g−1),42 JXNU-10(Y) (94.9 cm3 g−1),27 CuI@UIO-66-(COOH)2 (71 cm3 g−1),56 and MOF-OH (68.7 cm3 g−1).57 At 298 K, the adsorption amounts of C2H2 are 118 and 107 cm3 g−1 for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co), respectively, which are larger than those of the prominent MOFs FeNi-M′MOF (96 cm3 g−1),58 JNU-1 (63 cm3 g−1),59 NKMOF-1-Ni (61 cm3 g−1),60 UTSA-300a (69 cm3 g−1),47 and CPL-NH2 (41 cm3 g−1).39 In contrast, JXNU-11(Fe2M) showed low CO2 uptakes (Figures 2b and 2c). The high uptake of C2H2 is mainly caused by the polarized and fluorophilic microporous environments made of the cage-in-cage structures that are decorated with plenty of highly electronegative fluorine and carboxylate oxygen atoms that have high affinity for the acidic H atoms of C2H2 molecules but low affinity for CO2 molecules with two terminal electronegative O atoms. Such distinct gas sorption behaviors were confirmed by grand canonical Monte Carlo (GCMC) simulations. As depicted in Figures 3a and 3b, the distribution of C2H2 molecules in JXNU-11(Fe2M) is mainly within the octahedral cage and the space between the two cages, whereas much fewer CO2 molecules were distributed within the double-shell nest. Such a phenomenon suggests the double-shell nest is the desirable space for accommodating C2H2 molecules, in agreement with the experimental finding. The linear adsorption isotherms for CO2 further suggest the low affinity to CO2 for the host frameworks. The different adsorption behaviors between C2H2 and CO2 were also evidenced by the isosteric heats of adsorption (Qst) ( Supporting Information Figure S6). The obtained Qst of C2H2 were in the range of 31.6–29.7 kJ mol−1 for JXNU-11(Fe2M), which are notably higher than those of CO2 (16.0–19.6 kJ mol−1), further reflecting the strong affinity toward C2H2. The present Qst(C2H2) values are comparable to those of FJU-6-TATB (29 kJ mol−1),55 UTSA-74 (31 kJ mol−1),37 and JXNU-5a (32.9 kJ mol−1),35 but lower than those of CuI@UIO-66-(COOH)2 (74.5 kJ mol−1),56 JNU-1 (47.6 kJ mol−1),59 and NKMOF-1-Ni (53.9 kJ mol−1).60 The evident difference in adsorption enthalpy endows JXNU-11(Fe2M) with the thermodynamic separation possibility of a C2H2/CO2 mixture. Furthermore, JXNU-11(Fe2M) retained C2H2 storage ability after adsorption/desorption cycling experiments ( Supporting Information Figure S7), confirming their excellent recyclability for C2H2 adsorption. Figure 3 | Grand canonical Monte Carlo (GCMC) adsorption simulation. Computational C2H2 (a) and CO2 (b) distribution in the double-shell nest in JXNU-11(Fe2M) at 298 K and 1 atm. Download figure Download PowerPoint To evaluate the separation performance of JXNU-11(Fe2M), ideal adsorbed solution theory was used to calculate the separation selectivity of the C2H2/CO2 (v∶v = 50∶50) mixture. At 1 bar and 298 K, the adsorption selectivities for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) are 2.7 and 2.5, respectively ( Supporting Information Figure S8), which are comparable to the leading MOFs of UPC-200(Fe)-F-H2O (2.25)32 and FJU-6-TATB (3.1).55 The C2H2/CO2 uptake ratios for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) at 298 K and 1 bar are 2.1 and 1.9, respectively, which are larger than those of noted MOFs FJU-90a (1.75),36 SIFSIX-Cu-TPA (1.7),49 and ATU-Cu (1.2).33The high C2H2 adsorption capacity and moderate selectivity of JXNU-11(Fe2M) at ambient temperature reveal the greatly promising potential for C2H2/CO2 separation. To further investigate the C2H2/CO2 separation performance of JXNU-11(Fe2M), the breakthrough experiments were performed in practical separation processes for a C2H2/CO2 mixture (v∶v = 50∶50) at 298 K and 1 atm. The breakthrough experimental results show that CO2 was eluted first from the packed column, whereas C2H2 was retained in the bed for more than 120 and 108 min g−1 for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co) (Figure 4a), respectively. The C2H2 breakthrough times for JXNU-11(Fe2M) rank among the leading MOFs ( Supporting Information Table S4). Moreover, the significant roll-up phenomenon of CO2 was observed in the breakthrough experiments, indicating a large proportion of previously adsorbed CO2 molecules can be replaced by the later fed C2H2 molecules. Such a phenomenon means that the C2H2 molecule competes better than CO2 for the binding sites of JXNU-11(Fe2M), highlighting the excellent separation efficiency of JXNU-11(Fe2M) for a C2H2/CO2 mixture. The breakthrough experiments were terminated when the concentrations of effluent gases stabilized. Figure 4 | C2H2 and CO2 separation performances. (a) Breakthrough curves of JXNU-11(Fe2M) at the flow of 2 mL min−1, 298 K, and 1 atm. Breakthrough curves of JXNU-11(Fe2Ni) at different temperatures (b) and different total flow rates at 298 K (c). Download figure Download PowerPoint The captured amounts for C2H2 during the 0∼tbreak time under the dynamic conditions are 4.8 and 4.3 mmol g−1 for JXNU-11(Fe2Ni) and JXNU-11(Fe2Co), respectively. The C2H2 captured amounts are remarkable and outperform all other MOFs, including the top-performing ZJU-74a (3.64 mmol g−1),42 SNNU-45 (3.5 mmol g−1),41 BSF-3 (2.9 mmol g−1),38 JCM-1 (2.2 mmol g−1),61 and FJU-90a (1.87 mmol g−1)36 (Figure 5 and Supporting Information Table S4). In addition, a long interval time of 55 min g−1 between C2H2 and CO2 breakthrough for JXNU-10(Fe2Ni) was obtained. Such an interval time stands out amongst MOFs and is only shorter than those of SNNU-45 (79 min g−1)41 and SIFSIX-Cu-TPA (69 min g−1).49 Furthermore, their practical separation potentials under different temperatures were evaluated. The C2H2 breakthrough times increased significantly with a decrease of temperature (Figure 4b and Supporting information Figure S9). The breakthrough times of C2H2 reached 146 (283 K) and 165 min g−1 (273 K) with a gas flow of 2 mL min−1. Figure 5 | Experimental breakthrough performance. Comparison of MOFs with top-high C2H2/CO2 breakthrough performance at 298 K and 1 atm. Download figure Download PowerPoint The separation performance of JXNU-11(Fe2M) with different gas flow rates at 298 K were also studied. With the gas flow of 2 mL min−1, the longest C2H2 breakthrough time and the largest amount of the adsorbed C2H2 were achieved for JXNU-11(Fe2M) (Figure 4c and Supporting Information Figure S10). JXNU-11(Fe2M) exhibited no appreciable changes in breakthrough times after three cycles of dynamic breakthrough experiments ( Supporting Information Figure S11), indicative of an excellent recycling separation capability. With the advantage of the modest adsorption heat of C2H2, JXNU-11(Fe2M) materials can be easily regenerated through purging the column with He gas at ambient temperature, as evidenced by the overlapped breakthrough curves ( Supporting Information Figure S11). Such results further indicate JXNU-11(Fe2M) are highly desirable microporous materials for C2H2/CO2 separation. Finally, the desorption experiments were carried out on JXNU-11(Fe2Ni). After reaching the breakthrough equilibrium, the adsorption column was purged with He gas (flow rate: 4 mL min−1) under ambient pressure and 303 K. As depicted in Figure 6a, the adsorbed CO2 molecules were released from the adsorption bed quickly. In contrast, the desorption of C2H2 gas was much slower. The productivity of C2H2 with over 95% purity was estimated from the desorption curve to be 1.98 mmol for 1 g of JXNU-11(Fe2Ni) (Figure 6b). Figure 6 | Breakthrough experiments. (a) and (b) Breakthrough curves and desorption curves of JXNU-11(Fe2Ni) based on the equimolar C2H2/CO2 mixture. Download figure Download PowerPoint To fully understand the underlying mechanism of the preferential adsorption of C2H2 over CO2 in JXNU-11(Fe2M), we performed comprehensive GCMC simulations to investigate the gas adsorption sites for JXNU-11(Fe2M). The calculated preferential binding sites of C2H2 are shown in Figures 7a and 7b. Four primary binding sites (sites I–IV) for C2H2 were found. As expected, the primary binding sites for C2H2 molecules are F and carboxylate O sites. For site I, site II, and site III, strong C−H···F hydrogen bonding notably occur between C2H2 and F atoms (H···F = 2.49–3.01 Å), which confirms the strong binding affinity of JXNU-11(Fe2M) toward C2H2. Compared with these sites, site IV exhibits strong C−H⋯O hydrogen bonds with a pair of carboxylate O atoms with the distances of 2.61 and 2.60 Å. These H⋯F and H⋯O distances are shorter than the sums of the corresponding van der Waals radii of hydrogen and fluorine (2.67 Å) or hydrogen and oxygen (2.72 Å) atoms, suggesting the substantial interactions. Moreover, the C≡C group of C2H2 interacts with the adjacent open metal sites through π⋯M interactions, and weakly interacts with the neighboring aromatic ring units of the framework through π⋯π interactions. Thus these multiple interactions between C2H2 molecules and the framework synergistically result in the remarkable affinity for C2H2. Therefore, the nanosized double-shell nests with many highly electronegative F and O atoms provide the polarized and fluorophilic pore spaces to accommodate C2H2. In contrast, the negatively charged F and O sites are not desirable sites for CO2. The presence of the repulsive interactions between the F/O atoms and the O atoms of CO2 leads to a weaker interaction between CO2 and the framework, which is consistent with the experimental results. The computational results gave an average binding energy of 31.8 kJ mol−1 for C2H2, which is much higher than that of 20.9 kJ mol−1 for CO2, further corroborating the stronger affinity toward C2H2 in the present framework. These calculated results clearly reveal that the rich F and O atoms as well as the open metal sites in the cage-in-cage structures are synergistically responsible for the preferential adsorption of C2H2 over CO2, thus affording the top-level C2H2/CO2 separation performance of JXNU-11(Fe2M). Figure 7 | Mechanism study. adsorption sites of C2H2 in JXNU-11(Fe2M). (a) I and (b) and Ni F and O in and H in C2H2. The labeled is measured in Å. Download figure Download PowerPoint two MOFs M = Ni and featuring nanosized structures and highly polarized and fluorophilic pore Such a cage-in-cage structure formed is composed of an and an large The double-shell nested cages in the MOFs provide the fascinating and the optimized pore spaces for trapping C2H2. The abundant F and O atoms of the organic ligands in JXNU-11(Fe2M) JXNU-11(Fe2M) with multiple binding sites for C2H2. The F and O atoms possessing high in JXNU-11(Fe2M) the pore surfaces of the cages result in the highly dense on the pore which are the desirable sites for trapping C2H2 with charged H atoms. As a JXNU-11(Fe2M) shows a strong binding affinity to C2H2 over CO2. The nanosized double-shell cages in JXNU-11(Fe2M) afford the optimized pore spaces for trapping C2H2, leading to an efficient C2H2/CO2 separation. The practical breakthrough performance of JXNU-11(Fe2Ni) the breakthrough time for C2H2 gas and the C2H2-capture the highest observed at ambient confirming the C2H2/CO2 separation performance of JXNU-11(Fe2M) Thus we a strategy for the design of cage-in-cage structures in MOFs and develop an efficient for the of the highly polarized and fluorophilic pore cages in the fluorous Supporting Information Supporting Information is and detailed experimental crystal data for JXNU-11(Fe2M), adsorption column breakthrough experiments, and GCMC simulation. of is no of to Information This work is by the of and and the for and of of Jiangxi of and of Acetylene and by the of K. Chemistry and of Metal–Organic and of Frameworks to MOF and on Adsorption in a as the in the of the Chemistry of and Metal–Organic Frameworks and Zhang of Frameworks Featuring Separation Metal–Organic Frameworks and for of Surfaces by in Synthesis of in a

Water plays an important role in mediating protein-ligand interactions. Water rearrangement upon a ligand binding or modification can be very slow and beyond typical timescales used in molecular dynamics (MD) simulations. Thus, inadequate sampling of slow water motions in MD simulations often impairs the accuracy of the accuracy of ligand binding free energy calculations. Previous studies suggest grand canonical Monte Carlo (GCMC) outperforms normal MD simulations for water sampling, thus GCMC has been applied to help improve the accuracy of ligand binding free energy calculations. However, in prior work we observed protein and/or ligand motions impaired how well GCMC performs at water rehydration, suggesting more work is needed to improve this method to handle water sampling. In this work, we applied GCMC in 21 protein-ligand systems to assess the performance of GCMC for rehydrating buried water sites. While our results show that GCMC can rapidly rehydrate all selected water sites for most systems, it fails in five systems. In most failed systems, we observe protein/ligand motions, which occur in the absence of water, combine to close water sites and block instantaneous GCMC water insertion moves. For these five failed systems, we both extended our GCMC simulations and tested a new technique named grand canonical nonequilibrium candidate Monte Carlo (GCNCMC). GCNCMC combines GCMC with the nonequilibrium candidate Monte Carlo (NCMC) sampling technique to improve the probability of a successful water insertion/deletion. Our results show that GCNCMC and extended GCMC can rehydrate all target water sites for three of the five problematic systems and GCNCMC is more efficient than GCMC in two out of the three systems. In one system, only GCNCMC can rehydrate all target water sites, while GCMC fails. Both GCNCMC and GCMC fail in one system. This work suggests this new GCNCMC method is promising for water rehydration especially when protein/ligand motions may block water insertion/removal.

Pure and binary component sorption equilibria of light hydrocarbons in the zeolite silicalite from grand canonical Monte Carlo simulations

The adsorption of nitrogen, oxygen and argon has been studied in cadmium (II) ions exchanged zeolite A at 288.0 and 303.0 K. Experimentally measured adsorption isotherms are compared with theoretically calculated data using grand canonical Monte Carlo (GCMC) simulation. Nitrogen showed higher adsorption capacity and selectivity than oxygen and argon in these zeolite samples. The cadmium exchanged zeolite A showed increased adsorption capacity for nitrogen, oxygen, and argon with increase in cadmium (II) exchange levels. Isosteric heat of adsorption data showed stronger interactions of nitrogen molecules with cadmium (II) cations in zeolite samples. These observations have been explained in terms of higher electrostatic interaction of nitrogen with extra framework zeolite cations. The selectivity of oxygen over argon is explained in terms of its higher interaction of oxygen with cadmium exchanged zeolites than argon molecules. Heats of adsorption and adsorption isotherms were also calculated using grand canonical Monte Carlo simulation algorithm. Simulation studies expectedly show the proximity of nitrogen molecules to the locations of extra framework sodium and cadmium cations.

Grand canonical Monte Carlo (GCMC) simulation combined with ab initio quantum mechanics calculations were employed to study methane storage in homogeneous armchair open-ended single-walled silicon nanotubes (SWSiNTs), single-walled carbon nanotubes (SWCNTs), and single-walled silicon carbide nanotubes (SWSiCNTs) in triangular arrays. Two different groups of nanotubes were studied: the first were (12,12) SiNTs, (19,19) CNTs, and (15,15) SiCNTs and the second were (7,7) SiNTs, (11,11) CNTs, and (9,9) SiCNTs with the diameters of 26 and 15 A for the first and second groups, respectively. The simulations were carried out at different thermodynamic states. The potential energy functions were calculated using ab initio quantum mechanics and then fitted with (12,6) Lennard–Jones potential model as a bridge between first-principles calculations and GCMC simulations. The absolute, excess, and delivery adsorption isotherms of methane were calculated for two groups of nanotubes. The specific surface area and the isosteric heat of adsorption were computed. The radial distribution functions for the adsorbed molecules on different nanotubes were also calculated. Different isotherm models were fitted with the simulation adsorption data. According to the results, the excess uptake value of methane adsorption in (11,11) CNT array exceeded the US Department of Energy target (180 V/V at 298 K and 35 bar). The results also indicate that SiNTs and SiCNTs are not desirable materials compared with corresponding CNTs for natural gas storage.