TY - JOUR

T1 - Electron nuclear dynamics with plane wave basis sets

T2 - complete theory and formalism

AU - Teixeira, Erico S.

AU - Morales, Jorge A.

N1 - Publisher Copyright:
© 2020, Springer-Verlag GmbH Germany, part of Springer Nature.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.

PY - 2020/4/1

Y1 - 2020/4/1

N2 - Electron nuclear dynamics (END) is an ab initio quantum dynamics method that adopts a time-dependent, variational, direct, and non-adiabatic approach. The simplest-level (SL) END (SLEND) version employs a classical mechanics description for nuclei and a Thouless single-determinantal wave function for electrons. A higher-level END version, END/Kohn–Sham density functional theory, improves the electron correlation description of SLEND. While both versions can simulate various types of chemical reactions, they have difficulties to simulate scattering/capture of electrons to/from the continuum due to their reliance on localized Slater-type basis functions. To properly describe those processes, we formulate END with plane waves (PWs, END/PW), basis functions able to represent both bound and unbound electrons. As extra benefits, PWs also afford fast algorithms to simulate periodic systems, parametric independence from nuclear positions and momenta, and elimination of basis set linear dependencies and orthogonalization procedures. We obtain the END/PW formalism by extending the Thouless wave function and associated electron density to periodic systems, expressing the energy terms as functionals of the latter entities, and deriving the energy gradients with respect to nuclear and electronic variables. END/PW has a great potential to simulate electron processes in both periodic (crystal) and aperiodic (molecular) systems (the latter in a supercell approach). Following previous END studies, END/PW will be applied to electron scattering processes in proton cancer therapy reactions.

AB - Electron nuclear dynamics (END) is an ab initio quantum dynamics method that adopts a time-dependent, variational, direct, and non-adiabatic approach. The simplest-level (SL) END (SLEND) version employs a classical mechanics description for nuclei and a Thouless single-determinantal wave function for electrons. A higher-level END version, END/Kohn–Sham density functional theory, improves the electron correlation description of SLEND. While both versions can simulate various types of chemical reactions, they have difficulties to simulate scattering/capture of electrons to/from the continuum due to their reliance on localized Slater-type basis functions. To properly describe those processes, we formulate END with plane waves (PWs, END/PW), basis functions able to represent both bound and unbound electrons. As extra benefits, PWs also afford fast algorithms to simulate periodic systems, parametric independence from nuclear positions and momenta, and elimination of basis set linear dependencies and orthogonalization procedures. We obtain the END/PW formalism by extending the Thouless wave function and associated electron density to periodic systems, expressing the energy terms as functionals of the latter entities, and deriving the energy gradients with respect to nuclear and electronic variables. END/PW has a great potential to simulate electron processes in both periodic (crystal) and aperiodic (molecular) systems (the latter in a supercell approach). Following previous END studies, END/PW will be applied to electron scattering processes in proton cancer therapy reactions.

KW - Electron nuclear dynamics

KW - Plane wave basis sets

KW - Proton cancer therapy

KW - Thouless state

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

U2 - 10.1007/s00214-020-2578-z

DO - 10.1007/s00214-020-2578-z

M3 - Article

AN - SCOPUS:85083104875

VL - 139

JO - Theoretical Chemistry Accounts

JF - Theoretical Chemistry Accounts

SN - 1432-881X

IS - 4

M1 - 73

ER -