<div dir="ltr"><div>Dear Mejdeddine,</div><div><br></div><div>You are using a non-periodic Poisson solver, suitable for this system of isolated atoms. However, the cell dimensions and atomic coordinates have to be chosen in such a way that the electronic density on the edges of the unit cell is negligible. Keep in mind that the cubic unit cell is situated with one corner in location (0, 0, 0) and the opposite corner in (10, 10, 10) Center the atoms in the middle of the cell, either manually, or by switching on the <code class="gmail-upper"><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY/CENTER_COORDINATES.html">CENTER_COORDINATES</a> <span style="font-family:arial,sans-serif">statement in the</span><code class="gmail-upper"><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY.html"> TOPOLOGY</a></code></code><code> <span style="font-family:arial,sans-serif">section, which should center them in the middle of the unit cell by default.<br></span></code></div><div><code><span style="font-family:arial,sans-serif"><br></span></code></div><div><code><span style="font-family:arial,sans-serif">Apart from that, I don't see any other immediate issue. What binding energy plot do you get? What is strange about it? <br></span></code></div><div><code><span style="font-family:arial,sans-serif"><br></span></code></div><div><code><span style="font-family:arial,sans-serif">Yours sincerely,</span></code></div><div><code><span style="font-family:arial,sans-serif">Patrick Gono</span></code><br></div><div><code class="gmail-upper"><code class="gmail-upper"></code></code></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 6 Jul 2020 at 14:56, mejdeddine mokhtar <<a href="mailto:mejdit...@gmail.com">mejdit...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0798c015-5b94-4175-9179-71f3c2a5b5e3o%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/0798c015-5b94-4175-9179-71f3c2a5b5e3o%40googlegroups.com</a>.<br>
</blockquote></div>