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.02  18:57:09
 running on    3 total cores
 distrk:  each k-point on    3 cores,    1 groups
 distr:  one band on NCORE=   1 cores,    3 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 = 200
   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 = 1.11 0.32 0.73
   NPAR = 3

 POTCAR:    PAW_PBE Si 05Jan2001                  
 POTCAR:    PAW_PBE H 15Jun2001                   
 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 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 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 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         = 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         = 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         = 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 Si 05Jan2001                  :
 energy of atom  1       EATOM= -103.0669
 kinetic energy error for atom=    0.0012 (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 O 08Apr2002                   :
 energy of atom  3       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.977  0.087  0.601-  20 1.73
   2  0.050  0.262  0.659-  18 2.06
   3  0.924  0.498  0.509-  13 1.97
   4  0.539  0.048  0.196-   6 2.27  16 2.31
   5  0.556  0.953  0.524-  29 0.38
   6  0.352  0.134  0.244-   4 2.27
   7  0.591  0.396  0.465-  12 2.48
   8  0.858  0.370  0.740-  18 1.58
   9  0.783  0.729  0.377-  14 1.49  26 1.74
  10  0.242  0.979  0.151-  30 1.26
  11  0.400  0.180  0.606-
  12  0.301  0.358  0.391-   7 2.48
  13  0.108  0.568  0.501-  28 1.80   3 1.97
  14  0.703  0.795  0.419-   9 1.49  26 2.65
  15  0.891  0.593  0.270-  28 1.91
  16  0.722  0.959  0.245-   4 2.31
  17  0.942  0.671  0.769-
  18  0.062  0.360  0.725-   8 1.58   2 2.06
  19  0.756  0.207  0.026-  22 2.58
  20  0.138  0.058  0.502-  34 0.84   1 1.73
  21  0.040  0.205  0.429-
  22  0.976  0.108  0.004-  30 2.52  19 2.58
  23  0.979  0.433  0.181-  32 2.65
  24  0.025  0.736  0.044-  25 1.53
  25  0.107  0.770  0.932-  24 1.53
  26  0.968  0.711  0.464-   9 1.74  28 1.96  14 2.65
  27  0.882  0.807  0.638-
  28  0.058  0.630  0.383-  13 1.80  15 1.91  26 1.96
  29  0.543  0.935  0.514-   5 0.38
  30  0.199  0.021  0.069-  10 1.26  22 2.52
  31  0.552  0.606  0.623-
  32  0.221  0.517  0.265-  35 2.00  23 2.65
  33  0.454  0.754  0.860-
  34  0.149  0.099  0.484-  20 0.84
  35  0.406  0.572  0.182-  32 2.00
 
 -----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|     The distance between some ions is very small. Please check the          |
|     nearest-neighbor list in the OUTCAR file.                               |
|     I HOPE YOU KNOW WHAT YOU ARE DOING!                                     |
|                                                                             |
 -----------------------------------------------------------------------------

  LATTYP: Found a simple orthorhombic cell.
 ALAT       =     7.6631000000
 B/A-ratio  =     1.4142135689
 C/A-ratio  =     2.5659458966
  
  Lattice vectors:
  
 A1 = (  -7.6631000000,   0.0000000000,   0.0000000000)
 A2 = (   0.0000000000,   0.0000000000,  10.8372600000)
 A3 = (   0.0000000000,  19.6631000000,   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 orthorhombic supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of  8 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 orthorhombic supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of  8 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 orthorhombic supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of  8 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 :    1632.9616

  direct lattice vectors                    reciprocal lattice vectors
     7.663100000  0.000000000  0.000000000     0.130495491  0.000000000  0.000000000
     0.000000000 19.663100000  0.000000000     0.000000000  0.050856681  0.000000000
     0.000000000  0.000000000 10.837260000     0.000000000  0.000000000  0.092274246

  length of vectors
     7.663100000 19.663100000 10.837260000     0.130495491  0.050856681  0.092274246

  position of ions in fractional coordinates (direct lattice)
     0.976591790  0.086780240  0.600806050
     0.050350270  0.262156190  0.659128300
     0.924244260  0.498392990  0.509298650
     0.539370920  0.048074950  0.196081890
     0.555643030  0.952909820  0.524047320
     0.352368860  0.133754800  0.244098250
     0.591444270  0.395893690  0.464501480
     0.857793000  0.369647870  0.740209060
     0.782691520  0.729411700  0.377187140
     0.241723200  0.979232670  0.151433330
     0.399855240  0.179568000  0.605673910
     0.301223430  0.357642600  0.391302900
     0.107686950  0.568321640  0.500503210
     0.703151650  0.794556910  0.419348930
     0.891034310  0.593148910  0.270043730
     0.722303910  0.958677690  0.244938860
     0.942316810  0.671357730  0.769111960
     0.061909970  0.360439790  0.725270920
     0.756059160  0.206889280  0.025876600
     0.137925970  0.057547370  0.502280000
     0.040018910  0.205412460  0.429494320
     0.975699170  0.107974770  0.003503110
     0.979256270  0.432800560  0.181167090
     0.024576520  0.735772210  0.043950240
     0.107045360  0.770306820  0.931505630
     0.968221080  0.711497400  0.463644560
     0.881731980  0.806553250  0.638099050
     0.058407320  0.629545480  0.382570960
     0.543214730  0.934782430  0.514289320
     0.198827210  0.020903820  0.068652960
     0.552029530  0.606348530  0.623091630
     0.220908880  0.517483650  0.264984100
     0.454327150  0.754230550  0.860025470
     0.149047940  0.098770290  0.484219880
     0.406493100  0.571884380  0.181571440

  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

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

 
 
 KPOINTS: Automatic mesh                          

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

 Generating k-lattice:

  Cartesian coordinates                     Fractional coordinates (reciprocal lattice)
     0.065247746  0.000000000  0.000000000     0.500000000  0.000000000  0.000000000
     0.000000000  0.050856681  0.000000000     0.000000000  1.000000000  0.000000000
     0.000000000  0.000000000  0.046137123     0.000000000  0.000000000  0.500000000

  Length of vectors
     0.065247746  0.050856681  0.046137123

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

 
 Subroutine IBZKPT returns following result:
 ===========================================
 
 Found      4 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.000000  0.500000      1.000000
  0.500000  0.000000  0.500000      1.000000
 
 Following cartesian coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
  0.065248  0.000000  0.000000      1.000000
  0.000000  0.000000  0.046137      1.000000
  0.065248  0.000000  0.046137      1.000000
 


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




 Dimension of arrays:
   k-points           NKPTS =      4   k-points in BZ     NKDIM =      4   number of bands    NBANDS=    105
   number of dos      NEDOS =    301   number of ions     NIONS =     35
   non local maximal  LDIM  =      4   non local SUM 2l+1 LMDIM =      8
   total plane-waves  NPLWV = 207360
   max r-space proj   IRMAX =   1449   max aug-charges    IRDMAX=   4445
   dimension x,y,z NGX =    40 NGY =   96 NGZ =   54
   dimension x,y,z NGXF=    80 NGYF=  192 NGZF=  108
   support grid    NGXF=    80 NGYF=  192 NGZF=  108
   ions per type =              32   2   1
   NGX,Y,Z   is equivalent  to a cutoff of   8.68,  8.12,  8.28 a.u.
   NGXF,Y,Z  is equivalent  to a cutoff of  17.36, 16.23, 16.57 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.  12.50 32.07 17.67*2*pi/ulx,y,z
   ENINI  =  400.0     initial cutoff
   ENAUG  =  605.4 eV  augmentation charge cutoff
   NELM   =    200;   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
 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.134E-26a.u.
   SCALEE = 1.0000    scale energy and forces
   NPACO  =    256;   APACO  = 16.0  distance and # of slots for P.C.
   PSTRESS=    0.0 pullay stress

  Mass of Ions in am
   POMASS =  28.09  1.00 16.00
  Ionic Valenz
   ZVAL   =   4.00  1.00  6.00
  Atomic Wigner-Seitz radii
   RWIGS  =   1.11  0.32  0.73
  virtual crystal weights 
   VCA    =   1.00  1.00  1.00
   NELECT =     136.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.24E-07  absolut break condition
   DEPER  =   0.30     relativ break condition  

   TIME   =   0.40     timestep for ELM

  volume/ion in A,a.u.               =      46.66       314.85
  Fermi-wavevector in a.u.,A,eV,Ry     =   0.714928  1.351019  6.954243  0.511122
  Thomas-Fermi vector in A             =   1.802956
 
 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           37
 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 :     1632.96
      direct lattice vectors                 reciprocal lattice vectors
     7.663100000  0.000000000  0.000000000     0.130495491  0.000000000  0.000000000
     0.000000000 19.663100000  0.000000000     0.000000000  0.050856681  0.000000000
     0.000000000  0.000000000 10.837260000     0.000000000  0.000000000  0.092274246

  length of vectors
     7.663100000 19.663100000 10.837260000     0.130495491  0.050856681  0.092274246


 
 k-points in units of 2pi/SCALE and weight: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       0.250
   0.06524775  0.00000000  0.00000000       0.250
   0.00000000  0.00000000  0.04613712       0.250
   0.06524775  0.00000000  0.04613712       0.250
 
 k-points in reciprocal lattice and weights: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       0.250
   0.50000000  0.00000000  0.00000000       0.250
   0.00000000  0.00000000  0.50000000       0.250
   0.50000000  0.00000000  0.50000000       0.250
 
 position of ions in fractional coordinates (direct lattice) 
   0.97659179  0.08678024  0.60080605
   0.05035027  0.26215619  0.65912830
   0.92424426  0.49839299  0.50929865
   0.53937092  0.04807495  0.19608189
   0.55564303  0.95290982  0.52404732
   0.35236886  0.13375480  0.24409825
   0.59144427  0.39589369  0.46450148
   0.85779300  0.36964787  0.74020906
   0.78269152  0.72941170  0.37718714
   0.24172320  0.97923267  0.15143333
   0.39985524  0.17956800  0.60567391
   0.30122343  0.35764260  0.39130290
   0.10768695  0.56832164  0.50050321
   0.70315165  0.79455691  0.41934893
   0.89103431  0.59314891  0.27004373
   0.72230391  0.95867769  0.24493886
   0.94231681  0.67135773  0.76911196
   0.06190997  0.36043979  0.72527092
   0.75605916  0.20688928  0.02587660
   0.13792597  0.05754737  0.50228000
   0.04001891  0.20541246  0.42949432
   0.97569917  0.10797477  0.00350311
   0.97925627  0.43280056  0.18116709
   0.02457652  0.73577221  0.04395024
   0.10704536  0.77030682  0.93150563
   0.96822108  0.71149740  0.46364456
   0.88173198  0.80655325  0.63809905
   0.05840732  0.62954548  0.38257096
   0.54321473  0.93478243  0.51428932
   0.19882721  0.02090382  0.06865296
   0.55202953  0.60634853  0.62309163
   0.22090888  0.51748365  0.26498410
   0.45432715  0.75423055  0.86002547
   0.14904794  0.09877029  0.48421988
   0.40649310  0.57188438  0.18157144
 
 position of ions in cartesian coordinates  (Angst):
   7.48372055  1.70636854  6.51109137
   0.38583915  5.15480338  7.14314476
   7.08257619  9.79995120  5.51940189
   4.13325330  0.94530255  2.12499042
   4.25794810 18.73716108  5.67923706
   2.70023781  2.63003401  2.64535620
   4.53229659  7.78449722  5.03392331
   6.57335354  7.26842303  8.02183804
   5.99784339 14.34249520  4.08767510
   1.85234905 19.25474991  1.64112237
   3.06413069  3.53086354  6.56384564
   2.30830527  7.03236221  4.24065127
   0.82521587 11.17496524  5.42408342
   5.38832141 15.62345198  4.54459339
   6.82808502 11.66314633  2.92653411
   5.53508709 18.85057529  2.65446611
   7.22106795 13.20097418  8.33506628
   0.47442229  7.08736363  7.85994953
   5.79375695  4.06808460  0.28043144
   1.05694050  1.13155969  5.44333895
   0.30666891  4.03904574  4.65454161
   7.47688031  2.12311870  0.03796411
   7.50413872  8.51020069  1.96335486
   0.18833233 14.46756254  0.47630018
   0.82029930 15.14662003 10.09496870
   7.41957496 13.99024453  5.02463664
   6.75680034 15.85933721  6.91524531
   0.44758113 12.37881573  4.14602096
   4.16270880 18.38072040  5.57348708
   1.52363279  0.41103390  0.74400998
   4.23025749 11.92269178  6.75260600
   1.69284684 10.17533276  2.87170159
   3.48155438 14.83051073  9.32031963
   1.14216927  1.94213009  5.24761674
   3.11499727 11.24501975  1.96773690
 


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


 k-point  1 :   0.0000 0.0000 0.0000  plane waves:   29735
 k-point  2 :   0.5000 0.0000 0.0000  plane waves:   29568
 k-point  3 :   0.0000 0.0000 0.5000  plane waves:   29692
 k-point  4 :   0.5000 0.0000 0.5000  plane waves:   29632

 maximum and minimum number of plane-waves per node :     29735    29568

 maximum number of plane-waves:     29735
 maximum index in each direction: 
   IXMAX=   12   IYMAX=   32   IZMAX=   17
   IXMIN=  -12   IYMIN=  -32   IZMIN=  -18

 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    50 to avoid them
 WARNING: aliasing errors must be expected set NGY to   140 to avoid them
 WARNING: aliasing errors must be expected set NGZ to    72 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   254503. kBytes
=======================================================================

   base      :      30000. kBytes
   nonlr-proj:       3171. kBytes
   fftplans  :      19938. kBytes
   grid      :      66726. kBytes
   one-center:        215. kBytes
   wavefun   :     134453. kBytes
 
 Broyden mixing: mesh for mixing (old mesh)
   NGX = 25   NGY = 65   NGZ = 35
  (NGX  = 80   NGY  =192   NGZ  =108)
  gives a total of  56875 points

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


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


 Maximum index for non-local projection operator         1360
 Maximum index for augmentation-charges         1480 (set IRDMAX)


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


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



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



 eigenvalue-minimisations  :  1680
 total energy-change (2. order) : 0.9859646E+03  (-0.3820912E+04)
 number of electron     136.0000000 magnetization      35.0000000
 augmentation part      136.0000000 magnetization      35.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1730.95355737
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -290.06711855
  PAW double counting   =      3167.08814863    -2605.71786812
  entropy T*S    EENTRO =        -0.00678875
  eigenvalues    EBANDS =      -216.67222978
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       985.96460908 eV

  energy without entropy =      985.97139783  energy(sigma->0) =      985.96687200


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




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



 eigenvalue-minimisations  :  2280
 total energy-change (2. order) :-0.5745488E+03  (-0.5408739E+03)
 number of electron     136.0000000 magnetization      35.0000000
 augmentation part      136.0000000 magnetization      35.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1730.95355737
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -290.06711855
  PAW double counting   =      3167.08814863    -2605.71786812
  entropy T*S    EENTRO =        -0.00415363
  eigenvalues    EBANDS =      -791.22368195
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       411.41579202 eV

  energy without entropy =      411.41994566  energy(sigma->0) =      411.41717657


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




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



 eigenvalue-minimisations  :  2586
 total energy-change (2. order) :-0.1013562E+03  (-0.9239580E+02)
 number of electron     136.0000000 magnetization      35.0000000
 augmentation part      136.0000000 magnetization      35.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1730.95355737
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -290.06711855
  PAW double counting   =      3167.08814863    -2605.71786812
  entropy T*S    EENTRO =         0.02339725
  eigenvalues    EBANDS =      -892.60740043
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       310.05962443 eV

  energy without entropy =      310.03622718  energy(sigma->0) =      310.05182535


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




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



 eigenvalue-minimisations  :  2577
 total energy-change (2. order) :-0.1168687E+02  (-0.1109757E+02)
 number of electron     136.0000000 magnetization      35.0000000
 augmentation part      136.0000000 magnetization      35.0000000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1730.95355737
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -290.06711855
  PAW double counting   =      3167.08814863    -2605.71786812
  entropy T*S    EENTRO =         0.01105384
  eigenvalues    EBANDS =      -904.28192335
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       298.37275811 eV

  energy without entropy =      298.36170426  energy(sigma->0) =      298.36907349


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




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



 eigenvalue-minimisations  :  2280
 total energy-change (2. order) :-0.4056170E+00  (-0.4015197E+00)
 number of electron     136.0000005 magnetization      31.4184509
 augmentation part       -8.3316882 magnetization      20.4043409

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1730.95355737
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -290.06711855
  PAW double counting   =      3167.08814863    -2605.71786812
  entropy T*S    EENTRO =         0.00481601
  eigenvalues    EBANDS =      -904.68130250
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       297.96714112 eV

  energy without entropy =      297.96232511  energy(sigma->0) =      297.96553579


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




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



 eigenvalue-minimisations  :  2427
 total energy-change (2. order) : 0.1302349E+01  (-0.6783260E+02)
 number of electron     136.0000028 magnetization      29.5161277
 augmentation part       -8.3565492 magnetization      16.4116040

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

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1935.53720363
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -265.73014437
  PAW double counting   =   1016326.97031698 -1015770.54644138
  entropy T*S    EENTRO =         0.01905706
  eigenvalues    EBANDS =      -718.20011749
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       299.26949019 eV

  energy without entropy =      299.25043313  energy(sigma->0) =      299.26313784


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




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



 eigenvalue-minimisations  :  2466
 total energy-change (2. order) : 0.4366322E+02  (-0.2687182E+02)
 number of electron     136.0000006 magnetization      27.2043949
 augmentation part       -7.7643688 magnetization      15.4445372

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.9317
  1.3138  0.5497

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1846.15830000
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -269.32075102
  PAW double counting   =    544359.17277928  -543802.64245344
  entropy T*S    EENTRO =         0.02448767
  eigenvalues    EBANDS =      -760.43707134
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       342.93271418 eV

  energy without entropy =      342.90822650  energy(sigma->0) =      342.92455162


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




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



 eigenvalue-minimisations  :  2502
 total energy-change (2. order) : 0.2707736E+02  (-0.1840706E+02)
 number of electron     136.0000001 magnetization      27.1118737
 augmentation part       -7.5383255 magnetization      14.2423753

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6301
  1.2700  0.5564  0.0639

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1893.35009762
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -263.56512983
  PAW double counting   =    670146.82987369  -669591.37596400
  entropy T*S    EENTRO =         0.02782168
  eigenvalues    EBANDS =      -690.85045723
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       370.01006970 eV

  energy without entropy =      369.98224802  energy(sigma->0) =      370.00079581


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




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



 eigenvalue-minimisations  :  2388
 total energy-change (2. order) :-0.2985578E+01  (-0.4725996E+01)
 number of electron     136.0000001 magnetization      26.9970931
 augmentation part       -7.8422905 magnetization      13.2546533

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5157
  1.1901  0.5633  0.2262  0.0829

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1893.35399694
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -263.57347372
  PAW double counting   =    678211.55783590  -677656.12586628
  entropy T*S    EENTRO =         0.00801288
  eigenvalues    EBANDS =      -693.78204357
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       367.02449128 eV

  energy without entropy =      367.01647841  energy(sigma->0) =      367.02182032


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




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



 eigenvalue-minimisations  :  2361
 total energy-change (2. order) : 0.3945271E+01  (-0.5755133E+00)
 number of electron     136.0000001 magnetization      27.7623643
 augmentation part       -7.9093172 magnetization      14.0196501

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5371
  1.1139  0.5744  0.3865  0.3865  0.2244

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1887.11786987
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -262.73016163
  PAW double counting   =    666596.13996817  -666040.46154949
  entropy T*S    EENTRO =        -0.01179106
  eigenvalues    EBANDS =      -697.14285722
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       370.96976192 eV

  energy without entropy =      370.98155298  energy(sigma->0) =      370.97369228


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




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



 eigenvalue-minimisations  :  2073
 total energy-change (2. order) :-0.2553531E+01  (-0.1419522E+01)
 number of electron     136.0000002 magnetization      26.3202547
 augmentation part       -7.8354564 magnetization      12.4782384

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5861
  1.1552  0.5014  0.4838  0.4838  0.5695  0.3232

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1891.19537373
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -264.40124460
  PAW double counting   =    659738.58377592  -659183.36570197
  entropy T*S    EENTRO =        -0.02960458
  eigenvalues    EBANDS =      -693.46964307
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       368.41623097 eV

  energy without entropy =      368.44583555  energy(sigma->0) =      368.42609916


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




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



 eigenvalue-minimisations  :  2073
 total energy-change (2. order) :-0.9183411E+01  (-0.4038795E+01)
 number of electron     136.0000002 magnetization      23.5702495
 augmentation part       -7.9361159 magnetization      13.5846281

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6677
  0.6777  1.1735  0.7760  0.7760  0.5621  0.3543  0.3543

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1871.48786326
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -262.32048809
  PAW double counting   =    655755.64349942  -655199.53036895
  entropy T*S    EENTRO =        -0.03705527
  eigenvalues    EBANDS =      -725.32892683
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       359.23282004 eV

  energy without entropy =      359.26987531  energy(sigma->0) =      359.24517179


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




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



 eigenvalue-minimisations  :  1893
 total energy-change (2. order) : 0.3061168E+02  (-0.1545945E+02)
 number of electron     136.0000008 magnetization      16.1040496
 augmentation part       -8.2149838 magnetization       7.7733132

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7386
  0.7649  1.1675  1.1675  1.0989  0.5577  0.4099  0.4099  0.3329

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1853.38533180
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -247.74090330
  PAW double counting   =    647873.95735239  -647317.18810294
  entropy T*S    EENTRO =         0.01449678
  eigenvalues    EBANDS =      -728.10703544
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       389.84449871 eV

  energy without entropy =      389.83000193  energy(sigma->0) =      389.83966645


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




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



 eigenvalue-minimisations  :  2079
 total energy-change (2. order) : 0.3245388E+02  (-0.2108352E+02)
 number of electron     136.0000018 magnetization      15.0854061
 augmentation part       -9.1457820 magnetization      11.4681392

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6816
  0.7758  1.2261  1.2261  1.0812  0.5584  0.4145  0.4145  0.3457  0.0919

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1830.58452818
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -185.87131881
  PAW double counting   =    646355.90698044  -645795.90359179
  entropy T*S    EENTRO =         0.03886795
  eigenvalues    EBANDS =      -783.58205269
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       422.29837994 eV

  energy without entropy =      422.25951198  energy(sigma->0) =      422.28542395


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




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



 eigenvalue-minimisations  :  2568
 total energy-change (2. order) : 0.1379838E+01  (-0.3409705E+01)
 number of electron     136.0000018 magnetization      15.2752732
 augmentation part       -9.2146132 magnetization      12.7057076

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6293
  0.7786  1.2227  1.2227  1.0837  0.5622  0.4106  0.4106  0.3472  0.1273  0.1273

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1821.30436672
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -177.11710888
  PAW double counting   =    633592.14986600  -633033.49136080
  entropy T*S    EENTRO =         0.03953076
  eigenvalues    EBANDS =      -798.89236551
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       423.67821787 eV

  energy without entropy =      423.63868711  energy(sigma->0) =      423.66504095


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




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



 eigenvalue-minimisations  :  2145
 total energy-change (2. order) :-0.2103865E+00  (-0.3871790E+00)
 number of electron     136.0000017 magnetization      14.6466551
 augmentation part       -9.2335420 magnetization      12.4559059

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5971
  1.2208  1.2208  0.7815  1.0994  0.5613  0.4148  0.4148  0.3430  0.1848  0.1848
  0.1418

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1824.57164062
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -184.77311648
  PAW double counting   =    613098.87489042  -612539.62905755
  entropy T*S    EENTRO =         0.05396411
  eigenvalues    EBANDS =      -788.78123153
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       423.46783136 eV

  energy without entropy =      423.41386725  energy(sigma->0) =      423.44984333


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




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



 eigenvalue-minimisations  :  2028
 total energy-change (2. order) : 0.3558203E+01  (-0.2404395E+00)
 number of electron     136.0000017 magnetization      13.7651065
 augmentation part       -9.1675082 magnetization      11.7619648

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5857
  1.2252  1.2252  0.7886  1.1101  0.5599  0.3067  0.3067  0.4177  0.4177  0.3409
  0.1648  0.1648

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1821.31947044
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -178.34524492
  PAW double counting   =    611232.59045850  -610673.53686609
  entropy T*S    EENTRO =         0.05545045
  eigenvalues    EBANDS =      -794.71231655
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       427.02603398 eV

  energy without entropy =      426.97058352  energy(sigma->0) =      427.00755049


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




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



 eigenvalue-minimisations  :  2415
 total energy-change (2. order) :-0.1436326E+04  (-0.9506415E+03)
 number of electron     136.0000046 magnetization      13.7620001
 augmentation part       -9.8855518 magnetization      17.3229503

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5392
  1.2250  1.2250  1.1108  0.7883  0.5597  0.4180  0.4180  0.3410  0.3005  0.3005
  0.1607  0.1607  0.0016

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1818.65255778
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -169.32530074
  PAW double counting   =    602215.91342882  -601739.15266176
  entropy T*S    EENTRO =        -0.01362541
  eigenvalues    EBANDS =     -2160.36373492
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -1009.30042877 eV

  energy without entropy =    -1009.28680336  energy(sigma->0) =    -1009.29588697


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




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



 eigenvalue-minimisations  :  2376
 total energy-change (2. order) : 0.1205559E+04  (-0.5986034E+02)
 number of electron     136.0000047 magnetization      13.7619404
 augmentation part      -10.1120604 magnetization      10.3936964

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5006
  1.2251  1.2251  1.1108  0.7883  0.5598  0.4180  0.4180  0.3410  0.2999  0.2999
  0.1601  0.1601  0.0012  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1818.68954955
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -168.79723379
  PAW double counting   =    602113.39101336  -601572.81239991
  entropy T*S    EENTRO =         0.02923918
  eigenvalues    EBANDS =     -1019.15646861
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       196.25862369 eV

  energy without entropy =      196.22938451  energy(sigma->0) =      196.24887730


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




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



 eigenvalue-minimisations  :  2334
 total energy-change (2. order) :-0.1304765E+04  (-0.1188044E+03)
 number of electron     136.0000043 magnetization      13.6923132
 augmentation part      -10.1043022 magnetization       9.2554205

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4670
  1.2252  1.2252  1.1108  0.7882  0.5598  0.4181  0.4181  0.3410  0.2990  0.2990
  0.1593  0.1593  0.0005  0.0001  0.0011

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1818.70832337
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -168.46688156
  PAW double counting   =    602010.73441051  -601529.19780557
  entropy T*S    EENTRO =        -0.00164637
  eigenvalues    EBANDS =     -2265.15965673
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -1108.50588007 eV

  energy without entropy =    -1108.50423370  energy(sigma->0) =    -1108.50533128


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




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



 eigenvalue-minimisations  :  2445
 total energy-change (2. order) :-0.1215059E+05  (-0.9828564E+04)
 number of electron     136.0000013 magnetization      14.3191586
 augmentation part       -8.7867294 magnetization      13.9132959

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4375
  1.2249  1.2249  1.1072  0.7882  0.5602  0.4175  0.4175  0.3413  0.2955  0.2955
  0.1430  0.1430  0.0390  0.0004  0.0006  0.0006

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1819.23226868
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -167.06162278
  PAW double counting   =    600095.51355331  -599628.98734915
  entropy T*S    EENTRO =         0.00550022
  eigenvalues    EBANDS =    -14401.62341991
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -13259.09158397 eV

  energy without entropy =   -13259.09708419  energy(sigma->0) =   -13259.09341737


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




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



 eigenvalue-minimisations  :  1839
 total energy-change (2. order) : 0.1371915E+05  (-0.3482689E+03)
 number of electron     136.0000016 magnetization      13.1067940
 augmentation part       -8.9534301 magnetization      13.1825514

 Broyden mixing:
  rms(total) = 0.40478E+02    rms(broyden)= 0.40473E+02
  rms(prec ) = 0.40574E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4516
  1.2421  1.2421  0.7968  1.0764  0.4232  0.4232  0.5643  0.4095  0.4095  0.3515
  0.2193  0.2193  0.1810  0.1168  0.0004  0.0006  0.0006

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1819.12556196
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -166.46915587
  PAW double counting   =    591794.66241513  -591231.80477375
  entropy T*S    EENTRO =         0.03222059
  eigenvalues    EBANDS =      -779.53387089
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       460.05529627 eV

  energy without entropy =      460.02307567  energy(sigma->0) =      460.04455607


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




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



 eigenvalue-minimisations  :  2307
 total energy-change (2. order) :-0.1201995E+04  (-0.9076307E+03)
 number of electron     136.0000007 magnetization      13.1382945
 augmentation part      -10.0059816 magnetization      15.0570651

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4271
  1.2424  1.2424  0.7968  1.0741  0.4258  0.4258  0.5645  0.4087  0.4087  0.3526
  0.2218  0.2218  0.1807  0.1177  0.0006  0.0006  0.0004  0.0022

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1810.26480260
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -156.73253258
  PAW double counting   =    574375.30000613  -573884.24043002
  entropy T*S    EENTRO =         0.00851315
  eigenvalues    EBANDS =     -1928.30431481
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -741.93953772 eV

  energy without entropy =     -741.94805087  energy(sigma->0) =     -741.94237544


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




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



 eigenvalue-minimisations  :  2319
 total energy-change (2. order) :-0.2478651E+07  (-0.3488690E+06)
 number of electron     135.9980837 magnetization      13.1397810
 augmentation part      -12.0301259 magnetization      -3.9175217

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4046
  1.2424  1.2424  0.7968  1.0741  0.4258  0.4258  0.5645  0.4087  0.4087  0.3526
  0.2218  0.2218  0.1807  0.1177  0.0022  0.0006  0.0006  0.0004  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1810.43024016
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -159.12612741
  PAW double counting   =    574981.91436488  -573759.16095445
  entropy T*S    EENTRO =         0.02597351
  eigenvalues    EBANDS =  -2481308.03267694
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -2479392.51563755 eV

  energy without entropy = -2479392.54161107  energy(sigma->0) = -2479392.52429539


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




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



 eigenvalue-minimisations  :  2385
 total energy-change (2. order) : 0.2477119E+07  (-0.1953034E+04)
 number of electron     135.9999981 magnetization      13.0165028
 augmentation part      -10.1963600 magnetization      16.5094632

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3754
  1.2641  1.2641  0.7962  1.0347  0.5652  0.3269  0.3269  0.4026  0.4026  0.3634
  0.2312  0.2312  0.1943  0.0979  0.0045  0.0006  0.0006  0.0004  0.0002  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1810.47163607
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -155.13757511
  PAW double counting   =    574991.87940139  -573742.64647376
  entropy T*S    EENTRO =         0.05619933
  eigenvalues    EBANDS =     -4219.54776375
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -2273.57382497 eV

  energy without entropy =    -2273.63002430  energy(sigma->0) =    -2273.59255808


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




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



 eigenvalue-minimisations  :  2157
 total energy-change (2. order) : 0.2719557E+04  (-0.2663056E+03)
 number of electron     135.9999332 magnetization      12.4571295
 augmentation part       -9.3868664 magnetization      14.7182067

 Broyden mixing:
  rms(total) = 0.41357E+02    rms(broyden)= 0.41356E+02
  rms(prec ) = 0.41494E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3583
  1.2638  1.2638  0.7955  1.0345  0.5652  0.3170  0.3170  0.4028  0.4028  0.3637
  0.2252  0.2252  0.1939  0.1038  0.0471  0.0009  0.0006  0.0006  0.0001  0.0004
  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1811.78585572
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -153.13177333
  PAW double counting   =    578695.09463488  -578145.34724831
  entropy T*S    EENTRO =         0.01064156
  eigenvalues    EBANDS =      -801.15165057
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       445.98277153 eV

  energy without entropy =      445.97212997  energy(sigma->0) =      445.97922434


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




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



 eigenvalue-minimisations  :  2646
 total energy-change (2. order) :-0.1346045E+04  (-0.1128912E+04)
 number of electron     135.9997628 magnetization      10.9562897
 augmentation part       -8.6007360 magnetization      17.1388021

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3454
  1.2639  1.2639  1.0347  0.7957  0.5651  0.4029  0.4029  0.3127  0.3127  0.3637
  0.2212  0.2212  0.1934  0.1074  0.0676  0.0676  0.0010  0.0006  0.0006  0.0001
  0.0004  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1798.72817322
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -163.58535088
  PAW double counting   =    528015.42860290  -527471.43686178
  entropy T*S    EENTRO =         0.01905898
  eigenvalues    EBANDS =     -2144.05372171
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -900.06242270 eV

  energy without entropy =     -900.08148167  energy(sigma->0) =     -900.06877569


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




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



 eigenvalue-minimisations  :  2034
 total energy-change (2. order) : 0.1385433E+04  (-0.2084104E+03)
 number of electron     136.0000659 magnetization      10.4723526
 augmentation part       -8.8333473 magnetization      13.6491804

 Broyden mixing:
  rms(total) = 0.45070E+02    rms(broyden)= 0.45068E+02
  rms(prec ) = 0.45168E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3859
  1.2592  1.2592  0.8120  1.0190  0.6265  0.6265  0.5693  0.3361  0.3361  0.3969
  0.3969  0.3701  0.2246  0.2246  0.2102  0.1030  0.1030  0.0010  0.0006  0.0006
  0.0001  0.0004  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1797.65788989
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -135.34669179
  PAW double counting   =    521571.59469714  -521012.36712867
  entropy T*S    EENTRO =        -0.07252086
  eigenvalues    EBANDS =      -803.07405369
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       485.37043526 eV

  energy without entropy =      485.44295612  energy(sigma->0) =      485.39460888


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




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



 eigenvalue-minimisations  :  2166
 total energy-change (2. order) :-0.2426810E+02  (-0.5478975E+01)
 number of electron     136.0000206 magnetization      11.5722861
 augmentation part       -9.3741230 magnetization      11.8495469

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4270
  1.1902  1.1902  0.8250  1.0744  1.0744  0.9638  0.5821  0.4552  0.4552  0.3949
  0.3949  0.3645  0.2982  0.2982  0.2344  0.2344  0.1071  0.1071  0.0010  0.0006
  0.0006  0.0001  0.0002  0.0004

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1800.45795625
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -139.37691402
  PAW double counting   =    497432.66920592  -496877.21236279
  entropy T*S    EENTRO =        -0.06673469
  eigenvalues    EBANDS =      -816.74692762
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       461.10233356 eV

  energy without entropy =      461.16906825  energy(sigma->0) =      461.12457845


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




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



 eigenvalue-minimisations  :  1749
 total energy-change (2. order) : 0.5160533E+02  (-0.7742863E+03)
 number of electron     136.0000178 magnetization      11.4907224
 augmentation part       -9.4956602 magnetization      12.0706073

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4100
  1.2382  1.2382  0.8250  1.0258  1.0258  0.9539  0.5859  0.4424  0.4424  0.3897
  0.3897  0.3705  0.2972  0.2972  0.2336  0.2336  0.1071  0.1071  0.0439  0.0010
  0.0006  0.0006  0.0001  0.0004  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1801.71380646
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -153.00148806
  PAW double counting   =    450651.71739070  -450024.63151666
  entropy T*S    EENTRO =        -0.06860253
  eigenvalues    EBANDS =      -821.88833822
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       512.70766178 eV

  energy without entropy =      512.77626431  energy(sigma->0) =      512.73052929


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




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



 eigenvalue-minimisations  :  1893
 total energy-change (2. order) :-0.4331529E+02  (-0.1759317E+02)
 number of electron     136.0000020 magnetization      10.3915183
 augmentation part       -9.0013873 magnetization      12.0960836

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4057
  1.3724  1.3724  0.8264  0.9992  0.9992  0.9460  0.5881  0.4460  0.4460  0.3896
  0.3896  0.3728  0.3053  0.3053  0.2144  0.2144  0.1470  0.1021  0.1021  0.0083
  0.0010  0.0006  0.0006  0.0004  0.0002  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1801.56955565
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -151.65678095
  PAW double counting   =    454191.52946056  -453628.96029430
  entropy T*S    EENTRO =        -0.03819199
  eigenvalues    EBANDS =      -802.20628434
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       469.39237633 eV

  energy without entropy =      469.43056833  energy(sigma->0) =      469.40510700


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




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



 eigenvalue-minimisations  :  1929
 total energy-change (2. order) : 0.2429988E+02  (-0.7156877E+01)
 number of electron     135.9999688 magnetization      10.1153150
 augmentation part       -9.0909006 magnetization      13.4430425

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4383
  1.3233  1.3233  1.1522  0.8397  0.8844  0.8844  0.9379  0.5257  0.5257  0.5858
  0.4109  0.4109  0.3186  0.3186  0.3702  0.2657  0.2672  0.2672  0.1066  0.1066
  0.0060  0.0010  0.0006  0.0006  0.0001  0.0004  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1797.89343410
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -127.74068227
  PAW double counting   =    430095.25512857  -429533.50590832
  entropy T*S    EENTRO =        -0.00375477
  eigenvalues    EBANDS =      -804.71311552
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       493.69225661 eV

  energy without entropy =      493.69601138  energy(sigma->0) =      493.69350820


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




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



 eigenvalue-minimisations  :  2196
 total energy-change (2. order) : 0.1003896E+02  (-0.3222995E+01)
 number of electron     136.0000270 magnetization      10.9018543
 augmentation part       -8.9870205 magnetization      15.5522130

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4369
  1.4511  1.4511  1.1808  0.8376  0.8553  0.8553  0.8668  0.5477  0.5477  0.6135
  0.4083  0.4083  0.3217  0.3217  0.3756  0.2702  0.2702  0.2520  0.1762  0.1067
  0.1067  0.0060  0.0010  0.0006  0.0006  0.0001  0.0002  0.0004

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1796.74541478
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -122.22072452
  PAW double counting   =    428094.90638640  -427528.23862110
  entropy T*S    EENTRO =        -0.01046561
  eigenvalues    EBANDS =      -806.25397027
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       503.73121313 eV

  energy without entropy =      503.74167874  energy(sigma->0) =      503.73470167


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




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



 eigenvalue-minimisations  :  2133
 total energy-change (2. order) :-0.1947381E+06  (-0.4374004E+05)
 number of electron     136.0001381 magnetization      10.9090951
 augmentation part       -9.4429933 magnetization      10.0162260

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4218
  1.4616  1.4616  1.1588  0.8380  0.8529  0.8529  0.8614  0.5479  0.5479  0.6162
  0.4063  0.4063  0.3221  0.3221  0.3778  0.2707  0.2707  0.2480  0.1878  0.1067
  0.1067  0.0060  0.0010  0.0006  0.0006  0.0001  0.0002  0.0004  0.0003

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1790.72116475
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -140.21423355
  PAW double counting   =    396962.46359179  -395802.32942166
  entropy T*S    EENTRO =         0.01100730
  eigenvalues    EBANDS =   -196125.84754983
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =   -194234.34374769 eV

  energy without entropy =  -194234.35475500  energy(sigma->0) =  -194234.34741679


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




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



 eigenvalue-minimisations  :  2415
 total energy-change (2. order) : 0.1903873E+06  (-0.9989177E+04)
 number of electron     135.9997266 magnetization       9.1721887
 augmentation part       -9.7389068 magnetization      19.5461101

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4077
  1.4617  1.4617  1.1561  0.8380  0.8530  0.8530  0.8614  0.5477  0.5477  0.6162
  0.4061  0.4061  0.3221  0.3221  0.3779  0.2706  0.2706  0.2477  0.1881  0.1067
  0.1067  0.0060  0.0011  0.0010  0.0006  0.0006  0.0004  0.0001  0.0002  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1790.72611008
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -133.22921970
  PAW double counting   =    395671.09401717  -395416.37403553
  entropy T*S    EENTRO =         0.00872621
  eigenvalues    EBANDS =     -4840.14622494
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -3847.07882386 eV

  energy without entropy =    -3847.08755007  energy(sigma->0) =    -3847.08173259


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




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



 eigenvalue-minimisations  :  2613
 total energy-change (2. order) :-0.1126751E+07  (-0.2564385E+05)
 number of electron     136.2906045 magnetization       9.1521222
 augmentation part      -10.6343171 magnetization      -4.7852273

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3960
  1.4667  1.4667  1.1554  0.8380  0.8515  0.8515  0.8601  0.5480  0.5480  0.6174
  0.4061  0.4061  0.3221  0.3221  0.3779  0.2705  0.2705  0.2483  0.1856  0.1067
  0.1067  0.0200  0.0224  0.0060  0.0010  0.0006  0.0006  0.0001  0.0004  0.0002
  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1770.20744096
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -127.15585469
  PAW double counting   =    368050.15921397  -367574.19941299
  entropy T*S    EENTRO =        -0.01550806
  eigenvalues    EBANDS =  -1131839.00873132
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1130598.13371102 eV

  energy without entropy = -1130598.11820297  energy(sigma->0) = -1130598.12854167


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




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



 eigenvalue-minimisations  :  2124
 total energy-change (2. order) : 0.1130205E+07  (-0.1194415E+04)
 number of electron     141.5105504 magnetization       7.4777699
 augmentation part       -8.4558862 magnetization      15.0001114

 Broyden mixing:
  rms(total) = 0.54398E+02    rms(broyden)= 0.54392E+02
  rms(prec ) = 0.54462E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3910
  1.4907  1.4907  1.1423  0.8380  0.8390  0.8390  0.8737  0.5482  0.5482  0.6064
  0.4110  0.4110  0.3210  0.3210  0.3758  0.2681  0.2681  0.2574  0.1366  0.1366
  0.1636  0.1067  0.1067  0.0060  0.0013  0.0010  0.0006  0.0006  0.0001  0.0004
  0.0002  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1770.36590350
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -107.50783868
  PAW double counting   =    367690.39683023  -367129.89078804
  entropy T*S    EENTRO =        -0.00560903
  eigenvalues    EBANDS =     -1738.45351947
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -393.53280547 eV

  energy without entropy =     -393.52719644  energy(sigma->0) =     -393.53093580


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




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



 eigenvalue-minimisations  :  2085
 total energy-change (2. order) : 0.8644301E+03  (-0.4640491E+02)
 number of electron     139.5565958 magnetization       3.7414583
 augmentation part       -8.4820950 magnetization      10.2171811

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3972
  1.5071  1.5071  1.1381  0.8375  0.8131  0.8131  0.8950  0.5507  0.5507  0.6041
  0.3277  0.3277  0.4033  0.4033  0.3797  0.3220  0.3220  0.2807  0.2807  0.2419
  0.2419  0.1381  0.1067  0.1067  0.0060  0.0013  0.0010  0.0006  0.0006  0.0001
  0.0004  0.0002  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1782.78174945
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =      -123.23116483
  PAW double counting   =    442367.32195364  -441803.66996917
  entropy T*S    EENTRO =        -0.00587824
  eigenvalues    EBANDS =      -849.02992040
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       470.89729456 eV

  energy without entropy =      470.90317280  energy(sigma->0) =      470.89925398


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




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



 eigenvalue-minimisations  :  2511
 total energy-change (2. order) : 0.4874721E+02  (-0.1890946E+02)
 number of electron     137.7634985 magnetization       3.7367149
 augmentation part       -9.2019420 magnetization      11.9743522

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3857
  1.5072  1.5072  1.1280  0.8377  0.8141  0.8141  0.8935  0.5494  0.5494  0.6041
  0.3206  0.3206  0.4026  0.4026  0.3800  0.3219  0.3219  0.2810  0.2810  0.2435
  0.2435  0.1382  0.1067  0.1067  0.0284  0.0060  0.0013  0.0010  0.0006  0.0006
  0.0001  0.0004  0.0002  0.0001

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1789.94328336
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       -43.65207213
  PAW double counting   =    481009.62482725  -480457.54996198
  entropy T*S    EENTRO =        -0.01792212
  eigenvalues    EBANDS =      -861.11110503
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       519.64450565 eV

  energy without entropy =      519.66242777  energy(sigma->0) =      519.65047969


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




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



 eigenvalue-minimisations  :  2124
 total energy-change (2. order) :-0.1194776E+05  (-0.5068365E+04)
 number of electron     159.3711302 magnetization       3.7303533
 augmentation part      -15.9028927 magnetization      -2.9296827

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3738
  1.5087  1.5087  1.1170  0.8380  0.8140  0.8140  0.8931  0.5482  0.5482  0.6032
  0.3203  0.3203  0.4017  0.4017  0.3808  0.3218  0.3218  0.2796  0.2796  0.2420
  0.2420  0.1365  0.1067  0.1067  0.0170  0.0060  0.0013  0.0010  0.0006  0.0006
  0.0004  0.0002  0.0001  0.0001  0.0002

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1789.89414920
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       -41.44331281
  PAW double counting   =    481312.70552506  -480814.40772313
  entropy T*S    EENTRO =        -0.02385901
  eigenvalues    EBANDS =    -12757.34763576
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -11428.11713184 eV

  energy without entropy =   -11428.09327284  energy(sigma->0) =   -11428.10917884


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




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



 eigenvalue-minimisations  :  2112
 total energy-change (2. order) : 0.1074213E+05  (-0.1219100E+05)
 number of electron     153.5616622 magnetization       2.5242116
 augmentation part      -10.0307661 magnetization      18.3901514

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3633
  1.5088  1.5088  1.1177  0.8380  0.8139  0.8139  0.8931  0.5483  0.5483  0.6031
  0.3208  0.3208  0.4017  0.4017  0.3808  0.3218  0.3218  0.2794  0.2794  0.2416
  0.2416  0.1363  0.1067  0.1067  0.0121  0.0060  0.0013  0.0010  0.0006  0.0006
  0.0004  0.0004  0.0001  0.0002  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1791.43183684
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       -18.34958100
  PAW double counting   =    481770.70732176  -481142.98537141
  entropy T*S    EENTRO =        -0.02148422
  eigenvalues    EBANDS =     -2166.19786920
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -685.98479790 eV

  energy without entropy =     -685.96331367  energy(sigma->0) =     -685.97763649


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




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



 eigenvalue-minimisations  :  2292
 total energy-change (2. order) : 0.9430145E+03  (-0.3702112E+03)
 number of electron     148.7332916 magnetization       3.1730605
 augmentation part       -9.7961762 magnetization       5.8445546

 Broyden mixing:
  rms(total) = 0.65765E+02    rms(broyden)= 0.65762E+02
  rms(prec ) = 0.65915E+02
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3646
  1.5090  1.5090  1.1278  0.8376  0.8179  0.8179  0.8911  0.5461  0.5461  0.6040
  0.2914  0.2914  0.4019  0.4019  0.3812  0.3212  0.3212  0.2832  0.2832  0.2516
  0.2516  0.2206  0.2206  0.1405  0.1067  0.1067  0.0060  0.0013  0.0010  0.0006
  0.0006  0.0005  0.0004  0.0002  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1801.56614382
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        -5.72163058
  PAW double counting   =    531332.98527447  -530707.76189444
  entropy T*S    EENTRO =        -0.02864017
  eigenvalues    EBANDS =     -1223.17129287
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       257.02969560 eV

  energy without entropy =      257.05833577  energy(sigma->0) =      257.03924232


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




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



 eigenvalue-minimisations  :  2307
 total energy-change (2. order) :-0.1594520E+11  (-0.2940228E+10)
 number of electron     147.6144892 magnetization       3.1907799
 augmentation part       -8.5929976 magnetization      -1.1479338

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3552
  1.5097  1.5097  1.1261  0.8376  0.8910  0.8181  0.8181  0.6041  0.5460  0.5460
  0.4018  0.4018  0.2922  0.2922  0.3812  0.3212  0.3212  0.2832  0.2832  0.2518
  0.2518  0.2200  0.2200  0.1405  0.1067  0.1067  0.0060  0.0033  0.0013  0.0010
  0.0006  0.0006  0.0005  0.0004  0.0002  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1783.07165065
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       -58.45172100
  PAW double counting   =    577984.43027944  -577936.48789266
  entropy T*S    EENTRO =        -0.00245495
  eigenvalues    EBANDS = ******************
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  = ****************** eV

  energy without entropy =******************  energy(sigma->0) =******************


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




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



 eigenvalue-minimisations  :  2160
 total energy-change (2. order) : 0.1594385E+11  (-0.9966076E+06)
 number of electron     141.4278720 magnetization      -4.9759471
 augmentation part      -10.1265492 magnetization       7.4023099

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3496
  1.5099  1.5099  1.1270  0.8376  0.8186  0.8186  0.8913  0.5463  0.5463  0.6042
  0.3126  0.3126  0.4022  0.4022  0.3810  0.3210  0.3210  0.2840  0.2840  0.2546
  0.2546  0.2222  0.2222  0.0814  0.1412  0.1067  0.1067  0.0060  0.0021  0.0013
  0.0010  0.0006  0.0006  0.0005  0.0004  0.0002  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1784.50141139
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       -33.64384336
  PAW double counting   =    593806.88317085  -593291.64364410
  entropy T*S    EENTRO =        -0.03792425
  eigenvalues    EBANDS =  -1350436.11251897
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =  -1349076.76214820 eV

  energy without entropy = -1349076.72422394  energy(sigma->0) = -1349076.74950678


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




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



 eigenvalue-minimisations  :  1938
 total energy-change (2. order) : 0.1338549E+07  (-0.9720534E+04)
 number of electron     158.3283538 magnetization      -5.0833436
 augmentation part      -10.1842329 magnetization      23.1262313

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3421
  1.5099  1.5099  1.1254  0.8376  0.8187  0.8187  0.8903  0.5458  0.5458  0.6041
  0.3011  0.3011  0.4021  0.4021  0.3816  0.3209  0.3209  0.2837  0.2837  0.2548
  0.2548  0.2178  0.2178  0.1412  0.1067  0.1067  0.0845  0.0845  0.0060  0.0022
  0.0013  0.0010  0.0006  0.0006  0.0004  0.0005  0.0002  0.0001  0.0001  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1799.97987534
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       184.89220214
  PAW double counting   =   1019356.46431737 -1017602.91828758
  entropy T*S    EENTRO =         0.01178327
  eigenvalues    EBANDS =    -13328.38803680
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =    -10527.62387393 eV

  energy without entropy =   -10527.63565720  energy(sigma->0) =   -10527.62780168


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




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



 eigenvalue-minimisations  :  2778
 total energy-change (2. order) : 0.3048031E+04  (-0.8083870E+04)
 number of electron     149.0061959 magnetization      -5.2875461
 augmentation part       -8.9547359 magnetization       8.6740061

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

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3349
  1.5101  1.5101  1.1248  0.8377  0.8893  0.8147  0.8147  0.6030  0.5453  0.5453
  0.4012  0.4012  0.3830  0.3210  0.3210  0.2638  0.2638  0.2829  0.2829  0.2541
  0.2541  0.2010  0.2010  0.1673  0.1673  0.1412  0.1067  0.1067  0.0060  0.0031
  0.0022  0.0013  0.0010  0.0006  0.0006  0.0005  0.0004  0.0002  0.0001  0.0001
  0.0000

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =        24.76344702
  Ewald energy   TEWEN  =     -1117.79376468
  -Hartree energ DENC   =     -1800.52955310
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =       136.40119724
  PAW double counting   =   1024004.85061078 -1023468.06684744
  entropy T*S    EENTRO =         0.07400194
  eigenvalues    EBANDS =     -9014.61581490
  atomic energy  EATOM  =      3755.32434067
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =     -7479.59238247 eV

  energy without entropy =    -7479.66638440  energy(sigma->0) =    -7479.61704978


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




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


 -----------------------------------------------------------------------------
|                                                                             |
|     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 = 8 1 9                  |
|                                                                             |
|       ---->  I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <----       |
|                                                                             |
 -----------------------------------------------------------------------------