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.13  09:46:27
 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.693  0.622  0.446-
   2  0.712  0.661  0.645-
   3  0.593  0.388  0.423-
   4  0.593  0.477  0.281-
   5  0.420  0.534  0.413-
   6  0.492  0.571  0.676-
   7  0.485  0.025  0.876-
   8  0.270  0.349  0.780-
   9  0.396  0.392  0.247-
  10  0.249  0.311  0.326-
  11  0.297  0.386  0.904-
  12  0.533  0.359  0.328-  17 0.45
  13  0.472  0.976  0.556-
  14  0.855  0.597  0.276-
  15  0.345  0.686  0.920-
  16  0.967  0.267  0.765-
  17  0.526  0.388  0.331-  12 0.45
  18  0.040  0.346  0.292-
  19  0.040  0.385  0.613-
  20  0.547  0.115  0.408-
  21  0.478  0.412  0.156-
  22  0.571  0.762  0.252-
  23  0.144  0.136  0.451-
  24  0.269  0.483  0.359-
  25  0.466  0.561  0.577-  32 1.51
  26  0.480  0.582  0.283-
  27  0.163  0.348  0.861-
  28  0.469  0.251  0.496-  33 1.48
  29  0.395  0.113  0.353-
  30  0.334  0.549  0.189-
  31  0.964  0.974  0.763-
  32  0.534  0.539  0.540-  25 1.51
  33  0.420  0.324  0.507-  28 1.48
  34  0.832  0.127  0.228-
  35  0.687  0.222  0.609-
  36  0.661  0.407  0.937-
 
  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.692990530  0.622213090  0.446246620
     0.711976830  0.660814770  0.644620210
     0.593300380  0.387615420  0.422811770
     0.592702040  0.476877020  0.281490330
     0.419807370  0.534108500  0.413275690
     0.491971920  0.570936370  0.675705070
     0.484636190  0.025287050  0.875678470
     0.269628660  0.349119780  0.779705160
     0.395819100  0.392298890  0.246806360
     0.248743850  0.310891980  0.325554430
     0.296576220  0.386059690  0.903573650
     0.532821250  0.359255020  0.328265390
     0.471640270  0.975777780  0.556031220
     0.855147850  0.596612280  0.276260910
     0.345488300  0.686185780  0.919983090
     0.966690700  0.267412380  0.765200230
     0.526346300  0.387963290  0.330707990
     0.039732170  0.346332530  0.292197070
     0.039637810  0.385462160  0.613034670
     0.547039750  0.115008010  0.408428750
     0.477555750  0.411783910  0.155841150
     0.570976570  0.761720340  0.251877150
     0.144481890  0.135860040  0.450888410
     0.269357420  0.483196720  0.358941440
     0.465818180  0.561188960  0.577146810
     0.480238070  0.581983980  0.282817000
     0.163160440  0.347950090  0.861429250
     0.469108740  0.250575490  0.496193580
     0.394705590  0.112550720  0.353127760
     0.334466520  0.548769200  0.189107590
     0.963903540  0.974431650  0.762610870
     0.533723100  0.538709360  0.539539040
     0.420488520  0.323986320  0.507082510
     0.831993570  0.127005800  0.228241770
     0.686705570  0.222067820  0.608590310
     0.661190130  0.407261870  0.937226720

  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.69299053  0.62221309  0.44624662
   0.71197683  0.66081477  0.64462021
   0.59330038  0.38761542  0.42281177
   0.59270204  0.47687702  0.28149033
   0.41980737  0.53410850  0.41327569
   0.49197192  0.57093637  0.67570507
   0.48463619  0.02528705  0.87567847
   0.26962866  0.34911978  0.77970516
   0.39581910  0.39229889  0.24680636
   0.24874385  0.31089198  0.32555443
   0.29657622  0.38605969  0.90357365
   0.53282125  0.35925502  0.32826539
   0.47164027  0.97577778  0.55603122
   0.85514785  0.59661228  0.27626091
   0.34548830  0.68618578  0.91998309
   0.96669070  0.26741238  0.76520023
   0.52634630  0.38796329  0.33070799
   0.03973217  0.34633253  0.29219707
   0.03963781  0.38546216  0.61303467
   0.54703975  0.11500801  0.40842875
   0.47755575  0.41178391  0.15584115
   0.57097657  0.76172034  0.25187715
   0.14448189  0.13586004  0.45088841
   0.26935742  0.48319672  0.35894144
   0.46581818  0.56118896  0.57714681
   0.48023807  0.58198398  0.28281700
   0.16316044  0.34795009  0.86142925
   0.46910874  0.25057549  0.49619358
   0.39470559  0.11255072  0.35312776
   0.33446652  0.54876920  0.18910759
   0.96390354  0.97443165  0.76261087
   0.53372310  0.53870936  0.53953904
   0.42048852  0.32398632  0.50708251
   0.83199357  0.12700580  0.22824177
   0.68670557  0.22206782  0.60859031
   0.66119013  0.40726187  0.93722672
 
 position of ions in cartesian coordinates  (Angst):
  13.85981060  9.33319635  6.69369930
  14.23953660  9.91222155  9.66930315
  11.86600760  5.81423130  6.34217655
  11.85404080  7.15315530  4.22235495
   8.39614740  8.01162750  6.19913535
   9.83943840  8.56404555 10.13557605
   9.69272380  0.37930575 13.13517705
   5.39257320  5.23679670 11.69557740
   7.91638200  5.88448335  3.70209540
   4.97487700  4.66337970  4.88331645
   5.93152440  5.79089535 13.55360475
  10.65642500  5.38882530  4.92398085
   9.43280540 14.63666670  8.34046830
  17.10295700  8.94918420  4.14391365
   6.90976600 10.29278670 13.79974635
  19.33381400  4.01118570 11.47800345
  10.52692600  5.81944935  4.96061985
   0.79464340  5.19498795  4.38295605
   0.79275620  5.78193240  9.19552005
  10.94079500  1.72512015  6.12643125
   9.55111500  6.17675865  2.33761725
  11.41953140 11.42580510  3.77815725
   2.88963780  2.03790060  6.76332615
   5.38714840  7.24795080  5.38412160
   9.31636360  8.41783440  8.65720215
   9.60476140  8.72975970  4.24225500
   3.26320880  5.21925135 12.92143875
   9.38217480  3.75863235  7.44290370
   7.89411180  1.68826080  5.29691640
   6.68933040  8.23153800  2.83661385
  19.27807080 14.61647475 11.43916305
  10.67446200  8.08064040  8.09308560
   8.40977040  4.85979480  7.60623765
  16.63987140  1.90508700  3.42362655
  13.73411140  3.33101730  9.12885465
  13.22380260  6.10892805 14.05840080
 


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


 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  1277132. kBytes
=======================================================================

   base      :      30000. kBytes
   nonlr-proj:       4073. 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         2538
 Maximum index for augmentation-charges         2274 (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  :  2816
 total energy-change (2. order) : 0.5274549E+03  (-0.1851775E+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  =       239.48782979
  -Hartree energ DENC   =     -3114.85518880
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       230.50610591
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.00539561
  eigenvalues    EBANDS =      -360.84784259
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       527.45490252 eV

  energy without entropy =      527.44950691  energy(sigma->0) =      527.45310398


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




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



 eigenvalue-minimisations  :  4184
 total energy-change (2. order) :-0.3944313E+03  (-0.3716295E+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  =       239.48782979
  -Hartree energ DENC   =     -3114.85518880
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       230.50610591
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.00333890
  eigenvalues    EBANDS =      -755.27703952
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       133.02364888 eV

  energy without entropy =      133.02030998  energy(sigma->0) =      133.02253591


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




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



 eigenvalue-minimisations  :  3720
 total energy-change (2. order) :-0.1848826E+03  (-0.1653726E+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  =       239.48782979
  -Hartree energ DENC   =     -3114.85518880
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       230.50610591
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =        -0.03603918
  eigenvalues    EBANDS =      -940.12021274
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -51.85890243 eV

  energy without entropy =      -51.82286325  energy(sigma->0) =      -51.84688937


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




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



 eigenvalue-minimisations  :  3208
 total energy-change (2. order) :-0.1893965E+02  (-0.1813334E+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  =       239.48782979
  -Hartree energ DENC   =     -3114.85518880
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       230.50610591
  PAW double counting   =      2564.72054280    -2541.99018201
  entropy T*S    EENTRO =         0.01142478
  eigenvalues    EBANDS =      -959.10732264
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -70.79854836 eV

  energy without entropy =      -70.80997314  energy(sigma->0) =      -70.80235662


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




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



 eigenvalue-minimisations  :  3192
 total energy-change (2. order) :-0.5305993E+00  (-0.5206668E+00)
 number of electron      86.9999963 magnetization      34.1299503
 augmentation part        3.3551821 magnetization      11.7395661

 Broyden mixing:
  rms(total) = 0.35457E+01    rms(broyden)= 0.35442E+01
  rms(prec ) = 0.43346E+01
  weight for this iteration     100.00

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

  energy without entropy =      -71.34737534  energy(sigma->0) =      -71.33522355


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




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



 eigenvalue-minimisations  :  3016
 total energy-change (2. order) :-0.8011371E+02  (-0.2733057E+02)
 number of electron      86.9999950 magnetization      33.5030459
 augmentation part        0.3359919 magnetization      14.0809823

 Broyden mixing:
  rms(total) = 0.72678E+01    rms(broyden)= 0.72643E+01
  rms(prec ) = 0.92209E+01
  weight for this iteration     100.00

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3297.87875673
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       280.91559141
  PAW double counting   =      3104.21646812    -3081.63951709
  entropy T*S    EENTRO =        -0.03937911
  eigenvalues    EBANDS =      -906.93333503
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -151.44285684 eV

  energy without entropy =     -151.40347773  energy(sigma->0) =     -151.42973047


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




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



 eigenvalue-minimisations  :  3048
 total energy-change (2. order) : 0.1280022E+03  (-0.5973677E+01)
 number of electron      86.9999970 magnetization      33.0481443
 augmentation part        2.1638792 magnetization      10.9449877

 Broyden mixing:
  rms(total) = 0.31993E+01    rms(broyden)= 0.31945E+01
  rms(prec ) = 0.42099E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3543
  0.4937  0.2149

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3217.28809015
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       274.38144854
  PAW double counting   =      3030.20787637    -3007.36496166
  entropy T*S    EENTRO =        -0.04123171
  eigenvalues    EBANDS =      -853.25173218
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -23.44061920 eV

  energy without entropy =      -23.39938749  energy(sigma->0) =      -23.42687529


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




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



 eigenvalue-minimisations  :  2752
 total energy-change (2. order) : 0.8922248E+01  (-0.4664797E+01)
 number of electron      86.9999966 magnetization      33.1110212
 augmentation part        1.4283588 magnetization       2.9833134

 Broyden mixing:
  rms(total) = 0.35725E+01    rms(broyden)= 0.35703E+01
  rms(prec ) = 0.41903E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2605
  0.3279  0.3279  0.1257

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3208.45975748
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       271.99585866
  PAW double counting   =      3082.16824423    -3059.99866572
  entropy T*S    EENTRO =        -0.32373949
  eigenvalues    EBANDS =      -849.81638275
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -14.51837096 eV

  energy without entropy =      -14.19463146  energy(sigma->0) =      -14.41045779


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




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



 eigenvalue-minimisations  :  3096
 total energy-change (2. order) :-0.3697823E+01  (-0.5326505E+00)
 number of electron      86.9999938 magnetization      33.3453608
 augmentation part        1.5020817 magnetization      -4.1136897

 Broyden mixing:
  rms(total) = 0.34625E+01    rms(broyden)= 0.34613E+01
  rms(prec ) = 0.37165E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2180
  0.3491  0.2113  0.1559  0.1559

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3208.93106717
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       267.75569599
  PAW double counting   =      3092.15926968    -3069.86369289
  entropy T*S    EENTRO =        -0.10916911
  eigenvalues    EBANDS =      -849.14330246
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -18.21619437 eV

  energy without entropy =      -18.10702526  energy(sigma->0) =      -18.17980467


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




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



 eigenvalue-minimisations  :  3072
 total energy-change (2. order) :-0.1020543E+02  (-0.1009734E+01)
 number of electron      87.0000019 magnetization      33.4619411
 augmentation part        1.6415329 magnetization       3.0046600

 Broyden mixing:
  rms(total) = 0.35006E+01    rms(broyden)= 0.34989E+01
  rms(prec ) = 0.42243E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1987
  0.3492  0.1987  0.1987  0.1968  0.0498

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3205.96817192
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       262.84055755
  PAW double counting   =      3089.95447029    -3067.58279506
  entropy T*S    EENTRO =        -0.15820471
  eigenvalues    EBANDS =      -857.42355579
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -28.42162805 eV

  energy without entropy =      -28.26342334  energy(sigma->0) =      -28.36889315


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




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



 eigenvalue-minimisations  :  3312
 total energy-change (2. order) :-0.3245303E+01  (-0.1428522E+01)
 number of electron      86.9999920 magnetization      33.5109692
 augmentation part        1.8037659 magnetization       0.9848754

 Broyden mixing:
  rms(total) = 0.34967E+01    rms(broyden)= 0.34953E+01
  rms(prec ) = 0.40440E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1711
  0.3481  0.1976  0.1976  0.1955  0.0506  0.0373

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3204.23391597
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       260.62852009
  PAW double counting   =      3088.75232471    -3066.33404407
  entropy T*S    EENTRO =         0.06035730
  eigenvalues    EBANDS =      -860.45624443
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -31.66693078 eV

  energy without entropy =      -31.72728808  energy(sigma->0) =      -31.68704988


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




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



 eigenvalue-minimisations  :  3360
 total energy-change (2. order) : 0.1721415E+00  (-0.8851451E-01)
 number of electron      86.9999934 magnetization      33.6675262
 augmentation part        1.6446391 magnetization      -1.4048415

 Broyden mixing:
  rms(total) = 0.31759E+01    rms(broyden)= 0.31755E+01
  rms(prec ) = 0.35656E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1967
  0.2646  0.2646  0.3533  0.2367  0.1053  0.1053  0.0473

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3204.50619309
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       260.29154109
  PAW double counting   =      3089.41669891    -3066.99348395
  entropy T*S    EENTRO =        -0.11466935
  eigenvalues    EBANDS =      -859.50475450
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -31.49478930 eV

  energy without entropy =      -31.38011995  energy(sigma->0) =      -31.45656618


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




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



 eigenvalue-minimisations  :  2864
 total energy-change (2. order) :-0.4474921E+01  (-0.5115069E+00)
 number of electron      86.9999951 magnetization      33.7832532
 augmentation part        1.8671352 magnetization       1.5664768

 Broyden mixing:
  rms(total) = 0.29928E+01    rms(broyden)= 0.29917E+01
  rms(prec ) = 0.35094E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1955
  0.3137  0.3137  0.3602  0.2241  0.1207  0.1207  0.0662  0.0451

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3201.21151301
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       258.26932583
  PAW double counting   =      3087.81586915    -3065.39068461
  entropy T*S    EENTRO =        -0.09510591
  eigenvalues    EBANDS =      -865.27367312
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -35.96971007 eV

  energy without entropy =      -35.87460415  energy(sigma->0) =      -35.93800810


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




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



 eigenvalue-minimisations  :  2960
 total energy-change (2. order) :-0.1709282E+01  (-0.1280834E+00)
 number of electron      86.9999938 magnetization      33.7694728
 augmentation part        1.6905771 magnetization       0.4648648

 Broyden mixing:
  rms(total) = 0.32554E+01    rms(broyden)= 0.32551E+01
  rms(prec ) = 0.37893E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1776
  0.3077  0.3077  0.3631  0.2208  0.1209  0.1209  0.0729  0.0459  0.0384

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3201.45186708
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.91643474
  PAW double counting   =      3090.39387666    -3067.95231874
  entropy T*S    EENTRO =         0.02524340
  eigenvalues    EBANDS =      -865.52643254
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.67899197 eV

  energy without entropy =      -37.70423537  energy(sigma->0) =      -37.68740643


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




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



 eigenvalue-minimisations  :  3040
 total energy-change (2. order) : 0.9901889E-01  (-0.4397171E-02)
 number of electron      86.9999938 magnetization      33.7791830
 augmentation part        1.6715912 magnetization       0.4048140

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.1805
  0.3629  0.2937  0.2937  0.2152  0.1376  0.1376  0.1169  0.1169  0.0841  0.0465

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3201.63128463
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.00849416
  PAW double counting   =      3090.57869713    -3068.14100475
  entropy T*S    EENTRO =         0.02902593
  eigenvalues    EBANDS =      -865.33997250
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.57997307 eV

  energy without entropy =      -37.60899900  energy(sigma->0) =      -37.58964838


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




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



 eigenvalue-minimisations  :  2968
 total energy-change (2. order) : 0.1611578E-01  (-0.1001274E-01)
 number of electron      86.9999939 magnetization      33.7191281
 augmentation part        1.6676103 magnetization       0.2929507

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2057
  0.2911  0.2911  0.3626  0.2887  0.2887  0.2358  0.1354  0.1354  0.0934  0.0934
  0.0467

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3201.59971781
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.04628961
  PAW double counting   =      3091.43918407    -3068.99541113
  entropy T*S    EENTRO =         0.02246586
  eigenvalues    EBANDS =      -865.39273950
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.56385730 eV

  energy without entropy =      -37.58632316  energy(sigma->0) =      -37.57134592


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




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



 eigenvalue-minimisations  :  2848
 total energy-change (2. order) : 0.5829422E+00  (-0.3111070E-01)
 number of electron      86.9999935 magnetization      33.8441657
 augmentation part        1.6956418 magnetization       0.5535102

 Broyden mixing:
  rms(total) = 0.33225E+01    rms(broyden)= 0.33224E+01
  rms(prec ) = 0.38511E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2377
  0.4289  0.4289  0.3007  0.3666  0.2738  0.2738  0.2473  0.1437  0.1437  0.0991
  0.0991  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3201.55324694
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.52937652
  PAW double counting   =      3090.98266082    -3068.55919597
  entropy T*S    EENTRO =         0.03159423
  eigenvalues    EBANDS =      -865.32817537
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -36.98091511 eV

  energy without entropy =      -37.01250934  energy(sigma->0) =      -36.99144652


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




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



 eigenvalue-minimisations  :  2760
 total energy-change (2. order) :-0.7507548E+00  (-0.9093520E-01)
 number of electron      86.9999949 magnetization      33.8642729
 augmentation part        1.6834181 magnetization       0.5609600

 Broyden mixing:
  rms(total) = 0.31592E+01    rms(broyden)= 0.31587E+01
  rms(prec ) = 0.36922E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2382
  0.4563  0.4563  0.3222  0.3716  0.2964  0.2964  0.2300  0.1478  0.1478  0.1207
  0.1020  0.1020  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3199.95930580
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.78338525
  PAW double counting   =      3088.01384493    -3065.56903794
  entropy T*S    EENTRO =        -0.00181191
  eigenvalues    EBANDS =      -866.91481607
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.73166995 eV

  energy without entropy =      -37.72985804  energy(sigma->0) =      -37.73106598


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




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



 eigenvalue-minimisations  :  3088
 total energy-change (2. order) :-0.1904387E+00  (-0.6196458E-02)
 number of electron      86.9999949 magnetization      33.9005920
 augmentation part        1.6743304 magnetization       0.5291802

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2702
  0.5102  0.5102  0.3547  0.3547  0.3846  0.3296  0.3296  0.2253  0.2253  0.1552
  0.1552  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3200.01083757
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.62561792
  PAW double counting   =      3088.32148972    -3065.87519462
  entropy T*S    EENTRO =         0.01812006
  eigenvalues    EBANDS =      -866.91737570
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.92210860 eV

  energy without entropy =      -37.94022867  energy(sigma->0) =      -37.92814863


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




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



 eigenvalue-minimisations  :  3104
 total energy-change (2. order) :-0.1904252E+00  (-0.2503691E-01)
 number of electron      86.9999951 magnetization      33.9371263
 augmentation part        1.6513182 magnetization       0.4982574

 Broyden mixing:
  rms(total) = 0.31612E+01    rms(broyden)= 0.31611E+01
  rms(prec ) = 0.36847E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2679
  0.5366  0.5366  0.3902  0.3902  0.3891  0.3418  0.3418  0.2231  0.2231  0.1566
  0.1566  0.1006  0.1006  0.0466  0.0848

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3199.66289327
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.48761143
  PAW double counting   =      3085.95746196    -3063.49788921
  entropy T*S    EENTRO =         0.05556815
  eigenvalues    EBANDS =      -867.36846449
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -38.11253385 eV

  energy without entropy =      -38.16810200  energy(sigma->0) =      -38.13105657


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




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



 eigenvalue-minimisations  :  3200
 total energy-change (2. order) :-0.3473771E+00  (-0.1388265E-01)
 number of electron      86.9999950 magnetization      33.9493108
 augmentation part        1.6403353 magnetization       0.3678350

 Broyden mixing:
  rms(total) = 0.31624E+01    rms(broyden)= 0.31623E+01
  rms(prec ) = 0.36900E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2830
  0.5688  0.5688  0.3999  0.3999  0.3510  0.3510  0.3714  0.2604  0.2604  0.2209
  0.2209  0.1537  0.1537  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3199.31444521
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.28899270
  PAW double counting   =      3085.44414652    -3062.97984084
  entropy T*S    EENTRO =         0.05074169
  eigenvalues    EBANDS =      -867.86557742
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -38.45991097 eV

  energy without entropy =      -38.51065266  energy(sigma->0) =      -38.47682487


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




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



 eigenvalue-minimisations  :  3072
 total energy-change (2. order) :-0.1240597E+00  (-0.1892616E-01)
 number of electron      86.9999955 magnetization      34.1189772
 augmentation part        1.6550132 magnetization       0.2240593

 Broyden mixing:
  rms(total) = 0.30112E+01    rms(broyden)= 0.30111E+01
  rms(prec ) = 0.34995E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3353
  0.7338  0.7338  0.4129  0.4129  0.4898  0.4898  0.3885  0.3885  0.4043  0.2582
  0.2582  0.1542  0.1542  0.1743  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3200.07113898
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.10180277
  PAW double counting   =      3089.05134025    -3066.59069272
  entropy T*S    EENTRO =        -0.00830449
  eigenvalues    EBANDS =      -866.98304909
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -38.58397067 eV

  energy without entropy =      -38.57566618  energy(sigma->0) =      -38.58120251


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




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



 eigenvalue-minimisations  :  2856
 total energy-change (2. order) :-0.2953117E+01  (-0.1340943E+00)
 number of electron      86.9999953 magnetization      34.5414776
 augmentation part        1.5042107 magnetization      -0.7983714

 Broyden mixing:
  rms(total) = 0.31831E+01    rms(broyden)= 0.31821E+01
  rms(prec ) = 0.36893E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3423
  0.8058  0.8058  0.5112  0.5112  0.4130  0.4130  0.3816  0.3816  0.3290  0.3290
  0.2732  0.2732  0.1541  0.1541  0.1776  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3200.49493258
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       254.49122177
  PAW double counting   =      3094.95615685    -3072.48746851
  entropy T*S    EENTRO =        -0.01564805
  eigenvalues    EBANDS =      -867.90248835
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -41.53708728 eV

  energy without entropy =      -41.52143923  energy(sigma->0) =      -41.53187126


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




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



 eigenvalue-minimisations  :  2704
 total energy-change (2. order) :-0.1301737E+01  (-0.2380225E+00)
 number of electron      86.9999969 magnetization      35.0154799
 augmentation part        1.8144793 magnetization       0.9720777

 Broyden mixing:
  rms(total) = 0.26134E+01    rms(broyden)= 0.26121E+01
  rms(prec ) = 0.30466E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3662
  0.8772  0.8772  0.6049  0.6049  0.4153  0.4153  0.3994  0.3994  0.3944  0.3477
  0.3477  0.2707  0.2707  0.1542  0.1542  0.1775  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3195.75763722
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       253.33266123
  PAW double counting   =      3090.19576435    -3067.61720676
  entropy T*S    EENTRO =        -0.05008539
  eigenvalues    EBANDS =      -872.85839248
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -42.83882466 eV

  energy without entropy =      -42.78873928  energy(sigma->0) =      -42.82212953


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




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



 eigenvalue-minimisations  :  3008
 total energy-change (2. order) :-0.1228507E+01  (-0.5293689E+00)
 number of electron      86.9999949 magnetization      35.0069115
 augmentation part        1.9292678 magnetization      -0.0594905

 Broyden mixing:
  rms(total) = 0.21171E+01    rms(broyden)= 0.21147E+01
  rms(prec ) = 0.23733E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3573
  0.8719  0.8719  0.5801  0.5801  0.4166  0.4166  0.3753  0.3753  0.4114  0.3756
  0.3756  0.2633  0.2633  0.2334  0.1541  0.1541  0.1786  0.1005  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3191.50230427
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.19773920
  PAW double counting   =      3088.02533350    -3065.36350362
  entropy T*S    EENTRO =        -0.02616610
  eigenvalues    EBANDS =      -877.31450180
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -44.06733150 eV

  energy without entropy =      -44.04116539  energy(sigma->0) =      -44.05860946


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




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



 eigenvalue-minimisations  :  2976
 total energy-change (2. order) :-0.2008002E+01  (-0.2113088E+00)
 number of electron      86.9999949 magnetization      35.2023908
 augmentation part        1.8348989 magnetization      -0.7775539

 Broyden mixing:
  rms(total) = 0.20713E+01    rms(broyden)= 0.20707E+01
  rms(prec ) = 0.23532E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3799
  0.8732  0.8732  0.6880  0.6880  0.4184  0.4184  0.4358  0.4358  0.4493  0.4493
  0.3662  0.3662  0.2806  0.2806  0.2203  0.1541  0.1541  0.1790  0.1005  0.1005
  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3194.79815087
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       251.23749095
  PAW double counting   =      3109.34309475    -3086.74079398
  entropy T*S    EENTRO =        -0.07470862
  eigenvalues    EBANDS =      -874.95833776
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -46.07533393 eV

  energy without entropy =      -46.00062531  energy(sigma->0) =      -46.05043106


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




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



 eigenvalue-minimisations  :  2976
 total energy-change (2. order) :-0.2330091E+01  (-0.2768314E+00)
 number of electron      86.9999956 magnetization      35.6169332
 augmentation part        1.6908425 magnetization       0.4471139

 Broyden mixing:
  rms(total) = 0.24362E+01    rms(broyden)= 0.24354E+01
  rms(prec ) = 0.28557E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3946
  0.8649  0.8649  0.8173  0.8173  0.4189  0.4189  0.4650  0.4650  0.4402  0.4402
  0.4071  0.3683  0.3683  0.2725  0.2725  0.2445  0.1541  0.1541  0.1786  0.1005
  0.1005  0.0466

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3196.95055787
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       250.60677238
  PAW double counting   =      3125.48111134    -3102.84240683
  entropy T*S    EENTRO =        -0.09657135
  eigenvalues    EBANDS =      -874.51984453
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -48.40542525 eV

  energy without entropy =      -48.30885391  energy(sigma->0) =      -48.37323480


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




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



 eigenvalue-minimisations  :  3008
 total energy-change (2. order) : 0.1311182E+01  (-0.8925316E+00)
 number of electron      86.9999959 magnetization      36.2823389
 augmentation part        1.8532894 magnetization      -0.7903627

 Broyden mixing:
  rms(total) = 0.20047E+01    rms(broyden)= 0.20032E+01
  rms(prec ) = 0.22280E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4010
  0.9007  0.9007  0.8357  0.8357  0.5777  0.5777  0.4191  0.4191  0.4088  0.4088
  0.3780  0.3780  0.3191  0.3191  0.2724  0.2724  0.2658  0.1541  0.1541  0.0466
  0.1005  0.1005  0.1785

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3195.12966006
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       250.66619975
  PAW double counting   =      3136.08265630    -3113.37166434
  entropy T*S    EENTRO =        -0.03871616
  eigenvalues    EBANDS =      -875.21913044
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -47.09424335 eV

  energy without entropy =      -47.05552719  energy(sigma->0) =      -47.08133796


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




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



 eigenvalue-minimisations  :  2896
 total energy-change (2. order) : 0.1563843E+01  (-0.5642910E+00)
 number of electron      86.9999958 magnetization      37.2314594
 augmentation part        1.8362954 magnetization      -0.6822366

 Broyden mixing:
  rms(total) = 0.17477E+01    rms(broyden)= 0.17474E+01
  rms(prec ) = 0.19520E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4201
  0.8490  0.8490  0.8626  0.8626  0.7206  0.7206  0.4193  0.4193  0.4933  0.4933
  0.4101  0.4101  0.3963  0.3317  0.3317  0.2722  0.2722  0.2342  0.1541  0.1541
  0.0466  0.1005  0.1005  0.1786

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3194.66562511
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       252.27532105
  PAW double counting   =      3154.93276008    -3132.11455050
  entropy T*S    EENTRO =        -0.06411537
  eigenvalues    EBANDS =      -875.81026177
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -45.53040001 eV

  energy without entropy =      -45.46628464  energy(sigma->0) =      -45.50902822


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




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



 eigenvalue-minimisations  :  3024
 total energy-change (2. order) :-0.2870020E+00  (-0.3582353E+01)
 number of electron      86.9999967 magnetization      37.2520824
 augmentation part        2.0667114 magnetization       3.5423993

 Broyden mixing:
  rms(total) = 0.76711E+01    rms(broyden)= 0.76692E+01
  rms(prec ) = 0.78036E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4045
  0.8498  0.8498  0.8638  0.8638  0.7217  0.7217  0.4193  0.4193  0.4983  0.4983
  0.4105  0.4105  0.3997  0.3323  0.3323  0.2723  0.2723  0.2343  0.1541  0.1541
  0.1786  0.1005  0.1005  0.0466  0.0095

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3193.58389563
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       258.34430703
  PAW double counting   =      3192.73806892    -3170.53762877
  entropy T*S    EENTRO =        -0.06827680
  eigenvalues    EBANDS =      -882.62604840
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -45.81740206 eV

  energy without entropy =      -45.74912526  energy(sigma->0) =      -45.79464313


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




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



 eigenvalue-minimisations  :  3032
 total energy-change (2. order) : 0.1582430E+02  (-0.1493674E+03)
 number of electron      86.9999965 magnetization      37.3311816
 augmentation part        2.3079690 magnetization       0.4889197

 Broyden mixing:
  rms(total) = 0.16024E+01    rms(broyden)= 0.15964E+01
  rms(prec ) = 0.17322E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3939
  0.8540  0.8540  0.8763  0.8763  0.7283  0.7283  0.4193  0.4193  0.4859  0.4859
  0.4105  0.4105  0.4168  0.3216  0.3216  0.2719  0.2719  0.2315  0.1541  0.1541
  0.1786  0.1005  0.1005  0.0466  0.1137  0.0091

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3193.39275823
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       259.62301206
  PAW double counting   =      3185.69910260    -3158.05914721
  entropy T*S    EENTRO =        -0.07798263
  eigenvalues    EBANDS =      -873.70140138
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -29.99310321 eV

  energy without entropy =      -29.91512057  energy(sigma->0) =      -29.96710900


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




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



 eigenvalue-minimisations  :  3616
 total energy-change (2. order) :-0.6121261E+04  (-0.5454122E+04)
 number of electron      86.9999967 magnetization      37.2459561
 augmentation part        1.0300199 magnetization       1.3314286

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3795
  0.8535  0.8535  0.8744  0.8744  0.7291  0.7291  0.4193  0.4193  0.4859  0.4859
  0.4106  0.4106  0.4146  0.3219  0.3219  0.2718  0.2718  0.2311  0.1541  0.1541
  0.1787  0.1005  0.1005  0.0466  0.1257  0.0084  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3194.91815053
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       259.29142738
  PAW double counting   =      3197.45972517    -3169.03044810
  entropy T*S    EENTRO =        -0.03781334
  eigenvalues    EBANDS =     -6993.93517708
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -6151.25436491 eV

  energy without entropy =    -6151.21655157  energy(sigma->0) =    -6151.24176047


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




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



 eigenvalue-minimisations  :  3144
 total energy-change (2. order) : 0.6159796E+04  (-0.4143685E+03)
 number of electron      86.9999991 magnetization      38.0547007
 augmentation part        2.3053050 magnetization       1.6317056

 Broyden mixing:
  rms(total) = 0.27267E+01    rms(broyden)= 0.21953E+01
  rms(prec ) = 0.23257E+01
  weight for this iteration     100.00
  WARNING: grid for Broyden might be to small

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3713
  0.8548  0.8548  0.8933  0.8933  0.7196  0.7196  0.4193  0.4193  0.4834  0.4834
  0.4118  0.4118  0.4093  0.3267  0.3267  0.2715  0.2715  0.2298  0.1541  0.1541
  0.1005  0.1005  0.0466  0.1781  0.1813  0.0740  0.0020  0.0041

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3187.94047902
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       291.92234972
  PAW double counting   =      3194.19573007    -3168.89571255
  entropy T*S    EENTRO =        -0.04332168
  eigenvalues    EBANDS =      -870.61349588
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =         8.54114225 eV

  energy without entropy =        8.58446393  energy(sigma->0) =        8.55558281


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




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



 eigenvalue-minimisations  :  2928
 total energy-change (2. order) :-0.2121417E+02  (-0.7682315E+01)
 number of electron      86.9999993 magnetization      38.3659972
 augmentation part        2.3407372 magnetization      -1.1082246

 Broyden mixing:
  rms(total) = 0.24926E+01    rms(broyden)= 0.24697E+01
  rms(prec ) = 0.26924E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3729
  0.8536  0.8536  0.8856  0.8856  0.7140  0.7140  0.4193  0.4193  0.4785  0.4785
  0.4144  0.4144  0.4185  0.2186  0.3251  0.3251  0.2721  0.2721  0.2390  0.2390
  0.2314  0.1541  0.1541  0.1786  0.1005  0.1005  0.0466  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3177.48986720
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.83543368
  PAW double counting   =      3256.56202759    -3230.10755540
  entropy T*S    EENTRO =        -0.01050504
  eigenvalues    EBANDS =      -888.37863010
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -12.67302486 eV

  energy without entropy =      -12.66251983  energy(sigma->0) =      -12.66952318


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




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



 eigenvalue-minimisations  :  3096
 total energy-change (2. order) :-0.4198240E+01  (-0.8471284E+00)
 number of electron      86.9999987 magnetization      38.8157622
 augmentation part        2.3279420 magnetization      -0.7155110

 Broyden mixing:
  rms(total) = 0.22292E+01    rms(broyden)= 0.22274E+01
  rms(prec ) = 0.24352E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3729
  0.8507  0.8507  0.8827  0.8827  0.7078  0.7078  0.4193  0.4193  0.3079  0.4835
  0.4835  0.4149  0.4149  0.4212  0.3039  0.3039  0.3062  0.3062  0.2725  0.2725
  0.1541  0.1541  0.0466  0.1005  0.1005  0.2333  0.1786  0.2005  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3179.41256310
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       276.10151452
  PAW double counting   =      3260.64456126    -3236.22992137
  entropy T*S    EENTRO =        -0.04296905
  eigenvalues    EBANDS =      -887.84795869
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -16.87126484 eV

  energy without entropy =      -16.82829579  energy(sigma->0) =      -16.85694182


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




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



 eigenvalue-minimisations  :  3552
 total energy-change (2. order) :-0.1873325E+03  (-0.7804744E+02)
 number of electron      86.9999985 magnetization      38.3928610
 augmentation part        2.1352294 magnetization       4.8751201

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3642
  0.8513  0.8513  0.8883  0.8883  0.7073  0.7073  0.4193  0.4193  0.4837  0.4837
  0.4183  0.4183  0.4173  0.2538  0.3149  0.3149  0.2990  0.2990  0.2725  0.2725
  0.2341  0.1541  0.1541  0.2112  0.1786  0.1005  0.1005  0.0466  0.1219  0.0018
  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3179.62238759
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       280.78280792
  PAW double counting   =      3267.68655130    -3269.57288052
  entropy T*S    EENTRO =        -0.06802332
  eigenvalues    EBANDS =     -1053.32590257
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -204.20376318 eV

  energy without entropy =     -204.13573986  energy(sigma->0) =     -204.18108874


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




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



 eigenvalue-minimisations  :  3144
 total energy-change (2. order) : 0.2426545E+03  (-0.1798984E+04)
 number of electron      86.9999981 magnetization      37.1582920
 augmentation part        2.4562206 magnetization       0.2257416

 Broyden mixing:
  rms(total) = 0.26047E+01    rms(broyden)= 0.25971E+01
  rms(prec ) = 0.28500E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3765
  0.8419  0.8419  0.8296  0.8296  0.7169  0.7169  0.4522  0.4190  0.4190  0.4425
  0.4425  0.4087  0.4087  0.4627  0.4627  0.4603  0.3385  0.3385  0.2722  0.2722
  0.2696  0.2696  0.2311  0.1541  0.1541  0.1005  0.1005  0.0466  0.1786  0.1615
  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3179.10801629
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       287.03678339
  PAW double counting   =      3023.05150343    -2955.36923587
  entropy T*S    EENTRO =        -0.02250005
  eigenvalues    EBANDS =      -887.05385321
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        38.45075299 eV

  energy without entropy =       38.47325304  energy(sigma->0) =       38.45825301


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




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



 eigenvalue-minimisations  :  3344
 total energy-change (2. order) :-0.4727962E+03  (-0.3210266E+03)
 number of electron      86.9999965 magnetization      37.2862425
 augmentation part        2.1679526 magnetization       6.6105223

 Broyden mixing:
  rms(total) = 0.28668E+02    rms(broyden)= 0.28667E+02
  rms(prec ) = 0.28766E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3747
  0.8357  0.8357  0.8038  0.8038  0.7260  0.7260  0.4952  0.4190  0.4190  0.4872
  0.4872  0.4115  0.4115  0.4814  0.4814  0.4508  0.3412  0.3412  0.2767  0.2767
  0.2728  0.2728  0.2348  0.1541  0.1541  0.1005  0.1005  0.0466  0.1786  0.1993
  0.1328  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3188.86207289
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       265.06176632
  PAW double counting   =      3273.05379090    -3254.82583145
  entropy T*S    EENTRO =         0.00950971
  eigenvalues    EBANDS =     -1278.69873049
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -434.34549631 eV

  energy without entropy =     -434.35500602  energy(sigma->0) =     -434.34866622


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




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



 eigenvalue-minimisations  :  3192
 total energy-change (2. order) : 0.4097714E+03  (-0.8038451E+02)
 number of electron      86.9999949 magnetization      37.2601728
 augmentation part        2.2462612 magnetization      -2.9834762

 Broyden mixing:
  rms(total) = 0.22862E+01    rms(broyden)= 0.22793E+01
  rms(prec ) = 0.26107E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3909
  0.8266  0.8266  0.7991  0.7991  0.7262  0.7262  0.6799  0.6799  0.5878  0.4190
  0.4190  0.5277  0.5277  0.4163  0.4163  0.4128  0.4128  0.4384  0.3011  0.3011
  0.2739  0.2739  0.2477  0.2265  0.1541  0.1541  0.1786  0.1005  0.1005  0.0466
  0.1415  0.1415  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3188.48275123
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       268.64390132
  PAW double counting   =      3348.14635545    -3325.45186332
  entropy T*S    EENTRO =         0.00710254
  eigenvalues    EBANDS =      -877.35290902
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -24.57409265 eV

  energy without entropy =      -24.58119519  energy(sigma->0) =      -24.57646016


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




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



 eigenvalue-minimisations  :  3192
 total energy-change (2. order) :-0.1346342E+02  (-0.1679868E+01)
 number of electron      86.9999949 magnetization      37.2819167
 augmentation part        2.1025789 magnetization       1.4527816

 Broyden mixing:
  rms(total) = 0.20769E+01    rms(broyden)= 0.20745E+01
  rms(prec ) = 0.24236E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3950
  0.8411  0.8411  0.7694  0.7694  0.8199  0.8199  0.6195  0.7168  0.7168  0.4190
  0.4190  0.4958  0.4958  0.4168  0.4168  0.4479  0.4479  0.4388  0.3090  0.3090
  0.2734  0.2734  0.2614  0.2274  0.1541  0.1541  0.1925  0.1925  0.1786  0.1005
  0.1005  0.0466  0.1318  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3191.49575131
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       260.65470650
  PAW double counting   =      3397.27715513    -3374.74478505
  entropy T*S    EENTRO =        -0.02910331
  eigenvalues    EBANDS =      -879.61580917
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -38.03751560 eV

  energy without entropy =      -38.00841230  energy(sigma->0) =      -38.02781450


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




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



 eigenvalue-minimisations  :  3000
 total energy-change (2. order) : 0.7332167E+01  (-0.3893974E+02)
 number of electron      86.9999947 magnetization      36.9785514
 augmentation part        2.1385978 magnetization       2.1090320

 Broyden mixing:
  rms(total) = 0.11309E+02    rms(broyden)= 0.11308E+02
  rms(prec ) = 0.11344E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4153
  0.8965  0.8965  0.8987  0.8987  0.6538  0.8231  0.8231  0.7115  0.7115  0.4190
  0.4190  0.5113  0.5113  0.4153  0.4153  0.5010  0.5010  0.4368  0.4124  0.4124
  0.3056  0.3056  0.2746  0.2746  0.2458  0.1541  0.1541  0.2257  0.1005  0.1005
  0.0466  0.1784  0.1771  0.1335  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3189.99952331
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       262.37730398
  PAW double counting   =      3406.31119049    -3383.45858900
  entropy T*S    EENTRO =        -0.00014375
  eigenvalues    EBANDS =      -875.85165897
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -30.70534896 eV

  energy without entropy =      -30.70520522  energy(sigma->0) =      -30.70530105


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




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



 eigenvalue-minimisations  :  3072
 total energy-change (2. order) :-0.6313274E+01  (-0.2421666E+01)
 number of electron      86.9999964 magnetization      37.1444167
 augmentation part        2.2140322 magnetization       2.7740524

 Broyden mixing:
  rms(total) = 0.15008E+01    rms(broyden)= 0.14940E+01
  rms(prec ) = 0.16539E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4204
  0.9224  0.9224  0.9499  0.9499  0.6626  0.8553  0.8553  0.7023  0.7023  0.4190
  0.4190  0.5661  0.5661  0.4147  0.4147  0.4984  0.4984  0.4514  0.4135  0.4135
  0.3048  0.3048  0.2745  0.2745  0.2909  0.1541  0.1541  0.0466  0.1005  0.1005
  0.2310  0.2101  0.1787  0.1912  0.1332  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3190.68514664
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.50044848
  PAW double counting   =      3379.65473543    -3357.10922464
  entropy T*S    EENTRO =        -0.07509750
  eigenvalues    EBANDS =      -876.22040979
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.01862306 eV

  energy without entropy =      -36.94352556  energy(sigma->0) =      -36.99359056


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




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



 eigenvalue-minimisations  :  2576
 total energy-change (2. order) : 0.1017969E+01  (-0.2800105E+02)
 number of electron      86.9999970 magnetization      37.1641348
 augmentation part        2.2916431 magnetization       3.9942214

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4097
  0.9222  0.9222  0.9496  0.9496  0.6626  0.8566  0.8566  0.7025  0.7025  0.4190
  0.4190  0.5662  0.5662  0.4978  0.4978  0.4147  0.4147  0.4140  0.4140  0.4510
  0.3048  0.3048  0.2745  0.2745  0.2927  0.1541  0.1541  0.2313  0.1005  0.1005
  0.2118  0.1787  0.1900  0.0466  0.1332  0.0094  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3189.15110757
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.68777530
  PAW double counting   =      3380.59479829    -3358.04998630
  entropy T*S    EENTRO =        -0.05709048
  eigenvalues    EBANDS =      -876.94111491
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -36.00065407 eV

  energy without entropy =      -35.94356359  energy(sigma->0) =      -35.98162391


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




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



 eigenvalue-minimisations  :  2680
 total energy-change (2. order) : 0.2142273E+02  (-0.2004109E+03)
 number of electron      86.9999971 magnetization      37.1864270
 augmentation part        2.3750875 magnetization       1.7166411

 Broyden mixing:
  rms(total) = 0.15986E+01    rms(broyden)= 0.15980E+01
  rms(prec ) = 0.17002E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4047
  0.9239  0.9239  0.9512  0.9512  0.6633  0.8585  0.8585  0.7017  0.7017  0.4190
  0.4190  0.5637  0.5637  0.4146  0.4146  0.4970  0.4970  0.4480  0.4112  0.4112
  0.3047  0.3047  0.2745  0.2745  0.2914  0.1541  0.1541  0.2318  0.1005  0.1005
  0.0466  0.2130  0.1787  0.1879  0.1164  0.1164  0.1332  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3189.12169036
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       258.11255425
  PAW double counting   =      3394.52970860    -3345.25234928
  entropy T*S    EENTRO =        -0.04004767
  eigenvalues    EBANDS =      -882.72217202
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -14.57792489 eV

  energy without entropy =      -14.53787722  energy(sigma->0) =      -14.56457566


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




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



 eigenvalue-minimisations  :  3032
 total energy-change (2. order) :-0.2288747E+02  (-0.3994976E+00)
 number of electron      86.9999966 magnetization      37.1809656
 augmentation part        2.3618340 magnetization       1.6714089

 Broyden mixing:
  rms(total) = 0.15171E+01    rms(broyden)= 0.15165E+01
  rms(prec ) = 0.16143E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4055
  0.9220  0.9220  0.9525  0.9525  0.8654  0.8654  0.6614  0.7016  0.7016  0.4190
  0.4190  0.5655  0.5655  0.4973  0.4973  0.4147  0.4147  0.4447  0.4067  0.4067
  0.2591  0.3039  0.3039  0.2745  0.2745  0.2906  0.2068  0.2068  0.1541  0.1541
  0.0466  0.1005  0.1005  0.2328  0.2136  0.1787  0.1851  0.1333  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3189.24718550
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       256.96306318
  PAW double counting   =      3374.21363538    -3352.01926396
  entropy T*S    EENTRO =         0.01315479
  eigenvalues    EBANDS =      -877.30486661
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -37.46539113 eV

  energy without entropy =      -37.47854592  energy(sigma->0) =      -37.46977606


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




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



 eigenvalue-minimisations  :  3776
 total energy-change (2. order) :-0.1115115E+04  (-0.1036127E+04)
 number of electron      86.9999959 magnetization      39.6269132
 augmentation part        1.8159006 magnetization      12.0060659

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3952
  0.9221  0.9221  0.9526  0.9526  0.8589  0.8589  0.6615  0.7021  0.7021  0.4190
  0.4190  0.5647  0.5647  0.5001  0.5001  0.4147  0.4147  0.4433  0.4069  0.4069
  0.2544  0.3040  0.3040  0.2744  0.2744  0.2804  0.2092  0.2092  0.1541  0.1541
  0.0466  0.1005  0.1005  0.2313  0.2082  0.1787  0.1870  0.1333  0.0018  0.0050
  0.0042

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3189.07434028
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       257.01332831
  PAW double counting   =      3372.90143807    -3353.32294079
  entropy T*S    EENTRO =        -0.04108812
  eigenvalues    EBANDS =     -1989.97316576
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -1152.58069698 eV

  energy without entropy =    -1152.53960886  energy(sigma->0) =    -1152.56700094


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




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



 eigenvalue-minimisations  :  2904
 total energy-change (2. order) : 0.1221295E+04  (-0.4298319E+03)
 number of electron      86.9999990 magnetization      37.5213664
 augmentation part        2.9675099 magnetization       6.2086617

 Broyden mixing:
  rms(total) = 0.68874E+01    rms(broyden)= 0.67729E+01
  rms(prec ) = 0.70626E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3919
  0.9258  0.9258  0.9595  0.9595  0.6642  0.8647  0.8647  0.6991  0.6991  0.4190
  0.4190  0.5657  0.5657  0.4963  0.4963  0.4145  0.4145  0.4471  0.3975  0.3975
  0.2054  0.2054  0.3028  0.3028  0.2743  0.2743  0.2754  0.2754  0.2477  0.2250
  0.1541  0.1541  0.1790  0.1790  0.1005  0.1005  0.0466  0.1332  0.1100  0.1100
  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3197.04899766
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       350.60325041
  PAW double counting   =      3582.17205360    -3553.61184469
  entropy T*S    EENTRO =        -0.03407005
  eigenvalues    EBANDS =      -863.28237912
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        68.71408408 eV

  energy without entropy =       68.74815413  energy(sigma->0) =       68.72544076


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




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



 eigenvalue-minimisations  :  3224
 total energy-change (2. order) :-0.2099692E+04  (-0.1945074E+04)
 number of electron      86.9999992 magnetization      37.5578728
 augmentation part        2.6821679 magnetization       8.7893763

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3827
  0.9254  0.9254  0.9599  0.9599  0.8677  0.8677  0.6641  0.6993  0.6993  0.4190
  0.4190  0.5647  0.5647  0.4934  0.4934  0.4145  0.4145  0.4464  0.3973  0.3973
  0.2036  0.2036  0.3024  0.3024  0.2743  0.2743  0.2761  0.2761  0.2481  0.2251
  0.1541  0.1541  0.0466  0.1005  0.1005  0.1789  0.1789  0.1332  0.1107  0.1107
  0.0050  0.0018  0.0004

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3197.79157602
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       271.86738008
  PAW double counting   =      3687.99328338    -3661.35285109
  entropy T*S    EENTRO =        -0.02985483
  eigenvalues    EBANDS =     -2881.58006596
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -2030.97761284 eV

  energy without entropy =    -2030.94775801  energy(sigma->0) =    -2030.96766123


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




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



 eigenvalue-minimisations  :  3696
 total energy-change (2. order) :-0.4778854E+04  (-0.6254022E+04)
 number of electron      87.0000004 magnetization      36.3192187
 augmentation part        2.3697636 magnetization      10.6171527

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3746
  0.9254  0.9254  0.9572  0.9572  0.8742  0.8742  0.6643  0.6994  0.6994  0.4190
  0.4190  0.5637  0.5637  0.4857  0.4857  0.4145  0.4145  0.4464  0.4015  0.4015
  0.2234  0.2234  0.3043  0.3043  0.2739  0.2739  0.2532  0.2532  0.2523  0.2256
  0.1541  0.1541  0.1005  0.1005  0.1783  0.1783  0.0466  0.1331  0.1083  0.1083
  0.0330  0.0015  0.0018  0.0050

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3197.76266100
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       293.31488531
  PAW double counting   =      3710.90066173    -3689.26758567
  entropy T*S    EENTRO =        -0.03381574
  eigenvalues    EBANDS =     -7676.89901934
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -6809.83146311 eV

  energy without entropy =    -6809.79764738  energy(sigma->0) =    -6809.82019120


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




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



 eigenvalue-minimisations  :  3200
 total energy-change (2. order) : 0.4438811E+04  (-0.5420170E+03)
 number of electron      87.0000050 magnetization      36.6107490
 augmentation part        1.6701866 magnetization       7.9161098

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3554
  0.9094  0.9094  0.8591  0.8591  0.8276  0.7085  0.7085  0.4722  0.4722  0.3730
  0.5179  0.5179  0.3263  0.3263  0.4248  0.4248  0.1890  0.4465  0.3743  0.3743
  0.3442  0.2971  0.2971  0.1787  0.1787  0.2035  0.2035  0.2240  0.2240  0.1658
  0.1658  0.0020  0.0020  0.0016  0.0461  0.0526  0.1810  0.1578  0.1578  0.1135

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         6.34847707
  Ewald energy   TEWEN  =       239.48782979
  -Hartree energ DENC   =     -3229.48271449
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       348.65137822
  PAW double counting   =      3877.97876757    -3844.67758568
  entropy T*S    EENTRO =        -0.02287875
  eigenvalues    EBANDS =     -3273.38324706
  atomic energy  EATOM  =      3504.07976474
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -2371.02020859 eV

  energy without entropy =    -2370.99732984  energy(sigma->0) =    -2371.01258234


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




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


 -----------------------------------------------------------------------------
|                                                                             |
|     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!!! <----       |
|                                                                             |
 -----------------------------------------------------------------------------