<div dir="ltr"><div><div><div>Dear MP,<br><br></div>  My personal opinion: It depends on how large an error you can afford: If you need to get the effects close to the Dirac cone "correctly" (well, you might need a super cell of 60x60 or so), you cannot avoid the diagonalisation because the electronic structure in the ground state has a degeneracy exactly at the Fermi energy, and with OT you can only have one of those states occupied, breaking the symmetry. The density of states, however, is very small around the point K (or, the "Dirac point"), so a small error in the states around the K point does not cause large errors in many of the observables. If you have something that breaks the symmetry, for example an adsorbate, a small gap is usually opened and in principle OT might work. My general answer would be "it depends", one needs to test case by case, but in the particular case of graphene I would expect that the error made is not very large. Better to test 1) if you manage to converge the electronic structure and 2) how large the difference with OT and diagonalisation is in some cases.<br><br></div>    Greetings,<br><br></div>       apsi<br></div><div class="gmail_extra"><br><div class="gmail_quote">2014-11-18 12:37 GMT+01:00 mpol <span dir="ltr"><<a href="mailto:polynsk...@gmail.com" target="_blank">polynsk...@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2k developers and users,<div><br></div><div>Diagonalization instead of OT is suggested for KS wavefunction optimization in case of metals. But what about systems with small gap (semiconductors, transition metal clusters without ligands, etc.) or gapless systems (like e.g. graphene), can we use OT in these cases? Let's consider graphene, if one choses the supercell in a way that Dirac point is folded at gamma-point, the "gap" values may be small (~0.01 eV), even below the accuracy of currently available DFT methods. Can one use FULL_ALL preconditioner with ENERGY_GAP parameter equal to e.g. 0.0001 Ha in mentioned cases (small or zero gap)?</div><div><br></div><div>With best regards,</div><div>MP</div></div><span class="HOEnZb"><font color="#888888">

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br>  Ari Paavo Seitsonen / <a href="mailto:Ari.P.S...@iki.fi" target="_blank">Ari.P.S...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" target="_blank">http://www.iki.fi/~apsi/</a><br>  Institut für Chemie der Universität Zürich<br>  Tel: +41 44 63 55 44 97  /  Mobile: +41 79 71 90 935</div></div>
</div>