
<div dir="ltr"><div dir="ltr"><div dir="ltr">I have been running my optimization calculation on my WN slabÂ for 6 hours. From iteration 1 to iteration 500 my energy value has been changing to the 10^-2 or less. Do you think my energy value is going to converge but needs more time or is there a problem (and if so, how can I fix it)? This is the output file I got:<div><br></div><div><div>Message from routine get_command_line:</div><div>Â  Â  Â unexpected argument # 2 Â  Â  :-i</div><div><br></div><div>Â  Â  Â Program PWSCF v.6.1 (svn rev. 13369) starts on 30Jun2020 at Â 1:30:19Â </div><div><br></div><div>Â  Â  Â This program is part of the open-source Quantum ESPRESSO suite</div><div>Â  Â  Â for quantum simulation of materials; please cite</div><div>Â  Â  Â  Â  Â "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);</div><div>Â  Â  Â  Â  Â  URLÂ <a href="http://www.quantum-espresso.org/" target="_blank">http://www.quantum-espresso.org</a>",Â </div><div>Â  Â  Â in publications or presentations arising from this work. More details at</div><div>Â  Â  Â <a href="http://www.quantum-espresso.org/quote" target="_blank">http://www.quantum-espresso.org/quote</a></div><div><br></div><div>Â  Â  Â Parallel version (MPI), running on Â  192 processors</div><div>Â  Â  Â R & G space division: Â proc/nbgrp/npool/nimage = Â  Â  192</div><div>Â  Â  Â Waiting for input...</div><div>Â  Â  Â Reading input from standard input</div><div>Warning: card &CELL ignored</div><div>Warning: card / ignored</div><div><br></div><div>Â  Â  Â Current dimensions of program PWSCF are:</div><div>Â  Â  Â Max number of different atomic species (ntypx) = 10</div><div>Â  Â  Â Max number of k-points (npk) = Â 40000</div><div>Â  Â  Â Max angular momentum in pseudopotentials (lmaxx) = Â 3</div><div>Â  Â  Â  Â  Â  Â  Â  Â file N.pbe-n-kjpaw_psl.1.0.0.UPF: wavefunction(s) Â 2S renormalized</div><div>Â  Â  Â  Â  Â  Â  Â  Â file W.pbe-spn-kjpaw_psl.1.0.0.UPF: wavefunction(s) Â 6S 5P 5D renormalized</div><div><br></div><div>Â  Â  Â Subspace diagonalization in iterative solution of the eigenvalue problem:</div><div>Â  Â  Â one sub-group per band group will be used</div><div>Â  Â  Â scalapack distributed-memory algorithm (size of sub-group: Â 9* Â 9 procs)</div><div><br></div><div>Â  Â  Â Found symmetry operation: I + ( -0.5000 Â 0.0000 Â 0.0000)</div><div>Â  Â  Â This is a supercell, fractional translations are disabled</div><div>Â </div><div>Â  Â  Â Parallelization info</div><div>Â  Â  Â --------------------</div><div>Â  Â  Â sticks: Â  dense Â smooth Â  Â  PW Â  Â  G-vecs: Â  Â dense Â  smooth Â  Â  Â PW</div><div>Â  Â  Â Min Â  Â  Â  Â  Â 52 Â  Â  Â 22 Â  Â  Â 5 Â  Â  Â  Â  Â  Â  Â  Â 11415 Â  Â  3150 Â  Â  417</div><div>Â  Â  Â Max Â  Â  Â  Â  Â 53 Â  Â  Â 23 Â  Â  Â 6 Â  Â  Â  Â  Â  Â  Â  Â 11445 Â  Â  3168 Â  Â  449</div><div>Â  Â  Â Sum Â  Â  Â  10093 Â  Â 4291 Â  1151 Â  Â  Â  Â  Â  Â  Â 2194297 Â  606291 Â  84835</div><div>Â </div><div><br></div><div><br></div><div>Â  Â  Â bravais-lattice index Â  Â  = Â  Â  Â  Â  Â  14</div><div>Â  Â  Â lattice parameter (alat) Â = Â  Â  Â 15.7361 Â a.u.</div><div>Â  Â  Â unit-cell volume Â  Â  Â  Â  Â = Â  12543.4129 (a.u.)^3</div><div>Â  Â  Â number of atoms/cell Â  Â  Â = Â  Â  Â  Â  Â  48</div><div>Â  Â  Â number of atomic types Â  Â = Â  Â  Â  Â  Â  Â 2</div><div>Â  Â  Â number of electrons Â  Â  Â  = Â  Â  Â  456.00</div><div>Â  Â  Â number of Kohn-Sham states= Â  Â  Â  Â  Â 490</div><div>Â  Â  Â kinetic-energy cutoff Â  Â  = Â  Â  Â 50.3902 Â Ry</div><div>Â  Â  Â charge density cutoff Â  Â  = Â  Â  475.2210 Â Ry</div><div>Â  Â  Â convergence threshold Â  Â  = Â  Â  Â 1.0E-06</div><div>Â  Â  Â mixing beta Â  Â  Â  Â  Â  Â  Â  = Â  Â  Â  0.1581</div><div>Â  Â  Â number of iterations used = Â  Â  Â  Â  Â  Â 8 Â local-TF Â mixing</div><div>Â  Â  Â Exchange-correlation Â  Â  Â = Â SLA Â PW Â  PBX Â PBC ( 1 Â 4 Â 3 Â 4 0 0)</div><div>Â  Â  Â nstep Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  = Â  Â  Â  Â  Â 200</div><div><br></div><div><br></div><div>Â  Â  Â celldm(1)= Â 15.736052 Â celldm(2)= Â  1.079226 Â celldm(3)= Â  3.035453</div><div>Â  Â  Â celldm(4)= Â -0.185547 Â celldm(5)= Â  0.000000 Â celldm(6)= Â  0.000000</div><div><br></div><div>Â  Â  Â crystal axes: (cart. coord. in units of alat)</div><div>Â  Â  Â  Â  Â  Â  Â  Â a(1) = ( Â  1.000000 Â  0.000000 Â  0.000000 ) Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â a(2) = ( Â  0.000000 Â  1.079226 Â  0.000000 ) Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â a(3) = ( Â  0.000000 Â -0.563219 Â  2.982743 ) Â </div><div><br></div><div>Â  Â  Â reciprocal axes: (cart. coord. in units 2 pi/alat)</div><div>Â  Â  Â  Â  Â  Â  Â  Â b(1) = ( Â 1.000000 -0.000000 -0.000000 ) Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â b(2) = ( Â 0.000000 Â 0.926590 Â 0.174964 ) Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â b(3) = ( Â 0.000000 Â 0.000000 Â 0.335262 ) Â </div><div><br></div><div><br></div><div>Â  Â  Â PseudoPot. # 1 for Â N read from file:</div><div>Â  Â  Â /home/coralee/projects/def-jkopysci/coralee/.pseudopot/N.pbe-n-kjpaw_psl.1.0.0.UPF</div><div>Â  Â  Â MD5 check sum: e0e4e94a9c4025c5b51bd7d8793849bd</div><div>Â  Â  Â Pseudo is Projector augmented-wave + core cor, Zval = Â 5.0</div><div>Â  Â  Â Generated using "atomic" code by A. Dal Corso Â v.5.1.2</div><div>Â  Â  Â Shape of augmentation charge: PSQ</div><div>Â  Â  Â Using radial grid of 1085 points, Â 4 beta functions with:Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(1) = Â  0</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(2) = Â  0</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(3) = Â  1</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(4) = Â  1</div><div>Â  Â  Â Q(r) pseudized with 0 coefficientsÂ </div><div><br></div><div><br></div><div>Â  Â  Â PseudoPot. # 2 for Â W read from file:</div><div>Â  Â  Â /home/coralee/projects/def-jkopysci/coralee/.pseudopot/W.pbe-spn-kjpaw_psl.1.0.0.UPF</div><div>Â  Â  Â MD5 check sum: f3acacb803c85a3663896168a67a7ce2</div><div>Â  Â  Â Pseudo is Projector augmented-wave + core cor, Zval = 14.0</div><div>Â  Â  Â Generated using "atomic" code by A. Dal Corso Â v.5.1.2</div><div>Â  Â  Â Shape of augmentation charge: PSQ</div><div>Â  Â  Â Using radial grid of 1273 points, Â 6 beta functions with:Â </div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(1) = Â  0</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(2) = Â  0</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(3) = Â  1</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(4) = Â  1</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(5) = Â  2</div><div>Â  Â  Â  Â  Â  Â  Â  Â  l(6) = Â  2</div><div>Â  Â  Â Q(r) pseudized with 0 coefficientsÂ </div><div><br></div><div><br></div><div>Â  Â  Â atomic species Â  valence Â  Â mass Â  Â  pseudopotential</div><div>Â  Â  Â  Â  N Â  Â  Â  Â  Â  Â  Â 5.00 Â  Â 14.00674 Â  Â  Â N( 1.00)</div><div>Â  Â  Â  Â  W Â  Â  Â  Â  Â  Â  14.00 Â  183.84000 Â  Â  Â W( 1.00)</div><div><br></div><div>Â  Â  Â Starting magnetic structureÂ </div><div>Â  Â  Â atomic species Â  magnetization</div><div>Â  Â  Â  Â  N Â  Â  Â  Â  Â  Â 0.200</div><div>Â  Â  Â  Â  W Â  Â  Â  Â  Â  Â 0.200</div><div><br></div><div>Â  Â  Â No symmetry found</div><div><br></div><div><br></div><div><br></div><div>Â  Â Cartesian axes</div><div><br></div><div>Â  Â  Â site n. Â  Â  atom Â  Â  Â  Â  Â  Â  Â  Â  Â positions (alat units)</div><div>Â  Â  Â  Â  Â 1 Â  Â  Â  Â  Â  N Â  tau( Â  1) = ( Â  0.3852312 Â  0.3209930 Â  1.7818541 Â )</div><div>Â  Â  Â  Â  Â 2 Â  Â  Â  Â  Â  N Â  tau( Â  2) = ( Â  0.8852312 Â  0.3209930 Â  1.7818541 Â )</div><div>Â  Â  Â  Â  Â 3 Â  Â  Â  Â  Â  N Â  tau( Â  3) = ( Â  0.3852312 Â -0.2186200 Â  1.7818541 Â )</div><div>Â  Â  Â  Â  Â 4 Â  Â  Â  Â  Â  N Â  tau( Â  4) = ( Â  0.8852312 Â -0.2186200 Â  1.7818541 Â )</div><div>Â  Â  Â  Â  Â 5 Â  Â  Â  Â  Â  W Â  tau( Â  5) = ( Â  0.3852312 Â  0.5863722 Â  1.7711498 Â )</div><div>Â  Â  Â  Â  Â 6 Â  Â  Â  Â  Â  W Â  tau( Â  6) = ( Â  0.8852312 Â  0.5863722 Â  1.7711498 Â )</div><div>Â  Â  Â  Â  Â 7 Â  Â  Â  Â  Â  W Â  tau( Â  7) = ( Â  0.3852312 Â  0.0467592 Â  1.7711498 Â )</div><div>Â  Â  Â  Â  Â 8 Â  Â  Â  Â  Â  W Â  tau( Â  8) = ( Â  0.8852312 Â  0.0467592 Â  1.7711498 Â )</div><div>Â  Â  Â  Â  Â 9 Â  Â  Â  Â  Â  N Â  tau( Â  9) = ( Â  0.1352313 Â  0.6443508 Â  1.7374686 Â )</div><div>Â  Â  Â  Â  10 Â  Â  Â  Â  Â  N Â  tau( Â 10) = ( Â  0.6352312 Â  0.6443508 Â  1.7374686 Â )</div><div>Â  Â  Â  Â  11 Â  Â  Â  Â  Â  N Â  tau( Â 11) = ( Â  0.1352313 Â  0.1047377 Â  1.7374686 Â )</div><div>Â  Â  Â  Â  12 Â  Â  Â  Â  Â  N Â  tau( Â 12) = ( Â  0.6352312 Â  0.1047377 Â  1.7374686 Â )</div><div>Â  Â  Â  Â  13 Â  Â  Â  Â  Â  W Â  tau( Â 13) = ( Â  0.1352313 Â  0.3648474 Â  1.7311320 Â )</div><div>Â  Â  Â  Â  14 Â  Â  Â  Â  Â  W Â  tau( Â 14) = ( Â  0.6352312 Â  0.3648474 Â  1.7311320 Â )</div><div>Â  Â  Â  Â  15 Â  Â  Â  Â  Â  W Â  tau( Â 15) = ( Â  0.1352313 Â -0.1747656 Â  1.7311320 Â )</div><div>Â  Â  Â  Â  16 Â  Â  Â  Â  Â  W Â  tau( Â 16) = ( Â  0.6352312 Â -0.1747656 Â  1.7311320 Â )</div><div>Â  Â  Â  Â  17 Â  Â  Â  Â  Â  N Â  tau( Â 17) = ( Â  0.3852312 Â  0.6408628 Â  1.5167330 Â )</div><div>Â  Â  Â  Â  18 Â  Â  Â  Â  Â  N Â  tau( Â 18) = ( Â  0.8852312 Â  0.6408628 Â  1.5167330 Â )</div><div>Â  Â  Â  Â  19 Â  Â  Â  Â  Â  N Â  tau( Â 19) = ( Â  0.3852312 Â  0.1012498 Â  1.5167330 Â )</div><div>Â  Â  Â  Â  20 Â  Â  Â  Â  Â  N Â  tau( Â 20) = ( Â  0.8852312 Â  0.1012498 Â  1.5167330 Â )</div><div>Â  Â  Â  Â  21 Â  Â  Â  Â  Â  W Â  tau( Â 21) = ( Â  0.3852312 Â  0.3666288 Â  1.5060287 Â )</div><div>Â  Â  Â  Â  22 Â  Â  Â  Â  Â  W Â  tau( Â 22) = ( Â  0.8852312 Â  0.3666288 Â  1.5060287 Â )</div><div>Â  Â  Â  Â  23 Â  Â  Â  Â  Â  W Â  tau( Â 23) = ( Â  0.3852312 Â -0.1729842 Â  1.5060287 Â )</div><div>Â  Â  Â  Â  24 Â  Â  Â  Â  Â  W Â  tau( Â 24) = ( Â  0.8852312 Â -0.1729842 Â  1.5060287 Â )</div><div>Â  Â  Â  Â  25 Â  Â  Â  Â  Â  N Â  tau( Â 25) = ( Â  0.1352313 Â  0.4246075 Â  1.4723475 Â )</div><div>Â  Â  Â  Â  26 Â  Â  Â  Â  Â  N Â  tau( Â 26) = ( Â  0.6352312 Â  0.4246075 Â  1.4723475 Â )</div><div>Â  Â  Â  Â  27 Â  Â  Â  Â  Â  N Â  tau( Â 27) = ( Â  0.1352313 Â -0.1150055 Â  1.4723475 Â )</div><div>Â  Â  Â  Â  28 Â  Â  Â  Â  Â  N Â  tau( Â 28) = ( Â  0.6352312 Â -0.1150055 Â  1.4723475 Â )</div><div>Â  Â  Â  Â  29 Â  Â  Â  Â  Â  W Â  tau( Â 29) = ( Â  0.1352313 Â  0.6847172 Â  1.4660107 Â )</div><div>Â  Â  Â  Â  30 Â  Â  Â  Â  Â  W Â  tau( Â 30) = ( Â  0.6352312 Â  0.6847172 Â  1.4660107 Â )</div><div>Â  Â  Â  Â  31 Â  Â  Â  Â  Â  W Â  tau( Â 31) = ( Â  0.1352313 Â  0.1451041 Â  1.4660107 Â )</div><div>Â  Â  Â  Â  32 Â  Â  Â  Â  Â  W Â  tau( Â 32) = ( Â  0.6352312 Â  0.1451041 Â  1.4660107 Â )</div><div>Â  Â  Â  Â  33 Â  Â  Â  Â  Â  N Â  tau( Â 33) = ( Â  0.3852312 Â  0.4211194 Â  1.2516118 Â )</div><div>Â  Â  Â  Â  34 Â  Â  Â  Â  Â  N Â  tau( Â 34) = ( Â  0.8852312 Â  0.4211194 Â  1.2516118 Â )</div><div>Â  Â  Â  Â  35 Â  Â  Â  Â  Â  N Â  tau( Â 35) = ( Â  0.3852312 Â -0.1184936 Â  1.2516118 Â )</div><div>Â  Â  Â  Â  36 Â  Â  Â  Â  Â  N Â  tau( Â 36) = ( Â  0.8852312 Â -0.1184936 Â  1.2516118 Â )</div><div>Â  Â  Â  Â  37 Â  Â  Â  Â  Â  W Â  tau( Â 37) = ( Â  0.3852312 Â  0.6864986 Â  1.2409076 Â )</div><div>Â  Â  Â  Â  38 Â  Â  Â  Â  Â  W Â  tau( Â 38) = ( Â  0.8852312 Â  0.6864986 Â  1.2409076 Â )</div><div>Â  Â  Â  Â  39 Â  Â  Â  Â  Â  W Â  tau( Â 39) = ( Â  0.3852312 Â  0.1468856 Â  1.2409076 Â )</div><div>Â  Â  Â  Â  40 Â  Â  Â  Â  Â  W Â  tau( Â 40) = ( Â  0.8852312 Â  0.1468856 Â  1.2409076 Â )</div><div>Â  Â  Â  Â  41 Â  Â  Â  Â  Â  N Â  tau( Â 41) = ( Â  0.1352313 Â  0.7444773 Â  1.2072262 Â )</div><div>Â  Â  Â  Â  42 Â  Â  Â  Â  Â  N Â  tau( Â 42) = ( Â  0.6352312 Â  0.7444773 Â  1.2072262 Â )</div><div>Â  Â  Â  Â  43 Â  Â  Â  Â  Â  N Â  tau( Â 43) = ( Â  0.1352313 Â  0.2048643 Â  1.2072262 Â )</div><div>Â  Â  Â  Â  44 Â  Â  Â  Â  Â  N Â  tau( Â 44) = ( Â  0.6352312 Â  0.2048643 Â  1.2072262 Â )</div><div>Â  Â  Â  Â  45 Â  Â  Â  Â  Â  W Â  tau( Â 45) = ( Â  0.1352313 Â  0.4649739 Â  1.2008896 Â )</div><div>Â  Â  Â  Â  46 Â  Â  Â  Â  Â  W Â  tau( Â 46) = ( Â  0.6352312 Â  0.4649739 Â  1.2008896 Â )</div><div>Â  Â  Â  Â  47 Â  Â  Â  Â  Â  W Â  tau( Â 47) = ( Â  0.1352313 Â -0.0746391 Â  1.2008896 Â )</div><div>Â  Â  Â  Â  48 Â  Â  Â  Â  Â  W Â  tau( Â 48) = ( Â  0.6352312 Â -0.0746391 Â  1.2008896 Â )</div><div><br></div><div>Â  Â  Â number of k points= Â  Â  4 Â gaussian smearing, width (Ry)= Â 0.0200</div><div>Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â cart. coord. in units 2pi/alat</div><div>Â  Â  Â  Â  k( Â  Â 1) = ( Â  0.0000000 Â  0.0000000 Â  0.0000000), wk = Â  0.2500000</div><div>Â  Â  Â  Â  k( Â  Â 2) = ( Â  0.0000000 Â -0.4632949 Â -0.0874821), wk = Â  0.2500000</div><div>Â  Â  Â  Â  k( Â  Â 3) = ( Â -0.5000000 Â  0.0000000 Â  0.0000000), wk = Â  0.2500000</div><div>Â  Â  Â  Â  k( Â  Â 4) = ( Â -0.5000000 Â -0.4632949 Â -0.0874821), wk = Â  0.2500000</div><div><br></div><div>Â  Â  Â Dense Â grid: Â 2194297 G-vectors Â  Â  FFT dimensions: ( 120, 120, 360)</div><div><br></div><div>Â  Â  Â Smooth grid: Â  606291 G-vectors Â  Â  FFT dimensions: ( Â 72, Â 80, 216)</div><div><br></div><div>Â  Â  Â Estimated max dynamical RAM per process > Â  Â  Â 75.34MB</div><div><br></div><div>Â  Â  Â Estimated total allocated dynamical RAM > Â  14464.36MB</div><div>Â  Â  Â Generating pointlists ...</div><div>Â  Â  Â new r_m : Â  0.1068 (alat units) Â 1.6801 (a.u.) for type Â  Â 1</div><div>Â  Â  Â new r_m : Â  0.1068 (alat units) Â 1.6801 (a.u.) for type Â  Â 2</div><div><br></div><div>Â  Â  Â Initial potential from superposition of free atoms</div><div>Â  Â  Â Check: negative starting charge=(component1): Â  -0.001367</div><div>Â  Â  Â Check: negative starting charge=(component2): Â  -0.000912</div><div><br></div><div>Â  Â  Â starting charge Â 455.98646, renormalised to Â 456.00000</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 1.367E-03 9.115E-04</div><div>Â  Â  Â Starting wfc are Â 408 randomized atomic wfcs + Â  82 random wfc</div><div>Â  Â  Â Checking if some PAW data can be deallocated...Â </div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â  Â  11.8 secs</div><div><br></div><div>Â  Â  Â per-process dynamical memory: Â  111.1 Mb</div><div><br></div><div>Â  Â  Â Self-consistent Calculation</div><div><br></div><div>Â  Â  Â iteration # Â 1 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 6.4</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 1.815E-04 7.397E-05</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â  Â  91.7 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18922.62699425 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18923.50326373 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â 18.35022984 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â 17.52 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â 21.04 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration # Â 2 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 4.02E-03, Â avg # of iterations = Â 1.0</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 2.600E-01 2.782E-01</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â  Â 124.8 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18920.61522629 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.77912017 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â 11.59887509 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â 16.77 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â 20.11 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration # Â 3 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 2.54E-03, Â avg # of iterations = Â 1.4</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 1.835E-01 2.099E-01</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â  Â 157.9 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -18922.16475654 Ry</div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.19366580 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  1.17725139 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â 12.59 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â 13.91 Bohr mag/cell</div></div><div><br></div><div>.....</div><div><br></div><div><div style="font-family:arial,sans-serif;color:rgb(0,0,0)"><div>iteration #519 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 2.0</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.807E-02 8.736E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 25760.8 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -18922.33525065 Ry</div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33934098 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00032491 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  2.80 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.59 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #520 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 5.0</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.734E-02 8.750E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 25836.4 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -18922.33845071 Ry</div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33845795 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00029607 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  3.49 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.60 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #521 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 3.0</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.659E-02 8.729E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 25883.8 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18922.33852878 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33847917 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00032948 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  3.55 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.62 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #522 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 2.5</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.530E-02 8.190E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 25922.4 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18922.33802807 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33853778 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00039280 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  3.58 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.64 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #523 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 4.6</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.528E-02 8.146E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 25982.9 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18922.33916064 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33912989 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00026667 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  3.19 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.58 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #524 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div>Â  Â  Â ethr = Â 8.38E-10, Â avg # of iterations = Â 2.1</div><div><br></div><div>Â  Â  Â negative rho (up, down): Â 4.556E-02 8.074E-02</div><div><br></div><div>Â  Â  Â total cpu time spent up to now is Â  Â 26025.7 secs</div><div><br></div><div>Â  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-18922.33919903 Ry</span></div><div>Â  Â  Â Harris-Foulkes estimate Â  = Â -18922.33916652 Ry</div><div>Â  Â  Â estimated scf accuracy Â  Â < Â  Â  Â  0.00025649 Ry</div><div><br></div><div>Â  Â  Â total magnetization Â  Â  Â  = Â  Â  3.16 Bohr mag/cell</div><div>Â  Â  Â absolute magnetization Â  Â = Â  Â  5.58 Bohr mag/cell</div><div><br></div><div>Â  Â  Â iteration #525 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.16</div><div>Â  Â  Â Davidson diagonalization with overlap</div><div><br></div></div></div><div><br></div></div></div></div>