<div dir="ltr">On Fri, Sep 9, 2016 at 9:30 AM, Daniel Stoeffler <span dir="ltr"><<a href="mailto:daniel.stoeffler@ipcms.unistra.fr" target="_blank">daniel.stoeffler@ipcms.unistra.fr</a>></span> wrote:<br><div class="gmail_extra"><br><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">Surprisingly, the computation time is huge</blockquote><div><br></div><div>no sir: "unsurprisingly" the computation time is huge. The computational complexity of a calculation with hybrid functionals goes as (Nk*Nb)^2*Npw log(Npw) or so, where Nk=number of k-points, Nb=number of occupied bands, Npw=number of plane waves; with GGA functionals, as Nk*Nb*Npw log(Npw) or <br> Nk*Nb*Npw^2<br><br><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">could you try to explain me how to reduce the computation time<br clear="all"></blockquote></div></div><br></div><div class="gmail_extra">first of all you may try to reduce one of the sums over k-points using options nq1,nq2,nq3, as explained in the example for hybrid. For a metal, not sure it works well, though. Since you are using UltraSoft PseudoPotentials, try option "tqr=.true." (real-space treatment of augmentation functions): it should give a significant speedup. Another significant improvement is contained in the forthcoming 6.0 version of QE, but it doesn't work yet for USPP<br><br></div><div class="gmail_extra">Paolo<br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>