Repository logo
Log in(current)
Repository logoMIT Open ScholarshipDSpace@MIT
  1. Home
  2. MIT Open Access Articles
  3. MIT Open Access Articles
  4. Benchmarking binding energy calculations for organic structure-directing agents in pure-silica zeolites

Benchmarking binding energy calculations for organic structure-directing agents in pure-silica zeolites

Thumbnail Image
Download
Name

5.0044927.pdf

Description
Published version
Size

2.7 MB

Format

Adobe PDF

Checksum (MD5)

b4832e0f97eac09f8e468f217b8c8f72

sword-2022-05-13T15:34:14.original.xml (130 B)
Original SWORD entry document
Author(s)
Schwalbe-Koda, Daniel
•
Gómez-Bombarelli, Rafael
Date Issued
2021
Journal
The Journal of Chemical Physics
Publisher
AIP Publishing
Citation
Schwalbe-Koda, Daniel and Gómez-Bombarelli, Rafael. 2021. "Benchmarking binding energy calculations for organic structure-directing agents in pure-silica zeolites." The Journal of Chemical Physics, 154 (17).
Version
Final published version
Abstract
Molecular modeling plays an important role in the discovery of organic structure-directing agents (OSDAs) for zeolites. By quantifying the intensity of host-guest interactions, it is possible to select cost-effective molecules that maximize binding toward a given zeolite framework. Over the last few decades, a variety of methods and levels of theory have been used to calculate these binding energies. Nevertheless, there is no consensus on the best calculation strategy for high-throughput virtual screening undertakings. In this work, we compare binding affinities from density functional theory (DFT) and Dreiding force field calculations for 272 zeolite-OSDA pairs obtained from static and time-averaged simulations. Enabled by automation software, we show that Dreiding binding energies from the frozen pose method correlate best with DFT energies. They are also less sensitive to the choice of initial lattice parameters and optimization algorithms, as well as less computationally expensive than their time-averaged counterparts. Furthermore, we demonstrate that a broader exploration of the conformation space from molecular dynamics simulations does not provide significant improvements in binding energy trends over the frozen pose method despite being orders of magnitude more expensive. The code and benchmark data are open-sourced and provide robust and computationally efficient guidelines to calculating binding energies in zeolite-OSDA pairs.
MIT Department
Massachusetts Institute of Technology. Department of Materials Science and Engineering
Terms of Use
Creative Commons Attribution 4.0 International License
https://creativecommons.org/licenses/by/4.0
Persistent DSpace Link
https://hdl.handle.net/1721.1/142528
DOI of Published Version
https://doi.org/10.1063/5.0044927
Repository logo
PrivacyPermissionsAccessibilityContact us
Repository logo
Notify us about copyright concerns.