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">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">https://groups.google.com/d/msgid/cp2k/035d2574-d796-4112-bdd7-f8e4cad2ab01n%40googlegroups.com</a>.<br />