Multiscale coupling schemes spanning the quantum mechanical, atomistic forcefield, and continuum regimes

Roopam Khare, Steven L. Mielke, George C Schatz, Ted Belytschko

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

Strategies for coupling quantum mechanical (QM), molecular mechanical (MM), and continuum mechanical (CM) methods are described. For QM/MM coupling, we consider two overlapping domain schemes: the widely-used ONIOM method that involves full overlap, and a new minimal-overlap scheme denoted as the "quantum to molecular mechanical overlapping domain" (QtMMOD) method. In the QtMMOD method, the part of the region that is treated by quantum mechanics does not require MM calculations, which makes the method suitable for systems containing localized regions that cannot be modeled by available empirical potentials. We describe how the QtMMOD method can in turn be coupled to a finite element model by the bridging domain method, another overlapping domain method, yielding a QM/MM/CM model. We also show how the QtMMOD scheme can be used to couple a finite element model directly with a QM model. These coupling methods were used to calculate the fracture properties of graphene sheets containing defects. We obtained good agreement between stress-strain curves and fracture strengths calculated by the three methods and benchmark results of a strictly QM calculation. An extension of an atomistic strain calculation method based on a moving least squares scheme, which is applicable to the treatment of material interfaces, is given.

Original languageEnglish
Pages (from-to)3190-3202
Number of pages13
JournalComputer Methods in Applied Mechanics and Engineering
Volume197
Issue number41-42
DOIs
Publication statusPublished - Jul 1 2008

Fingerprint

continuums
Quantum theory
Stress-strain curves
Graphene
Fracture toughness
Defects
fracture strength
quantum mechanics
graphene
defects
curves

Keywords

  • Fracture
  • Graphene sheets
  • Molecular mechanics
  • Multiscale
  • Nanofracture
  • Quantum mechanics

ASJC Scopus subject areas

  • Computer Science Applications
  • Computational Mechanics

Cite this

Multiscale coupling schemes spanning the quantum mechanical, atomistic forcefield, and continuum regimes. / Khare, Roopam; Mielke, Steven L.; Schatz, George C; Belytschko, Ted.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 197, No. 41-42, 01.07.2008, p. 3190-3202.

Research output: Contribution to journalArticle

@article{04056f4e94e849bd9475b49e4fdc41fd,
title = "Multiscale coupling schemes spanning the quantum mechanical, atomistic forcefield, and continuum regimes",
abstract = "Strategies for coupling quantum mechanical (QM), molecular mechanical (MM), and continuum mechanical (CM) methods are described. For QM/MM coupling, we consider two overlapping domain schemes: the widely-used ONIOM method that involves full overlap, and a new minimal-overlap scheme denoted as the {"}quantum to molecular mechanical overlapping domain{"} (QtMMOD) method. In the QtMMOD method, the part of the region that is treated by quantum mechanics does not require MM calculations, which makes the method suitable for systems containing localized regions that cannot be modeled by available empirical potentials. We describe how the QtMMOD method can in turn be coupled to a finite element model by the bridging domain method, another overlapping domain method, yielding a QM/MM/CM model. We also show how the QtMMOD scheme can be used to couple a finite element model directly with a QM model. These coupling methods were used to calculate the fracture properties of graphene sheets containing defects. We obtained good agreement between stress-strain curves and fracture strengths calculated by the three methods and benchmark results of a strictly QM calculation. An extension of an atomistic strain calculation method based on a moving least squares scheme, which is applicable to the treatment of material interfaces, is given.",
keywords = "Fracture, Graphene sheets, Molecular mechanics, Multiscale, Nanofracture, Quantum mechanics",
author = "Roopam Khare and Mielke, {Steven L.} and Schatz, {George C} and Ted Belytschko",
year = "2008",
month = "7",
day = "1",
doi = "10.1016/j.cma.2007.11.029",
language = "English",
volume = "197",
pages = "3190--3202",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0374-2830",
publisher = "Elsevier",
number = "41-42",

}

TY - JOUR

T1 - Multiscale coupling schemes spanning the quantum mechanical, atomistic forcefield, and continuum regimes

AU - Khare, Roopam

AU - Mielke, Steven L.

AU - Schatz, George C

AU - Belytschko, Ted

PY - 2008/7/1

Y1 - 2008/7/1

N2 - Strategies for coupling quantum mechanical (QM), molecular mechanical (MM), and continuum mechanical (CM) methods are described. For QM/MM coupling, we consider two overlapping domain schemes: the widely-used ONIOM method that involves full overlap, and a new minimal-overlap scheme denoted as the "quantum to molecular mechanical overlapping domain" (QtMMOD) method. In the QtMMOD method, the part of the region that is treated by quantum mechanics does not require MM calculations, which makes the method suitable for systems containing localized regions that cannot be modeled by available empirical potentials. We describe how the QtMMOD method can in turn be coupled to a finite element model by the bridging domain method, another overlapping domain method, yielding a QM/MM/CM model. We also show how the QtMMOD scheme can be used to couple a finite element model directly with a QM model. These coupling methods were used to calculate the fracture properties of graphene sheets containing defects. We obtained good agreement between stress-strain curves and fracture strengths calculated by the three methods and benchmark results of a strictly QM calculation. An extension of an atomistic strain calculation method based on a moving least squares scheme, which is applicable to the treatment of material interfaces, is given.

AB - Strategies for coupling quantum mechanical (QM), molecular mechanical (MM), and continuum mechanical (CM) methods are described. For QM/MM coupling, we consider two overlapping domain schemes: the widely-used ONIOM method that involves full overlap, and a new minimal-overlap scheme denoted as the "quantum to molecular mechanical overlapping domain" (QtMMOD) method. In the QtMMOD method, the part of the region that is treated by quantum mechanics does not require MM calculations, which makes the method suitable for systems containing localized regions that cannot be modeled by available empirical potentials. We describe how the QtMMOD method can in turn be coupled to a finite element model by the bridging domain method, another overlapping domain method, yielding a QM/MM/CM model. We also show how the QtMMOD scheme can be used to couple a finite element model directly with a QM model. These coupling methods were used to calculate the fracture properties of graphene sheets containing defects. We obtained good agreement between stress-strain curves and fracture strengths calculated by the three methods and benchmark results of a strictly QM calculation. An extension of an atomistic strain calculation method based on a moving least squares scheme, which is applicable to the treatment of material interfaces, is given.

KW - Fracture

KW - Graphene sheets

KW - Molecular mechanics

KW - Multiscale

KW - Nanofracture

KW - Quantum mechanics

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

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

U2 - 10.1016/j.cma.2007.11.029

DO - 10.1016/j.cma.2007.11.029

M3 - Article

VL - 197

SP - 3190

EP - 3202

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0374-2830

IS - 41-42

ER -