<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>