<div><font face="Arial">Dear CP2K users,<br /><br />I am calculating the energy of a slab (periodic in XY and non-periodic in Z) using CP2K. I have included my original CP2K input file for reference.<br /><br />The problem is that I cannot achieve convergence with UKS. Perhaps it's due to my impatience, but the convergence cannot go below 0.01 after a few hundred iterations.<br /><br />SCAN (including R2SCAN) is known for its difficulty in converging. However, for silicon, it converges well with the same settings, though it doesn't require UKS (this was just for experimentation).<br /><br />After changing the POISSON_SOLVER to IMPLICIT, I observed a slight improvement in the convergence profile. However, it still did not fully converge.<br /><br />Upon adding different options like the ones below, I saw a significant improvement<br /><br /> &POISSON<br />   PERIODIC XY<br />   POISSON_SOLVER IMPLICIT<br />   &IMPLICIT<br />    BOUNDARY_CONDITIONS NEUMANN<br />    NEUMANN_DIRECTIONS Z<br />   &END IMPLICIT<br />  &END POISSON<br /><br />However, there is a discrepancy between the energy scale where the original script struggled (around -17794 Hartree) and the one shown by this modification (-18097 Hartree). The simulation is still running but seems promising, at least in terms of convergence trend through iterations.<br />My questions are:<br /><br />(1) Are the boundary conditions shown above appropriate for solid metal slab calculations? The energy discrepancy induced by adding options makes me suspicious of using an improper method. Or could it be that we've finally found the correct energy level?<br />(2) Does my original code have an intrinsic problem that prevents convergence, or is it simply a matter of needing more patience to allow it to converge?<br /><br />Thanks in advance,<br /><br />Ji Woong</font><br /></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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d664aa41-d272-4e00-96de-70cccaf63af4n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d664aa41-d272-4e00-96de-70cccaf63af4n%40googlegroups.com</a>.<br />