Electric field force on nuclei

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
aohara
Newbie
Newbie
Posts: 13
Joined: Sat Apr 16, 2016 5:54 am
License Nr.: 5-225

Electric field force on nuclei

#1 Post by aohara » Thu Jul 11, 2019 6:44 pm

A colleague and I were recently discussing some calculations on relaxation and molecular dynamics of slabs under the use of an electric field (using the saw-tooth implementation of EFIELD / LDIPOL). From my reading through of the manual and also the code (specifically in dipol.F), the only contribution to the forces is through the dipole term added to the electronic contribution in the Hellman-Feynman forces. However, as far as I'm aware, there is no inclusion of the force due to the direct interaction of the electric field with the pseudo-nuclei. Is there a reasonable straightforward way to include this interaction even at a semi-classical level? Otherwise, is it correct to assume that nuclei only move due to their following the perturbation of the electron charge density due to the field? Thanks for any advice!

Post Reply