./iterations/neb1_max2_image01_iter91_OUTCAR.out output for 41: 3.3A-->1.3A (Job 121) O-H...N hydrogen bond by TSA-OH, TSA-NH2 (TSS)

Status: finished
 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.26  23:18:15
 running on    1 total cores
 distrk:  each k-point on    1 cores,    1 groups
 distr:  one band on NCORE=   1 cores,    1 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 = 1
   INIWAV = 1
   ISTART = 0
   ICHARG = 2
   LWAVE = .FALSE.
   LCHARG = .FALSE.
   ADDGRID = .FALSE.
   ISMEAR = 1
   SIGMA = 0.2
   LREAL = Auto
   LSCALAPACK = .FALSE.
   RWIGS = 0.75 1.11 0.73 0.32

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

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

   l    n(q)    QCUT    max X(q) W(low)/X(q) W(high)/X(q)  e(spline) 
   0     10    10.053    79.467    0.76E-04    0.72E-04    0.56E-06
   0     10    10.053    66.151    0.76E-04    0.72E-04    0.55E-06
   1     10    10.053     8.350    0.25E-03    0.92E-03    0.41E-05
   1     10    10.053     5.531    0.27E-03    0.10E-02    0.45E-05
 Optimization of the real space projectors (new method)

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

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

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

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

 maximal supplied QI-value         = 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
  PAW_PBE N 08Apr2002                   :
 energy of atom  1       EATOM= -264.5486
 kinetic energy error for atom=    0.0736 (will be added to EATOM!!)
  PAW_PBE Si 05Jan2001                  :
 energy of atom  2       EATOM= -103.0669
 kinetic energy error for atom=    0.0012 (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!!)
  PAW_PBE H 15Jun2001                   :
 energy of atom  4       EATOM=  -12.4884
 kinetic energy error for atom=    0.0098 (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.077  0.981  0.642-
   2  0.960  0.656  0.781-
   3  0.475  0.415  0.435-
   4  0.087  0.314  0.771-
   5  0.253  0.061  0.491-
   6  0.435  0.120  0.836-
   7  0.190  0.011  0.817-
   8  0.796  0.117  0.585-  20 1.25
   9  0.028  0.760  0.331-
  10  0.146  0.441  0.390-
  11  0.092  0.625  0.305-
  12  0.508  0.656  0.176-
  13  0.342  0.135  0.398-
  14  0.843  0.896  0.813-
  15  0.114  0.805  0.059-
  16  0.275  0.744  0.339-
  17  0.299  0.551  0.737-
  18  0.673  0.758  0.833-
  19  0.997  0.317  0.536-
  20  0.830  0.198  0.635-   8 1.25
  21  0.668  0.315  0.465-
  22  0.783  0.183  0.370-
  23  0.030  0.907  0.752-
  24  0.819  0.352  0.779-
  25  0.886  0.781  0.689-
  26  0.496  0.814  0.740-
  27  0.457  0.660  0.737-
  28  0.852  0.010  0.350-
  29  0.555  0.206  0.314-
 
  LATTYP: Found a simple tetragonal cell.
 ALAT       =    12.0000000000
 C/A-ratio  =     1.2500000000
  
  Lattice vectors:
  
 A1 = (   0.0000000000,  12.0000000000,   0.0000000000)
 A2 = (   0.0000000000,   0.0000000000,  12.0000000000)
 A3 = (  15.0000000000,   0.0000000000,   0.0000000000)


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

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


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


The static configuration has the point symmetry C_1 .


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

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


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


The dynamic configuration has the point symmetry C_1 .


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


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

                                     Primitive cell                                     

  volume of cell :    2160.0000

  direct lattice vectors                    reciprocal lattice vectors
    15.000000000  0.000000000  0.000000000     0.066666667  0.000000000  0.000000000
     0.000000000 12.000000000  0.000000000     0.000000000  0.083333333  0.000000000
     0.000000000  0.000000000 12.000000000     0.000000000  0.000000000  0.083333333

  length of vectors
    15.000000000 12.000000000 12.000000000     0.066666667  0.083333333  0.083333333

  position of ions in fractional coordinates (direct lattice)
     0.077384000  0.981015660  0.641587860
     0.960488550  0.655947870  0.781445430
     0.474935080  0.415230110  0.435127510
     0.086712430  0.313739310  0.770961970
     0.252675700  0.060902090  0.491479830
     0.434568490  0.120253990  0.836428790
     0.190196280  0.011340450  0.816869980
     0.796250100  0.116809790  0.585341970
     0.027537300  0.760438450  0.330808530
     0.145916520  0.441359460  0.389845980
     0.091706120  0.625093730  0.305094510
     0.507848560  0.655567590  0.176232010
     0.342370380  0.134785040  0.398459760
     0.842900500  0.895968920  0.812767320
     0.114163000  0.805374750  0.058811960
     0.274716210  0.744169120  0.338712200
     0.298771170  0.551202930  0.737279580
     0.673072430  0.757869870  0.833189050
     0.997337580  0.316870180  0.536177160
     0.829590130  0.198357380  0.634930560
     0.667585380  0.315329120  0.464784290
     0.782711250  0.183465170  0.369778010
     0.030108540  0.906680730  0.751905750
     0.818968720  0.351918540  0.778857480
     0.886183420  0.780666600  0.688619350
     0.496228220  0.814036410  0.739556270
     0.457380360  0.659840360  0.736583570
     0.851851740  0.010383050  0.350294310
     0.555127330  0.206146210  0.314481020

  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

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

 
 
 KPOINTS: Automatic mesh                          

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

 Generating k-lattice:

  Cartesian coordinates                     Fractional coordinates (reciprocal lattice)
     0.066666667  0.000000000  0.000000000     1.000000000  0.000000000  0.000000000
     0.000000000  0.041666667  0.000000000     0.000000000  0.500000000  0.000000000
     0.000000000  0.000000000  0.041666667     0.000000000  0.000000000  0.500000000

  Length of vectors
     0.066666667  0.041666667  0.041666667

  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.000000  0.500000  0.000000      1.000000
  0.000000  0.000000  0.500000      1.000000
  0.000000  0.500000  0.500000      1.000000
 
 Following cartesian coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000
  0.000000  0.041667  0.000000      1.000000
  0.000000  0.000000  0.041667      1.000000
  0.000000  0.041667  0.041667      1.000000
 


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




 Dimension of arrays:
   k-points           NKPTS =      4   k-points in BZ     NKDIM =      4   number of bands    NBANDS=     47
   number of dos      NEDOS =    301   number of ions     NIONS =     29
   non local maximal  LDIM  =      4   non local SUM 2l+1 LMDIM =      8
   total plane-waves  NPLWV = 288000
   max r-space proj   IRMAX =   2551   max aug-charges    IRDMAX=   4657
   dimension x,y,z NGX =    80 NGY =   60 NGZ =   60
   dimension x,y,z NGXF=   160 NGYF=  120 NGZF=  120
   support grid    NGXF=   160 NGYF=  120 NGZF=  120
   ions per type =               3   6   1  19
   NGX,Y,Z   is equivalent  to a cutoff of   8.87,  8.31,  8.31 a.u.
   NGXF,Y,Z  is equivalent  to a cutoff of  17.73, 16.62, 16.62 a.u.

 SYSTEM =  No title                                
 POSCAR =  No title                                

 Startparameter for this run:
   NWRITE =      1    write-flag & timer
   PREC   = normal    normal or accurate (medium, high low for compatibility)
   ISTART =      0    job   : 0-new  1-cont  2-samecut
   ICHARG =      2    charge: 1-file 2-atom 10-const
   ISPIN  =      1    spin polarized calculation?
   LNONCOLLINEAR =      F non collinear calculations
   LSORBIT =      F    spin-orbit coupling
   INIWAV =      1    electr: 0-lowe 1-rand  2-diag
   LASPH  =      F    aspherical Exc in radial PAW
 Electronic Relaxation 1
   ENCUT  =  400.0 eV  29.40 Ry    5.42 a.u.  24.46 19.57 19.57*2*pi/ulx,y,z
   ENINI  =  400.0     initial cutoff
   ENAUG  =  627.1 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  -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.514E-26a.u.
   SCALEE = 1.0000    scale energy and forces
   NPACO  =    256;   APACO  = 16.0  distance and # of slots for P.C.
   PSTRESS=    0.0 pullay stress

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

   TIME   =   0.40     timestep for ELM

  volume/ion in A,a.u.               =      74.48       502.63
  Fermi-wavevector in a.u.,A,eV,Ry     =   0.506583  0.957302  3.491607  0.256626
  Thomas-Fermi vector in A             =   1.517676
 
 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
 non-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           15
 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 :     2160.00
      direct lattice vectors                 reciprocal lattice vectors
    15.000000000  0.000000000  0.000000000     0.066666667  0.000000000  0.000000000
     0.000000000 12.000000000  0.000000000     0.000000000  0.083333333  0.000000000
     0.000000000  0.000000000 12.000000000     0.000000000  0.000000000  0.083333333

  length of vectors
    15.000000000 12.000000000 12.000000000     0.066666667  0.083333333  0.083333333


 
 k-points in units of 2pi/SCALE and weight: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       0.250
   0.00000000  0.04166667  0.00000000       0.250
   0.00000000  0.00000000  0.04166667       0.250
   0.00000000  0.04166667  0.04166667       0.250
 
 k-points in reciprocal lattice and weights: Automatic mesh                          
   0.00000000  0.00000000  0.00000000       0.250
   0.00000000  0.50000000  0.00000000       0.250
   0.00000000  0.00000000  0.50000000       0.250
   0.00000000  0.50000000  0.50000000       0.250
 
 position of ions in fractional coordinates (direct lattice) 
   0.07738400  0.98101566  0.64158786
   0.96048855  0.65594787  0.78144543
   0.47493508  0.41523011  0.43512751
   0.08671243  0.31373931  0.77096197
   0.25267570  0.06090209  0.49147983
   0.43456849  0.12025399  0.83642879
   0.19019628  0.01134045  0.81686998
   0.79625010  0.11680979  0.58534197
   0.02753730  0.76043845  0.33080853
   0.14591652  0.44135946  0.38984598
   0.09170612  0.62509373  0.30509451
   0.50784856  0.65556759  0.17623201
   0.34237038  0.13478504  0.39845976
   0.84290050  0.89596892  0.81276732
   0.11416300  0.80537475  0.05881196
   0.27471621  0.74416912  0.33871220
   0.29877117  0.55120293  0.73727958
   0.67307243  0.75786987  0.83318905
   0.99733758  0.31687018  0.53617716
   0.82959013  0.19835738  0.63493056
   0.66758538  0.31532912  0.46478429
   0.78271125  0.18346517  0.36977801
   0.03010854  0.90668073  0.75190575
   0.81896872  0.35191854  0.77885748
   0.88618342  0.78066660  0.68861935
   0.49622822  0.81403641  0.73955627
   0.45738036  0.65984036  0.73658357
   0.85185174  0.01038305  0.35029431
   0.55512733  0.20614621  0.31448102
 
 position of ions in cartesian coordinates  (Angst):
   1.16076000 11.77218792  7.69905432
  14.40732825  7.87137444  9.37734516
   7.12402620  4.98276132  5.22153012
   1.30068645  3.76487172  9.25154364
   3.79013550  0.73082508  5.89775796
   6.51852735  1.44304788 10.03714548
   2.85294420  0.13608540  9.80243976
  11.94375150  1.40171748  7.02410364
   0.41305950  9.12526140  3.96970236
   2.18874780  5.29631352  4.67815176
   1.37559180  7.50112476  3.66113412
   7.61772840  7.86681108  2.11478412
   5.13555570  1.61742048  4.78151712
  12.64350750 10.75162704  9.75320784
   1.71244500  9.66449700  0.70574352
   4.12074315  8.93002944  4.06454640
   4.48156755  6.61443516  8.84735496
  10.09608645  9.09443844  9.99826860
  14.96006370  3.80244216  6.43412592
  12.44385195  2.38028856  7.61916672
  10.01378070  3.78394944  5.57741148
  11.74066875  2.20158204  4.43733612
   0.45162810 10.88016876  9.02286900
  12.28453080  4.22302248  9.34628976
  13.29275130  9.36799920  8.26343220
   7.44342330  9.76843692  8.87467524
   6.86070540  7.91808432  8.83900284
  12.77777610  0.12459660  4.20353172
   8.32690995  2.47375452  3.77377224
 


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


 k-point  1 :   0.0000 0.0000 0.0000  plane waves:   39297
 k-point  2 :   0.0000 0.5000 0.0000  plane waves:   39224
 k-point  3 :   0.0000 0.0000 0.5000  plane waves:   39224
 k-point  4 :   0.0000 0.5000 0.5000  plane waves:   39160

 maximum and minimum number of plane-waves per node :     39297    39160

 maximum number of plane-waves:     39297
 maximum index in each direction: 
   IXMAX=   24   IYMAX=   19   IZMAX=   19
   IXMIN=  -24   IYMIN=  -20   IZMIN=  -20

 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    98 to avoid them
 WARNING: aliasing errors must be expected set NGY to    80 to avoid them
 WARNING: aliasing errors must be expected set NGZ to    80 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   411352. kBytes
=======================================================================

   base      :      30000. kBytes
   nonlr-proj:       3185. kBytes
   fftplans  :      79712. kBytes
   grid      :     179635. kBytes
   one-center:         89. kBytes
   wavefun   :     118731. kBytes
 
 Broyden mixing: mesh for mixing (old mesh)
   NGX = 49   NGY = 39   NGZ = 39
  (NGX  =160   NGY  =120   NGZ  =120)
  gives a total of  74529 points

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


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


 Maximum index for non-local projection operator         2428
 Maximum index for augmentation-charges         4333 (set IRDMAX)


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


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



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



 eigenvalue-minimisations  :   376
 total energy-change (2. order) : 0.2834000E+03  (-0.1424555E+04)
 number of electron      64.0000000 magnetization 
 augmentation part       64.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1265.48989666
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        51.51743548
  PAW double counting   =      1530.63699070    -1430.61760059
  entropy T*S    EENTRO =         0.00679079
  eigenvalues    EBANDS =      -229.05340311
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       283.39997027 eV

  energy without entropy =      283.39317948  energy(sigma->0) =      283.39770667


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




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



 eigenvalue-minimisations  :   528
 total energy-change (2. order) :-0.2527173E+03  (-0.2445830E+03)
 number of electron      64.0000000 magnetization 
 augmentation part       64.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1265.48989666
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        51.51743548
  PAW double counting   =      1530.63699070    -1430.61760059
  entropy T*S    EENTRO =        -0.00992067
  eigenvalues    EBANDS =      -481.75400401
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =        30.68265791 eV

  energy without entropy =       30.69257858  energy(sigma->0) =       30.68596480


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




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



 eigenvalue-minimisations  :   540
 total energy-change (2. order) :-0.7483470E+02  (-0.6351863E+02)
 number of electron      64.0000000 magnetization 
 augmentation part       64.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1265.48989666
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        51.51743548
  PAW double counting   =      1530.63699070    -1430.61760059
  entropy T*S    EENTRO =         0.04872886
  eigenvalues    EBANDS =      -556.64735337
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -44.15204192 eV

  energy without entropy =      -44.20077078  energy(sigma->0) =      -44.16828487


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




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



 eigenvalue-minimisations  :   568
 total energy-change (2. order) :-0.1324907E+02  (-0.9949800E+01)
 number of electron      64.0000000 magnetization 
 augmentation part       64.0000000 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1265.48989666
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        51.51743548
  PAW double counting   =      1530.63699070    -1430.61760059
  entropy T*S    EENTRO =        -0.17214864
  eigenvalues    EBANDS =      -569.67554740
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -57.40111344 eV

  energy without entropy =      -57.22896481  energy(sigma->0) =      -57.34373056


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




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



 eigenvalue-minimisations  :   568
 total energy-change (2. order) :-0.1509849E+01  (-0.1426861E+01)
 number of electron      64.0000075 magnetization 
 augmentation part        0.7940245 magnetization 

 Broyden mixing:
  rms(total) = 0.26209E+01    rms(broyden)= 0.26196E+01
  rms(prec ) = 0.37316E+01
  weight for this iteration     100.00

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1265.48989666
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        51.51743548
  PAW double counting   =      1530.63699070    -1430.61760059
  entropy T*S    EENTRO =        -0.15455359
  eigenvalues    EBANDS =      -571.20299137
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -58.91096237 eV

  energy without entropy =      -58.75640878  energy(sigma->0) =      -58.85944450


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




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



 eigenvalue-minimisations  :   656
 total energy-change (2. order) :-0.5809294E+02  (-0.6701566E+02)
 number of electron      63.9999981 magnetization 
 augmentation part       -2.1246998 magnetization 

 Broyden mixing:
  rms(total) = 0.47748E+01    rms(broyden)= 0.47722E+01
  rms(prec ) = 0.64845E+01
  weight for this iteration     100.00

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1411.36087290
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        59.67840482
  PAW double counting   =      1866.45721903    -1767.20835446
  entropy T*S    EENTRO =        -0.00305978
  eigenvalues    EBANDS =      -490.96689028
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =      -117.00389991 eV

  energy without entropy =     -117.00084013  energy(sigma->0) =     -117.00287998


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




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



 eigenvalue-minimisations  :   536
 total energy-change (2. order) : 0.6243879E+02  (-0.6480918E+01)
 number of electron      64.0000055 magnetization 
 augmentation part       -1.0472173 magnetization 

 Broyden mixing:
  rms(total) = 0.28177E+01    rms(broyden)= 0.28168E+01
  rms(prec ) = 0.38524E+01
  weight for this iteration     100.00

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

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1311.30346401
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        55.03100894
  PAW double counting   =      1723.48914680    -1623.52618971
  entropy T*S    EENTRO =         0.11375883
  eigenvalues    EBANDS =      -524.76902276
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -54.56510824 eV

  energy without entropy =      -54.67886708  energy(sigma->0) =      -54.60302786


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




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



 eigenvalue-minimisations  :   576
 total energy-change (2. order) : 0.7483642E+01  (-0.3074578E+01)
 number of electron      64.0000055 magnetization 
 augmentation part       -0.8535999 magnetization 

 Broyden mixing:
  rms(total) = 0.24704E+01    rms(broyden)= 0.24700E+01
  rms(prec ) = 0.34031E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.2961
  0.3618  0.2632  0.2632

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1306.77548267
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        54.76942183
  PAW double counting   =      1721.19165276    -1621.19064570
  entropy T*S    EENTRO =         0.01504808
  eigenvalues    EBANDS =      -521.49111424
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -47.08146629 eV

  energy without entropy =      -47.09651436  energy(sigma->0) =      -47.08648231


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




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



 eigenvalue-minimisations  :   656
 total energy-change (2. order) : 0.3193344E+01  (-0.5522635E+00)
 number of electron      64.0000044 magnetization 
 augmentation part       -1.1419808 magnetization 

 Broyden mixing:
  rms(total) = 0.23491E+01    rms(broyden)= 0.23489E+01
  rms(prec ) = 0.32239E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3826
  0.4500  0.4500  0.3152  0.3152

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1309.23130007
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        54.85640211
  PAW double counting   =      1700.10008479    -1600.13928666
  entropy T*S    EENTRO =        -0.02538918
  eigenvalues    EBANDS =      -515.84828670
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -43.88812206 eV

  energy without entropy =      -43.86273288  energy(sigma->0) =      -43.87965900


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




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



 eigenvalue-minimisations  :   681
 total energy-change (2. order) :-0.1706523E+01  (-0.9630954E+01)
 number of electron      63.9999988 magnetization 
 augmentation part       -0.9511756 magnetization 

 Broyden mixing:
  rms(total) = 0.31125E+01    rms(broyden)= 0.31098E+01
  rms(prec ) = 0.42116E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3776
  0.4916  0.4916  0.3418  0.3418  0.2212

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1297.50138640
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        54.19920480
  PAW double counting   =      1622.93972916    -1522.90569181
  entropy T*S    EENTRO =         0.00910572
  eigenvalues    EBANDS =      -528.73526045
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -45.59464532 eV

  energy without entropy =      -45.60375104  energy(sigma->0) =      -45.59768056


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




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



 eigenvalue-minimisations  :   608
 total energy-change (2. order) :-0.2240069E+01  (-0.2172228E+01)
 number of electron      63.9999995 magnetization 
 augmentation part       -0.8239209 magnetization 

 Broyden mixing:
  rms(total) = 0.34458E+01    rms(broyden)= 0.34453E+01
  rms(prec ) = 0.47780E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3806
  0.6031  0.6031  0.3324  0.3324  0.2062  0.2062

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1286.60323673
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.62076849
  PAW double counting   =      1599.60062402    -1499.40948542
  entropy T*S    EENTRO =         0.04472592
  eigenvalues    EBANDS =      -541.48776428
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -47.83471434 eV

  energy without entropy =      -47.87944026  energy(sigma->0) =      -47.84962298


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




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



 eigenvalue-minimisations  :   656
 total energy-change (2. order) : 0.3549776E+01  (-0.4455013E+01)
 number of electron      64.0000062 magnetization 
 augmentation part        0.3131217 magnetization 

 Broyden mixing:
  rms(total) = 0.23132E+01    rms(broyden)= 0.23117E+01
  rms(prec ) = 0.32662E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3665
  0.8777  0.3334  0.3334  0.3334  0.3334  0.1770  0.1770

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1281.52605823
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.18451643
  PAW double counting   =      1640.20165739    -1539.96920460
  entropy T*S    EENTRO =        -0.01206283
  eigenvalues    EBANDS =      -542.56343985
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -44.28493804 eV

  energy without entropy =      -44.27287521  energy(sigma->0) =      -44.28091709


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




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



 eigenvalue-minimisations  :   600
 total energy-change (2. order) : 0.3327604E+01  (-0.1430306E+01)
 number of electron      64.0000055 magnetization 
 augmentation part        0.1832186 magnetization 

 Broyden mixing:
  rms(total) = 0.18458E+01    rms(broyden)= 0.18453E+01
  rms(prec ) = 0.26237E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3809
  1.0880  0.3467  0.3467  0.3314  0.3314  0.2150  0.2150  0.1730

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1289.53970117
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.50006576
  PAW double counting   =      1668.91200884    -1568.75129157
  entropy T*S    EENTRO =        -0.11766762
  eigenvalues    EBANDS =      -531.36040242
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -40.95733453 eV

  energy without entropy =      -40.83966691  energy(sigma->0) =      -40.91811199


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




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



 eigenvalue-minimisations  :   624
 total energy-change (2. order) : 0.5865524E+00  (-0.7030772E+00)
 number of electron      63.9999985 magnetization 
 augmentation part       -1.0416248 magnetization 

 Broyden mixing:
  rms(total) = 0.28507E+01    rms(broyden)= 0.28492E+01
  rms(prec ) = 0.38442E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3627
  1.1874  0.3617  0.3617  0.3362  0.3362  0.2637  0.1460  0.1460  0.1253

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1297.36730525
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.84661185
  PAW double counting   =      1690.08298626    -1589.98023696
  entropy T*S    EENTRO =         0.07391015
  eigenvalues    EBANDS =      -523.42640183
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -40.37078213 eV

  energy without entropy =      -40.44469228  energy(sigma->0) =      -40.39541885


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




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



 eigenvalue-minimisations  :   658
 total energy-change (2. order) : 0.8155649E+00  (-0.3954370E+00)
 number of electron      64.0000037 magnetization 
 augmentation part       -0.9274258 magnetization 

 Broyden mixing:
  rms(total) = 0.14024E+01    rms(broyden)= 0.13979E+01
  rms(prec ) = 0.19577E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3615
  1.2819  0.3570  0.3570  0.3537  0.3537  0.3404  0.1951  0.1569  0.1569  0.0620

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1298.36545714
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.85579380
  PAW double counting   =      1707.75303445    -1607.63876566
  entropy T*S    EENTRO =        -0.14468859
  eigenvalues    EBANDS =      -521.41478771
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.55521720 eV

  energy without entropy =      -39.41052861  energy(sigma->0) =      -39.50698767


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




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



 eigenvalue-minimisations  :   672
 total energy-change (2. order) : 0.7564372E-01  (-0.5552102E+00)
 number of electron      64.0000028 magnetization 
 augmentation part       -0.7727276 magnetization 

 Broyden mixing:
  rms(total) = 0.14842E+01    rms(broyden)= 0.14822E+01
  rms(prec ) = 0.20471E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3694
  1.2712  0.6943  0.3484  0.3484  0.3415  0.2876  0.2876  0.1558  0.1558  0.1252
  0.0478

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1297.10009062
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.62025962
  PAW double counting   =      1707.15583262    -1607.02871049
  entropy T*S    EENTRO =        -0.13435711
  eigenvalues    EBANDS =      -522.39216115
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.47957347 eV

  energy without entropy =      -39.34521637  energy(sigma->0) =      -39.43478777


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




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



 eigenvalue-minimisations  :   674
 total energy-change (2. order) : 0.1581647E+00  (-0.3318569E+00)
 number of electron      64.0000059 magnetization 
 augmentation part       -0.2950889 magnetization 

 Broyden mixing:
  rms(total) = 0.11190E+01    rms(broyden)= 0.11166E+01
  rms(prec ) = 0.15848E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3806
  1.3486  0.9178  0.3486  0.3486  0.3546  0.3116  0.3116  0.1511  0.1511  0.1714
  0.1067  0.0450

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1294.56011436
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.40848473
  PAW double counting   =      1708.36273152    -1608.20722577
  entropy T*S    EENTRO =        -0.12353318
  eigenvalues    EBANDS =      -524.60140537
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.32140877 eV

  energy without entropy =      -39.19787560  energy(sigma->0) =      -39.28023105


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




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



 eigenvalue-minimisations  :   560
 total energy-change (2. order) :-0.2607445E+00  (-0.1985526E+00)
 number of electron      64.0000060 magnetization 
 augmentation part       -0.4976483 magnetization 

 Broyden mixing:
  rms(total) = 0.14406E+01    rms(broyden)= 0.14395E+01
  rms(prec ) = 0.19721E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.3995
  1.6251  0.8622  0.5261  0.3498  0.3498  0.3112  0.3112  0.2294  0.1538  0.1538
  0.1647  0.1110  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1293.71602466
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.23130062
  PAW double counting   =      1715.76991956    -1615.59441561
  entropy T*S    EENTRO =        -0.01785202
  eigenvalues    EBANDS =      -525.65473486
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.58215331 eV

  energy without entropy =      -39.56430129  energy(sigma->0) =      -39.57620263


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




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



 eigenvalue-minimisations  :   656
 total energy-change (2. order) : 0.2959790E+00  (-0.7141547E-01)
 number of electron      64.0000063 magnetization 
 augmentation part       -0.3580523 magnetization 

 Broyden mixing:
  rms(total) = 0.11265E+01    rms(broyden)= 0.11262E+01
  rms(prec ) = 0.15443E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4060
  1.5694  1.0600  0.6062  0.3517  0.3517  0.3115  0.3115  0.3264  0.1522  0.1522
  0.1821  0.1562  0.1074  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1293.05889108
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.13522765
  PAW double counting   =      1712.14989374    -1611.96564991
  entropy T*S    EENTRO =        -0.14052082
  eigenvalues    EBANDS =      -525.80588752
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.28617428 eV

  energy without entropy =      -39.14565347  energy(sigma->0) =      -39.23933401


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




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



 eigenvalue-minimisations  :   576
 total energy-change (2. order) : 0.6565395E-01  (-0.1648023E+00)
 number of electron      64.0000032 magnetization 
 augmentation part       -0.6224695 magnetization 

 Broyden mixing:
  rms(total) = 0.11889E+01    rms(broyden)= 0.11876E+01
  rms(prec ) = 0.16856E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4240
  1.7090  1.4648  0.3513  0.3513  0.4765  0.4765  0.3114  0.3114  0.1555  0.1555
  0.2036  0.1560  0.1074  0.0841  0.0449

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1294.11459641
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        53.11898518
  PAW double counting   =      1706.02626649    -1605.86414078
  entropy T*S    EENTRO =        -0.24332775
  eigenvalues    EBANDS =      -524.54336072
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.22052034 eV

  energy without entropy =      -38.97719259  energy(sigma->0) =      -39.13941109


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




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



 eigenvalue-minimisations  :   632
 total energy-change (2. order) :-0.1872504E+00  (-0.2883841E+00)
 number of electron      64.0000062 magnetization 
 augmentation part       -0.0986561 magnetization 

 Broyden mixing:
  rms(total) = 0.12765E+01    rms(broyden)= 0.12758E+01
  rms(prec ) = 0.17841E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4186
  1.8814  1.4122  0.3510  0.3510  0.4894  0.4894  0.3118  0.3118  0.2606  0.1553
  0.1553  0.1491  0.1491  0.1165  0.0449  0.0689

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1291.37990219
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.95082393
  PAW double counting   =      1702.93349358    -1602.74849732
  entropy T*S    EENTRO =        -0.07869845
  eigenvalues    EBANDS =      -527.48464389
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.40777069 eV

  energy without entropy =      -39.32907224  energy(sigma->0) =      -39.38153787


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




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



 eigenvalue-minimisations  :   548
 total energy-change (2. order) : 0.2564371E+00  (-0.1843052E-01)
 number of electron      64.0000059 magnetization 
 augmentation part       -0.2386592 magnetization 

 Broyden mixing:
  rms(total) = 0.95478E+00    rms(broyden)= 0.95465E+00
  rms(prec ) = 0.13259E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4310
  1.9944  1.4563  0.5260  0.5260  0.3508  0.3508  0.4138  0.3197  0.3197  0.1549
  0.1549  0.1964  0.1660  0.1660  0.1132  0.0449  0.0734

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.17152229
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.94817156
  PAW double counting   =      1702.82783019    -1602.65085717
  entropy T*S    EENTRO =        -0.17942297
  eigenvalues    EBANDS =      -526.32518656
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.15133359 eV

  energy without entropy =      -38.97191061  energy(sigma->0) =      -39.09152593


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




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



 eigenvalue-minimisations  :   584
 total energy-change (2. order) : 0.1243151E+00  (-0.2213547E-01)
 number of electron      64.0000049 magnetization 
 augmentation part       -0.4665617 magnetization 

 Broyden mixing:
  rms(total) = 0.45931E+00    rms(broyden)= 0.45887E+00
  rms(prec ) = 0.64098E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4372
  2.0736  1.4191  0.6470  0.6470  0.3508  0.3508  0.3206  0.3206  0.3378  0.3378
  0.1552  0.1552  0.2065  0.1586  0.1586  0.1130  0.0449  0.0727

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.99889701
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.95288504
  PAW double counting   =      1703.40320669    -1603.23463173
  entropy T*S    EENTRO =        -0.27982668
  eigenvalues    EBANDS =      -525.26940846
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.02701850 eV

  energy without entropy =      -38.74719181  energy(sigma->0) =      -38.93374294


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




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



 eigenvalue-minimisations  :   649
 total energy-change (2. order) :-0.7312321E-01  (-0.1774455E-01)
 number of electron      64.0000057 magnetization 
 augmentation part       -0.3716884 magnetization 

 Broyden mixing:
  rms(total) = 0.74725E+00    rms(broyden)= 0.74711E+00
  rms(prec ) = 0.10153E+01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4470
  1.9878  1.6421  0.7381  0.7381  0.3507  0.3507  0.4050  0.4050  0.3174  0.3174
  0.2306  0.1550  0.1550  0.1656  0.1517  0.1517  0.1133  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.37522794
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.87039457
  PAW double counting   =      1703.10391305    -1602.92846772
  entropy T*S    EENTRO =        -0.23266565
  eigenvalues    EBANDS =      -525.93774167
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.10014171 eV

  energy without entropy =      -38.86747606  energy(sigma->0) =      -39.02258649


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




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



 eigenvalue-minimisations  :   656
 total energy-change (2. order) : 0.4934639E-01  (-0.8281462E-02)
 number of electron      64.0000050 magnetization 
 augmentation part       -0.4276183 magnetization 

 Broyden mixing:
  rms(total) = 0.42288E+00    rms(broyden)= 0.42280E+00
  rms(prec ) = 0.58326E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4564
  1.9874  1.8393  0.8059  0.8059  0.3507  0.3507  0.4548  0.4548  0.3157  0.3157
  0.2590  0.1550  0.1550  0.1992  0.1548  0.1548  0.1386  0.1132  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.50831060
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.84654071
  PAW double counting   =      1701.78210727    -1601.60993808
  entropy T*S    EENTRO =        -0.26991464
  eigenvalues    EBANDS =      -525.69093365
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05079532 eV

  energy without entropy =      -38.78088069  energy(sigma->0) =      -38.96082378


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




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



 eigenvalue-minimisations  :   579
 total energy-change (2. order) : 0.9220002E-02  (-0.2097975E-01)
 number of electron      64.0000037 magnetization 
 augmentation part       -0.5994445 magnetization 

 Broyden mixing:
  rms(total) = 0.24002E+00    rms(broyden)= 0.23904E+00
  rms(prec ) = 0.31966E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4447
  2.0822  1.7220  0.8314  0.8314  0.3507  0.3507  0.4297  0.4297  0.3181  0.3181
  0.3353  0.2206  0.1550  0.1550  0.1548  0.1548  0.1576  0.1139  0.1104  0.0449
  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.82285994
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.83161063
  PAW double counting   =      1700.52132177    -1600.35224892
  entropy T*S    EENTRO =        -0.33826538
  eigenvalues    EBANDS =      -525.28078714
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.04157532 eV

  energy without entropy =      -38.70330994  energy(sigma->0) =      -38.92882019


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




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



 eigenvalue-minimisations  :   560
 total energy-change (2. order) :-0.3103654E-01  (-0.5383138E-02)
 number of electron      64.0000029 magnetization 
 augmentation part       -0.6733165 magnetization 

 Broyden mixing:
  rms(total) = 0.62869E+00    rms(broyden)= 0.62847E+00
  rms(prec ) = 0.84837E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4247
  2.0816  1.7172  0.8324  0.8324  0.3507  0.3507  0.4254  0.4254  0.3184  0.3184
  0.3390  0.2215  0.1549  0.1549  0.1547  0.1547  0.1567  0.1137  0.1089  0.0728
  0.0449  0.0141

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.90896391
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.83028739
  PAW double counting   =      1699.69140183    -1599.52208706
  entropy T*S    EENTRO =        -0.34823765
  eigenvalues    EBANDS =      -525.21466612
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.07261187 eV

  energy without entropy =      -38.72437422  energy(sigma->0) =      -38.95653265


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




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



 eigenvalue-minimisations  :   404
 total energy-change (2. order) : 0.1472419E-01  (-0.1410295E-03)
 number of electron      64.0000030 magnetization 
 augmentation part       -0.6624628 magnetization 

 Broyden mixing:
  rms(total) = 0.59750E+00    rms(broyden)= 0.59749E+00
  rms(prec ) = 0.80685E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4456
  2.0800  1.6685  0.9022  0.9022  0.3507  0.3507  0.4679  0.4679  0.4308  0.4308
  0.3125  0.3125  0.1550  0.1550  0.2195  0.2195  0.1571  0.1571  0.1603  0.1132
  0.1172  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.88119767
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.83791057
  PAW double counting   =      1699.65833395    -1599.48811776
  entropy T*S    EENTRO =        -0.34692499
  eigenvalues    EBANDS =      -525.23754544
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05788768 eV

  energy without entropy =      -38.71096269  energy(sigma->0) =      -38.94224602


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




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



 eigenvalue-minimisations  :   592
 total energy-change (2. order) :-0.3424371E-01  (-0.4773007E-01)
 number of electron      64.0000053 magnetization 
 augmentation part       -0.4569830 magnetization 

 Broyden mixing:
  rms(total) = 0.71037E+00    rms(broyden)= 0.70917E+00
  rms(prec ) = 0.96730E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.4821
  2.2387  1.7694  1.0292  1.0292  0.7824  0.3507  0.3507  0.5414  0.4451  0.4451
  0.4213  0.3128  0.3128  0.1550  0.1550  0.2061  0.2061  0.1568  0.1568  0.1557
  0.1131  0.1186  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.35295306
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.79638931
  PAW double counting   =      1700.99575648    -1600.81327413
  entropy T*S    EENTRO =        -0.22465216
  eigenvalues    EBANDS =      -525.89305149
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.09213139 eV

  energy without entropy =      -38.86747923  energy(sigma->0) =      -39.01724733


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




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



 eigenvalue-minimisations  :   432
 total energy-change (2. order) : 0.3510988E-01  (-0.5688119E-02)
 number of electron      64.0000048 magnetization 
 augmentation part       -0.5106900 magnetization 

 Broyden mixing:
  rms(total) = 0.40010E+00    rms(broyden)= 0.40007E+00
  rms(prec ) = 0.54302E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5084
  2.4219  2.2717  1.0447  0.8983  0.8983  0.3507  0.3507  0.4600  0.4600  0.4977
  0.4977  0.3905  0.3133  0.3133  0.1550  0.1550  0.2063  0.2063  0.1568  0.1568
  0.1565  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.38994917
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.76622086
  PAW double counting   =      1699.55807716    -1599.37538139
  entropy T*S    EENTRO =        -0.27770818
  eigenvalues    EBANDS =      -525.73793445
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05702151 eV

  energy without entropy =      -38.77931333  energy(sigma->0) =      -38.96445212


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




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



 eigenvalue-minimisations  :   576
 total energy-change (2. order) : 0.1195022E-01  (-0.6784587E-02)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5544460 magnetization 

 Broyden mixing:
  rms(total) = 0.58071E-01    rms(broyden)= 0.56137E-01
  rms(prec ) = 0.79072E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5270
  2.6705  2.3350  1.0537  1.0537  0.6834  0.6454  0.6454  0.3507  0.3507  0.4730
  0.4730  0.3991  0.3991  0.3132  0.3132  0.1550  0.1550  0.2069  0.2069  0.1568
  0.1568  0.1562  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.36348183
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.75452107
  PAW double counting   =      1698.16513456    -1597.98598630
  entropy T*S    EENTRO =        -0.32115290
  eigenvalues    EBANDS =      -525.69375955
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.04507129 eV

  energy without entropy =      -38.72391840  energy(sigma->0) =      -38.93802033


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




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



 eigenvalue-minimisations  :   422
 total energy-change (2. order) :-0.1048906E-01  (-0.5799321E-03)
 number of electron      64.0000042 magnetization 
 augmentation part       -0.5409419 magnetization 

 Broyden mixing:
  rms(total) = 0.61788E-01    rms(broyden)= 0.61756E-01
  rms(prec ) = 0.83795E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5727
  3.4384  2.2120  1.2001  1.2001  0.8771  0.8771  0.3507  0.3507  0.5541  0.5541
  0.4668  0.4668  0.3732  0.3732  0.3132  0.3132  0.1550  0.1550  0.2067  0.2067
  0.1568  0.1568  0.1562  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.30271656
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.74000692
  PAW double counting   =      1698.50520838    -1598.32449120
  entropy T*S    EENTRO =        -0.31372057
  eigenvalues    EBANDS =      -525.75950098
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05556035 eV

  energy without entropy =      -38.74183978  energy(sigma->0) =      -38.95098683


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




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



 eigenvalue-minimisations  :   482
 total energy-change (2. order) :-0.1371477E-02  (-0.1131787E-03)
 number of electron      64.0000042 magnetization 
 augmentation part       -0.5404206 magnetization 

 Broyden mixing:
  rms(total) = 0.82919E-01    rms(broyden)= 0.82908E-01
  rms(prec ) = 0.11270E+00
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.5922
  3.7131  2.2258  1.6725  0.9798  0.9798  0.7890  0.3507  0.3507  0.5906  0.5906
  0.4600  0.4600  0.5018  0.3132  0.3132  0.3741  0.3741  0.1550  0.1550  0.2068
  0.2068  0.1568  0.1568  0.1562  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.26157147
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73211970
  PAW double counting   =      1699.32124616    -1599.13951584
  entropy T*S    EENTRO =        -0.31246263
  eigenvalues    EBANDS =      -525.79640140
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05693183 eV

  energy without entropy =      -38.74446920  energy(sigma->0) =      -38.95277762


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




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



 eigenvalue-minimisations  :   626
 total energy-change (2. order) :-0.3420609E-03  (-0.2554919E-03)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5650533 magnetization 

 Broyden mixing:
  rms(total) = 0.51496E-01    rms(broyden)= 0.51455E-01
  rms(prec ) = 0.71708E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6235
  4.0733  2.4437  1.7587  0.9901  0.9901  0.8436  0.7333  0.7333  0.3507  0.3507
  0.4602  0.4602  0.4874  0.4874  0.3132  0.3132  0.3748  0.3748  0.1550  0.1550
  0.2068  0.2068  0.1568  0.1568  0.1562  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.35570428
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73673929
  PAW double counting   =      1700.12241511    -1599.94162957
  entropy T*S    EENTRO =        -0.31934693
  eigenvalues    EBANDS =      -525.69940116
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05727389 eV

  energy without entropy =      -38.73792696  energy(sigma->0) =      -38.95082491


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




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



 eigenvalue-minimisations  :   506
 total energy-change (2. order) :-0.6666321E-03  (-0.2460895E-04)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5651439 magnetization 

 Broyden mixing:
  rms(total) = 0.45930E-01    rms(broyden)= 0.45927E-01
  rms(prec ) = 0.64033E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6475
  4.0817  2.6835  1.9631  0.9837  0.9837  0.9662  0.9662  0.7091  0.3507  0.3507
  0.4612  0.4612  0.5571  0.5571  0.3132  0.3132  0.4335  0.3732  0.3732  0.1550
  0.1550  0.2068  0.2068  0.1568  0.1568  0.1562  0.1131  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.38240465
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73874207
  PAW double counting   =      1700.74641449    -1600.56618947
  entropy T*S    EENTRO =        -0.31930063
  eigenvalues    EBANDS =      -525.67485597
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05794052 eV

  energy without entropy =      -38.73863989  energy(sigma->0) =      -38.95150698


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




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



 eigenvalue-minimisations  :   538
 total energy-change (2. order) :-0.2782437E-03  (-0.2088159E-03)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5448595 magnetization 

 Broyden mixing:
  rms(total) = 0.40503E-01    rms(broyden)= 0.40387E-01
  rms(prec ) = 0.56673E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6493
  4.0992  2.7954  1.9669  1.2691  0.9242  0.9242  0.8922  0.6726  0.6726  0.3507
  0.3507  0.4606  0.4606  0.4932  0.4932  0.3132  0.3132  0.3779  0.3779  0.3768
  0.1550  0.1550  0.2068  0.2068  0.1568  0.1568  0.1562  0.1131  0.1185  0.0449
  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.29687468
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73591734
  PAW double counting   =      1700.74454848    -1600.56342100
  entropy T*S    EENTRO =        -0.31890430
  eigenvalues    EBANDS =      -525.75913825
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05821876 eV

  energy without entropy =      -38.73931446  energy(sigma->0) =      -38.95191733


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




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



 eigenvalue-minimisations  :   440
 total energy-change (2. order) :-0.4150547E-03  (-0.1769505E-03)
 number of electron      64.0000040 magnetization 
 augmentation part       -0.5623856 magnetization 

 Broyden mixing:
  rms(total) = 0.70505E-01    rms(broyden)= 0.70471E-01
  rms(prec ) = 0.95910E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6679
  4.3274  2.8625  2.0318  1.4588  1.0021  0.9301  0.9301  0.7007  0.7007  0.3507
  0.3507  0.4612  0.4612  0.5303  0.5303  0.3132  0.3132  0.4176  0.4176  0.3696
  0.3696  0.1550  0.1550  0.2068  0.2068  0.1568  0.1568  0.1562  0.1131  0.1185
  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.35168686
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73805393
  PAW double counting   =      1700.76744705    -1600.58697133
  entropy T*S    EENTRO =        -0.32505397
  eigenvalues    EBANDS =      -525.70007628
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05863382 eV

  energy without entropy =      -38.73357985  energy(sigma->0) =      -38.95028249


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




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



 eigenvalue-minimisations  :   395
 total energy-change (2. order) : 0.3752947E-03  (-0.1983625E-04)
 number of electron      64.0000040 magnetization 
 augmentation part       -0.5670562 magnetization 

 Broyden mixing:
  rms(total) = 0.54019E-01    rms(broyden)= 0.54014E-01
  rms(prec ) = 0.73245E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.6909
  4.5525  3.0610  2.0499  1.6176  1.0415  1.0415  0.7271  0.7271  0.7597  0.7597
  0.3507  0.3507  0.4615  0.4615  0.5420  0.5420  0.3132  0.3132  0.4382  0.3813
  0.3813  0.3823  0.1550  0.1550  0.2068  0.2068  0.1568  0.1568  0.1562  0.1131
  0.1185  0.0449  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.38084845
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.74018947
  PAW double counting   =      1700.85729028    -1600.67679009
  entropy T*S    EENTRO =        -0.32516514
  eigenvalues    EBANDS =      -525.67258826
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05825852 eV

  energy without entropy =      -38.73309339  energy(sigma->0) =      -38.94987014


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




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



 eigenvalue-minimisations  :   389
 total energy-change (2. order) : 0.4461640E-04  (-0.2359024E-03)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5506213 magnetization 

 Broyden mixing:
  rms(total) = 0.24764E-01    rms(broyden)= 0.24640E-01
  rms(prec ) = 0.33546E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7200
  4.9813  3.0175  2.0765  2.0765  1.1212  1.1212  0.9065  0.9065  0.6500  0.6500
  0.3507  0.3507  0.5851  0.5851  0.4615  0.4615  0.4467  0.4467  0.3132  0.3132
  0.3763  0.3763  0.3633  0.1550  0.1550  0.2068  0.2068  0.1568  0.1568  0.1562
  0.0449  0.1131  0.1185  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.34429683
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73874950
  PAW double counting   =      1700.88441139    -1600.70317000
  entropy T*S    EENTRO =        -0.31851879
  eigenvalues    EBANDS =      -525.71504284
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05821391 eV

  energy without entropy =      -38.73969512  energy(sigma->0) =      -38.95204098


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




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



 eigenvalue-minimisations  :   333
 total energy-change (2. order) :-0.1968928E-03  (-0.4042694E-04)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5557843 magnetization 

 Broyden mixing:
  rms(total) = 0.89771E-02    rms(broyden)= 0.89606E-02
  rms(prec ) = 0.12442E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7208
  5.0724  3.0635  2.0750  2.0750  1.1190  1.1190  0.9241  0.9241  0.3507  0.3507
  0.6701  0.6701  0.4615  0.4615  0.5772  0.5772  0.5079  0.5079  0.3132  0.3132
  0.4253  0.3803  0.3803  0.3650  0.1550  0.1550  0.2068  0.2068  0.1568  0.1568
  0.1562  0.0449  0.1131  0.1185  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.34990719
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73799796
  PAW double counting   =      1700.71820465    -1600.53686828
  entropy T*S    EENTRO =        -0.32099749
  eigenvalues    EBANDS =      -525.70649411
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05841080 eV

  energy without entropy =      -38.73741331  energy(sigma->0) =      -38.95141164


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




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



 eigenvalue-minimisations  :   230
 total energy-change (2. order) :-0.1368468E-04  (-0.1226113E-05)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5550477 magnetization 

 Broyden mixing:
  rms(total) = 0.78148E-02    rms(broyden)= 0.78141E-02
  rms(prec ) = 0.10945E-01
  weight for this iteration     100.00

 eigenvalues of (default mixing * dielectric matrix)
  average eigenvalue GAMMA=   0.7382
  5.2276  3.1618  2.1456  2.1456  1.1095  1.1095  0.9761  0.9761  0.6784  0.6784
  0.7230  0.7230  0.3507  0.3507  0.4615  0.4615  0.5719  0.5719  0.3132  0.3132
  0.4328  0.4328  0.3762  0.3762  0.3671  0.1550  0.1550  0.2068  0.2068  0.1568
  0.1568  0.1562  0.0449  0.1131  0.1185  0.0728

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.34826105
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73805743
  PAW double counting   =      1700.65537841    -1600.47387894
  entropy T*S    EENTRO =        -0.32077724
  eigenvalues    EBANDS =      -525.70859675
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05842448 eV

  energy without entropy =      -38.73764724  energy(sigma->0) =      -38.95149874


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




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



 eigenvalue-minimisations  :   230
 total energy-change (2. order) :-0.7448833E-05  (-0.5463241E-05)
 number of electron      64.0000041 magnetization 
 augmentation part       -0.5550477 magnetization 

 Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         5.46929015
  Ewald energy   TEWEN  =      -460.24633778
  -Hartree energ DENC   =     -1292.35733546
  -exchange      EXHF   =         0.00000000
  -V(xc)+E(xc)   XCENC  =        52.73834161
  PAW double counting   =      1700.65316055    -1600.47168608
  entropy T*S    EENTRO =        -0.32028591
  eigenvalues    EBANDS =      -525.70028029
  atomic energy  EATOM  =      2081.17670128
  Solvation  Ediel_sol  =         0.00000000
  ---------------------------------------------------
  free energy    TOTEN  =       -39.05843193 eV

  energy without entropy =      -38.73814602  energy(sigma->0) =      -38.95166996


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




 average (electrostatic) potential at core
  the test charge radii are     0.7089  0.9892  0.7215  0.5201
  (the norm of the test charge is              1.0000)
       1 -73.8436       2 -74.0318       3 -74.1483       4 -96.2341       5 -95.8685
       6 -96.0849       7 -95.5560       8 -94.7963       9 -95.7409      10 -79.0099
      11 -39.8927      12 -40.6886      13 -39.8692      14 -40.7882      15 -40.0290
      16 -40.3449      17 -40.3533      18 -40.8554      19 -40.3631      20 -42.8837
      21 -40.7718      22 -40.9490      23 -40.9447      24 -40.3256      25 -41.0075
      26 -40.7523      27 -41.0344      28 -40.6254      29 -41.0151
 
 
 
 E-fermi :  -4.8666     XC(G=0):  -3.1465     alpha+bet : -1.0645


 k-point     1 :       0.0000    0.0000    0.0000
  band No.  band energies     occupation 
      1     -19.7478      2.00000
      2     -16.5884      2.00000
      3     -16.3700      2.00000
      4     -16.3472      2.00000
      5     -12.9918      2.00000
      6     -11.7470      2.00000
      7     -11.6443      2.00000
      8     -11.5804      2.00000
      9     -11.5414      2.00000
     10     -10.9905      2.00000
     11      -7.2871      2.00000
     12      -7.1909      2.00000
     13      -6.5785      2.00000
     14      -6.5029      2.00000
     15      -6.4242      2.00000
     16      -6.3460      2.00000
     17      -5.9790      2.00000
     18      -5.6997      2.00000
     19      -5.5811      2.00001
     20      -5.4139      2.00076
     21      -5.3372      2.00436
     22      -5.2228      2.03033
     23      -5.1542      2.06060
     24      -5.0672      2.05080
     25      -5.0448      2.01962
     26      -5.0010      1.89922
     27      -4.9849      1.83244
     28      -4.9664      1.73910
     29      -4.9319      1.52149
     30      -4.9183      1.42187
     31      -4.8930      1.22079
     32      -4.8873      1.17427
     33      -4.8565      0.91417
     34      -4.8443      0.81216
     35      -4.8289      0.68727
     36      -4.8049      0.50444
     37      -4.7711      0.28463
     38      -4.7483      0.16767
     39      -4.6515     -0.06260
     40      -4.6174     -0.07078
     41      -4.5503     -0.04783
     42      -4.4362     -0.00949
     43      -4.3266     -0.00090
     44      -4.2459     -0.00010
     45      -4.1820     -0.00001
     46      -4.0797     -0.00000
     47      -3.9550     -0.00000

 k-point     2 :       0.0000    0.5000    0.0000
  band No.  band energies     occupation 
      1     -19.7477      2.00000
      2     -16.5884      2.00000
      3     -16.3700      2.00000
      4     -16.3471      2.00000
      5     -12.9917      2.00000
      6     -11.7470      2.00000
      7     -11.6443      2.00000
      8     -11.5804      2.00000
      9     -11.5413      2.00000
     10     -10.9905      2.00000
     11      -7.2880      2.00000
     12      -7.1911      2.00000
     13      -6.5789      2.00000
     14      -6.5008      2.00000
     15      -6.4262      2.00000
     16      -6.3439      2.00000
     17      -5.9823      2.00000
     18      -5.6942      2.00000
     19      -5.5866      2.00000
     20      -5.4033      2.00098
     21      -5.3554      2.00297
     22      -5.1998      2.04009
     23      -5.1563      2.05977
     24      -5.0693      2.05293
     25      -5.0408      2.01209
     26      -4.9958      1.87923
     27      -4.9897      1.85359
     28      -4.9524      1.65706
     29      -4.9160      1.40443
     30      -4.9112      1.36710
     31      -4.9006      1.28332
     32      -4.8688      1.01884
     33      -4.8645      0.98188
     34      -4.8561      0.91093
     35      -4.8285      0.68363
     36      -4.8059      0.51165
     37      -4.7701      0.27916
     38      -4.7513      0.18145
     39      -4.6741     -0.04165
     40      -4.6299     -0.07036
     41      -4.5701     -0.05685
     42      -4.4330     -0.00895
     43      -4.3267     -0.00091
     44      -4.2431     -0.00010
     45      -4.1769     -0.00001
     46      -4.0823     -0.00000
     47      -3.9651     -0.00000

 k-point     3 :       0.0000    0.0000    0.5000
  band No.  band energies     occupation 
      1     -19.7477      2.00000
      2     -16.5884      2.00000
      3     -16.3700      2.00000
      4     -16.3471      2.00000
      5     -12.9917      2.00000
      6     -11.7468      2.00000
      7     -11.6444      2.00000
      8     -11.5804      2.00000
      9     -11.5414      2.00000
     10     -10.9905      2.00000
     11      -7.2869      2.00000
     12      -7.1909      2.00000
     13      -6.5779      2.00000
     14      -6.5031      2.00000
     15      -6.4241      2.00000
     16      -6.3455      2.00000
     17      -5.9788      2.00000
     18      -5.7008      2.00000
     19      -5.5816      2.00001
     20      -5.4145      2.00074
     21      -5.3377      2.00431
     22      -5.2031      2.03863
     23      -5.1271      2.06925
     24      -5.0980      2.06936
     25      -5.0653      2.04890
     26      -5.0221      1.96805
     27      -4.9794      1.80622
     28      -4.9474      1.62572
     29      -4.9379      1.56316
     30      -4.9296      1.50470
     31      -4.9047      1.31599
     32      -4.8763      1.08207
     33      -4.8507      0.86591
     34      -4.8365      0.74813
     35      -4.8265      0.66800
     36      -4.7861      0.37579
     37      -4.7664      0.25853
     38      -4.7528      0.18874
     39      -4.6495     -0.06375
     40      -4.6187     -0.07085
     41      -4.5562     -0.05058
     42      -4.4381     -0.00981
     43      -4.3294     -0.00097
     44      -4.2463     -0.00010
     45      -4.1829     -0.00001
     46      -4.0826     -0.00000
     47      -3.9643     -0.00000

 k-point     4 :       0.0000    0.5000    0.5000
  band No.  band energies     occupation 
      1     -19.7476      2.00000
      2     -16.5883      2.00000
      3     -16.3700      2.00000
      4     -16.3471      2.00000
      5     -12.9917      2.00000
      6     -11.7468      2.00000
      7     -11.6445      2.00000
      8     -11.5804      2.00000
      9     -11.5413      2.00000
     10     -10.9905      2.00000
     11      -7.2878      2.00000
     12      -7.1911      2.00000
     13      -6.5782      2.00000
     14      -6.5010      2.00000
     15      -6.4261      2.00000
     16      -6.3434      2.00000
     17      -5.9821      2.00000
     18      -5.6950      2.00000
     19      -5.5871      2.00000
     20      -5.4022      2.00101
     21      -5.3555      2.00296
     22      -5.1616      2.05751
     23      -5.1535      2.06091
     24      -5.1031      2.07033
     25      -5.0683      2.05191
     26      -5.0173      1.95419
     27      -4.9607      1.70669
     28      -4.9500      1.64232
     29      -4.9186      1.42376
     30      -4.9144      1.39210
     31      -4.8991      1.27084
     32      -4.8758      1.07790
     33      -4.8591      0.93678
     34      -4.8475      0.83901
     35      -4.8231      0.64135
     36      -4.8000      0.46923
     37      -4.7628      0.23896
     38      -4.7515      0.18264
     39      -4.6726     -0.04343
     40      -4.6307     -0.07025
     41      -4.5755     -0.05917
     42      -4.4339     -0.00909
     43      -4.3297     -0.00098
     44      -4.2431     -0.00010
     45      -4.1777     -0.00001
     46      -4.0863     -0.00000
     47      -3.9557     -0.00000


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


 soft charge-density along one line, spin component           1
         0         1         2         3         4         5         6         7         8         9
 total charge-density along one line
 
 pseudopotential strength for first ion, spin component:           1
 13.920  16.716  -0.000   0.001   0.000   0.000  -0.002   0.000
 16.716  20.075  -0.000   0.001   0.000   0.001  -0.002   0.000
 -0.000  -0.000  -7.355   0.002   0.001 -10.202   0.003   0.002
  0.001   0.001   0.002  -7.356   0.005   0.003 -10.203   0.008
  0.000   0.000   0.001   0.005  -7.338   0.002   0.008 -10.176
  0.000   0.001 -10.202   0.003   0.002 -13.504   0.004   0.003
 -0.002  -0.002   0.003 -10.203   0.008   0.004 -13.505   0.012
  0.000   0.000   0.002   0.008 -10.176   0.003   0.012 -13.463
 total augmentation occupancy for first ion, spin component:           1
  2.805  -0.461  -0.069   0.242  -0.002   0.012  -0.033  -0.001
 -0.461   0.187   0.089  -0.320   0.006  -0.011   0.032  -0.001
 -0.069   0.089   1.093  -0.040   0.063  -0.048   0.027  -0.023
  0.242  -0.320  -0.040   1.179   0.029   0.027  -0.091   0.009
 -0.002   0.006   0.063   0.029   1.287  -0.023   0.009  -0.060
  0.012  -0.011  -0.048   0.027  -0.023   0.003  -0.003   0.002
 -0.033   0.032   0.027  -0.091   0.009  -0.003   0.009  -0.001
 -0.001  -0.001  -0.023   0.009  -0.060   0.002  -0.001   0.004


------------------------ aborting loop because EDIFF is reached ----------------------------------------



  FORCE on cell =-STRESS in cart. coord.  units (eV):
  Direction    XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------------
  Alpha Z     5.46929     5.46929     5.46929
  Ewald    -229.32337    90.57282  -321.49787     6.34005   -54.72698     6.16154
  Hartree   362.43046   654.66472   275.51678    -3.48277   -55.44843    -3.79211
  E(xc)    -195.13087  -195.32328  -195.10324    -0.09887    -0.33217    -0.04272
  Local    -766.80443 -1373.94128  -584.08438    -8.32608   108.64082    -7.65063
  n-local   165.64893   171.90177   164.18534     3.15031     4.61524     1.38825
  augment   -33.34960   -34.10767   -32.88063    -0.62842    -0.51191    -0.15404
  Kinetic   664.83486   655.68347   661.17933     7.45526     4.31552     4.92359
  Fock        0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  -------------------------------------------------------------------------------------
  Total      -26.2247351    -25.0801653    -27.2153732      4.4094874      6.5520856      0.8338842
  in kB      -19.4521648    -18.6031816    -20.1869694      3.2707318      4.8600014      0.6185326
  external PRESSURE =     -19.4141053 kB  Pullay stress =       0.0000000 kB


 VOLUME and BASIS-vectors are now :
 -----------------------------------------------------------------------------
  energy-cutoff  :      400.00
  volume of cell :     2160.00
      direct lattice vectors                 reciprocal lattice vectors
    15.000000000  0.000000000  0.000000000     0.066666667  0.000000000  0.000000000
     0.000000000 12.000000000  0.000000000     0.000000000  0.083333333  0.000000000
     0.000000000  0.000000000 12.000000000     0.000000000  0.000000000  0.083333333

  length of vectors
    15.000000000 12.000000000 12.000000000     0.066666667  0.083333333  0.083333333


 FORCES acting on ions
    electron-ion (+dipol)            ewald-force                    non-local-force                 convergence-correction
 -----------------------------------------------------------------------------------------------
   0.152E+02 0.361E+01 0.196E+02   -.146E+02 -.122E+01 -.285E+02   -.638E-01 -.360E+01 0.136E+02   -.150E-02 -.631E-02 0.887E-02
   0.174E+02 0.210E+02 -.232E+02   -.155E+02 -.235E+02 0.256E+02   -.295E+01 0.402E+01 -.353E+01   -.102E-01 0.637E-02 0.503E-02
   -.179E+02 -.135E+02 0.828E+01   0.173E+02 0.144E+02 -.790E+01   0.117E+01 -.164E+01 -.753E+00   0.175E-01 0.185E-02 0.834E-02
   0.712E+00 -.147E+02 -.245E+02   -.599E-01 0.153E+02 0.264E+02   -.123E+01 -.996E+00 -.313E+01   -.143E-02 -.156E-02 0.992E-02
   -.266E+02 0.480E+01 0.222E+02   0.229E+02 -.696E+01 -.204E+02   0.511E+01 0.311E+01 -.281E+01   0.164E-01 -.306E-02 0.332E-02
   -.235E+02 -.403E+01 -.166E+02   0.250E+02 0.473E+01 0.174E+02   -.228E+01 -.108E+01 -.112E+01   0.614E-02 0.321E-02 -.467E-02
   -.297E+02 0.167E+01 -.445E+02   0.302E+02 -.224E+01 0.476E+02   -.175E+01 0.685E+00 -.536E+01   0.388E-02 -.228E-02 0.125E-02
   0.372E+02 0.256E+02 0.165E+02   -.393E+02 -.288E+02 -.133E+02   -.164E+01 -.296E+01 -.785E+01   -.120E-01 -.226E-03 0.170E-02
   0.207E+02 -.939E+01 0.193E+02   -.234E+02 0.120E+02 -.180E+02   0.462E+01 -.466E+01 -.213E+01   -.101E-01 0.180E-01 -.189E-02
   -.686E+01 0.712E+01 0.246E+02   0.740E+01 -.720E+01 -.248E+02   -.795E+00 0.177E+00 0.209E+00   0.254E-01 -.181E-02 0.602E-02
   -.261E+01 0.470E+01 0.105E+02   0.221E+01 -.419E+01 -.102E+02   -.820E+00 0.131E+01 0.231E+00   -.403E-02 0.141E-01 -.642E-02
   -.318E+01 -.165E+01 0.270E+01   0.312E+01 0.160E+01 -.261E+01   0.279E-02 0.141E-01 -.759E-02   0.891E-03 -.295E-02 0.372E-02
   -.136E+02 -.306E+01 0.128E+02   0.129E+02 0.259E+01 -.121E+02   -.944E+00 -.661E+00 0.797E+00   0.675E-02 0.277E-02 0.136E-02
   0.980E+01 -.344E+01 -.911E+01   -.991E+01 0.266E+01 0.856E+01   0.795E-02 -.698E-01 -.150E+00   -.436E-03 -.165E-02 0.505E-04
   -.192E+01 0.421E+00 -.681E+00   0.179E+01 -.409E+00 0.100E+01   -.366E-01 0.444E-01 0.791E-01   -.395E-03 0.128E-02 -.236E-02
   -.861E+01 -.117E+01 0.404E+01   0.831E+01 0.127E+01 -.405E+01   -.858E-01 0.109E+00 0.237E-01   0.575E-02 0.267E-02 0.460E-03
   -.417E+01 -.626E+00 -.440E+01   0.446E+01 0.737E+00 0.438E+01   0.305E-01 -.210E-01 0.296E-01   0.474E-02 -.275E-03 -.745E-03
   0.357E+01 0.192E+01 -.390E+01   -.338E+01 -.162E+01 0.380E+01   0.991E-02 0.650E-01 -.561E-02   0.143E-03 -.323E-03 -.740E-03
   0.744E+01 -.156E+01 0.456E+01   -.686E+01 0.164E+01 -.380E+01   0.776E-01 -.649E-01 0.238E+00   -.795E-02 -.438E-02 0.609E-02
   -.201E+01 -.178E+02 -.122E+02   0.734E+01 0.281E+02 0.180E+02   -.185E+01 -.376E+01 -.236E+01   -.345E-02 -.274E-03 0.658E-03
   -.799E+00 -.400E+01 0.165E+01   0.884E+00 0.348E+01 -.189E+01   0.819E-01 -.144E+00 0.288E-01   -.115E-01 0.493E-02 0.497E-04
   0.354E+01 -.320E+01 0.121E+02   -.339E+01 0.261E+01 -.114E+02   0.698E-01 -.142E+00 0.303E+00   -.398E-02 0.175E-03 0.213E-03
   0.126E+02 0.118E+02 -.208E+02   -.116E+02 -.105E+02 0.187E+02   0.468E+00 0.562E+00 -.668E+00   -.175E-02 -.237E-02 0.170E-02
   0.724E+01 -.307E+01 -.692E+01   -.696E+01 0.307E+01 0.688E+01   0.467E-01 -.794E-01 -.109E-01   -.245E-02 0.824E-03 0.122E-03
   0.150E+02 -.313E+01 0.830E+01   -.143E+02 0.236E+01 -.679E+01   0.165E+00 -.130E+00 0.272E+00   0.264E-03 -.258E-02 0.380E-02
   -.413E+01 -.213E+01 -.941E+00   0.376E+01 0.403E+00 0.101E+01   -.120E+00 -.304E+00 0.213E-01   0.849E-03 -.868E-03 -.105E-02
   -.314E+01 0.428E+01 -.262E+01   0.327E+01 -.273E+01 0.264E+01   0.571E-01 0.357E+00 0.149E-01   0.179E-02 -.244E-02 -.304E-02
   0.397E+01 0.385E+01 0.919E+01   -.432E+01 -.295E+01 -.890E+01   -.349E-01 0.135E+00 0.890E-01   -.258E-02 -.159E-03 -.276E-03
   -.282E+01 0.532E+01 0.773E+01   0.269E+01 -.463E+01 -.724E+01   -.269E-01 0.543E-01 0.451E-01   -.357E-02 0.858E-02 0.535E-02
 -----------------------------------------------------------------------------------------------
   0.269E+01 0.964E+01 0.138E+02   0.799E-14 -.284E-13 -.266E-14   -.270E+01 -.967E+01 -.139E+02   0.133E-01 0.312E-01 0.468E-01
 
 
 POSITION                                       TOTAL-FORCE (eV/Angst)
 -----------------------------------------------------------------------------------
      1.16076     11.77219      7.69905         0.514126     -1.207016      4.771586
     14.40733      7.87137      9.37735        -1.062221      1.492314     -1.108441
      7.12403      4.98276      5.22153         0.575838     -0.711638     -0.365125
      1.30069      3.76487      9.25154        -0.577789     -0.369196     -1.227667
      3.79014      0.73083      5.89776         1.385615      0.953312     -1.012689
      6.51853      1.44305     10.03715        -0.765812     -0.372093     -0.337545
      2.85294      0.13609      9.80244        -1.276759      0.114372     -2.188942
     11.94375      1.40172      7.02410        -3.710642     -6.148942     -4.673706
      0.41306      9.12526      3.96970         1.928946     -2.036982     -0.865931
      2.18875      5.29631      4.67815        -0.222232      0.100720      0.073748
      1.37559      7.50112      3.66113        -1.228418      1.840475      0.518042
      7.61773      7.86681      2.11478        -0.054898     -0.042931      0.087358
      5.13556      1.61742      4.78152        -1.624985     -1.126200      1.465300
     12.64351     10.75163      9.75321        -0.098576     -0.847923     -0.699883
      1.71244      9.66450      0.70574        -0.163992      0.058332      0.400701
      4.12074      8.93003      4.06455        -0.377647      0.211106      0.015661
      4.48157      6.61444      8.84735         0.325317      0.090128      0.010627
     10.09609      9.09444      9.99827         0.206266      0.368845     -0.108040
     14.96006      3.80244      6.43413         0.656501      0.004998      1.004013
     12.44385      2.38029      7.61917         3.482057      6.468876      3.432508
     10.01378      3.78395      5.57741         0.155594     -0.656560     -0.209479
     11.74067      2.20158      4.43734         0.209663     -0.735625      0.988343
      0.45163     10.88017      9.02287         1.389975      1.865158     -2.749623
     12.28453      4.22302      9.34629         0.324558     -0.072910     -0.049688
     13.29275      9.36800      8.26343         0.872803     -0.895333      1.783349
      7.44342      9.76844      8.87468        -0.493980     -2.031425      0.094157
      6.86071      7.91808      8.83900         0.181439      1.908065      0.035029
     12.77778      0.12460      4.20353        -0.386454      1.030293      0.379514
      8.32691      2.47375      3.77377        -0.164291      0.747783      0.536824
 -----------------------------------------------------------------------------------
    total drift:                                0.005093      0.003793     -0.009285


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



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

  energy  without entropy=      -38.7381460215  energy(sigma->0) =      -38.95166996
 


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




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



volume of typ            1:     0.2 %

volume of typ            2:     1.6 %

volume of typ            3:     0.1 %

volume of typ            4:     0.1 %
 
total charge     
 
# of ion       s       p       d       tot
------------------------------------------
    1        1.123   1.595   0.001   2.719
    2        1.137   1.566   0.000   2.703
    3        1.142   1.554   0.000   2.695
    4        0.959   0.423   0.001   1.383
    5        0.947   0.446   0.010   1.403
    6        0.948   0.468   0.000   1.417
    7        0.938   0.446   0.004   1.388
    8        0.922   0.658   0.058   1.639
    9        0.939   0.459   0.009   1.407
   10        1.324   2.530   0.000   3.854
   11        0.111   0.000   0.000   0.112
   12        0.097   0.000   0.000   0.097
   13        0.116   0.000   0.000   0.116
   14        0.094   0.000   0.000   0.094
   15        0.100   0.000   0.000   0.100
   16        0.098   0.000   0.000   0.098
   17        0.098   0.000   0.000   0.098
   18        0.094   0.000   0.000   0.094
   19        0.096   0.000   0.000   0.096
   20        0.188   0.001   0.000   0.190
   21        0.093   0.000   0.000   0.093
   22        0.092   0.000   0.000   0.092
   23        0.094   0.000   0.000   0.094
   24        0.100   0.000   0.000   0.100
   25        0.088   0.000   0.000   0.088
   26        0.094   0.000   0.000   0.094
   27        0.088   0.000   0.000   0.088
   28        0.097   0.000   0.000   0.097
   29        0.093   0.000   0.000   0.093
--------------------------------------------------
tot          12.31   10.15    0.08   22.54
 

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

   base      :      30000. kBytes
   nonlr-proj:       3185. kBytes
   fftplans  :      79712. kBytes
   grid      :     179635. kBytes
   one-center:         89. kBytes
   wavefun   :     118731. kBytes
 
  
  
 General timing and accounting informations for this job:
 ========================================================
  
                  Total CPU time used (sec):       39.148
                            User time (sec):       37.469
                          System time (sec):        1.679
                         Elapsed time (sec):       39.261
  
                   Maximum memory used (kb):     1233064.
                   Average memory used (kb):          N/A
  
                          Minor page faults:       218781
                          Major page faults:            0
                 Voluntary context switches:          507