vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex                        
  
 MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign.
 com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1
 6685
 
 This VASP executable licensed from Materials Design, Inc.
 
 executed on                        Lin64 date 2024.08.19  09:35:48
 running on    2 total cores
 distrk:  each k-point on    2 cores,    1 groups
 distr:  one band on NCORE=   1 cores,    2 groups


--------------------------------------------------------------------------------------------------------


 INCAR:
   SYSTEM = No title
   PREC = Normal
   ENCUT = 400.000
   IBRION = -1
   NSW = 0
   ISIF = 2
   NELMIN = 2
   EDIFF = 1.0e-05
   EDIFFG = -0.02
   VOSKOWN = 1
   NBLOCK = 1
   NWRITE = 1
   NELM = 95
   LHFCALC = .TRUE.
   HFSCREEN = 0.2
   PRECFOCK = Normal
   ALGO = Damped
   TIME = 0.4
   LMAXFOCK = 4
   NKREDX = 1
   NKREDY = 1
   NKREDZ = 1
   ISPIN = 1
   INIWAV = 1
   ISTART = 0
   ICHARG = 2
   LWAVE = .TRUE.
   LCHARG = .FALSE.
   ADDGRID = .FALSE.
   ISMEAR = 1
   SIGMA = 0.2
   LREAL = Auto
   LSCALAPACK = .FALSE.
   RWIGS = 1.11 0.73 0.77 0.32 0.75
   NPAR = 2

 POTCAR:    PAW_PBE Si 05Jan2001                  
 POTCAR:    PAW_PBE O 08Apr2002                   
 POTCAR:    PAW_PBE C 08Apr2002                   
 POTCAR:    PAW_PBE H 15Jun2001                   
 POTCAR:    PAW_PBE N 08Apr2002                   
 POTCAR:    PAW_PBE Si 05Jan2001                  
  local pseudopotential read in
  partial core-charges read in
  partial kinetic energy density read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in
 
   number of l-projection  operators is LMAX  =           4
   number of lm-projection operators is LMMAX =           8
 
 POTCAR:    PAW_PBE O 08Apr2002                   
  local pseudopotential read in
  partial core-charges read in
  partial kinetic energy density read in
  kinetic energy density of atom read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in
 
   number of l-projection  operators is LMAX  =           4
   number of lm-projection operators is LMMAX =           8
 
 POTCAR:    PAW_PBE C 08Apr2002                   
  local pseudopotential read in
  partial core-charges read in
  partial kinetic energy density read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in
 
   number of l-projection  operators is LMAX  =           4
   number of lm-projection operators is LMMAX =           8
 
 POTCAR:    PAW_PBE H 15Jun2001                   
  local pseudopotential read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in
 
   number of l-projection  operators is LMAX  =           3
   number of lm-projection operators is LMMAX =           5
 
 POTCAR:    PAW_PBE N 08Apr2002                   
  local pseudopotential read in
  partial core-charges read in
  partial kinetic energy density read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in
 
   number of l-projection  operators is LMAX  =           4
   number of lm-projection operators is LMMAX =           8
 
 Optimization of the real space projectors (new method)

 maximal supplied QI-value         = 19.84
 optimisation between [QCUT,QGAM] = [ 10.12, 20.44] = [ 28.68,116.96] Ry 
 Optimized for a Real-space Cutoff    1.23 Angstroem

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0      7    10.119   159.560    0.56E-04    0.22E-03    0.45E-07
   0      7    10.119   115.863    0.56E-04    0.21E-03    0.45E-07
   1      7    10.119    88.339    0.34E-03    0.49E-03    0.11E-06
   1      7    10.119    48.592    0.33E-03    0.48E-03    0.11E-06
 Optimization of the real space projectors (new method)

 maximal supplied QI-value         = 24.76
 optimisation between [QCUT,QGAM] = [ 10.15, 20.30] = [ 28.85,115.39] Ry 
 Optimized for a Real-space Cutoff    1.38 Angstroem

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0      8    10.150    20.381    0.22E-03    0.32E-03    0.29E-06
   0      8    10.150    15.268    0.23E-03    0.35E-03    0.30E-06
   1      8    10.150     5.964    0.46E-03    0.53E-03    0.21E-06
   1      8    10.150     5.382    0.38E-03    0.45E-03    0.19E-06
 Optimization of the real space projectors (new method)

 maximal supplied QI-value         = 25.13
 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry 
 Optimized for a Real-space Cutoff    1.30 Angstroem

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0      8    10.053   115.676    0.49E-03    0.72E-03    0.18E-06
   0      8    10.053    87.132    0.49E-03    0.71E-03    0.18E-06
   1      7    10.053     4.429    0.32E-03    0.31E-03    0.18E-06
   1      7    10.053     2.733    0.23E-03    0.19E-03    0.20E-06
 Optimization of the real space projectors (new method)

 maximal supplied QI-value         = 34.20
 optimisation between [QCUT,QGAM] = [  9.92, 20.18] = [ 27.55,114.04] Ry 
 Optimized for a Real-space Cutoff    1.26 Angstroem

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0      8     9.919    19.460    0.50E-03    0.23E-03    0.29E-06
   0      8     9.919    12.209    0.48E-03    0.23E-03    0.28E-06
   1      7     9.919     4.655    0.17E-03    0.75E-03    0.30E-06
 Optimization of the real space projectors (new method)

 maximal supplied QI-value         = 25.13
 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry 
 Optimized for a Real-space Cutoff    1.65 Angstroem

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0     10    10.053    79.467    0.76E-04    0.72E-04    0.56E-06
   0     10    10.053    66.151    0.76E-04    0.72E-04    0.55E-06
   1     10    10.053     8.350    0.25E-03    0.92E-03    0.41E-05
   1     10    10.053     5.531    0.27E-03    0.10E-02    0.45E-05
  PAW_PBE Si 05Jan2001                  :
 energy of atom  1       EATOM= -103.0669
 kinetic energy error for atom=    0.0012 (will be added to EATOM!!)
  PAW_PBE O 08Apr2002                   :
 energy of atom  2       EATOM= -432.3788
 kinetic energy error for atom=    0.1156 (will be added to EATOM!!)
  PAW_PBE C 08Apr2002                   :
 energy of atom  3       EATOM= -147.1560
 kinetic energy error for atom=    0.0288 (will be added to EATOM!!)
  PAW_PBE H 15Jun2001                   :
 energy of atom  4       EATOM=  -12.4884
 kinetic energy error for atom=    0.0098 (will be added to EATOM!!)
  PAW_PBE N 08Apr2002                   :
 energy of atom  5       EATOM= -264.5486
 kinetic energy error for atom=    0.0736 (will be added to EATOM!!)
 
 
 POSCAR: No title
  positions in direct lattice
  No initial velocities read in
 exchange correlation table for  LEXCH =        8
   RHO(1)=    0.500       N(1)  =     2000
   RHO(2)=  100.500       N(2)  =     4000
 


--------------------------------------------------------------------------------------------------------


 ion  position               nearest neighbor table
   1  0.429  0.450  0.483-   2 1.61   4 1.85   3 1.94   5 2.12
   2  0.490  0.478  0.548-  12 1.26   1 1.61
   3  0.382  0.339  0.458-   6 1.56   1 1.94
   4  0.472  0.475  0.378-   8 1.07   9 1.09  10 1.23   1 1.85
   5  0.370  0.519  0.578-   7 1.82   1 2.12
   6  0.377  0.255  0.519-  14 1.09  15 1.09  16 1.09   3 1.56
   7  0.326  0.526  0.472-  19 0.98  17 1.16  18 1.21   5 1.82
   8  0.466  0.423  0.331-   4 1.07
   9  0.451  0.540  0.361-   4 1.09
  10  0.532  0.490  0.386-   4 1.23
  11  0.329  0.473  0.646-
  12  0.461  0.443  0.614-   2 1.26
  13  0.423  0.585  0.582-
  14  0.336  0.222  0.485-   6 1.09
  15  0.426  0.225  0.512-   6 1.09
  16  0.364  0.272  0.588-   6 1.09
  17  0.336  0.601  0.459-   7 1.16
  18  0.331  0.497  0.396-   7 1.21
  19  0.287  0.491  0.490-   7 0.98
  20  0.509  0.639  0.444-  22 1.09
  21  0.545  0.647  0.542-  22 1.13
  22  0.497  0.669  0.509-  20 1.09  21 1.13
 
  LATTYP: Found a simple tetragonal cell.
 ALAT       =    15.0000000000
 C/A-ratio  =     1.3333333333
  
  Lattice vectors:
  
 A1 = (   0.0000000000,  15.0000000000,   0.0000000000)
 A2 = (   0.0000000000,   0.0000000000,  15.0000000000)
 A3 = (  20.0000000000,   0.0000000000,   0.0000000000)


Analysis of symmetry for initial positions (statically):
=====================================================================
 Subroutine PRICEL returns:
 Original cell was already a primitive cell.
 

 Routine SETGRP: Setting up the symmetry group for a 
 simple tetragonal supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of 16 trial point group operations.


The static configuration has the point symmetry C_1 .


Analysis of symmetry for dynamics (positions and initial velocities):
=====================================================================
 Subroutine PRICEL returns:
 Original cell was already a primitive cell.
 

 Routine SETGRP: Setting up the symmetry group for a 
 simple tetragonal supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of 16 trial point group operations.


The dynamic configuration has the point symmetry C_1 .


 Subroutine INISYM returns: Found  1 space group operations
 (whereof  1 operations are pure point group operations),
 and found     1 'primitive' translations


----------------------------------------------------------------------------------------

                                     Primitive cell                                     

  volume of cell :    4500.0000

  direct lattice vectors                    reciprocal lattice vectors
    20.000000000  0.000000000  0.000000000     0.050000000  0.000000000  0.000000000
     0.000000000 15.000000000  0.000000000     0.000000000  0.066666667  0.000000000
     0.000000000  0.000000000 15.000000000     0.000000000  0.000000000  0.066666667

  length of vectors
    20.000000000 15.000000000 15.000000000     0.050000000  0.066666667  0.066666667

  position of ions in fractional coordinates (direct lattice)
     0.428635980  0.450034510  0.483076360
     0.489635420  0.477628430  0.547671860
     0.382398530  0.338992240  0.457839930
     0.472230380  0.475416970  0.377616250
     0.369587990  0.518720520  0.577991790
     0.376631060  0.255265060  0.519227650
     0.326373890  0.526029960  0.471547290
     0.465848670  0.422814580  0.330559210
     0.450896820  0.540313280  0.361126020
     0.532239990  0.489772260  0.385842310
     0.328921550  0.473206780  0.646296680
     0.461249380  0.443053530  0.614000640
     0.422748950  0.585221850  0.581639770
     0.335904320  0.221664790  0.484509280
     0.425780560  0.224502040  0.512093300
     0.364313180  0.272124500  0.588191660
     0.336077250  0.601416300  0.458810840
     0.330721230  0.497344750  0.396358830
     0.287318720  0.490943100  0.490332380
     0.508812760  0.639266020  0.444376610
     0.545307940  0.647181110  0.542242190
     0.497471120  0.669079150  0.508812980

  ion indices of the primitive-cell ions
   primitive index   ion index
                 1           1
                 2           2
                 3           3
                 4           4
                 5           5
                 6           6
                 7           7
                 8           8
                 9           9
                10          10
                11          11
                12          12
                13          13
                14          14
                15          15
                16          16
                17          17
                18          18
                19          19
                20          20
                21          21
                22          22

----------------------------------------------------------------------------------------

 
 
 KPOINTS: Automatic mesh                          

Automatic generation of k-mesh.
 Grid dimensions read from file:
 generate k-points for:    1    1    1

 Generating k-lattice:

  Cartesian coordinates                     Fractional coordinates (reciprocal lattice)
     0.050000000  0.000000000  0.000000000     1.000000000  0.000000000  0.000000000
     0.000000000  0.066666667  0.000000000     0.000000000  1.000000000  0.000000000
     0.000000000  0.000000000  0.066666667     0.000000000  0.000000000  1.000000000

  Length of vectors
     0.050000000  0.066666667  0.066666667

  Shift w.r.t. Gamma in fractional coordinates (k-lattice)
     0.000000000  0.000000000  0.000000000

 
 Subroutine IBZKPT returns following result:
 ===========================================
 
 Found      1 irreducible k-points:
 
 Following reciprocal coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
 
 Following cartesian coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
 
 
 Subroutine IBZKPT_HF returns following result:
 ==============================================
 
 Found      1 k-points in 1st BZ
 the following      1 k-points will be used (e.g. in the exchange kernel)
 Following reciprocal coordinates:   # in IRBZ
  0.000000  0.000000  0.000000    1.00000000   1 t-inv F


--------------------------------------------------------------------------------------------------------




 Dimension of arrays:
   k-points           NKPTS =      1   k-points in BZ     NKDIM =      1   number of bands    NBANDS=     38
   number of dos      NEDOS =    301   number of ions     NIONS =     22
   non local maximal  LDIM  =      4   non local SUM 2l+1 LMDIM =      8
   total plane-waves  NPLWV = 627200
   max r-space proj   IRMAX =   2665   max aug-charges    IRDMAX=   4859
   dimension x,y,z NGX =    98 NGY =   80 NGZ =   80
   dimension x,y,z NGXF=   196 NGYF=  160 NGZF=  160
   support grid    NGXF=   196 NGYF=  160 NGZF=  160
   ions per type =               1   2   4  14   1
   NGX,Y,Z   is equivalent  to a cutoff of   8.15,  8.87,  8.87 a.u.
   NGXF,Y,Z  is equivalent  to a cutoff of  16.29, 17.73, 17.73 a.u.

 SYSTEM =  No title                                
 POSCAR =  No title                                

 Startparameter for this run:
   NWRITE =      1    write-flag & timer
   PREC   = normal    normal or accurate (medium, high low for compatibility)
   ISTART =      0    job   : 0-new  1-cont  2-samecut
   ICHARG =      2    charge: 1-file 2-atom 10-const
   ISPIN  =      1    spin polarized calculation?
   LNONCOLLINEAR =      F non collinear calculations
   LSORBIT =      F    spin-orbit coupling
   INIWAV =      1    electr: 0-lowe 1-rand  2-diag
   LASPH  =      F    aspherical Exc in radial PAW
 Electronic Relaxation 1
   ENCUT  =  400.0 eV  29.40 Ry    5.42 a.u.  32.61 24.46 24.46*2*pi/ulx,y,z
   ENINI  =  400.0     initial cutoff
   ENAUG  =  644.9 eV  augmentation charge cutoff
   NELM   =     95;   NELMIN=  2; NELMDL= -5     # of ELM steps 
   EDIFF  = 0.1E-04   stopping-criterion for ELM
   LREAL  =      T    real-space projection
   NLSPLINE    = F    spline interpolate recip. space projectors
   LCOMPAT=      F    compatible to vasp.4.4
   GGA_COMPAT  = T    GGA compatible to vasp.4.4-vasp.4.6
   LMAXPAW     = -100 max onsite density
   LMAXMIX     =    2 max onsite mixed and CHGCAR
   VOSKOWN=      1    Vosko Wilk Nusair interpolation
   ROPT   =   -0.00050  -0.00050  -0.00050  -0.00050
   ROPT   =   -0.00050
 Ionic relaxation
   EDIFFG = -.2E-01   stopping-criterion for IOM
   NSW    =      0    number of steps for IOM
   NBLOCK =      1;   KBLOCK =      1    inner block; outer block 
   IBRION =     -1    ionic relax: 0-MD 1-quasi-New 2-CG
   NFREE  =      0    steps in history (QN), initial steepest desc. (CG)
   ISIF   =      2    stress and relaxation
   IWAVPR =     10    prediction:  0-non 1-charg 2-wave 3-comb
   ISYM   =      3    0-nonsym 1-usesym 2-fastsym
   LCORR  =      T    Harris-Foulkes like correction to forces

   POTIM  = 0.5000    time-step for ionic-motion
   TEIN   =    0.0    initial temperature
   TEBEG  =    0.0;   TEEND  =   0.0 temperature during run
   SMASS  =  -3.00    Nose mass-parameter (am)
   estimated Nose-frequenzy (Omega)   =  0.10E-29 period in steps = 0.13E+47 mass=  -0.914E-26a.u.
   SCALEE = 1.0000    scale energy and forces
   NPACO  =    256;   APACO  = 16.0  distance and # of slots for P.C.
   PSTRESS=    0.0 pullay stress

  Mass of Ions in am
   POMASS =  28.09 16.00 12.01  1.00 14.00
  Ionic Valenz
   ZVAL   =   4.00  6.00  4.00  1.00  5.00
  Atomic Wigner-Seitz radii
   RWIGS  =   1.11  0.73  0.77  0.32  0.75
  virtual crystal weights 
   VCA    =   1.00  1.00  1.00  1.00  1.00
   NELECT =      51.0000    total number of electrons
   NUPDOWN=      -1.0000    fix difference up-down

 DOS related values:
   EMIN   =  10.00;   EMAX   =-10.00  energy-range for DOS
   EFERMI =   0.00
   ISMEAR =     1;   SIGMA  =   0.20  broadening in eV -4-tet -1-fermi 0-gaus

 Electronic relaxation 2 (details)
   IALGO  =     53    algorithm
   LDIAG  =      T    sub-space diagonalisation (order eigenvalues)
   LSUBROT=      F    optimize rotation matrix (better conditioning)
   TURBO    =      0    0=normal 1=particle mesh
   IRESTART =      0    0=no restart 2=restart with 2 vectors
   NREBOOT  =      0    no. of reboots
   NMIN     =      0    reboot dimension
   EREF     =   0.00    reference energy to select bands
   IMIX   =      4    mixing-type and parameters
     AMIX     =   0.40;   BMIX     =  1.00
     AMIX_MAG =   1.60;   BMIX_MAG =  1.00
     AMIN     =   0.10
     WC   =   100.;   INIMIX=   1;  MIXPRE=   1;  MAXMIX= -45

 Intra band minimization:
   WEIMIN = 0.0000     energy-eigenvalue tresh-hold
   EBREAK =  0.66E-07  absolut break condition
   DEPER  =   0.30     relativ break condition  

   TIME   =   0.40     timestep for ELM

  volume/ion in A,a.u.               =     204.55      1380.34
  Fermi-wavevector in a.u.,A,eV,Ry     =   0.367729  0.694906  1.839838  0.135224
  Thomas-Fermi vector in A             =   1.293057
 
 Write flags
   LWAVE        =      T    write WAVECAR
   LDOWNSAMPLE  =      F    k-point downsampling of WAVECAR
   LCHARG       =      F    write CHGCAR
   LVTOT        =      F    write LOCPOT, total local potential
   LVHAR        =      F    write LOCPOT, Hartree potential only
   LELF         =      F    write electronic localiz. function (ELF)
   LORBIT       =      0    0 simple, 1 ext, 2 COOP (PROOUT), +10 PAW based schemes


 Dipole corrections
   LMONO  =      F    monopole corrections only (constant potential shift)
   LDIPOL =      F    correct potential (dipole corrections)
   IDIPOL =      0    1-x, 2-y, 3-z, 4-all directions 
   EPSILON=  1.0000000 bulk dielectric constant

 Exchange correlation treatment:
   GGA     =    --    GGA type
   LEXCH   =     8    internal setting for exchange type
   VOSKOWN=      1    Vosko Wilk Nusair interpolation
   EXXOEP  =     0    0=HF, 1=EXX-LHF (local Hartree Fock) 2=EXX OEP
   LHFCALC =     T    Hartree Fock is set to
   LSYMGRAD=     F    symmetrize gradient (conserves proper symmetry)
   PRECFOCK=normal    Normal, Fast or Accurate (Low or Medium for compatibility)
   LRHFCALC=     F    long range Hartree Fock
   LRSCOR  =     F    long range correlation only (use DFT for short range part)
   LTHOMAS =     F    Thomas Fermi screening in HF
   LMODELHF=     F    short range full HF, long range fraction AEXX
   ENCUT4O =  -1.0   cutoff for four orbital integrals eV
   LMAXFOCK=     4    L truncation for augmentation on plane wave grid
   LMAXFOCKAE=  -1    L truncation for all-electron charge restoration on plane wave grid
   NMAXFOCKAE=   1    number of basis functions for all-electron charge restoration
   LFOCKAEDFT=     F  apply the AE augmentation even for DFT
   NKREDX  =     1    reduce k-point grid by
   NKREDY  =     1    reduce k-point grid by
   NKREDZ  =     1    reduce k-point grid by
   SHIFTRED=     F    shift reduced grid of Gamma
   HFKIDENT=     F    idential grid for each k-point
   ODDONLY =     F    use only odd q-grid points
   EVENONLY=     F    use only even q-grid points
   HFALPHA =   -1.0000 decay constant for conv. correction
   MCALPHA =    0.0000 extent of test-charge in conv. correction in multipole expansion
   AEXX    =    0.2500 exact exchange contribution
   HFSCREEN=    0.2000 screening length (either q_TF or 0.3 A-1)
   HFSCREENC=   0.2000 screening length for correlation (either q_TF or 0.3 A-1)
   HFRCUT  =    0.0000 spherical cutoff for potential kernel
   ALDAX   =    0.7500 LDA exchange part
   AGGAX   =    0.7500 GGA exchange part
   ALDAC   =    1.0000 LDA correlation
   AGGAC   =    1.0000 GGA correlation
   NBANDSGWLOW=     1    first orbital included in HF term
   ENCUTFOCK=  -1.0 apply spherical cutoff to Coloumb kernel
 Linear response parameters
   LEPSILON=     F    determine dielectric tensor
   LRPA    =     F    only Hartree local field effects (RPA)
   LNABLA  =     F    use nabla operator in PAW spheres
   LVEL    =     F    velocity operator in full k-point grid
   LINTERFAST=   F  fast interpolation
   KINTER  =     0    interpolate to denser k-point grid
   CSHIFT  =0.1000    complex shift for real part using Kramers Kronig
   OMEGAMAX=  -1.0    maximum frequency
   DEG_THRESHOLD= 0.2000000E-02 threshold for treating states as degnerate
   RTIME   =   -0.100 relaxation time in fs
  (WPLASMAI=    0.000 imaginary part of plasma frequency in eV, 0.658/RTIME)
   DFIELD  = 0.0000000 0.0000000 0.0000000 field for delta impulse in time
 
 Orbital magnetization related:
   ORBITALMAG=     F  switch on orbital magnetization
   LCHIMAG   =     F  perturbation theory with respect to B field
   DQ        =  0.001000  dq finite difference perturbation B field
   LLRAUG    =     F  two centre corrections for induced B field



--------------------------------------------------------------------------------------------------------


 Static calculation
 charge density and potential will be updated during run
 non-spin polarized calculation
 Conjugate gradient for all bands (Freysoldt, et al. PRB 79, 241103 (2009))
 perform sub-space diagonalisation
    before iterative eigenvector-optimisation
 modified Broyden-mixing scheme, WC =      100.0
 initial mixing is a Kerker type mixing with AMIX =  0.4000 and BMIX =      1.0000
 Hartree-type preconditioning will be used
 using additional bands           12
 real space projection scheme for non local part
 use partial core corrections
 no Harris-corrections to forces 
 use gradient corrections 
 use of overlap-Matrix (Vanderbilt PP)
 Methfessel and Paxton  Order N= 1 SIGMA  =   0.20


--------------------------------------------------------------------------------------------------------


  energy-cutoff  :      400.00
  volume of cell :     4500.00
      direct lattice vectors                 reciprocal lattice vectors
    20.000000000  0.000000000  0.000000000     0.050000000  0.000000000  0.000000000
     0.000000000 15.000000000  0.000000000     0.000000000  0.066666667  0.000000000
     0.000000000  0.000000000 15.000000000     0.000000000  0.000000000  0.066666667

  length of vectors
    20.000000000 15.000000000 15.000000000     0.050000000  0.066666667  0.066666667


 
 k-points in units of 2pi/SCALE and weight: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       1.000
 
 k-points in reciprocal lattice and weights: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       1.000
 
 position of ions in fractional coordinates (direct lattice) 
   0.42863598  0.45003451  0.48307636
   0.48963542  0.47762843  0.54767186
   0.38239853  0.33899224  0.45783993
   0.47223038  0.47541697  0.37761625
   0.36958799  0.51872052  0.57799179
   0.37663106  0.25526506  0.51922765
   0.32637389  0.52602996  0.47154729
   0.46584867  0.42281458  0.33055921
   0.45089682  0.54031328  0.36112602
   0.53223999  0.48977226  0.38584231
   0.32892155  0.47320678  0.64629668
   0.46124938  0.44305353  0.61400064
   0.42274895  0.58522185  0.58163977
   0.33590432  0.22166479  0.48450928
   0.42578056  0.22450204  0.51209330
   0.36431318  0.27212450  0.58819166
   0.33607725  0.60141630  0.45881084
   0.33072123  0.49734475  0.39635883
   0.28731872  0.49094310  0.49033238
   0.50881276  0.63926602  0.44437661
   0.54530794  0.64718111  0.54224219
   0.49747112  0.66907915  0.50881298
 
 position of ions in cartesian coordinates  (Angst):
   8.57271960  6.75051765  7.24614540
   9.79270840  7.16442645  8.21507790
   7.64797060  5.08488360  6.86759895
   9.44460760  7.13125455  5.66424375
   7.39175980  7.78080780  8.66987685
   7.53262120  3.82897590  7.78841475
   6.52747780  7.89044940  7.07320935
   9.31697340  6.34221870  4.95838815
   9.01793640  8.10469920  5.41689030
  10.64479980  7.34658390  5.78763465
   6.57843100  7.09810170  9.69445020
   9.22498760  6.64580295  9.21000960
   8.45497900  8.77832775  8.72459655
   6.71808640  3.32497185  7.26763920
   8.51561120  3.36753060  7.68139950
   7.28626360  4.08186750  8.82287490
   6.72154500  9.02124450  6.88216260
   6.61442460  7.46017125  5.94538245
   5.74637440  7.36414650  7.35498570
  10.17625520  9.58899030  6.66564915
  10.90615880  9.70771665  8.13363285
   9.94942240 10.03618725  7.63219470
 


--------------------------------------------------------------------------------------------------------


 k-point  1 :   0.0000 0.0000 0.0000  plane waves:   81909

 maximum and minimum number of plane-waves per node :     81909    81909

 maximum number of plane-waves:     81909
 maximum index in each direction: 
   IXMAX=   32   IYMAX=   24   IZMAX=   24
   IXMIN=  -32   IYMIN=  -24   IZMIN=  -24

 The following grids will avoid any aliasing or wrap around errors in the Hartre
 e energy
  - symmetry arguments have not been applied
  - exchange correlation energies might require even more grid points
  - we recommend to set PREC=Normal or Accurate and rely on VASP defaults
 WARNING: aliasing errors must be expected set NGX to   140 to avoid them
 WARNING: aliasing errors must be expected set NGY to    98 to avoid them
 WARNING: aliasing errors must be expected set NGZ to    98 to avoid them

 serial   3D FFT for wavefunctions
 parallel 3D FFT for charge:
    minimum data exchange during FFTs selected (reduces bandwidth)

 
 Radii for the augmentation spheres in the non-local exchange
 for species   1 augmentation radius   1.237 (default was   0.989)
       energy cutoff for augmentation   1600.0
 for species   2 augmentation radius   0.902 (default was   0.722)
       energy cutoff for augmentation   1600.0
 for species   3 augmentation radius   0.874 (default was   0.699)
       energy cutoff for augmentation   1600.0
 for species   4 augmentation radius   0.650 (default was   0.520)
       energy cutoff for augmentation   1600.0
 for species   5 augmentation radius   0.886 (default was   0.709)
       energy cutoff for augmentation   1600.0
  SETUP_FOCK is finished

 total amount of memory used by VASP MPI-rank0   423957. kBytes
=======================================================================

   base      :      30000. kBytes
   nonlr-proj:       3414. kBytes
   fftplans  :      87706. kBytes
   grid      :     227924. kBytes
   one-center:         67. kBytes
   HF        :         23. kBytes
   wavefun   :      74823. kBytes
 
 Broyden mixing: mesh for mixing (old mesh)
   NGX = 65   NGY = 49   NGZ = 49
  (NGX  =196   NGY  =160   NGZ  =160)
  gives a total of 156065 points

 initial charge density was supplied:
 charge density of overlapping atoms calculated
 number of electron      51.0000000 magnetization 
 keeping initial charge density in first step


--------------------------------------------------------------------------------------------------------


 Maximum index for non-local projection operator         2542
 Maximum index for augmentation-charges         2251 (set IRDMAX)


--------------------------------------------------------------------------------------------------------


 First call to EWALD:  gamma=   0.107
 Maximum number of real-space cells 2x 3x 3
 Maximum number of reciprocal cells 3x 3x 3



--------------------------------------- Iteration      1(   1)  ---------------------------------------



 eigenvalue-minimisations  :    90
 total energy-change (2. order) : 0.4746116E+03  (-0.7093404E+03)
 number of electron      51.0000000 magnetization 
 augmentation part       51.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4448.50954483
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       147.58523416
  PAW double counting   =      1450.12005463    -1743.20965345
  entropy T*S    EENTRO =        -0.02120880
  eigenvalues    EBANDS =       125.31609624
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =       474.61161430 eV

  energy without entropy =      474.63282310  energy(sigma->0) =      474.61868390
  exchange ACFDT corr.  =         0.00000000  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   2)  ---------------------------------------



 eigenvalue-minimisations  :   152
 total energy-change (2. order) :-0.3308501E+03  (-0.3089064E+03)
 number of electron      51.0000000 magnetization 
 augmentation part       51.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4448.50954483
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       147.58523416
  PAW double counting   =      1450.12005463    -1743.20965345
  entropy T*S    EENTRO =        -0.01889495
  eigenvalues    EBANDS =      -205.53633961
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =       143.76149229 eV

  energy without entropy =      143.78038725  energy(sigma->0) =      143.76779061
  exchange ACFDT corr.  =         0.00000000  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   3)  ---------------------------------------



 eigenvalue-minimisations  :   152
 total energy-change (2. order) :-0.3442338E+03  (-0.3183133E+03)
 number of electron      51.0000000 magnetization 
 augmentation part       51.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4448.50954483
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       147.58523416
  PAW double counting   =      1450.12005463    -1743.20965345
  entropy T*S    EENTRO =        -0.00736882
  eigenvalues    EBANDS =      -549.78171299
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -200.47235496 eV

  energy without entropy =     -200.46498614  energy(sigma->0) =     -200.46989869
  exchange ACFDT corr.  =         0.00000000  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   4)  ---------------------------------------



 eigenvalue-minimisations  :   152
 total energy-change (2. order) :-0.1818614E+03  (-0.1777132E+03)
 number of electron      51.0000000 magnetization 
 augmentation part       51.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4448.50954483
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       147.58523416
  PAW double counting   =      1450.12005463    -1743.20965345
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -731.59409010
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -382.33378220 eV

  energy without entropy =     -382.27736325  energy(sigma->0) =     -382.31497588
  exchange ACFDT corr.  =         0.00000000  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   5)  ---------------------------------------



 eigenvalue-minimisations  :   152
 total energy-change (2. order) : 0.2611881E+03  (-0.6110664E+02)
 number of electron      51.0000060 magnetization 
 augmentation part        2.3583376 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4448.50954483
  -exchange      EXHF   =       110.69855731
  -V(xc)+E(xc)   XCENC  =       113.54986398
  PAW double counting   =      1479.12216159    -1469.61099418
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -849.66994295
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -121.14568170 eV

  energy without entropy =     -121.08926274  energy(sigma->0) =     -121.12687538
  exchange ACFDT corr.  =         0.00000000  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   6)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.3485513E+02  (-0.6644606E+02)
 number of electron      51.0000045 magnetization 
 augmentation part        1.7684765 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4730.72629009
  -exchange      EXHF   =       122.63579738
  -V(xc)+E(xc)   XCENC  =       127.93818385
  PAW double counting   =      4817.25532001    -4810.08425893
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -556.58352080
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =       -86.29055120 eV

  energy without entropy =      -86.23413224  energy(sigma->0) =      -86.27174488
  exchange ACFDT corr.  =        -0.25691140  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   7)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2573542E+02  (-0.1287391E+02)
 number of electron      51.0000044 magnetization 
 augmentation part        1.6533216 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4505.70001297
  -exchange      EXHF   =       110.70467039
  -V(xc)+E(xc)   XCENC  =       117.03568762
  PAW double counting   =      4751.02225193    -4742.29868107
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -786.06410270
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -112.02596942 eV

  energy without entropy =     -111.96955046  energy(sigma->0) =     -112.00716310
  exchange ACFDT corr.  =        -0.30792875  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   8)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.5999784E+01  (-0.8022728E+01)
 number of electron      51.0000047 magnetization 
 augmentation part        1.7269622 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4584.51211935
  -exchange      EXHF   =       112.16184598
  -V(xc)+E(xc)   XCENC  =       118.22837812
  PAW double counting   =      4629.54433037    -4620.37422027
  entropy T*S    EENTRO =        -0.05641881
  eigenvalues    EBANDS =      -716.34818530
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.02575306 eV

  energy without entropy =     -117.96933426  energy(sigma->0) =     -118.00694680
  exchange ACFDT corr.  =        -0.27355837  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(   9)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3493570E+01  (-0.2177860E+01)
 number of electron      51.0000048 magnetization 
 augmentation part        1.7692100 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4668.76547847
  -exchange      EXHF   =       115.66503581
  -V(xc)+E(xc)   XCENC  =       121.36615911
  PAW double counting   =      4444.10293290    -4434.97805618
  entropy T*S    EENTRO =        -0.05641351
  eigenvalues    EBANDS =      -642.18413331
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -121.51932260 eV

  energy without entropy =     -121.46290909  energy(sigma->0) =     -121.50051810
  exchange ACFDT corr.  =        -0.33179798  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  10)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.1202157E+01  (-0.1167060E+01)
 number of electron      51.0000048 magnetization 
 augmentation part        1.7650036 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4671.41159424
  -exchange      EXHF   =       116.22344972
  -V(xc)+E(xc)   XCENC  =       121.83608067
  PAW double counting   =      4275.64636393    -4266.63653877
  entropy T*S    EENTRO =        -0.05638208
  eigenvalues    EBANDS =      -641.65346345
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.72147932 eV

  energy without entropy =     -122.66509724  energy(sigma->0) =     -122.70268529
  exchange ACFDT corr.  =        -0.42452290  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  11)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.5815284E+00  (-0.4302868E+00)
 number of electron      51.0000048 magnetization 
 augmentation part        1.7686746 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4649.24412005
  -exchange      EXHF   =       115.49861142
  -V(xc)+E(xc)   XCENC  =       121.15009586
  PAW double counting   =      4158.19906657    -4149.19868358
  entropy T*S    EENTRO =        -0.05512963
  eigenvalues    EBANDS =      -662.98223219
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.30300771 eV

  energy without entropy =     -123.24787808  energy(sigma->0) =     -123.28463116
  exchange ACFDT corr.  =        -0.50412398  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  12)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2503584E+00  (-0.1719650E+00)
 number of electron      51.0000049 magnetization 
 augmentation part        1.7739971 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4647.07387059
  -exchange      EXHF   =       115.42158604
  -V(xc)+E(xc)   XCENC  =       121.10051301
  PAW double counting   =      4107.45424631    -4098.44694308
  entropy T*S    EENTRO =        -0.04928295
  eigenvalues    EBANDS =      -665.28440450
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.55336608 eV

  energy without entropy =     -123.50408313  energy(sigma->0) =     -123.53693843
  exchange ACFDT corr.  =        -0.48249866  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  13)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.9861188E-01  (-0.7691961E-01)
 number of electron      51.0000049 magnetization 
 augmentation part        1.7721657 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4654.34847757
  -exchange      EXHF   =       115.61428609
  -V(xc)+E(xc)   XCENC  =       121.28305017
  PAW double counting   =      4093.45766387    -4084.45901769
  entropy T*S    EENTRO =        -0.04205700
  eigenvalues    EBANDS =      -658.48083623
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.65197796 eV

  energy without entropy =     -123.60992097  energy(sigma->0) =     -123.63795896
  exchange ACFDT corr.  =        -0.43527175  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  14)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4338561E-01  (-0.5203161E-01)
 number of electron      51.0000048 magnetization 
 augmentation part        1.7706206 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4656.68206005
  -exchange      EXHF   =       115.64344723
  -V(xc)+E(xc)   XCENC  =       121.29620256
  PAW double counting   =      4090.73783269    -4081.74298595
  entropy T*S    EENTRO =        -0.03266655
  eigenvalues    EBANDS =      -656.23637940
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.69536357 eV

  energy without entropy =     -123.66269702  energy(sigma->0) =     -123.68447472
  exchange ACFDT corr.  =        -0.38866732  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  15)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.1626875E-01  (-0.9079161E-01)
 number of electron      51.0000049 magnetization 
 augmentation part        1.7674754 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4652.02948877
  -exchange      EXHF   =       115.58032880
  -V(xc)+E(xc)   XCENC  =       121.24508096
  PAW double counting   =      4091.85886208    -4082.85617269
  entropy T*S    EENTRO =        -0.03401810
  eigenvalues    EBANDS =      -660.80821251
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.71163232 eV

  energy without entropy =     -123.67761423  energy(sigma->0) =     -123.70029296
  exchange ACFDT corr.  =        -0.25425776  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  16)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.1692022E-01  (-0.2009590E+00)
 number of electron      51.0000047 magnetization 
 augmentation part        1.7753956 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4654.62113607
  -exchange      EXHF   =       115.56952333
  -V(xc)+E(xc)   XCENC  =       121.19763640
  PAW double counting   =      4094.33594568    -4085.33424333
  entropy T*S    EENTRO =        -0.03199379
  eigenvalues    EBANDS =      -658.13905636
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.69471210 eV

  energy without entropy =     -123.66271831  energy(sigma->0) =     -123.68404751
  exchange ACFDT corr.  =        -0.33233094  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  17)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2548403E+00  (-0.2955805E+01)
 number of electron      51.0000050 magnetization 
 augmentation part        1.7643698 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4646.37780393
  -exchange      EXHF   =       115.60840570
  -V(xc)+E(xc)   XCENC  =       121.28764676
  PAW double counting   =      4099.67260701    -4090.66119717
  entropy T*S    EENTRO =        -0.05437377
  eigenvalues    EBANDS =      -666.26817273
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.43987182 eV

  energy without entropy =     -123.38549804  energy(sigma->0) =     -123.42174723
  exchange ACFDT corr.  =        -0.17224929  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  18)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.1689353E+00  (-0.7477349E+00)
 number of electron      51.0000047 magnetization 
 augmentation part        1.7764241 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4660.78236098
  -exchange      EXHF   =       115.60019839
  -V(xc)+E(xc)   XCENC  =       121.16914728
  PAW double counting   =      4101.09588729    -4092.10354541
  entropy T*S    EENTRO =        -0.03493244
  eigenvalues    EBANDS =      -651.52652568
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.27093654 eV

  energy without entropy =     -123.23600410  energy(sigma->0) =     -123.25929239
  exchange ACFDT corr.  =        -0.63199166  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  19)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.1153658E+00  (-0.5208551E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7640160 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4643.86365623
  -exchange      EXHF   =       115.62000410
  -V(xc)+E(xc)   XCENC  =       121.30012336
  PAW double counting   =      4106.60908212    -4097.59520844
  entropy T*S    EENTRO =        -0.05529971
  eigenvalues    EBANDS =      -668.52161956
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.15557075 eV

  energy without entropy =     -123.10027104  energy(sigma->0) =     -123.13713751
  exchange ACFDT corr.  =        -0.18858228  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  20)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.1183241E+00  (-0.5781214E+00)
 number of electron      51.0000048 magnetization 
 augmentation part        1.7727824 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.22412603
  -exchange      EXHF   =       115.62057290
  -V(xc)+E(xc)   XCENC  =       121.17743263
  PAW double counting   =      4106.29639693    -4097.30419450
  entropy T*S    EENTRO =        -0.03695482
  eigenvalues    EBANDS =      -651.11531343
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.27389486 eV

  energy without entropy =     -123.23694004  energy(sigma->0) =     -123.26157659
  exchange ACFDT corr.  =        -0.65780673  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  21)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3400524E+00  (-0.7242193E+00)
 number of electron      51.0000050 magnetization 
 augmentation part        1.7679210 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4648.06811510
  -exchange      EXHF   =       115.57621565
  -V(xc)+E(xc)   XCENC  =       121.25166405
  PAW double counting   =      4108.25780436    -4099.24751445
  entropy T*S    EENTRO =        -0.03145025
  eigenvalues    EBANDS =      -664.67768327
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.61394724 eV

  energy without entropy =     -123.58249699  energy(sigma->0) =     -123.60346382
  exchange ACFDT corr.  =        -0.27985701  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  22)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.6343459E-02  (-0.5503706E+00)
 number of electron      51.0000046 magnetization 
 augmentation part        1.7817164 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4657.67276202
  -exchange      EXHF   =       115.63730240
  -V(xc)+E(xc)   XCENC  =       121.21531958
  PAW double counting   =      4108.41107497    -4099.41452750
  entropy T*S    EENTRO =        -0.03708090
  eigenvalues    EBANDS =      -655.08319730
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.60760378 eV

  energy without entropy =     -123.57052289  energy(sigma->0) =     -123.59524348
  exchange ACFDT corr.  =        -0.36773989  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  23)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.1360949E+01  (-0.1572474E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7615449 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4639.22113775
  -exchange      EXHF   =       115.71861006
  -V(xc)+E(xc)   XCENC  =       121.35956096
  PAW double counting   =      4110.69537157    -4101.67825952
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -672.41435540
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.24665463 eV

  energy without entropy =     -122.19023570  energy(sigma->0) =     -122.22784832
  exchange ACFDT corr.  =        -0.19030147  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  24)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.5026405E+00  (-0.8834272E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7640966 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4664.66912282
  -exchange      EXHF   =       115.62807665
  -V(xc)+E(xc)   XCENC  =       121.15854814
  PAW double counting   =      4106.87673309    -4097.88935383
  entropy T*S    EENTRO =        -0.04633153
  eigenvalues    EBANDS =      -647.12839375
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.74929511 eV

  energy without entropy =     -122.70296358  energy(sigma->0) =     -122.73385127
  exchange ACFDT corr.  =        -0.68222017  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  25)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.6829653E+00  (-0.4106834E+00)
 number of electron      51.0000043 magnetization 
 augmentation part        1.7905610 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4659.40883974
  -exchange      EXHF   =       115.62193112
  -V(xc)+E(xc)   XCENC  =       121.17623765
  PAW double counting   =      4110.00002928    -4101.00636307
  entropy T*S    EENTRO =        -0.04492599
  eigenvalues    EBANDS =      -653.09956044
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.43226038 eV

  energy without entropy =     -123.38733439  energy(sigma->0) =     -123.41728505
  exchange ACFDT corr.  =        -0.59122062  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  26)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4525662E+01  (-0.3440518E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7549094 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4631.03841515
  -exchange      EXHF   =       116.05827365
  -V(xc)+E(xc)   XCENC  =       121.50219825
  PAW double counting   =      4118.50570818    -4109.48574035
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -677.73433329
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.90659837 eV

  energy without entropy =     -118.85017941  energy(sigma->0) =     -118.88779205
  exchange ACFDT corr.  =        -0.30797114  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  27)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3399576E+01  (-0.9358738E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7598273 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4667.52806158
  -exchange      EXHF   =       115.65980309
  -V(xc)+E(xc)   XCENC  =       121.16963017
  PAW double counting   =      4105.94902540    -4096.95261744
  entropy T*S    EENTRO =        -0.05641879
  eigenvalues    EBANDS =      -643.87817177
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.30617474 eV

  energy without entropy =     -122.24975595  energy(sigma->0) =     -122.28736848
  exchange ACFDT corr.  =        -0.66801836  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  28)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7520397E+00  (-0.9849103E+00)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7845410 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4663.89765645
  -exchange      EXHF   =       115.67503011
  -V(xc)+E(xc)   XCENC  =       121.19755596
  PAW double counting   =      4111.67647606    -4102.68340863
  entropy T*S    EENTRO =        -0.04844700
  eigenvalues    EBANDS =      -648.30042904
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.05821444 eV

  energy without entropy =     -123.00976744  energy(sigma->0) =     -123.04206544
  exchange ACFDT corr.  =        -0.69735117  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  29)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2438615E+01  (-0.2708289E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7574482 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4633.95788017
  -exchange      EXHF   =       115.84028445
  -V(xc)+E(xc)   XCENC  =       121.41339718
  PAW double counting   =      4118.63979990    -4109.61831262
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -676.21907752
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.61959944 eV

  energy without entropy =     -120.56318048  energy(sigma->0) =     -120.60079312
  exchange ACFDT corr.  =        -0.36626209  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  30)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2233266E+01  (-0.4931346E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7610129 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4664.85076841
  -exchange      EXHF   =       115.67348950
  -V(xc)+E(xc)   XCENC  =       121.18538730
  PAW double counting   =      4109.75825393    -4100.75764873
  entropy T*S    EENTRO =        -0.05609599
  eigenvalues    EBANDS =      -647.13579609
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.85286511 eV

  energy without entropy =     -122.79676912  energy(sigma->0) =     -122.83416644
  exchange ACFDT corr.  =        -0.68479511  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  31)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4162693E+00  (-0.5731775E+00)
 number of electron      51.0000043 magnetization 
 augmentation part        1.7893001 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.04650348
  -exchange      EXHF   =       115.63843741
  -V(xc)+E(xc)   XCENC  =       121.16989811
  PAW double counting   =      4111.63080901    -4102.63157028
  entropy T*S    EENTRO =        -0.03158837
  eigenvalues    EBANDS =      -651.30474555
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.26913443 eV

  energy without entropy =     -123.23754606  energy(sigma->0) =     -123.25860497
  exchange ACFDT corr.  =        -0.70378091  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  32)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4563383E+01  (-0.1616146E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7519608 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.23935370
  -exchange      EXHF   =       116.05982575
  -V(xc)+E(xc)   XCENC  =       121.46609528
  PAW double counting   =      4120.36171905    -4111.33931255
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.31377356
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.70575176 eV

  energy without entropy =     -118.64933280  energy(sigma->0) =     -118.68694544
  exchange ACFDT corr.  =        -0.18727298  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  33)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3544264E+01  (-0.8843922E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7558559 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4667.76557682
  -exchange      EXHF   =       115.66508995
  -V(xc)+E(xc)   XCENC  =       121.17247017
  PAW double counting   =      4105.37870391    -4096.37529985
  entropy T*S    EENTRO =        -0.05641893
  eigenvalues    EBANDS =      -643.59962026
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.25001551 eV

  energy without entropy =     -122.19359657  energy(sigma->0) =     -122.23120920
  exchange ACFDT corr.  =        -0.66402339  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  34)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7200692E+00  (-0.9996615E+00)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7820385 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.68759439
  -exchange      EXHF   =       115.61108381
  -V(xc)+E(xc)   XCENC  =       121.13877956
  PAW double counting   =      4109.23290963    -4100.23207641
  entropy T*S    EENTRO =        -0.05333332
  eigenvalues    EBANDS =      -649.30740432
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.97008471 eV

  energy without entropy =     -122.91675139  energy(sigma->0) =     -122.95230694
  exchange ACFDT corr.  =        -0.69990893  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  35)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2385036E+01  (-0.2751150E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7548472 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4631.56410333
  -exchange      EXHF   =       115.76474138
  -V(xc)+E(xc)   XCENC  =       121.33995856
  PAW double counting   =      4116.44539436    -4107.42006311
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.42827931
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.58504842 eV

  energy without entropy =     -120.52862946  energy(sigma->0) =     -120.56624210
  exchange ACFDT corr.  =        -0.42988219  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  36)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2216611E+01  (-0.5156227E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7594038 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4664.84095005
  -exchange      EXHF   =       115.65612262
  -V(xc)+E(xc)   XCENC  =       121.17414256
  PAW double counting   =      4107.77162582    -4098.76491010
  entropy T*S    EENTRO =        -0.05562822
  eigenvalues    EBANDS =      -647.07190811
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.80165989 eV

  energy without entropy =     -122.74603167  energy(sigma->0) =     -122.78311715
  exchange ACFDT corr.  =        -0.68922979  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  37)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4322301E+00  (-0.5356017E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7899353 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.06421148
  -exchange      EXHF   =       115.62124355
  -V(xc)+E(xc)   XCENC  =       121.15836183
  PAW double counting   =      4110.39014806    -4101.38829833
  entropy T*S    EENTRO =        -0.03275817
  eigenvalues    EBANDS =      -651.22614174
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.23388999 eV

  energy without entropy =     -123.20113182  energy(sigma->0) =     -123.22297060
  exchange ACFDT corr.  =        -0.70730935  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  38)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4727517E+01  (-0.1632700E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7521223 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4628.71462746
  -exchange      EXHF   =       116.07871569
  -V(xc)+E(xc)   XCENC  =       121.46521200
  PAW double counting   =      4120.55023816    -4111.53098387
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.65280519
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.50637250 eV

  energy without entropy =     -118.44995354  energy(sigma->0) =     -118.48756618
  exchange ACFDT corr.  =        -0.20964622  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  39)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3712065E+01  (-0.8801668E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7552728 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4668.10214804
  -exchange      EXHF   =       115.67676255
  -V(xc)+E(xc)   XCENC  =       121.18466465
  PAW double counting   =      4105.51394479    -4096.50923059
  entropy T*S    EENTRO =        -0.05641892
  eigenvalues    EBANDS =      -643.25664848
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.21843774 eV

  energy without entropy =     -122.16201881  energy(sigma->0) =     -122.19963143
  exchange ACFDT corr.  =        -0.66596544  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  40)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7162849E+00  (-0.1003712E+01)
 number of electron      51.0000043 magnetization 
 augmentation part        1.7821427 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.47072403
  -exchange      EXHF   =       115.59697488
  -V(xc)+E(xc)   XCENC  =       121.12703462
  PAW double counting   =      4108.98301928    -4099.98105753
  entropy T*S    EENTRO =        -0.05048265
  eigenvalues    EBANDS =      -649.46418728
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.93472264 eV

  energy without entropy =     -122.88423999  energy(sigma->0) =     -122.91789509
  exchange ACFDT corr.  =        -0.70113890  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  41)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2682274E+01  (-0.2846737E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7540373 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4630.05782127
  -exchange      EXHF   =       115.78590068
  -V(xc)+E(xc)   XCENC  =       121.33048693
  PAW double counting   =      4116.71358332    -4107.68845823
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.61629351
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.25244840 eV

  energy without entropy =     -120.19602944  energy(sigma->0) =     -120.23364208
  exchange ACFDT corr.  =        -0.41603809  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  42)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2450644E+01  (-0.5736473E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7583306 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.49635464
  -exchange      EXHF   =       115.66311581
  -V(xc)+E(xc)   XCENC  =       121.17944411
  PAW double counting   =      4107.74065096    -4098.73260387
  entropy T*S    EENTRO =        -0.05612709
  eigenvalues    EBANDS =      -646.33156216
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.70309241 eV

  energy without entropy =     -122.64696531  energy(sigma->0) =     -122.68438337
  exchange ACFDT corr.  =        -0.68560210  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  43)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4788245E+00  (-0.6310733E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7890651 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.36225106
  -exchange      EXHF   =       115.61602519
  -V(xc)+E(xc)   XCENC  =       121.15396038
  PAW double counting   =      4110.29268363    -4101.29048277
  entropy T*S    EENTRO =        -0.03186409
  eigenvalues    EBANDS =      -650.86636152
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.18191690 eV

  energy without entropy =     -123.15005281  energy(sigma->0) =     -123.17129554
  exchange ACFDT corr.  =        -0.71055073  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  44)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4657180E+01  (-0.3140789E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7520668 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4628.45772163
  -exchange      EXHF   =       116.07041271
  -V(xc)+E(xc)   XCENC  =       121.45495314
  PAW double counting   =      4120.39215709    -4111.37331864
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.90999150
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.52473657 eV

  energy without entropy =     -118.46831761  energy(sigma->0) =     -118.50593025
  exchange ACFDT corr.  =        -0.19820341  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  45)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3719019E+01  (-0.8656517E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7551687 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4668.38077754
  -exchange      EXHF   =       115.69615935
  -V(xc)+E(xc)   XCENC  =       121.20187802
  PAW double counting   =      4105.83978501    -4096.83312283
  entropy T*S    EENTRO =        -0.05641891
  eigenvalues    EBANDS =      -643.04189458
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.24375519 eV

  energy without entropy =     -122.18733627  energy(sigma->0) =     -122.22494888
  exchange ACFDT corr.  =        -0.66524662  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  46)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7018925E+00  (-0.9714644E+00)
 number of electron      51.0000043 magnetization 
 augmentation part        1.7820537 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.47095513
  -exchange      EXHF   =       115.60175565
  -V(xc)+E(xc)   XCENC  =       121.13174998
  PAW double counting   =      4109.07515552    -4100.07263672
  entropy T*S    EENTRO =        -0.05030922
  eigenvalues    EBANDS =      -649.48493443
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.94564770 eV

  energy without entropy =     -122.89533848  energy(sigma->0) =     -122.92887796
  exchange ACFDT corr.  =        -0.70099273  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  47)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2714711E+01  (-0.2850673E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7538845 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.99939264
  -exchange      EXHF   =       115.78898833
  -V(xc)+E(xc)   XCENC  =       121.33095647
  PAW double counting   =      4116.84412035    -4107.81985984
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.65607627
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.23093647 eV

  energy without entropy =     -120.17451751  energy(sigma->0) =     -120.21213015
  exchange ACFDT corr.  =        -0.41429924  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  48)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2462185E+01  (-0.5832940E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7581800 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.81955271
  -exchange      EXHF   =       115.67771663
  -V(xc)+E(xc)   XCENC  =       121.19099712
  PAW double counting   =      4107.93353855    -4098.92448343
  entropy T*S    EENTRO =        -0.05616256
  eigenvalues    EBANDS =      -646.02555508
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.69312153 eV

  energy without entropy =     -122.63695897  energy(sigma->0) =     -122.67440068
  exchange ACFDT corr.  =        -0.68252530  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  49)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4849021E+00  (-0.6360316E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7890365 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.46775023
  -exchange      EXHF   =       115.62053991
  -V(xc)+E(xc)   XCENC  =       121.15764417
  PAW double counting   =      4110.43955530    -4101.43745143
  entropy T*S    EENTRO =        -0.03159356
  eigenvalues    EBANDS =      -650.76503510
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.17802360 eV

  energy without entropy =     -123.14643004  energy(sigma->0) =     -123.16749242
  exchange ACFDT corr.  =        -0.70948808  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  50)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4643955E+01  (-0.4034129E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7521550 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4628.52688317
  -exchange      EXHF   =       116.07308296
  -V(xc)+E(xc)   XCENC  =       121.45684683
  PAW double counting   =      4120.59040450    -4111.57249829
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.85406390
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.53406828 eV

  energy without entropy =     -118.47764932  energy(sigma->0) =     -118.51526196
  exchange ACFDT corr.  =        -0.19320040  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  51)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3712806E+01  (-0.8770604E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7552031 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4668.63005507
  -exchange      EXHF   =       115.71315595
  -V(xc)+E(xc)   XCENC  =       121.21456385
  PAW double counting   =      4106.03439200    -4097.02705081
  entropy T*S    EENTRO =        -0.05641892
  eigenvalues    EBANDS =      -642.82609735
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.24687404 eV

  energy without entropy =     -122.19045512  energy(sigma->0) =     -122.22806773
  exchange ACFDT corr.  =        -0.66155780  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  52)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7065489E+00  (-0.9640606E+00)
 number of electron      51.0000043 magnetization 
 augmentation part        1.7818744 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.45285601
  -exchange      EXHF   =       115.60674646
  -V(xc)+E(xc)   XCENC  =       121.13496960
  PAW double counting   =      4109.18835769    -4100.18569245
  entropy T*S    EENTRO =        -0.05255366
  eigenvalues    EBANDS =      -649.51916563
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.95342291 eV

  energy without entropy =     -122.90086925  energy(sigma->0) =     -122.93590502
  exchange ACFDT corr.  =        -0.69918336  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  53)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2526843E+01  (-0.2805989E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7543109 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4630.60492397
  -exchange      EXHF   =       115.76803571
  -V(xc)+E(xc)   XCENC  =       121.32754913
  PAW double counting   =      4117.00800118    -4107.98492013
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.21840498
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.42658038 eV

  energy without entropy =     -120.37016143  energy(sigma->0) =     -120.40777407
  exchange ACFDT corr.  =        -0.42887576  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  54)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2310595E+01  (-0.5648117E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7588780 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.73527702
  -exchange      EXHF   =       115.68832924
  -V(xc)+E(xc)   XCENC  =       121.19874459
  PAW double counting   =      4108.46425052    -4099.45511127
  entropy T*S    EENTRO =        -0.05590975
  eigenvalues    EBANDS =      -646.17232921
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.73717577 eV

  energy without entropy =     -122.68126603  energy(sigma->0) =     -122.71853919
  exchange ACFDT corr.  =        -0.68089733  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  55)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4679362E+00  (-0.5840140E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7902155 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.40266280
  -exchange      EXHF   =       115.62854667
  -V(xc)+E(xc)   XCENC  =       121.16465528
  PAW double counting   =      4111.02660454    -4102.02543477
  entropy T*S    EENTRO =        -0.03231228
  eigenvalues    EBANDS =      -650.87154749
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.20511197 eV

  energy without entropy =     -123.17279969  energy(sigma->0) =     -123.19434121
  exchange ACFDT corr.  =        -0.70725282  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  56)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4699354E+01  (-0.2487140E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7524755 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4628.86415314
  -exchange      EXHF   =       116.08874955
  -V(xc)+E(xc)   XCENC  =       121.47089579
  PAW double counting   =      4121.54717184    -4112.53202061
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.51472568
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.50575819 eV

  energy without entropy =     -118.44933923  energy(sigma->0) =     -118.48695187
  exchange ACFDT corr.  =        -0.20046866  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  57)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3705574E+01  (-0.9278295E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7554315 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4668.90524694
  -exchange      EXHF   =       115.72571584
  -V(xc)+E(xc)   XCENC  =       121.22178246
  PAW double counting   =      4106.38550976    -4097.37898906
  entropy T*S    EENTRO =        -0.05641893
  eigenvalues    EBANDS =      -642.53432119
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.21133175 eV

  energy without entropy =     -122.15491282  energy(sigma->0) =     -122.19252544
  exchange ACFDT corr.  =        -0.65622842  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  58)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7407709E+00  (-0.9848515E+00)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7821766 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.43413636
  -exchange      EXHF   =       115.61005441
  -V(xc)+E(xc)   XCENC  =       121.13574996
  PAW double counting   =      4109.56388393    -4100.56173078
  entropy T*S    EENTRO =        -0.05466307
  eigenvalues    EBANDS =      -649.54014122
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.95210264 eV

  energy without entropy =     -122.89743957  energy(sigma->0) =     -122.93388162
  exchange ACFDT corr.  =        -0.69693172  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  59)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2368197E+01  (-0.2764290E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7549062 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4631.21234714
  -exchange      EXHF   =       115.75298995
  -V(xc)+E(xc)   XCENC  =       121.32770477
  PAW double counting   =      4117.59320191    -4108.57187227
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.74955578
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.58390528 eV

  energy without entropy =     -120.52748632  energy(sigma->0) =     -120.56509896
  exchange ACFDT corr.  =        -0.44332485  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  60)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2191998E+01  (-0.5493880E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7596764 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.60138110
  -exchange      EXHF   =       115.69415394
  -V(xc)+E(xc)   XCENC  =       121.20293128
  PAW double counting   =      4109.19682375    -4100.18838680
  entropy T*S    EENTRO =        -0.05554435
  eigenvalues    EBANDS =      -646.35426202
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.77590356 eV

  energy without entropy =     -122.72035921  energy(sigma->0) =     -122.75738878
  exchange ACFDT corr.  =        -0.68067862  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  61)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4528852E+00  (-0.5271805E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7914535 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.34917892
  -exchange      EXHF   =       115.63529541
  -V(xc)+E(xc)   XCENC  =       121.17073267
  PAW double counting   =      4111.78235587    -4102.78225483
  entropy T*S    EENTRO =        -0.03303197
  eigenvalues    EBANDS =      -650.96083100
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.22878880 eV

  energy without entropy =     -123.19575684  energy(sigma->0) =     -123.21777815
  exchange ACFDT corr.  =        -0.70532767  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  62)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4762657E+01  (-0.1494073E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7525633 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.18563092
  -exchange      EXHF   =       116.10464821
  -V(xc)+E(xc)   XCENC  =       121.48520640
  PAW double counting   =      4122.54966505    -4113.53710861
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -679.18051641
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.46613189 eV

  energy without entropy =     -118.40971294  energy(sigma->0) =     -118.44732557
  exchange ACFDT corr.  =        -0.20778170  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  63)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3700252E+01  (-0.9851334E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555016 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4669.05054636
  -exchange      EXHF   =       115.73026707
  -V(xc)+E(xc)   XCENC  =       121.22223783
  PAW double counting   =      4106.59759522    -4097.59188123
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -642.34827412
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.16638420 eV

  energy without entropy =     -122.10996527  energy(sigma->0) =     -122.14757789
  exchange ACFDT corr.  =        -0.65140373  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  64)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.7805674E+00  (-0.1011054E+01)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7825263 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.39956525
  -exchange      EXHF   =       115.61071093
  -V(xc)+E(xc)   XCENC  =       121.13448816
  PAW double counting   =      4109.84351505    -4100.84144512
  entropy T*S    EENTRO =        -0.05610794
  eigenvalues    EBANDS =      -649.56887280
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.94695162 eV

  energy without entropy =     -122.89084368  energy(sigma->0) =     -122.92824897
  exchange ACFDT corr.  =        -0.69516228  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  65)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2267873E+01  (-0.2735311E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7552820 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4631.63964760
  -exchange      EXHF   =       115.74437747
  -V(xc)+E(xc)   XCENC  =       121.32899025
  PAW double counting   =      4118.11323098    -4109.09332771
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.40723056
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.67907878 eV

  energy without entropy =     -120.62265982  energy(sigma->0) =     -120.66027246
  exchange ACFDT corr.  =        -0.45365138  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  66)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2121433E+01  (-0.5399338E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7602191 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.48456633
  -exchange      EXHF   =       115.69622877
  -V(xc)+E(xc)   XCENC  =       121.20432818
  PAW double counting   =      4109.70911937    -4100.70116372
  entropy T*S    EENTRO =        -0.05520165
  eigenvalues    EBANDS =      -646.49867588
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.80051222 eV

  energy without entropy =     -122.74531057  energy(sigma->0) =     -122.78211167
  exchange ACFDT corr.  =        -0.68093071  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  67)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4428077E+00  (-0.4824279E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7923790 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.30893518
  -exchange      EXHF   =       115.63901247
  -V(xc)+E(xc)   XCENC  =       121.17415957
  PAW double counting   =      4112.31426464    -4103.31493245
  entropy T*S    EENTRO =        -0.03357934
  eigenvalues    EBANDS =      -651.02232364
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.24331991 eV

  energy without entropy =     -123.20974056  energy(sigma->0) =     -123.23212679
  exchange ACFDT corr.  =        -0.70380390  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  68)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4809471E+01  (-0.9060686E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7526397 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.40506009
  -exchange      EXHF   =       116.11594800
  -V(xc)+E(xc)   XCENC  =       121.49529553
  PAW double counting   =      4123.26334621    -4114.25273529
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.94770051
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.43384916 eV

  energy without entropy =     -118.37743020  energy(sigma->0) =     -118.41504284
  exchange ACFDT corr.  =        -0.21332356  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  69)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3699750E+01  (-0.1027574E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555481 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4669.16860647
  -exchange      EXHF   =       115.73321924
  -V(xc)+E(xc)   XCENC  =       121.22282081
  PAW double counting   =      4106.74780558    -4097.74258959
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -642.20046650
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.13359954 eV

  energy without entropy =     -122.07718060  energy(sigma->0) =     -122.11479323
  exchange ACFDT corr.  =        -0.64870221  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  70)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.8096465E+00  (-0.1029006E+01)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7827453 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.35931300
  -exchange      EXHF   =       115.61015173
  -V(xc)+E(xc)   XCENC  =       121.13292871
  PAW double counting   =      4110.03456932    -4101.03256972
  entropy T*S    EENTRO =        -0.05700951
  eigenvalues    EBANDS =      -649.60323050
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.94324607 eV

  energy without entropy =     -122.88623655  energy(sigma->0) =     -122.92424289
  exchange ACFDT corr.  =        -0.69430554  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  71)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2205079E+01  (-0.2712511E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555038 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4631.89337515
  -exchange      EXHF   =       115.73826406
  -V(xc)+E(xc)   XCENC  =       121.32912892
  PAW double counting   =      4118.45915916    -4109.44012009
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.20485064
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.73816692 eV

  energy without entropy =     -120.68174796  energy(sigma->0) =     -120.71936060
  exchange ACFDT corr.  =        -0.46005533  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  72)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2077249E+01  (-0.5350710E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7605293 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.40878259
  -exchange      EXHF   =       115.69723187
  -V(xc)+E(xc)   XCENC  =       121.20497097
  PAW double counting   =      4110.04776247    -4101.04008309
  entropy T*S    EENTRO =        -0.05491301
  eigenvalues    EBANDS =      -646.59073268
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.81541567 eV

  energy without entropy =     -122.76050266  energy(sigma->0) =     -122.79711133
  exchange ACFDT corr.  =        -0.68115028  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  73)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4367977E+00  (-0.4490074E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7929369 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.29807044
  -exchange      EXHF   =       115.64184556
  -V(xc)+E(xc)   XCENC  =       121.17670798
  PAW double counting   =      4112.63727297    -4103.63835650
  entropy T*S    EENTRO =        -0.03402659
  eigenvalues    EBANDS =      -651.04733628
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.25221337 eV

  energy without entropy =     -123.21818679  energy(sigma->0) =     -123.24087118
  exchange ACFDT corr.  =        -0.70237723  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  74)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4837908E+01  (-0.5068840E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7526643 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.55412777
  -exchange      EXHF   =       116.12303864
  -V(xc)+E(xc)   XCENC  =       121.50188704
  PAW double counting   =      4123.68825425    -4114.67881777
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.79114962
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.41430548 eV

  energy without entropy =     -118.35788652  energy(sigma->0) =     -118.39549916
  exchange ACFDT corr.  =        -0.21785178  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  75)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3698661E+01  (-0.1055428E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555575 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4669.24833622
  -exchange      EXHF   =       115.73548326
  -V(xc)+E(xc)   XCENC  =       121.22351242
  PAW double counting   =      4106.80572235    -4097.80073592
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -642.10282966
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.11296638 eV

  energy without entropy =     -122.05654744  energy(sigma->0) =     -122.09416007
  exchange ACFDT corr.  =        -0.64695978  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  76)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.8284306E+00  (-0.1039795E+01)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7828724 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.32409585
  -exchange      EXHF   =       115.60964028
  -V(xc)+E(xc)   XCENC  =       121.13178291
  PAW double counting   =      4110.12457593    -4101.12254232
  entropy T*S    EENTRO =        -0.05758359
  eigenvalues    EBANDS =      -649.63497533
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.94139697 eV

  energy without entropy =     -122.88381338  energy(sigma->0) =     -122.92220244
  exchange ACFDT corr.  =        -0.69374810  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  77)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2165204E+01  (-0.2696058E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7556477 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4632.05427956
  -exchange      EXHF   =       115.73414505
  -V(xc)+E(xc)   XCENC  =       121.32904230
  PAW double counting   =      4118.66691284    -4109.64837716
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.07668879
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.77619257 eV

  energy without entropy =     -120.71977361  energy(sigma->0) =     -120.75738625
  exchange ACFDT corr.  =        -0.46412344  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  78)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2048699E+01  (-0.5324604E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7607212 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.36751574
  -exchange      EXHF   =       115.69811499
  -V(xc)+E(xc)   XCENC  =       121.20561180
  PAW double counting   =      4110.26018038    -4101.25265265
  entropy T*S    EENTRO =        -0.05469876
  eigenvalues    EBANDS =      -646.64284733
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.82489116 eV

  energy without entropy =     -122.77019240  energy(sigma->0) =     -122.80665824
  exchange ACFDT corr.  =        -0.68126128  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  79)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4330510E+00  (-0.4255442E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7932880 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.30469715
  -exchange      EXHF   =       115.64416116
  -V(xc)+E(xc)   XCENC  =       121.17878482
  PAW double counting   =      4112.82814802    -4103.82949175
  entropy T*S    EENTRO =        -0.03431767
  eigenvalues    EBANDS =      -651.05078482
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.25794214 eV

  energy without entropy =     -123.22362447  energy(sigma->0) =     -123.24650292
  exchange ACFDT corr.  =        -0.70124899  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  80)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4857020E+01  (-0.2872029E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7526358 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.65842503
  -exchange      EXHF   =       116.12790315
  -V(xc)+E(xc)   XCENC  =       121.50645179
  PAW double counting   =      4123.93771217    -4114.92900806
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.68187492
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.40092224 eV

  energy without entropy =     -118.34450328  energy(sigma->0) =     -118.38211592
  exchange ACFDT corr.  =        -0.22075546  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  81)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3697653E+01  (-0.1074822E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555300 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4669.29272788
  -exchange      EXHF   =       115.73662901
  -V(xc)+E(xc)   XCENC  =       121.22361760
  PAW double counting   =      4106.79493409    -4097.79002843
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -642.04521732
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.09857556 eV

  energy without entropy =     -122.04215662  energy(sigma->0) =     -122.07976925
  exchange ACFDT corr.  =        -0.64572828  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  82)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.8414706E+00  (-0.1047154E+01)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7829434 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.29243006
  -exchange      EXHF   =       115.60897190
  -V(xc)+E(xc)   XCENC  =       121.13072651
  PAW double counting   =      4110.14770983    -4101.14556939
  entropy T*S    EENTRO =        -0.05796294
  eigenvalues    EBANDS =      -649.66367230
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.94004612 eV

  energy without entropy =     -122.88208318  energy(sigma->0) =     -122.92072514
  exchange ACFDT corr.  =        -0.69335716  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  83)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.2138486E+01  (-0.2684176E+02)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7557393 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4632.16168185
  -exchange      EXHF   =       115.73131099
  -V(xc)+E(xc)   XCENC  =       121.32892115
  PAW double counting   =      4118.78518583    -4109.76694863
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -677.99102053
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -120.80155964 eV

  energy without entropy =     -120.74514068  energy(sigma->0) =     -120.78275332
  exchange ACFDT corr.  =        -0.46681040  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  84)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.2029477E+01  (-0.5311301E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7608442 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4665.34479862
  -exchange      EXHF   =       115.69876004
  -V(xc)+E(xc)   XCENC  =       121.20612590
  PAW double counting   =      4110.38997249    -4101.38255184
  entropy T*S    EENTRO =        -0.05453655
  eigenvalues    EBANDS =      -646.67276156
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.83103620 eV

  energy without entropy =     -122.77649965  energy(sigma->0) =     -122.81285735
  exchange ACFDT corr.  =        -0.68134545  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  85)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.4307143E+00  (-0.4085397E+00)
 number of electron      51.0000042 magnetization 
 augmentation part        1.7935180 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4661.31772491
  -exchange      EXHF   =       115.64598630
  -V(xc)+E(xc)   XCENC  =       121.18042988
  PAW double counting   =      4112.94166368    -4103.94316969
  entropy T*S    EENTRO =        -0.03453524
  eigenvalues    EBANDS =      -651.04503560
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -123.26175054 eV

  energy without entropy =     -123.22721530  energy(sigma->0) =     -123.25023879
  exchange ACFDT corr.  =        -0.70036904  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  86)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) : 0.4869853E+01  (-0.1363324E+00)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7525995 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4629.73282469
  -exchange      EXHF   =       116.13123999
  -V(xc)+E(xc)   XCENC  =       121.50963992
  PAW double counting   =      4124.08686604    -4115.07862287
  entropy T*S    EENTRO =        -0.05641896
  eigenvalues    EBANDS =      -678.60429720
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -118.39189772 eV

  energy without entropy =     -118.33547876  energy(sigma->0) =     -118.37309140
  exchange ACFDT corr.  =        -0.22292963  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  87)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.3696962E+01  (-0.1088043E+01)
 number of electron      51.0000051 magnetization 
 augmentation part        1.7555000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4669.32247898
  -exchange      EXHF   =       115.73738945
  -V(xc)+E(xc)   XCENC  =       121.22367611
  PAW double counting   =      4106.76547758    -4097.76057953
  entropy T*S    EENTRO =        -0.05641894
  eigenvalues    EBANDS =      -642.00656128
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.08885927 eV

  energy without entropy =     -122.03244033  energy(sigma->0) =     -122.07005296
  exchange ACFDT corr.  =        -0.64487946  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  88)  ---------------------------------------



 eigenvalue-minimisations  :     0
 total energy-change (2. order) :-0.8503072E+00  (-0.1052066E+01)
 number of electron      51.0000044 magnetization 
 augmentation part        1.7829850 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18207684
  Ewald energy   TEWEN  =      2969.37918699
  -Hartree energ DENC   =     -4662.26773138
  -exchange      EXHF   =       115.60837700
  -V(xc)+E(xc)   XCENC  =       121.12989031
  PAW double counting   =      4110.14633288    -4101.14409650
  entropy T*S    EENTRO =        -0.05821977
  eigenvalues    EBANDS =      -649.68615614
  atomic energy  EATOM  =      1971.76937252
  ---------------------------------------------------
  free energy    TOTEN  =      -122.93916643 eV

  energy without entropy =     -122.88094666  energy(sigma->0) =     -122.91975984
  exchange ACFDT corr.  =        -0.69308078  see jH, gK, PRB 81, 115126


--------------------------------------------------------------------------------------------------------




--------------------------------------- Iteration      1(  89)  ---------------------------------------