Understanding the kinetic and thermodynamic origins of xylene separation in UiO-66(Zr) via molecular simulation

Xylene isomers are precursors in many important chemical processes, yet their separation via crystallization or distillation is energy intensive. Adsorption presents an attractive, lower-energy alternative and the discovery of adsorbents which outperform the current state-of-the-art zeolitic materials represents one of the key challenges in materials design, with metal-organic frameworks receiving particular attention. One of the most well-studied systems in this context is UiO-66(Zr), which selectively adsorbs ortho-xylene over the other C8 alkylaromatics. The mechanism behind this separation has remained unclear, however. In this work, we employ a wide range of computational techniques to explore both the equilibrium and dynamic behavior of the xylene isomers in UiO-66(Zr). In addition to correctly predicting the experimentally-observed ortho-selectivity, we demonstrate that the equilibrium selectivity is based upon the complete encapsulation of ortho-xylene within the pores of the framework. Furthermore the flexible nature of the adsorbent is crucial in facilitating xylene diffusion and our simulations reveal for the first time significant differences between the intracrystalline diffusion mechanisms of the three isomers resulting in a kinetic contribution to the selectivity. Consequently it is important to include both equilibrium and kinetic effects when screening MOFs for xylene separations

Developing and validating Fuzzy-Border continuum solvation model with POlarizable Simulations Second order Interaction Model (POSSIM) force field for proteins

The accurate, fast and low cost computational tools are indispensable for studying the structure and dynamics of biological macromolecules in aqueous solution. The goal of this thesis is development and validation of continuum Fuzzy-Border (FB) solvation model to work with the Polarizable Simulations Second-order Interaction Model (POSSIM) force field for proteins developed by Professor G A Kaminski. The implicit FB model has advantages over the popularly used Poisson Boltzmann (PB) solvation model. The FB continuum model attenuates the noise and convergence issues commonly present in numerical treatments of the PB model by employing fixed position cubic grid to compute interactions. It also uses either second or first-order approximation for the solvent polarization which is similar to the second-order explicit polarization applied in POSSIM force field. The FB model was first developed and parameterized with nonpolarizable OPLS-AA force field for small molecules which are not only important in themselves but also building blocks of proteins and peptide side chains. The hydration parameters are fitted to reproduce the experimental or quantum mechanical hydration energies of the molecules with the overall average unsigned error of ca. 0.076kcal/mol. It was further validated by computing the absolute pKa values of 11 substituted phenols with the average unsigned error of 0.41pH units in comparison with the quantum mechanical error of 0.38pH units for this set of molecules. There was a good transferability of hydration parameters and the results were produced only with fitting of the specific atoms to the hydration energy and pKa targets. This clearly demonstrates the numerical and physical basis of the model is good enough and with proper fitting can reproduce the acidity constants for other systems as well. After the successful development of FB model with the fixed charges OPLS-AA force field, it was expanded to permit simulations with Polarizable Simulations Second-order Interaction Model (POSSIM) force field. The hydration parameters of the small molecules representing analogues of protein side chains were fitted to their solvation energies at 298.15K with an average error of ca.0.136kcal/mol. Second, the resulting parameters were used to reproduce the pKa values of the reference systems and the carboxylic (Asp7, Glu10, Glu19, Asp27 and Glu43) and basic residues (Lys13, Lys29, Lys34, His52 and Lys55) of the turkey ovomucoid third domain (OMTKY3) protein. The overall average unsigned error in the pKa values of the acid residues was found to be 0.37pH units and the basic residues was 0.38 pH units compared to 0.58pH units and 0.72 pH units calculated previously using polarizable force field (PFF) and Poisson Boltzmann formalism (PBF) continuum solvation model. These results are produced with fitting of specific atoms of the reference systems and carboxylic and basic residues of the OMTKY3 protein. Since FB model has produced improved pKa shifts of carboxylic residues and basic protein residues in OMTKY3 protein compared to PBF/PFF, it suggests the methodology of first-order FB continuum solvation model works well in such calculations. In this study the importance of explicit treatment of the electrostatic polarization in calculating pKa of both acid and basic protein residues is also emphasized. Moreover, the presented results demonstrate not only the consistently good degree of accuracy of protein pKa calculations with the second-degree POSSIM approximation of the polarizable calculations and the first-order approximation used in the Fuzzy-Border model for the continuum solvation energy, but also a high degree of transferability of both the POSSIM and continuum solvent Fuzzy Border parameters. Therefore, the FB model of solvation combined with the POSSIM force field can be successfully applied to study the protein and protein-ligand systems in water

The right isotherms for the right reasons?: validation of generic force fields for prediction of methane adsorption in metal-organic frameworks

In recent years, the use of computational tools to aid in the evaluation, understanding and design of advanced porous materials for gas storage and separation processes has become ever-more widespread. High-performance computing facilities have become more powerful and more accessible and molecular simulation of gas adsorption has become routine, often involving the use of a number of default and commonly-used parameters as a result. In this work, we consider the application of molecular simulation in one particular field of adsorption – the prediction of methane adsorption in metal-organic frameworks in the low-loading regime – and employ a range of computational techniques to evaluate the appropriateness of many commonly chosen simulation parameters to these systems. In addition to confirming the power of relatively simple generic force fields to quickly and accurately predict methane adsorption isotherms in a range of MOFs, we demonstrate that these force fields are capable of providing detailed molecular-level information which is in very good agreement with quantum chemical predictions. We highlight a number of chemical systems in which molecular-level insight from generic force fields should be approached with a degree of caution and provide some general recommendations for best-practice in simulations of CH4 adsorption in MOFs

Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly

Quantum Mechanical effects in n-alkane droplets

According to quantum chemical (QC) calculations of a series of n-alkanes (CnH2n+2 at n = 1-8) the gap ΔHL between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) decreases with increasing methylene group (-CH2-) number in the n-alkane chain. While the LUMO level appears to be relatively unchanged, in good agreement with experimental results, the HOMO level is unstabilized and ΔHL decreases from approximately 11 eV (methane) to 6.5 eV (n-octane). Comparative ab initio calculations show not only the quantum confinement effects (QCE) in n-alkanes in the gas phase due to the reduction of the ΔHL value, but also enabled selection of cost-effective methodologies for modelling long-chain n-alkanes (n > 8). The used methodologies include ‘Local Spin Density Approximation’, combining special exchange functional with suitable correlation functional. Electronic structures and energetics of n-pentane (C5H12) in the evaporation/condensation processes are studied to understand the molecular mechanism for these processes. Two main step processes, liquid (L) 1 liquid-gas (L-G) interface 2 gas (G), are analysed. While phase transformations between L and L-G is estimated as activation processes, it is found that there is no energy barrier in the transitions between L-G and

Accuracy of quantum mechanically derived force-fields parameterized from dispersion-corrected DFT data: the benzene dimer as a prototype for aromatic interactions

A multilevel approach is presented to assess the ability of several popular dispersion corrected density functionals (M06-2X, CAM-B3LYP-D3, BLYP-D3, and B3LYP-D3) to reliably describe two-body interaction potential energy surfaces (IPESs). To this end, the automated Picky procedure (Cacelli et al. J. Comput. Chem. 2012, 33, 1055) was exploited, which consists in parametrizing specific intermolecular force fields through an iterative approach, based on the comparison with quantum mechanical data. For each of the tested functionals, the resulting force field was employed in classical Monte Carlo and Molecular Dynamics simulations, performed on systems of up to 1000 molecules in ambient conditions, to calculate a number of condensed phase properties. The comparison of the resulting structural and dynamic properties with experimental data allows us to assess the quality of each IPES and, consequently, even the quality of the DFT functionals. The methodology is tested on the benzene dimer, commonly used as a benchmark molecule, a prototype of aromatic interactions. The best results were obtained with the CAM-B3LYP-D3 functional. Besides assessing the reliability of DFT functionals in describing aromatic IPESs, this work provides a further step toward a robust protocol for the derivation of sound force field parameters from quantum mechanical data. This method can be relevant in all those cases where standard force fields fail in giving accurate predictions

A molecular interpretation of the dynamics of diffusive mass transport of water within a glassy polyetherimide

The diffusion process of water molecules within a polyetherimide (PEI) glassy matrix has been analyzed by combining the experimental analysis of water sorption kinetics performed by FTIR spectroscopy with theoretical information gathered from Molecular Dynamics simulations and with the expression of water chemical potential provided by a non‐equilibrium lattice fluid model able to describe the thermodynamics of glassy polymers. This approach allowed us to construct a convincing description of the diffusion mechanism of water in PEI providing molecular details of the process related to the effects of the cross‐ and self‐hydrogen bonding established in the system on the dynamics of water mass transport