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.

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 Zn²⁺ 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.

Strongly interacting chiral additives play a complex role during the crystallisation of chiral isoreticular MOF-74. They can coordinate to the open sites of the metal ions and induce strain on the framework, leading to a local twisting of the organic linker.

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 (). Efficient bookkeeping is implemented so that the quantity 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.

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.

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.

Multiple Free Energy Calculations from Single State Point Continuous Fractional Component Monte Carlo Simulation Using Umbrella Sampling

We introduce an alternative method to perform free energy calculations for mixtures at multiple temperatures and pressures from a single simulation, by combining umbrella sampling and the continuous fractional component Monte Carlo method. One can perform a simulation of a mixture at a certain pressure and temperature and accurately compute the chemical potential at other pressures and temperatures close to the simulation conditions. This method has the following advantages: (1) Accurate estimates of the chemical potential as a function of pressure and temperature are obtained from a single state simulation without additional postprocessing. This can potentially reduce the number of simulations of a system for free energy calculations for a specific temperature and/or pressure range. (2) Partial molar volumes and enthalpies are obtained directly from the estimated chemical potentials. We tested our method for a Lennard-Jones system, aqueous mixtures of methanol at T = 298 K and P = 1 bar, and a mixture of ammonia, nitrogen, and hydrogen at T = 573 K and P = 800 bar. For pure methanol (N = 410 molecules), we observed that the estimated chemical potentials from umbrella sampling are in excellent agreement with the reference values obtained from independent simulations, for ΔT = ±15 K and ΔP = 100 bar (with respect to the simulated system). For larger systems, this range becomes smaller because the relative fluctuations of energy and volume become smaller. Without sufficient overlap, the performance of the method will become poor especially for nonlinear variations of the chemical potential.

Water-Ethanol and Methanol-Ethanol Separations Using in Situ Confined Polymer Chains in a Metal-Organic Framework

This study presents a straightforward approach for the in situ polymerization of poly(N-isopropylacrylamide) (PNIPAM) chains within the one-dimensional (1D) pores of the five-coordinated zinc-based metal-organic framework DMOF in order to obtain new MOF-based composites. The loading amount of PNIPAM within DMOF ⊃ PNIPAM composites can be tuned by changing the initial weight ratio between NIPAM, which is the monomer of PNIPAM, and DMOF. The guest PNIPAM chains in the composites block partially the 1D pores of DMOF, thus leading to a narrowed nanospace. The water adsorption studies reveal that the water uptake increased by increasing the loading of PNIPAM in the final DMOF ⊃ PNIPAM composites, indicating that the exposed amide groups of PNIPAM gradually alter the hydrophobicity of pristine DMOF and lead to hydrophilic DMOF ⊃ PNIPAM composites. The composite with the highest loading of PNIPAM displays a selective adsorption for water and methanol over ethanol when using equimolar mixtures of methanol-ethanol and water-ethanol. This is confirmed by the single-component adsorption measurements as well as ideal adsorbed solution theory molecular simulations. Additionally, the water stability of pristine DMOF has been greatly improved after the incorporation of PNIPAM in its pores. PNIPAM can undergo a phase transition between hydrophobic and hydrophilic phases in response to a low temperature change. This property is used in order to control the desorption of water and methanol molecules, thus enabling an efficient and cost-effective regeneration process.

Efficient Separation of Ethanol-Methanol and Ethanol-Water Mixtures Using ZIF-8 Supported on a Hierarchical Porous Mixed-Oxide Substrate

This work reports a new approach for the synthesis of a zeolitic imidazolate framework (ZIF-8) composite. It employs the direct growth of the crystalline ZIF-8 on a mixed-metal oxide support TiO₂-SiO₂ (TSO), which mimics the porous structure of Populus nigra. Using the natural leaf as a template, the TSO support was prepared using a sol-gel method. The growth of the ZIF-8 layer on the TSO support was carried out by the seeds and second growth method. This method facilitates the homogeneous dispersion of ZIF-8 crystals at the surface of the TSO composite. The ZIF-8@TSO composite adsorbs methanol selectively, mainly due to the hierarchical porous structure of the mixed oxide support. As compared with the as-synthesized ZIF-8, a 50% methanol uptake is achieved in the ZIF-8@TSO composite, with only 25 wt % ZIF-8 loading. IAST simulations show that the ZIF-8@TSO composite has a preferential adsorption toward methanol when using an equimolar methanol-ethanol mixture. An opposite behavior is observed for the as-synthesized ZIF-8. The results show that combining MOFs and mixed-oxide supports with bioinspired structures opens opportunities for synthesizing new materials with unique and enhanced adsorption and separation properties.

Polarizable Force Field for CO2 in M-MOF-74 Derived from Quantum Mechanics

On the short term, carbon capture is a viable solution to reduce human-induced CO₂ emissions, which requires an energy efficient separation of CO₂. Metal-organic frameworks (MOFs) may offer opportunities for carbon capture and other industrially relevant separations. Especially, MOFs with embedded open metal sites have been shown to be promising. Molecular simulation is a useful tool to predict the performance of MOFs even before the synthesis of the material. This reduces the experimental effort, and the selection process of the most suitable MOF for a particular application can be accelerated. To describe the interactions between open metal sites and guest molecules in molecular simulation is challenging. Polarizable force fields have potential to improve the description of such specific interactions. Previously, we tested the applicability of polarizable force fields for CO₂ in M-MOF-74 by verifying the ability to reproduce experimental measurements. Here, we develop a predictive polarizable force field for CO₂ in M-MOF-74 (M = Co, Fe, Mg, Mn, Ni, Zn) without the requirement of experimental data. The force field is derived from energies predicted from quantum mechanics. The procedure is easily transferable to other MOFs. To incorporate explicit polarization, the induced dipole method is applied between the framework and the guest molecule. Atomic polarizabilities are assigned according to the literature. Only the Lennard-Jones parameters of the open metal sites are parameterized to reproduce energies from quantum mechanics. The created polarizable force field for CO₂ in M-MOF-74 can describe the adsorption well and even better than that in our previous work.

Mechanical Properties in Metal-Organic Frameworks: Emerging Opportunities and Challenges for Device Functionality and Technological Applications

Some of the most remarkable recent developments in metal-organic framework (MOF) performance properties can only be rationalized by the mechanical properties endowed by their hybrid inorganic-organic nanoporous structures. While these characteristics create intriguing application prospects, the same attributes also present challenges that will need to be overcome to enable the integration of MOFs with technologies where these promising traits can be exploited. In this review, emerging opportunities and challenges are identified for MOF-enabled device functionality and technological applications that arise from their fascinating mechanical properties. This is discussed not only in the context of their more well-studied gas storage and separation applications, but also for instances where MOFs serve as components of functional nanodevices. Recent advances in understanding MOF mechanical structure-property relationships due to attributes such as defects and interpenetration are highlighted, and open questions related to state-of-the-art computational approaches for quantifying their mechanical properties are critically discussed.

On flexible force fields for metal-organic frameworks: Recent developments and future prospects

Classical force field simulations can be used to study structural, diffusion, and adsorption properties of metal-organic frameworks (MOFs). To account for the dynamic behavior of the material, parameterization schemes have been developed to derive force constants and the associated reference values by fitting on ab initio energies, vibrational frequencies, and elastic constants. Here, we review recent developments in flexible force field models for MOFs. Existing flexible force field models are generally able to reproduce the majority of experimentally observed structural and dynamic properties of MOFs. The lack of efficient sampling schemes for capturing stimuli-driven phase transitions, however, currently limits the full predictive potential of existing flexible force fields from being realized. This article is categorized under: Structure and Mechanism > Computational Materials ScienceMolecular and Statistical Mechanics > Molecular Mechanics.

Elucidating the Variable-Temperature Mechanical Properties of a Negative Thermal Expansion Metal-Organic Framework

We report the first experimental study into the thermomechanical and viscoelastic properties of a metal-organic framework (MOF) material. Nanoindentations show a decrease in the Young's modulus, consistent with classical molecular dynamics simulations, and hardness of HKUST-1 with increasing temperature over the 25-100 °C range. Variable-temperature dynamic mechanical analysis reveals significant creep behavior, with a reduction of 56% and 88% of the hardness over 10 min at 25 and 100 °C, respectively. This result suggests that, despite the increased density that results from increasing temperature in the negative thermal expansion MOF, the thermally induced softening due to vibrational and entropic contributions plays a more dominant role in dictating the material's temperature-dependent mechanical behavior.

Absorption Refrigeration Cycles with Ammonia-Ionic Liquid Working Pairs Studied by Molecular Simulation

For absorption refrigeration, it has been shown that ionic liquids have the potential to replace conventional working pairs. Due to the huge number of possibilities, conducting lab experiments to find the optimal ionic liquid is infeasible. Here, we provide a proof-of-principle study of an alternative computational approach. The required thermodynamic properties, i.e., solubility, heat capacity, and heat of absorption, are determined via molecular simulations. These properties are used in a model of the absorption refrigeration cycle to estimate the circulation ratio and the coefficient of performance. We selected two ionic liquids as absorbents: [emim][Tf₂N], and [emim][SCN]. As refrigerant NH₃ was chosen due to its favorable operating range. The results are compared to the traditional approach in which parameters of a thermodynamic model are fitted to reproduce experimental data. The work shows that simulations can be used to predict the required thermodynamic properties to estimate the performance of absorption refrigeration cycles. However, high-quality force fields are required to accurately predict the cycle performance.

Adsorption equilibrium of nitrogen dioxide in porous materials

The effect of confinement on the equilibrium reactive system containing nitrogen dioxide and dinitrogen tetroxide is studied by molecular simulation and the reactive Monte Carlo (RxMC) approach. The bulk-phase reaction was successfully reproduced and five all-silica zeolites (i.e. FAU, FER, MFI, MOR, and TON) with different topologies were selected to study their adoption behavior. Dinitrogen tetroxide showed a stronger affinity than nitrogen dioxide in all the zeolites due to size effects, but exclusive adsorption sites in MOR allowed the adsorption of nitrogen dioxide with no competition at these sites. From the study of the adsorption isotherms and isobars of the reacting mixture, confinement enhanced the formation of dimers over the full range of pressure and temperature, finding the largest deviations from bulk fractions at low temperature and high pressure. The channel size and shape of the zeolite have a noticeable influence on the dinitrogen tetroxide formation, being more important in MFI, closely followed by TON and MOR, and finally FER and FAU. Preferential adsorption sites in MOR lead to an unusually strong selective adsorption towards nitrogen dioxide, demonstrating that the topological structure has a crucial influence on the composition of the mixture and must be carefully considered in systems containing nitrogen dioxide.

Selective CO2 adsorption in water-stable alkaline-earth based metal–organic frameworks

Four novel metal–organic frameworks (MOFs) built from alkaline-earth metal ions and the flexible ligand tetrakis[4-(carboxyphenyl)oxamethyl]methane acid (H₄L) were synthesized. Among them, Ca-based MOFs have high stability under humidity conditions and show good CO₂ uptake.

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

Including explicit polarization significantly improves the description of the adsorption in comparison to non-polarizable generic force fields.

Efficient Application of Continuous Fractional Component Monte Carlo in the Reaction Ensemble

A new formulation of the Reaction Ensemble Monte Carlo technique (RxMC) combined with the Continuous Fractional Component Monte Carlo method is presented. This method is denoted by serial Rx/CFC. The key ingredient is that fractional molecules of either reactants or reaction products are present and that chemical reactions always involve fractional molecules. Serial Rx/CFC has the following advantages compared to other approaches: (1) One directly obtains chemical potentials of all reactants and reaction products. Obtained chemical potentials can be used directly as an independent check to ensure that chemical equilibrium is achieved. (2) Independent biasing is applied to the fractional molecules of reactants and reaction products. Therefore, the efficiency of the algorithm is significantly increased, compared to the other approaches. (3) Changes in the maximum scaling parameter of intermolecular interactions can be chosen differently for reactants and reaction products. (4) The number of fractional molecules is reduced. As a proof of principle, our method is tested for Lennard-Jones systems at various pressures and for various chemical reactions. Excellent agreement was found both for average densities and equilibrium mixture compositions computed using serial Rx/CFC, RxMC/CFCMC previously introduced by Rosch and Maginn (Journal of Chemical Theory and Computation, 2011, 7, 269-279), and the conventional RxMC approach. The serial Rx/CFC approach is also tested for the reaction of ammonia synthesis at various temperatures and pressures. Excellent agreement was found between results obtained from serial Rx/CFC, experimental results from literature, and thermodynamic modeling using the Peng-Robinson equation of state. The efficiency of reaction trial moves is improved by a factor of 2 to 3 (depending on the system) compared to the RxMC/CFCMC formulation by Rosch and Maginn.

Flexible Force Field Parameterization through Fitting on the Ab Initio-Derived Elastic Tensor

Constructing functional forms and their corresponding force field parameters for the metal–linker interface of metal–organic frameworks is challenging. We propose fitting these parameters on the elastic tensor, computed from ab initio density functional theory calculations. The advantage of this top-down approach is that it becomes evident if functional forms are missing when components of the elastic tensor are off. As a proof-of-concept, a new flexible force field for MIL-47(V) is derived. Negative thermal expansion is observed and framework flexibility has a negligible effect on adsorption and transport properties for small guest molecules. We believe that this force field parametrization approach can serve as a useful tool for developing accurate flexible force field models that capture the correct mechanical behavior of the full periodic structure.

Behavior of the Enthalpy of Adsorption in Nanoporous Materials Close to Saturation Conditions

Many important industrial separation processes based on adsorption operate close to saturation. In this regime, the underlying adsorption processes are mostly driven by entropic forces. At equilibrium, the entropy of adsorption is closely related to the enthalpy of adsorption. Thus, studying the behavior of the enthalpy of adsorption as a function of loading is fundamental to understanding separation processes. Unfortunately, close to saturation, the enthalpy of adsorption is hard to measure experimentally and hard to compute in simulations. In simulations, the enthalpy of adsorption is usually obtained from energy/particle fluctuations in the grand-canonical ensemble, but this methodology is hampered by vanishing insertions/deletions at high loading. To investigate the fundamental behavior of the enthalpy and entropy of adsorption at high loading, we develop a simplistic model of adsorption in a channel and show that at saturation the enthalpy of adsorption diverges to large positive values due to repulsive intermolecular interactions. However, there are many systems that can avoid repulsive intermolecular interactions and hence do not show this drastic increase in enthalpy of adsorption close to saturation. We find that the conventional grand-canonical Monte Carlo method is incapable of determining the enthalpy of adsorption from energy/particle fluctuations at high loading. Here, we show that by using the continuous fractional component Monte Carlo, the enthalpy of adsorption close to saturation conditions can be reliably obtained from the energy/particle fluctuations in the grand-canonical ensemble. The best method to study properties at saturation is the NVT energy (local-) slope methodology.

Polarizable Force Fields for CO2 and CH4 Adsorption in M-MOF-74

The family of M-MOF-74, with M = Co, Cr, Cu, Fe, Mg, Mn, Ni, Ti, V, and Zn, provides opportunities for numerous energy related gas separation applications. The pore structure of M-MOF-74 exhibits a high internal surface area and an exceptionally large adsorption capacity. The chemical environment of the adsorbate molecule in M-MOF-74 can be tuned by exchanging the metal ion incorporated in the structure. To optimize materials for a given separation process, insights into how the choice of the metal ion affects the interaction strength with adsorbate molecules and how to model these interactions are essential. Here, we quantitatively highlight the importance of polarization by comparing the proposed polarizable force field to orbital interaction energies from DFT calculations. Adsorption isotherms and heats of adsorption are computed for CO2, CH4, and their mixtures in M-MOF-74 with all 10 metal ions. The results are compared to experimental data, and to previous simulation results using nonpolarizable force fields derived from quantum mechanics. To the best of our knowledge, the developed polarizable force field is the only one so far trying to cover such a large set of possible metal ions. For the majority of metal ions, our simulations are in good agreement with experiments, demonstrating the effectiveness of our polarizable potential and the transferability of the adopted approach.

Controlling Thermal Expansion: A Metal-Organic Frameworks Route

Controlling thermal expansion is an important, not yet resolved, and challenging problem in materials research. A conceptual design is introduced here, for the first time, for the use of metal-organic frameworks (MOFs) as platforms for controlling thermal expansion devices that can operate in the negative, zero, and positive expansion regimes. A detailed computer simulation study, based on molecular dynamics, is presented to support the targeted application. MOF-5 has been selected as model material, along with three molecules of similar size and known differences in terms of the nature of host-guest interactions. It has been shown that adsorbate molecules can control, in a colligative way, the thermal expansion of the solid, so that changing the adsorbate molecules induces the solid to display positive, zero, or negative thermal expansion. We analyze in depth the distortion mechanisms, beyond the ligand metal junction, to cover the ligand distortions, and the energetic and entropic effect on the thermo-structural behavior. We provide an unprecedented atomistic insight on the effect of adsorbates on the thermal expansion of MOFs as a basic tool toward controlling the thermal expansion.

A Simulation Study of Alkanes in Linde Type A Zeolites

Monte Carlo simulations were performed to study the adsorption and diffusion of small hydrocarbons in Linde Type A zeolites as a function of their calcium/sodium ratio. The diffusion studies were focused on methane whereas the adsorption simulations were performed from methane up to pentane. The results obtained showed that an increase in the number of cations in the structure (exchange of univalent sodium ions by divalent calcium ions) led to an increase in the adsorption of linear alkanes at low and medium pressure, but caused a decrease in adsorption at the highest pressures. An increase in the amount of cations favours molecular attraction and hence results in lower mobility. At higher cation loading the ions block the windows interconnecting the LTA cages, leading to a further decrease in diffusion. Methane self-diffusion coefficients obtained from our simulations were twice as high for the Linde Type 5A zeolite as for the Linde Type 4A zeolite. These results are consistent with previous experimental studies and provide a molecular picture of the influence of the zeolite type, the amount of cations contained and their location in the structure.

Understanding and solving disorder in the substitution pattern of amino functionalized MIL-47(V)

Electronic energies and elastic constants of four amino functionalized MIL-47(V) supercells were computed using the plane wave density functional theory to determine the influence of the substituent positions on the organic linker.

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.

Computing the Heat of Adsorption using Molecular Simulations: The Effect of Strong Coulombic Interactions

Molecular simulations are an important tool for the study of adsorption of hydrocarbons in nanoporous materials such as zeolites. The heat of adsorption is an important thermodynamic quantity that can be measured both in experiments and molecular simulations, and therefore it is often used to investigate the quality of a force field for a certain guest-host (g - h) system. In molecular simulations, the heat of adsorption in zeolites is often computed using either of the following methods: (1) using the Clausius-Clapeyron equation, which requires the partial derivative of the pressure with respect to temperature at constant loading, (2) using the energy difference between the host with and without a single guest molecule present, and (3) from energy/particle fluctuations in the grand-canonical ensemble. To calculate the heat of adsorption from experiments (besides direct calorimetry), only the first method is usually applicable. Although the computation of the heat of adsorption is straightforward for all-silica zeolites, severe difficulties arise when applying the conventional methods to systems with nonframework cations present. The reason for this is that these nonframework cations have very strong Coulombic interactions with the zeolite. We will present an alternative method based on biased interactions of guest molecules that suffers less from these difficulties. This method requires only a single simulation of the host structure. In addition, we will review some of the other important issues concerning the handling of these strong Coulombic interactions in simulating the adsorption of guest molecules. It turns out that the recently proposed Wolf method ( J. Chem. Phys. 1999, 110 , 8254 ) performs poorly for zeolites as a large cutoff radius is needed for convergence.

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.