[CP2K-user] [CP2K:20427] UKS + BCC Iron ( setup and convergence )
Ji Woong Yu
jiwoongs0510 at gmail.com
Sat Jul 6 13:39:54 UTC 2024
Dear CP2K users,
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.
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.
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).
After changing the POISSON_SOLVER to IMPLICIT, I observed a slight
improvement in the convergence profile. However, it still did not fully
converge.
Upon adding different options like the ones below, I saw a significant
improvement
&POISSON
PERIODIC XY
POISSON_SOLVER IMPLICIT
&IMPLICIT
BOUNDARY_CONDITIONS NEUMANN
NEUMANN_DIRECTIONS Z
&END IMPLICIT
&END POISSON
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.
My questions are:
(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?
(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?
Thanks in advance,
Ji Woong
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/d664aa41-d272-4e00-96de-70cccaf63af4n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240706/95f9565c/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cp2k.inp
Type: chemical/x-gamess-input
Size: 11659 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240706/95f9565c/attachment.inp>
More information about the CP2K-user
mailing list