Stable absorbing boundary conditions for molecular dynamics in general domains

Xiaojie Wu, Xiantao Li

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

A systematic strategy of absorbing boundary conditions for molecular dynamics simulations are presented. The exact boundary conditions for crystalline solids with harmonic approximation are expressed as a dynamic Dirichlet-to-Neumann (DtN) map. It connects the displacement of the atoms at the boundary to the traction on these atoms. The DtN map is valid for a domain with general geometry. To avoid evaluating the time convolution of the dynamic DtN map, we approximate the associated kernel function by rational functions in the Laplace domain. The parameters in the approximations are determined by interpolations. The explicit forms of the zeroth, first, and second order approximations will be presented. The stability of the molecular dynamics model, supplemented with these absorbing boundary conditions is established. Two numerical simulations are performed to demonstrate the effectiveness of the methods.

Original languageEnglish (US)
Pages (from-to)1259-1272
Number of pages14
JournalComputational Mechanics
Volume62
Issue number6
DOIs
StatePublished - Dec 1 2018

Fingerprint

Dirichlet-to-Neumann Map
Absorbing Boundary Conditions
Molecular Dynamics
Molecular dynamics
Boundary conditions
Atoms
Second-order Approximation
Rational functions
Zeroth
Computer simulation
Approximation
Kernel Function
Laplace
Convolution
Rational function
Molecular Dynamics Simulation
Dynamic models
Dynamic Model
Interpolation
Harmonic

All Science Journal Classification (ASJC) codes

  • Computational Mechanics
  • Ocean Engineering
  • Mechanical Engineering
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics

Cite this

@article{7af1c1b2c1e4466291c411eaf8c9c53e,
title = "Stable absorbing boundary conditions for molecular dynamics in general domains",
abstract = "A systematic strategy of absorbing boundary conditions for molecular dynamics simulations are presented. The exact boundary conditions for crystalline solids with harmonic approximation are expressed as a dynamic Dirichlet-to-Neumann (DtN) map. It connects the displacement of the atoms at the boundary to the traction on these atoms. The DtN map is valid for a domain with general geometry. To avoid evaluating the time convolution of the dynamic DtN map, we approximate the associated kernel function by rational functions in the Laplace domain. The parameters in the approximations are determined by interpolations. The explicit forms of the zeroth, first, and second order approximations will be presented. The stability of the molecular dynamics model, supplemented with these absorbing boundary conditions is established. Two numerical simulations are performed to demonstrate the effectiveness of the methods.",
author = "Xiaojie Wu and Xiantao Li",
year = "2018",
month = "12",
day = "1",
doi = "10.1007/s00466-018-1562-0",
language = "English (US)",
volume = "62",
pages = "1259--1272",
journal = "Computational Mechanics",
issn = "0178-7675",
publisher = "Springer Verlag",
number = "6",

}

Stable absorbing boundary conditions for molecular dynamics in general domains. / Wu, Xiaojie; Li, Xiantao.

In: Computational Mechanics, Vol. 62, No. 6, 01.12.2018, p. 1259-1272.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Stable absorbing boundary conditions for molecular dynamics in general domains

AU - Wu, Xiaojie

AU - Li, Xiantao

PY - 2018/12/1

Y1 - 2018/12/1

N2 - A systematic strategy of absorbing boundary conditions for molecular dynamics simulations are presented. The exact boundary conditions for crystalline solids with harmonic approximation are expressed as a dynamic Dirichlet-to-Neumann (DtN) map. It connects the displacement of the atoms at the boundary to the traction on these atoms. The DtN map is valid for a domain with general geometry. To avoid evaluating the time convolution of the dynamic DtN map, we approximate the associated kernel function by rational functions in the Laplace domain. The parameters in the approximations are determined by interpolations. The explicit forms of the zeroth, first, and second order approximations will be presented. The stability of the molecular dynamics model, supplemented with these absorbing boundary conditions is established. Two numerical simulations are performed to demonstrate the effectiveness of the methods.

AB - A systematic strategy of absorbing boundary conditions for molecular dynamics simulations are presented. The exact boundary conditions for crystalline solids with harmonic approximation are expressed as a dynamic Dirichlet-to-Neumann (DtN) map. It connects the displacement of the atoms at the boundary to the traction on these atoms. The DtN map is valid for a domain with general geometry. To avoid evaluating the time convolution of the dynamic DtN map, we approximate the associated kernel function by rational functions in the Laplace domain. The parameters in the approximations are determined by interpolations. The explicit forms of the zeroth, first, and second order approximations will be presented. The stability of the molecular dynamics model, supplemented with these absorbing boundary conditions is established. Two numerical simulations are performed to demonstrate the effectiveness of the methods.

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

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

U2 - 10.1007/s00466-018-1562-0

DO - 10.1007/s00466-018-1562-0

M3 - Article

VL - 62

SP - 1259

EP - 1272

JO - Computational Mechanics

JF - Computational Mechanics

SN - 0178-7675

IS - 6

ER -