Procedure for building a consistent embedding at the QM-CM interface

Aditi Mallik, Decarlos E. Taylor, Keith A Runge, James W. Dufty, H. P. Cheng

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

We propose and test a method for a consistent embedding of a domain treated with detailed quantum mechanical methods (QM) inside a domain treated using classical mechanical (CM) potentials. The physical context of this embedding is the response of a system to mechanical strain which leads to fracture. To provide a quantitative test of qualitative ideas, a model system capable of being treated by QM in its entirety is chosen: a silica nano-rod, comprised of 108 atoms. The embedding is constructed so that the CM description yields the same linear response to mechanical strain as the QM description to within a few percent. An acceptable composite representation of the full system requires (1) a CM potential for the classical domain with appropriate linear mechanical response, (2) pseudo-atoms for the termination of dangling bonds at the QM/CM interface, and (3) a dipole description for the polarization of the remainder of the CM region. A key test for the fidelity of the modeled QM domain is accurate forces and electronic charge densities. We show that the composite has the small strain behavior of the full QM treatment of the nano-rod. We find similar success in the application of this method to two other silica model systems: the notched rod, and a nano-ring. Both transfer Hamiltonian (TH) and density functional theory (DFT) are used for the underlying QM.

Original languageEnglish (US)
Pages (from-to)45-60
Number of pages16
JournalJournal of Computer-Aided Materials Design
Volume13
Issue number1-3
DOIs
StatePublished - Oct 2006
Externally publishedYes

Fingerprint

Silicon Dioxide
Silica
Hamiltonians
Atoms
Dangling bonds
Composite materials
Charge density
Density functional theory
Polarization

Keywords

  • Classical pair-potentials
  • Embedding
  • Multi-scale
  • QM/CM
  • Silica

ASJC Scopus subject areas

  • Materials Science(all)
  • Computer Science Applications
  • Computational Theory and Mathematics

Cite this

Procedure for building a consistent embedding at the QM-CM interface. / Mallik, Aditi; Taylor, Decarlos E.; Runge, Keith A; Dufty, James W.; Cheng, H. P.

In: Journal of Computer-Aided Materials Design, Vol. 13, No. 1-3, 10.2006, p. 45-60.

Research output: Contribution to journalArticle

Mallik, Aditi ; Taylor, Decarlos E. ; Runge, Keith A ; Dufty, James W. ; Cheng, H. P. / Procedure for building a consistent embedding at the QM-CM interface. In: Journal of Computer-Aided Materials Design. 2006 ; Vol. 13, No. 1-3. pp. 45-60.
@article{421e25aa73a4431397d04c7cd867937a,
title = "Procedure for building a consistent embedding at the QM-CM interface",
abstract = "We propose and test a method for a consistent embedding of a domain treated with detailed quantum mechanical methods (QM) inside a domain treated using classical mechanical (CM) potentials. The physical context of this embedding is the response of a system to mechanical strain which leads to fracture. To provide a quantitative test of qualitative ideas, a model system capable of being treated by QM in its entirety is chosen: a silica nano-rod, comprised of 108 atoms. The embedding is constructed so that the CM description yields the same linear response to mechanical strain as the QM description to within a few percent. An acceptable composite representation of the full system requires (1) a CM potential for the classical domain with appropriate linear mechanical response, (2) pseudo-atoms for the termination of dangling bonds at the QM/CM interface, and (3) a dipole description for the polarization of the remainder of the CM region. A key test for the fidelity of the modeled QM domain is accurate forces and electronic charge densities. We show that the composite has the small strain behavior of the full QM treatment of the nano-rod. We find similar success in the application of this method to two other silica model systems: the notched rod, and a nano-ring. Both transfer Hamiltonian (TH) and density functional theory (DFT) are used for the underlying QM.",
keywords = "Classical pair-potentials, Embedding, Multi-scale, QM/CM, Silica",
author = "Aditi Mallik and Taylor, {Decarlos E.} and Runge, {Keith A} and Dufty, {James W.} and Cheng, {H. P.}",
year = "2006",
month = "10",
doi = "10.1007/s10820-006-9014-0",
language = "English (US)",
volume = "13",
pages = "45--60",
journal = "Journal of Computer-Aided Materials Design",
issn = "0928-1045",
publisher = "ESCOM Science Publishers",
number = "1-3",

}

TY - JOUR

T1 - Procedure for building a consistent embedding at the QM-CM interface

AU - Mallik, Aditi

AU - Taylor, Decarlos E.

AU - Runge, Keith A

AU - Dufty, James W.

AU - Cheng, H. P.

PY - 2006/10

Y1 - 2006/10

N2 - We propose and test a method for a consistent embedding of a domain treated with detailed quantum mechanical methods (QM) inside a domain treated using classical mechanical (CM) potentials. The physical context of this embedding is the response of a system to mechanical strain which leads to fracture. To provide a quantitative test of qualitative ideas, a model system capable of being treated by QM in its entirety is chosen: a silica nano-rod, comprised of 108 atoms. The embedding is constructed so that the CM description yields the same linear response to mechanical strain as the QM description to within a few percent. An acceptable composite representation of the full system requires (1) a CM potential for the classical domain with appropriate linear mechanical response, (2) pseudo-atoms for the termination of dangling bonds at the QM/CM interface, and (3) a dipole description for the polarization of the remainder of the CM region. A key test for the fidelity of the modeled QM domain is accurate forces and electronic charge densities. We show that the composite has the small strain behavior of the full QM treatment of the nano-rod. We find similar success in the application of this method to two other silica model systems: the notched rod, and a nano-ring. Both transfer Hamiltonian (TH) and density functional theory (DFT) are used for the underlying QM.

AB - We propose and test a method for a consistent embedding of a domain treated with detailed quantum mechanical methods (QM) inside a domain treated using classical mechanical (CM) potentials. The physical context of this embedding is the response of a system to mechanical strain which leads to fracture. To provide a quantitative test of qualitative ideas, a model system capable of being treated by QM in its entirety is chosen: a silica nano-rod, comprised of 108 atoms. The embedding is constructed so that the CM description yields the same linear response to mechanical strain as the QM description to within a few percent. An acceptable composite representation of the full system requires (1) a CM potential for the classical domain with appropriate linear mechanical response, (2) pseudo-atoms for the termination of dangling bonds at the QM/CM interface, and (3) a dipole description for the polarization of the remainder of the CM region. A key test for the fidelity of the modeled QM domain is accurate forces and electronic charge densities. We show that the composite has the small strain behavior of the full QM treatment of the nano-rod. We find similar success in the application of this method to two other silica model systems: the notched rod, and a nano-ring. Both transfer Hamiltonian (TH) and density functional theory (DFT) are used for the underlying QM.

KW - Classical pair-potentials

KW - Embedding

KW - Multi-scale

KW - QM/CM

KW - Silica

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

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

U2 - 10.1007/s10820-006-9014-0

DO - 10.1007/s10820-006-9014-0

M3 - Article

AN - SCOPUS:33749004293

VL - 13

SP - 45

EP - 60

JO - Journal of Computer-Aided Materials Design

JF - Journal of Computer-Aided Materials Design

SN - 0928-1045

IS - 1-3

ER -