<div dir="ltr"><div>Dear Dorothea,</div><div>To break a complex problem down: Should I consider it a bug, that I get non-constant potentials in the plot of the hartree energies, when using MT or implicit BCs and that potential in the metal does not match EXT_POTENTIAL? (see plot below - also posted in the original question) Or is there a good reason for this behavior?</div><div><br></div><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="margin-left:40px">red:      IMPLICIT poissn solver with Neumann BC at z=0 and z=60 and homogeneous Dirichlet BC at z=30</div><div style="margin-left:40px">green:  MT (Martyna-Tuckerman) poisson solver (shifted by 0.4V)</div><div style="margin-left:40px">blue:    normal Ewald summation (shifted by 0.4V)<br></div><div><br></div><p style="text-align:center;clear:both"><img style="margin-left:1em;margin-right:1em" src="https://groups.google.com/group/cp2k/attach/13361e99fe472/tmp3.png?part=0.3&view=1&authuser=0" alt="tmp3.png" width="400" height="337"></p><p style="text-align:center;clear:both"><br></p><br></div></blockquote></div></blockquote><div>Thank you for your help!</div><div>Best regards,</div><div>Katharina</div><div><br></div></div>