<div dir="ltr"><div><div><div><div><div><div><div><div>Dear all<br></div>For this simple case of extended in the XY plane system the procedure of 2d basis reduction<br></div>controlled by parameters ewind and epsproj should not be very efficient, in fact it is better<br>
</div>if you use the full 2d basis, which is in this case n2d=271... You can do that specifying<br></div>very big value for ewind, say ewind=400.<br></div>I would also reduce delgep parameter to say delgep=5.d-10. It needs to be set to bigger<br>
</div>values if the code crashes...<br></div>good luck,<br></div>Alex<br><div><div><div><br><div> <br></div></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-07-18 16:39 GMT+02:00 Sclauzero  Gabriele <span dir="ltr"><<a href="mailto:gabriele.sclauzero@mat.ethz.ch" target="_blank">gabriele.sclauzero@mat.ethz.ch</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Vladislav,<br>
<br>
    I mostly agree with what Nicki said here below. I also believe that in your case the problem is the ewind value, which is too small and can lead to some weirdnesses in the 2D-basis set reduction. You should try at least 4.0 and check convergence as suggested by Nicki.<br>

In my experience, the choice of epsproj depends a lot on the system. Even values around 1.d-5 could assure a decent convergence, unless the basis set reduction went bananas (e.g., because of a too small ewind).<br>
Also, if you push epsproj very low, but you keep ewind very small, then you could also get into trouble.<br>
In normal calculations, I never found the need to increase the nz1 from its default.<br>
<br>
Finally, a word about the crystal structure. Are you sure that’s what you want? The neighbors of the Co atoms seem to have different bond distances. Maybe you need to choose another value for celldm(3)?<br>
<br>
HTH<br>
<br>
<br>
GS<br>
<br>
<br>
<br>
<br>
<br>
Il giorno 18/lug/2014, alle ore 12:53, nicvok <<a href="mailto:nicvok@freenet.de">nicvok@freenet.de</a>> ha scritto:<br>
<div><div class="h5"><br>
> The parameters epsproj and ewind control the amount of basis functions<br>
> used in the underlying complex<br>
> bandstructure (CBS) calculation in your case 121; cf. Choi&Ihm (<a href="http://journals.aps.org/prb/pdf/10.1103/PhysRevB.59.2267" target="_blank">http://journals.aps.org/prb/pdf/10.1103/PhysRevB.59.2267</a><br>

> ) and<br>
> DalCorso et al. (<a href="http://journals.aps.org/prb/abstract/10.1103/PhysRevB.74.045429" target="_blank">http://journals.aps.org/prb/abstract/10.1103/PhysRevB.74.045429</a><br>
> ). So one should take a<br>
> look at the CBS before and check for reasonableness and convergence.<br>
> Then adapt step by step epsproj and ewind (even more nz1) and check<br>
> for convergence of the CBS.<br>
> In your case ewind=2 might be okay, nevertheless I would slightly<br>
> increase this (3-5). Did you ever tried<br>
> some meaningful epsproj< 1.d-7, which should be a good value? The<br>
> number of integral points with 22 seems to<br>
> be okay. That it worked for other systems could be by accident. Again,<br>
> as always, one finds the optimal<br>
> value of the parameters, by checking for convergence of the underlying<br>
> physics.<br>
><br>
> bests Nicki<br>
><br>
> -------------------------------------------------------------<br>
> Nicki Frank Hinsche, Dr. rer. nat.<br>
> Institute of physics - Theoretical physics,<br>
> Martin-Luther-University Halle-Wittenberg,<br>
> Von-Seckendorff-Platz 1, Room 1.07<br>
> D-06120 Halle/Saale, Germany<br>
> Tel.: ++49 345 5525460<br>
> -------------------------------------------------------------<br>
> _______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
</div></div>Dr. Gabriele Sclauzero<br>
Materials Theory - ETHZ<br>
ETH Zurich, HIT G 43.2<br>
Wolfgang-Pauli-Str. 27<br>
8093 Zürich, Switzerland<br>
<br>
Phone <a href="tel:%2B41%2044%20633%2094%2010" value="+41446339410">+41 44 633 94 10</a><br>
Fax <a href="tel:%2B41%2044%20633%2014%2059" value="+41446331459">+41 44 633 14 59</a><br>
<a href="mailto:gabriele.sclauzero@mat.ethz.ch">gabriele.sclauzero@mat.ethz.ch</a><br>
<a href="http://www.theory.mat.ethz.ch" target="_blank">www.theory.mat.ethz.ch</a><br>
<div class="HOEnZb"><div class="h5"><br>
<br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div>