<div dir="ltr"><div dir="ltr">Dear QE Community,<br></div>Could someone clarify if the following is correct: assume_isolated='2D' implies that the system (2D slab) should be centered around z=0 and the total thickness of the system should be less than half of the Z dimension of the simulation cell. In other words, if Z_dim is the height of the simulation cell, all atoms should be within [-Z_dim/4, Z_dim/4]. <div>Thanks in advance,</div><div><br></div><div>Michal Krompiec</div><div>Merck KGaA & University of Southampton<br></div><div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 12 Apr 2019 at 11:07, Thomas Brumme <<a href="mailto:thomas.brumme@uni-leipzig.de">thomas.brumme@uni-leipzig.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div bgcolor="#FFFFFF">
Dear Julien,<br>
<br>
it's actually described indirectly in the paper - in all figures you
can see that the<br>
system is centered a z = 0. Otherwise, the cutoff will cut off some
of the real<br>
Coulomb interactions. Probably this should be added to the
description, but this<br>
should be done by the original author. If he doesn't answer here, I
will contact<br>
him.<br>
<br>
The way you center should be more like (zh+zl)=0 - I don't
understand why<br>
you want to divide by 2 (zero divided by 2?!). That the slab isn't
centered should<br>
not be problematic. From your input I also see that the vacuum
region should be<br>
large enough but maybe you can try increasing it a bit more - maybe
increase the<br>
cell dimension in z to 45 or 50 Angstrom?<br>
<br>
Cheerio<br>
<br>
Thomas<br>
<br><br></div></blockquote></div></div></div></div>