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.18  20:16:42
 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
   ALGO = Normal (blocked Davidson)
   ISPIN = 2
   INIWAV = 1
   ISTART = 0
   ICHARG = 2
   LWAVE = .FALSE.
   LCHARG = .FALSE.
   ADDGRID = .FALSE.
   ISMEAR = 1
   SIGMA = 0.2
   LREAL = Auto
   LSCALAPACK = .FALSE.
   RWIGS = 0.75 0.32 0.77 1.11 0.73
   NPAR = 2

 POTCAR:    PAW_PBE N 08Apr2002                   
 POTCAR:    PAW_PBE H 15Jun2001                   
 POTCAR:    PAW_PBE C 08Apr2002                   
 POTCAR:    PAW_PBE Si 05Jan2001                  
 POTCAR:    PAW_PBE O 08Apr2002                   
 -----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|     You use a magnetic or noncollinear calculation, but did not specify     |
|     the initial magnetic moment with the MAGMOM tag. Note that a            |
|     default of 1 will be used for all atoms. This ferromagnetic setup       |
|     may break the symmetry of the crystal, in particular it may rule        |
|     out finding an antiferromagnetic solution. Thence, we recommend         |
|     setting the initial magnetic moment manually or verifying carefully     |
|     that this magnetic setup is desired.                                    |
|                                                                             |
 -----------------------------------------------------------------------------

 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
 
 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 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 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
 
 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
 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.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         = 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
  PAW_PBE N 08Apr2002                   :
 energy of atom  1       EATOM= -264.5486
 kinetic energy error for atom=    0.0736 (will be added to EATOM!!)
  PAW_PBE H 15Jun2001                   :
 energy of atom  2       EATOM=  -12.4884
 kinetic energy error for atom=    0.0098 (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 Si 05Jan2001                  :
 energy of atom  4       EATOM= -103.0669
 kinetic energy error for atom=    0.0012 (will be added to EATOM!!)
  PAW_PBE O 08Apr2002                   :
 energy of atom  5       EATOM= -432.3788
 kinetic energy error for atom=    0.1156 (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.719  0.440  0.592-   3 1.04   2 1.05
   2  0.742  0.452  0.654-   1 1.05
   3  0.760  0.442  0.549-   1 1.04
   4  0.149  0.518  0.366-  25 1.10
   5  0.205  0.612  0.369-  25 1.10
   6  0.220  0.522  0.293-  25 1.10
   7  0.109  0.550  0.657-  26 1.10
   8  0.149  0.639  0.599-  26 1.11
   9  0.105  0.554  0.539-  26 1.10
  10  0.418  0.432  0.295-  27 1.10
  11  0.417  0.342  0.372-  27 1.10
  12  0.344  0.410  0.359-  27 1.10
  13  0.532  0.276  0.574-  28 1.11
  14  0.532  0.349  0.478-  28 1.10
  15  0.555  0.390  0.587-  28 1.12
  16  0.246  0.691  0.505-  29 1.10
  17  0.303  0.667  0.592-  29 1.10
  18  0.328  0.656  0.478-  29 1.10
  19  0.440  0.620  0.559-  30 1.10
  20  0.449  0.557  0.659-  30 1.10
  21  0.512  0.552  0.575-  30 1.10
  22  0.323  0.446  0.641-  24 1.10
  23  0.308  0.378  0.548-  24 1.10
  24  0.325  0.445  0.568-  22 1.10  23 1.10  32 1.86  31 1.87
  25  0.202  0.538  0.361-   6 1.10   4 1.10   5 1.10  33 1.42
  26  0.137  0.567  0.597-   7 1.10   9 1.10   8 1.11  34 1.42
  27  0.400  0.410  0.360-  10 1.10  12 1.10  11 1.10  35 1.43
  28  0.521  0.344  0.550-  14 1.10  13 1.11  15 1.12  36 1.44
  29  0.287  0.646  0.525-  18 1.10  17 1.10  16 1.10  32 1.88
  30  0.457  0.557  0.586-  19 1.10  21 1.10  20 1.10  31 1.86
  31  0.415  0.459  0.535-  35 1.66  36 1.68  30 1.86  24 1.87
  32  0.261  0.525  0.526-  34 1.66  33 1.66  24 1.86  29 1.88
  33  0.243  0.492  0.423-  25 1.42  32 1.66
  34  0.196  0.513  0.595-  26 1.42  32 1.66
  35  0.424  0.471  0.425-  27 1.43  31 1.66
  36  0.453  0.366  0.571-  28 1.44  31 1.68
 
  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 .


Analysis of structural, dynamic, and magnetic symmetry:
=====================================================================
 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 overall 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.719219550  0.440279480  0.592010130
     0.741651470  0.452456050  0.654305900
     0.760066500  0.441862240  0.549292160
     0.149447780  0.517673390  0.366479290
     0.204864450  0.611562930  0.368793330
     0.219663260  0.522085130  0.292950080
     0.108790890  0.550254020  0.657298620
     0.148770500  0.638859070  0.598508400
     0.104905910  0.553632120  0.538516390
     0.418326670  0.431793390  0.294545690
     0.417182140  0.341604460  0.372396240
     0.344444250  0.410495570  0.358513530
     0.531596330  0.276320590  0.574045900
     0.532389210  0.349225920  0.477829450
     0.555302930  0.389884020  0.586662650
     0.245725090  0.690830040  0.504641100
     0.302975510  0.667381920  0.591632880
     0.328320350  0.656491570  0.478371580
     0.439569640  0.620395330  0.558623980
     0.448638690  0.556588320  0.658778350
     0.511537450  0.552111200  0.575051740
     0.323040750  0.446281240  0.641238590
     0.307923040  0.378375820  0.547693610
     0.324909900  0.445482540  0.567782310
     0.202220940  0.538432490  0.360652820
     0.137247050  0.566830580  0.596693550
     0.399537640  0.410164790  0.359948630
     0.520801270  0.344325030  0.549592700
     0.286574950  0.645938730  0.524947740
     0.457322470  0.557004790  0.586357550
     0.414673890  0.459325010  0.534687840
     0.260752580  0.525486680  0.526494320
     0.242921300  0.492251970  0.423346220
     0.195999920  0.513373940  0.594883480
     0.424194720  0.471121930  0.425301900
     0.452747030  0.366144590  0.571041000

  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
                23          23
                24          24
                25          25
                26          26
                27          27
                28          28
                29          29
                30          30
                31          31
                32          32
                33          33
                34          34
                35          35
                36          36

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

 
 
 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
 


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




 Dimension of arrays:
   k-points           NKPTS =      1   k-points in BZ     NKDIM =      1   number of bands    NBANDS=     80
   number of dos      NEDOS =    301   number of ions     NIONS =     36
   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  22   7   2   4
   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  =      2    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   =      2    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 =  14.00  1.00 12.01 28.09 16.00
  Ionic Valenz
   ZVAL   =   5.00  1.00  4.00  4.00  6.00
  Atomic Wigner-Seitz radii
   RWIGS  =   0.75  0.32  0.77  1.11  0.73
  virtual crystal weights 
   VCA    =   1.00  1.00  1.00  1.00  1.00
   NELECT =      87.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  =     38    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.31E-07  absolut break condition
   DEPER  =   0.30     relativ break condition  

   TIME   =   0.40     timestep for ELM

  volume/ion in A,a.u.               =     125.00       843.54
  Fermi-wavevector in a.u.,A,eV,Ry     =   0.439383  0.830314  2.626712  0.193058
  Thomas-Fermi vector in A             =   1.413435
 
 Write flags
   LWAVE        =      F    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
   LHFCALC =     F    Hartree Fock is set to
   LHFONE  =     F    Hartree Fock one center treatment
   AEXX    =    0.0000 exact exchange contribution

 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
 spin polarized calculation
 Variant of blocked Davidson
 Davidson routine will perform the subspace rotation
 perform sub-space diagonalisation
    after 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           36
 real space projection scheme for non local part
 use partial core corrections
 calculate Harris-corrections to forces 
   (improved forces if not selfconsistent)
 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.71921955  0.44027948  0.59201013
   0.74165147  0.45245605  0.65430590
   0.76006650  0.44186224  0.54929216
   0.14944778  0.51767339  0.36647929
   0.20486445  0.61156293  0.36879333
   0.21966326  0.52208513  0.29295008
   0.10879089  0.55025402  0.65729862
   0.14877050  0.63885907  0.59850840
   0.10490591  0.55363212  0.53851639
   0.41832667  0.43179339  0.29454569
   0.41718214  0.34160446  0.37239624
   0.34444425  0.41049557  0.35851353
   0.53159633  0.27632059  0.57404590
   0.53238921  0.34922592  0.47782945
   0.55530293  0.38988402  0.58666265
   0.24572509  0.69083004  0.50464110
   0.30297551  0.66738192  0.59163288
   0.32832035  0.65649157  0.47837158
   0.43956964  0.62039533  0.55862398
   0.44863869  0.55658832  0.65877835
   0.51153745  0.55211120  0.57505174
   0.32304075  0.44628124  0.64123859
   0.30792304  0.37837582  0.54769361
   0.32490990  0.44548254  0.56778231
   0.20222094  0.53843249  0.36065282
   0.13724705  0.56683058  0.59669355
   0.39953764  0.41016479  0.35994863
   0.52080127  0.34432503  0.54959270
   0.28657495  0.64593873  0.52494774
   0.45732247  0.55700479  0.58635755
   0.41467389  0.45932501  0.53468784
   0.26075258  0.52548668  0.52649432
   0.24292130  0.49225197  0.42334622
   0.19599992  0.51337394  0.59488348
   0.42419472  0.47112193  0.42530190
   0.45274703  0.36614459  0.57104100
 
 position of ions in cartesian coordinates  (Angst):
  14.38439100  6.60419220  8.88015195
  14.83302940  6.78684075  9.81458850
  15.20133000  6.62793360  8.23938240
   2.98895560  7.76510085  5.49718935
   4.09728900  9.17344395  5.53189995
   4.39326520  7.83127695  4.39425120
   2.17581780  8.25381030  9.85947930
   2.97541000  9.58288605  8.97762600
   2.09811820  8.30448180  8.07774585
   8.36653340  6.47690085  4.41818535
   8.34364280  5.12406690  5.58594360
   6.88888500  6.15743355  5.37770295
  10.63192660  4.14480885  8.61068850
  10.64778420  5.23838880  7.16744175
  11.10605860  5.84826030  8.79993975
   4.91450180 10.36245060  7.56961650
   6.05951020 10.01072880  8.87449320
   6.56640700  9.84737355  7.17557370
   8.79139280  9.30592995  8.37935970
   8.97277380  8.34882480  9.88167525
  10.23074900  8.28166800  8.62577610
   6.46081500  6.69421860  9.61857885
   6.15846080  5.67563730  8.21540415
   6.49819800  6.68223810  8.51673465
   4.04441880  8.07648735  5.40979230
   2.74494100  8.50245870  8.95040325
   7.99075280  6.15247185  5.39922945
  10.41602540  5.16487545  8.24389050
   5.73149900  9.68908095  7.87421610
   9.14644940  8.35507185  8.79536325
   8.29347780  6.88987515  8.02031760
   5.21505160  7.88230020  7.89741480
   4.85842600  7.38377955  6.35019330
   3.91999840  7.70060910  8.92325220
   8.48389440  7.06682895  6.37952850
   9.05494060  5.49216885  8.56561500
 


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


 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)


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

   base      :      30000. kBytes
   nonlr-proj:       4077. kBytes
   fftplans  :      87706. kBytes
   grid      :     314101. kBytes
   one-center:        221. kBytes
   wavefun   :     105128. 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      87.0000000 magnetization      36.0000000
 keeping initial charge density in first step


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


 Maximum index for non-local projection operator         2542
 Maximum index for augmentation-charges         2263 (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  :   384
 total energy-change (2. order) : 0.6880075E+03  (-0.2067122E+04)
 number of electron      87.0000000 magnetization      36.0000000
 augmentation part       87.0000000 magnetization      36.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8222.29457434
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       263.93852034
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.01798507
  eigenvalues    EBANDS =      -541.98458791
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       688.00745458 eV

  energy without entropy =      688.02543965  energy(sigma->0) =      688.01344960


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




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



 eigenvalue-minimisations  :   544
 total energy-change (2. order) :-0.5837745E+03  (-0.5269159E+03)
 number of electron      87.0000000 magnetization      36.0000000
 augmentation part       87.0000000 magnetization      36.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8222.29457434
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       263.93852034
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.01465677
  eigenvalues    EBANDS =     -1125.79176406
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       104.23292028 eV

  energy without entropy =      104.21826350  energy(sigma->0) =      104.22803468


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.2505412E+03  (-0.2489141E+03)
 number of electron      87.0000000 magnetization      36.0000000
 augmentation part       87.0000000 magnetization      36.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8222.29457434
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       263.93852034
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.01371475
  eigenvalues    EBANDS =     -1376.33198317
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -146.30824086 eV

  energy without entropy =     -146.32195561  energy(sigma->0) =     -146.31281244


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.2284357E+02  (-0.2272304E+02)
 number of electron      87.0000000 magnetization      36.0000000
 augmentation part       87.0000000 magnetization      36.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8222.29457434
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       263.93852034
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.01642370
  eigenvalues    EBANDS =     -1399.17825719
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -169.15180593 eV

  energy without entropy =     -169.16822963  energy(sigma->0) =     -169.15728049


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) :-0.4960715E+00  (-0.4947756E+00)
 number of electron      87.0000001 magnetization      30.9695110
 augmentation part        4.2557863 magnetization      30.1098455

 Broyden mixing:
  rms(total) = 0.41763E+01    rms(broyden)= 0.41735E+01
  rms(prec ) = 0.43484E+01
  weight for this iteration     100.00

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8222.29457434
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       263.93852034
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.01657284
  eigenvalues    EBANDS =     -1399.67447783
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -169.64787743 eV

  energy without entropy =     -169.66445027  energy(sigma->0) =     -169.65340171


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




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



 eigenvalue-minimisations  :   376
 total energy-change (2. order) : 0.9912806E+02  (-0.2422271E+02)
 number of electron      87.0000001 magnetization      26.4554883
 augmentation part        3.8582514 magnetization      25.1472797

 Broyden mixing:
  rms(total) = 0.20303E+01    rms(broyden)= 0.20292E+01
  rms(prec ) = 0.21145E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9139
  0.9139

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8424.60246969
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       360.79651518
  PAW double counting   =      4227.23018571    -4206.74280432
  entropy T*S    EENTRO =         0.01542672
  eigenvalues    EBANDS =     -1192.85239391
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -70.51981954 eV

  energy without entropy =      -70.53524626  energy(sigma->0) =      -70.52496178


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




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



 eigenvalue-minimisations  :   416
 total energy-change (2. order) :-0.4993400E+02  (-0.4013871E+01)
 number of electron      87.0000001 magnetization      22.4586694
 augmentation part        3.5807524 magnetization      21.2188663

 Broyden mixing:
  rms(total) = 0.12854E+01    rms(broyden)= 0.12851E+01
  rms(prec ) = 0.13238E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9149
  1.0481  0.7817

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8530.04726537
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       325.17984352
  PAW double counting   =      5726.88981682    -5707.14615113
  entropy T*S    EENTRO =         0.01879009
  eigenvalues    EBANDS =     -1100.98457398
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -120.45381928 eV

  energy without entropy =     -120.47260937  energy(sigma->0) =     -120.46008264


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.1513649E+02  (-0.6523909E+00)
 number of electron      87.0000001 magnetization      16.8998539
 augmentation part        3.5642440 magnetization      15.7348888

 Broyden mixing:
  rms(total) = 0.82471E+00    rms(broyden)= 0.82464E+00
  rms(prec ) = 0.84864E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.1638
  1.7128  1.1217  0.6570

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8570.01856866
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       316.22442731
  PAW double counting   =      6628.74360242    -6608.87585474
  entropy T*S    EENTRO =         0.02120045
  eigenvalues    EBANDS =     -1067.32083370
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -135.59030617 eV

  energy without entropy =     -135.61150661  energy(sigma->0) =     -135.59737298


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




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



 eigenvalue-minimisations  :   384
 total energy-change (2. order) :-0.2916595E+02  (-0.1270059E+01)
 number of electron      87.0000000 magnetization      13.7521954
 augmentation part        3.5803874 magnetization      12.6474147

 Broyden mixing:
  rms(total) = 0.40927E+00    rms(broyden)= 0.40912E+00
  rms(prec ) = 0.41968E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2406
  2.2088  1.3137  0.7589  0.6808

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8606.85173336
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       295.63473093
  PAW double counting   =      7442.67634016    -7422.63253003
  entropy T*S    EENTRO =         0.02675549
  eigenvalues    EBANDS =     -1039.24553864
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -164.75625468 eV

  energy without entropy =     -164.78301017  energy(sigma->0) =     -164.76517318


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




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



 eigenvalue-minimisations  :   392
 total energy-change (2. order) :-0.1180811E+02  (-0.5675901E+00)
 number of electron      87.0000000 magnetization      10.1901269
 augmentation part        3.5534857 magnetization       9.1037224

 Broyden mixing:
  rms(total) = 0.31045E+00    rms(broyden)= 0.31033E+00
  rms(prec ) = 0.31950E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2366
  1.9324  1.9324  0.8480  0.8480  0.6222

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8621.52982885
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       287.87013879
  PAW double counting   =      7596.68431448    -7576.59228453
  entropy T*S    EENTRO =         0.03261687
  eigenvalues    EBANDS =     -1028.66504283
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -176.56436531 eV

  energy without entropy =     -176.59698217  energy(sigma->0) =     -176.57523759


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




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



 eigenvalue-minimisations  :   384
 total energy-change (2. order) :-0.7443507E+01  (-0.3849339E+00)
 number of electron      87.0000001 magnetization       6.9867331
 augmentation part        3.5244124 magnetization       5.9626138

 Broyden mixing:
  rms(total) = 0.21763E+00    rms(broyden)= 0.21745E+00
  rms(prec ) = 0.22461E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.3281
  2.3430  2.3430  0.9795  0.9795  0.7070  0.6163

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8623.75529834
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       282.13348394
  PAW double counting   =      7493.94667488    -7473.78044686
  entropy T*S    EENTRO =         0.02997629
  eigenvalues    EBANDS =     -1028.21798269
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -184.00787200 eV

  energy without entropy =     -184.03784829  energy(sigma->0) =     -184.01786410


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




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



 eigenvalue-minimisations  :   368
 total energy-change (2. order) :-0.3905416E+01  (-0.1697272E+00)
 number of electron      87.0000004 magnetization       5.0979774
 augmentation part        3.5239946 magnetization       4.1426225

 Broyden mixing:
  rms(total) = 0.14714E+00    rms(broyden)= 0.14707E+00
  rms(prec ) = 0.15111E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.3685
  3.1359  2.4460  1.2519  0.7800  0.6731  0.6462  0.6462

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.98900196
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       278.89528875
  PAW double counting   =      7425.73613416    -7405.49131747
  entropy T*S    EENTRO =         0.02995750
  eigenvalues    EBANDS =     -1027.73006978
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -187.91328803 eV

  energy without entropy =     -187.94324553  energy(sigma->0) =     -187.92327386


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.1334212E+01  (-0.4117379E-01)
 number of electron      87.0000001 magnetization       3.7351673
 augmentation part        3.5203724 magnetization       2.7625534

 Broyden mixing:
  rms(total) = 0.97869E-01    rms(broyden)= 0.97798E-01
  rms(prec ) = 0.10228E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.4345
  4.0563  2.3193  1.2980  0.7747  0.7747  0.6157  0.8530  0.7841

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.67435884
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.74012754
  PAW double counting   =      7393.17940664    -7372.92349998
  entropy T*S    EENTRO =         0.03275860
  eigenvalues    EBANDS =     -1026.23765430
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -189.24749956 eV

  energy without entropy =     -189.28025816  energy(sigma->0) =     -189.25841909


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




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



 eigenvalue-minimisations  :   392
 total energy-change (2. order) :-0.9114362E+00  (-0.1809668E-01)
 number of electron      87.0000003 magnetization       2.6956551
 augmentation part        3.5244948 magnetization       1.7381537

 Broyden mixing:
  rms(total) = 0.69831E-01    rms(broyden)= 0.69756E-01
  rms(prec ) = 0.74073E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.4747
  4.7762  2.1926  1.6746  1.0474  0.9083  0.7160  0.6149  0.6712  0.6712

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.79339755
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.82680008
  PAW double counting   =      7357.28482738    -7337.02594930
  entropy T*S    EENTRO =         0.03529354
  eigenvalues    EBANDS =     -1026.12223072
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -190.15893579 eV

  energy without entropy =     -190.19422933  energy(sigma->0) =     -190.17070030


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.6061863E+00  (-0.8134935E-02)
 number of electron      87.0000002 magnetization       1.9663977
 augmentation part        3.5242699 magnetization       1.0113092

 Broyden mixing:
  rms(total) = 0.42653E-01    rms(broyden)= 0.42602E-01
  rms(prec ) = 0.45070E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5030
  5.2976  2.0075  2.0075  1.1592  0.9777  0.9777  0.6917  0.6169  0.6472  0.6472

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8627.49958181
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.22292727
  PAW double counting   =      7340.13017036    -7319.88778343
  entropy T*S    EENTRO =         0.03500249
  eigenvalues    EBANDS =     -1025.40157770
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -190.76512206 eV

  energy without entropy =     -190.80012455  energy(sigma->0) =     -190.77678956


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




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



 eigenvalue-minimisations  :   368
 total energy-change (2. order) :-0.5489484E+00  (-0.3062953E-02)
 number of electron      86.9999999 magnetization       1.4613916
 augmentation part        3.5207704 magnetization       0.5056513

 Broyden mixing:
  rms(total) = 0.27146E-01    rms(broyden)= 0.27133E-01
  rms(prec ) = 0.28858E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5710
  5.8309  2.3421  2.3421  1.4416  1.0678  0.6514  0.6514  0.8213  0.8213  0.6921
  0.6194

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8627.45621217
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       275.64718850
  PAW double counting   =      7328.85344914    -7308.61838672
  entropy T*S    EENTRO =         0.03480283
  eigenvalues    EBANDS =     -1025.41063278
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -191.31407044 eV

  energy without entropy =     -191.34887326  energy(sigma->0) =     -191.32567138


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.4151448E+00  (-0.1387034E-02)
 number of electron      87.0000000 magnetization       1.2520696
 augmentation part        3.5184544 magnetization       0.3029057

 Broyden mixing:
  rms(total) = 0.16419E-01    rms(broyden)= 0.16402E-01
  rms(prec ) = 0.17854E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5928
  6.2626  2.5750  2.5750  1.4098  1.0571  1.0571  0.6279  0.6279  0.8205  0.8205
  0.6557  0.6242

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8627.41637964
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       275.23810945
  PAW double counting   =      7334.53052739    -7314.29859891
  entropy T*S    EENTRO =         0.03525597
  eigenvalues    EBANDS =     -1025.45385030
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -191.72921528 eV

  energy without entropy =     -191.76447124  energy(sigma->0) =     -191.74096726


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.1897411E+00  (-0.6221971E-03)
 number of electron      87.0000000 magnetization       1.1585055
 augmentation part        3.5204414 magnetization       0.2087929

 Broyden mixing:
  rms(total) = 0.10206E-01    rms(broyden)= 0.10190E-01
  rms(prec ) = 0.11065E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5703
  6.4800  2.8659  2.3975  1.5000  1.1499  1.1499  0.6312  0.6312  0.8184  0.8184
  0.6062  0.6829  0.6829

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.92192916
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       275.02251304
  PAW double counting   =      7336.37734996    -7316.14132619
  entropy T*S    EENTRO =         0.03508359
  eigenvalues    EBANDS =     -1025.92636841
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -191.91895641 eV

  energy without entropy =     -191.95404000  energy(sigma->0) =     -191.93065094


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.9642129E-01  ( 0.1034552E-02)
 number of electron      87.0000001 magnetization       1.1064743
 augmentation part        3.5212768 magnetization       0.1576533

 Broyden mixing:
  rms(total) = 0.66892E-02    rms(broyden)= 0.66841E-02
  rms(prec ) = 0.73957E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5723
  6.6601  3.3494  2.2098  1.7569  1.4055  0.8656  0.8656  0.8980  0.8980  0.6241
  0.6241  0.6710  0.6212  0.5630

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.56709623
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.90948405
  PAW double counting   =      7339.17162953    -7318.93120570
  entropy T*S    EENTRO =         0.03503198
  eigenvalues    EBANDS =     -1026.26894209
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.01537769 eV

  energy without entropy =     -192.05040968  energy(sigma->0) =     -192.02705502


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.7325479E-01  ( 0.1637904E-02)
 number of electron      87.0000000 magnetization       1.0856973
 augmentation part        3.5209334 magnetization       0.1372757

 Broyden mixing:
  rms(total) = 0.39404E-02    rms(broyden)= 0.39333E-02
  rms(prec ) = 0.45063E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5618
  6.7610  3.5872  2.0864  2.0864  1.5036  0.8782  0.8782  0.8475  0.8475  0.6402
  0.6402  0.7111  0.7111  0.6246  0.6246

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.49335105
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.83841246
  PAW double counting   =      7343.29161815    -7323.05040979
  entropy T*S    EENTRO =         0.03503046
  eigenvalues    EBANDS =     -1026.34565348
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.08863248 eV

  energy without entropy =     -192.12366295  energy(sigma->0) =     -192.10030930


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.2591472E-01  ( 0.1912349E-02)
 number of electron      87.0000000 magnetization       1.0783064
 augmentation part        3.5204713 magnetization       0.1295785

 Broyden mixing:
  rms(total) = 0.30302E-02    rms(broyden)= 0.30193E-02
  rms(prec ) = 0.34404E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.5172
  6.8143  3.6641  2.1545  2.1545  1.5313  0.7231  0.7231  0.8679  0.8679  0.8684
  0.8684  0.6276  0.6276  0.6202  0.6084  0.5542

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.45332627
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.81256595
  PAW double counting   =      7344.53190274    -7324.29110583
  entropy T*S    EENTRO =         0.03502726
  eigenvalues    EBANDS =     -1026.38533182
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.11454720 eV

  energy without entropy =     -192.14957446  energy(sigma->0) =     -192.12622295


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.5477434E-02  ( 0.2311450E-02)
 number of electron      87.0000000 magnetization       1.0729905
 augmentation part        3.5203710 magnetization       0.1241251

 Broyden mixing:
  rms(total) = 0.27006E-02    rms(broyden)= 0.26862E-02
  rms(prec ) = 0.30465E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.4918
  6.8767  3.8495  2.2720  2.2720  1.4948  0.9880  0.9880  0.8564  0.8564  0.6509
  0.6509  0.5700  0.5700  0.6492  0.6492  0.6242  0.5414

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.41175425
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80671677
  PAW double counting   =      7345.44396182    -7325.20300791
  entropy T*S    EENTRO =         0.03504116
  eigenvalues    EBANDS =     -1026.42670300
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.12002464 eV

  energy without entropy =     -192.15506579  energy(sigma->0) =     -192.13170502


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




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



 eigenvalue-minimisations  :   384
 total energy-change (2. order) :-0.3171197E-02  ( 0.2832057E-02)
 number of electron      87.0000000 magnetization       1.0722698
 augmentation part        3.5203358 magnetization       0.1231350

 Broyden mixing:
  rms(total) = 0.29013E-02    rms(broyden)= 0.28888E-02
  rms(prec ) = 0.31395E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.4317
  6.8857  3.8826  2.2976  2.2976  1.4807  1.0473  1.0473  0.8395  0.8395  0.6379
  0.6379  0.6821  0.6821  0.6257  0.4639  0.4639  0.5323  0.4264

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.39978732
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80425591
  PAW double counting   =      7345.60032996    -7325.35961096
  entropy T*S    EENTRO =         0.03506629
  eigenvalues    EBANDS =     -1026.43917047
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.12319583 eV

  energy without entropy =     -192.15826212  energy(sigma->0) =     -192.13488460


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) :-0.1386313E+02  (-0.1367320E+02)
 number of electron      87.0000000 magnetization       1.0722450
 augmentation part        3.4093207 magnetization       0.9052907

 Broyden mixing:
  rms(total) = 0.52826E+00    rms(broyden)= 0.52672E+00
  rms(prec ) = 0.55793E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.3563
  6.8856  3.8826  2.2975  2.2975  1.4800  1.0479  1.0479  0.8387  0.8387  0.6378
  0.6378  0.6824  0.6824  0.4623  0.4623  0.6258  0.5308  0.4314  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.39890494
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80381534
  PAW double counting   =      7345.44587739    -7325.20532872
  entropy T*S    EENTRO =        -0.02455941
  eigenvalues    EBANDS =     -1040.24294274
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -205.98632231 eV

  energy without entropy =     -205.96176290  energy(sigma->0) =     -205.97813584


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




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



 eigenvalue-minimisations  :   368
 total energy-change (2. order) : 0.1405844E+02  (-0.1797165E+02)
 number of electron      87.0000004 magnetization       1.0708095
 augmentation part        3.5124273 magnetization       0.1284787

 Broyden mixing:
  rms(total) = 0.40460E-01    rms(broyden)= 0.25046E-01
  rms(prec ) = 0.26023E-01
  weight for this iteration     100.00
  WARNING: grid for Broyden might be to small

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2892
  6.8880  3.8802  2.3050  2.3050  1.4669  1.0514  1.0514  0.8298  0.8298  0.6375
  0.6375  0.6908  0.6908  0.4633  0.4633  0.6261  0.5366  0.4240  0.0061  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.39683479
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.83398697
  PAW double counting   =      7345.45577317    -7325.21521880
  entropy T*S    EENTRO =         0.03511018
  eigenvalues    EBANDS =     -1026.27642436
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -191.92788687 eV

  energy without entropy =     -191.96299705  energy(sigma->0) =     -191.93959026


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.1864302E+00  (-0.1358243E+00)
 number of electron      87.0000000 magnetization       1.0723617
 augmentation part        3.5188878 magnetization       0.1238883

 Broyden mixing:
  rms(total) = 0.88947E-02    rms(broyden)= 0.55260E-02
  rms(prec ) = 0.57597E-02
  weight for this iteration     100.00
  WARNING: grid for Broyden might be to small

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2402
  6.8903  3.8950  2.3076  2.3076  1.4371  1.0762  1.0762  0.8223  0.8223  0.6321
  0.6321  0.6911  0.6911  0.6289  0.5132  0.5132  0.4654  0.4654  0.0886  0.0886
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.34125711
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80245893
  PAW double counting   =      7344.88178164    -7324.63984961
  entropy T*S    EENTRO =         0.03505155
  eigenvalues    EBANDS =     -1026.48822325
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.11431709 eV

  energy without entropy =     -192.14936864  energy(sigma->0) =     -192.12600094


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




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



 eigenvalue-minimisations  :   384
 total energy-change (2. order) :-0.1078813E-01  (-0.1183085E-02)
 number of electron      87.0000000 magnetization       1.0719819
 augmentation part        3.5198260 magnetization       0.1213373

 Broyden mixing:
  rms(total) = 0.40412E-02    rms(broyden)= 0.37116E-02
  rms(prec ) = 0.38708E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2413
  6.9111  3.9803  2.3352  2.3352  1.2867  1.2867  1.0892  0.7863  0.7863  0.8033
  0.8033  0.5935  0.5935  0.6027  0.6027  0.6008  0.6008  0.5096  0.5096  0.0000
  0.1463  0.1463

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.38222584
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80146074
  PAW double counting   =      7342.98601361    -7322.74535777
  entropy T*S    EENTRO =         0.03504496
  eigenvalues    EBANDS =     -1026.45576169
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.12510522 eV

  energy without entropy =     -192.16015018  energy(sigma->0) =     -192.13678687


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.2039778E-02  ( 0.3634414E-02)
 number of electron      87.0000000 magnetization       1.0718509
 augmentation part        3.5203357 magnetization       0.1206211

 Broyden mixing:
  rms(total) = 0.31162E-02    rms(broyden)= 0.30850E-02
  rms(prec ) = 0.32704E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.2111
  6.9146  3.9985  2.3545  2.3545  1.2825  1.2825  1.1460  0.7089  0.7089  0.8123
  0.8123  0.6692  0.6692  0.6460  0.6304  0.6203  0.6203  0.4952  0.4952  0.3529
  0.1400  0.1400  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.37160449
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80072566
  PAW double counting   =      7343.20314752    -7322.96224326
  entropy T*S    EENTRO =         0.03502993
  eigenvalues    EBANDS =     -1026.46792113
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.12714500 eV

  energy without entropy =     -192.16217493  energy(sigma->0) =     -192.13882164


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.4073567E-02  ( 0.9301484E-03)
 number of electron      87.0000000 magnetization       1.0716277
 augmentation part        3.5203546 magnetization       0.1205170

 Broyden mixing:
  rms(total) = 0.40907E-02    rms(broyden)= 0.40760E-02
  rms(prec ) = 0.42598E-02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.1727
  6.9196  4.0069  2.3699  2.3699  1.3257  1.3257  1.0189  0.9060  0.6871  0.6871
  0.7424  0.7424  0.7458  0.6242  0.6242  0.6525  0.6187  0.4994  0.4994  0.3479
  0.0000  0.1442  0.1442  0.1415

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.35854949
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.80024536
  PAW double counting   =      7343.48072068    -7323.23969209
  entropy T*S    EENTRO =         0.03503568
  eigenvalues    EBANDS =     -1026.48469948
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -192.13121857 eV

  energy without entropy =     -192.16625425  energy(sigma->0) =     -192.14289713


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) :-0.5889675E+05  (-0.5888835E+05)
 number of electron      87.0186072 magnetization       1.0709308
 augmentation part        2.7370887 magnetization       4.0670869

 Broyden mixing:
  rms(total) = 0.10356E+02    rms(broyden)= 0.10355E+02
  rms(prec ) = 0.10405E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.1257
  6.9196  4.0069  2.3699  2.3699  1.3252  1.3252  1.0170  0.9086  0.6869  0.6869
  0.7423  0.7423  0.7458  0.6241  0.6241  0.6529  0.6188  0.4993  0.4993  0.3481
  0.1442  0.1442  0.1415  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.35344325
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.79985242
  PAW double counting   =      7343.48697299    -7323.24592987
  entropy T*S    EENTRO =        -0.04482854
  eigenvalues    EBANDS =    -59923.16374992
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -59088.88540540 eV

  energy without entropy =   -59088.84057686  energy(sigma->0) =   -59088.87046255


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




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



 eigenvalue-minimisations  :   416
 total energy-change (2. order) : 0.3916787E+05  (-0.1995116E+05)
 number of electron      86.9834044 magnetization       1.0695659
 augmentation part        3.0297216 magnetization       7.0452887

 Broyden mixing:
  rms(total) = 0.14513E+02    rms(broyden)= 0.14507E+02
  rms(prec ) = 0.14562E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.0827
  6.9199  4.0084  2.3685  2.3685  1.3253  1.3253  1.0214  0.6882  0.6882  0.8961
  0.7447  0.7447  0.7552  0.6242  0.6242  0.6571  0.6135  0.4999  0.4999  0.3471
  0.1442  0.1442  0.1421  0.0000  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.33401731
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       275.05716788
  PAW double counting   =      7344.78712298    -7324.54602316
  entropy T*S    EENTRO =        -0.01530695
  eigenvalues    EBANDS =    -20755.60236045
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -19921.01769624 eV

  energy without entropy =   -19921.00238929  energy(sigma->0) =   -19921.01259392


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




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



 eigenvalue-minimisations  :   416
 total energy-change (2. order) : 0.1661267E+05  (-0.2552466E+04)
 number of electron      87.7845713 magnetization       1.0258813
 augmentation part        2.6249496 magnetization      15.5137134

 Broyden mixing:
  rms(total) = 0.69636E+01    rms(broyden)= 0.69592E+01
  rms(prec ) = 0.72418E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.0425
  6.9215  4.0071  2.3699  2.3699  1.3229  1.3229  1.0086  0.6799  0.6799  0.8935
  0.7500  0.7500  0.7607  0.6225  0.6225  0.6573  0.6141  0.5065  0.5065  0.3550
  0.1438  0.1438  0.1388  0.0001  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.35975223
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       275.51061278
  PAW double counting   =      7346.10818329    -7325.86681265
  entropy T*S    EENTRO =         0.00261144
  eigenvalues    EBANDS =     -4143.37717091
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -3308.34660752 eV

  energy without entropy =    -3308.34921896  energy(sigma->0) =    -3308.34747800


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




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



 eigenvalue-minimisations  :   496
 total energy-change (2. order) :-0.1198596E+04  (-0.2541383E+04)
 number of electron      87.8003259 magnetization       1.0267511
 augmentation part        2.7573931 magnetization      24.8093704

 Broyden mixing:
  rms(total) = 0.10631E+02    rms(broyden)= 0.10626E+02
  rms(prec ) = 0.10994E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   1.0056
  6.9216  4.0102  2.3722  2.3722  1.3239  1.3239  1.0134  0.9000  0.6757  0.6757
  0.7533  0.7533  0.7446  0.6616  0.6154  0.6218  0.6218  0.5053  0.5053  0.3585
  0.1439  0.1439  0.1393  0.0004  0.0004  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.57450941
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.20190754
  PAW double counting   =      7349.18684683    -7328.94458975
  entropy T*S    EENTRO =         0.01216545
  eigenvalues    EBANDS =     -5343.45969312
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -4506.94215170 eV

  energy without entropy =    -4506.95431715  energy(sigma->0) =    -4506.94620685


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




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



 eigenvalue-minimisations  :   496
 total energy-change (2. order) :-0.9309382E+04  (-0.1349251E+05)
 number of electron      87.0498665 magnetization       1.0274431
 augmentation part        2.0603727 magnetization      23.0065109

 Broyden mixing:
  rms(total) = 0.10665E+02    rms(broyden)= 0.10662E+02
  rms(prec ) = 0.11054E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9710
  6.9217  4.0109  2.3723  2.3723  1.3218  1.3218  1.0043  0.9112  0.6779  0.6779
  0.7531  0.7531  0.7427  0.6225  0.6225  0.6624  0.6144  0.5061  0.5061  0.3565
  0.1439  0.1439  0.1391  0.0004  0.0005  0.0005  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.65197611
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.21975276
  PAW double counting   =      7347.58163098    -7327.33795842
  entropy T*S    EENTRO =        -0.01242076
  eigenvalues    EBANDS =    -14653.75905388
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -13816.32430467 eV

  energy without entropy =   -13816.31188391  energy(sigma->0) =   -13816.32016442


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) : 0.7976608E+04  (-0.2356302E+04)
 number of electron      89.5038384 magnetization       1.0249192
 augmentation part        2.5962832 magnetization      17.4860650

 Broyden mixing:
  rms(total) = 0.86298E+01    rms(broyden)= 0.86287E+01
  rms(prec ) = 0.89059E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9387
  6.9225  4.0119  2.3738  2.3738  1.3218  1.3218  0.9796  0.9334  0.6723  0.6723
  0.7485  0.7485  0.7368  0.6720  0.6152  0.6219  0.6219  0.5170  0.5170  0.3522
  0.1437  0.1437  0.1367  0.0010  0.0010  0.0009  0.0007  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.38899132
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.29983060
  PAW double counting   =      7348.58700203    -7328.34133071
  entropy T*S    EENTRO =         0.01681140
  eigenvalues    EBANDS =     -6676.52551576
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -5839.71647300 eV

  energy without entropy =    -5839.73328440  energy(sigma->0) =    -5839.72207680


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) :-0.1670760E+04  (-0.5860210E+04)
 number of electron      89.2422890 magnetization       1.0184039
 augmentation part        2.4557385 magnetization      24.2969836

 Broyden mixing:
  rms(total) = 0.92122E+01    rms(broyden)= 0.92119E+01
  rms(prec ) = 0.95344E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9085
  6.9225  4.0118  2.3739  2.3739  1.3230  1.3230  0.9723  0.9395  0.6714  0.6714
  0.7487  0.7487  0.7359  0.6716  0.6153  0.6218  0.6218  0.5173  0.5173  0.3529
  0.1437  0.1437  0.1367  0.0011  0.0011  0.0008  0.0007  0.0005  0.0002  0.0002
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.33395288
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.68194861
  PAW double counting   =      7349.25377149    -7329.00718656
  entropy T*S    EENTRO =        -0.01320064
  eigenvalues    EBANDS =     -8348.69320009
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -7510.47609930 eV

  energy without entropy =    -7510.46289866  energy(sigma->0) =    -7510.47169909


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) :-0.6995325E+04  (-0.4718079E+04)
 number of electron      89.2555877 magnetization       1.0170536
 augmentation part        2.2789534 magnetization      26.0117633

 Broyden mixing:
  rms(total) = 0.10712E+02    rms(broyden)= 0.10712E+02
  rms(prec ) = 0.11077E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.8801
  6.9226  4.0117  2.3740  2.3740  1.3212  1.3212  0.9571  0.9571  0.6714  0.6714
  0.7486  0.7486  0.7379  0.6723  0.6151  0.6218  0.6218  0.5157  0.5157  0.3531
  0.1437  0.1437  0.1366  0.0013  0.0013  0.0010  0.0010  0.0009  0.0009  0.0002
  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.33374399
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.58252849
  PAW double counting   =      7350.09270051    -7329.84620084
  entropy T*S    EENTRO =        -0.03547984
  eigenvalues    EBANDS =    -15343.89698334
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -14505.80145824 eV

  energy without entropy =   -14505.76597839  energy(sigma->0) =   -14505.78963162


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) :-0.8166329E+05  (-0.9176932E+05)
 number of electron      88.9699127 magnetization       1.0142079
 augmentation part        1.8760563 magnetization      31.7379419

 Broyden mixing:
  rms(total) = 0.12980E+02    rms(broyden)= 0.12979E+02
  rms(prec ) = 0.13375E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.8534
  6.9236  4.0103  2.3732  2.3732  1.3201  1.3201  0.9510  0.9510  0.6705  0.6705
  0.7550  0.7550  0.7417  0.6704  0.6141  0.6243  0.6243  0.5132  0.5132  0.3467
  0.1445  0.1445  0.1430  0.0013  0.0013  0.0011  0.0011  0.0010  0.0010  0.0010
  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.33664175
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.60361872
  PAW double counting   =      7350.43913234    -7330.19260422
  entropy T*S    EENTRO =        -0.02654594
  eigenvalues    EBANDS =    -97007.21124287
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -96169.08856296 eV

  energy without entropy =   -96169.06201702  energy(sigma->0) =   -96169.07971431


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.5647303E+05  (-0.1497631E+06)
 number of electron      86.9864411 magnetization       1.0073937
 augmentation part        1.6179822 magnetization      35.8063946

 Broyden mixing:
  rms(total) = 0.16884E+02    rms(broyden)= 0.16883E+02
  rms(prec ) = 0.17360E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.8282
  6.9235  4.0107  2.3725  2.3725  1.3207  1.3207  0.9505  0.9505  0.6702  0.6702
  0.7544  0.7544  0.7453  0.6244  0.6244  0.6697  0.6132  0.5122  0.5122  0.3473
  0.1445  0.1445  0.1433  0.0011  0.0011  0.0012  0.0012  0.0010  0.0010  0.0010
  0.0002  0.0002  0.0000  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.29346562
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.23691010
  PAW double counting   =      7350.75548920    -7330.50854994
  entropy T*S    EENTRO =        -0.02931349
  eigenvalues    EBANDS =   -153479.91107183
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -152642.11428081 eV

  energy without entropy =  -152642.08496732  energy(sigma->0) =  -152642.10450964


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




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



 eigenvalue-minimisations  :   504
 total energy-change (2. order) :-0.9970883E+05  (-0.2499365E+06)
 number of electron      87.6509800 magnetization       1.0165901
 augmentation part        1.8904744 magnetization      46.4847796

 Broyden mixing:
  rms(total) = 0.28238E+02    rms(broyden)= 0.28235E+02
  rms(prec ) = 0.28575E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.8049
  6.9231  4.0109  2.3695  2.3695  1.3119  1.3119  1.0481  0.6750  0.6750  0.8480
  0.7580  0.7580  0.7697  0.6262  0.6262  0.6630  0.6118  0.5158  0.5158  0.3406
  0.1446  0.1446  0.1441  0.0013  0.0013  0.0011  0.0011  0.0011  0.0011  0.0009
  0.0009  0.0004  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.20301676
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.41772154
  PAW double counting   =      7349.03046462    -7328.78323430
  entropy T*S    EENTRO =         0.00477472
  eigenvalues    EBANDS =   -253188.04453362
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -252350.94210303 eV

  energy without entropy =  -252350.94687775  energy(sigma->0) =  -252350.94369460


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) :-0.1179533E+06  (-0.3702322E+06)
 number of electron      88.1396525 magnetization       1.0165153
 augmentation part        1.5348088 magnetization      45.6779238

 Broyden mixing:
  rms(total) = 0.41470E+02    rms(broyden)= 0.41470E+02
  rms(prec ) = 0.41681E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7826
  6.9232  4.0109  2.3698  2.3698  1.3122  1.3122  1.0484  0.6740  0.6740  0.8474
  0.7585  0.7585  0.7708  0.6263  0.6263  0.6617  0.6114  0.5162  0.5162  0.3423
  0.1445  0.1445  0.1430  0.0013  0.0013  0.0012  0.0012  0.0011  0.0011  0.0009
  0.0009  0.0004  0.0002  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.01322148
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       278.15167043
  PAW double counting   =      7345.47177020    -7325.22312746
  entropy T*S    EENTRO =        -0.00140420
  eigenvalues    EBANDS =   -371143.21950181
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -370304.19809355 eV

  energy without entropy =  -370304.19668935  energy(sigma->0) =  -370304.19762548


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) : 0.1273801E+06  (-0.2439078E+06)
 number of electron      88.8094444 magnetization       1.0234718
 augmentation part        1.8069698 magnetization      31.2042414

 Broyden mixing:
  rms(total) = 0.39582E+02    rms(broyden)= 0.39582E+02
  rms(prec ) = 0.39725E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7614
  6.9230  4.0119  2.3697  2.3697  1.3145  1.3145  1.0419  0.6740  0.6740  0.8562
  0.7575  0.7575  0.7674  0.6262  0.6262  0.6617  0.6105  0.5152  0.5152  0.3435
  0.1444  0.1444  0.1423  0.0013  0.0013  0.0011  0.0011  0.0011  0.0011  0.0009
  0.0009  0.0003  0.0003  0.0004  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.00162816
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.59834157
  PAW double counting   =      7346.36712830    -7326.11740125
  entropy T*S    EENTRO =        -0.02440928
  eigenvalues    EBANDS =   -243762.53748242
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -242924.07973047 eV

  energy without entropy =  -242924.05532119  energy(sigma->0) =  -242924.07159405


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




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



 eigenvalue-minimisations  :   424
 total energy-change (2. order) :-0.2498284E+06  (-0.3377935E+06)
 number of electron      88.2734059 magnetization       1.0233212
 augmentation part        1.8734648 magnetization      23.9123083

 Broyden mixing:
  rms(total) = 0.42352E+02    rms(broyden)= 0.42352E+02
  rms(prec ) = 0.42447E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7414
  6.9230  4.0118  2.3698  2.3698  1.3146  1.3146  1.0416  0.6740  0.6740  0.8560
  0.7575  0.7575  0.7678  0.6262  0.6262  0.6617  0.6105  0.5152  0.5152  0.3434
  0.1444  0.1444  0.1423  0.0013  0.0013  0.0011  0.0011  0.0011  0.0011  0.0009
  0.0009  0.0004  0.0003  0.0003  0.0002  0.0002  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.98861238
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.41507098
  PAW double counting   =      7346.53617226    -7326.28651567
  entropy T*S    EENTRO =        -0.01714156
  eigenvalues    EBANDS =   -493590.73606356
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -492752.44136916 eV

  energy without entropy =  -492752.42422760  energy(sigma->0) =  -492752.43565531


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




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



 eigenvalue-minimisations  :   448
 total energy-change (2. order) :-0.1504754E+07  (-0.1443891E+07)
 number of electron      87.6016652 magnetization       1.0210219
 augmentation part        1.4742603 magnetization      31.1252297

 Broyden mixing:
  rms(total) = 0.53890E+02    rms(broyden)= 0.53890E+02
  rms(prec ) = 0.54005E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7222
  6.9229  4.0112  2.3680  2.3680  1.3134  1.3134  1.0592  0.6732  0.6732  0.7587
  0.7587  0.8161  0.7911  0.6265  0.6265  0.6610  0.6106  0.5144  0.5144  0.3419
  0.1445  0.1445  0.1429  0.0013  0.0013  0.0012  0.0012  0.0010  0.0010  0.0009
  0.0009  0.0005  0.0005  0.0004  0.0003  0.0003  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.98588750
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.84530549
  PAW double counting   =      7346.72618918    -7326.47659201
  entropy T*S    EENTRO =        -0.02779941
  eigenvalues    EBANDS =  -1998344.00306333
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1997506.28612681 eV

  energy without entropy = -1997506.25832741  energy(sigma->0) = -1997506.27686034


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




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



 eigenvalue-minimisations  :   448
 total energy-change (2. order) :-0.2423676E+07  (-0.1520936E+07)
 number of electron      87.2640174 magnetization       1.0205661
 augmentation part        1.2326479 magnetization      27.8822841

 Broyden mixing:
  rms(total) = 0.54576E+02    rms(broyden)= 0.54576E+02
  rms(prec ) = 0.54699E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7044
  6.9228  4.0126  2.3691  2.3691  1.3127  1.3127  1.0642  0.6688  0.6688  0.7581
  0.7581  0.8081  0.8081  0.6249  0.6249  0.6589  0.6096  0.5157  0.5157  0.3457
  0.1447  0.1447  0.1449  0.0012  0.0012  0.0012  0.0012  0.0010  0.0010  0.0010
  0.0010  0.0004  0.0004  0.0004  0.0004  0.0004  0.0003  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.05809027
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.73696350
  PAW double counting   =      7347.24994696    -7327.00042433
  entropy T*S    EENTRO =        -0.01372577
  eigenvalues    EBANDS =  -4422019.83480353
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -4421182.28441268 eV

  energy without entropy = -4421182.27068691  energy(sigma->0) = -4421182.27983742


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




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



 eigenvalue-minimisations  :   480
 total energy-change (2. order) :-0.2420628E+06  (-0.3236910E+07)
 number of electron      87.2847323 magnetization       1.0220552
 augmentation part        0.7907386 magnetization      47.1953017

 Broyden mixing:
  rms(total) = 0.74870E+02    rms(broyden)= 0.74870E+02
  rms(prec ) = 0.75017E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6874
  6.9231  4.0120  2.3692  2.3692  1.3116  1.3116  1.0678  0.6727  0.6727  0.7581
  0.7581  0.8066  0.8066  0.6252  0.6252  0.6595  0.6103  0.5168  0.5168  0.3429
  0.1446  0.1446  0.1443  0.0013  0.0013  0.0011  0.0011  0.0010  0.0010  0.0010
  0.0010  0.0004  0.0004  0.0004  0.0004  0.0004  0.0004  0.0002  0.0002  0.0000
  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.95413404
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.23346245
  PAW double counting   =      7346.63018885    -7326.38024344
  entropy T*S    EENTRO =        -0.01484895
  eigenvalues    EBANDS =  -4664082.22726673
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -4663245.07712109 eV

  energy without entropy = -4663245.06227215  energy(sigma->0) = -4663245.07217144


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




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



 eigenvalue-minimisations  :   488
 total energy-change (2. order) : 0.7082052E+06  (-0.3861455E+07)
 number of electron      87.8626110 magnetization       1.0252651
 augmentation part        0.6029404 magnetization      44.2510737

 Broyden mixing:
  rms(total) = 0.82443E+02    rms(broyden)= 0.82443E+02
  rms(prec ) = 0.82565E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6711
  6.9232  4.0120  2.3692  2.3692  1.3147  1.3147  1.0646  0.6751  0.6751  0.7594
  0.7594  0.8113  0.8005  0.6251  0.6251  0.6602  0.6104  0.5146  0.5146  0.3416
  0.1447  0.1447  0.1453  0.0013  0.0013  0.0011  0.0011  0.0010  0.0010  0.0010
  0.0010  0.0004  0.0004  0.0004  0.0004  0.0004  0.0004  0.0001  0.0002  0.0002
  0.0000  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.96512105
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.48090801
  PAW double counting   =      7346.37431108    -7326.12442325
  entropy T*S    EENTRO =         0.00408024
  eigenvalues    EBANDS =  -3955877.24638215
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3955039.84090636 eV

  energy without entropy = -3955039.84498660  energy(sigma->0) = -3955039.84226644


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.2159477E+07  (-0.5960746E+07)
 number of electron      87.4669778 magnetization       1.0280071
 augmentation part        1.3005786 magnetization      25.3980051

 Broyden mixing:
  rms(total) = 0.10011E+03    rms(broyden)= 0.10011E+03
  rms(prec ) = 0.10018E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6556
  6.9238  4.0119  2.3714  2.3714  1.3215  1.3215  1.0498  0.6742  0.6742  0.7549
  0.7549  0.8191  0.7966  0.6263  0.6263  0.6589  0.6108  0.5192  0.5192  0.3349
  0.1448  0.1448  0.1466  0.0013  0.0013  0.0011  0.0011  0.0010  0.0010  0.0010
  0.0007  0.0007  0.0005  0.0005  0.0004  0.0004  0.0002  0.0002  0.0003  0.0002
  0.0002  0.0003  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.97458033
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.81820638
  PAW double counting   =      7345.68991172    -7325.43902695
  entropy T*S    EENTRO =        -0.04123136
  eigenvalues    EBANDS =  -6115354.63933196
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -6114516.95033174 eV

  energy without entropy = -6114516.90910038  energy(sigma->0) = -6114516.93658795


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




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



 eigenvalue-minimisations  :   424
 total energy-change (2. order) : 0.4641614E+07  (-0.8256241E+06)
 number of electron      88.3095059 magnetization       1.1542934
 augmentation part        1.3085936 magnetization      23.4921809

 Broyden mixing:
  rms(total) = 0.63277E+02    rms(broyden)= 0.63277E+02
  rms(prec ) = 0.63366E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6408
  6.9239  4.0124  2.3720  2.3720  1.3243  1.3243  1.0317  0.6766  0.6766  0.8475
  0.7520  0.7520  0.7821  0.6265  0.6265  0.6599  0.6104  0.5199  0.5199  0.3325
  0.1449  0.1449  0.1475  0.0013  0.0013  0.0011  0.0011  0.0009  0.0009  0.0009
  0.0009  0.0005  0.0005  0.0005  0.0005  0.0005  0.0005  0.0003  0.0003  0.0003
  0.0003  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.13459250
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.49365882
  PAW double counting   =      7347.01148637    -7326.76190548
  entropy T*S    EENTRO =        -0.03645014
  eigenvalues    EBANDS =  -1473740.43875947
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1472903.23084162 eV

  energy without entropy = -1472903.19439148  energy(sigma->0) = -1472903.21869157


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




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



 eigenvalue-minimisations  :   424
 total energy-change (2. order) :-0.8699356E+07  (-0.9590734E+07)
 number of electron      88.7989401 magnetization       1.0728992
 augmentation part        1.1432146 magnetization      28.0104817

 Broyden mixing:
  rms(total) = 0.65240E+02    rms(broyden)= 0.65240E+02
  rms(prec ) = 0.65329E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5466
  6.8440  3.9621  2.3030  1.1377  1.1377  0.6706  0.6706  0.6929  0.6929  0.6373
  0.6420  0.6420  0.5224  0.5224  0.3322  0.1382  0.1382  0.1602  0.0015  0.0015
  0.0013  0.0013  0.0013  0.0013  0.0012  0.0012  0.0008  0.0008  0.0006  0.0006
  0.0006  0.0006  0.0004  0.0004  0.0003  0.0003  0.0003  0.0000  0.0001  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8622.32357968
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.61299159
  PAW double counting   =      7313.65293251    -7293.36758169
  entropy T*S    EENTRO =         0.00753657
  eigenvalues    EBANDS = -10173099.21071387
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  = -10172258.99269380 eV

  energy without entropy =-10172259.00023037  energy(sigma->0) =-10172258.99520599


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




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



 eigenvalue-minimisations  :   448
 total energy-change (2. order) : 0.9381613E+07  (-0.7618652E+06)
 number of electron      89.3372918 magnetization       1.0701270
 augmentation part        1.2542964 magnetization      20.5704430

 Broyden mixing:
  rms(total) = 0.54423E+02    rms(broyden)= 0.54423E+02
  rms(prec ) = 0.54488E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5332
  6.8439  3.9622  2.3033  1.1362  1.1362  0.6703  0.6703  0.6924  0.6924  0.6374
  0.6423  0.6423  0.5220  0.5220  0.3319  0.1385  0.1385  0.1601  0.0016  0.0016
  0.0016  0.0016  0.0015  0.0015  0.0015  0.0015  0.0009  0.0009  0.0006  0.0006
  0.0006  0.0004  0.0004  0.0004  0.0004  0.0003  0.0003  0.0003  0.0003  0.0000
  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.46706437
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.17924422
  PAW double counting   =      7333.88215283    -7313.59870460
  entropy T*S    EENTRO =         0.00161986
  eigenvalues    EBANDS =   -791485.07811320
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -790646.44514449 eV

  energy without entropy =  -790646.44676435  energy(sigma->0) =  -790646.44568444


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.1172220E+07  (-0.1835683E+07)
 number of electron      88.8687931 magnetization       1.1391734
 augmentation part        1.1826017 magnetization      33.4956157

 Broyden mixing:
  rms(total) = 0.77487E+02    rms(broyden)= 0.77487E+02
  rms(prec ) = 0.77562E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5206
  6.8435  3.9647  2.3029  1.1394  1.1394  0.6658  0.6658  0.6923  0.6923  0.6360
  0.6391  0.6391  0.5246  0.5246  0.3336  0.1392  0.1392  0.1625  0.0017  0.0017
  0.0016  0.0016  0.0016  0.0016  0.0015  0.0015  0.0015  0.0015  0.0010  0.0010
  0.0006  0.0006  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0002  0.0002
  0.0003  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.49031205
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.41446260
  PAW double counting   =      7333.73430601    -7313.45183301
  entropy T*S    EENTRO =         0.01247510
  eigenvalues    EBANDS =  -1963704.93000249
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1962866.07518307 eV

  energy without entropy = -1962866.08765817  energy(sigma->0) = -1962866.07934143


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) : 0.1172034E+07  (-0.7586666E+06)
 number of electron      89.2399836 magnetization       1.2133420
 augmentation part        1.8136205 magnetization      14.6186542

 Broyden mixing:
  rms(total) = 0.50360E+02    rms(broyden)= 0.50359E+02
  rms(prec ) = 0.50421E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5085
  6.8434  3.9652  2.3032  1.1415  1.1415  0.6638  0.6638  0.6900  0.6900  0.6349
  0.6383  0.6383  0.5241  0.5241  0.3378  0.1387  0.1387  0.1604  0.0019  0.0019
  0.0018  0.0018  0.0019  0.0019  0.0017  0.0017  0.0016  0.0016  0.0014  0.0010
  0.0010  0.0006  0.0006  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0003
  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.53606393
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.17818247
  PAW double counting   =      7336.05283694    -7315.76100443
  entropy T*S    EENTRO =         0.00888556
  eigenvalues    EBANDS =   -791671.01487401
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -790832.43631664 eV

  energy without entropy =  -790832.44520220  energy(sigma->0) =  -790832.43927849


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.3047611E+06  (-0.5465346E+06)
 number of electron      88.7018543 magnetization       1.2109000
 augmentation part        1.3144336 magnetization      33.3117141

 Broyden mixing:
  rms(total) = 0.69793E+02    rms(broyden)= 0.69793E+02
  rms(prec ) = 0.69872E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4972
  6.8450  3.9652  2.3032  1.1420  1.1420  0.6652  0.6652  0.6346  0.6852  0.6852
  0.6421  0.6421  0.5266  0.5266  0.3352  0.1391  0.1391  0.1588  0.0036  0.0032
  0.0032  0.0025  0.0025  0.0025  0.0020  0.0020  0.0017  0.0017  0.0018  0.0018
  0.0011  0.0011  0.0006  0.0005  0.0005  0.0004  0.0004  0.0003  0.0003  0.0003
  0.0003  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.67918970
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.38581946
  PAW double counting   =      7333.50752356    -7313.19424110
  entropy T*S    EENTRO =        -0.00063510
  eigenvalues    EBANDS =  -1096431.16834823
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1095593.51335034 eV

  energy without entropy = -1095593.51271525  energy(sigma->0) = -1095593.51313865


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.6393761E+06  (-0.8649519E+06)
 number of electron      87.3785996 magnetization       1.0446288
 augmentation part        1.1340257 magnetization      20.4957447

 Broyden mixing:
  rms(total) = 0.52329E+02    rms(broyden)= 0.52329E+02
  rms(prec ) = 0.52430E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3191
  4.9354  2.1970  1.0795  0.8375  0.6217  0.6217  0.6129  0.5104  0.5104  0.3390
  0.1760  0.1446  0.1446  0.0040  0.0034  0.0034  0.0035  0.0023  0.0023  0.0015
  0.0015  0.0015  0.0015  0.0016  0.0014  0.0014  0.0007  0.0007  0.0006  0.0005
  0.0005  0.0004  0.0004  0.0003  0.0003  0.0003  0.0002  0.0002  0.0001  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.59844131
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.99515094
  PAW double counting   =      7332.57802930    -7312.26266951
  entropy T*S    EENTRO =        -0.00515454
  eigenvalues    EBANDS =  -1735807.00049034
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1734969.65785470 eV

  energy without entropy = -1734969.65270016  energy(sigma->0) = -1734969.65613652


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) :-0.2198261E+07  (-0.3194263E+06)
 number of electron      87.6714375 magnetization       1.0662044
 augmentation part        1.1619078 magnetization      25.0120894

 Broyden mixing:
  rms(total) = 0.53309E+02    rms(broyden)= 0.53309E+02
  rms(prec ) = 0.53428E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3117
  4.9361  2.1969  1.0759  0.8358  0.6240  0.6240  0.6130  0.5101  0.5101  0.3415
  0.1761  0.1442  0.1442  0.0052  0.0041  0.0041  0.0034  0.0035  0.0035  0.0029
  0.0029  0.0027  0.0023  0.0023  0.0016  0.0014  0.0010  0.0010  0.0006  0.0005
  0.0005  0.0004  0.0004  0.0003  0.0002  0.0003  0.0003  0.0003  0.0000  0.0001
  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.39613139
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.21370483
  PAW double counting   =      7335.83063777    -7315.48417942
  entropy T*S    EENTRO =        -0.03782420
  eigenvalues    EBANDS =  -3934066.56029408
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3933230.79836574 eV

  energy without entropy = -3933230.76054154  energy(sigma->0) = -3933230.78575768


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) : 0.2277868E+07  (-0.9265083E+06)
 number of electron      88.0810398 magnetization       1.0761337
 augmentation part        1.0263481 magnetization      26.9921684

 Broyden mixing:
  rms(total) = 0.55100E+02    rms(broyden)= 0.55100E+02
  rms(prec ) = 0.55217E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3046
  4.9370  2.1976  1.0775  0.8347  0.6240  0.6240  0.6140  0.5119  0.5119  0.3417
  0.1737  0.1443  0.1443  0.0059  0.0059  0.0062  0.0048  0.0048  0.0040  0.0034
  0.0024  0.0024  0.0020  0.0020  0.0016  0.0016  0.0016  0.0016  0.0010  0.0010
  0.0005  0.0005  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0003  0.0002
  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8626.31617130
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.37620879
  PAW double counting   =      7343.21693378    -7322.86894236
  entropy T*S    EENTRO =         0.02425980
  eigenvalues    EBANDS =  -1656198.53000339
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1655362.46199391 eV

  energy without entropy = -1655362.48625371  energy(sigma->0) = -1655362.47008051


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




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



 eigenvalue-minimisations  :   440
 total energy-change (2. order) : 0.3675954E+06  (-0.1724037E+06)
 number of electron      88.2609914 magnetization       1.0777169
 augmentation part        1.0828287 magnetization      27.1484965

 Broyden mixing:
  rms(total) = 0.53211E+02    rms(broyden)= 0.53211E+02
  rms(prec ) = 0.53326E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2979
  4.9370  2.1992  1.0803  0.8295  0.6246  0.6246  0.6140  0.5120  0.5120  0.3435
  0.1724  0.1444  0.1444  0.0088  0.0080  0.0080  0.0071  0.0071  0.0034  0.0030
  0.0030  0.0030  0.0025  0.0025  0.0023  0.0023  0.0017  0.0017  0.0014  0.0010
  0.0010  0.0006  0.0006  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0003
  0.0000  0.0002  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8624.89801751
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.60606308
  PAW double counting   =      7338.44292200    -7318.08865981
  entropy T*S    EENTRO =         0.02420277
  eigenvalues    EBANDS =  -1288604.79821887
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1287767.07598757 eV

  energy without entropy = -1287767.10019034  energy(sigma->0) = -1287767.08405516


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




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



 eigenvalue-minimisations  :   440
 total energy-change (2. order) :-0.7842092E+07  (-0.3507071E+06)
 number of electron      87.1717529 magnetization       1.0430729
 augmentation part        1.1433831 magnetization      24.0184851

 Broyden mixing:
  rms(total) = 0.44974E+02    rms(broyden)= 0.44974E+02
  rms(prec ) = 0.45118E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2916
  4.9439  2.1992  1.0864  0.8310  0.6247  0.6247  0.6145  0.5131  0.5131  0.3387
  0.1725  0.1452  0.1452  0.0086  0.0086  0.0089  0.0089  0.0075  0.0065  0.0037
  0.0037  0.0031  0.0031  0.0025  0.0017  0.0017  0.0015  0.0015  0.0015  0.0014
  0.0009  0.0009  0.0007  0.0005  0.0005  0.0005  0.0004  0.0004  0.0003  0.0003
  0.0003  0.0000  0.0001  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8625.10785293
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.66111726
  PAW double counting   =      7337.94569704    -7317.59330606
  entropy T*S    EENTRO =         0.00619321
  eigenvalues    EBANDS =  -9130696.34195354
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -9129858.79438426 eV

  energy without entropy = -9129858.80057746  energy(sigma->0) = -9129858.79644866


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




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



 eigenvalue-minimisations  :   440
 total energy-change (2. order) : 0.2786367E+07  (-0.4716028E+06)
 number of electron      87.4479782 magnetization       1.1615112
 augmentation part        1.2494407 magnetization      15.1504694

 Broyden mixing:
  rms(total) = 0.38419E+02    rms(broyden)= 0.38419E+02
  rms(prec ) = 0.38567E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0952
  1.3650  0.5515  0.5515  0.4926  0.3245  0.1504  0.1504  0.1170  0.0111  0.0111
  0.0112  0.0112  0.0081  0.0081  0.0071  0.0071  0.0056  0.0031  0.0023  0.0023
  0.0021  0.0021  0.0018  0.0016  0.0014  0.0014  0.0010  0.0010  0.0005  0.0005
  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0000  0.0001  0.0004  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8628.07858177
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.52825012
  PAW double counting   =      7374.16093966    -7353.79844689
  entropy T*S    EENTRO =         0.03091405
  eigenvalues    EBANDS =  -6344325.77501606
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -6343491.29622014 eV

  energy without entropy = -6343491.32713419  energy(sigma->0) = -6343491.30652482


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




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



 eigenvalue-minimisations  :   448
 total energy-change (2. order) : 0.1374632E+07  (-0.1639442E+07)
 number of electron      86.7895560 magnetization       1.1399713
 augmentation part        0.7798385 magnetization      41.2151824

 Broyden mixing:
  rms(total) = 0.47992E+02    rms(broyden)= 0.47991E+02
  rms(prec ) = 0.48229E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0936
  1.3651  0.5508  0.5508  0.4801  0.3319  0.1526  0.1526  0.1118  0.0194  0.0194
  0.0176  0.0132  0.0132  0.0094  0.0094  0.0062  0.0054  0.0037  0.0032  0.0032
  0.0022  0.0022  0.0018  0.0018  0.0018  0.0016  0.0009  0.0009  0.0008  0.0006
  0.0006  0.0006  0.0003  0.0003  0.0004  0.0004  0.0003  0.0003  0.0000  0.0002
  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8628.27624140
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.52892972
  PAW double counting   =      7379.32013222    -7358.96998384
  entropy T*S    EENTRO =        -0.01369819
  eigenvalues    EBANDS =  -4969693.26818073
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -4968859.04332146 eV

  energy without entropy = -4968859.02962326  energy(sigma->0) = -4968859.03875539


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




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



 eigenvalue-minimisations  :   480
 total energy-change (2. order) : 0.2927833E+07  (-0.1347232E+07)
 number of electron      87.8831515 magnetization       1.1253072
 augmentation part        0.3793223 magnetization      53.8810708

 Broyden mixing:
  rms(total) = 0.63255E+02    rms(broyden)= 0.63254E+02
  rms(prec ) = 0.63466E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0919
  1.3640  0.5539  0.5539  0.4457  0.3603  0.1528  0.1528  0.1100  0.0200  0.0200
  0.0152  0.0152  0.0150  0.0150  0.0147  0.0147  0.0062  0.0036  0.0031  0.0023
  0.0023  0.0021  0.0021  0.0020  0.0020  0.0015  0.0015  0.0014  0.0011  0.0011
  0.0005  0.0005  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0002  0.0002
  0.0003  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8628.54936409
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.41342200
  PAW double counting   =      7388.30281310    -7367.96716922
  entropy T*S    EENTRO =         0.02844220
  eigenvalues    EBANDS =  -2041860.29292913
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -2041026.42906439 eV

  energy without entropy = -2041026.45750658  energy(sigma->0) = -2041026.43854512


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.5037820E+06  (-0.1641092E+07)
 number of electron      87.5721792 magnetization       1.1249921
 augmentation part        1.1235697 magnetization      24.6622343

 Broyden mixing:
  rms(total) = 0.55854E+02    rms(broyden)= 0.55854E+02
  rms(prec ) = 0.55983E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0900
  1.3639  0.5533  0.5533  0.4504  0.3557  0.1519  0.1519  0.1121  0.0176  0.0176
  0.0167  0.0167  0.0150  0.0150  0.0151  0.0151  0.0102  0.0062  0.0033  0.0030
  0.0030  0.0026  0.0026  0.0024  0.0023  0.0023  0.0015  0.0013  0.0012  0.0010
  0.0010  0.0006  0.0006  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003  0.0003
  0.0002  0.0002  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8627.83722376
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.99701313
  PAW double counting   =      7383.27442542    -7362.92166355
  entropy T*S    EENTRO =        -0.00988874
  eigenvalues    EBANDS =  -2545643.53535124
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -2544808.39696798 eV

  energy without entropy = -2544808.38707924  energy(sigma->0) = -2544808.39367173


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) : 0.6411407E+06  (-0.1510251E+07)
 number of electron      87.6832435 magnetization       1.1683911
 augmentation part        0.8948316 magnetization      42.0140573

 Broyden mixing:
  rms(total) = 0.66512E+02    rms(broyden)= 0.66512E+02
  rms(prec ) = 0.66673E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0888
  1.3532  0.5427  0.5427  0.5283  0.3265  0.1523  0.1523  0.1104  0.0187  0.0187
  0.0175  0.0175  0.0148  0.0148  0.0148  0.0148  0.0139  0.0139  0.0060  0.0043
  0.0042  0.0026  0.0026  0.0019  0.0019  0.0018  0.0018  0.0019  0.0019  0.0015
  0.0011  0.0011  0.0005  0.0005  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003
  0.0000  0.0002  0.0002  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8628.62341109
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.77185620
  PAW double counting   =      7391.91261525    -7371.56521640
  entropy T*S    EENTRO =        -0.02128880
  eigenvalues    EBANDS =  -1904501.76738671
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1903667.65711077 eV

  energy without entropy = -1903667.63582197  energy(sigma->0) = -1903667.65001450


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




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



 eigenvalue-minimisations  :   400
 total energy-change (2. order) :-0.1844966E+07  (-0.6494092E+06)
 number of electron      87.6062210 magnetization       1.2449467
 augmentation part        2.0094446 magnetization       8.1014692

 Broyden mixing:
  rms(total) = 0.30541E+02    rms(broyden)= 0.30540E+02
  rms(prec ) = 0.30657E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0309
  0.3929  0.3929  0.1090  0.1090  0.0310  0.0228  0.0228  0.0174  0.0174  0.0172
  0.0172  0.0158  0.0158  0.0101  0.0101  0.0064  0.0045  0.0034  0.0034  0.0026
  0.0020  0.0020  0.0015  0.0015  0.0014  0.0014  0.0011  0.0011  0.0006  0.0005
  0.0005  0.0004  0.0003  0.0003  0.0003  0.0003  0.0002  0.0002  0.0001  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8628.77212079
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.14536484
  PAW double counting   =      7395.00195178    -7374.65692788
  entropy T*S    EENTRO =        -0.03574825
  eigenvalues    EBANDS =  -3749468.43718898
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3748634.11894852 eV

  energy without entropy = -3748634.08320027  energy(sigma->0) = -3748634.10703244


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




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



 eigenvalue-minimisations  :   448
 total energy-change (2. order) :-0.4554157E+07  (-0.3026648E+06)
 number of electron      86.7068978 magnetization       1.2624029
 augmentation part        0.9505512 magnetization      30.3402662

 Broyden mixing:
  rms(total) = 0.53411E+02    rms(broyden)= 0.53411E+02
  rms(prec ) = 0.53587E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0312
  0.3934  0.3934  0.1090  0.1090  0.0320  0.0320  0.0170  0.0209  0.0209  0.0189
  0.0189  0.0181  0.0181  0.0159  0.0159  0.0064  0.0059  0.0059  0.0045  0.0026
  0.0026  0.0027  0.0025  0.0025  0.0019  0.0019  0.0010  0.0008  0.0008  0.0008
  0.0006  0.0006  0.0005  0.0002  0.0001  0.0003  0.0003  0.0003  0.0003  0.0001
  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8622.14969891
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.66406173
  PAW double counting   =      7396.19181777    -7375.93875004
  entropy T*S    EENTRO =        -0.04680638
  eigenvalues    EBANDS =  -8303630.97833841
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -8302790.62199347 eV

  energy without entropy = -8302790.57518709  energy(sigma->0) = -8302790.60639134


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) : 0.4981315E+07  (-0.2811808E+07)
 number of electron      88.0222399 magnetization       1.2654751
 augmentation part        0.9464015 magnetization      21.6855612

 Broyden mixing:
  rms(total) = 0.56687E+02    rms(broyden)= 0.56687E+02
  rms(prec ) = 0.56811E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0316
  0.3917  0.3917  0.1091  0.1091  0.0408  0.0408  0.0235  0.0235  0.0196  0.0196
  0.0200  0.0200  0.0168  0.0168  0.0160  0.0160  0.0125  0.0052  0.0060  0.0047
  0.0047  0.0021  0.0021  0.0019  0.0017  0.0017  0.0014  0.0014  0.0012  0.0009
  0.0009  0.0006  0.0006  0.0006  0.0002  0.0003  0.0003  0.0003  0.0003  0.0003
  0.0001  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8621.03900966
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.58435959
  PAW double counting   =      7387.33423890    -7367.08984502
  entropy T*S    EENTRO =        -0.01697248
  eigenvalues    EBANDS =  -3322316.83991134
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3321475.43141924 eV

  energy without entropy = -3321475.41444676  energy(sigma->0) = -3321475.42576175


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




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



 eigenvalue-minimisations  :   496
 total energy-change (2. order) :-0.6765139E+06  (-0.3963921E+07)
 number of electron      88.4267571 magnetization       1.2651742
 augmentation part        0.1140304 magnetization      50.0497853

 Broyden mixing:
  rms(total) = 0.80552E+02    rms(broyden)= 0.80552E+02
  rms(prec ) = 0.80688E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0322
  0.3928  0.3928  0.1088  0.1088  0.0453  0.0453  0.0316  0.0316  0.0287  0.0287
  0.0216  0.0202  0.0202  0.0158  0.0158  0.0153  0.0153  0.0109  0.0063  0.0038
  0.0026  0.0026  0.0024  0.0024  0.0025  0.0025  0.0019  0.0019  0.0013  0.0010
  0.0010  0.0005  0.0005  0.0005  0.0005  0.0004  0.0004  0.0003  0.0003  0.0002
  0.0002  0.0001  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8621.57727152
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.12539089
  PAW double counting   =      7399.79956086    -7379.57406309
  entropy T*S    EENTRO =        -0.07717864
  eigenvalues    EBANDS =  -3998830.63167734
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3997989.29951806 eV

  energy without entropy = -3997989.22233943  energy(sigma->0) = -3997989.27379185


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) : 0.1069677E+07  (-0.2803737E+07)
 number of electron      88.7468380 magnetization       1.2627148
 augmentation part        0.9343885 magnetization      16.7179370

 Broyden mixing:
  rms(total) = 0.58381E+02    rms(broyden)= 0.58380E+02
  rms(prec ) = 0.58460E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0313
  0.3944  0.3944  0.1087  0.1087  0.0453  0.0453  0.0305  0.0305  0.0270  0.0270
  0.0199  0.0199  0.0155  0.0155  0.0154  0.0154  0.0099  0.0099  0.0091  0.0063
  0.0027  0.0027  0.0026  0.0026  0.0025  0.0025  0.0020  0.0020  0.0013  0.0010
  0.0010  0.0006  0.0006  0.0006  0.0006  0.0003  0.0004  0.0004  0.0003  0.0003
  0.0000  0.0001  0.0001  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8621.58488590
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.18856745
  PAW double counting   =      7399.77609863    -7379.55057006
  entropy T*S    EENTRO =         0.04199900
  eigenvalues    EBANDS =  -2929153.64690313
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -2928312.13997324 eV

  energy without entropy = -2928312.18197224  energy(sigma->0) = -2928312.15397291


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.6587579E+06  (-0.2730202E+07)
 number of electron      87.9686900 magnetization       1.1476730
 augmentation part        0.4561342 magnetization      38.9871301

 Broyden mixing:
  rms(total) = 0.96501E+02    rms(broyden)= 0.96501E+02
  rms(prec ) = 0.96595E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0102
  0.0486  0.0486  0.0398  0.0398  0.0312  0.0312  0.0206  0.0211  0.0211  0.0155
  0.0155  0.0150  0.0150  0.0060  0.0049  0.0049  0.0038  0.0038  0.0023  0.0023
  0.0021  0.0021  0.0017  0.0017  0.0015  0.0015  0.0009  0.0009  0.0008  0.0008
  0.0005  0.0006  0.0004  0.0004  0.0004  0.0004  0.0001  0.0001  0.0003  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8621.36572910
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       277.13440526
  PAW double counting   =      7401.12667661    -7380.89451773
  entropy T*S    EENTRO =         0.00650527
  eigenvalues    EBANDS =  -3587911.73193585
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3587070.08887476 eV

  energy without entropy = -3587070.09538004  energy(sigma->0) = -3587070.09104319


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.3769664E+07  (-0.2153136E+07)
 number of electron      87.6500596 magnetization       1.1791582
 augmentation part        1.7814299 magnetization      11.8976445

 Broyden mixing:
  rms(total) = 0.62801E+02    rms(broyden)= 0.62800E+02
  rms(prec ) = 0.62854E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0153
  0.2131  0.0527  0.0527  0.0387  0.0387  0.0311  0.0311  0.0198  0.0213  0.0213
  0.0156  0.0156  0.0150  0.0150  0.0060  0.0046  0.0046  0.0038  0.0038  0.0024
  0.0024  0.0021  0.0021  0.0017  0.0017  0.0015  0.0015  0.0009  0.0009  0.0010
  0.0010  0.0005  0.0006  0.0004  0.0004  0.0004  0.0004  0.0001  0.0001  0.0003
  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8588.65183069
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       283.37643574
  PAW double counting   =      7306.24803172    -7285.57844354
  entropy T*S    EENTRO =         0.00080010
  eigenvalues    EBANDS =  -7357615.12952759
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -7356734.09881349 eV

  energy without entropy = -7356734.09961359  energy(sigma->0) = -7356734.09908019


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




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



 eigenvalue-minimisations  :   480
 total energy-change (2. order) : 0.1458600E+07  (-0.3970129E+07)
 number of electron      86.6179641 magnetization       1.1846069
 augmentation part        0.8891903 magnetization      41.3620130

 Broyden mixing:
  rms(total) = 0.91431E+02    rms(broyden)= 0.91430E+02
  rms(prec ) = 0.91546E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0161
  0.1533  0.1077  0.0533  0.0533  0.0385  0.0385  0.0311  0.0311  0.0198  0.0218
  0.0218  0.0151  0.0151  0.0150  0.0150  0.0060  0.0051  0.0043  0.0037  0.0037
  0.0026  0.0026  0.0021  0.0021  0.0017  0.0017  0.0015  0.0015  0.0009  0.0009
  0.0010  0.0010  0.0005  0.0006  0.0004  0.0004  0.0004  0.0004  0.0001  0.0001
  0.0003  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8577.08267828
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       282.45620009
  PAW double counting   =      7324.14419201    -7303.41889343
  entropy T*S    EENTRO =        -0.00913116
  eigenvalues    EBANDS =  -5899025.81160827
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -5898134.08619827 eV

  energy without entropy = -5898134.07706712  energy(sigma->0) = -5898134.08315455


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




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



 eigenvalue-minimisations  :   496
 total energy-change (2. order) : 0.2497572E+07  (-0.2056314E+07)
 number of electron      87.5599593 magnetization       1.1415337
 augmentation part        1.0616216 magnetization      25.1453528

 Broyden mixing:
  rms(total) = 0.97125E+02    rms(broyden)= 0.97125E+02
  rms(prec ) = 0.97196E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0179
  0.1866  0.1133  0.0558  0.0558  0.0551  0.0385  0.0385  0.0309  0.0309  0.0212
  0.0212  0.0163  0.0154  0.0154  0.0148  0.0148  0.0061  0.0048  0.0048  0.0039
  0.0039  0.0024  0.0024  0.0021  0.0021  0.0017  0.0017  0.0015  0.0015  0.0009
  0.0009  0.0009  0.0009  0.0006  0.0005  0.0004  0.0004  0.0004  0.0004  0.0001
  0.0001  0.0003  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8571.82613974
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       282.03167891
  PAW double counting   =      7282.51551309    -7261.84531087
  entropy T*S    EENTRO =         0.02730566
  eigenvalues    EBANDS =  -3401458.97396632
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -3400562.43519851 eV

  energy without entropy = -3400562.46250417  energy(sigma->0) = -3400562.44430039


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




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



 eigenvalue-minimisations  :   504
 total energy-change (2. order) :-0.1584830E+07  (-0.2503545E+07)
 number of electron      86.9189209 magnetization       1.1277995
 augmentation part        0.1501826 magnetization      60.4558498

 Broyden mixing:
  rms(total) = 0.15954E+03    rms(broyden)= 0.15954E+03
  rms(prec ) = 0.15963E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0203
  0.1663  0.1663  0.0733  0.0733  0.0525  0.0525  0.0393  0.0393  0.0310  0.0310
  0.0219  0.0211  0.0211  0.0153  0.0153  0.0147  0.0147  0.0060  0.0046  0.0046
  0.0036  0.0036  0.0024  0.0024  0.0021  0.0021  0.0017  0.0017  0.0015  0.0015
  0.0009  0.0009  0.0009  0.0009  0.0006  0.0004  0.0004  0.0004  0.0004  0.0004
  0.0001  0.0001  0.0003  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8567.98251959
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       281.52755053
  PAW double counting   =      7260.32002863    -7239.69426939
  entropy T*S    EENTRO =         0.00469862
  eigenvalues    EBANDS =  -4986292.10139854
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -4985392.29018898 eV

  energy without entropy = -4985392.29488760  energy(sigma->0) = -4985392.29175518


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




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



 eigenvalue-minimisations  :   488
 total energy-change (2. order) :-0.1021072E+07  (-0.6000056E+07)
 number of electron      86.9483105 magnetization       1.0053747
 augmentation part        0.2435265 magnetization      58.4138980

 Broyden mixing:
  rms(total) = 0.19264E+03    rms(broyden)= 0.19264E+03
  rms(prec ) = 0.19272E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0196
  0.1650  0.1650  0.0521  0.0521  0.0414  0.0440  0.0440  0.0244  0.0197  0.0197
  0.0168  0.0168  0.0167  0.0167  0.0136  0.0126  0.0126  0.0089  0.0080  0.0080
  0.0046  0.0046  0.0028  0.0024  0.0021  0.0021  0.0014  0.0008  0.0008  0.0008
  0.0006  0.0006  0.0005  0.0004  0.0002  0.0002  0.0002  0.0002  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8582.22781686
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       280.69108014
  PAW double counting   =      7294.12257677    -7273.59465756
  entropy T*S    EENTRO =         0.00622912
  eigenvalues    EBANDS =  -6007348.87578769
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -6006464.24265531 eV

  energy without entropy = -6006464.24888443  energy(sigma->0) = -6006464.24473169


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




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



 eigenvalue-minimisations  :   480
 total energy-change (2. order) :-0.1134712E+07  (-0.7130017E+07)
 number of electron      87.3013230 magnetization       1.1210446
 augmentation part        0.8157848 magnetization      58.5246424

 Broyden mixing:
  rms(total) = 0.25368E+03    rms(broyden)= 0.25368E+03
  rms(prec ) = 0.25373E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0219
  0.1661  0.1661  0.1114  0.0547  0.0512  0.0512  0.0401  0.0401  0.0262  0.0171
  0.0171  0.0168  0.0168  0.0167  0.0167  0.0139  0.0128  0.0128  0.0100  0.0078
  0.0078  0.0043  0.0043  0.0032  0.0024  0.0023  0.0023  0.0015  0.0008  0.0008
  0.0007  0.0007  0.0006  0.0005  0.0005  0.0002  0.0002  0.0002  0.0002  0.0001
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8543.24051530
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       294.34191703
  PAW double counting   =      7670.51858480    -7627.37553595
  entropy T*S    EENTRO =         0.02946513
  eigenvalues    EBANDS =  -7142135.94381847
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -7141176.03418198 eV

  energy without entropy = -7141176.06364712  energy(sigma->0) = -7141176.04400370


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




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



 eigenvalue-minimisations  :   496
 total energy-change (2. order) : 0.4326378E+07  (-0.2792225E+07)
 number of electron      87.4481454 magnetization       1.1472190
 augmentation part        0.8542773 magnetization      47.5994815

 Broyden mixing:
  rms(total) = 0.33604E+03    rms(broyden)= 0.33604E+03
  rms(prec ) = 0.33608E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0223
  0.1660  0.1660  0.1047  0.0637  0.0637  0.0518  0.0434  0.0434  0.0212  0.0212
  0.0195  0.0172  0.0172  0.0155  0.0155  0.0144  0.0144  0.0108  0.0108  0.0103
  0.0103  0.0099  0.0048  0.0048  0.0022  0.0022  0.0019  0.0015  0.0015  0.0009
  0.0009  0.0006  0.0004  0.0004  0.0004  0.0004  0.0003  0.0003  0.0002  0.0002
  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8547.19820646
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       293.66236343
  PAW double counting   =      7786.34045017    -7765.00583576
  entropy T*S    EENTRO =        -0.04049397
  eigenvalues    EBANDS =  -2815731.61380845
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -2814798.21981027 eV

  energy without entropy = -2814798.17931630  energy(sigma->0) = -2814798.20631228


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




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



 eigenvalue-minimisations  :   504
 total energy-change (2. order) : 0.1810523E+07  (-0.9503342E+06)
 number of electron      88.0663995 magnetization       1.1457981
 augmentation part        1.6014528 magnetization      28.8967663

 Broyden mixing:
  rms(total) = 0.31103E+03    rms(broyden)= 0.31103E+03
  rms(prec ) = 0.31105E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0228
  0.1779  0.1779  0.1053  0.0584  0.0584  0.0473  0.0410  0.0410  0.0257  0.0257
  0.0254  0.0177  0.0177  0.0172  0.0172  0.0174  0.0174  0.0144  0.0144  0.0115
  0.0088  0.0088  0.0082  0.0051  0.0051  0.0021  0.0021  0.0017  0.0017  0.0014
  0.0012  0.0012  0.0006  0.0006  0.0006  0.0004  0.0004  0.0002  0.0002  0.0002
  0.0002  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8550.07024078
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       293.56246023
  PAW double counting   =      7783.13591569    -7761.77422522
  entropy T*S    EENTRO =        -0.01851536
  eigenvalues    EBANDS =  -1005205.20932121
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1004274.73820589 eV

  energy without entropy = -1004274.71969053  energy(sigma->0) = -1004274.73203410


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) : 0.1937921E+05  (-0.9884912E+06)
 number of electron      88.2342179 magnetization       1.0048660
 augmentation part        2.0948460 magnetization      19.5482181

 Broyden mixing:
  rms(total) = 0.25009E+03    rms(broyden)= 0.25009E+03
  rms(prec ) = 0.25011E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0240
  0.1822  0.1822  0.0985  0.0560  0.0560  0.0481  0.0472  0.0472  0.0414  0.0414
  0.0270  0.0270  0.0209  0.0202  0.0202  0.0175  0.0175  0.0146  0.0146  0.0145
  0.0145  0.0082  0.0084  0.0084  0.0049  0.0049  0.0021  0.0018  0.0018  0.0013
  0.0013  0.0011  0.0009  0.0009  0.0005  0.0005  0.0005  0.0003  0.0003  0.0002
  0.0002  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8549.49295961
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       293.71180607
  PAW double counting   =      7797.46702133    -7776.09270098
  entropy T*S    EENTRO =        -0.01649207
  eigenvalues    EBANDS =   -985826.73870430
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -984895.52630880 eV

  energy without entropy =  -984895.50981672  energy(sigma->0) =  -984895.52081144


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




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



 eigenvalue-minimisations  :   360
 total energy-change (2. order) : 0.7769425E+06  (-0.2014474E+06)
 number of electron      88.1202984 magnetization       1.1640182
 augmentation part        3.3610738 magnetization       2.9166597

 Broyden mixing:
  rms(total) = 0.10488E+03    rms(broyden)= 0.10488E+03
  rms(prec ) = 0.10488E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0263
  0.1585  0.1585  0.1589  0.0571  0.0571  0.0486  0.0486  0.0408  0.0408  0.0335
  0.0335  0.0247  0.0247  0.0250  0.0250  0.0174  0.0137  0.0137  0.0122  0.0098
  0.0098  0.0057  0.0057  0.0054  0.0054  0.0058  0.0046  0.0044  0.0009  0.0009
  0.0008  0.0005  0.0005  0.0002  0.0001  0.0001  0.0001  0.0001  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8551.99242319
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       293.49583272
  PAW double counting   =      7716.17695492    -7694.69297146
  entropy T*S    EENTRO =        -0.03286847
  eigenvalues    EBANDS =   -208881.58086674
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -207952.99062146 eV

  energy without entropy =  -207952.95775299  energy(sigma->0) =  -207952.97966530


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




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



 eigenvalue-minimisations  :   368
 total energy-change (2. order) : 0.6196408E+05  (-0.1469216E+06)
 number of electron      88.1157544 magnetization       1.2978650
 augmentation part        3.1642803 magnetization       7.5114361

 Broyden mixing:
  rms(total) = 0.11525E+03    rms(broyden)= 0.11525E+03
  rms(prec ) = 0.11526E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0354
  0.3464  0.2474  0.1390  0.1390  0.0632  0.0537  0.0537  0.0562  0.0337  0.0337
  0.0332  0.0332  0.0275  0.0275  0.0251  0.0251  0.0171  0.0124  0.0124  0.0098
  0.0098  0.0075  0.0075  0.0057  0.0057  0.0049  0.0049  0.0049  0.0049  0.0010
  0.0010  0.0008  0.0004  0.0004  0.0002  0.0001  0.0001  0.0001  0.0001  0.0001
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8476.75200554
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       303.47927064
  PAW double counting   =      7586.57913298    -7564.87949571
  entropy T*S    EENTRO =         0.05574612
  eigenvalues    EBANDS =   -147003.02790839
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -145988.90953915 eV

  energy without entropy =  -145988.96528526  energy(sigma->0) =  -145988.92812118


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) : 0.1420255E+06  (-0.4288032E+04)
 number of electron      87.5373680 magnetization       1.3298859
 augmentation part        3.3703103 magnetization       5.2880008

 Broyden mixing:
  rms(total) = 0.12873E+03    rms(broyden)= 0.12873E+03
  rms(prec ) = 0.12874E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0393
  0.3071  0.3071  0.1610  0.1458  0.1458  0.0788  0.0626  0.0626  0.0357  0.0357
  0.0310  0.0310  0.0303  0.0272  0.0272  0.0233  0.0233  0.0168  0.0150  0.0150
  0.0114  0.0099  0.0099  0.0056  0.0056  0.0051  0.0051  0.0044  0.0044  0.0036
  0.0012  0.0012  0.0007  0.0005  0.0005  0.0002  0.0001  0.0001  0.0001  0.0001
  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8481.11893495
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       299.27108569
  PAW double counting   =      7738.56855332    -7717.19775032
  entropy T*S    EENTRO =         0.03465782
  eigenvalues    EBANDS =     -4968.57361061
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -3963.38027828 eV

  energy without entropy =    -3963.41493610  energy(sigma->0) =    -3963.39183088


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




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



 eigenvalue-minimisations  :   520
 total energy-change (2. order) :-0.4756438E+05  (-0.6525061E+04)
 number of electron      87.3539767 magnetization       1.3367346
 augmentation part        1.9676211 magnetization      35.6162465

 Broyden mixing:
  rms(total) = 0.15313E+03    rms(broyden)= 0.15313E+03
  rms(prec ) = 0.15317E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0408
  0.3236  0.3236  0.1618  0.1552  0.1552  0.0813  0.0653  0.0653  0.0361  0.0361
  0.0354  0.0354  0.0313  0.0313  0.0272  0.0272  0.0259  0.0259  0.0162  0.0125
  0.0125  0.0095  0.0094  0.0094  0.0070  0.0070  0.0057  0.0057  0.0041  0.0041
  0.0017  0.0017  0.0011  0.0010  0.0003  0.0003  0.0002  0.0002  0.0002  0.0001
  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8479.97161251
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       298.42399606
  PAW double counting   =      7734.04362781    -7712.67714975
  entropy T*S    EENTRO =        -0.00685267
  eigenvalues    EBANDS =    -52533.20645787
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -51527.75872816 eV

  energy without entropy =   -51527.75187549  energy(sigma->0) =   -51527.75644394


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




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



 eigenvalue-minimisations  :   504
 total energy-change (2. order) :-0.4906396E+05  (-0.1001356E+06)
 number of electron      87.2442409 magnetization       1.2643468
 augmentation part        1.1958259 magnetization      43.2353966

 Broyden mixing:
  rms(total) = 0.18502E+03    rms(broyden)= 0.18502E+03
  rms(prec ) = 0.18505E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0401
  0.3112  0.3112  0.1600  0.1527  0.1527  0.0868  0.0631  0.0631  0.0569  0.0367
  0.0369  0.0369  0.0301  0.0301  0.0313  0.0313  0.0256  0.0256  0.0166  0.0121
  0.0119  0.0119  0.0086  0.0086  0.0088  0.0088  0.0087  0.0057  0.0057  0.0046
  0.0035  0.0035  0.0011  0.0011  0.0008  0.0004  0.0004  0.0002  0.0001  0.0001
  0.0000  0.0000  0.0001  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8483.53356323
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       298.47500938
  PAW double counting   =      7749.77761179    -7728.42490648
  entropy T*S    EENTRO =        -0.04154144
  eigenvalues    EBANDS =   -101593.60398988
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -100591.71565908 eV

  energy without entropy =  -100591.67411764  energy(sigma->0) =  -100591.70181194


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) : 0.8375656E+05  (-0.1627959E+05)
 number of electron      88.6901378 magnetization       1.1242683
 augmentation part        1.8751213 magnetization      28.0421565

 Broyden mixing:
  rms(total) = 0.14217E+03    rms(broyden)= 0.14217E+03
  rms(prec ) = 0.14220E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0381
  0.2575  0.2575  0.2198  0.0977  0.0977  0.0795  0.0619  0.0619  0.0436  0.0436
  0.0313  0.0313  0.0290  0.0290  0.0216  0.0234  0.0173  0.0173  0.0157  0.0157
  0.0096  0.0092  0.0092  0.0081  0.0081  0.0061  0.0061  0.0039  0.0031  0.0031
  0.0024  0.0023  0.0004  0.0003  0.0003  0.0002  0.0002  0.0001  0.0000  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8481.65007880
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       295.61304331
  PAW double counting   =      7732.72380371    -7711.35461944
  entropy T*S    EENTRO =         0.00030007
  eigenvalues    EBANDS =    -17836.12156394
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -16835.15339433 eV

  energy without entropy =   -16835.15369439  energy(sigma->0) =   -16835.15349435


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




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



 eigenvalue-minimisations  :   408
 total energy-change (2. order) :-0.1062869E+05  (-0.2822708E+05)
 number of electron      89.1142455 magnetization       1.0823738
 augmentation part        2.6044074 magnetization      14.2306998

 Broyden mixing:
  rms(total) = 0.10049E+03    rms(broyden)= 0.10049E+03
  rms(prec ) = 0.10051E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0445
  0.2736  0.2736  0.2239  0.2239  0.1152  0.1152  0.0722  0.0722  0.0770  0.0361
  0.0361  0.0331  0.0331  0.0241  0.0241  0.0218  0.0218  0.0197  0.0184  0.0167
  0.0167  0.0114  0.0076  0.0077  0.0077  0.0082  0.0082  0.0055  0.0055  0.0027
  0.0026  0.0026  0.0026  0.0005  0.0004  0.0004  0.0002  0.0001  0.0001  0.0000
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8472.46114259
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       292.75378271
  PAW double counting   =      7703.34328934    -7681.97942411
  entropy T*S    EENTRO =         0.01317425
  eigenvalues    EBANDS =    -28471.14805972
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -27463.84265936 eV

  energy without entropy =   -27463.85583361  energy(sigma->0) =   -27463.84705078


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




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



 eigenvalue-minimisations  :   456
 total energy-change (2. order) : 0.6930171E+04  (-0.1336240E+05)
 number of electron      88.7849758 magnetization       1.2350143
 augmentation part        2.2199062 magnetization      23.1279258

 Broyden mixing:
  rms(total) = 0.10763E+03    rms(broyden)= 0.10763E+03
  rms(prec ) = 0.10765E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0532
  0.3364  0.3364  0.2699  0.2346  0.2346  0.1096  0.1096  0.0742  0.0742  0.0761
  0.0370  0.0370  0.0294  0.0294  0.0273  0.0273  0.0232  0.0232  0.0187  0.0210
  0.0155  0.0155  0.0124  0.0088  0.0088  0.0084  0.0084  0.0055  0.0055  0.0046
  0.0028  0.0027  0.0027  0.0024  0.0005  0.0004  0.0004  0.0002  0.0001  0.0001
  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8485.48343307
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       289.06792716
  PAW double counting   =      8083.03860900    -8062.18286331
  entropy T*S    EENTRO =        -0.01234584
  eigenvalues    EBANDS =    -21523.73549439
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -20533.67187968 eV

  energy without entropy =   -20533.65953385  energy(sigma->0) =   -20533.66776440


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) :-0.1002652E+06  (-0.1422741E+05)
 number of electron      87.9753178 magnetization       1.1495271
 augmentation part        2.1914688 magnetization      24.4230155

 Broyden mixing:
  rms(total) = 0.13821E+03    rms(broyden)= 0.13821E+03
  rms(prec ) = 0.13823E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0548
  0.3563  0.3563  0.2899  0.2265  0.2265  0.0954  0.0954  0.1100  0.1100  0.0730
  0.0467  0.0336  0.0336  0.0254  0.0254  0.0279  0.0279  0.0193  0.0193  0.0209
  0.0209  0.0169  0.0165  0.0167  0.0097  0.0097  0.0083  0.0083  0.0061  0.0061
  0.0050  0.0037  0.0028  0.0025  0.0025  0.0005  0.0002  0.0002  0.0002  0.0002
  0.0002  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8505.96869834
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       287.38403844
  PAW double counting   =      8192.97509557    -8172.53320599
  entropy T*S    EENTRO =         0.05368561
  eigenvalues    EBANDS =   -121766.43477915
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -120798.88814311 eV

  energy without entropy =  -120798.94182872  energy(sigma->0) =  -120798.90603831


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




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



 eigenvalue-minimisations  :   520
 total energy-change (2. order) :-0.1014515E+07  (-0.5126666E+06)
 number of electron      86.9423068 magnetization       1.2009118
 augmentation part        1.4662322 magnetization      41.4387605

 Broyden mixing:
  rms(total) = 0.28090E+03    rms(broyden)= 0.28090E+03
  rms(prec ) = 0.28093E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0551
  0.3592  0.3592  0.3046  0.2160  0.2160  0.1064  0.1064  0.0765  0.0765  0.0776
  0.0656  0.0656  0.0412  0.0412  0.0320  0.0320  0.0278  0.0278  0.0214  0.0208
  0.0208  0.0195  0.0195  0.0171  0.0123  0.0123  0.0106  0.0083  0.0083  0.0049
  0.0049  0.0023  0.0023  0.0027  0.0027  0.0024  0.0005  0.0004  0.0004  0.0002
  0.0002  0.0001  0.0000  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8503.52571458
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       286.40249878
  PAW double counting   =      8086.21081933    -8065.56962051
  entropy T*S    EENTRO =        -0.00552618
  eigenvalues    EBANDS =  -1136282.54430556
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1135313.39612795 eV

  energy without entropy = -1135313.39060177  energy(sigma->0) = -1135313.39428589


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




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



 eigenvalue-minimisations  :   504
 total energy-change (2. order) : 0.3684212E+06  (-0.5603447E+06)
 number of electron      87.7057260 magnetization       1.5771730
 augmentation part        0.8575080 magnetization      52.0433058

 Broyden mixing:
  rms(total) = 0.44966E+03    rms(broyden)= 0.44966E+03
  rms(prec ) = 0.44968E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0539
  0.3165  0.3165  0.2838  0.1940  0.1258  0.1258  0.1093  0.1093  0.0960  0.0606
  0.0606  0.0339  0.0339  0.0403  0.0403  0.0279  0.0279  0.0229  0.0229  0.0210
  0.0137  0.0137  0.0125  0.0108  0.0108  0.0038  0.0038  0.0036  0.0036  0.0020
  0.0021  0.0021  0.0016  0.0007  0.0007  0.0007  0.0005  0.0000  0.0001  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8508.07707024
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       285.94023385
  PAW double counting   =      8104.11347358    -8083.47928705
  entropy T*S    EENTRO =         0.02909998
  eigenvalues    EBANDS =   -767856.40000560
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -766892.23783472 eV

  energy without entropy =  -766892.26693470  energy(sigma->0) =  -766892.24753471


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) :-0.1230653E+07  (-0.4821958E+06)
 number of electron      87.1866122 magnetization       1.5747920
 augmentation part        0.5922439 magnetization      52.9658328

 Broyden mixing:
  rms(total) = 0.49766E+03    rms(broyden)= 0.49766E+03
  rms(prec ) = 0.49767E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0550
  0.3228  0.3228  0.2835  0.1991  0.1450  0.1450  0.1274  0.1274  0.0973  0.0449
  0.0449  0.0439  0.0439  0.0344  0.0344  0.0281  0.0281  0.0260  0.0260  0.0164
  0.0164  0.0147  0.0147  0.0137  0.0137  0.0097  0.0097  0.0039  0.0025  0.0025
  0.0022  0.0022  0.0019  0.0015  0.0007  0.0006  0.0006  0.0006  0.0000  0.0001
  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8500.68394206
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       288.95366208
  PAW double counting   =      7758.88425505    -7737.56925771
  entropy T*S    EENTRO =        -0.05958347
  eigenvalues    EBANDS =  -1998520.36694607
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1997545.20609142 eV

  energy without entropy = -1997545.14650795  energy(sigma->0) = -1997545.18623026


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




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



 eigenvalue-minimisations  :   488
 total energy-change (2. order) :-0.3851615E+07  (-0.3973139E+07)
 number of electron      87.1461975 magnetization       1.6219334
 augmentation part        1.5807750 magnetization      21.0557841

 Broyden mixing:
  rms(total) = 0.50139E+03    rms(broyden)= 0.50139E+03
  rms(prec ) = 0.50140E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0564
  0.2943  0.2943  0.2874  0.2276  0.2276  0.1434  0.1434  0.1192  0.1192  0.0512
  0.0479  0.0479  0.0342  0.0342  0.0335  0.0335  0.0341  0.0341  0.0326  0.0204
  0.0141  0.0141  0.0132  0.0118  0.0109  0.0109  0.0078  0.0078  0.0038  0.0029
  0.0029  0.0020  0.0020  0.0020  0.0007  0.0011  0.0006  0.0006  0.0007  0.0001
  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8500.68439989
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       288.83730407
  PAW double counting   =      7785.90075286    -7764.63271532
  entropy T*S    EENTRO =        -0.05947497
  eigenvalues    EBANDS =  -5850134.96712652
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -5849159.96993900 eV

  energy without entropy = -5849159.91046403  energy(sigma->0) = -5849159.95011401


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




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



 eigenvalue-minimisations  :   464
 total energy-change (2. order) : 0.4163453E+07  (-0.6409688E+06)
 number of electron      87.1107119 magnetization       1.6175606
 augmentation part        0.7057708 magnetization      36.7839261

 Broyden mixing:
  rms(total) = 0.68581E+03    rms(broyden)= 0.68581E+03
  rms(prec ) = 0.68582E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0581
  0.3132  0.3132  0.2829  0.2448  0.2448  0.1469  0.1469  0.1087  0.1087  0.0897
  0.0501  0.0501  0.0309  0.0309  0.0340  0.0340  0.0328  0.0328  0.0337  0.0337
  0.0206  0.0178  0.0178  0.0127  0.0127  0.0117  0.0117  0.0047  0.0047  0.0041
  0.0026  0.0026  0.0023  0.0017  0.0018  0.0018  0.0007  0.0005  0.0005  0.0004
  0.0000  0.0001  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8489.60987574
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       290.47070231
  PAW double counting   =      7782.83811690    -7761.65680047
  entropy T*S    EENTRO =        -0.01454365
  eigenvalues    EBANDS =  -1686694.25056733
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1685706.58724723 eV

  energy without entropy = -1685706.57270357  energy(sigma->0) = -1685706.58239934


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




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



 eigenvalue-minimisations  :   472
 total energy-change (2. order) :-0.3330779E+07  (-0.1471510E+07)
 number of electron      88.2169963 magnetization       2.0121448
 augmentation part        1.2117075 magnetization      20.1745056

 Broyden mixing:
  rms(total) = 0.43096E+03    rms(broyden)= 0.43096E+03
  rms(prec ) = 0.43097E+03
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.0587
  0.3220  0.3220  0.2776  0.2465  0.2465  0.1267  0.1267  0.1190  0.1190  0.0905
  0.0802  0.0802  0.0338  0.0338  0.0358  0.0358  0.0309  0.0309  0.0276  0.0276
  0.0244  0.0244  0.0179  0.0179  0.0142  0.0194  0.0183  0.0089  0.0036  0.0037
  0.0037  0.0024  0.0024  0.0021  0.0021  0.0017  0.0005  0.0006  0.0006  0.0005
  0.0001  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8488.57534784
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       290.11390295
  PAW double counting   =      7784.27389440    -7763.09212245
  entropy T*S    EENTRO =        -0.01432439
  eigenvalues    EBANDS =  -5017474.20505674
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -5016485.86333330 eV

  energy without entropy = -5016485.84900891  energy(sigma->0) = -5016485.85855850


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




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



 eigenvalue-minimisations  :   488
 total energy-change (2. order) :-0.7197263E+07  (-0.1042155E+08)
 number of electron      88.2169963 magnetization       2.0121448
 augmentation part        1.2117075 magnetization      20.1745056

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      5655.20747896
  -Hartree energ DENC   =     -8441.07069835
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       296.31021889
  PAW double counting   =      7543.58940070    -7521.86330634
  entropy T*S    EENTRO =        -0.07488351
  eigenvalues    EBANDS = -12214791.34746988
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  = -12213748.82101773 eV

  energy without entropy =-12213748.74613422  energy(sigma->0) =-12213748.79605656


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




 average (electrostatic) potential at core
  the test charge radii are     0.7089  0.5201  0.6991  0.9892  0.7215
  (the norm of the test charge is              1.0000)
       1 -75.4224       2 -42.5141       3 -42.4886       4 -44.0514       5 -43.3525
       6 -43.3260       7 -40.5798       8 -41.5573       9 -41.7792      10 -42.4807
      11 -43.1381      12 -43.1879      13 -41.7724      14 -42.4555      15 -42.8095
      16 -43.7020      17 -43.4973      18 -43.8019      19 -41.3759      20 -41.3686
      21 -42.7851      22 -42.7247      23 -44.4140      24 -61.5574      25 -61.8531
      26 -55.7614      27 -62.6195      28 -60.6475      29 -62.1209      30 -59.4838
      31 -94.3859      32 -93.9590      33 -85.0601      34 -81.8892      35 -83.9168
      36 -83.1743
 
 
 
 E-fermi : -13.1943     XC(G=0):  -0.7356     alpha+bet : -0.4834


 spin component 1

 k-point     1 :       0.0000    0.0000    0.0000
  band No.  band energies     occupation 
      1     -28.4798      1.00000
      2     -28.3165      1.00000
      3     -27.3720      1.00000
      4     -25.6205      1.00000
      5     -22.0543      1.00000
      6     -20.5024      1.00000
      7     -20.0929      1.00000
      8     -19.3280      1.00000
      9     -19.0736      1.00000
     10     -18.7598      1.00000
     11     -18.0620      1.00000
     12     -16.5903      1.00000
     13     -15.1455      1.00000
     14     -14.9698      1.00000
     15     -13.6896      1.00129
     16     -13.4782      1.03103
     17     -13.2794      0.82633
     18     -13.1698      0.39701
     19     -13.1644      0.37515
     20     -13.0266      0.00068
     21     -12.7546     -0.00400
     22     -12.5351     -0.00002
     23     -12.4772     -0.00000
     24     -12.4125     -0.00000
     25     -12.2601     -0.00000
     26     -12.1098     -0.00000
     27     -11.8740     -0.00000
     28     -11.4265     -0.00000
     29     -11.2258     -0.00000
     30     -11.1245     -0.00000
     31     -11.1178     -0.00000
     32     -10.6882     -0.00000
     33     -10.5230     -0.00000
     34     -10.2157     -0.00000
     35     -10.1199     -0.00000
     36      -9.9187     -0.00000
     37      -9.6279     -0.00000
     38      -9.5885     -0.00000
     39      -9.4717     -0.00000
     40      -9.1730     -0.00000
     41      -8.7839     -0.00000
     42      -8.0918     -0.00000
     43      -7.9243     -0.00000
     44      -7.4980      0.00000
     45      -0.0188      0.00000
     46       0.1889      0.00000
     47       0.2373      0.00000
     48       0.3169      0.00000
     49       0.7674      0.00000
     50       0.8086      0.00000
     51       0.8651      0.00000
     52       0.9053      0.00000
     53       0.9556      0.00000
     54       0.9864      0.00000
     55       1.0640      0.00000
     56       1.2219      0.00000
     57       1.2523      0.00000
     58       1.4055      0.00000
     59       1.4798      0.00000
     60       1.5397      0.00000
     61       1.5473      0.00000
     62       1.6134      0.00000
     63       1.6575      0.00000
     64       1.6582      0.00000
     65       1.7123      0.00000
     66       1.7378      0.00000
     67       1.8022      0.00000
     68       1.8430      0.00000
     69       1.8488      0.00000
     70       1.9157      0.00000
     71       1.9719      0.00000
     72       2.0677      0.00000
     73       2.1480      0.00000
     74       2.1541      0.00000
     75       2.3132      0.00000
     76       2.3527      0.00000
     77       2.4459      0.00000
     78       2.5555      0.00000
     79       2.6095      0.00000
     80       2.7577      0.00000

 spin component 2

 k-point     1 :       0.0000    0.0000    0.0000
  band No.  band energies     occupation 
      1   **********      1.00000
      2   **********      1.00000
      3   **********      1.00000
      4   **********      1.00000
      5   **********      1.00000
      6   **********      1.00000
      7   **********      1.00000
      8   **********      1.00000
      9   **********      1.00000
     10   **********      1.00000
     11   **********      1.00000
     12   **********      1.00000
     13   **********      1.00000
     14   **********      1.00000
     15   **********      1.00000
     16   **********      1.00000
     17   **********      1.00000
     18   **********      1.00000
     19   **********      1.00000
     20   **********      1.00000
     21   **********      1.00000
     22   **********      1.00000
     23   **********      1.00000
     24   **********      1.00000
     25   **********      1.00000
     26   **********      1.00000
     27   -5900.0794      1.00000
     28   -5796.0001      1.00000
     29   -5563.5786      1.00000
     30   -4819.3258      1.00000
     31   -4438.7816      1.00000
     32   -3550.5418      1.00000
     33   -3333.0364      1.00000
     34   -2921.1506      1.00000
     35   -2819.5475      1.00000
     36   -2358.3386      1.00000
     37   -1827.6189      1.00000
     38   -1784.2302      1.00000
     39   -1741.2648      1.00000
     40   -1717.9806      1.00000
     41   -1477.5876      1.00000
     42   -1122.7814      1.00000
     43    -932.5320      1.00000
     44    -900.1694      1.00000
     45    -752.7113      1.00000
     46    -729.2241      1.00000
     47    -660.7988      1.00000
     48    -635.6463      1.00000
     49    -305.7613      1.00000
     50    -284.1936      1.00000
     51    -234.4116      1.00000
     52    -212.9054      1.00000
     53    -179.5007      1.00000
     54    -162.6481      1.00000
     55    -158.0091      1.00000
     56    -130.9625      1.00000
     57    -127.6812      1.00000
     58     -50.6748      1.00000
     59     -44.7383      1.00000
     60     -37.4152      1.00000
     61     -28.0358      1.00000
     62     -26.4470      1.00000
     63     -23.2625      1.00000
     64     -21.3406      1.00000
     65     -19.9044      1.00000
     66     -18.1455      1.00000
     67     -17.1799      1.00000
     68     -15.9975      1.00000
     69     -14.2036      1.00000
     70     -13.1638      0.37252
     71     -11.2767     -0.00000
     72      -8.3371     -0.00000
     73      -5.5207      0.00000
     74      -4.4258      0.00000
     75      -1.7395      0.00000
     76      -0.9174      0.00000
     77       3.5655      0.00000
     78      14.3781      0.00000
     79     446.9351      0.00000
     80    1634.4199      0.00000


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


 soft charge-density along one line, spin component           1
         0         1         2         3         4         5         6         7         8         9
 total charge-density along one line
 
 soft charge-density along one line, spin component           2
         0         1         2         3         4         5         6         7         8         9
 total charge-density along one line
 
 pseudopotential strength for first ion, spin component:           1
 14.414  17.318  -0.005  -0.004  -0.004  -0.009  -0.007  -0.009
 17.318  20.810  -0.006  -0.004  -0.005  -0.011  -0.008  -0.012
 -0.005  -0.006  -7.667   0.033   0.023 -10.698   0.052   0.037
 -0.004  -0.004   0.033  -7.714  -0.004   0.052 -10.775  -0.006
 -0.004  -0.005   0.023  -0.004  -7.618   0.037  -0.006 -10.620
 -0.009  -0.011 -10.698   0.052   0.037 -14.288   0.083   0.059
 -0.007  -0.008   0.052 -10.775  -0.006   0.083 -14.412  -0.008
 -0.009  -0.012   0.037  -0.006 -10.620   0.059  -0.008 -14.165
 pseudopotential strength for first ion, spin component:           2
 13.871  16.654  -0.007   0.002   0.003  -0.015   0.002  -0.005
 16.654  19.997  -0.009   0.002   0.003  -0.018   0.002  -0.005
 -0.007  -0.009  -7.348  -0.044  -0.032 -10.204  -0.071  -0.052
  0.002   0.002  -0.044  -7.260  -0.008  -0.071 -10.064  -0.012
  0.003   0.003  -0.032  -0.008  -7.417  -0.052  -0.012 -10.317
 -0.015  -0.018 -10.204  -0.071  -0.052 -13.529  -0.114  -0.082
  0.002   0.002  -0.071 -10.064  -0.012  -0.114 -13.306  -0.017
 -0.005  -0.005  -0.052  -0.012 -10.317  -0.082  -0.017 -13.709
 total augmentation occupancy for first ion, spin component:           1
******* *******  -2.090  -4.036   1.321   4.709   4.002  -1.096
******* 723.206   1.880   3.421  -0.604  -4.008  -3.388   0.740
 -2.090   1.880   6.862  -6.542   5.648  -4.674   3.973  -3.576
 -4.036   3.421  -6.542  10.153  -6.268   4.586  -5.690   3.775
  1.321  -0.604   5.648  -6.268   6.385  -3.985   3.670  -3.917
  4.709  -4.008  -4.674   4.586  -3.985   3.425  -2.698   2.507
  4.002  -3.388   3.973  -5.690   3.670  -2.698   3.528  -2.246
 -1.096   0.740  -3.576   3.775  -3.917   2.507  -2.246   2.539
 total augmentation occupancy for first ion, spin component:           2
******* 863.345   2.299   4.239  -0.405  -4.765  -4.063   0.808
863.345 *******  -1.956  -3.504   0.281   4.028   3.411  -0.638
  2.299  -1.956  -6.841   6.717  -5.648   4.669  -4.001   3.570
  4.239  -3.504   6.717  -7.943   5.809  -4.625   5.342  -3.696
 -0.405   0.281  -5.648   5.809  -6.113   3.983  -3.603   3.845
 -4.765   4.028   4.669  -4.625   3.983  -3.424   2.705  -2.505
 -4.063   3.411  -4.001   5.342  -3.603   2.705  -3.473   2.236
  0.808  -0.638   3.570  -3.696   3.845  -2.505   2.236  -2.519


------------------------ aborting loop EDIFF was not reached (unconverged)  ----------------------------


 -----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|     The electronic self-consistency was not achieved in the given           |
|     number of steps (NELM). The forces and other quantities evaluated       |
|     might not be reliable so examine the results carefully. If you find     |
|     spurious results, we suggest increasing NELM, if you were close to      |
|     convergence or switching to a different ALGO or adjusting the           |
|     density mixing parameters otherwise.                                    |
|                                                                             |
 -----------------------------------------------------------------------------


  FORCE on cell =-STRESS in cart. coord.  units (eV):
  Direction    XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------------
  Alpha Z     6.34848     6.34848     6.34848
  Ewald    2965.19875   960.59903  1729.40690 -1102.41236   188.83938   234.00523
  Hartree  1662.89370   941.34841  1187.27822  -412.03015    32.12259    68.03771
  E(xc)    -340.66486  -342.13411  -341.29019    -0.72120     0.29202     0.29237
  Local   -5003.72492 -2577.86588 -3443.41248  1377.54137  -181.45492  -276.56498
  n-local 69567.06184 65671.50104111135.23507-15581.28656  9095.31128 41625.85575
  augment     2.60976     0.53371     0.28944    -5.43478     4.79903     4.35079
  Kinetic  2750.92000  3069.30108  2813.45108    93.23951   -63.35128   -92.74729
  Fock        0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  -------------------------------------------------------------------------------------
  Total    71610.6427449  67729.6317571 113087.3065175 -15631.1041812   9076.5581115  41563.2295988
  in kB    25496.2107539  24114.4179045  40263.5375140  -5565.2890582   3231.6123646  14798.1476055
  external PRESSURE =   29958.0553908 kB  Pullay stress =       0.0000000 kB


 VOLUME and BASIS-vectors are now :
 -----------------------------------------------------------------------------
  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


 FORCES acting on ions
    electron-ion (+dipol)            ewald-force                    non-local-force                 convergence-correction
 -----------------------------------------------------------------------------------------------
   0.325E+02 0.201E+02 0.307E+00   0.157E-01 -.103E+02 -.223E+01   -.604E+05 0.163E+06 0.871E+05   -.290E+02 -.710E+01 -.162E+01
   -.179E+02 -.528E+01 -.373E+02   0.374E+02 0.121E+02 0.670E+02   0.335E+04 0.513E+04 -.114E+05   -.110E+02 -.339E+01 -.138E+02
   -.329E+02 -.746E-01 0.260E+02   0.636E+02 0.919E+00 -.457E+02   0.102E+03 -.712E+03 -.344E+03   -.157E+02 -.940E+00 0.882E+01
   0.594E+02 0.799E+01 0.197E+02   -.901E+02 -.148E+02 -.261E+02   -.780E+03 0.114E+04 0.280E+04   0.176E+02 0.332E+01 0.759E+01
   0.160E+02 -.528E+02 0.214E+02   -.224E+02 0.811E+02 -.300E+02   0.235E+04 0.849E+04 -.406E+04   0.738E+01 -.140E+02 0.948E+01
   0.204E+01 0.375E+01 0.611E+02   -.197E+01 -.721E+01 -.921E+02   0.525E+04 0.477E+03 -.535E+04   0.387E+01 0.101E+01 0.178E+02
   0.410E+02 0.648E+00 -.408E+02   -.644E+02 -.644E+01 0.684E+02   0.313E+04 0.172E+04 -.237E+04   0.159E+02 0.197E+01 -.153E+02
   0.181E+02 -.486E+02 -.132E+02   -.241E+02 0.825E+02 0.183E+02   -.122E+04 -.328E+04 -.444E+03   0.868E+01 -.189E+02 -.451E+01
   0.509E+02 -.214E+01 0.199E+02   -.792E+02 -.191E+01 -.401E+02   -.849E+03 0.139E+04 -.654E+04   0.195E+02 0.880E+00 0.835E+01
   -.188E+02 -.249E+01 0.559E+02   0.296E+02 0.724E+01 -.927E+02   0.201E+04 0.147E+04 -.563E+03   -.516E+01 -.334E+00 0.219E+02
   -.160E+02 0.523E+02 0.197E+02   0.230E+02 -.857E+02 -.267E+02   0.859E+03 0.124E+04 -.332E+04   -.297E+01 0.194E+02 0.890E+01
   0.287E+02 0.180E+02 0.299E+02   -.565E+02 -.266E+02 -.438E+02   0.480E+04 0.952E+03 -.258E+04   0.124E+02 0.787E+01 0.127E+02
   -.245E+02 0.452E+02 -.145E+02   0.379E+02 -.814E+02 0.253E+02   0.276E+04 0.366E+04 -.396E+04   -.978E+01 0.207E+02 -.534E+01
   -.339E+02 0.131E+02 0.372E+02   0.516E+02 -.197E+02 -.660E+02   0.443E+04 0.170E+04 -.646E+03   -.138E+02 0.691E+01 0.150E+02
   -.419E+02 -.105E+02 -.239E+02   0.708E+02 0.210E+02 0.383E+02   -.418E+04 -.793E+03 0.465E+04   -.186E+02 -.191E+01 -.807E+01
   0.190E+02 -.457E+02 0.715E+01   -.445E+02 0.772E+02 -.160E+02   -.881E+04 -.585E+04 0.180E+05   0.119E+02 -.200E+02 0.461E+01
   -.133E+02 -.333E+02 -.358E+02   0.184E+02 0.549E+02 0.662E+02   0.169E+04 0.729E+03 0.158E+03   0.430E-02 -.162E+02 -.151E+02
   -.271E+02 -.324E+02 0.187E+02   0.470E+02 0.517E+02 -.417E+02   -.278E+04 -.810E+03 0.195E+04   -.659E+01 -.156E+02 0.121E+02
   -.943E+01 -.420E+02 -.115E+01   0.785E+01 0.818E+02 -.873E+01   0.162E+04 0.359E+04 -.310E+04   -.311E+01 -.238E+02 0.309E+01
   -.964E+01 -.161E+02 -.418E+02   0.105E+02 0.222E+02 0.848E+02   0.163E+04 0.873E+03 -.682E+04   -.289E+01 -.620E+01 -.252E+02
   -.426E+02 -.178E+02 -.642E+01   0.839E+02 0.250E+02 0.471E+01   -.148E+04 -.154E+04 0.442E+03   -.244E+02 -.786E+01 -.779E+00
   0.105E+01 0.905E+01 -.527E+02   -.748E+01 -.155E+02 0.924E+02   0.138E+02 0.547E+03 -.166E+04   0.500E+01 0.641E+01 -.231E+02
   0.155E+02 0.512E+02 -.584E+01   -.314E+02 -.873E+02 0.211E+01   -.189E+04 -.232E+04 0.420E+04   0.105E+02 0.223E+02 -.325E+00
   0.269E+01 0.360E+02 -.800E+02   -.277E+02 -.750E+02 0.100E+03   -.290E+05 -.362E+05 0.605E+05   0.214E+02 0.354E+02 -.227E+02
   0.114E+03 -.592E+02 0.157E+03   -.156E+03 0.695E+02 -.202E+03   0.187E+05 -.830E+04 -.442E+05   0.404E+02 -.851E+01 0.466E+02
   0.169E+03 -.796E+02 -.542E+02   -.227E+03 0.915E+02 0.722E+02   0.149E+05 -.314E+05 -.356E+04   0.616E+02 -.183E+02 -.184E+02
   -.978E+01 0.891E+02 0.166E+03   0.666E+00 -.125E+03 -.220E+03   -.471E+02 -.629E+05 -.624E+05   0.639E+01 0.337E+02 0.558E+02
   -.159E+03 0.851E+02 -.131E+02   0.211E+03 -.120E+03 0.191E+01   -.699E+05 0.244E+04 0.681E+05   -.515E+02 0.355E+02 0.519E+01
   -.337E+02 -.117E+03 -.216E+02   0.848E+01 0.175E+03 0.176E+02   -.219E+05 0.175E+05 0.389E+05   0.188E+02 -.583E+02 0.429E+01
   -.688E+02 -.899E+02 -.653E+02   0.103E+03 0.124E+03 0.843E+02   -.721E+04 0.205E+04 0.428E+05   -.319E+02 -.342E+02 -.221E+02
   -.169E+02 -.531E+02 -.719E+02   0.215E+02 0.335E+02 0.731E+02   0.366E+05 -.121E+05 -.321E+04   -.398E+01 0.953E+01 -.602E+01
   -.325E+02 -.339E+02 -.536E+02   0.298E+01 0.346E+02 0.477E+02   0.218E+05 -.422E+04 -.168E+05   0.226E+02 -.182E+01 0.269E+01
   0.247E+02 0.122E+03 0.848E+02   -.509E+02 -.166E+03 -.150E+03   0.782E+05 0.130E+06 -.101E+05   0.371E+02 0.323E+02 0.671E+02
   0.642E+02 0.835E+02 -.166E+03   -.117E+03 -.139E+03 0.223E+03   -.778E+05 -.443E+06 0.737E+05   0.477E+02 0.440E+02 -.535E+02
   -.161E+03 -.134E+03 0.118E+03   0.184E+03 0.158E+03 -.171E+03   -.355E+04 0.907E+04 0.671E+05   -.236E+02 -.919E+01 0.549E+02
   0.186E+01 0.163E+03 -.145E+03   -.127E+02 -.222E+03 0.188E+03   0.175E+05 0.409E+05 0.162E+06   0.576E+01 0.567E+02 -.283E+02
 -----------------------------------------------------------------------------------------------
   -.110E+03 -.752E+02 -.101E+03   0.320E-13 0.284E-13 0.284E-13   -.701E+05 -.215E+06 0.439E+06   0.121E+03 0.712E+02 0.103E+03
 
 
 POSITION                                       TOTAL-FORCE (eV/Angst)
 -----------------------------------------------------------------------------------
     14.38439      6.60419      8.88015    -58475.332123 169205.010683  74913.949697
     14.83303      6.78684      9.81459      5302.242768  11094.728256 -23596.015123
     15.20133      6.62793      8.23938      2063.454896   5252.737649 -12543.374143
      2.98896      7.76510      5.49719      1153.693885   7102.009210  -9389.491994
      4.09729      9.17344      5.53190      4296.194578  14465.996111 -16247.987405
      4.39327      7.83128      4.39425      7203.683070   6439.654622 -17552.498065
      2.17582      8.25381      9.85948      5071.037866   7684.991367 -14550.805804
      2.97541      9.58289      8.97763       727.051523   2700.143063 -12632.315706
      2.09812      8.30448      8.07775      1088.915225   7353.561778 -18736.661463
      8.36653      6.47690      4.41819      3960.217528   7441.299107 -12766.512220
      8.34364      5.12407      5.58594      2809.849480   7188.870858 -15509.223332
      6.88889      6.15743      5.37770      6728.694898   6915.978236 -14771.714084
     10.63193      4.14481      8.61069      4711.640998   9608.022228 -16141.183262
     10.64778      5.23839      7.16744      6381.017961   7662.062752 -12848.757464
     11.10606      5.84826      8.79994     -2221.774434   5180.855767  -7530.016333
      4.91450     10.36245      7.56962     -6877.970911    130.923090   5842.080621
      6.05951     10.01073      8.87449      3643.737935   6699.006264 -12015.141931
      6.56641      9.84737      7.17557      -819.014654   5158.265371 -10249.353715
      8.79139      9.30593      8.37936      3565.880834   9567.755419 -15295.567853
      8.97277      8.34882      9.88168      3570.751805   6837.642235 -18989.446711
     10.23075      8.28167      8.62578       485.166467   4421.105285 -11748.950739
      6.46081      6.69422      9.61858      1959.087476   6511.727899 -13830.224921
      6.15846      5.67564      8.21540        52.006189   3629.475301  -7989.283173
      6.49820      6.68224      8.51673    -27087.590997 -30197.454587  48309.087467
      4.04442      8.07649      5.40979     20627.898801  -2333.530093 -56428.986030
      2.74494      8.50246      8.95040     16832.369409 -25444.322844 -15750.800354
      7.99075      6.15247      5.39923      1896.925304 -56970.400963 -74626.411527
     10.41603      5.16488      8.24389    -67957.550334   8409.956519  55904.354040
      5.73150      9.68908      7.87422    -19910.471551  23417.358680  26690.084068
      9.14645      8.35507      8.79536     -5262.587347   8018.702200  30625.620290
      8.29348      6.88988      8.02032     38542.813532  -6188.891944 -15401.840872
      5.21505      7.88230      7.89741     23697.337923   1744.555452 -29015.402725
      4.85843      7.38378      6.35019     80198.571168 136169.844471 -22272.091136
      3.92000      7.70061      8.92325    -75809.781948-436790.822451  61467.651262
      8.48389      7.06683      6.37953     -1600.229770  15047.603917  54880.771349
      9.05494      5.49217      8.56562     19452.062553  46865.579091 149796.459293
 -----------------------------------------------------------------------------------
    total drift:                           -70081.587277-214734.578371 438797.240793


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



  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy   TOTEN  = -12213748.8210177291 eV

  energy  without entropy=-12213748.7461342197  energy(sigma->0) =-12213748.79605656
 


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




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



volume of typ            1:     0.0 %

volume of typ            2:     0.1 %

volume of typ            3:     0.3 %

volume of typ            4:     0.3 %

volume of typ            5:     0.1 %
 
total charge     
 
# of ion       s       p       d       tot
------------------------------------------
    1        0.894   0.632   0.017   1.543
    2        0.103  -0.008   0.000   0.095
    3        0.098  -0.007   0.000   0.091
    4        0.105  -0.009   0.000   0.096
    5        0.096   0.006   0.000   0.102
    6        0.099  -0.010   0.000   0.090
    7        0.076  -0.008   0.000   0.068
    8        0.079  -0.007   0.001   0.073
    9        0.077  -0.008   0.000   0.068
   10        0.078  -0.010   0.000   0.068
   11        0.095  -0.009   0.000   0.086
   12        0.114  -0.008   0.000   0.106
   13        0.067  -0.010   0.000   0.056
   14        0.101  -0.009   0.000   0.092
   15        0.081  -0.010   0.000   0.072
   16        0.080   0.107   0.000   0.187
   17        0.078  -0.010   0.000   0.068
   18        0.073  -0.010   0.000   0.064
   19        0.069  -0.010   0.000   0.059
   20        0.083  -0.009   0.000   0.074
   21        0.076  -0.009   0.000   0.067
   22        0.089  -0.009   0.000   0.080
   23        0.094  -0.009   0.000   0.086
   24        0.606   0.584   0.026   1.216
   25        0.600   0.906   0.046   1.552
   26        0.535   1.078   0.049   1.662
   27        0.592   0.819   0.043   1.453
   28        0.559   0.751   0.044   1.354
   29        0.601   0.530   0.026   1.157
   30        0.505   0.463   0.019   0.987
   31        0.572   0.663   0.289   1.524
   32        0.517   0.591   0.277   1.385
   33        1.258   1.770   0.040   3.068
   34        1.326   1.201   0.044   2.571
   35        1.161   1.540   0.032   2.733
   36        1.223   1.535   0.035   2.792
--------------------------------------------------
tot          12.86   13.00    0.99   26.84
 
magnetization (x)
 
# of ion       s       p       d       tot
------------------------------------------
    1        0.084  -0.596  -0.016  -0.528
    2       -0.060   0.009  -0.000  -0.051
    3       -0.056   0.009  -0.000  -0.048
    4       -0.019   0.011  -0.000  -0.009
    5       -0.022  -0.005  -0.000  -0.027
    6       -0.020   0.011  -0.000  -0.009
    7       -0.041   0.009  -0.000  -0.032
    8       -0.048   0.007  -0.001  -0.041
    9       -0.041   0.009  -0.000  -0.032
   10       -0.037   0.011  -0.000  -0.026
   11       -0.013   0.011  -0.000  -0.002
   12       -0.035   0.010  -0.000  -0.025
   13       -0.037   0.011  -0.000  -0.025
   14       -0.053   0.010  -0.000  -0.043
   15       -0.026   0.011  -0.000  -0.015
   16       -0.038  -0.106  -0.000  -0.144
   17       -0.030   0.011  -0.000  -0.019
   18       -0.037   0.011  -0.000  -0.026
   19       -0.033   0.011  -0.000  -0.022
   20       -0.045   0.010  -0.000  -0.036
   21       -0.039   0.010  -0.000  -0.029
   22       -0.043   0.010  -0.000  -0.033
   23       -0.039   0.010  -0.000  -0.029
   24        0.097  -0.411  -0.024  -0.338
   25        0.083  -0.164  -0.029  -0.110
   26        0.061  -1.019  -0.044  -1.002
   27        0.112  -0.201  -0.028  -0.117
   28        0.117  -0.483  -0.034  -0.400
   29        0.113  -0.415  -0.025  -0.327
   30        0.190  -0.429  -0.018  -0.257
   31       -0.111  -0.319  -0.192  -0.622
   32       -0.107  -0.317  -0.191  -0.615
   33        0.003  -0.218  -0.034  -0.248
   34       -0.096  -0.918  -0.042  -1.056
   35        0.077  -0.301  -0.025  -0.248
   36        0.014  -0.756  -0.030  -0.773
--------------------------------------------------
tot          -0.17   -6.46   -0.74   -7.37
 

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

   base      :      30000. kBytes
   nonlr-proj:       4077. kBytes
   fftplans  :      87706. kBytes
   grid      :     314101. kBytes
   one-center:        221. kBytes
   wavefun   :     105128. kBytes
 
  
  
 General timing and accounting informations for this job:
 ========================================================
  
                  Total CPU time used (sec):      341.914
                            User time (sec):      320.294
                          System time (sec):       21.621
                         Elapsed time (sec):      344.213
  
                   Maximum memory used (kb):     1619716.
                   Average memory used (kb):          N/A
  
                          Minor page faults:       373100
                          Major page faults:            0
                 Voluntary context switches:        19344