Problem with using EFIELD_PEAD for MD and relaxation

[huangj3]

With EFIELD_PEAD given in the INCAR, the first SCF loop looks ok. It shows the SDA and DMP lines. However, starting from the 2nd ionic step, the electric field is probably not enforced. By doing an ionic relaxation with and without the electric field, I saw that the first ionic step energies are different but the following steps are the same, indicating that the EFIELD_PEAD wasn't enforced in an ionic relaxation from the 2nd step. This is also true for an MD simulation. Am I doing the calculations correctly or EFIELD_PEAD is only meant for static calculation?

Thanks a lot.

Huangj3

[Li Huang]

I have exactly the same question. Any reply would be highly appreciated!