Hi, CP2K developers and users,<br><br> I think a bug likely exits in "efield_utils.F" that handles the finite external electric fields. In the subroutine "make_field", the statement (line 152 of the latest CVS version) of <br>
<br> "field = 0.0_dp" should be changed to "field = field + 0.0_dp" (or simply be removed)<br><br> The bugfix allows the appropriate handling of multiple sources of external electric field, at least one of which is with constant envelope.<br>
<br> Please correct me if I am wrong.<br><br> Thanks.<br><br>Hanning Chen, Ph.D.,<br>Department of Chemistry<br>Northwestern University<br>Evanston, IL 60091<br>U.S.A<br> <br><br><br><br>