Tim M. Becker

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.

Investigating polarization effects of CO2 adsorption in MgMOF-74

Abstract MgMOF-74 is a promising candidate for a variety of gas separation applications, e.g., carbon capture and natural gas sweetening due to its high CO2 uptake capacity and its favorable selectivity toward CO2. Motivated by its promising properties, MgMOF-74 has been extensively studied both experimentally and computationally. Experimentally determined adsorption isotherms show an inflection at a loading of approximately one CO2 molecule per magnesium ion due to strong adsorption sites close to the ions. It is a great challenge to accurately reproduce this behavior in molecular simulations. In this study, we explicitly consider polarization between the adsorbed CO2 molecules and the framework of MgMOF-74 via the induced point dipole method. Back-polarization is neglected to achieve reasonable simulation times. To account for implicitly incorporated polarization, we rescale the Lennard–Jones energy parameters with respect to the atomic polarizabilities. A series of Monte Carlo simulations of CO2 in MgMOF-74 is conducted. The computed CO2 adsorption isotherm is in good agreement with experimental measurements and previous simulation results using a DFT-derived force field. This indicates that polarization is important for describing the adsorption of CO2 in MgMOF-74. The direct inclusion of polarization will lead to force fields with better physical justification and transferability.

Finite-size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard–Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell–Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer (J. Phys. Chem. B, 2004, 108, 15873−15879). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard–Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell–Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell–Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

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][Tf2N], and [emim][SCN]. As refrigerant NH3 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.

Gibbs ensemble Monte Carlo simulations of multicomponent natural gas mixtures

Abstract Vapour–liquid equilibrium (VLE) and volumetric data of multicomponent mixtures are extremely important for natural gas production and processing, but it is time consuming and challenging to experimentally obtain these properties. An alternative tool is provided by means of molecular simulation. Here, Monte Carlo (MC) simulations in the Gibbs ensemble are used to compute the VLE of multicomponent natural gas mixtures. Two multicomponent systems, one containing a mixture of six components (, , , S, and ), and the other containing a mixture of nine components (, , , S, , , , and ) are simulated. The computed VLE from the MC simulations is in good agreement with available experimental data and the GERG-2008 equation of state modelling. The results show that molecular simulation can be used to predict properties of multicomponent systems relevant for the natural gas industry. Guidelines are provided to setup Gibbs ensemble simulations for multicomponent systems, which is a challenging task due to the increased number of degrees of freedom.

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 CO2 emissions, which requires an energy efficient separation of CO2. 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 CO2 in M-MOF-74 by verifying the ability to reproduce experimental measurements. Here, we develop a predictive polarizable force field for CO2 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 CO2 in M-MOF-74 can describe the adsorption well and even better than that in our previous work.

Thermodynamic and Transport Properties of Crown-Ethers: Force Field Development and Molecular Simulations

Crown-ethers have recently been used to assemble porous liquids (PLs), which are liquids with permanent porosity formed by mixing bulky solvent molecules (e.g., 15-crown-5 ether) with solvent-inaccessible organic cages. PLs and crown-ethers belong to a novel class of materials, which can potentially be used for gas separation and storage, but their performance for this purpose needs to be assessed thoroughly. Here, we use molecular simulations to study the gas separation performance of crown-ethers as the solvent of porous liquids. The TraPPE force field for linear ether molecules has been adjusted by fitting a new set of torsional potentials to accurately describe cyclic crown-ether molecules. Molecular dynamics (MD) simulations have been used to compute densities, shear viscosities, and self-diffusion coefficients of 12-crown-4, 15-crown-5, and 18-crown-6 ethers. In addition, Monte Carlo (MC) simulations have been used to compute the solubility of the gases CO2, CH4, and N2 in 12-crown-4 and 15-crown-5 ether. The computed properties are compared with available experimental data of crown-ethers and their linear counterparts, i.e., polyethylene glycol dimethyl ethers.