<div dir="ltr">Dear Thomas,<div><br></div><div>thank you for your explanation. I am now curious why this does not seem to effect my VASP calculation but I guess I have to ascribed it to different implementations of either the dipole correction or how the background is treated... AS VASP only allows to charge "cubic cells" I guess they implement the m-p scheme:</div><div><span style="color:rgb(0,0,0);font-family:sans-serif;font-size:12.7px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:sans-serif;font-size:12.7px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:sans-serif;font-size:12.7px">I will try assume_isolated='2D' and 'mp' and see if I can reproduce the VASP results!</span></div><div><br></div><div>Thank you again for your help!</div><div><br></div><div>Chris </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 28, 2018 at 7:04 PM, Dr. Thomas Brumme <span dir="ltr"><<a href="mailto:thomas.brumme@uni-leipzig.de" target="_blank">thomas.brumme@uni-leipzig.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Chris,<br>
<br>
The potential shows the typical quadratic dependence on z since you're calculating a charged system - there is a homogeneous background charge since the 3D pbc system is assumed to be neutral. This has nothing to do with the dipole correction. Depending on what you want to do next it might be useful to set the flag assume_isolated='2D'<br>
<br>
Regards<br>
<br>
Thomas<br>
<br>
Zitat von Christoph Wolf <wolf.christoph@qns.science>:<div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear all,<br>
<br>
I am still observing something strange in my slab + dipole correction<br>
calculation that I do not fully understand.<br>
<br>
When using dipfield+tefield (eopreg and emaxpos well within the vacuum<br>
region) I encounter a "saggy" electrostatic potential (plot_num=11) despite<br>
the sawtooth efield potential (plot_num=12) looking as usual. Maybe someone<br>
can give it a look and confirm if this is due to the excess charge in the<br>
system (this does not happen when running the same system in VASP)?<br>
<br>
I attach input and the plot of the potential for 2 and 4 layers of vacuum<br>
(more vacuum does seem to improve the situation).<br>
<br>
Thank you in advance for your time and assistance!<br>
<br>
Chris<br>
<br>
PS: I am using a QE version with dipole bug fix:<br>
<br>
     Program PWSCF v.6.2.2 starts on 25Jun2018 at 10:39:24<br>
<br>
--<br>
Postdoctoral Researcher<br>
Center for Quantum Nanoscience, Institute for Basic Science<br>
Ewha Womans University, Seoul, South Korea<br>
</blockquote>
<br></div></div>
--<br>
Dr. rer. nat. Thomas Brumme<br>
Wilhelm-Ostwald-Institute for Physical and Theoretical Chemistry<br>
Leipzig University<br>
Phillipp-Rosenthal-Strasse 31<br>
04103 Leipzig<br>
<br>
Tel:  +49 (0)341 97 36456<br>
<br>
email: <a href="mailto:thomas.brumme@uni-leipzig.de" target="_blank">thomas.brumme@uni-leipzig.de</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div>
</div>