4.6 Article

Improved Force-Field Parameters for QM/MM Simulations of the Energies of Adsorption for Molecules in Zeolites and a Free Rotor Correction to the Rigid Rotor Harmonic Oscillator Model for Adsorption Enthalpies

Journal

JOURNAL OF PHYSICAL CHEMISTRY C
Volume 119, Issue 4, Pages 1840-1850

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/jp509921r

Keywords

-

Funding

  1. XC2 program - BP
  2. National Science Foundation [CHE-1265731]
  3. Direct For Mathematical & Physical Scien
  4. Division Of Chemistry [1265731] Funding Source: National Science Foundation

Ask authors/readers for more resources

Quantum mechanics/molecular mechanics (QM/MM) simulations provide an efficient avenue for studying reactions catalyzed in zeolite systems; however, the accuracy of such calculations is highly dependent on the zeolite MM parameters used. Previously reported parameters (P1), which were chosen to minimize the root mean square (RMS) deviations of adsorption energies compared with full QM omega B97X-D/6-31+G** adsorption energies, are shown to overestimate binding energies compared with experimental values, particularly for larger substrates. To address this issue, a new parameter set (P2) is derived by rescaling the previously reported characteristic energies of the Lennard-Jones potential in P1. The accuracy of the thermal correction for adsorption enthalpies determined by the rigid rotor-harmonic oscillator approximation (RRHO) is examined and shown to be improved by treating low-lying vibrational modes as free translational and rotational modes via a quasi-RRHO model. With P2 and quasi-RRHO, adsorption energies calculated with QM/MM agree with experimental values with an RMS error of 1.8 kcal/mol for both nonpolar and polar molecules adsorbed in MFI, H-MFI, and H-BEA. By contrast, the RMS error for the same test sets obtained using parameter set P1 is 8.3 kcal/mol. Glucose-fructose isomerization catalyzed by Sn-BEA is taken as an example to demonstrate that improved values for apparent activation energies can be obtained using the methodology reported here. With parameter set P2, the apparent activation energy calculated with QM/MM reproduces the experimental value to within 1 kcal/mol. By contrast, using parameter set P1, the error is -12.9 kcal/mol.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.6
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available