Understanding shape selectivity effects of hydroisomerization using a reaction equilibrium model.

We study important aspects of shape selectivity effects of zeolites for hydroisomerization of linear alkanes, which produces a myriad of isomers, particularly for long chain hydrocarbons. To investigate the conditions for achieving an optimal yield of branched hydrocarbons, it is important to understand the role of chemical equilibrium in these reversible reactions. We conduct an extensive analysis of shape selectivity effects of different zeolites for the hydroisomerization of C7 and C8 isomers at chemical reaction equilibrium conditions. The reaction ensemble Monte Carlo method, coupled with grand-canonical Monte Carlo simulations, is commonly used for computing reaction equilibrium of heterogeneous reactions. The computational demands become prohibitive for a large number of reactions. We used a faster alternative in which reaction equilibrium is obtained by imposing chemical equilibrium in the gas phase and phase equilibrium between the gas phase components and the adsorbed phase counterparts. This effectively mimics the chemical equilibrium distribution in the adsorbed phase. Using Henry's law at infinite dilution and mixture adsorption isotherm models at elevated pressures, we calculate the adsorbed loadings in the zeolites. This study shows that zeolites with cage or channel-like structures exhibit significant differences in selectivity for alkane isomers. We also observe a minimal impact of pressure on the gas-phase equilibrium of these reactions at typical experimental reaction temperatures 400-700K. This study marks initial strides in understanding the reaction product distribution for long-chain alkanes.

The Influence of UiO-66 Metal-Organic Framework Structural Defects on Adsorption and Separation of Hexane Isomers.

In this work, adsorption properties of the UiO-66 metal-organic framework were investigated, with particular emphasis on the influence of structural defects. A series of UiO-66 samples were synthesized and characterized using a wide range of experimental techniques. Type I adsorption isotherms for low-temperature adsorption of N2 and Ar showed that micropore volume and specific surface area significantly increase with the number of defects. Adsorption of hexane isomers in UiO-66 was studied by means of quasi-equilibrated temperature-programmed desorption and adsorption (QE-TPDA) experimental and Monte Carlo simulation techniques. QE-TPDA profiles revealed that only defect-free UiO-66 exhibits distinct two adsorption states. This technique also yielded high-quality adsorption isobars that were successfully recreated using Grand-Canonical Monte Carlo molecular simulations, which, however, required refinement of the existing force fields. The calculations demonstrated the detailed mechanism of adsorption and separation of hexane isomers in the UiO-66 structure. The preferred tetrahedral cages provide suitable voids for bulky molecules, which is the reason for unusual "reverse" selectivity of UiO-66 towards di-branched alkanes. Interconnection of the tetrahedral cavities due to missing organic linkers greatly reduces the selectivity of the defected material.

New Features of the Open Source Monte Carlo Software Brick-CFCMC: Thermodynamic Integration and Hybrid Trial Moves.

We present several new major features added to the Monte Carlo (MC) simulation code Brick-CFCMC for phase- and reaction equilibria calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one is thermodynamic integration for the computation of excess chemical potentials (µex). For this purpose, we implemented the computation of the ensemble average of the derivative of the potential energy with respect to the scaling factor for intermolecular interactions (⟨∂U∂λ⟩). Efficient bookkeeping is implemented so that the quantity ∂U∂λ is updated after every MC trial move with negligible computational cost. We demonstrate the accuracy and reliability of the calculation of µex for sodium chloride in water. Second, we implemented hybrid MC/MD translation and rotation trial moves to increase the efficiency of sampling of the configuration space. In these trial moves, short Molecular Dynamics (MD) trajectories are performed to collectively displace or rotate all molecules in the system. These trajectories are accepted or rejected based on the total energy drift. The efficiency of these trial moves can be tuned by changing the time step and the trajectory length. The new trial moves are demonstrated using MC simulations of a viscous fluid (deep eutectic solvent).

Effect of Water Content on Thermodynamic Properties of Compressed Hydrogen.

Force field-based molecular simulations were used to calculate thermal expansivities, heat capacities, and Joule-Thomson coefficients of binary (standard) hydrogen-water mixtures for temperatures between 366.15 and 423.15 K and pressures between 50 and 1000 bar. The mole fraction of water in saturated hydrogen-water mixtures in the gas phase ranges from 0.004 to 0.138. The same properties were calculated for pure hydrogen at 323.15 K and pressures between 100 and 1000 bar. Simulations were performed using the TIP3P and a modified TIP4P force field for water and the Marx, Vrabec, Cracknell, Buch, and Hirschfelder force fields for hydrogen. The vapor-liquid equilibria of hydrogen-water mixtures were calculated along the melting line of ice Ih, corresponding to temperatures between 264.21 and 272.4 K, using the TIP3P force field for water and the Marx force field for hydrogen. In this temperature range, the solubilities and the chemical potentials of hydrogen and water were obtained. Based on the computed solubility data of hydrogen in water, the freezing-point depression of water was computed ranging from 264.21 to 272.4 K. The modified TIP4P and Marx force fields were used to improve the solubility calculations of hydrogen-water mixtures reported in our previous study [Rahbari A.;J. Chem. Eng. Data2019, 64, 4103-4115] for temperatures between 323 and 423 K and pressures ranging from 100 to 1000 bar. The chemical potentials of ice Ih were calculated as a function of pressure between 100 and 1000 bar, along the melting line for temperatures between 264.21 and 272.4 K, using the IAPWS equation of state for ice Ih. We show that at low pressures, the presence of water has a large effect on the thermodynamic properties of compressed hydrogen. Our conclusions may have consequences for the energetics of a hydrogen refueling station using electrochemical hydrogen compressors.

Synthesis of Chiral MOF-74 Frameworks by Post-Synthetic Modification by Using an Amino Acid.

The synthesis of chiral metal-organic frameworks (MOFs) is highly relevant for asymmetric heterogenous catalysis, yet very challenging. Chiral MOFs with MOF-74 topology were synthesised by using post-synthetic modification with proline. Vibrational circular dichroism studies demonstrate that proline is the source of chirality. The solvents used in the synthesis play a key role in tuning the loading of proline and its interaction with the MOF-74 framework. In N,N'-dimethylformamide, proline coordinates monodentate to the Zn2+ ions within the MOF-74 framework, whereas it is only weakly bound to the framework when using methanol as solvent. Introducing chirality within the MOF-74 framework also leads to the formation of defects, with both the organic linker and metal ions missing from the framework. The formation of defects combined with the coordination of DMF and proline within the framework leads to a pore blocking effect. This is confirmed by adsorption studies and testing of the chiral MOFs in the asymmetric aldol reaction between acetone and para-nitrobenzaldehyde.

Brick-CFCMC: Open Source Software for Monte Carlo Simulations of Phase and Reaction Equilibria Using the Continuous Fractional Component Method.

We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.

Potential of polarizable force fields for predicting the separation performance of small hydrocarbons in M-MOF-74.

The separation of light olefins from paraffins via cryogenic distillation is a very energy intensive process. Solid adsorbents and especially metal-organic frameworks with open metal sites have the potential to significantly lower the required energy. Specifically, M-MOF-74 has drawn considerable attention for application in olefin/paraffin separation. To investigate how the separation proceeds on a molecular level and to design better materials, molecular simulation can be a useful tool. Unfortunately, it is still a challenge to model the adsorption behavior of many adsorbates in metal-organic frameworks with open metal sites. Previously, the inclusion of explicit polarization has been suggested to improve the quality of classical force fields for such systems. Here, the potential of polarizable force fields for the description of olefins and paraffins in metal-organic frameworks with open metal sites is investigated. In particular, heats of adsorption, binding geometries, and adsorption isotherms are calculated for C2H4, C2H6, C3H6, and C3H8 in M-MOF-74 (with M = Co, Mn, Fe, and Ni). In this study, no force field parameters are adjusted to improve the model. The results show that including explicit polarization significantly improves the description of the adsorption in comparison to non-polarizable generic force fields which do not consider explicit polarization. The study also reveals that simulation predictions are sensitive to the assigned repulsive potential and framework charges. A fully re-parametrized polarizable force field may have the capability to improve the predictions even further.

Predicting Multicomponent Adsorption Isotherms in Open-Metal Site Materials Using Force Field Calculations Based on Energy Decomposed Density Functional Theory.

For the design of adsorptive-separation units, knowledge is required of the multicomponent adsorption behavior. Ideal adsorbed solution theory (IAST) breaks down for olefin adsorption in open-metal site (OMS) materials due to non-ideal donor-acceptor interactions. Using a density-function-theory-based energy decomposition scheme, we develop a physically justifiable classical force field that incorporates the missing orbital interactions using an appropriate functional form. Our first-principles derived force field shows greatly improved quantitative agreement with the inflection points, initial uptake, saturation capacity, and enthalpies of adsorption obtained from our in-house adsorption experiments. While IAST fails to make accurate predictions, our improved force field model is able to correctly predict the multicomponent behavior. Our approach is also transferable to other OMS structures, allowing the accurate study of their separation performances for olefins/paraffins and further mixtures involving complex donor-acceptor interactions.

Aqueous Solutions of Ionic Liquids: Microscopic Assembly.

Aqueous solutions of ionic liquids are of special interest, due to the distinctive properties of ionic liquids, in particular, their amphiphilic character. A better understanding of the structure-property relationships of such systems is hence desirable. One of the crucial molecular-level interactions that influences the macroscopic behavior is hydrogen bonding. In this work, we conduct molecular dynamics simulations to investigate the effects of ionic liquids on the hydrogen-bond network of water in dilute aqueous solutions of ionic liquids with various combinations of cations and anions. Calculations are performed for imidazolium-based cations with alkyl chains of different lengths and for a variety of anions, namely, [Br](-), [NO3](-), [SCN](-) [BF4](-), [PF6](-), and [Tf2N](-). The structure of water and the water-ionic liquid interactions involved in the formation of a heterogeneous network are analyzed by using radial distribution functions and hydrogen-bond statistics. To this end, we employ the geometric criterion of the hydrogen-bond definition and it is shown that the structure of water is sensitive to the amount of ionic liquid and to the anion type. In particular, [SCN](-) and [Tf2N](-) were found to be the most hydrophilic and hydrophobic anions, respectively. Conversely, the cation chain length did not influence the results.

Exploiting Large-Pore Metal-Organic Frameworks for Separations through Entropic Molecular Mechanisms.

We review the molecular mechanisms behind adsorption and the separations of mixtures in metal-organic frameworks and zeolites. Separation mechanisms can be based on differences in the affinity of the adsorbate with the framework and on entropic effects. To develop next-generation adsorbents, the separation efficiency of the materials needs to be improved. The performance under industrially relevant conditions largely depends on two factors: 1) the separation selectivity and 2) the pore volume capacity of the material. Enthalpic mechanisms can lead to increased selectivities, but these are mostly restricted to the low loading regime, and hence these mechanisms are unable to make use of all of the large-pore volume that a metal-organic framework can provide. Industrial processes routinely operate in the pore saturation regime. In this Review, we focus on entropic molecular separation mechanisms that are effective under these conditions and, in particular, on a recent methodology to obtain high selectivities at high pore loading.

Entropic separations of mixtures of aromatics by selective face-to-face molecular stacking in one-dimensional channels of metal-organic frameworks and zeolites.

Separation of challenging mixtures using metal-organic frameworks can be achieved by an entropy-driven mechanism, where one of the components can arrange into a "face-to-face" stacking, thus reducing its "footprint" and reaching a higher saturation loading.

Separation of Amyl Alcohol Isomers in ZIF-77.

The separation of pentanol isomer mixtures is shown to be very efficient using the nanoporous adsorbent zeolitic imidazolate framework ZIF-77. Through molecular simulations, we demonstrate that this material achieves a complete separation of linear from monobranched-and these from dibranched-isomers. Remarkably, the adsorption and diffusion behaviors follow the same decreasing trend, produced by the channel size of ZIF-77 and the guest shape. This separation based on molecular branching applies to alkanes and alcohols and promises to encompass numerous other functional groups.

Entropic separation of styrene/ethylbenzene mixtures by exploitation of subtle differences in molecular configurations in ordered crystalline nanoporous adsorbents.

The separation of styrene/ethylbenzene mixture is of great importance in the petrochemical industry. Current technology uses distillation; this separation is difficult because of the small, 9 K, difference in the boiling points. An alternative separation method uses selective adsorption in nanoporous materials such as zeolites and metal-organic frameworks. Here we present a simulation screening study for the separation of styrene/ethylbenzene mixture by adsorptive means in nanoporous materials near pore saturation conditions. Under these conditions, different entropic mechanisms can dictate the separation process. Commensurate stacking has the best trade-off between selectivity and saturation capacity and offers a geometrical solution to the separation problem. MIL-47 has the right channel size and topology for styrene to exhibit commensurate stacking offering high capacity and selectivity for styrene over ethylbenzene. Out of all the screened structures, MIL-47 was found to be the best candidate for the separation of styrene/ethylbenzene mixture.

Separating xylene isomers by commensurate stacking of p-xylene within channels of MAF-X8.

The development of energy-efficient processes for selective separation of p-xylene from mixtures with its isomers is of vital importance in the petrochemical industries. Current industrial practice uses BaX zeolite that has high adsorption selectivity for p-xylene. Finding para-selective structures is challenging. With state-of-the-art simulation methodologies we systematically screened a wide variety of zeolites and metal-organic frameworks (MOFs). Our investigations highlight the crucial importance of the channel dimension on the separation. MAF-X8 is particularly noteworthy because the channel dimensions and geometry allow "commensurate stacking" which we exploit as a separation mechanism at saturation conditions. Due to a significantly improved capacity compared to BaX, the cycle times for p-xylene with MAF-X8 are found to be about a factor of 4.5 longer. This is expected to result in significant process improvements.

Molecular-level insight into unusual low pressure CO2 affinity in pillared metal-organic frameworks.

Fundamental insight into how low pressure adsorption properties are affected by chemical functionalization is critical to the development of next-generation porous materials for postcombustion CO2 capture. In this work, we present a systematic approach to understanding low pressure CO2 affinity in isostructural metal-organic frameworks (MOFs) using molecular simulations and apply it to obtain quantitative, molecular-level insight into interesting experimental low pressure adsorption trends in a series of pillared MOFs. Our experimental results show that increasing the number of nonpolar functional groups on the benzene dicarboxylate (BDC) linker in the pillared DMOF-1 [Zn2(BDC)2(DABCO)] structure is an effective way to tune the CO2 Henry's coefficient in this isostructural series. These findings are contrary to the common scenario where polar functional groups induce the greatest increase in low pressure affinity through polarization of the CO2 molecule. Instead, MOFs in this isostructural series containing nitro, hydroxyl, fluorine, chlorine, and bromine functional groups result in little increase to the low pressure CO2 affinity. Strong agreement between simulated and experimental Henry's coefficient values is obtained from simulations on representative structures, and a powerful yet simple approach involving the analysis of the simulated heats of adsorption, adsorbate density distributions, and minimum energy 0 K binding sites is presented to elucidate the intermolecular interactions governing these interesting trends. Through a combined experimental and simulation approach, we demonstrate how subtle, structure-specific differences in CO2 affinity induced by functionalization can be understood at the molecular-level through classical simulations. This work also illustrates how structure-property relationships resulting from chemical functionalization can be very specific to the topology and electrostatic environment in the structure of interest. Given the excellent agreement between experiments and simulation, predicted CO2 selectivities over N2, CH4, and CO are also investigated to demonstrate that methyl groups also provide the greatest increase in CO2 selectivity relative to the other functional groups. These results indicate that methyl ligand functionalization may be a promising approach for creating both water stable and CO2 selective variations of other MOFs for various industrial applications.

Electro-osmotic Drag and Thermodynamic Properties of Water in Hydrated Nafion Membranes from Molecular Dynamics.

One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Molecular simulation investigation into the performance of Cu-BTC metal-organic frameworks for carbon dioxide-methane separations.

We report a molecular simulation study for Cu-BTC metal-organic frameworks as carbon dioxide-methane separation devices. For this study we have computed adsorption and diffusion of methane and carbon dioxide in the structure, both as pure components and mixtures over the full range of bulk gas compositions. From the single component isotherms, mixture adsorption is predicted using the ideal adsorbed solution theory. These predictions are in very good agreement with our computed mixture isotherms and with previously reported data. Adsorption and diffusion selectivities and preferential sitings are also discussed with the aim to provide new molecular level information for all studied systems.

Competitive Adsorption of Xylenes at Chemical Equilibrium in Zeolites.

The separation of xylenes is one of the most important processes in the petrochemical industry. In this article, the competitive adsorption from a fluid-phase mixture of xylenes in zeolites is studied. Adsorption from both vapor and liquid phases is considered. Computations of adsorption of pure xylenes and a mixture of xylenes at chemical equilibrium in several zeolite types at 250 °C are performed by Monte Carlo simulations. It is observed that shape and size selectivity entropic effects are predominant for small one-dimensional systems. Entropic effects due to the efficient arrangement of xylenes become relevant for large one-dimensional systems. For zeolites with two intersecting channels, the selectivity is determined by a competition between enthalpic and entropic effects. Such effects are related to the orientation of the methyl groups of the xylenes. m-Xylene is preferentially adsorbed if xylenes fit tightly in the intersection of the channels. If the intersection is much larger than the adsorbed molecules, p-xylene is preferentially adsorbed. This study provides insight into how the zeolite topology can influence the competitive adsorption and selectivity of xylenes at reaction conditions. Different selectivities are observed when a vapor phase is adsorbed compared to the adsorption from a liquid phase. These insight have a direct impact on the design criteria for future applications of zeolites in the industry. MRE-type and AFI-type zeolites exclusively adsorb p-xylene and o-xylene from the mixture of xylenes in the liquid phase, respectively. These zeolite types show potential to be used as high-performing molecular sieves for xylene separation and catalysis.

Correction: Role of additives and solvents in the synthesis of chiral isoreticular MOF-74 topologies.

Correction for 'Role of additives and solvents in the synthesis of chiral isoreticular MOF-74 topologies' by Andreea Gheorghe et al., Dalton Trans., 2021, DOI: 10.1039/D1DT01945G.

Role of additives and solvents in the synthesis of chiral isoreticular MOF-74 topologies.

Chiral induction is a simple and inexpensive approach to synthesise chiral metal-organic frameworks, even when using achiral building-blocks. The challenge lies in selecting the proper chiral inductor. This can only be achieved upon understanding the mechanism behind the chirality transfer from the chiral guest to the achiral MOF. In this work, the role of two types of chiral additives and different solvents was investigated in the crystallization of isoreticular MOF-74. We show that pyrrolidone-based solvents can interact with the framework walls and influence the thermal stability of the MOF. The role of the different chiral additives is related to the strength of their interaction with the MOF. Unlike cinchona alkaloids that have weak interactions with the framework, L- or D-trans-4-hydroxyproline (L- or D-Hyp) can strongly bind to the Zn2+ metal centres and cause the twisting of the organic linker. Moreover, L- and D-Hyp additives can affect the IRMOF-74 nucleation process depending on their concentration and handedness.