Nature of noncovalent interactions in catenane supramolecular complexes: Calibrating the MM3 force field with ab initio, DFT, and SAPT methods

Tomekia M. Simeon, Mark A Ratner, George C Schatz

Research output: Contribution to journalArticle

21 Citations (Scopus)

Abstract

The design and assembly of mechanically interlocked molecules, such as catenanes and rotaxanes, are dictated by various types of noncovalent interactions. In particular, [C-H···O] hydrogen-bonding and π-π stacking interactions in these supramolecular complexes have been identified as important noncovalent interactions. With this in mind, we examined the [3]catenane 2·4PF6 using molecular mechanics (MM3), ab initio methods (HF, MP2), several versions of density functional theory (DFT) (B3LYP, M0X), and the dispersion-corrected method DFT-D3. Symmetry adapted perturbation theory (DFT-SAPT) provides the highest level of theory considered, and we use the DFT-SAPT results both to calibrate the other electronic structure methods, and the empirical potential MM3 force field that is often used to describe larger catenane and rotaxane structures where [C- H···O] hydrogen-bonding and π-π stacking interactions play a role. Our results indicate that the MM3 calculated complexation energies agree qualitatively with the energetic ordering from DFT-SAPT calculations with an aug-cc-pVTZ basis, both for structures dominated by [C- H···O] hydrogen-bonding and π-π stacking interactions. When the DFT-SAPT energies are decomposed into components, we find that electrostatic interactions dominate the [C-H···O] hydrogen-bonding interactions, while dispersion makes a significant contribution to π-π stacking. Another important conclusion is that DFT-D3 based on M06 or M06-2X provides interaction energies that are in near-quantitative agreement with DFT-SAPT. DFT results without the D3 correction have important differences compared to DFT-SAPT, while HF and even MP2 results are in poor agreement with DFT-SAPT.

Original languageEnglish
Pages (from-to)7918-7927
Number of pages10
JournalJournal of Physical Chemistry A
Volume117
Issue number33
DOIs
Publication statusPublished - Aug 22 2013

Fingerprint

calibrating
field theory (physics)
Density functional theory
density functional theory
interactions
Hydrogen bonds
Rotaxanes
hydrogen
Catenanes
catenane
Molecular mechanics
Coulomb interactions
Complexation
Electronic structure
energy
assembly
perturbation theory
electrostatics
electronic structure
Molecules

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

@article{ad585cd946fd45bb9d55b5104bc45296,
title = "Nature of noncovalent interactions in catenane supramolecular complexes: Calibrating the MM3 force field with ab initio, DFT, and SAPT methods",
abstract = "The design and assembly of mechanically interlocked molecules, such as catenanes and rotaxanes, are dictated by various types of noncovalent interactions. In particular, [C-H···O] hydrogen-bonding and π-π stacking interactions in these supramolecular complexes have been identified as important noncovalent interactions. With this in mind, we examined the [3]catenane 2·4PF6 using molecular mechanics (MM3), ab initio methods (HF, MP2), several versions of density functional theory (DFT) (B3LYP, M0X), and the dispersion-corrected method DFT-D3. Symmetry adapted perturbation theory (DFT-SAPT) provides the highest level of theory considered, and we use the DFT-SAPT results both to calibrate the other electronic structure methods, and the empirical potential MM3 force field that is often used to describe larger catenane and rotaxane structures where [C- H···O] hydrogen-bonding and π-π stacking interactions play a role. Our results indicate that the MM3 calculated complexation energies agree qualitatively with the energetic ordering from DFT-SAPT calculations with an aug-cc-pVTZ basis, both for structures dominated by [C- H···O] hydrogen-bonding and π-π stacking interactions. When the DFT-SAPT energies are decomposed into components, we find that electrostatic interactions dominate the [C-H···O] hydrogen-bonding interactions, while dispersion makes a significant contribution to π-π stacking. Another important conclusion is that DFT-D3 based on M06 or M06-2X provides interaction energies that are in near-quantitative agreement with DFT-SAPT. DFT results without the D3 correction have important differences compared to DFT-SAPT, while HF and even MP2 results are in poor agreement with DFT-SAPT.",
author = "Simeon, {Tomekia M.} and Ratner, {Mark A} and Schatz, {George C}",
year = "2013",
month = "8",
day = "22",
doi = "10.1021/jp400051b",
language = "English",
volume = "117",
pages = "7918--7927",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "33",

}

TY - JOUR

T1 - Nature of noncovalent interactions in catenane supramolecular complexes

T2 - Calibrating the MM3 force field with ab initio, DFT, and SAPT methods

AU - Simeon, Tomekia M.

AU - Ratner, Mark A

AU - Schatz, George C

PY - 2013/8/22

Y1 - 2013/8/22

N2 - The design and assembly of mechanically interlocked molecules, such as catenanes and rotaxanes, are dictated by various types of noncovalent interactions. In particular, [C-H···O] hydrogen-bonding and π-π stacking interactions in these supramolecular complexes have been identified as important noncovalent interactions. With this in mind, we examined the [3]catenane 2·4PF6 using molecular mechanics (MM3), ab initio methods (HF, MP2), several versions of density functional theory (DFT) (B3LYP, M0X), and the dispersion-corrected method DFT-D3. Symmetry adapted perturbation theory (DFT-SAPT) provides the highest level of theory considered, and we use the DFT-SAPT results both to calibrate the other electronic structure methods, and the empirical potential MM3 force field that is often used to describe larger catenane and rotaxane structures where [C- H···O] hydrogen-bonding and π-π stacking interactions play a role. Our results indicate that the MM3 calculated complexation energies agree qualitatively with the energetic ordering from DFT-SAPT calculations with an aug-cc-pVTZ basis, both for structures dominated by [C- H···O] hydrogen-bonding and π-π stacking interactions. When the DFT-SAPT energies are decomposed into components, we find that electrostatic interactions dominate the [C-H···O] hydrogen-bonding interactions, while dispersion makes a significant contribution to π-π stacking. Another important conclusion is that DFT-D3 based on M06 or M06-2X provides interaction energies that are in near-quantitative agreement with DFT-SAPT. DFT results without the D3 correction have important differences compared to DFT-SAPT, while HF and even MP2 results are in poor agreement with DFT-SAPT.

AB - The design and assembly of mechanically interlocked molecules, such as catenanes and rotaxanes, are dictated by various types of noncovalent interactions. In particular, [C-H···O] hydrogen-bonding and π-π stacking interactions in these supramolecular complexes have been identified as important noncovalent interactions. With this in mind, we examined the [3]catenane 2·4PF6 using molecular mechanics (MM3), ab initio methods (HF, MP2), several versions of density functional theory (DFT) (B3LYP, M0X), and the dispersion-corrected method DFT-D3. Symmetry adapted perturbation theory (DFT-SAPT) provides the highest level of theory considered, and we use the DFT-SAPT results both to calibrate the other electronic structure methods, and the empirical potential MM3 force field that is often used to describe larger catenane and rotaxane structures where [C- H···O] hydrogen-bonding and π-π stacking interactions play a role. Our results indicate that the MM3 calculated complexation energies agree qualitatively with the energetic ordering from DFT-SAPT calculations with an aug-cc-pVTZ basis, both for structures dominated by [C- H···O] hydrogen-bonding and π-π stacking interactions. When the DFT-SAPT energies are decomposed into components, we find that electrostatic interactions dominate the [C-H···O] hydrogen-bonding interactions, while dispersion makes a significant contribution to π-π stacking. Another important conclusion is that DFT-D3 based on M06 or M06-2X provides interaction energies that are in near-quantitative agreement with DFT-SAPT. DFT results without the D3 correction have important differences compared to DFT-SAPT, while HF and even MP2 results are in poor agreement with DFT-SAPT.

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

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

U2 - 10.1021/jp400051b

DO - 10.1021/jp400051b

M3 - Article

C2 - 23941280

AN - SCOPUS:84883171167

VL - 117

SP - 7918

EP - 7927

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 33

ER -