<div dir="ltr">Dear CP2K developers and users,<div><br></div><div>I've tried to compute the binding energy of Na _ Cl (singlet state so the Multiplicity should be 1)but I got a strange curve at large distances. I don't have a lot o experience in using cp2k so I attached the script that I've used to calculate the potential energy at varying distance. Could someone check if I did something wrong? any help would be appreciated. </div><div><br></div><div>&GLOBAL</div><div>    PROJECT NaCl</div><div>    RUN_TYPE ENERGY             </div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div><br></div><div>  METHOD Quickstep</div><div>  &DFT</div><div>      BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>      POTENTIAL_FILE_NAME GTH_POTENTIALS            </div><div>      CHARGE 0</div><div>      MULTIPLICITY 1</div><div><br></div><div>    &MGRID</div><div>       CUTOFF [Ry] 400 </div><div>    &END</div><div><br></div><div>    &QS</div><div>       METHOD GPW </div><div>       EPS_DEFAULT 1.0E-6 </div><div>    &END</div><div><br></div><div>    &POISSON                      # POISSON solver for non-periodic calculation</div><div>       PERIODIC NONE</div><div>       PSOLVER WAVELET</div><div>    &END</div><div>    &SCF                              </div><div>      SCF_GUESS ATOMIC ! can be used to RESTART an interrupted calculation</div><div>      MAX_SCF 300</div><div>      EPS_SCF 1.0E-6 ! accuracy of the SCF procedure typically 1.0E-6 - 1.0E-7</div><div>      &OT</div><div>        PRECONDITIONER FULL_SINGLE_INVERSE</div><div>        MINIMIZER DIIS</div><div>      &END OT</div><div>      &MIXING</div><div>         METHOD BROYDEN_MIXING</div><div>         ALPHA  0.4</div><div>         BETA   1.5</div><div>         NBROYDEN 8</div><div>      &END MIXING</div><div>      &OUTER_SCF ! repeat the inner SCF cycle 10 times<br></div><div>        MAX_SCF 100</div><div>        EPS_SCF 1.0E-5 ! must match the above</div><div>        OPTIMIZER  DIIS</div><div>      &END</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE<br></div><div>      &END XC_FUNCTIONAL </div><div>    &END XC</div><div>  &END DFT</div><div>   &SUBSYS</div><div>    &CELL </div><div>      ABC [angstrom] 10.00 10.00 10.00<br></div><div>    &END CELL</div><div>    &COORD</div><div>Na    0.0 0.0 0.0 </div><div>Cl    0.0 0.0 MYDIST</div><div>    &END COORD</div><div>    &TOPOLOGY</div><div>      CONNECTIVITY GENERATE</div><div>      &GENERATE</div><div>        BONDLENGTH_MAX 9</div><div>      &END</div><div>    &END</div><div>    &KIND Na                              </div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH         </div><div>      POTENTIAL GTH-PBE          </div><div>    &END KIND</div><div>    &KIND Cl</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><div><br></div><div><br></div></div><div>The bash script used to get the calculation running is the following: </div><div><br></div><div><div>rm -f ener_profile</div><div><br></div><div>#for dist in `seq 2.3 0.1 6.0`</div><div>for dist in 4.0 5.0 6.0 7.0 8.0 10.0; do</div><div><br></div><div>  #</div><div>  # compute energy </div><div>  # </div><div>  echo "Computing potential energy for distance $dist"</div><div>  sed "s/MYDIST/${dist}/g" mode1.inp > inp</div><div>  cp2k.popt inp > out </div><div>  ener=`grep ' ENERGY| Total FORCE_EVAL' out | awk '{print $NF}'`</div><div>  echo $dist $ener >> ener_profile</div><div><br></div><div>done</div></div><div><br></div></div>