<div dir="ltr">Dear Jing,<div><br></div><div>generally speaking, if you treat the system as a semiconductor (no occupations keywork, no nbnd specified), the default number of bands will be the number of electrons divided by two.<br></div><div>For metallic systems or systems that are treated as they were metallic (occupations = 'smearing') the systems increases this number of bands by a factor (if I remember well by 20%), because</div><div>this is needed in the convergence of the scf cycles when a subset ob bands may cross the Fermi level.</div><div><br></div><div>So, if QE automatically sets the number of bands to 960, it means that you have 960/1.2 = 800 KS states. It is clear that if your system is metallic or if it has a very small gap, the scf cycle will never </div><div>(or quite hardly) converge if you set nbnd = 800 (nor the code will output an error, that is instead issued if you set nbnd < nelec/2, for scf/relax).</div><div><br></div><div>Concerning the time needed to complete the calculation, this depends on a number of factors: number of electrons (high in your case), number of k-points (doubled for spin-polarized calculations</div><div>as yours), starting geometry (far or close to the equilibrium one?), number of needed plane waves (your cutoff is not that small). Moreover, it can significantly depend on the number of processes </div><div>that you allocate on your supercomputer and on the way you use the parallel features of QE.</div><div><br></div><div>So, the answers to your questions:</div><div>- And how should I determine an appropriate value of nbnd before starting any calculation?</div><div>The easiest way (for relax/scf) is to let QE choose the default, bearing in mind that if your system has small gaps, is metallic or is spin polarized you need occupations = 'smearing'.</div><div><br></div><div>- Is this a normal timescale for a system with >100 atoms? <br></div><div>It can be normal, but there is not sufficient information to establish whether or not you're efficiently exploring the parallelization features of QE.</div><div><br></div><div>Giovanni</div><div><br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno sab 3 giu 2023 alle ore 22:18 Jing Lian Ng <<a href="mailto:jinglian@utexas.edu">jinglian@utexas.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello all,<div><br></div><div>I am trying to do a vc-relax calculation on y-Fe2O3 with 160 atoms. Currently I am facing convergence issue with the input script below. With nbnd = 800, the calculation usually finishes in ~ 1 hr on a supercomputer, but fails to reach convergence even on the first BFGS step. When I removed the nbnd parameter, QE automatically set the number of Kohn-Sham states to 960 and the calculation was able to proceed. However, the relaxation calculation is taking more than 72 hours and hasn't reached completion (It is still running now, on 15th BFGS cycle). Is this a normal timescale for a system with >100 atoms? And how should I determine an appropriate value of nbnd before starting any calculation?</div><div><br></div><div>############ Input Script ####################</div><div>&SYSTEM<br> ibrav = 14<br> A = 8.48998 <br> B = 8.48998 <br> C = 25.46994 <br> cosAB = 0 <br> cosAC = 0 <br> cosBC = 0 <br> nat = 160<br> ntyp = 3<br> ecutwfc = 65<br> ecutrho = 780<br> nbnd = 800<br> occupations = 'smearing' <br> degauss = 0.001<br> smearing = 'mv' <br> nspin = 2 <br> starting_magnetization(1) = 0.5<br> starting_magnetization(2) = -0.5<br> starting_magnetization(3) = 0<br> /<br> &ELECTRONS<br> conv_thr = 0.00001<br> mixing_mode = 'plain' <br> mixing_beta = 0.3<br> electron_maxstep = 200<br> /<br> &IONS<br> /<br> &CELL<br> !cell_dofree = 'volume' <br> /<br>ATOMIC_SPECIES<br> FeO 55.84500 Fe.pbesol-spn-kjpaw_psl.0.2.1.UPF <br> FeT 55.84500 Fe.pbesol-spn-kjpaw_psl.0.2.1.UPF <br> O 15.99900 O.pbesol-n-kjpaw_psl.0.1.UPF <br></div><div>##########################################################</div><div><br></div><div><div>Thanks,</div><div>Jing Lian Ng</div><div>Graduate student at CHE, University of Texas at Austin</div></div><div><br></div></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>