<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><br>&nbsp;&nbsp;&nbsp; Gentlemen,<br><br>My structures do not see to be converging
with increased ecutwfc. &nbsp; I'm setting ecutrho to ecutwfc * 12 and using
US PP. &nbsp;&nbsp; below is a table of ecut values vs energies vs difference
from previous steps followed by my input file followed by part of the output
file.&nbsp;&nbsp;&nbsp; It appears to begin to converge and 120 ecut but then jumps at
130 and then more at 140. &nbsp;&nbsp; I understood that convergence should be
pretty well achieved by 30 ecut for US PP.&nbsp; The only thing I change
between runs are the two ecut parameters.&nbsp; Any ideas as to what I'm
doing wrong or what I should be expecting?&nbsp; <br><br>A follow up
question in that I'm trying to do vc-relax based on this scf and the
forces decrease until a value of about 0.0003 and then oscillate around
this value.&nbsp;&nbsp;&nbsp; could that be related to the ecut convergence problem
I'm having?<br><br>thanks,<br><br>Tim Mason,&nbsp;&nbsp; University of Missouri St. Louis Department of Physics and astronomy<br>------------------------------------------------------------------------------------<br>ecutwfc&nbsp; total energy&nbsp;&nbsp;&nbsp; difference from previous <br><br>100&nbsp;&nbsp;&nbsp; -48.21607621<br>105&nbsp;&nbsp;&nbsp; -48.21651331&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0004371<br>110&nbsp;&nbsp;&nbsp; -48.21672787&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00021456<br>120&nbsp;&nbsp;&nbsp; -48.21686396&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00013609<br>130&nbsp;&nbsp;&nbsp; -48.2172448&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00038084<br>140&nbsp;&nbsp;&nbsp; -48.21819748&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00095268<br><br><br>----------------------------------------------------------------------<br>&amp;CONTROL<br>calculation = "scf",<br>prefix = 'n4_scf.out',<br>tprnfor = .TRUE.<br>pseudo_dir =
 '/home/thmmqc/pwscf/psps/US_GGA',<br>outdir='/home/thmmqc/tmp'<br>etot_conv_thr = 1.d-8<br>/<br><br>&amp;SYSTEM<br>ibrav =
 0, nat = 8, ntyp = 3,<br>ecutwfc = 100<br>ecutrho = 1200<br>celldm(1) = 1.889725989<br><br>/<br><br>&amp;ELECTRONS<br>diagonalization = 'cg',<br>mixing_mode = 'plain',<br>mixing_beta = 0.7,<br>conv_thr = 1.0D-8<br>/<br><br><br>ATOMIC_SPECIES<br>Li 6.941000 Li.pbe-n-van.UPF<br>N 14.006740 N.pbe-van_ak.UPF<br>H 1.007940 H.pbe-van_ak.UPF<br><br>CELL_PARAMETERS<br>10.52658029 -1.08441956 -0.28831668 <br>-2.96581684 6.61198966 -0.14661465 <br>1.15738394 0.29429438 3.51384172 <br><br>ATOMIC_POSITIONS {crystal}<br>Li 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_0">403906747000</span> 0.720101910000 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_1">473406799000</span> <br>Li 0.228004145000 0.310445951000 0.053578743000 <br>&nbsp;N 0.230067712000 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;"
 class="yshortcuts" id="lw_1225985948_2">303154666000</span> 0.552891390000 <br>&nbsp;N 0.403264134000 0.730145360000 -0.024840278000 <br>&nbsp;H 0.185705821000 0.137557379000 0.553855070000 <br>&nbsp;H 0.132830606000 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_3">318409519000</span> 0.624406047000 <br>&nbsp;H 0.486061391000 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_4">685122024000</span> -0.095482152000 <br>&nbsp;H 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_5">467101514000</span> 0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_6">897467040000</span> -0.<span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts"
 id="lw_1225985948_7">051109646000</span>
 <br><br>K_POINTS automatic<br>2 2 2 0 0 0 <br><br>---------------------------------------------------------------------------------------------------------------------<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Program PWSCF&nbsp;&nbsp;&nbsp;&nbsp; v.4.0.2&nbsp; starts ...<br>&nbsp;&nbsp;&nbsp;&nbsp; Today is&nbsp; 6Nov2008 at&nbsp; 8:30: 1 <br><br>&nbsp;&nbsp;&nbsp;&nbsp; Parallel version (MPI)<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Number of processors in use:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br><br>&nbsp;&nbsp;&nbsp;&nbsp; For Norm-Conserving or Ultrasoft (Vanderbilt) Pseudopotentials or PAW<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Current dimensions of program pwscf are:<br>&nbsp;&nbsp;&nbsp;&nbsp; Max number of different atomic species (ntypx) = 10<br>&nbsp;&nbsp;&nbsp;&nbsp; Max number of k-points (npk) =&nbsp; 40000<br>&nbsp;&nbsp;&nbsp;&nbsp; Max <span class="yshortcuts" id="lw_1225985948_8">angular momentum</span> in pseudopotentials (lmaxx) =&nbsp;
 3<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Iterative solution of the
 eigenvalue problem<br>&nbsp;&nbsp;&nbsp;&nbsp; Too few procs for parallel algorithm<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; we need at least 4 procs per pool<br><br>&nbsp;&nbsp;&nbsp;&nbsp; a serial algorithm will be used<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; Planes per process (thick) : nr3 = 96 npp =&nbsp; 96 ncplane =51840<br>&nbsp;&nbsp;&nbsp;&nbsp; Planes per process (smooth): nr3s= 54 npps=&nbsp; 54 ncplanes=17280<br>&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp; Proc/&nbsp; planes cols&nbsp;&nbsp;&nbsp;&nbsp; G&nbsp;&nbsp;&nbsp; planes cols&nbsp;&nbsp;&nbsp; G&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; columns&nbsp; G<br>&nbsp;&nbsp;&nbsp;&nbsp; Pool&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (dense grid)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (smooth grid)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (wavefct grid)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 96&nbsp; 32843&nbsp; 1954809&nbsp;&nbsp; 54&nbsp; 10579&nbsp;&nbsp; 357021&nbsp;&nbsp; 2893&nbsp;&nbsp;&nbsp;
 51177<br>&nbsp;<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; bravais-lattice index&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp;&nbsp;&nbsp;&nbsp; lattice parameter (a_0)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.8897&nbsp; a.u.<br>&nbsp;&nbsp;&nbsp;&nbsp; unit-cell volume&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; 1595.0644 (a.u.)^3<br>&nbsp;&nbsp;&nbsp;&nbsp; number of atoms/cell&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8<br>&nbsp;&nbsp;&nbsp;&nbsp; number of atomic types&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3<br>&nbsp;&nbsp;&nbsp;&nbsp; <span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_9">number of electrons</span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 16.00<br>&nbsp;&nbsp;&nbsp;&nbsp; number of Kohn-Sham
 states=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8<br>&nbsp;&nbsp;&nbsp;&nbsp; kinetic-energy cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 140.0000&nbsp; <span style="border-bottom: 1px dashed rgb(0, 102, 204); cursor: pointer;" class="yshortcuts" id="lw_1225985948_10">Ry</span><br>&nbsp;&nbsp;&nbsp;&nbsp; <span class="yshortcuts" id="lw_1225985948_11">charge density</span> cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; 1740.0000&nbsp; Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; convergence threshold&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0E-08<br>&nbsp;&nbsp;&nbsp;&nbsp; mixing beta&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.7000<br>&nbsp;&nbsp;&nbsp;&nbsp; number of iterations used =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp; plain&nbsp;&nbsp;&nbsp;&nbsp; mixing<br>&nbsp;&nbsp;&nbsp;&nbsp;
 Exchange-correlation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; SLA&nbsp; PW&nbsp;&nbsp; PBE&nbsp; PBE (1434)<br><br>&nbsp;&nbsp;&nbsp;&nbsp; celldm(1)=&nbsp;&nbsp; 1.889726&nbsp;
 celldm(2)=&nbsp;&nbsp; 0.000000&nbsp; celldm(3)=&nbsp;&nbsp; 0.000000<br>&nbsp;&nbsp;&nbsp;&nbsp; celldm(4)=&nbsp;&nbsp; 0.000000&nbsp; celldm(5)=&nbsp;&nbsp; 0.000000&nbsp; celldm(6)=&nbsp;&nbsp; 0.000000<br><br>&nbsp;&nbsp;&nbsp;&nbsp; crystal axes: (cart. coord. in units of a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(1) = ( 10.526580 -1.084420 -0.288317 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = ( -2.965817&nbsp; 6.611990 -0.146615 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = (&nbsp; 1.157384&nbsp; 0.294294&nbsp; 3.513842 )&nbsp; <br><br>&nbsp;&nbsp;&nbsp;&nbsp; reciprocal axes: (cart. coord. in units 2 pi/a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(1) = (&nbsp; 0.098478&nbsp; 0.043373 -0.036069 )&nbsp;
 <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(2) = (&nbsp; 0.015762&nbsp; 0.157902 -0.018417 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(3) = (&nbsp; 0.008738&nbsp; 0.010147&nbsp; 0.280861 )&nbsp; <br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; PseudoPot. # 1 for Li read from file Li.pbe-n-van.UPF<br>&nbsp;&nbsp;&nbsp;&nbsp; Pseudo is Ultrasoft + core correction, Zval =&nbsp; 1.0<br>&nbsp;&nbsp;&nbsp;&nbsp; Generated by new atomic code, or converted to UPF format<br>&nbsp;&nbsp;&nbsp;&nbsp; Using radial grid of&nbsp; 751 points,&nbsp; 2 beta functions with: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(1) =&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(2) =&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; Q(r) pseudized with&nbsp; 8 coefficients,&nbsp; rinner
 =&nbsp;&nbsp;&nbsp; 1.000&nbsp;&nbsp; 1.000&nbsp;&nbsp; 1.000<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; PseudoPot. # 2 for N&nbsp; read from file N.pbe-van_ak.UPF<br>&nbsp;&nbsp;&nbsp;&nbsp; Pseudo is Ultrasoft, Zval =&nbsp; 5.0<br>&nbsp;&nbsp;&nbsp;&nbsp; Generated by new atomic code, or converted to UPF format<br>&nbsp;&nbsp;&nbsp;&nbsp; Using radial grid of&nbsp; 729 points,&nbsp; 4 beta functions with: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(1) =&nbsp;&nbsp; 0<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(2) =&nbsp;&nbsp; 0<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(3) =&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(4) =&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; Q(r) pseudized with&nbsp; 8 coefficients,&nbsp; rinner =&nbsp;&nbsp;&nbsp;
 0.800&nbsp;&nbsp; 0.800&nbsp;&nbsp; 0.800<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; PseudoPot. # 3 for H&nbsp; read from file H.pbe-van_ak.UPF<br>&nbsp;&nbsp;&nbsp;&nbsp; Pseudo is Ultrasoft, Zval =&nbsp; 1.0<br>&nbsp;&nbsp;&nbsp;&nbsp; Generated by new atomic code, or converted to UPF format<br>&nbsp;&nbsp;&nbsp;&nbsp; Using radial grid of&nbsp; 615 points,&nbsp; 1 beta functions with: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; l(1) =&nbsp;&nbsp; 0<br>&nbsp;&nbsp;&nbsp;&nbsp; Q(r) pseudized with&nbsp; 8 coefficients,&nbsp; rinner =&nbsp;&nbsp;&nbsp; 0.800<br><br>&nbsp;&nbsp;&nbsp;&nbsp; atomic species&nbsp;&nbsp; valence&nbsp;&nbsp;&nbsp; mass&nbsp;&nbsp;&nbsp;&nbsp; pseudopotential<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Li&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.00&nbsp;&nbsp;&nbsp;&nbsp; 6.94100&nbsp;&nbsp;&nbsp;&nbsp; Li(
 1.00)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.00&nbsp;&nbsp;&nbsp; 14.00674&nbsp;&nbsp;&nbsp;&nbsp; N ( 1.00)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.00&nbsp;&nbsp;&nbsp;&nbsp; 1.00794&nbsp;&nbsp;&nbsp;&nbsp; H ( 1.00)<br><br>&nbsp;&nbsp;&nbsp;&nbsp; No symmetry!<br><br>&nbsp;&nbsp; Cartesian axes<br><br>&nbsp;&nbsp;&nbsp;&nbsp; site n.&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; positions (a_0 units)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Li&nbsp; tau(&nbsp; 1) = (&nbsp;&nbsp; 2.6639799&nbsp;&nbsp; 4.4626230&nbsp;&nbsp; 1.4414460&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Li&nbsp; tau(&nbsp; 2) = (&nbsp;&nbsp; 1.5413893&nbsp;&nbsp; 1.8211812&nbsp;&nbsp; 0.0770139&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp; tau(&nbsp; 3) = (&nbsp;&nbsp; 2.1626326&nbsp;&nbsp; 1.9176784&nbsp;&nbsp; 1.8319936&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp; tau(&nbsp; 4) = (&nbsp;&nbsp; 2.0507651&nbsp;&nbsp; 4.3830957&nbsp; -0.3106026&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; tau(&nbsp; 5) = (&nbsp;&nbsp; 2.1879002&nbsp;&nbsp; 0.8711414&nbsp;&nbsp; 1.8724490&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; tau(&nbsp; 6) =
 (&nbsp;&nbsp; 1.1765853&nbsp;&nbsp; 2.1450355&nbsp;&nbsp; 2.1090832&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; tau(&nbsp; 7) = (&nbsp;&nbsp; 2.9741083&nbsp;&nbsp; 3.9748254&nbsp; -0.5760977&nbsp; )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; tau(&nbsp; 8) = (&nbsp;&nbsp; 2.1961052&nbsp;&nbsp; 5.4124675&nbsp; -0.4458462&nbsp; )<br><br>&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp; Self-consistent Calculation<br><br>&nbsp;&nbsp;&nbsp;&nbsp; iteration #&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; ecut=&nbsp;&nbsp; 140.00 Ry&nbsp;&nbsp;&nbsp;&nbsp; beta=0.70<br>&nbsp;&nbsp;&nbsp;&nbsp; CG style diagonalization<br>&nbsp;&nbsp;&nbsp;&nbsp; ethr =&nbsp; 1.00E-02,&nbsp; avg #
 of iterations =&nbsp; 3.3<br><br>&nbsp;&nbsp;&nbsp;&nbsp; negative rho (up, down):&nbsp; 0.143E-04 0.000E+00<br><br>&nbsp;&nbsp;&nbsp; <br></td></tr></table>