<div dir="ltr">Dear Otani,<div>I have used a proper supercell with lattice constant near to the experimental value. Is the method only for charged surfaces? Also, when I applied pbc with esm, I got good results. The problem is only with bc1 boundary condition.</div><div><br></div><div>Thank you</div><div><br></div><div>Rajesh</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jul 24, 2017 at 9:11 AM, 大谷実 <span dir="ltr"><<a href="mailto:minoru.otani@aist.go.jp" target="_blank">minoru.otani@aist.go.jp</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div id="m_3696643046896593414divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
Hi Rajesh,
<div><br>
</div>
<div>I could not find a problem from your email. But the most probable error is the position of atoms. You need to carefully prepare the input file for the spacial arrangement of the unit cell and atoms. Please look at the original paper (Fig.1 in PRB 73, 115407
 (2006)). As you can see, the center of mass of the slab should locate around the origin (z=0).</div>
<div><br>
</div>
<div>Best regards,</div>
<div>Minoru</div>
<div><br>
</div>
<div>
<div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------</div>
<div>Minoru Otani</div>
<div>National Institute of Advanced Industrial Science and Technology</div>
<div>Research Centre for Computational Design of Advanced Functional Materials</div>
<div>email : <a href="mailto:minoru.otani@aist.go.jp" target="_blank">minoru.otani@aist.go.jp</a><br>
<div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>---------------</div>
</div>
</div>
<div><br>
<div class="m_3696643046896593414AppleOriginalContents" style="direction:ltr">
<blockquote type="cite"><div><div class="h5">
<div>On Jul 24, 2017, at 11:40, Rajesh <<a href="mailto:creativeidleminds@gmail.com" target="_blank">creativeidleminds@gmail.com</a>> wrote:</div>
<br class="m_3696643046896593414Apple-interchange-newline">
</div></div><div><div><div class="h5">
<div dir="ltr">Dear Users
<div>I am new to qe. I would like to study a 2D material, graphene, and want to introduce vacuum in the direction perpendicular to nanosheet plane. Is esm with bc1 boundary conditions is right? I tried with this option but energy of the sheet is coming
 out to be positive which is not possible. Can anyone guide me?</div>
<div><br>
</div>
<div><br>
</div>
<div>Thank you.</div>
<div>Rajesh</div>
</div></div></div>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a></div>
</blockquote>
</div>
<br>
</div>
</div>
</div>

<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>