<div dir="auto"><div>Continuing the thread for clarification; please correct my errors:<div dir="auto"><br></div><div dir="auto">If I remember right, lmaxv & lmaxe change how the LM expansion is done inside the spheres, so need to be larger with higher L. The CPU scales something like a square, and they are CPU expensive.</div><div dir="auto"><br></div><div dir="auto">GMAX is a plane wave expansion, so (like lapw0) would need to be larger for small RMT such as H or odd cases. Again something like a square dependence, although I am not sure how CPU expensive except for a surface.</div><div dir="auto"><br></div><div dir="auto">The CPU time will scale something like the product of the outer k-mesh and the inner (for which a reduced mesh can be used). Standard accuracy issue.</div><div dir="auto"><br></div><div dir="auto">NBand in case.inhf is I think linear in CPU; I am not exactly sure why it needs to be larger than the number of occupied bands (forgotten).</div><div dir="auto"><br></div><div dir="auto">RKMAX is standard, with a 2-3 power dependence of CPU time.</div><div dir="auto"><br></div><div dir="auto">N.B., I think there is also a dependence upon how different the straight lapw1 & hf are, for instance the band gap in lapw1. This matters for how large the number of states used in lapw1 has to be.</div><div dir="auto"><br></div><div dir="auto"><br><br><div data-smartmail="gmail_signature" dir="auto">_____<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 15, 2020, 08:38 Tran, Fabien <<a href="mailto:fabien.tran@tuwien.ac.at">fabien.tran@tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">That's difficult since it may depend on the system, RMT, property, and required accuracy, but roughly (concerning  the parameters in case.inhf):<br>
<br>
-For geometry: none of the parameters really need to be increased compared to default. In particular, nband can be set to minimum required (number of occupied bands plus one). Maybe gmax (from case.inhf) should be tested if RMT are chosen clearly smaller than default. The option "-nonself" could also be used for non-magnetic systems.<br>
<br>
-For cohesive energy: preferably test nband, gmax, lmaxe and lmaxv. Especially the latter two for systems with an open  f-shell.<br>
<br>
-For electronic structure: test lmaxe and lmaxv for systems with an open  f-shell. For the band gap of non-magnetic systems, the option "-diaghf" (<a href="https://urldefense.com/v3/__https://doi.org/10.1016/j.physleta.2012.01.022__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rjxGqJUjQ$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__https://doi.org/10.1016/j.physleta.2012.01.022__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rjxGqJUjQ$</a> ) can be used to reduce the computational time by two orders of magnitude.<br>
<br>
-For electron density (e.g., for EFG): Test carefully all parameters.<br>
<br>
-k-mesh: In general, the convergence of properties with HF/hybrids is slower than with LDA/GGA/MGGA (fortunately only slightly slower with screened hybrids like HSE06).<br>
<br>
-RKmax: the convergence is the same as for LDA/GGA/MGGA.<br>
<br>
-Computational time: strongly influenced by k-mesh, RKmax, nband, lmaxe and lmaxv. Less by gmax (from case.inhf).<br>
<br>
<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Laurence Marks <<a href="mailto:laurence.marks@gmail.com" target="_blank" rel="noreferrer">laurence.marks@gmail.com</a>><br>
Sent: Saturday, August 15, 2020 2:40 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] Wien post from <a href="mailto:pascal.boulet@univ-amu.fr" target="_blank" rel="noreferrer">pascal.boulet@univ-amu.fr</a> (Errors with DGEMM and hybrid calculations)<br>
  <br>
It would be good to have a brief summary of what matters with hf:<br>
a) For speed<br>
b) For accuracy<br>
<br>
It is probably somewhere in the docu, but another cite for the list would be useful (to me as well as others).<br>
_____<br>
Professor Laurence Marks<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br>
<a href="http://www.numis.northwestern.edu" rel="noreferrer noreferrer" target="_blank">http://www.numis.northwestern.edu</a><br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rgXzLvd8g$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rgXzLvd8g$</a> <br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rjETTfF-Q$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!AN3JSEcq-L36m1Ickxc1xjixv8KB9lE8fbCvyAvgSj7ySnH9QhVXJB8tDGGm-rjETTfF-Q$</a> <br>
</blockquote></div></div></div>