Ab initio molecular dynamics with discrete variable representation basis sets

Techniques and application to liquid water

Hee Seung Lee, Mark Tuckerman

Research output: Contribution to journalArticle

Abstract

Finite temperature ab initio molecular dynamics (AIMD), in which forces are obtained from "on-the-fly" electronic structure calculations, is a widely used technique for studying structural and dynamical properties of chemically active systems. Recently, we introduced an AIMD scheme based on discrete variable representation (DVR) basis sets, which was shown to have improved convergence properties over the conventional plane wave (PW) basis set [Liu,Y.; et al. Phys. Rev. B 2003, 68, 125110]. In the present work, the numerical algorithms for the DVR based AIMD scheme (DVR/AIMD) are provided in detail, and the latest developments of the approach are presented. The accuracy and stability of the current implementation of the DVR/AIMD scheme are tested by performing a simulation of liquid water at ambient conditions. The structural information obtained from the present work is in good agreement with the result of recent AIMD simulations with a PW basis set (PW/AIMD). Advantages of using the DVR/AIMD scheme over the PW/ AIMD method are discussed. In particular, it is shown that a DVR/AIMD simulation of liquid water in the complete basis set limit is possible with a relatively small number of grid points.

Original languageEnglish (US)
Pages (from-to)5549-5560
Number of pages12
JournalJournal of Physical Chemistry A
Volume110
Issue number16
DOIs
StatePublished - Apr 27 2006

Fingerprint

Molecular dynamics
molecular dynamics
Water
Liquids
liquids
water
plane waves
simulation
Computer simulation
Electronic structure
grids
electronic structure

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Ab initio molecular dynamics with discrete variable representation basis sets : Techniques and application to liquid water. / Lee, Hee Seung; Tuckerman, Mark.

In: Journal of Physical Chemistry A, Vol. 110, No. 16, 27.04.2006, p. 5549-5560.

Research output: Contribution to journalArticle

@article{28bcf253c22a4bf6bf49eafc01a6425a,
title = "Ab initio molecular dynamics with discrete variable representation basis sets: Techniques and application to liquid water",
abstract = "Finite temperature ab initio molecular dynamics (AIMD), in which forces are obtained from {"}on-the-fly{"} electronic structure calculations, is a widely used technique for studying structural and dynamical properties of chemically active systems. Recently, we introduced an AIMD scheme based on discrete variable representation (DVR) basis sets, which was shown to have improved convergence properties over the conventional plane wave (PW) basis set [Liu,Y.; et al. Phys. Rev. B 2003, 68, 125110]. In the present work, the numerical algorithms for the DVR based AIMD scheme (DVR/AIMD) are provided in detail, and the latest developments of the approach are presented. The accuracy and stability of the current implementation of the DVR/AIMD scheme are tested by performing a simulation of liquid water at ambient conditions. The structural information obtained from the present work is in good agreement with the result of recent AIMD simulations with a PW basis set (PW/AIMD). Advantages of using the DVR/AIMD scheme over the PW/ AIMD method are discussed. In particular, it is shown that a DVR/AIMD simulation of liquid water in the complete basis set limit is possible with a relatively small number of grid points.",
author = "Lee, {Hee Seung} and Mark Tuckerman",
year = "2006",
month = "4",
day = "27",
doi = "10.1021/jp0570770",
language = "English (US)",
volume = "110",
pages = "5549--5560",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "16",

}

TY - JOUR

T1 - Ab initio molecular dynamics with discrete variable representation basis sets

T2 - Techniques and application to liquid water

AU - Lee, Hee Seung

AU - Tuckerman, Mark

PY - 2006/4/27

Y1 - 2006/4/27

N2 - Finite temperature ab initio molecular dynamics (AIMD), in which forces are obtained from "on-the-fly" electronic structure calculations, is a widely used technique for studying structural and dynamical properties of chemically active systems. Recently, we introduced an AIMD scheme based on discrete variable representation (DVR) basis sets, which was shown to have improved convergence properties over the conventional plane wave (PW) basis set [Liu,Y.; et al. Phys. Rev. B 2003, 68, 125110]. In the present work, the numerical algorithms for the DVR based AIMD scheme (DVR/AIMD) are provided in detail, and the latest developments of the approach are presented. The accuracy and stability of the current implementation of the DVR/AIMD scheme are tested by performing a simulation of liquid water at ambient conditions. The structural information obtained from the present work is in good agreement with the result of recent AIMD simulations with a PW basis set (PW/AIMD). Advantages of using the DVR/AIMD scheme over the PW/ AIMD method are discussed. In particular, it is shown that a DVR/AIMD simulation of liquid water in the complete basis set limit is possible with a relatively small number of grid points.

AB - Finite temperature ab initio molecular dynamics (AIMD), in which forces are obtained from "on-the-fly" electronic structure calculations, is a widely used technique for studying structural and dynamical properties of chemically active systems. Recently, we introduced an AIMD scheme based on discrete variable representation (DVR) basis sets, which was shown to have improved convergence properties over the conventional plane wave (PW) basis set [Liu,Y.; et al. Phys. Rev. B 2003, 68, 125110]. In the present work, the numerical algorithms for the DVR based AIMD scheme (DVR/AIMD) are provided in detail, and the latest developments of the approach are presented. The accuracy and stability of the current implementation of the DVR/AIMD scheme are tested by performing a simulation of liquid water at ambient conditions. The structural information obtained from the present work is in good agreement with the result of recent AIMD simulations with a PW basis set (PW/AIMD). Advantages of using the DVR/AIMD scheme over the PW/ AIMD method are discussed. In particular, it is shown that a DVR/AIMD simulation of liquid water in the complete basis set limit is possible with a relatively small number of grid points.

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

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

U2 - 10.1021/jp0570770

DO - 10.1021/jp0570770

M3 - Article

VL - 110

SP - 5549

EP - 5560

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 16

ER -