Rice University logo

Ab initio free energy calculations for adsorption and reactions in nanoporous systems



By: Joachim Sauer
Professor of Chemistry
From: Humboldt University Berlin
When: Wednesday, April 26, 2017
4:00 PM - 5:00 PM
Where: Dell Butcher Hall
Abstract: The ab initio prediction of reaction rate or equilibrium constants for systems with hundreds of atoms with an accuracy that is comparable to experiment is a challenge for computational quantum chemistry. We present a divide-and-conquer strategy that departs from a potential energy surface obtained by standard density functional theory with inclusion of dispersion. The energies of the reactant and transition structures are refined by wave-function-type electron correlation calculations for the reaction site.[1] Thermal effects and entropies are calculated from vibrational partition functions. Anharmonic frequencies are calculated for each vibrational mode separately.[2] For a key reaction of an industrially relevant catalytic process, the methylation of small alkenes over zeolites, we obtain results that agree within chemical accuracy limits with experiment. Deviations are within one order of magnitude for rate constants (free energies), and pre-exponential factors (entropies), and within 4 kJ/mol for enthalpy barriers.[3] This methodology also yields chemically accurate (± 4 kJ/mol or less) free energies of adsorption (Henry constants) for small molecules on Brønsted sites in zeolites,[4] or on metal ion and linker sites on the internal surfaces of metal organic frameworks (MOF).[5] The accurate ab initio prediction of adsorption isotherms and selectivities with no other input than the atomic positions is prerequisite to a rational design of materials for gas storage, e.g. for energy carriers such as H2, and separation, e.g. removal of CO2 from CH4. This we achieve with Grand Canonical Monte Carlo (GCMC) simulations on a lattice of adsorption sites. The Hamiltonian is defined by Gibbs free energies of adsorption on individual sites and lateral interaction energies (adsorbate-adsorbate) calculated ab initio. [1] C. Tuma, J. Sauer, Phys. Chem. Chem. Phys. 2006, 8, 3955-3965. [2] G. Piccini, J. Sauer, J. Chem. Theory Comput. 2014, 10, 2479-2487. [3] G. Piccini, M. Alessio, J. Sauer, Angew. Chem., Int. Ed. 2016, 55, 5235-5237. [4] G. Piccini, M. Alessio, J. Sauer, Y. Zhi, Y. Liu, R. Kolvenbach, A. Jentys, J. A. Lercher, J. Phys. Chem. C 2015, 119, 6128-6137. [5] A. Kundu, G. Piccini, K. Sillar, J. Sauer, J. Am. Chem. Soc. 2016, 138, 14047-14056.