<div dir="ltr">Hi,<div><br></div><div>About the vacuum level, when a charged system is calculated, the charges in the slab system interacts with the replica in the "non-periodic" direction, creating a electric field itself and making the Hartree potential non plane in the vacuum zone... this also happens sometimes for neutral slab systems with dipoles into the slab or on the surfaces which do not cancel each other. In VASP there are dipoles and electric field corrections to accomplish this task without this effect maintaining the 3D periodicity... I guess you are using those in your VASP calculations. In CP2K AFAIK there is dipole corrections, but do not remember if there is electric field correction, but here there is the possibility to turn off the periodicity of the cell in the "non-periodical" direction of the slab, maybe this is the way you need to go.</div><div><br></div><div>Regards - Lucas</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El vie, 7 abr 2023 a las 18:38, ling chen (<<a href="mailto:ling.chen@adelaide.edu.au">ling.chen@adelaide.edu.au</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Experts,<div><br></div><div>I am currently working on calculating the work function of a
system and have encountered a few challenges that I would appreciate your
insights on.<br></div><div>

<p>For the zero-charge system, my current method involves
calculating the Hartree energy, taking the value in the vacuum region (usually
a horizontal line, as shown in the attached figure), and subtracting the Fermi
energy obtained from the .out file. This approach yields results that agree
well with those obtained using VASP. </p>

<p> However, when dealing with charged systems, such as
charge = -2, the Hartree energy in the vacuum region becomes a diagonal line
(as shown in the attached figure). Using the average value and subtracting the
Fermi energy results in a significant discrepancy compared to the VASP calculations.
It has been suggested that I use the SCCS implicit water model to stabilize the
vacuum energy level, which does indeed restore the horizontal line in the
Hartree energy. However, the calculated work function deviates even further
from the VASP results.</p>

<p>I am unsure about the cause of this discrepancy and would
like to inquire if you have any recommendations for calculating the work
function in charged systems. Additionally, I have noticed that the work
function calculation is significantly influenced by the CUTOFF and REL_CUTOFF
parameters in MGRID. What values would you suggest to be appropriate for these
parameters? </p>

<p> I would be grateful for your expertise in addressing
these concerns. Thank you in advance for your time and assistance. </p><p>Kind Regards,</p><p><img alt="Screenshot 2023-04-08 073641.png" width="460px" height="271px" src="cid:8548158c-287c-446d-bfb5-0bfc9da8be41"><img alt="Screenshot 2023-04-08 073716.png" width="466px" height="253px" src="cid:a48cb721-da69-4e14-975c-65fbef7de58d">Ling Chen</p></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" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/035d2574-d796-4112-bdd7-f8e4cad2ab01n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/035d2574-d796-4112-bdd7-f8e4cad2ab01n%40googlegroups.com</a>.<br>
</blockquote></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/CAOFT4PJGPUdHgkqtmo_kGnJAWqDuSVVwCr1P1UYop4BNtHe0NQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAOFT4PJGPUdHgkqtmo_kGnJAWqDuSVVwCr1P1UYop4BNtHe0NQ%40mail.gmail.com</a>.<br />