Tailoring Parameters for QM/MM Simulations: Accurate Modeling of Adsorption and Catalysis in Zirconium-Based Metal-Organic Frameworks

Yu-Chi Kao,Yi-Ming Wang,Jyun-Yi Yeh,Shih-Cheng Li,Kevin C.-W. Wu,Li-Chiang Lin,Yi-Pei Li
DOI: https://doi.org/10.26434/chemrxiv-2024-8md2k
2024-03-05
Abstract:Quantum mechanics/molecular mechanics (QM/MM) simulations offer an efficient way to model reactions occurring in complex environments. This study introduces a specialized set of charge and Lennard-Jones parameters tailored for electrostatically embedded QM/MM calculations, aiming to accurately model both adsorption processes and catalytic reactions in zirconium-based metal-organic frameworks (Zr-MOFs). To validate our approach, we compare adsorption energies derived from QM/MM simulations against experimental results and Monte Carlo simulation outcomes. The developed parameters showcase the ability of QM/MM simulations to represent long-range electrostatic and van der Waals interactions faithfully. This capability is evidenced by the prediction of adsorption energies with a low root mean square error of 1.1 kcal/mol across a wide range of adsorbates. The practical applicability of our QM/MM model is further illustrated through the study of glucose isomerization and epimerization reactions catalyzed by two structurally distinct Zr-MOF catalysts, UiO-66 and MOF-808. Our QM/MM calculations closely align with experimental activation energies. Importantly, the parameter set introduced here is shown to be compatible with the widely used universal force field (UFF). Moreover, we thoroughly explore how the size of the cluster model and the choice of density functional theory (DFT) methodologies influence the simulation outcomes. This work provides an accurate and computationally efficient framework for modeling complex catalytic reactions within Zr-MOFs, contributing valuable insights into their mechanistic behaviors and facilitating further advancements in this dynamic area of research.
Chemistry
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a set of charge and Lennard - Jones parameters specifically for zirconium - based metal - organic frameworks (Zr - MOFs) in order to improve the accuracy of quantum mechanics/molecular mechanics (QM/MM) simulations in adsorption processes and catalytic reactions. Specifically, the objectives of the paper include: 1. **Parameter Optimization**: By adjusting the charge and Lennard - Jones parameters, enable QM/MM simulations to more accurately describe the long - range electrostatic interactions and van der Waals forces in Zr - MOFs. This includes verifying the effectiveness of these parameters through experimental data and Monte Carlo simulation results. 2. **Model Validation**: By calculating the adsorption energies of different adsorbates in Zr - MOFs and comparing them with experimental results, verify the accuracy of the developed parameter set. The paper shows that these parameters have a low root - mean - square error (RMSE) of 1.1 kcal/mol and a mean absolute error (MAE) of 1.0 kcal/mol in predicting adsorption energies. 3. **Application Example**: By studying the catalytic mechanisms of glucose isomerization and epimerization reactions on two structurally different Zr - MOF catalysts (UiO - 66 and MOF - 808), demonstrate the practical application value of the developed QM/MM model. These calculation results are highly consistent with the experimental activation energies, further proving the effectiveness of the model. 4. **Method Compatibility**: Explore the influence of different density functional theory (DFT) methods on the simulation results, and ensure that the developed parameter set is compatible with multiple commonly used functional functions (such as ωB97X - D, ωB97X - V, B97 - D, and M062X), thereby providing an accurate and computationally efficient framework. In summary, this paper aims to provide a reliable method for modeling complex catalytic reactions in Zr - MOFs by optimizing and validating the parameters used in QM/MM simulations, thereby promoting further development in this field.