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 E. 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 - Jan 1 2008

ASJC Scopus subject areas

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

Fingerprint Dive into the research topics of 'Efficient solution of Poisson's equation using discrete variable representation basis sets for Car-Parrinello ab initio molecular dynamics simulations with cluster boundary conditions'. Together they form a unique fingerprint.

  • Cite this