Efficient solution of Poisson's equation using discrete variable representation basis sets for Car-Parrinello ab initio molecular dynamics simulations with cluster boundary conditions

Hee Seung Lee, Mark Tuckerman

Research output: Contribution to journalArticle

Abstract

An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu, Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O (N 43)] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N2 and a water dimer.

Original languageEnglish (US)
Article number224108
JournalJournal of Chemical Physics
Volume129
Issue number22
DOIs
StatePublished - 2008

Fingerprint

Poisson equation
Molecular dynamics
Railroad cars
Boundary conditions
boundary conditions
molecular dynamics
Computer simulation
simulation
Dimers
Carrier concentration
free boundaries
Water
grids
dimers
Costs
formalism
costs
scaling
orbitals
cells

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

@article{e6fb481d8ea645ada0418affe34541ee,
title = "Efficient solution of Poisson's equation using discrete variable representation basis sets for Car-Parrinello ab initio molecular dynamics simulations with cluster boundary conditions",
abstract = "An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu, Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O (N 43)] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N2 and a water dimer.",
author = "Lee, {Hee Seung} and Mark Tuckerman",
year = "2008",
doi = "10.1063/1.3036423",
language = "English (US)",
volume = "129",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "22",

}

TY - JOUR

T1 - Efficient solution of Poisson's equation using discrete variable representation basis sets for Car-Parrinello ab initio molecular dynamics simulations with cluster boundary conditions

AU - Lee, Hee Seung

AU - Tuckerman, Mark

PY - 2008

Y1 - 2008

N2 - An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu, Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O (N 43)] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N2 and a water dimer.

AB - An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu, Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O (N 43)] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N2 and a water dimer.

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

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

U2 - 10.1063/1.3036423

DO - 10.1063/1.3036423

M3 - Article

C2 - 19071908

AN - SCOPUS:57649195670

VL - 129

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 22

M1 - 224108

ER -