Dear cp2k users/developers,<div><br></div><div>I am wondering if the function of applying electric field through "periodic_efield" was designed to also work with DFTB.  I simply took the H2O-efiled.inp from tests/QS/regtest-p-efield and modified it to use DFTB (parameters from data/DFTB/scc).  I was using version 7.1.</div><div><br></div><div>I am comparing the Mulliken charges from cp2k-DFT, cp2k-DFTB, and also from DFTB+.  Apparently, results from cp2k-DFTB under the field are a little odd.</div><div><br></div><div>Thanks, Wei</div><div><br></div><div>efield:                                                               -0.01                  0                         0.01</div><div>Mulliken charge (cp2k-DFT):                 -0.389706          -0.426154           -0.461808<br></div><div>Mulliken charge (cp2k-DFTB):               0.060747          -0.547117          -1.089736<br></div><div>Mulliken charge (DFTB+):                       -0.525336         -0.547117          -0.568806 </div><div><br></div><div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    &PRINT</div><div>      &MOMENTS</div><div>        PERIODIC T</div><div>      &END</div><div>    &END</div><div>    SURFACE_DIPOLE_CORRECTION</div><div>    SURF_DIP_DIR Z</div><div>    &PERIODIC_EFIELD</div><div>       POLARISATION 0 0 1</div><div>       INTENSITY -0.01</div><div>    &END    </div><div>    &QS</div><div>      METHOD DFTB</div><div>      &DFTB</div><div>        SELF_CONSISTENT    T</div><div>        DO_EWALD           T</div><div>        DISPERSION         F</div><div>        &PARAMETER</div><div>            SK_FILE H H     ./H-H.skf</div><div>            SK_FILE H O     ./H-O.skf</div><div>            SK_FILE O H     ./O-H.skf</div><div>            SK_FILE O O     ./O-O.skf</div><div>        &END PARAMETER</div><div>      &END DFTB</div><div>    &END QS</div><div>    &SCF</div><div>      EPS_SCF 1.0E-6</div><div>      SCF_GUESS MOPAC</div><div>      MAX_SCF 100</div><div>      &OT</div><div>        PRECONDITIONER FULL_ALL</div><div>      &END</div><div>    &END SCF</div><div>    &POISSON</div><div>      PERIODIC XYZ</div><div>      &EWALD</div><div>        EWALD_TYPE SPME</div><div>        ALPHA  .35</div><div>        GMAX   41</div><div>      &END EWALD        </div><div>    &END POISSON</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &TOPOLOGY</div><div>    &CENTER_COORDINATES</div><div>    &END</div><div>    &END</div><div>    &CELL</div><div>      ABC 16.0 16.0 16.0</div><div>    &END CELL</div><div>    &COORD</div><div>    O   0.000000    0.000000   -0.065587 H2O</div><div>    H   0.000000   -0.757136    0.520545 H2O</div><div>    H   0.000000    0.757136    0.520545 H2O</div><div>    &END COORD</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div>  PROJECT H2O-field</div><div>  PRINT_LEVEL MEDIUM</div><div>  RUN_TYPE ENERGY</div><div>&END GLOBAL</div></div><div><br></div><div><br></div><div><br></div>