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.09.12  18:30:28
 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 = 60
   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.555  0.179  0.914-
   2  0.640  0.502  0.950-
   3  0.262  0.652  0.491-
   4  0.745  0.741  0.829-
   5  0.360  0.143  0.733-
   6  0.578  0.322  0.588-
   7  0.439  0.753  0.267-
   8  0.300  0.106  0.935-
   9  0.047  0.260  0.276-
  10  0.438  0.126  0.495-
  11  0.984  0.806  0.373-
  12  0.633  0.794  0.873-  24 0.51
  13  0.226  0.429  0.415-
  14  0.051  0.388  0.168-
  15  0.087  0.671  0.061-
  16  0.877  0.104  0.233-
  17  0.227  0.483  0.328-
  18  0.544  0.823  0.585-
  19  0.981  0.294  0.656-
  20  0.389  0.982  0.007-
  21  0.233  0.308  0.691-
  22  0.911  0.059  0.275-
  23  0.141  0.132  0.781-
  24  0.657  0.803  0.868-  12 0.51
  25  0.309  0.799  0.056-
  26  0.690  0.833  0.422-
  27  0.493  0.494  0.064-
  28  0.772  0.979  0.887-
  29  0.261  0.655  0.676-
  30  0.783  0.428  0.032-
  31  0.818  0.552  0.871-
  32  0.102  0.215  0.884-
  33  0.401  0.238  0.891-
  34  0.831  0.962  0.608-
  35  0.629  0.553  0.655-
  36  0.998  0.059  0.227-
 
 -----------------------------------------------------------------------------
|                                                                             |
|           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 distance between some ions is very small. Please check the          |
|     nearest-neighbor list in the OUTCAR file.                               |
|     I HOPE YOU KNOW WHAT YOU ARE DOING!                                     |
|                                                                             |
 -----------------------------------------------------------------------------

  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.555423050  0.178776220  0.913678730
     0.640043090  0.502488520  0.950038570
     0.262367520  0.652172480  0.490856970
     0.744942140  0.740959720  0.829449130
     0.360285520  0.143211400  0.733267670
     0.577723330  0.322385970  0.588313590
     0.438976010  0.752928500  0.267029760
     0.300143550  0.105797980  0.935095940
     0.046528080  0.259760190  0.275863640
     0.438134180  0.125942060  0.494854900
     0.984115430  0.805823880  0.373415430
     0.632900420  0.794296480  0.872785170
     0.225578080  0.429074350  0.415319740
     0.050818460  0.388229780  0.168493340
     0.086600430  0.670854350  0.061328580
     0.877305680  0.104321740  0.232693220
     0.226776090  0.483364800  0.328113810
     0.544448970  0.822511100  0.584927940
     0.981319490  0.293932140  0.656392690
     0.388550520  0.981951780  0.006817820
     0.233448830  0.307979880  0.691040230
     0.910873810  0.058684300  0.274810220
     0.140848170  0.131946800  0.780728180
     0.657229510  0.802895920  0.867836960
     0.308693610  0.799121610  0.055700590
     0.689597350  0.832506650  0.422180440
     0.493210320  0.493600980  0.063819320
     0.772255030  0.978838990  0.886570600
     0.260681630  0.655498930  0.676299390
     0.782753040  0.427849830  0.032018100
     0.817746160  0.552033540  0.871495700
     0.101523580  0.215279910  0.883595320
     0.401429740  0.237642060  0.890839510
     0.830621690  0.961737870  0.608135020
     0.628941460  0.553117960  0.655155730
     0.997760520  0.059420350  0.227258390

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

 Generating k-lattice:

  Cartesian coordinates                     Fractional coordinates (reciprocal lattice)
     0.025000000  0.000000000  0.000000000     0.500000000  0.000000000  0.000000000
     0.000000000  0.033333333  0.000000000     0.000000000  0.500000000  0.000000000
     0.000000000  0.000000000  0.033333333     0.000000000  0.000000000  0.500000000

  Length of vectors
     0.025000000  0.033333333  0.033333333

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

 
 Subroutine IBZKPT returns following result:
 ===========================================
 
 Found      8 irreducible k-points:
 
 Following reciprocal coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
  0.500000  0.000000  0.000000      1.000000
  0.000000  0.500000  0.000000      1.000000
  0.000000  0.000000  0.500000      1.000000
  0.500000  0.500000  0.000000      1.000000
  0.500000  0.000000  0.500000      1.000000
  0.000000  0.500000  0.500000      1.000000
  0.500000  0.500000  0.500000      1.000000
 
 Following cartesian coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
  0.025000  0.000000  0.000000      1.000000
  0.000000  0.033333  0.000000      1.000000
  0.000000  0.000000  0.033333      1.000000
  0.025000  0.033333  0.000000      1.000000
  0.025000  0.000000  0.033333      1.000000
  0.000000  0.033333  0.033333      1.000000
  0.025000  0.033333  0.033333      1.000000
 


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




 Dimension of arrays:
   k-points           NKPTS =      8   k-points in BZ     NKDIM =      8   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   =     60;   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       0.125
   0.02500000  0.00000000  0.00000000       0.125
   0.00000000  0.03333333  0.00000000       0.125
   0.00000000  0.00000000  0.03333333       0.125
   0.02500000  0.03333333  0.00000000       0.125
   0.02500000  0.00000000  0.03333333       0.125
   0.00000000  0.03333333  0.03333333       0.125
   0.02500000  0.03333333  0.03333333       0.125
 
 k-points in reciprocal lattice and weights: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       0.125
   0.50000000  0.00000000  0.00000000       0.125
   0.00000000  0.50000000  0.00000000       0.125
   0.00000000  0.00000000  0.50000000       0.125
   0.50000000  0.50000000  0.00000000       0.125
   0.50000000  0.00000000  0.50000000       0.125
   0.00000000  0.50000000  0.50000000       0.125
   0.50000000  0.50000000  0.50000000       0.125
 
 position of ions in fractional coordinates (direct lattice) 
   0.55542305  0.17877622  0.91367873
   0.64004309  0.50248852  0.95003857
   0.26236752  0.65217248  0.49085697
   0.74494214  0.74095972  0.82944913
   0.36028552  0.14321140  0.73326767
   0.57772333  0.32238597  0.58831359
   0.43897601  0.75292850  0.26702976
   0.30014355  0.10579798  0.93509594
   0.04652808  0.25976019  0.27586364
   0.43813418  0.12594206  0.49485490
   0.98411543  0.80582388  0.37341543
   0.63290042  0.79429648  0.87278517
   0.22557808  0.42907435  0.41531974
   0.05081846  0.38822978  0.16849334
   0.08660043  0.67085435  0.06132858
   0.87730568  0.10432174  0.23269322
   0.22677609  0.48336480  0.32811381
   0.54444897  0.82251110  0.58492794
   0.98131949  0.29393214  0.65639269
   0.38855052  0.98195178  0.00681782
   0.23344883  0.30797988  0.69104023
   0.91087381  0.05868430  0.27481022
   0.14084817  0.13194680  0.78072818
   0.65722951  0.80289592  0.86783696
   0.30869361  0.79912161  0.05570059
   0.68959735  0.83250665  0.42218044
   0.49321032  0.49360098  0.06381932
   0.77225503  0.97883899  0.88657060
   0.26068163  0.65549893  0.67629939
   0.78275304  0.42784983  0.03201810
   0.81774616  0.55203354  0.87149570
   0.10152358  0.21527991  0.88359532
   0.40142974  0.23764206  0.89083951
   0.83062169  0.96173787  0.60813502
   0.62894146  0.55311796  0.65515573
   0.99776052  0.05942035  0.22725839
 
 position of ions in cartesian coordinates  (Angst):
  11.10846100  2.68164330 13.70518095
  12.80086180  7.53732780 14.25057855
   5.24735040  9.78258720  7.36285455
  14.89884280 11.11439580 12.44173695
   7.20571040  2.14817100 10.99901505
  11.55446660  4.83578955  8.82470385
   8.77952020 11.29392750  4.00544640
   6.00287100  1.58696970 14.02643910
   0.93056160  3.89640285  4.13795460
   8.76268360  1.88913090  7.42282350
  19.68230860 12.08735820  5.60123145
  12.65800840 11.91444720 13.09177755
   4.51156160  6.43611525  6.22979610
   1.01636920  5.82344670  2.52740010
   1.73200860 10.06281525  0.91992870
  17.54611360  1.56482610  3.49039830
   4.53552180  7.25047200  4.92170715
  10.88897940 12.33766650  8.77391910
  19.62638980  4.40898210  9.84589035
   7.77101040 14.72927670  0.10226730
   4.66897660  4.61969820 10.36560345
  18.21747620  0.88026450  4.12215330
   2.81696340  1.97920200 11.71092270
  13.14459020 12.04343880 13.01755440
   6.17387220 11.98682415  0.83550885
  13.79194700 12.48759975  6.33270660
   9.86420640  7.40401470  0.95728980
  15.44510060 14.68258485 13.29855900
   5.21363260  9.83248395 10.14449085
  15.65506080  6.41774745  0.48027150
  16.35492320  8.28050310 13.07243550
   2.03047160  3.22919865 13.25392980
   8.02859480  3.56463090 13.36259265
  16.61243380 14.42606805  9.12202530
  12.57882920  8.29676940  9.82733595
  19.95521040  0.89130525  3.40887585
 


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


 k-point  1 :   0.0000 0.0000 0.0000  plane waves:   81909
 k-point  2 :   0.5000 0.0000 0.0000  plane waves:   81810
 k-point  3 :   0.0000 0.5000 0.0000  plane waves:   81748
 k-point  4 :   0.0000 0.0000 0.5000  plane waves:   81748
 k-point  5 :   0.5000 0.5000 0.0000  plane waves:   81700
 k-point  6 :   0.5000 0.0000 0.5000  plane waves:   81700
 k-point  7 :   0.0000 0.5000 0.5000  plane waves:   81676
 k-point  8 :   0.5000 0.5000 0.5000  plane waves:   81648

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

 maximum number of plane-waves:     81909
 maximum index in each direction: 
   IXMAX=   32   IYMAX=   24   IZMAX=   24
   IXMIN=  -33   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  1277134. kBytes
=======================================================================

   base      :      30000. kBytes
   nonlr-proj:       4075. kBytes
   fftplans  :      87706. kBytes
   grid      :     314101. kBytes
   one-center:        221. kBytes
   wavefun   :     841031. 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         2539
 Maximum index for augmentation-charges         2241 (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  :  2808
 total energy-change (2. order) : 0.5319703E+03  (-0.1843998E+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  =      -600.23300183
  -Hartree energ DENC   =     -2265.73674255
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       227.93777744
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.00528887
  eigenvalues    EBANDS =      -363.15106451
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       531.97028228 eV

  energy without entropy =      531.97557116  energy(sigma->0) =      531.97204524


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




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



 eigenvalue-minimisations  :  4112
 total energy-change (2. order) :-0.3756463E+03  (-0.3534416E+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  =      -600.23300183
  -Hartree energ DENC   =     -2265.73674255
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       227.93777744
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.00588952
  eigenvalues    EBANDS =      -738.79678305
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       156.32396309 eV

  energy without entropy =      156.32985261  energy(sigma->0) =      156.32592626


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




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



 eigenvalue-minimisations  :  3792
 total energy-change (2. order) :-0.1868730E+03  (-0.1650388E+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  =      -600.23300183
  -Hartree energ DENC   =     -2265.73674255
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       227.93777744
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.11876749
  eigenvalues    EBANDS =      -925.55688571
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -30.54901754 eV

  energy without entropy =      -30.43025005  energy(sigma->0) =      -30.50942838


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




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



 eigenvalue-minimisations  :  2936
 total energy-change (2. order) :-0.2213406E+02  (-0.2019458E+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  =      -600.23300183
  -Hartree energ DENC   =     -2265.73674255
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       227.93777744
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.16972941
  eigenvalues    EBANDS =      -947.63998105
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -52.68307479 eV

  energy without entropy =      -52.51334538  energy(sigma->0) =      -52.62649832


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




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



 eigenvalue-minimisations  :  3264
 total energy-change (2. order) :-0.6344771E+00  (-0.6279978E+00)
 number of electron      87.0000038 magnetization      34.1220739
 augmentation part        3.7148435 magnetization      11.7855444

 Broyden mixing:
  rms(total) = 0.41757E+01    rms(broyden)= 0.41737E+01
  rms(prec ) = 0.50449E+01
  weight for this iteration     100.00

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

  energy without entropy =      -53.13881115  energy(sigma->0) =      -53.25797162


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




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



 eigenvalue-minimisations  :  2872
 total energy-change (2. order) :-0.1034348E+03  (-0.2835205E+02)
 number of electron      87.0000011 magnetization      33.5368440
 augmentation part       -0.3175464 magnetization      12.5679223

 Broyden mixing:
  rms(total) = 0.81649E+01    rms(broyden)= 0.81607E+01
  rms(prec ) = 0.99771E+01
  weight for this iteration     100.00

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2503.77759969
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       284.35877531
  PAW double counting   =      3104.32298754    -3082.84871485
  entropy T*S    EENTRO =        -0.04127227
  eigenvalues    EBANDS =      -868.96174708
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -156.75233106 eV

  energy without entropy =     -156.71105880  energy(sigma->0) =     -156.73857364


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




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



 eigenvalue-minimisations  :  2856
 total energy-change (2. order) : 0.1485998E+03  (-0.9010369E+01)
 number of electron      87.0000041 magnetization      33.1148856
 augmentation part        2.6245632 magnetization      10.5931540

 Broyden mixing:
  rms(total) = 0.37649E+01    rms(broyden)= 0.37595E+01
  rms(prec ) = 0.47311E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3436
  0.4630  0.2243

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2376.36132748
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       272.58680171
  PAW double counting   =      2898.54642888    -2876.37509459
  entropy T*S    EENTRO =         0.04878376
  eigenvalues    EBANDS =      -836.79337430
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -8.15254205 eV

  energy without entropy =       -8.20132581  energy(sigma->0) =       -8.16880330


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




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



 eigenvalue-minimisations  :  2832
 total energy-change (2. order) : 0.1787053E+02  (-0.2326002E+01)
 number of electron      87.0000030 magnetization      33.1504173
 augmentation part        1.4206326 magnetization       2.3821080

 Broyden mixing:
  rms(total) = 0.45002E+01    rms(broyden)= 0.44981E+01
  rms(prec ) = 0.54183E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2564
  0.3253  0.3253  0.1187

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2378.61038873
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       268.93718710
  PAW double counting   =      2903.04431697    -2881.21360357
  entropy T*S    EENTRO =        -0.20490126
  eigenvalues    EBANDS =      -812.42986277
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =         9.71798773 eV

  energy without entropy =        9.92288898  energy(sigma->0) =        9.78628815


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




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



 eigenvalue-minimisations  :  3176
 total energy-change (2. order) :-0.3934504E+00  (-0.2062185E+01)
 number of electron      87.0000053 magnetization      33.2154731
 augmentation part        1.8888410 magnetization      -1.0192665

 Broyden mixing:
  rms(total) = 0.38957E+01    rms(broyden)= 0.38928E+01
  rms(prec ) = 0.45111E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2080
  0.3146  0.3146  0.1567  0.0461

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2376.51184537
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       266.32505262
  PAW double counting   =      2902.43164264    -2880.55060680
  entropy T*S    EENTRO =        -0.07669640
  eigenvalues    EBANDS =      -812.48824936
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =         9.32453731 eV

  energy without entropy =        9.40123371  energy(sigma->0) =        9.35010278


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




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



 eigenvalue-minimisations  :  3152
 total energy-change (2. order) :-0.1154070E+01  (-0.1162094E+01)
 number of electron      87.0000023 magnetization      33.3094998
 augmentation part        1.9029504 magnetization      -2.6849261

 Broyden mixing:
  rms(total) = 0.36347E+01    rms(broyden)= 0.36330E+01
  rms(prec ) = 0.40765E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1699
  0.2864  0.2864  0.1625  0.0571  0.0571

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2375.88435454
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       264.26103270
  PAW double counting   =      2902.77242790    -2880.85032668
  entropy T*S    EENTRO =        -0.10042648
  eigenvalues    EBANDS =      -812.22312584
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =         8.17046704 eV

  energy without entropy =        8.27089352  energy(sigma->0) =        8.20394253


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




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



 eigenvalue-minimisations  :  3176
 total energy-change (2. order) :-0.2779768E+01  (-0.2606593E+00)
 number of electron      87.0000050 magnetization      33.4767021
 augmentation part        1.6277274 magnetization      -4.4058497

 Broyden mixing:
  rms(total) = 0.40098E+01    rms(broyden)= 0.40068E+01
  rms(prec ) = 0.44878E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1693
  0.2777  0.2777  0.1816  0.1816  0.0611  0.0362

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2375.48597207
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       261.90875568
  PAW double counting   =      2903.11190878    -2881.14905739
  entropy T*S    EENTRO =         0.07116476
  eigenvalues    EBANDS =      -813.26134054
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =         5.39069919 eV

  energy without entropy =        5.31953443  energy(sigma->0) =        5.36697760


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




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



 eigenvalue-minimisations  :  3288
 total energy-change (2. order) :-0.8315098E+01  (-0.1119954E+01)
 number of electron      87.0000040 magnetization      33.6233132
 augmentation part        1.7599923 magnetization      -3.6915422

 Broyden mixing:
  rms(total) = 0.38431E+01    rms(broyden)= 0.38407E+01
  rms(prec ) = 0.43639E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1626
  0.3280  0.2135  0.2135  0.2407  0.0676  0.0374  0.0374

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.58137150
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       258.50381695
  PAW double counting   =      2900.00983055    -2877.99082891
  entropy T*S    EENTRO =        -0.04338954
  eigenvalues    EBANDS =      -820.01769672
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -2.92439919 eV

  energy without entropy =       -2.88100965  energy(sigma->0) =       -2.90993601


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




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



 eigenvalue-minimisations  :  3072
 total energy-change (2. order) :-0.1958730E+01  (-0.1888661E+00)
 number of electron      87.0000050 magnetization      33.6417312
 augmentation part        1.9683885 magnetization      -2.8992188

 Broyden mixing:
  rms(total) = 0.33461E+01    rms(broyden)= 0.33450E+01
  rms(prec ) = 0.37685E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1545
  0.3310  0.2125  0.2125  0.2375  0.0743  0.0743  0.0501  0.0440

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.80189347
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.32494568
  PAW double counting   =      2898.03325308    -2875.99384362
  entropy T*S    EENTRO =        -0.02115452
  eigenvalues    EBANDS =      -821.61967609
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -4.88312896 eV

  energy without entropy =       -4.86197444  energy(sigma->0) =       -4.87607745


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




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



 eigenvalue-minimisations  :  3088
 total energy-change (2. order) :-0.2377297E+00  (-0.7164365E-02)
 number of electron      87.0000050 magnetization      33.5296508
 augmentation part        1.9438505 magnetization      -3.3653651

 Broyden mixing:
  rms(total) = 0.33676E+01    rms(broyden)= 0.33676E+01
  rms(prec ) = 0.37942E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1885
  0.3519  0.2850  0.2850  0.2349  0.2352  0.1071  0.1071  0.0451  0.0451

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.35249444
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.23788428
  PAW double counting   =      2898.96984006    -2876.93703151
  entropy T*S    EENTRO =        -0.00481390
  eigenvalues    EBANDS =      -821.22948313
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -5.12085866 eV

  energy without entropy =       -5.11604476  energy(sigma->0) =       -5.11925403


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




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



 eigenvalue-minimisations  :  2856
 total energy-change (2. order) : 0.1148412E+01  (-0.3550980E-01)
 number of electron      87.0000042 magnetization      33.5386569
 augmentation part        1.8633949 magnetization      -3.7491030

 Broyden mixing:
  rms(total) = 0.36666E+01    rms(broyden)= 0.36661E+01
  rms(prec ) = 0.41287E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1722
  0.3586  0.2805  0.2805  0.2324  0.1302  0.1302  0.1232  0.0965  0.0450  0.0450

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.70754451
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.71472477
  PAW double counting   =      2898.12301445    -2876.09419574
  entropy T*S    EENTRO =        -0.00395726
  eigenvalues    EBANDS =      -821.19972787
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -3.97244619 eV

  energy without entropy =       -3.96848893  energy(sigma->0) =       -3.97112710


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




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



 eigenvalue-minimisations  :  3000
 total energy-change (2. order) :-0.2575134E+00  (-0.3426539E-02)
 number of electron      87.0000043 magnetization      33.5728255
 augmentation part        1.8737446 magnetization      -3.4120053

 Broyden mixing:
  rms(total) = 0.36442E+01    rms(broyden)= 0.36442E+01
  rms(prec ) = 0.41153E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2118
  0.3384  0.3384  0.3115  0.3115  0.3449  0.2218  0.1441  0.1441  0.0852  0.0449
  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.77508341
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.58294969
  PAW double counting   =      2898.19932474    -2876.17289016
  entropy T*S    EENTRO =        -0.02171543
  eigenvalues    EBANDS =      -821.23778502
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -4.22995962 eV

  energy without entropy =       -4.20824419  energy(sigma->0) =       -4.22272114


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




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



 eigenvalue-minimisations  :  2880
 total energy-change (2. order) : 0.1770083E+00  (-0.4526010E-01)
 number of electron      87.0000044 magnetization      33.7017836
 augmentation part        1.8078722 magnetization      -4.2893903

 Broyden mixing:
  rms(total) = 0.37442E+01    rms(broyden)= 0.37440E+01
  rms(prec ) = 0.41634E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2295
  0.4453  0.4453  0.3377  0.3377  0.3422  0.2577  0.1460  0.1460  0.1184  0.0875
  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.63104123
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.35392926
  PAW double counting   =      2898.66984214    -2876.63623379
  entropy T*S    EENTRO =         0.08110183
  eigenvalues    EBANDS =      -821.08578953
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -4.05295135 eV

  energy without entropy =       -4.13405318  energy(sigma->0) =       -4.07998529


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




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



 eigenvalue-minimisations  :  3024
 total energy-change (2. order) :-0.1390976E+01  (-0.6110206E-01)
 number of electron      87.0000053 magnetization      33.6958298
 augmentation part        1.8719030 magnetization      -4.1676882

 Broyden mixing:
  rms(total) = 0.36478E+01    rms(broyden)= 0.36474E+01
  rms(prec ) = 0.40637E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2309
  0.4843  0.4843  0.3381  0.3381  0.3524  0.2253  0.1564  0.1564  0.1450  0.1450
  0.0860  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.34992887
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.93189641
  PAW double counting   =      2900.62339935    -2878.59464058
  entropy T*S    EENTRO =         0.13532531
  eigenvalues    EBANDS =      -821.38521922
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -5.44392763 eV

  energy without entropy =       -5.57925294  energy(sigma->0) =       -5.48903607


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




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



 eigenvalue-minimisations  :  3056
 total energy-change (2. order) : 0.7181497E-01  (-0.2808477E-02)
 number of electron      87.0000053 magnetization      33.8862566
 augmentation part        1.8817785 magnetization      -3.9781324

 Broyden mixing:
  rms(total) = 0.36515E+01    rms(broyden)= 0.36515E+01
  rms(prec ) = 0.40659E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2819
  0.6388  0.6388  0.3805  0.3805  0.3071  0.3071  0.3176  0.3176  0.1777  0.1777
  0.1273  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.71851814
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.93518438
  PAW double counting   =      2899.43112324    -2877.39728490
  entropy T*S    EENTRO =         0.13653590
  eigenvalues    EBANDS =      -821.95439311
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -5.37211266 eV

  energy without entropy =       -5.50864856  energy(sigma->0) =       -5.41762463


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




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



 eigenvalue-minimisations  :  2888
 total energy-change (2. order) :-0.2409772E+01  (-0.3846579E-01)
 number of electron      87.0000047 magnetization      33.9099388
 augmentation part        1.8320007 magnetization      -3.5380362

 Broyden mixing:
  rms(total) = 0.35822E+01    rms(broyden)= 0.35817E+01
  rms(prec ) = 0.40361E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2812
  0.6759  0.6759  0.3850  0.3850  0.3223  0.3223  0.3068  0.2411  0.2411  0.1794
  0.1794  0.1285  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.50006657
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.79406372
  PAW double counting   =      2899.78789168    -2877.76014061
  entropy T*S    EENTRO =         0.14981046
  eigenvalues    EBANDS =      -821.44868368
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.78188502 eV

  energy without entropy =       -7.93169548  energy(sigma->0) =       -7.83182184


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




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



 eigenvalue-minimisations  :  3280
 total energy-change (2. order) :-0.1309229E+00  (-0.5679685E-02)
 number of electron      87.0000046 magnetization      33.8764172
 augmentation part        1.8191295 magnetization      -3.5199104

 Broyden mixing:
  rms(total) = 0.35288E+01    rms(broyden)= 0.35287E+01
  rms(prec ) = 0.39738E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2698
  0.6747  0.6747  0.3899  0.3899  0.3285  0.3285  0.3219  0.2481  0.2481  0.1823
  0.1823  0.1281  0.0864  0.0449  0.0449  0.0433

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.57361080
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.64846668
  PAW double counting   =      2900.03281455    -2878.00215499
  entropy T*S    EENTRO =         0.12067187
  eigenvalues    EBANDS =      -821.33423518
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.91280789 eV

  energy without entropy =       -8.03347976  energy(sigma->0) =       -7.95303185


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




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



 eigenvalue-minimisations  :  2896
 total energy-change (2. order) : 0.1615353E+00  (-0.1000061E+00)
 number of electron      87.0000046 magnetization      33.7810194
 augmentation part        1.7872331 magnetization      -2.9067381

 Broyden mixing:
  rms(total) = 0.34062E+01    rms(broyden)= 0.34059E+01
  rms(prec ) = 0.38154E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2877
  0.6479  0.6479  0.4043  0.4043  0.3432  0.3432  0.2542  0.2542  0.3438  0.2878
  0.2878  0.1840  0.1840  0.1278  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.76517226
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.73797141
  PAW double counting   =      2900.81900275    -2878.80149174
  entropy T*S    EENTRO =         0.04875234
  eigenvalues    EBANDS =      -820.98557507
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.75127260 eV

  energy without entropy =       -7.80002494  energy(sigma->0) =       -7.76752338


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




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



 eigenvalue-minimisations  :  3008
 total energy-change (2. order) : 0.2448146E+00  (-0.1796339E-01)
 number of electron      87.0000050 magnetization      33.7642652
 augmentation part        1.8455213 magnetization      -2.3401659

 Broyden mixing:
  rms(total) = 0.35136E+01    rms(broyden)= 0.35134E+01
  rms(prec ) = 0.39551E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3056
  0.6424  0.6424  0.4167  0.4167  0.3477  0.3477  0.3625  0.3625  0.3471  0.3471
  0.2981  0.2981  0.1841  0.1841  0.1277  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.57974725
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.83689330
  PAW double counting   =      2899.33299136    -2877.30912270
  entropy T*S    EENTRO =         0.09302993
  eigenvalues    EBANDS =      -822.07574262
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.50645800 eV

  energy without entropy =       -7.59948793  energy(sigma->0) =       -7.53746798


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




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



 eigenvalue-minimisations  :  3312
 total energy-change (2. order) : 0.1576451E+00  (-0.1831912E-02)
 number of electron      87.0000050 magnetization      33.6595084
 augmentation part        1.8462020 magnetization      -2.4634053

 Broyden mixing:
  rms(total) = 0.35208E+01    rms(broyden)= 0.35208E+01
  rms(prec ) = 0.39619E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3305
  0.5736  0.5736  0.5383  0.5383  0.4388  0.4388  0.3981  0.3981  0.3447  0.3447
  0.3642  0.3279  0.3279  0.1840  0.1840  0.1277  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.67118197
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.91243546
  PAW double counting   =      2899.57867967    -2877.55694079
  entropy T*S    EENTRO =         0.09943320
  eigenvalues    EBANDS =      -821.90647848
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.34881294 eV

  energy without entropy =       -7.44824614  energy(sigma->0) =       -7.38195734


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




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



 eigenvalue-minimisations  :  2880
 total energy-change (2. order) : 0.6833455E+00  (-0.4823571E-02)
 number of electron      87.0000049 magnetization      34.0839672
 augmentation part        1.8409595 magnetization      -2.2852513

 Broyden mixing:
  rms(total) = 0.35581E+01    rms(broyden)= 0.35580E+01
  rms(prec ) = 0.39951E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3494
  0.7081  0.7081  0.4446  0.4231  0.4231  0.4669  0.4669  0.4189  0.4189  0.4009
  0.3709  0.3709  0.3473  0.3473  0.1840  0.1840  0.1277  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.76852043
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.24752071
  PAW double counting   =      2901.28219705    -2879.26513730
  entropy T*S    EENTRO =         0.08955867
  eigenvalues    EBANDS =      -821.44632608
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -6.66546740 eV

  energy without entropy =       -6.75502607  energy(sigma->0) =       -6.69532029


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




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



 eigenvalue-minimisations  :  2784
 total energy-change (2. order) :-0.2134356E+01  (-0.7414044E-01)
 number of electron      87.0000042 magnetization      34.3008302
 augmentation part        1.6938429 magnetization      -3.5965721

 Broyden mixing:
  rms(total) = 0.33185E+01    rms(broyden)= 0.33178E+01
  rms(prec ) = 0.37072E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3715
  0.7963  0.7963  0.5274  0.4759  0.4759  0.5035  0.5035  0.4294  0.4294  0.3935
  0.3935  0.3746  0.3746  0.3275  0.3275  0.1840  0.1840  0.1277  0.0864  0.0449
  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2375.72839128
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.09013179
  PAW double counting   =      2901.74560613    -2879.71245765
  entropy T*S    EENTRO =        -0.02010462
  eigenvalues    EBANDS =      -820.36984798
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -8.79982363 eV

  energy without entropy =       -8.77971902  energy(sigma->0) =       -8.79312209


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




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



 eigenvalue-minimisations  :  2880
 total energy-change (2. order) :-0.1626346E+01  (-0.8590433E-01)
 number of electron      87.0000041 magnetization      34.2703330
 augmentation part        1.7423347 magnetization      -2.8485982

 Broyden mixing:
  rms(total) = 0.30388E+01    rms(broyden)= 0.30379E+01
  rms(prec ) = 0.33803E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3613
  0.7773  0.7773  0.4968  0.4663  0.4663  0.5030  0.5030  0.4343  0.4343  0.4051
  0.4051  0.3920  0.3920  0.3310  0.3310  0.1840  0.1840  0.1277  0.0864  0.0449
  0.0449  0.1614

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.65879494
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.91796407
  PAW double counting   =      2898.24129697    -2876.15414772
  entropy T*S    EENTRO =        -0.00576177
  eigenvalues    EBANDS =      -823.96196582
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -10.42616923 eV

  energy without entropy =      -10.42040746  energy(sigma->0) =      -10.42424864


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




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



 eigenvalue-minimisations  :  3032
 total energy-change (2. order) : 0.7440001E-01  (-0.1595700E-01)
 number of electron      87.0000042 magnetization      34.4477428
 augmentation part        1.7534256 magnetization      -2.6970438

 Broyden mixing:
  rms(total) = 0.30297E+01    rms(broyden)= 0.30297E+01
  rms(prec ) = 0.33607E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3709
  0.8044  0.8044  0.4625  0.4625  0.5239  0.5239  0.4605  0.4605  0.4356  0.4356
  0.4048  0.4048  0.3663  0.3663  0.3284  0.3284  0.2858  0.1840  0.1840  0.1277
  0.0864  0.0449  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.63627741
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.90294289
  PAW double counting   =      2898.61761429    -2876.53207392
  entropy T*S    EENTRO =         0.00562527
  eigenvalues    EBANDS =      -823.90484033
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -10.35176923 eV

  energy without entropy =      -10.35739450  energy(sigma->0) =      -10.35364432


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




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



 eigenvalue-minimisations  :  2952
 total energy-change (2. order) :-0.1598374E+01  (-0.7376409E-01)
 number of electron      87.0000047 magnetization      34.4586582
 augmentation part        1.5775264 magnetization      -2.8201805

 Broyden mixing:
  rms(total) = 0.29337E+01    rms(broyden)= 0.29324E+01
  rms(prec ) = 0.33246E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3604
  0.8093  0.8093  0.5039  0.5242  0.5242  0.4609  0.4609  0.4329  0.4329  0.4144
  0.4144  0.3839  0.3839  0.3306  0.3306  0.2802  0.1840  0.1840  0.1277  0.0864
  0.0449  0.0449  0.2410  0.2410

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.80554449
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.49686012
  PAW double counting   =      2899.36217002    -2877.26145591
  entropy T*S    EENTRO =        -0.14263743
  eigenvalues    EBANDS =      -823.79477553
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -11.95014325 eV

  energy without entropy =      -11.80750582  energy(sigma->0) =      -11.90259744


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




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



 eigenvalue-minimisations  :  2904
 total energy-change (2. order) : 0.5189613E-01  (-0.1837029E-01)
 number of electron      87.0000050 magnetization      34.7082402
 augmentation part        1.5376136 magnetization      -2.4167511

 Broyden mixing:
  rms(total) = 0.31192E+01    rms(broyden)= 0.31189E+01
  rms(prec ) = 0.35720E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3556
  0.7964  0.7964  0.4543  0.4543  0.4398  0.4398  0.4802  0.4802  0.4064  0.4064
  0.4183  0.4183  0.4076  0.4076  0.3329  0.3329  0.2732  0.2732  0.1840  0.1840
  0.0449  0.0449  0.0864  0.1277  0.1996

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.88625151
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.63506057
  PAW double counting   =      2898.84868517    -2876.75152409
  entropy T*S    EENTRO =        -0.14447153
  eigenvalues    EBANDS =      -823.79498571
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -11.89824712 eV

  energy without entropy =      -11.75377559  energy(sigma->0) =      -11.85008994


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




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



 eigenvalue-minimisations  :  3096
 total energy-change (2. order) :-0.4271950E+01  (-0.5034043E+00)
 number of electron      87.0000056 magnetization      34.7263293
 augmentation part        1.5547519 magnetization      -1.0061691

 Broyden mixing:
  rms(total) = 0.32189E+01    rms(broyden)= 0.32177E+01
  rms(prec ) = 0.37721E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3518
  0.7906  0.7906  0.4719  0.4719  0.4168  0.4168  0.4596  0.4596  0.3686  0.3686
  0.4057  0.4057  0.4217  0.4217  0.3989  0.3989  0.3322  0.3322  0.2635  0.1840
  0.1840  0.1277  0.0449  0.0449  0.0864  0.0791

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.50273887
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.40792129
  PAW double counting   =      2894.23570356    -2872.12648749
  entropy T*S    EENTRO =        -0.09523679
  eigenvalues    EBANDS =      -827.28459880
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -16.17019713 eV

  energy without entropy =      -16.07496034  energy(sigma->0) =      -16.13845153


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




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



 eigenvalue-minimisations  :  3088
 total energy-change (2. order) :-0.3072277E+00  (-0.2171359E-01)
 number of electron      87.0000057 magnetization      34.9336996
 augmentation part        1.5498831 magnetization      -0.1164497

 Broyden mixing:
  rms(total) = 0.32861E+01    rms(broyden)= 0.32859E+01
  rms(prec ) = 0.38553E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3577
  0.7920  0.7920  0.4736  0.4736  0.4399  0.4399  0.4190  0.4583  0.4583  0.3243
  0.3243  0.4136  0.4136  0.4262  0.4262  0.4030  0.4030  0.3340  0.3340  0.3184
  0.1840  0.1840  0.1277  0.0864  0.0449  0.0449  0.1181

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.43393492
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.51921563
  PAW double counting   =      2893.84986958    -2871.73437212
  entropy T*S    EENTRO =        -0.03750867
  eigenvalues    EBANDS =      -827.83593435
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -16.47742487 eV

  energy without entropy =      -16.43991621  energy(sigma->0) =      -16.46492198


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




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



 eigenvalue-minimisations  :  2824
 total energy-change (2. order) :-0.1406830E+01  (-0.1171907E+00)
 number of electron      87.0000063 magnetization      34.9234889
 augmentation part        1.4986129 magnetization      -1.2101057

 Broyden mixing:
  rms(total) = 0.34414E+01    rms(broyden)= 0.34407E+01
  rms(prec ) = 0.40637E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3559
  0.8184  0.8184  0.4711  0.4416  0.4416  0.4560  0.4560  0.4490  0.4490  0.4172
  0.4172  0.4354  0.4354  0.4076  0.4076  0.3803  0.3351  0.3351  0.3054  0.2463
  0.2463  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.1220

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.89276655
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.51321999
  PAW double counting   =      2891.99514074    -2869.85629393
  entropy T*S    EENTRO =        -0.00675225
  eigenvalues    EBANDS =      -828.83204305
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -17.88425507 eV

  energy without entropy =      -17.87750282  energy(sigma->0) =      -17.88200432


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




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



 eigenvalue-minimisations  :  2744
 total energy-change (2. order) :-0.7797443E+00  (-0.7506135E-01)
 number of electron      87.0000062 magnetization      35.7568772
 augmentation part        1.4437527 magnetization       0.6971697

 Broyden mixing:
  rms(total) = 0.35794E+01    rms(broyden)= 0.35789E+01
  rms(prec ) = 0.42789E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3720
  0.9107  0.9107  0.5381  0.5408  0.5408  0.4701  0.4701  0.4391  0.4391  0.4383
  0.4383  0.4521  0.4521  0.4279  0.4279  0.3761  0.3370  0.3370  0.2997  0.2997
  0.2924  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.1524  0.1255

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.80938531
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.14845578
  PAW double counting   =      2891.40509659    -2869.27141560
  entropy T*S    EENTRO =        -0.03778427
  eigenvalues    EBANDS =      -829.29420659
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -18.66399942 eV

  energy without entropy =      -18.62621515  energy(sigma->0) =      -18.65140466


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




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



 eigenvalue-minimisations  :  2800
 total energy-change (2. order) :-0.3498622E+01  (-0.1017807E+01)
 number of electron      87.0000063 magnetization      36.0232013
 augmentation part        1.2887304 magnetization       0.5476663

 Broyden mixing:
  rms(total) = 0.36125E+01    rms(broyden)= 0.36100E+01
  rms(prec ) = 0.43989E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3693
  0.9397  0.9397  0.5418  0.5782  0.5782  0.4716  0.4716  0.4421  0.4421  0.4113
  0.4113  0.4479  0.4479  0.3046  0.3046  0.2951  0.4208  0.4208  0.3935  0.3371
  0.3371  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.1736  0.1736  0.1225

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.20931172
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       251.82025740
  PAW double counting   =      2882.72862732    -2860.47143491
  entropy T*S    EENTRO =        -0.08079320
  eigenvalues    EBANDS =      -832.14520631
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -22.16262145 eV

  energy without entropy =      -22.08182825  energy(sigma->0) =      -22.13569038


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




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



 eigenvalue-minimisations  :  3200
 total energy-change (2. order) :-0.1808900E+00  (-0.3152899E+00)
 number of electron      87.0000062 magnetization      36.0478021
 augmentation part        1.4608245 magnetization       0.7251332

 Broyden mixing:
  rms(total) = 0.32455E+01    rms(broyden)= 0.32450E+01
  rms(prec ) = 0.39494E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3704
  0.9610  0.9610  0.5555  0.6069  0.6069  0.4752  0.4752  0.3121  0.3121  0.2857
  0.4369  0.4369  0.3856  0.3856  0.4276  0.4276  0.4044  0.4044  0.3825  0.3350
  0.3350  0.3044  0.3044  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.1641
  0.1230

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.85470526
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       251.64241893
  PAW double counting   =      2883.41040001    -2861.12076567
  entropy T*S    EENTRO =        -0.08139105
  eigenvalues    EBANDS =      -832.53470842
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -22.34351149 eV

  energy without entropy =      -22.26212043  energy(sigma->0) =      -22.31638114


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




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



 eigenvalue-minimisations  :  3264
 total energy-change (2. order) : 0.4787147E-01  (-0.2142976E+00)
 number of electron      87.0000057 magnetization      36.1973257
 augmentation part        1.5676865 magnetization       0.1125958

 Broyden mixing:
  rms(total) = 0.27679E+01    rms(broyden)= 0.27664E+01
  rms(prec ) = 0.33037E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3610
  0.9650  0.9650  0.5576  0.6085  0.6085  0.4759  0.4759  0.4438  0.4438  0.3124
  0.3124  0.2840  0.4275  0.4275  0.3716  0.3716  0.4077  0.4077  0.3870  0.3353
  0.3353  0.3095  0.3095  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.1570
  0.1232  0.0576

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.30802113
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       250.97412031
  PAW double counting   =      2886.43823935    -2864.15854009
  entropy T*S    EENTRO =        -0.09698666
  eigenvalues    EBANDS =      -831.33969178
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -22.29564002 eV

  energy without entropy =      -22.19865336  energy(sigma->0) =      -22.26331113


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




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



 eigenvalue-minimisations  :  2992
 total energy-change (2. order) : 0.6618733E+00  (-0.2670258E+00)
 number of electron      87.0000057 magnetization      36.2927249
 augmentation part        1.6565850 magnetization       1.0159462

 Broyden mixing:
  rms(total) = 0.26618E+01    rms(broyden)= 0.26608E+01
  rms(prec ) = 0.31679E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3679
  0.9740  0.9740  0.5573  0.6000  0.6000  0.4758  0.4758  0.4460  0.4460  0.4242
  0.4242  0.4346  0.4346  0.3144  0.3144  0.2853  0.4059  0.4059  0.3769  0.3360
  0.3360  0.2816  0.2816  0.2808  0.2808  0.1840  0.1840  0.0449  0.0449  0.0864
  0.1277  0.1808  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.80340910
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       251.54303615
  PAW double counting   =      2882.45355871    -2860.14604068
  entropy T*S    EENTRO =        -0.03803985
  eigenvalues    EBANDS =      -831.83811190
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -21.63376669 eV

  energy without entropy =      -21.59572684  energy(sigma->0) =      -21.62108674


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




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



 eigenvalue-minimisations  :  3464
 total energy-change (2. order) : 0.1004730E+01  (-0.1177899E+00)
 number of electron      87.0000056 magnetization      36.2961764
 augmentation part        1.7098671 magnetization       0.5244201

 Broyden mixing:
  rms(total) = 0.25061E+01    rms(broyden)= 0.25054E+01
  rms(prec ) = 0.29847E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3645
  0.9726  0.9726  0.5569  0.6113  0.6113  0.4763  0.4763  0.4560  0.4560  0.4253
  0.4253  0.4361  0.4361  0.3141  0.3141  0.2858  0.4083  0.4083  0.3661  0.3376
  0.3376  0.2990  0.2990  0.2846  0.2846  0.1840  0.1840  0.0449  0.0449  0.0864
  0.1277  0.1894  0.1568  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.61742191
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.32517098
  PAW double counting   =      2876.74291776    -2854.40459219
  entropy T*S    EENTRO =         0.02246451
  eigenvalues    EBANDS =      -831.89281629
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -20.62903717 eV

  energy without entropy =      -20.65150168  energy(sigma->0) =      -20.63652534


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




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



 eigenvalue-minimisations  :  3088
 total energy-change (2. order) :-0.1729601E-02  (-0.6142730E-02)
 number of electron      87.0000056 magnetization      36.3486665
 augmentation part        1.7186225 magnetization       0.5671463

 Broyden mixing:
  rms(total) = 0.24870E+01    rms(broyden)= 0.24870E+01
  rms(prec ) = 0.29627E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3827
  0.9590  0.9590  0.6480  0.6480  0.5357  0.5357  0.4778  0.4778  0.4800  0.4800
  0.4556  0.4556  0.4377  0.4377  0.3153  0.3153  0.2821  0.3406  0.3406  0.4145
  0.4145  0.4018  0.3618  0.3618  0.3354  0.3354  0.1840  0.1840  0.0449  0.0449
  0.0864  0.1277  0.1972  0.1972  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2373.86083517
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.44918337
  PAW double counting   =      2877.38758124    -2855.05333078
  entropy T*S    EENTRO =         0.02250782
  eigenvalues    EBANDS =      -831.77111322
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -20.63076677 eV

  energy without entropy =      -20.65327459  energy(sigma->0) =      -20.63826938


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




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



 eigenvalue-minimisations  :  2688
 total energy-change (2. order) :-0.1760941E+00  (-0.1695098E+00)
 number of electron      87.0000056 magnetization      36.5094765
 augmentation part        1.6920624 magnetization       0.2133238

 Broyden mixing:
  rms(total) = 0.25125E+01    rms(broyden)= 0.25123E+01
  rms(prec ) = 0.29960E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4023
  0.9554  0.9554  0.7504  0.7155  0.7155  0.5947  0.4792  0.4792  0.5260  0.5260
  0.4799  0.4799  0.3155  0.3155  0.2826  0.3514  0.3514  0.4284  0.4284  0.4191
  0.4191  0.4239  0.4239  0.3942  0.3942  0.3369  0.3369  0.1840  0.1840  0.0449
  0.0449  0.0864  0.1277  0.2054  0.2054  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.27811650
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.41225952
  PAW double counting   =      2874.39285109    -2852.01914604
  entropy T*S    EENTRO =         0.00416675
  eigenvalues    EBANDS =      -833.51411562
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -20.80686083 eV

  energy without entropy =      -20.81102758  energy(sigma->0) =      -20.80824975


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




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



 eigenvalue-minimisations  :  3000
 total energy-change (2. order) :-0.9154263E+00  (-0.1028583E+01)
 number of electron      87.0000057 magnetization      36.6456309
 augmentation part        1.6667368 magnetization       0.2103301

 Broyden mixing:
  rms(total) = 0.26978E+01    rms(broyden)= 0.26950E+01
  rms(prec ) = 0.32680E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4137
  1.0190  1.0190  0.7953  0.7531  0.7531  0.5859  0.5747  0.5747  0.4787  0.4787
  0.4735  0.4735  0.3157  0.3157  0.2828  0.4626  0.4626  0.4362  0.4362  0.3488
  0.3488  0.4285  0.4285  0.4176  0.3831  0.3831  0.3371  0.3371  0.1840  0.1840
  0.0449  0.0449  0.0864  0.1277  0.2035  0.2035  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2370.41573506
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.39551790
  PAW double counting   =      2868.55760211    -2846.14212694
  entropy T*S    EENTRO =         0.01528716
  eigenvalues    EBANDS =      -837.32807232
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -21.72228718 eV

  energy without entropy =      -21.73757434  energy(sigma->0) =      -21.72738290


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




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



 eigenvalue-minimisations  :  2928
 total energy-change (2. order) :-0.2913292E+01  (-0.4719936E+00)
 number of electron      87.0000049 magnetization      37.6760462
 augmentation part        1.6723210 magnetization       0.9226950

 Broyden mixing:
  rms(total) = 0.26635E+01    rms(broyden)= 0.26629E+01
  rms(prec ) = 0.32275E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4373
  1.0667  1.0667  0.8404  0.8030  0.8030  0.5804  0.7433  0.7433  0.4779  0.4779
  0.4795  0.4795  0.5259  0.5259  0.3157  0.3157  0.2829  0.3484  0.3484  0.4644
  0.4644  0.4379  0.4379  0.4439  0.4439  0.4113  0.4113  0.3369  0.3369  0.1840
  0.1840  0.0449  0.0449  0.0864  0.1277  0.2039  0.2039  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.28879567
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       251.74554873
  PAW double counting   =      2877.40516307    -2854.99320343
  entropy T*S    EENTRO =         0.04296338
  eigenvalues    EBANDS =      -837.74249493
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -24.63557888 eV

  energy without entropy =      -24.67854226  energy(sigma->0) =      -24.64990000


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




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



 eigenvalue-minimisations  :  2944
 total energy-change (2. order) : 0.3037038E+01  (-0.1051004E+01)
 number of electron      87.0000028 magnetization      38.3037695
 augmentation part        1.4310444 magnetization      -0.1505868

 Broyden mixing:
  rms(total) = 0.27086E+01    rms(broyden)= 0.27066E+01
  rms(prec ) = 0.32505E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4362
  1.0302  1.0302  0.8624  0.8097  0.8097  0.8115  0.8115  0.5792  0.4778  0.4778
  0.5246  0.5246  0.4860  0.4860  0.3158  0.3158  0.2829  0.3484  0.3484  0.4452
  0.4452  0.4498  0.4498  0.4051  0.4051  0.4064  0.3979  0.3979  0.3370  0.3370
  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.2038  0.2038  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2374.95867925
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       255.87306782
  PAW double counting   =      2881.14050441    -2858.66715261
  entropy T*S    EENTRO =        -0.03476179
  eigenvalues    EBANDS =      -835.14675971
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -21.59854115 eV

  energy without entropy =      -21.56377936  energy(sigma->0) =      -21.58695389


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




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



 eigenvalue-minimisations  :  3144
 total energy-change (2. order) : 0.1056671E+02  (-0.5911203E+01)
 number of electron      87.0000028 magnetization      38.6736448
 augmentation part        1.6777325 magnetization      -0.3158655

 Broyden mixing:
  rms(total) = 0.22412E+01    rms(broyden)= 0.22371E+01
  rms(prec ) = 0.26399E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4296
  1.0341  1.0341  0.8661  0.8151  0.8151  0.7944  0.7944  0.5790  0.4778  0.4778
  0.4898  0.4898  0.5112  0.5112  0.3158  0.3158  0.2829  0.3483  0.3483  0.4369
  0.4369  0.4386  0.4386  0.4097  0.4097  0.4102  0.4102  0.3372  0.3372  0.3067
  0.3067  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.2039  0.2039  0.1229

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2372.32076329
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       261.68736836
  PAW double counting   =      2872.37010259    -2849.81619053
  entropy T*S    EENTRO =        -0.00189994
  eigenvalues    EBANDS =      -833.14568633
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -11.03182918 eV

  energy without entropy =      -11.02992924  energy(sigma->0) =      -11.03119586


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




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



 eigenvalue-minimisations  :  3168
 total energy-change (2. order) : 0.3344224E+01  (-0.1169167E+01)
 number of electron      87.0000030 magnetization      38.6841272
 augmentation part        1.7566169 magnetization       0.2748933

 Broyden mixing:
  rms(total) = 0.30150E+01    rms(broyden)= 0.30146E+01
  rms(prec ) = 0.33123E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4191
  1.0336  1.0336  0.8661  0.8153  0.8153  0.7945  0.7945  0.5790  0.4778  0.4778
  0.4898  0.4898  0.5108  0.5108  0.3158  0.3158  0.2829  0.3483  0.3483  0.4369
  0.4369  0.4378  0.4378  0.4105  0.4105  0.4119  0.4119  0.3372  0.3372  0.3015
  0.3015  0.1840  0.1840  0.0449  0.0449  0.0864  0.1277  0.2039  0.2039  0.1229
  0.0104

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2370.81350100
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       266.70757953
  PAW double counting   =      2858.35507326    -2836.12873519
  entropy T*S    EENTRO =         0.04167839
  eigenvalues    EBANDS =      -836.04494062
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -7.68760566 eV

  energy without entropy =       -7.72928405  energy(sigma->0) =       -7.70149846


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




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



 eigenvalue-minimisations  :  3232
 total energy-change (2. order) :-0.3560053E+00  (-0.3703892E-01)
 number of electron      87.0000037 magnetization      38.7085048
 augmentation part        1.7952135 magnetization       1.9091627

 Broyden mixing:
  rms(total) = 0.38404E+01    rms(broyden)= 0.38396E+01
  rms(prec ) = 0.41397E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4085
  1.0192  1.0192  0.8656  0.8167  0.8167  0.7991  0.7991  0.5789  0.4778  0.4778
  0.4900  0.4900  0.5068  0.5068  0.3158  0.3158  0.2829  0.3482  0.3482  0.4367
  0.4367  0.4386  0.4386  0.3969  0.3969  0.3980  0.3980  0.3376  0.3376  0.0449
  0.0449  0.2442  0.1840  0.1840  0.0864  0.2035  0.2035  0.1277  0.2043  0.2043
  0.1229  0.0133

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2370.81914964
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       266.95314533
  PAW double counting   =      2860.43599250    -2838.26789726
  entropy T*S    EENTRO =         0.03960872
  eigenvalues    EBANDS =      -836.58055064
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -8.04361100 eV

  energy without entropy =       -8.08321972  energy(sigma->0) =       -8.05681391


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




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



 eigenvalue-minimisations  :  3344
 total energy-change (2. order) :-0.1869713E+02  (-0.4006675E+01)
 number of electron      87.0000038 magnetization      38.7856673
 augmentation part        1.8653820 magnetization       3.5993838

 Broyden mixing:
  rms(total) = 0.23184E+02    rms(broyden)= 0.23183E+02
  rms(prec ) = 0.23266E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4067
  1.0184  1.0184  0.8667  0.8190  0.8190  0.7860  0.7860  0.5788  0.4778  0.4778
  0.4913  0.4913  0.4940  0.4940  0.3158  0.3158  0.2829  0.3481  0.3481  0.4364
  0.4364  0.4397  0.4397  0.3968  0.3968  0.3969  0.3969  0.3818  0.3371  0.3371
  0.3227  0.3227  0.1840  0.1840  0.0449  0.0449  0.0864  0.2038  0.2038  0.1277
  0.1229  0.0068  0.0068

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2370.92247757
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       267.36813038
  PAW double counting   =      2858.51943288    -2836.94905100
  entropy T*S    EENTRO =         0.03074256
  eigenvalues    EBANDS =      -854.98275525
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -26.74073802 eV

  energy without entropy =      -26.77148059  energy(sigma->0) =      -26.75098554


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




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



 eigenvalue-minimisations  :  3040
 total energy-change (2. order) : 0.1264909E+02  (-0.4539294E+01)
 number of electron      87.0000039 magnetization      38.8646165
 augmentation part        1.8900656 magnetization       3.2608649

 Broyden mixing:
  rms(total) = 0.32740E+01    rms(broyden)= 0.32735E+01
  rms(prec ) = 0.37783E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3986
  1.0214  1.0214  0.8671  0.8194  0.8194  0.7806  0.7806  0.5788  0.4778  0.4778
  0.4924  0.4924  0.4909  0.4909  0.3158  0.3158  0.2829  0.3481  0.3481  0.4357
  0.4357  0.4415  0.4415  0.3956  0.3956  0.3981  0.3981  0.3906  0.3371  0.3371
  0.3381  0.3381  0.1840  0.1840  0.2039  0.2039  0.1277  0.0864  0.0449  0.0449
  0.1229  0.0171  0.0071  0.0071

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.06968492
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       269.07922896
  PAW double counting   =      2857.65250167    -2833.61553717
  entropy T*S    EENTRO =         0.02185707
  eigenvalues    EBANDS =      -846.35525466
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -14.09164906 eV

  energy without entropy =      -14.11350613  energy(sigma->0) =      -14.09893475


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




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



 eigenvalue-minimisations  :  3584
 total energy-change (2. order) :-0.1619611E+04  (-0.8373762E+03)
 number of electron      87.0000021 magnetization      38.8072542
 augmentation part        1.7742057 magnetization       4.1606161

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3664
  0.8736  0.7829  0.7829  0.8078  0.6013  0.6013  0.3974  0.6050  0.6050  0.4505
  0.4505  0.3053  0.3053  0.2482  0.4960  0.4960  0.3479  0.3479  0.4564  0.4564
  0.3780  0.3780  0.4448  0.4448  0.3916  0.3916  0.3190  0.3029  0.3029  0.1506
  0.1506  0.1265  0.1265  0.1383  0.0072  0.0072  0.0015  0.0437  0.0437  0.0883

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.16854334
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       270.79589203
  PAW double counting   =      2837.77031952    -2823.94187569
  entropy T*S    EENTRO =         0.06866636
  eigenvalues    EBANDS =     -2457.42200295
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -1633.70230410 eV

  energy without entropy =    -1633.77097046  energy(sigma->0) =    -1633.72519289


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




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



 eigenvalue-minimisations  :  2992
 total energy-change (2. order) : 0.1658950E+04  (-0.3174886E+03)
 number of electron      87.0000045 magnetization      39.6214039
 augmentation part        2.2818783 magnetization       3.9698245

 Broyden mixing:
  rms(total) = 0.33599E+01    rms(broyden)= 0.33437E+01
  rms(prec ) = 0.36431E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3803
  0.8753  0.7936  0.7936  0.6522  0.6522  0.6910  0.6910  0.4458  0.5074  0.5074
  0.2935  0.2935  0.2736  0.4957  0.4957  0.3331  0.3331  0.5516  0.4855  0.4855
  0.3755  0.3755  0.4756  0.4756  0.3980  0.3980  0.3549  0.3549  0.3636  0.3636
  0.2025  0.2025  0.0072  0.0072  0.0018  0.0433  0.1007  0.1007  0.1263  0.1263
  0.0882

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.78920453
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       280.11859070
  PAW double counting   =      2869.29387888    -2837.37484699
  entropy T*S    EENTRO =         0.00339943
  eigenvalues    EBANDS =      -825.19927907
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        25.24777840 eV

  energy without entropy =       25.24437896  energy(sigma->0) =       25.24664525


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




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



 eigenvalue-minimisations  :  2816
 total energy-change (2. order) : 0.4789134E+02  (-0.6289653E+02)
 number of electron      87.0000060 magnetization      40.0392976
 augmentation part        2.1619682 magnetization      -4.1789740

 Broyden mixing:
  rms(total) = 0.36319E+01    rms(broyden)= 0.36247E+01
  rms(prec ) = 0.39437E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3961
  0.8616  0.8939  0.8939  0.7070  0.7070  0.5021  0.6967  0.6967  0.7326  0.5192
  0.5192  0.2912  0.2912  0.2754  0.4921  0.4921  0.3442  0.3442  0.5044  0.5044
  0.3818  0.3818  0.4736  0.4736  0.4117  0.4117  0.3942  0.3942  0.3768  0.3768
  0.3124  0.1971  0.1971  0.0071  0.0071  0.0018  0.0435  0.0995  0.0995  0.1196
  0.1196  0.0882

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.81330908
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       312.23342289
  PAW double counting   =      3277.34704100    -3230.10203365
  entropy T*S    EENTRO =         0.01372551
  eigenvalues    EBANDS =      -824.73496397
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        73.13912267 eV

  energy without entropy =       73.12539717  energy(sigma->0) =       73.13454751


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




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



 eigenvalue-minimisations  :  3200
 total energy-change (2. order) :-0.5724564E+03  (-0.4680999E+03)
 number of electron      86.9999987 magnetization      39.9263099
 augmentation part        1.8174155 magnetization       0.8836877

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3872
  0.8641  0.8961  0.8961  0.7082  0.7082  0.4947  0.7037  0.7037  0.7319  0.5190
  0.5190  0.2905  0.2905  0.2737  0.4918  0.4918  0.3441  0.3441  0.5049  0.5049
  0.3815  0.3815  0.4734  0.4734  0.4189  0.4189  0.3910  0.3910  0.3751  0.3751
  0.3081  0.1963  0.1963  0.1215  0.1215  0.0981  0.0981  0.0435  0.0882  0.0071
  0.0071  0.0018  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2369.95550444
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       288.12163701
  PAW double counting   =      3320.91614193    -3294.58714890
  entropy T*S    EENTRO =         0.02236156
  eigenvalues    EBANDS =     -1354.02996017
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -499.31723304 eV

  energy without entropy =     -499.33959459  energy(sigma->0) =     -499.32468689


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




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



 eigenvalue-minimisations  :  3368
 total energy-change (2. order) :-0.1173093E+05  (-0.2828376E+04)
 number of electron      86.9999980 magnetization      39.7798622
 augmentation part        1.6165437 magnetization       2.2502114

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3786
  0.8642  0.8975  0.8975  0.7067  0.7067  0.4952  0.7051  0.7051  0.7325  0.5199
  0.5199  0.2903  0.2903  0.2751  0.4921  0.4921  0.3440  0.3440  0.5047  0.5047
  0.3813  0.3813  0.4732  0.4732  0.4192  0.4192  0.3904  0.3904  0.3750  0.3750
  0.3081  0.1966  0.1966  0.1218  0.1218  0.0971  0.0971  0.0435  0.0897  0.0071
  0.0071  0.0038  0.0018  0.0007

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2369.66001700
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       299.59018416
  PAW double counting   =      3351.91806105    -3327.01892629
  entropy T*S    EENTRO =         0.12626383
  eigenvalues    EBANDS =    -13095.39771914
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -12230.24691341 eV

  energy without entropy =   -12230.37317724  energy(sigma->0) =   -12230.28900136


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




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



 eigenvalue-minimisations  :  2728
 total energy-change (2. order) : 0.1232414E+05  (-0.4353683E+03)
 number of electron      87.0000042 magnetization      42.1952171
 augmentation part        2.2110303 magnetization      -2.1055836

 Broyden mixing:
  rms(total) = 0.49659E+01    rms(broyden)= 0.45722E+01
  rms(prec ) = 0.47870E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3541
  0.7033  0.7693  0.7693  0.5393  0.5393  0.4490  0.5145  0.5145  0.6683  0.5995
  0.5995  0.5699  0.5699  0.2495  0.2524  0.2524  0.4195  0.4195  0.2617  0.2617
  0.4069  0.4069  0.4174  0.4174  0.3933  0.3933  0.2456  0.2456  0.2691  0.2691
  0.2502  0.1350  0.1140  0.1140  0.1014  0.0468  0.0073  0.0073  0.0018  0.0008

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2379.41462237
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       335.52447310
  PAW double counting   =      2739.48260384    -2698.45299467
  entropy T*S    EENTRO =         0.00838871
  eigenvalues    EBANDS =      -813.45141416
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        93.89167443 eV

  energy without entropy =       93.88328572  energy(sigma->0) =       93.88887820


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




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



 eigenvalue-minimisations  :  3288
 total energy-change (2. order) :-0.3088206E+04  (-0.2462265E+04)
 number of electron      87.0000046 magnetization      41.3767862
 augmentation part        2.3150309 magnetization       8.8309399

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3656
  0.7035  0.8693  0.8693  0.5875  0.5875  0.4326  0.5480  0.5480  0.6484  0.6484
  0.2702  0.2702  0.2378  0.4170  0.4170  0.5898  0.5898  0.5276  0.2252  0.2252
  0.3949  0.3949  0.4429  0.4429  0.3245  0.3245  0.3963  0.3963  0.2775  0.2775
  0.2955  0.2955  0.1338  0.1160  0.1160  0.0838  0.0478  0.0074  0.0074  0.0018
  0.0008

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2359.76173582
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       290.82413229
  PAW double counting   =      3090.12099587    -3432.54779659
  entropy T*S    EENTRO =        -0.01396131
  eigenvalues    EBANDS =     -3493.13084558
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -2994.31397117 eV

  energy without entropy =    -2994.30000986  energy(sigma->0) =    -2994.30931740


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




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



 eigenvalue-minimisations  :  2856
 total energy-change (2. order) : 0.2988294E+04  (-0.2425451E+03)
 number of electron      87.0000029 magnetization      44.7443372
 augmentation part        2.1440199 magnetization       4.2298361

 Broyden mixing:
  rms(total) = 0.17430E+02    rms(broyden)= 0.17428E+02
  rms(prec ) = 0.17463E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3722
  0.8913  0.8913  0.6945  0.5668  0.5668  0.4672  0.7040  0.7040  0.5678  0.5678
  0.6382  0.2480  0.2623  0.2623  0.4254  0.4254  0.2504  0.2504  0.4735  0.4735
  0.4728  0.4502  0.4502  0.3948  0.3948  0.4064  0.4064  0.3301  0.3301  0.2990
  0.2990  0.2716  0.2716  0.1220  0.1220  0.1246  0.0878  0.0517  0.0074  0.0074
  0.0018  0.0008

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2365.96201637
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       286.48270041
  PAW double counting   =      3422.17890416    -3391.19670388
  entropy T*S    EENTRO =        -0.02485686
  eigenvalues    EBANDS =      -867.69365801
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        -6.02039058 eV

  energy without entropy =       -5.99553373  energy(sigma->0) =       -6.01210496


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




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



 eigenvalue-minimisations  :  3264
 total energy-change (2. order) :-0.6900365E+01  (-0.6014273E+02)
 number of electron      87.0000031 magnetization      44.1179582
 augmentation part        2.3644597 magnetization      14.6461578

 Broyden mixing:
  rms(total) = 0.11175E+02    rms(broyden)= 0.11169E+02
  rms(prec ) = 0.11466E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3674
  0.8736  0.8736  0.6945  0.5607  0.5607  0.7453  0.7453  0.4383  0.6198  0.6198
  0.2488  0.2596  0.2596  0.4233  0.4233  0.5509  0.5509  0.4641  0.4641  0.2364
  0.2364  0.3994  0.3994  0.4272  0.4272  0.4082  0.4082  0.3273  0.3273  0.2938
  0.2938  0.2736  0.2736  0.1489  0.1420  0.1420  0.0945  0.0945  0.0529  0.0074
  0.0074  0.0018  0.0007

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2371.19572765
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       321.99801166
  PAW double counting   =      3515.97787694    -3476.41924489
  entropy T*S    EENTRO =        -0.01770158
  eigenvalues    EBANDS =      -913.45920970
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -12.92075524 eV

  energy without entropy =      -12.90305366  energy(sigma->0) =      -12.91485471


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




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



 eigenvalue-minimisations  :  3376
 total energy-change (2. order) : 0.6522373E+02  (-0.8103107E+03)
 number of electron      87.0000025 magnetization      45.1227880
 augmentation part        1.6027445 magnetization       9.3695714

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3607
  0.8881  0.8881  0.6955  0.7407  0.7407  0.5568  0.5568  0.6146  0.6146  0.4421
  0.4253  0.4253  0.5518  0.5518  0.2591  0.2591  0.2510  0.4644  0.4644  0.3994
  0.3994  0.4277  0.4277  0.4087  0.4087  0.2372  0.2372  0.3253  0.3253  0.2916
  0.2916  0.2763  0.2763  0.1658  0.1658  0.1339  0.0594  0.0594  0.0742  0.0742
  0.0073  0.0073  0.0018  0.0008

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =      -600.23300183
  -Hartree energ DENC   =     -2375.53430447
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       325.79327619
  PAW double counting   =      3491.46471624    -3377.37256491
  entropy T*S    EENTRO =         0.02526897
  eigenvalues    EBANDS =      -922.26865770
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        52.30297429 eV

  energy without entropy =       52.27770532  energy(sigma->0) =       52.29455130


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




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


 -----------------------------------------------------------------------------
|                                                                             |
|     EEEEEEE  RRRRRR   RRRRRR   OOOOOOO  RRRRRR      ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     EEEEE    RRRRRR   RRRRRR   O     O  RRRRRR       #       #       #      |
|     E        R   R    R   R    O     O  R   R                               |
|     E        R    R   R    R   O     O  R    R      ###     ###     ###     |
|     EEEEEEE  R     R  R     R  OOOOOOO  R     R     ###     ###     ###     |
|                                                                             |
|     Error EDDDAV: Call to ZHEGV failed. Returncode = 7 1 8                  |
|                                                                             |
|       ---->  I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <----       |
|                                                                             |
 -----------------------------------------------------------------------------