Path integral molecular dynamics methods

Pierre A Deymier, Keith A Runge, Ki Dong Oh, G. E. Jabbour

Research output: Chapter in Book/Report/Conference proceedingChapter

Abstract

We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.

Original languageEnglish (US)
Title of host publicationSpringer Series in Materials Science
PublisherSpringer Verlag
Pages13-106
Number of pages94
Volume226
DOIs
StatePublished - 2016

Publication series

NameSpringer Series in Materials Science
Volume226
ISSN (Print)0933033X

Fingerprint

Molecular dynamics
Alkali Metals
Electrons
Alkali metals
Metals
Ions
Quantum theory
Potential energy
Kinetic energy
Polymers
Monte Carlo methods
Phase transitions
Thermodynamics
Plasmas
Heating
Temperature

ASJC Scopus subject areas

  • Materials Science(all)

Cite this

Deymier, P. A., Runge, K. A., Oh, K. D., & Jabbour, G. E. (2016). Path integral molecular dynamics methods. In Springer Series in Materials Science (Vol. 226, pp. 13-106). (Springer Series in Materials Science; Vol. 226). Springer Verlag. https://doi.org/10.1007/978-3-319-24529-4_2

Path integral molecular dynamics methods. / Deymier, Pierre A; Runge, Keith A; Oh, Ki Dong; Jabbour, G. E.

Springer Series in Materials Science. Vol. 226 Springer Verlag, 2016. p. 13-106 (Springer Series in Materials Science; Vol. 226).

Research output: Chapter in Book/Report/Conference proceedingChapter

Deymier, PA, Runge, KA, Oh, KD & Jabbour, GE 2016, Path integral molecular dynamics methods. in Springer Series in Materials Science. vol. 226, Springer Series in Materials Science, vol. 226, Springer Verlag, pp. 13-106. https://doi.org/10.1007/978-3-319-24529-4_2
Deymier PA, Runge KA, Oh KD, Jabbour GE. Path integral molecular dynamics methods. In Springer Series in Materials Science. Vol. 226. Springer Verlag. 2016. p. 13-106. (Springer Series in Materials Science). https://doi.org/10.1007/978-3-319-24529-4_2
Deymier, Pierre A ; Runge, Keith A ; Oh, Ki Dong ; Jabbour, G. E. / Path integral molecular dynamics methods. Springer Series in Materials Science. Vol. 226 Springer Verlag, 2016. pp. 13-106 (Springer Series in Materials Science).
@inbook{b0c8d28cec2c4901be271c78927460f5,
title = "Path integral molecular dynamics methods",
abstract = "We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.",
author = "Deymier, {Pierre A} and Runge, {Keith A} and Oh, {Ki Dong} and Jabbour, {G. E.}",
year = "2016",
doi = "10.1007/978-3-319-24529-4_2",
language = "English (US)",
volume = "226",
series = "Springer Series in Materials Science",
publisher = "Springer Verlag",
pages = "13--106",
booktitle = "Springer Series in Materials Science",

}

TY - CHAP

T1 - Path integral molecular dynamics methods

AU - Deymier, Pierre A

AU - Runge, Keith A

AU - Oh, Ki Dong

AU - Jabbour, G. E.

PY - 2016

Y1 - 2016

N2 - We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.

AB - We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.

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

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

U2 - 10.1007/978-3-319-24529-4_2

DO - 10.1007/978-3-319-24529-4_2

M3 - Chapter

VL - 226

T3 - Springer Series in Materials Science

SP - 13

EP - 106

BT - Springer Series in Materials Science

PB - Springer Verlag

ER -