Combined quantum mechanical and molecular mechanical method for metal-organic frameworks

Proton topologies of NU-1000

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Metal-organic frameworks (MOFs) are materials with applications in catalysis, gas separations, and storage. Quantum mechanical (QM) calculations can provide valuable guidance to understand and predict their properties. In order to make the calculations faster, rather than modeling these materials as periodic (infinite) systems, it is useful to construct finite models (called cluster models) and use subsystem methods such as fragment methods or combined quantum mechanical and molecular mechanical (QM/MM) methods. Here we employ a QM/MM methodology to study one particular MOF that has been of widespread interest because of its wide pores and good solvent and thermal stability, namely NU-1000, which contains hexanuclear zirconium nodes and 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy4-) linkers. A modified version of the Bristow-Tiana-Walsh transferable force field has been developed to allow QM/MM calculations on NU-1000; we call the new parametrization the NU1T force field. We consider isomeric structures corresponding to various proton topologies of the [Zr63-O)8O8H16]8+ node of NU-1000, and we compute their relative energies using a QM/MM scheme designed for the present kind of problem. We compared the results to full quantum mechanical (QM) energy calculations and found that the QM/MM models can reproduce the full QM relative energetics (which span a range of 334 kJ mol-1) with a mean unsigned deviation (MUD) of only 2 kJ mol-1. Furthermore, we found that the structures optimized by QM/MM are nearly identical to their full QM optimized counterparts.

Original languageEnglish (US)
Pages (from-to)1778-1786
Number of pages9
JournalPhysical Chemistry Chemical Physics
Volume20
Issue number3
DOIs
StatePublished - Jan 1 2018

Fingerprint

Protons
topology
Metals
Topology
protons
metals
field theory (physics)
Benzoic Acid
benzoic acid
pyrenes
Catalysis
catalysis
Thermodynamic stability
thermal stability
Gases
fragments
methodology
porosity
deviation
energy

PubMed: MeSH publication types

  • Journal Article

Cite this

@article{efb5943792cd41f4b3ab4088e78d1793,
title = "Combined quantum mechanical and molecular mechanical method for metal-organic frameworks: Proton topologies of NU-1000",
abstract = "Metal-organic frameworks (MOFs) are materials with applications in catalysis, gas separations, and storage. Quantum mechanical (QM) calculations can provide valuable guidance to understand and predict their properties. In order to make the calculations faster, rather than modeling these materials as periodic (infinite) systems, it is useful to construct finite models (called cluster models) and use subsystem methods such as fragment methods or combined quantum mechanical and molecular mechanical (QM/MM) methods. Here we employ a QM/MM methodology to study one particular MOF that has been of widespread interest because of its wide pores and good solvent and thermal stability, namely NU-1000, which contains hexanuclear zirconium nodes and 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy4-) linkers. A modified version of the Bristow-Tiana-Walsh transferable force field has been developed to allow QM/MM calculations on NU-1000; we call the new parametrization the NU1T force field. We consider isomeric structures corresponding to various proton topologies of the [Zr6(μ3-O)8O8H16]8+ node of NU-1000, and we compute their relative energies using a QM/MM scheme designed for the present kind of problem. We compared the results to full quantum mechanical (QM) energy calculations and found that the QM/MM models can reproduce the full QM relative energetics (which span a range of 334 kJ mol-1) with a mean unsigned deviation (MUD) of only 2 kJ mol-1. Furthermore, we found that the structures optimized by QM/MM are nearly identical to their full QM optimized counterparts.",
author = "Xinping Wu and Laura Gagliardi and Truhlar, {Donald G}",
year = "2018",
month = "1",
day = "1",
doi = "10.1039/c7cp06751h",
language = "English (US)",
volume = "20",
pages = "1778--1786",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "3",

}

TY - JOUR

T1 - Combined quantum mechanical and molecular mechanical method for metal-organic frameworks

T2 - Proton topologies of NU-1000

AU - Wu, Xinping

AU - Gagliardi, Laura

AU - Truhlar, Donald G

PY - 2018/1/1

Y1 - 2018/1/1

N2 - Metal-organic frameworks (MOFs) are materials with applications in catalysis, gas separations, and storage. Quantum mechanical (QM) calculations can provide valuable guidance to understand and predict their properties. In order to make the calculations faster, rather than modeling these materials as periodic (infinite) systems, it is useful to construct finite models (called cluster models) and use subsystem methods such as fragment methods or combined quantum mechanical and molecular mechanical (QM/MM) methods. Here we employ a QM/MM methodology to study one particular MOF that has been of widespread interest because of its wide pores and good solvent and thermal stability, namely NU-1000, which contains hexanuclear zirconium nodes and 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy4-) linkers. A modified version of the Bristow-Tiana-Walsh transferable force field has been developed to allow QM/MM calculations on NU-1000; we call the new parametrization the NU1T force field. We consider isomeric structures corresponding to various proton topologies of the [Zr6(μ3-O)8O8H16]8+ node of NU-1000, and we compute their relative energies using a QM/MM scheme designed for the present kind of problem. We compared the results to full quantum mechanical (QM) energy calculations and found that the QM/MM models can reproduce the full QM relative energetics (which span a range of 334 kJ mol-1) with a mean unsigned deviation (MUD) of only 2 kJ mol-1. Furthermore, we found that the structures optimized by QM/MM are nearly identical to their full QM optimized counterparts.

AB - Metal-organic frameworks (MOFs) are materials with applications in catalysis, gas separations, and storage. Quantum mechanical (QM) calculations can provide valuable guidance to understand and predict their properties. In order to make the calculations faster, rather than modeling these materials as periodic (infinite) systems, it is useful to construct finite models (called cluster models) and use subsystem methods such as fragment methods or combined quantum mechanical and molecular mechanical (QM/MM) methods. Here we employ a QM/MM methodology to study one particular MOF that has been of widespread interest because of its wide pores and good solvent and thermal stability, namely NU-1000, which contains hexanuclear zirconium nodes and 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy4-) linkers. A modified version of the Bristow-Tiana-Walsh transferable force field has been developed to allow QM/MM calculations on NU-1000; we call the new parametrization the NU1T force field. We consider isomeric structures corresponding to various proton topologies of the [Zr6(μ3-O)8O8H16]8+ node of NU-1000, and we compute their relative energies using a QM/MM scheme designed for the present kind of problem. We compared the results to full quantum mechanical (QM) energy calculations and found that the QM/MM models can reproduce the full QM relative energetics (which span a range of 334 kJ mol-1) with a mean unsigned deviation (MUD) of only 2 kJ mol-1. Furthermore, we found that the structures optimized by QM/MM are nearly identical to their full QM optimized counterparts.

UR - http://www.scopus.com/inward/record.url?scp=85040947492&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85040947492&partnerID=8YFLogxK

U2 - 10.1039/c7cp06751h

DO - 10.1039/c7cp06751h

M3 - Article

VL - 20

SP - 1778

EP - 1786

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 3

ER -