vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.08.31 02:42:54 running on 2 total cores distrk: each k-point on 2 cores, 1 groups distr: one band on NCORE= 1 cores, 2 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 200 ALGO = Normal (blocked Davidson) ISPIN = 2 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 0.73 1.11 0.32 0.75 NPAR = 2 POTCAR: PAW_PBE O 08Apr2002 POTCAR: PAW_PBE Si 05Jan2001 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE N 08Apr2002 ----------------------------------------------------------------------------- | | | W W AA RRRRR N N II N N GGGG !!! | | W W A A R R NN N II NN N G G !!! | | W W A A R R N N N II N N N G !!! | | W WW W AAAAAA RRRRR N N N II N N N G GGG ! | | WW WW A A R R N NN II N NN G G | | W W A A R R N N II N N GGGG !!! | | | | You use a magnetic or noncollinear calculation, but did not specify | | the initial magnetic moment with the MAGMOM tag. Note that a | | default of 1 will be used for all atoms. This ferromagnetic setup | | may break the symmetry of the crystal, in particular it may rule | | out finding an antiferromagnetic solution. Thence, we recommend | | setting the initial magnetic moment manually or verifying carefully | | that this magnetic setup is desired. | | | ----------------------------------------------------------------------------- POTCAR: PAW_PBE 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 Si 05Jan2001 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 POTCAR: PAW_PBE 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 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 = 19.84 optimisation between [QCUT,QGAM] = [ 10.12, 20.44] = [ 28.68,116.96] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 7 10.119 159.560 0.56E-04 0.22E-03 0.45E-07 0 7 10.119 115.863 0.56E-04 0.21E-03 0.45E-07 1 7 10.119 88.339 0.34E-03 0.49E-03 0.11E-06 1 7 10.119 48.592 0.33E-03 0.48E-03 0.11E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 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 PAW_PBE O 08Apr2002 : energy of atom 1 EATOM= -432.3788 kinetic energy error for atom= 0.1156 (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 H 15Jun2001 : energy of atom 3 EATOM= -12.4884 kinetic energy error for atom= 0.0098 (will be added to EATOM!!) PAW_PBE N 08Apr2002 : energy of atom 4 EATOM= -264.5486 kinetic energy error for atom= 0.0736 (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.267 0.049 0.055- 2 0.066 0.156 0.521- 3 0.395 0.140 0.721- 4 0.335 0.189 0.555- 5 0.917 0.097 0.389- 6 0.361 0.947 0.403- 7 0.434 0.461 0.722- 8 0.455 0.380 0.594- 9 0.672 0.311 0.832- 10 0.825 0.848 0.847- 11 0.788 0.594 0.335- 12 0.722 0.266 0.522- LATTYP: Found a simple cubic cell. ALAT = 10.0000000000 Lattice vectors: A1 = ( 10.0000000000, 0.0000000000, 0.0000000000) A2 = ( 0.0000000000, 10.0000000000, 0.0000000000) A3 = ( 0.0000000000, 0.0000000000, 10.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 cubic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 48 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 cubic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 48 trial point group operations. The dynamic configuration has the point symmetry C_1 . Analysis of constrained symmetry for selective dynamics: ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple cubic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 48 trial point group operations. The constrained configuration has the point symmetry C_1 . Analysis of structural, dynamic, and magnetic symmetry: ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple cubic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 48 trial point group operations. The overall configuration has the point symmetry C_1 . Subroutine INISYM returns: Found 1 space group operations (whereof 1 operations are pure point group operations), and found 1 'primitive' translations ---------------------------------------------------------------------------------------- Primitive cell volume of cell : 1000.0000 direct lattice vectors reciprocal lattice vectors 10.000000000 0.000000000 0.000000000 0.100000000 0.000000000 0.000000000 0.000000000 10.000000000 0.000000000 0.000000000 0.100000000 0.000000000 0.000000000 0.000000000 10.000000000 0.000000000 0.000000000 0.100000000 length of vectors 10.000000000 10.000000000 10.000000000 0.100000000 0.100000000 0.100000000 position of ions in fractional coordinates (direct lattice) 0.267301560 0.049403020 0.054937530 0.066028380 0.155839470 0.521159270 0.395211750 0.140125940 0.720536630 0.335371870 0.188510980 0.555287870 0.916538830 0.096963150 0.388780050 0.361391410 0.947493740 0.403139610 0.434455570 0.460647040 0.722108890 0.454743330 0.380032410 0.593817340 0.671585560 0.311445340 0.832473460 0.824525060 0.848416220 0.847359760 0.787928270 0.593699900 0.335250820 0.722268350 0.265942330 0.522266720 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 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 2 2 2 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.050000000 0.000000000 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.500000000 Length of vectors 0.050000000 0.050000000 0.050000000 Shift w.r.t. Gamma in fractional coordinates (k-lattice) 0.000000000 0.000000000 0.000000000 Subroutine IBZKPT returns following result: =========================================== Found 8 irreducible k-points: Following reciprocal coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.500000 0.000000 1.000000 0.000000 0.000000 0.500000 1.000000 0.500000 0.500000 0.000000 1.000000 0.000000 0.500000 0.500000 1.000000 0.500000 0.000000 0.500000 1.000000 0.500000 0.500000 0.500000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.050000 0.000000 0.000000 1.000000 0.000000 0.050000 0.000000 1.000000 0.000000 0.000000 0.050000 1.000000 0.050000 0.050000 0.000000 1.000000 0.000000 0.050000 0.050000 1.000000 0.050000 0.000000 0.050000 1.000000 0.050000 0.050000 0.050000 1.000000 -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 8 k-points in BZ NKDIM = 8 number of bands NBANDS= 26 number of dos NEDOS = 301 number of ions NIONS = 12 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 125000 max r-space proj IRMAX = 2395 max aug-charges IRDMAX= 4378 dimension x,y,z NGX = 50 NGY = 50 NGZ = 50 dimension x,y,z NGXF= 100 NGYF= 100 NGZF= 100 support grid NGXF= 100 NGYF= 100 NGZF= 100 ions per type = 1 2 8 1 NGX,Y,Z is equivalent to a cutoff of 8.31, 8.31, 8.31 a.u. NGXF,Y,Z is equivalent to a cutoff of 16.62, 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 = 2 spin polarized calculation? LNONCOLLINEAR = F non collinear calculations LSORBIT = F spin-orbit coupling INIWAV = 1 electr: 0-lowe 1-rand 2-diag LASPH = F aspherical Exc in radial PAW Electronic Relaxation 1 ENCUT = 400.0 eV 29.40 Ry 5.42 a.u. 16.31 16.31 16.31*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.229E-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 = 16.00 28.09 1.00 14.00 Ionic Valenz ZVAL = 6.00 4.00 1.00 5.00 Atomic Wigner-Seitz radii RWIGS = 0.73 1.11 0.32 0.75 virtual crystal weights VCA = 1.00 1.00 1.00 1.00 NELECT = 27.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.96E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 83.33 562.36 Fermi-wavevector in a.u.,A,eV,Ry = 0.491130 0.928100 3.281838 0.241208 Thomas-Fermi vector in A = 1.494349 Write flags LWAVE = F write WAVECAR LDOWNSAMPLE = F k-point downsampling of WAVECAR LCHARG = F write CHGCAR LVTOT = F write LOCPOT, total local potential LVHAR = F write LOCPOT, Hartree potential only LELF = F write electronic localiz. function (ELF) LORBIT = 0 0 simple, 1 ext, 2 COOP (PROOUT), +10 PAW based schemes Dipole corrections LMONO = F monopole corrections only (constant potential shift) LDIPOL = F correct potential (dipole corrections) IDIPOL = 0 1-x, 2-y, 3-z, 4-all directions EPSILON= 1.0000000 bulk dielectric constant Exchange correlation treatment: GGA = -- GGA type LEXCH = 8 internal setting for exchange type VOSKOWN= 1 Vosko Wilk Nusair interpolation LHFCALC = F Hartree Fock is set to LHFONE = F Hartree Fock one center treatment AEXX = 0.0000 exact exchange contribution Linear response parameters LEPSILON= F determine dielectric tensor LRPA = F only Hartree local field effects (RPA) LNABLA = F use nabla operator in PAW spheres LVEL = F velocity operator in full k-point grid LINTERFAST= F fast interpolation KINTER = 0 interpolate to denser k-point grid CSHIFT =0.1000 complex shift for real part using Kramers Kronig OMEGAMAX= -1.0 maximum frequency DEG_THRESHOLD= 0.2000000E-02 threshold for treating states as degnerate RTIME = -0.100 relaxation time in fs (WPLASMAI= 0.000 imaginary part of plasma frequency in eV, 0.658/RTIME) DFIELD = 0.0000000 0.0000000 0.0000000 field for delta impulse in time Orbital magnetization related: ORBITALMAG= F switch on orbital magnetization LCHIMAG = F perturbation theory with respect to B field DQ = 0.001000 dq finite difference perturbation B field LLRAUG = F two centre corrections for induced B field -------------------------------------------------------------------------------------------------------- Static calculation charge density and potential will be updated during run spin polarized calculation Variant of blocked Davidson Davidson routine will perform the subspace rotation perform sub-space diagonalisation after iterative eigenvector-optimisation modified Broyden-mixing scheme, WC = 100.0 initial mixing is a Kerker type mixing with AMIX = 0.4000 and BMIX = 1.0000 Hartree-type preconditioning will be used using additional bands 12 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 : 1000.00 direct lattice vectors reciprocal lattice vectors 10.000000000 0.000000000 0.000000000 0.100000000 0.000000000 0.000000000 0.000000000 10.000000000 0.000000000 0.000000000 0.100000000 0.000000000 0.000000000 0.000000000 10.000000000 0.000000000 0.000000000 0.100000000 length of vectors 10.000000000 10.000000000 10.000000000 0.100000000 0.100000000 0.100000000 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 0.125 0.05000000 0.00000000 0.00000000 0.125 0.00000000 0.05000000 0.00000000 0.125 0.00000000 0.00000000 0.05000000 0.125 0.05000000 0.05000000 0.00000000 0.125 0.00000000 0.05000000 0.05000000 0.125 0.05000000 0.00000000 0.05000000 0.125 0.05000000 0.05000000 0.05000000 0.125 k-points in reciprocal lattice and weights: Automatic mesh 0.00000000 0.00000000 0.00000000 0.125 0.50000000 0.00000000 0.00000000 0.125 0.00000000 0.50000000 0.00000000 0.125 0.00000000 0.00000000 0.50000000 0.125 0.50000000 0.50000000 0.00000000 0.125 0.00000000 0.50000000 0.50000000 0.125 0.50000000 0.00000000 0.50000000 0.125 0.50000000 0.50000000 0.50000000 0.125 position of ions in fractional coordinates (direct lattice) 0.26730156 0.04940302 0.05493753 0.06602838 0.15583947 0.52115927 0.39521175 0.14012594 0.72053663 0.33537187 0.18851098 0.55528787 0.91653883 0.09696315 0.38878005 0.36139141 0.94749374 0.40313961 0.43445557 0.46064704 0.72210889 0.45474333 0.38003241 0.59381734 0.67158556 0.31144534 0.83247346 0.82452506 0.84841622 0.84735976 0.78792827 0.59369990 0.33525082 0.72226835 0.26594233 0.52226672 position of ions in cartesian coordinates (Angst): 2.67301560 0.49403020 0.54937530 0.66028380 1.55839470 5.21159270 3.95211750 1.40125940 7.20536630 3.35371870 1.88510980 5.55287870 9.16538830 0.96963150 3.88780050 3.61391410 9.47493740 4.03139610 4.34455570 4.60647040 7.22108890 4.54743330 3.80032410 5.93817340 6.71585560 3.11445340 8.32473460 8.24525060 8.48416220 8.47359760 7.87928270 5.93699900 3.35250820 7.22268350 2.65942330 5.22266720 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 18037 k-point 2 : 0.5000 0.0000 0.0000 plane waves: 18184 k-point 3 : 0.0000 0.5000 0.0000 plane waves: 18184 k-point 4 : 0.0000 0.0000 0.5000 plane waves: 18184 k-point 5 : 0.5000 0.5000 0.0000 plane waves: 18176 k-point 6 : 0.0000 0.5000 0.5000 plane waves: 18176 k-point 7 : 0.5000 0.0000 0.5000 plane waves: 18176 k-point 8 : 0.5000 0.5000 0.5000 plane waves: 18048 maximum and minimum number of plane-waves per node : 18184 18037 maximum number of plane-waves: 18184 maximum index in each direction: IXMAX= 16 IYMAX= 16 IZMAX= 16 IXMIN= -16 IYMIN= -16 IZMIN= -16 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 70 to avoid them WARNING: aliasing errors must be expected set NGY to 70 to avoid them WARNING: aliasing errors must be expected set NGZ to 70 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 168555. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 1227. kBytes fftplans : 17620. kBytes grid : 58880. kBytes one-center: 73. kBytes wavefun : 60755. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 33 NGY = 33 NGZ = 33 (NGX =100 NGY =100 NGZ =100) gives a total of 35937 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 27.0000000 magnetization 12.0000000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 2280 Maximum index for augmentation-charges 2030 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.177 Maximum number of real-space cells 3x 3x 3 Maximum number of reciprocal cells 3x 3x 3 ----------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 848 total energy-change (2. order) : 0.1374167E+03 (-0.6045117E+03) number of electron 27.0000000 magnetization 12.0000000 augmentation part 27.0000000 magnetization 12.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -609.76393957 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 41.14871112 PAW double counting = 741.62526794 -709.47216318 entropy T*S EENTRO = -0.00228510 eigenvalues EBANDS = -102.43817702 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 137.41674827 eV energy without entropy = 137.41903337 energy(sigma->0) = 137.41750997 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 1304 total energy-change (2. order) :-0.1170354E+03 (-0.1096718E+03) number of electron 27.0000000 magnetization 12.0000000 augmentation part 27.0000000 magnetization 12.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -609.76393957 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 41.14871112 PAW double counting = 741.62526794 -709.47216318 entropy T*S EENTRO = 0.00251994 eigenvalues EBANDS = -219.47836656 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 20.38136378 eV energy without entropy = 20.37884383 energy(sigma->0) = 20.38052380 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 1026 total energy-change (2. order) :-0.4268069E+02 (-0.3827661E+02) number of electron 27.0000000 magnetization 12.0000000 augmentation part 27.0000000 magnetization 12.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -609.76393957 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 41.14871112 PAW double counting = 741.62526794 -709.47216318 entropy T*S EENTRO = 0.03297827 eigenvalues EBANDS = -262.18951650 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.29932783 eV energy without entropy = -22.33230610 energy(sigma->0) = -22.31032059 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 988 total energy-change (2. order) :-0.2569820E+01 (-0.2557717E+01) number of electron 27.0000000 magnetization 12.0000000 augmentation part 27.0000000 magnetization 12.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -609.76393957 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 41.14871112 PAW double counting = 741.62526794 -709.47216318 entropy T*S EENTRO = 0.01880181 eigenvalues EBANDS = -264.74516016 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.86914796 eV energy without entropy = -24.88794977 energy(sigma->0) = -24.87541523 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 1102 total energy-change (2. order) :-0.4583860E-01 (-0.4580133E-01) number of electron 26.9999966 magnetization 9.7355222 augmentation part 0.9576957 magnetization 2.8676005 Broyden mixing: rms(total) = 0.20679E+01 rms(broyden)= 0.20672E+01 rms(prec ) = 0.34698E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -609.76393957 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 41.14871112 PAW double counting = 741.62526794 -709.47216318 entropy T*S EENTRO = 0.01867375 eigenvalues EBANDS = -264.79087071 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.91498656 eV energy without entropy = -24.93366031 energy(sigma->0) = -24.92121114 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 948 total energy-change (2. order) :-0.2189718E+02 (-0.2769922E+02) number of electron 26.9999997 magnetization 9.0023896 augmentation part -1.0347725 magnetization 4.1273739 Broyden mixing: rms(total) = 0.39091E+01 rms(broyden)= 0.39073E+01 rms(prec ) = 0.65418E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3628 0.3628 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -690.56748297 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 54.16412781 PAW double counting = 914.00612067 -882.22749769 entropy T*S EENTRO = 0.00753405 eigenvalues EBANDS = -218.51430555 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -46.81216959 eV energy without entropy = -46.81970364 energy(sigma->0) = -46.81468094 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 1054 total energy-change (2. order) : 0.2812111E+02 (-0.1261610E+02) number of electron 26.9999975 magnetization 8.8988129 augmentation part 0.2343700 magnetization 1.9612217 Broyden mixing: rms(total) = 0.21619E+01 rms(broyden)= 0.21605E+01 rms(prec ) = 0.37103E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2846 0.3440 0.2252 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -636.79972715 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 50.11511134 PAW double counting = 837.01289170 -804.95144546 entropy T*S EENTRO = -0.01661681 eigenvalues EBANDS = -240.37060350 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -18.69105578 eV energy without entropy = -18.67443897 energy(sigma->0) = -18.68551684 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 932 total energy-change (2. order) : 0.4988824E+01 (-0.3488175E+00) number of electron 26.9999974 magnetization 8.8074136 augmentation part 0.2415288 magnetization 1.8718942 Broyden mixing: rms(total) = 0.18738E+01 rms(broyden)= 0.18737E+01 rms(prec ) = 0.32458E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3274 0.3691 0.3065 0.3065 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -637.35027860 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 49.34958531 PAW double counting = 838.17413408 -806.14456734 entropy T*S EENTRO = -0.04735519 eigenvalues EBANDS = -234.00308402 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -13.70223167 eV energy without entropy = -13.65487648 energy(sigma->0) = -13.68644661 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 962 total energy-change (2. order) : 0.1608180E+01 (-0.1935452E+00) number of electron 26.9999972 magnetization 8.7476908 augmentation part 0.3470054 magnetization 1.9170691 Broyden mixing: rms(total) = 0.14809E+01 rms(broyden)= 0.14806E+01 rms(prec ) = 0.24288E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3233 0.3497 0.3497 0.2970 0.2970 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -637.75691684 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 48.74363392 PAW double counting = 837.41575072 -805.37965677 entropy T*S EENTRO = -0.01345380 eigenvalues EBANDS = -231.42274285 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -12.09405153 eV energy without entropy = -12.08059773 energy(sigma->0) = -12.08956693 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 868 total energy-change (2. order) :-0.1250562E+01 (-0.3846069E-01) number of electron 26.9999972 magnetization 8.6856963 augmentation part 0.3062380 magnetization 1.8271895 Broyden mixing: rms(total) = 0.15311E+01 rms(broyden)= 0.15311E+01 rms(prec ) = 0.26172E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3809 0.4548 0.4548 0.3360 0.3360 0.3229 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -638.67915024 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 48.49960648 PAW double counting = 839.04776256 -807.00579842 entropy T*S EENTRO = -0.01525880 eigenvalues EBANDS = -231.51110910 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -13.34461344 eV energy without entropy = -13.32935464 energy(sigma->0) = -13.33952717 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 1022 total energy-change (2. order) :-0.3726959E+00 (-0.7838620E-02) number of electron 26.9999971 magnetization 8.5613445 augmentation part 0.3063763 magnetization 1.7070432 Broyden mixing: rms(total) = 0.14917E+01 rms(broyden)= 0.14917E+01 rms(prec ) = 0.25542E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5038 0.7116 0.7116 0.4141 0.4141 0.4457 0.3260 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -639.41479168 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 48.25148386 PAW double counting = 838.88060180 -806.83207583 entropy T*S EENTRO = -0.00678188 eigenvalues EBANDS = -230.91507974 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -13.71730938 eV energy without entropy = -13.71052750 energy(sigma->0) = -13.71504875 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 972 total energy-change (2. order) :-0.7953469E+00 (-0.2924708E-01) number of electron 26.9999971 magnetization 8.3104079 augmentation part 0.3015599 magnetization 1.5295788 Broyden mixing: rms(total) = 0.14211E+01 rms(broyden)= 0.14210E+01 rms(prec ) = 0.24566E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6558 1.1211 1.1211 0.4615 0.4615 0.5476 0.5476 0.3300 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -641.12512036 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 47.57301930 PAW double counting = 838.71713187 -806.67232264 entropy T*S EENTRO = 0.02622149 eigenvalues EBANDS = -229.35092005 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -14.51265632 eV energy without entropy = -14.53887781 energy(sigma->0) = -14.52139681 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 1016 total energy-change (2. order) :-0.1146915E+01 (-0.9986351E-01) number of electron 26.9999966 magnetization 8.3392517 augmentation part -0.2290118 magnetization -0.0326567 Broyden mixing: rms(total) = 0.15672E+01 rms(broyden)= 0.15643E+01 rms(prec ) = 0.23646E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5770 1.0654 1.0654 0.4587 0.4587 0.5694 0.5694 0.3230 0.1061 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -643.70806236 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 46.74032027 PAW double counting = 838.10842919 -806.04155727 entropy T*S EENTRO = -0.08137993 eigenvalues EBANDS = -226.99665543 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -15.65957144 eV energy without entropy = -15.57819152 energy(sigma->0) = -15.63244480 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 904 total energy-change (2. order) :-0.1600361E+01 (-0.1657025E+00) number of electron 26.9999972 magnetization 8.5061286 augmentation part 0.0324180 magnetization 1.2875633 Broyden mixing: rms(total) = 0.13156E+01 rms(broyden)= 0.13148E+01 rms(prec ) = 0.15483E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5632 1.1697 1.1697 0.4620 0.4620 0.5489 0.5489 0.3173 0.1949 0.1951 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -640.31808831 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.40720402 PAW double counting = 824.71941882 -792.56572613 entropy T*S EENTRO = -0.03160911 eigenvalues EBANDS = -230.79046618 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -17.25993281 eV energy without entropy = -17.22832369 energy(sigma->0) = -17.24939644 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 846 total energy-change (2. order) : 0.2106174E+00 (-0.1563783E+00) number of electron 26.9999973 magnetization 8.5586603 augmentation part -0.0606803 magnetization -1.4383324 Broyden mixing: rms(total) = 0.15855E+01 rms(broyden)= 0.15828E+01 rms(prec ) = 0.19696E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5192 1.1336 1.1336 0.4716 0.4716 0.5167 0.5167 0.3315 0.2436 0.2436 0.1295 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -639.81965914 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.16796821 PAW double counting = 825.05033350 -792.90739734 entropy T*S EENTRO = -0.01364004 eigenvalues EBANDS = -230.84625464 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -17.04931537 eV energy without entropy = -17.03567533 energy(sigma->0) = -17.04476869 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 860 total energy-change (2. order) :-0.2137643E+01 (-0.2663901E+00) number of electron 26.9999977 magnetization 8.5226763 augmentation part 0.1131034 magnetization -0.8084862 Broyden mixing: rms(total) = 0.10607E+01 rms(broyden)= 0.10600E+01 rms(prec ) = 0.12244E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5041 1.1027 1.1027 0.4731 0.4731 0.5198 0.5198 0.3416 0.3416 0.3267 0.2263 0.1180 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -639.82495965 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.28421538 PAW double counting = 821.93978093 -789.77179095 entropy T*S EENTRO = -0.01548445 eigenvalues EBANDS = -232.11805400 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -19.18695866 eV energy without entropy = -19.17147421 energy(sigma->0) = -19.18179718 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 1040 total energy-change (2. order) :-0.1127611E+01 (-0.1583029E+00) number of electron 26.9999976 magnetization 8.5243566 augmentation part -0.0955282 magnetization -0.9614437 Broyden mixing: rms(total) = 0.13289E+01 rms(broyden)= 0.13279E+01 rms(prec ) = 0.18194E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4755 1.1027 1.1027 0.4812 0.4812 0.5187 0.5187 0.3384 0.3384 0.3281 0.2464 0.1453 0.1040 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -641.58643993 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.13561393 PAW double counting = 823.83839756 -791.66781887 entropy T*S EENTRO = -0.01256415 eigenvalues EBANDS = -231.34109260 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.31456998 eV energy without entropy = -20.30200583 energy(sigma->0) = -20.31038193 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 1016 total energy-change (2. order) :-0.7651561E-01 (-0.1075904E-01) number of electron 26.9999978 magnetization 8.5355464 augmentation part -0.1104537 magnetization -0.8853535 Broyden mixing: rms(total) = 0.13341E+01 rms(broyden)= 0.13339E+01 rms(prec ) = 0.18670E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4852 1.1110 1.1110 0.5260 0.5260 0.4755 0.4755 0.3637 0.3637 0.3135 0.3135 0.3224 0.2901 0.1152 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -641.26798256 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.16536040 PAW double counting = 822.98311102 -790.80519309 entropy T*S EENTRO = -0.01310313 eigenvalues EBANDS = -231.77261231 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.39108559 eV energy without entropy = -20.37798246 energy(sigma->0) = -20.38671788 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 924 total energy-change (2. order) : 0.2439077E-02 (-0.7714146E-02) number of electron 26.9999979 magnetization 8.5362498 augmentation part -0.1115840 magnetization -0.5333525 Broyden mixing: rms(total) = 0.13115E+01 rms(broyden)= 0.13113E+01 rms(prec ) = 0.18593E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4694 1.1068 1.1068 0.5279 0.5279 0.4720 0.4720 0.3973 0.3973 0.3150 0.3150 0.3053 0.3053 0.2079 0.1146 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -640.89032513 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.18897521 PAW double counting = 821.84509021 -789.66300541 entropy T*S EENTRO = -0.00344069 eigenvalues EBANDS = -232.18527479 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.38864652 eV energy without entropy = -20.38520583 energy(sigma->0) = -20.38749962 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 1006 total energy-change (2. order) : 0.3955475E-01 (-0.3807757E-03) number of electron 26.9999978 magnetization 8.5498816 augmentation part -0.1112030 magnetization -0.5284776 Broyden mixing: rms(total) = 0.13122E+01 rms(broyden)= 0.13122E+01 rms(prec ) = 0.18604E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4517 1.1067 1.1067 0.4119 0.4119 0.4732 0.4732 0.5238 0.5238 0.3125 0.3125 0.3129 0.3129 0.2790 0.1148 0.0993 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -640.89039096 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.22184220 PAW double counting = 821.84920298 -789.66714418 entropy T*S EENTRO = -0.00348786 eigenvalues EBANDS = -232.17844802 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.34909176 eV energy without entropy = -20.34560390 energy(sigma->0) = -20.34792914 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 926 total energy-change (2. order) :-0.2397393E+00 (-0.9203588E-02) number of electron 26.9999979 magnetization 8.7188311 augmentation part -0.1246391 magnetization -0.5972611 Broyden mixing: rms(total) = 0.13337E+01 rms(broyden)= 0.13337E+01 rms(prec ) = 0.19186E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4620 1.1153 1.1153 0.4974 0.4974 0.4753 0.4753 0.5170 0.5170 0.3366 0.3273 0.3273 0.3132 0.3132 0.2248 0.2248 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -641.76380119 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.27695560 PAW double counting = 826.90793829 -794.73300507 entropy T*S EENTRO = -0.01593605 eigenvalues EBANDS = -231.58031678 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.58883111 eV energy without entropy = -20.57289506 energy(sigma->0) = -20.58351909 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 884 total energy-change (2. order) :-0.6883668E+00 (-0.5302593E-01) number of electron 26.9999980 magnetization 8.7473400 augmentation part -0.1246073 magnetization -1.7024089 Broyden mixing: rms(total) = 0.13359E+01 rms(broyden)= 0.13356E+01 rms(prec ) = 0.19312E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4506 1.1151 1.1151 0.5094 0.5094 0.4752 0.4752 0.5202 0.5202 0.3538 0.3305 0.3305 0.3164 0.3164 0.2523 0.2523 0.1149 0.1537 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -642.63345534 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.20382618 PAW double counting = 835.71857157 -803.57564975 entropy T*S EENTRO = 0.04880761 eigenvalues EBANDS = -231.35863230 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.27719794 eV energy without entropy = -21.32600555 energy(sigma->0) = -21.29346714 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 1084 total energy-change (2. order) :-0.1152888E-01 (-0.1759614E-02) number of electron 26.9999980 magnetization 8.6363118 augmentation part -0.1230952 magnetization -1.8652731 Broyden mixing: rms(total) = 0.13293E+01 rms(broyden)= 0.13293E+01 rms(prec ) = 0.19242E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5021 1.1074 1.1074 0.6807 0.6878 0.6878 0.4772 0.4772 0.5101 0.5101 0.4343 0.4343 0.3160 0.3160 0.3176 0.2841 0.2841 0.2910 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -642.43129573 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.16094075 PAW double counting = 835.20258655 -803.05797455 entropy T*S EENTRO = 0.04965384 eigenvalues EBANDS = -231.53197177 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.28872682 eV energy without entropy = -21.33838066 energy(sigma->0) = -21.30527810 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 956 total energy-change (2. order) :-0.9688820E-01 (-0.5912325E-02) number of electron 26.9999980 magnetization 8.0173749 augmentation part -0.1424000 magnetization -2.3512687 Broyden mixing: rms(total) = 0.13679E+01 rms(broyden)= 0.13678E+01 rms(prec ) = 0.20099E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5818 1.0541 1.1758 1.1758 0.9872 0.9872 0.6228 0.6228 0.4718 0.4718 0.2969 0.5534 0.5534 0.3157 0.3157 0.4252 0.2929 0.2929 0.3236 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -643.15931757 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.29836709 PAW double counting = 839.71798367 -807.56379481 entropy T*S EENTRO = 0.03272210 eigenvalues EBANDS = -231.03090958 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.38561502 eV energy without entropy = -21.41833712 energy(sigma->0) = -21.39652239 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 1008 total energy-change (2. order) :-0.2692911E+01 (-0.1140854E+01) number of electron 26.9999987 magnetization 7.9494496 augmentation part -0.3255607 magnetization -0.3780602 Broyden mixing: rms(total) = 0.18587E+01 rms(broyden)= 0.18550E+01 rms(prec ) = 0.29727E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5656 1.0512 1.1785 1.1785 0.9715 0.9715 0.6047 0.6047 0.2966 0.4700 0.4700 0.5257 0.5257 0.4077 0.4077 0.3176 0.3176 0.3250 0.2861 0.2861 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -647.45449661 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.98274493 PAW double counting = 880.53429162 -848.38874018 entropy T*S EENTRO = 0.01254949 eigenvalues EBANDS = -230.08420912 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.07852580 eV energy without entropy = -24.09107528 energy(sigma->0) = -24.08270896 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 1048 total energy-change (2. order) :-0.1879432E+00 (-0.5298016E-01) number of electron 26.9999988 magnetization 8.0304787 augmentation part -0.3027020 magnetization -0.6181568 Broyden mixing: rms(total) = 0.17845E+01 rms(broyden)= 0.17841E+01 rms(prec ) = 0.28510E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5950 1.0833 1.2832 1.2832 0.9809 0.9809 0.6020 0.6020 0.6102 0.6102 0.2963 0.4706 0.4706 0.5385 0.5385 0.4894 0.3174 0.3174 0.3231 0.2909 0.2909 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -647.19404143 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.80585228 PAW double counting = 885.42545450 -853.27257758 entropy T*S EENTRO = 0.04101520 eigenvalues EBANDS = -230.39150605 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.26646899 eV energy without entropy = -24.30748418 energy(sigma->0) = -24.28014072 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 934 total energy-change (2. order) :-0.2639828E-01 (-0.2118018E-01) number of electron 26.9999989 magnetization 8.0048714 augmentation part -0.3027972 magnetization -0.6350322 Broyden mixing: rms(total) = 0.17929E+01 rms(broyden)= 0.17928E+01 rms(prec ) = 0.28560E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5988 1.1039 1.3493 1.3493 1.0042 1.0042 0.6654 0.6654 0.6097 0.6097 0.2962 0.4713 0.4713 0.5431 0.5431 0.3171 0.3171 0.4177 0.4177 0.3250 0.2887 0.2887 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -648.23056860 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.78854908 PAW double counting = 896.41696244 -864.34012807 entropy T*S EENTRO = 0.04878320 eigenvalues EBANDS = -229.29579940 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.29286727 eV energy without entropy = -24.34165046 energy(sigma->0) = -24.30912833 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 1038 total energy-change (2. order) :-0.1807681E+00 (-0.1849878E-01) number of electron 26.9999990 magnetization 7.8963664 augmentation part -0.3033863 magnetization -0.7230804 Broyden mixing: rms(total) = 0.18221E+01 rms(broyden)= 0.18221E+01 rms(prec ) = 0.28841E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6350 1.1313 1.0584 1.0584 1.2149 1.2149 0.9104 0.9104 0.6844 0.6844 0.2962 0.4700 0.4700 0.6223 0.6223 0.5497 0.5497 0.5041 0.3174 0.3174 0.2901 0.2901 0.3236 0.1149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -649.41646825 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.04900797 PAW double counting = 914.99213557 -882.94838426 entropy T*S EENTRO = 0.04149276 eigenvalues EBANDS = -228.51075325 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -24.47363538 eV energy without entropy = -24.51512814 energy(sigma->0) = -24.48746630 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 882 total energy-change (2. order) : 0.6755546E+00 (-0.1338790E+00) number of electron 26.9999994 magnetization 7.8957166 augmentation part -0.3025166 magnetization -0.8115463 Broyden mixing: rms(total) = 0.18537E+01 rms(broyden)= 0.18533E+01 rms(prec ) = 0.29068E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6143 1.1345 1.0678 1.0678 1.2044 1.2044 0.9370 0.9370 0.6777 0.6777 0.2962 0.4699 0.4699 0.6131 0.6131 0.5499 0.5499 0.4967 0.3174 0.3174 0.3235 0.2901 0.2901 0.1149 0.1235 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -648.30392758 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.18175683 PAW double counting = 933.63397960 -901.61676765 entropy T*S EENTRO = 0.04062064 eigenvalues EBANDS = -229.05307668 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.79808074 eV energy without entropy = -23.83870139 energy(sigma->0) = -23.81162096 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 1042 total energy-change (2. order) :-0.3353593E-02 (-0.7237686E-02) number of electron 26.9999994 magnetization 7.8960023 augmentation part -0.2996372 magnetization -0.8152695 Broyden mixing: rms(total) = 0.18509E+01 rms(broyden)= 0.18509E+01 rms(prec ) = 0.29038E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5959 1.1365 1.2155 1.2155 1.0673 1.0673 0.9470 0.9470 0.6750 0.6750 0.2962 0.4699 0.4699 0.6094 0.6094 0.5485 0.5485 0.4963 0.3174 0.3174 0.3236 0.2901 0.2901 0.1149 0.1247 0.1247 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -648.33048190 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.20082371 PAW double counting = 933.67787176 -901.66048719 entropy T*S EENTRO = 0.04142757 eigenvalues EBANDS = -229.04992238 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.80143434 eV energy without entropy = -23.84286190 energy(sigma->0) = -23.81524353 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 888 total energy-change (2. order) : 0.1252317E+00 (-0.2018153E-02) number of electron 26.9999995 magnetization 7.8830456 augmentation part -0.2943095 magnetization -0.7714033 Broyden mixing: rms(total) = 0.18374E+01 rms(broyden)= 0.18373E+01 rms(prec ) = 0.28774E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5996 1.1430 1.0672 1.0672 1.1939 1.1939 0.8903 0.8903 0.6741 0.6741 0.2962 0.6209 0.6209 0.4701 0.4701 0.4463 0.4463 0.5465 0.5465 0.4967 0.3174 0.3174 0.3236 0.2901 0.2901 0.1149 0.1822 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -648.29540354 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.27417376 PAW double counting = 935.97279977 -903.95955697 entropy T*S EENTRO = 0.03037661 eigenvalues EBANDS = -229.01792631 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.67620260 eV energy without entropy = -23.70657921 energy(sigma->0) = -23.68632814 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 930 total energy-change (2. order) : 0.2257678E+00 (-0.6178044E-02) number of electron 26.9999994 magnetization 7.9542569 augmentation part -0.2571979 magnetization -0.3414389 Broyden mixing: rms(total) = 0.17004E+01 rms(broyden)= 0.17001E+01 rms(prec ) = 0.26681E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6113 1.1488 1.0599 1.0599 1.0746 1.0746 0.9534 0.9534 0.7053 0.7053 0.7192 0.7192 0.2962 0.6329 0.6329 0.4701 0.4701 0.5507 0.5507 0.1149 0.4848 0.3174 0.3174 0.3239 0.2885 0.2885 0.2964 0.2964 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -647.76704442 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.30143716 PAW double counting = 932.69547866 -900.66047391 entropy T*S EENTRO = -0.01760897 eigenvalues EBANDS = -229.32155742 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.45043482 eV energy without entropy = -23.43282585 energy(sigma->0) = -23.44456516 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 956 total energy-change (2. order) : 0.2965312E+00 (-0.2359506E-01) number of electron 26.9999987 magnetization 7.9979052 augmentation part -0.1139920 magnetization 1.0436207 Broyden mixing: rms(total) = 0.12618E+01 rms(broyden)= 0.12594E+01 rms(prec ) = 0.19374E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6073 1.1405 1.2116 1.2116 1.0265 1.0265 1.0748 1.0748 0.6838 0.6838 0.4616 0.4616 0.2962 0.6384 0.6384 0.4701 0.4701 0.5328 0.5328 0.4930 0.4930 0.4847 0.3174 0.3174 0.1149 0.3234 0.2901 0.2901 0.2452 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -645.98938241 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.92262273 PAW double counting = 920.09923868 -888.05637907 entropy T*S EENTRO = -0.03972465 eigenvalues EBANDS = -230.40961301 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.15390365 eV energy without entropy = -23.11417900 energy(sigma->0) = -23.14066210 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 922 total energy-change (2. order) : 0.7351345E-01 (-0.1902810E-01) number of electron 26.9999985 magnetization 7.9845347 augmentation part -0.0901509 magnetization 1.2250211 Broyden mixing: rms(total) = 0.12572E+01 rms(broyden)= 0.12570E+01 rms(prec ) = 0.18813E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6084 1.1325 1.3132 1.3132 1.0127 1.0127 1.0259 1.0259 0.6867 0.6867 0.4462 0.4462 0.2962 0.6243 0.6243 0.6152 0.6152 0.4701 0.4701 0.5455 0.5455 0.1149 0.3174 0.3174 0.4738 0.2901 0.2901 0.3481 0.3230 0.2597 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.75524534 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.24903504 PAW double counting = 941.94731138 -909.96475736 entropy T*S EENTRO = -0.02911093 eigenvalues EBANDS = -229.84695709 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.08039021 eV energy without entropy = -23.05127927 energy(sigma->0) = -23.07068656 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 1080 total energy-change (2. order) : 0.4178772E-01 (-0.7409599E-02) number of electron 26.9999985 magnetization 7.9880324 augmentation part -0.0860320 magnetization 1.2331831 Broyden mixing: rms(total) = 0.12650E+01 rms(broyden)= 0.12650E+01 rms(prec ) = 0.18831E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5892 1.1327 1.3122 1.3122 1.0130 1.0130 1.0257 1.0257 0.6871 0.6871 0.4479 0.4479 0.2962 0.6259 0.6259 0.6150 0.6150 0.4701 0.4701 0.5457 0.5457 0.1149 0.3174 0.3174 0.4736 0.2901 0.2901 0.3548 0.3231 0.2586 0.0225 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -647.08747570 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.40539792 PAW double counting = 950.98118775 -919.02237834 entropy T*S EENTRO = -0.02646903 eigenvalues EBANDS = -229.60819918 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.03860249 eV energy without entropy = -23.01213346 energy(sigma->0) = -23.02977948 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 1030 total energy-change (2. order) :-0.2356222E-01 (-0.1001663E-02) number of electron 26.9999986 magnetization 7.9214178 augmentation part -0.0958092 magnetization 1.0724960 Broyden mixing: rms(total) = 0.12876E+01 rms(broyden)= 0.12876E+01 rms(prec ) = 0.19270E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6258 1.1168 1.2379 1.2379 0.7836 1.0228 1.0228 0.9221 0.9221 0.9191 0.9191 0.6851 0.6851 0.6672 0.6672 0.4701 0.4701 0.5512 0.5512 0.5494 0.5494 0.4720 0.4720 0.3174 0.3174 0.3234 0.2901 0.2901 0.1149 0.2962 0.3053 0.2515 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -647.17217659 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.42522909 PAW double counting = 950.21133995 -918.24672947 entropy T*S EENTRO = -0.03163624 eigenvalues EBANDS = -229.56752554 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23.06216471 eV energy without entropy = -23.03052847 energy(sigma->0) = -23.05161930 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 950 total energy-change (2. order) : 0.1343067E+00 (-0.7485102E-02) number of electron 26.9999985 magnetization 7.8674824 augmentation part -0.0769149 magnetization 1.1855642 Broyden mixing: rms(total) = 0.12289E+01 rms(broyden)= 0.12289E+01 rms(prec ) = 0.18330E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6473 1.1039 0.9220 1.1488 1.1488 1.1360 1.1360 1.0414 1.0414 1.0051 1.0051 0.7038 0.7038 0.7507 0.7507 0.4701 0.4701 0.5830 0.5830 0.5591 0.5591 0.3174 0.3174 0.4700 0.4613 0.4613 0.3234 0.2901 0.2901 0.1149 0.2962 0.2970 0.2512 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.65422481 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.38621144 PAW double counting = 943.45512352 -911.46922400 entropy T*S EENTRO = -0.02508072 eigenvalues EBANDS = -229.93999749 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.92785798 eV energy without entropy = -22.90277726 energy(sigma->0) = -22.91949774 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 896 total energy-change (2. order) : 0.1679645E-01 (-0.1902384E-02) number of electron 26.9999985 magnetization 7.9634065 augmentation part -0.0775564 magnetization 1.2663879 Broyden mixing: rms(total) = 0.12133E+01 rms(broyden)= 0.12132E+01 rms(prec ) = 0.18231E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6636 1.4218 1.4218 1.1108 0.8874 0.8874 1.0039 1.0039 0.9995 0.9995 0.9737 0.9737 0.7114 0.7114 0.7883 0.7883 0.4701 0.4701 0.5754 0.5754 0.5735 0.5735 0.3174 0.3174 0.5015 0.5015 0.4802 0.3234 0.2901 0.2901 0.1149 0.2962 0.2947 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.25742167 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.31881652 PAW double counting = 936.04666276 -904.04067327 entropy T*S EENTRO = -0.02873036 eigenvalues EBANDS = -230.26904959 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.91106153 eV energy without entropy = -22.88233117 energy(sigma->0) = -22.90148475 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 904 total energy-change (2. order) :-0.3748361E-01 (-0.5400562E-02) number of electron 26.9999985 magnetization 8.0791465 augmentation part -0.0820861 magnetization 1.3522540 Broyden mixing: rms(total) = 0.12462E+01 rms(broyden)= 0.12461E+01 rms(prec ) = 0.18556E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6854 1.6024 1.6024 1.0826 1.0765 1.0765 1.0131 1.0131 1.0638 1.0638 0.9781 0.9781 0.7137 0.7137 0.8047 0.8047 0.4701 0.4701 0.5814 0.5814 0.5715 0.5715 0.3174 0.3174 0.5090 0.5090 0.4789 0.4789 0.3234 0.2901 0.2901 0.1149 0.2962 0.2952 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.86144993 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.40097919 PAW double counting = 943.93035677 -911.94794008 entropy T*S EENTRO = -0.03224525 eigenvalues EBANDS = -229.75757992 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.94854514 eV energy without entropy = -22.91629989 energy(sigma->0) = -22.93779673 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 1078 total energy-change (2. order) : 0.4148083E-01 (-0.5530295E-02) number of electron 26.9999986 magnetization 8.1173032 augmentation part -0.0886116 magnetization 1.2798489 Broyden mixing: rms(total) = 0.12576E+01 rms(broyden)= 0.12575E+01 rms(prec ) = 0.18723E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6705 1.6024 1.6024 1.0913 1.0913 1.0716 1.0142 1.0142 1.0756 1.0756 0.9828 0.9828 0.7134 0.7134 0.7951 0.7951 0.4701 0.4701 0.5851 0.5851 0.1149 0.5591 0.5591 0.3174 0.3174 0.4891 0.4891 0.4749 0.4491 0.3234 0.2901 0.2901 0.2962 0.2952 0.2514 0.2215 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.83544415 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.31150596 PAW double counting = 943.35564176 -911.37717881 entropy T*S EENTRO = -0.04434361 eigenvalues EBANDS = -229.63657955 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.90706431 eV energy without entropy = -22.86272071 energy(sigma->0) = -22.89228311 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 1042 total energy-change (2. order) : 0.2502743E-01 (-0.3947083E-02) number of electron 26.9999986 magnetization 8.0346901 augmentation part -0.0771946 magnetization 1.2377664 Broyden mixing: rms(total) = 0.12241E+01 rms(broyden)= 0.12241E+01 rms(prec ) = 0.18109E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7271 1.6278 1.6278 1.5119 1.5119 1.1260 1.1260 0.9880 1.0138 1.0138 1.0066 1.0066 0.8859 0.8859 0.7183 0.7183 0.7707 0.7707 0.4701 0.4701 0.5704 0.5704 0.5979 0.5979 0.3174 0.3174 0.5322 0.5322 0.5359 0.4921 0.3234 0.2901 0.2901 0.1149 0.2962 0.2952 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -646.20751350 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.12818174 PAW double counting = 937.76723266 -905.77904199 entropy T*S EENTRO = -0.04083978 eigenvalues EBANDS = -230.06939010 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.88203689 eV energy without entropy = -22.84119711 energy(sigma->0) = -22.86842363 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 898 total energy-change (2. order) : 0.2587940E+00 (-0.3087474E-01) number of electron 26.9999985 magnetization 7.9970738 augmentation part -0.0755790 magnetization 1.2611499 Broyden mixing: rms(total) = 0.11779E+01 rms(broyden)= 0.11776E+01 rms(prec ) = 0.17838E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7071 1.6217 1.6217 1.5011 1.5011 1.1204 1.1204 0.9877 1.0139 1.0139 1.0070 1.0070 0.8844 0.8844 0.7183 0.7183 0.7701 0.7701 0.4701 0.4701 0.5701 0.5701 0.5976 0.5976 0.3174 0.3174 0.5324 0.5324 0.5365 0.4921 0.3234 0.2901 0.2901 0.1149 0.2962 0.2952 0.2513 0.0375 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -644.15039636 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.97763636 PAW double counting = 928.85834264 -896.81149468 entropy T*S EENTRO = -0.03854530 eigenvalues EBANDS = -231.77811960 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.62324287 eV energy without entropy = -22.58469757 energy(sigma->0) = -22.61039444 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 888 total energy-change (2. order) : 0.1712487E-01 (-0.4972584E-02) number of electron 26.9999985 magnetization 7.9997696 augmentation part -0.0742016 magnetization 1.2731106 Broyden mixing: rms(total) = 0.11714E+01 rms(broyden)= 0.11713E+01 rms(prec ) = 0.17804E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7268 1.6424 1.6424 1.4587 1.4587 1.1481 1.1213 1.1213 0.9960 1.0141 1.0141 1.0014 1.0014 0.8941 0.8941 0.7178 0.7178 0.7727 0.7727 0.4701 0.4701 0.2962 0.2952 0.1149 0.5698 0.5698 0.6009 0.6009 0.3174 0.3174 0.5327 0.5327 0.5247 0.4926 0.3234 0.2901 0.2901 0.3692 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -644.45534422 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.07839635 PAW double counting = 927.02505132 -894.97600095 entropy T*S EENTRO = -0.03944507 eigenvalues EBANDS = -231.55810952 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.60611800 eV energy without entropy = -22.56667293 energy(sigma->0) = -22.59296964 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 1032 total energy-change (2. order) : 0.4306739E-01 (-0.4972660E-03) number of electron 26.9999985 magnetization 8.2185361 augmentation part -0.0742472 magnetization 1.4876454 Broyden mixing: rms(total) = 0.11753E+01 rms(broyden)= 0.11753E+01 rms(prec ) = 0.17832E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7847 1.9551 1.9551 1.9253 1.9253 1.2606 1.2606 0.9898 1.0122 1.0122 1.0949 1.0949 1.0032 1.0032 0.9172 0.9172 0.7170 0.7170 0.7548 0.7548 0.4701 0.4701 0.5661 0.5661 0.5907 0.5907 0.3174 0.3174 0.5348 0.5348 0.5133 0.5133 0.4883 0.3234 0.2901 0.2901 0.1149 0.2962 0.2952 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -644.41358273 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.12568395 PAW double counting = 928.74739613 -896.69732457 entropy T*S EENTRO = -0.03879765 eigenvalues EBANDS = -231.60575982 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -22.56305061 eV energy without entropy = -22.52425296 energy(sigma->0) = -22.55011806 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 904 total energy-change (2. order) : 0.1156124E+01 (-0.2863833E+00) number of electron 26.9999984 magnetization 8.3341465 augmentation part -0.0593679 magnetization 1.2465495 Broyden mixing: rms(total) = 0.11486E+01 rms(broyden)= 0.11485E+01 rms(prec ) = 0.16931E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7868 2.0080 2.0080 2.0451 2.0451 1.2729 1.2729 0.9901 1.0124 1.0124 1.0838 1.0838 1.0286 1.0286 0.9401 0.9401 0.7171 0.7171 0.7804 0.7804 0.4701 0.4701 0.2962 0.2952 0.1149 0.5676 0.5676 0.5963 0.5963 0.3174 0.3174 0.5274 0.5274 0.5228 0.4889 0.3234 0.2901 0.2901 0.4383 0.4383 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -640.73009510 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.94738135 PAW double counting = 923.69589367 -891.60257429 entropy T*S EENTRO = -0.05595219 eigenvalues EBANDS = -233.98091433 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.40692680 eV energy without entropy = -21.35097462 energy(sigma->0) = -21.38827607 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 1092 total energy-change (2. order) : 0.5776145E+00 (-0.5972950E-01) number of electron 26.9999984 magnetization 8.1316605 augmentation part 0.0068308 magnetization 1.0910762 Broyden mixing: rms(total) = 0.10967E+01 rms(broyden)= 0.10964E+01 rms(prec ) = 0.15619E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7701 2.0693 2.0693 1.8516 1.8516 1.2691 1.2691 0.9897 1.0124 1.0124 1.1310 1.1310 0.9037 0.9037 0.9675 0.9675 0.7172 0.7172 0.7499 0.7499 0.2962 0.2952 0.1149 0.4701 0.4701 0.5645 0.5645 0.6115 0.6115 0.3174 0.3174 0.5429 0.5429 0.5226 0.5226 0.4881 0.3234 0.2901 0.2901 0.4170 0.4170 0.2513 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -638.96110685 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.82811374 PAW double counting = 918.30782206 -886.18739017 entropy T*S EENTRO = -0.09386249 eigenvalues EBANDS = -235.04222266 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.82931229 eV energy without entropy = -20.73544979 energy(sigma->0) = -20.79802479 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 842 total energy-change (2. order) :-0.2829443E+00 (-0.9524348E-01) number of electron 26.9999985 magnetization 8.3785333 augmentation part 0.0667382 magnetization 2.5488920 Broyden mixing: rms(total) = 0.13903E+01 rms(broyden)= 0.13888E+01 rms(prec ) = 0.19889E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7764 2.0548 2.0548 2.0621 2.0621 1.2757 1.2757 0.9904 1.0126 1.0126 1.1005 1.1005 1.0557 1.0557 0.8601 0.8601 0.7173 0.7173 0.8200 0.8200 0.4701 0.4701 0.6256 0.6256 0.5700 0.5700 0.3174 0.3174 0.5738 0.5738 0.5348 0.5348 0.5331 0.5331 0.4830 0.3234 0.2901 0.2901 0.2962 0.2952 0.1149 0.2513 0.1065 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -639.29177984 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.88502403 PAW double counting = 914.91704847 -882.79280590 entropy T*S EENTRO = -0.03628660 eigenvalues EBANDS = -235.11279081 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.11225656 eV energy without entropy = -21.07596996 energy(sigma->0) = -21.10016102 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 854 total energy-change (2. order) : 0.5222733E+00 (-0.2366415E+00) number of electron 26.9999976 magnetization 8.4163127 augmentation part 0.2091045 magnetization 1.2260172 Broyden mixing: rms(total) = 0.88598E+00 rms(broyden)= 0.87894E+00 rms(prec ) = 0.10768E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7698 2.0424 2.0424 2.0888 2.0888 1.2772 1.2772 0.9903 1.0125 1.0125 1.2048 1.2048 1.0679 1.0679 0.9134 0.9134 0.7172 0.7172 0.8186 0.8186 0.4701 0.4701 0.6376 0.6376 0.5695 0.5695 0.3174 0.3174 0.5579 0.5579 0.5527 0.5527 0.5281 0.5281 0.4855 0.3234 0.2901 0.2901 0.2962 0.2952 0.1149 0.2513 0.1055 0.1055 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -636.56925840 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.72690536 PAW double counting = 907.93587542 -875.78171030 entropy T*S EENTRO = -0.04674880 eigenvalues EBANDS = -237.17438058 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.58998322 eV energy without entropy = -20.54323443 energy(sigma->0) = -20.57440029 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 924 total energy-change (2. order) :-0.4844723E-03 (-0.1270131E-01) number of electron 26.9999977 magnetization 8.4300767 augmentation part 0.2163845 magnetization 1.4651644 Broyden mixing: rms(total) = 0.94149E+00 rms(broyden)= 0.94114E+00 rms(prec ) = 0.11631E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7524 2.0422 2.0422 2.0886 2.0886 1.2773 1.2773 0.9903 1.0125 1.0125 1.2054 1.2054 1.0675 1.0675 0.9135 0.9135 0.7172 0.7172 0.8188 0.8188 0.4701 0.4701 0.6377 0.6377 0.5694 0.5694 0.5579 0.5579 0.5527 0.5527 0.5282 0.5282 0.4855 0.3174 0.3174 0.3234 0.2901 0.2901 0.2962 0.2952 0.2513 0.1149 0.1070 0.1070 0.0040 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -636.43526214 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.79408376 PAW double counting = 916.66881555 -884.53179966 entropy T*S EENTRO = -0.04207095 eigenvalues EBANDS = -237.36356832 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.59046770 eV energy without entropy = -20.54839674 energy(sigma->0) = -20.57644404 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 998 total energy-change (2. order) :-0.1667384E-01 (-0.2222472E-02) number of electron 26.9999977 magnetization 9.3853754 augmentation part 0.2128749 magnetization 2.4073140 Broyden mixing: rms(total) = 0.93121E+00 rms(broyden)= 0.93117E+00 rms(prec ) = 0.11473E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7372 1.9476 1.9476 1.9161 1.9161 1.0190 1.1580 1.1580 1.2808 1.2808 0.8486 0.8486 1.0337 0.7603 0.7603 0.3184 0.2600 0.2381 0.8107 0.8107 0.1407 0.1047 0.7686 0.7686 0.7172 0.6373 0.6373 0.3996 0.3996 0.4644 0.4644 0.1045 0.2538 0.2538 0.5551 0.5551 0.3846 0.3846 0.5352 0.4653 0.1801 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -636.48434349 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.78756248 PAW double counting = 917.63223946 -885.49738902 entropy T*S EENTRO = -0.04357499 eigenvalues EBANDS = -237.32097006 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.60714153 eV energy without entropy = -20.56356654 energy(sigma->0) = -20.59261654 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 902 total energy-change (2. order) : 0.7291342E+01 (-0.2730102E+01) number of electron 26.9999963 magnetization 9.5129907 augmentation part 0.4365770 magnetization 2.6359154 Broyden mixing: rms(total) = 0.12837E+01 rms(broyden)= 0.12803E+01 rms(prec ) = 0.16057E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7161 1.9744 1.9744 1.6530 1.6530 1.0185 1.1603 1.1603 1.3233 1.3233 0.7854 0.7854 0.8578 0.8578 1.0429 0.8343 0.8343 0.3163 0.2721 0.2325 0.1408 0.1037 0.7789 0.7789 0.3313 0.3313 0.0922 0.4823 0.4823 0.6164 0.6164 0.2691 0.2691 0.6290 0.4156 0.4156 0.5647 0.5647 0.5509 0.4817 0.1930 0.1946 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -629.24194410 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 48.36108310 PAW double counting = 902.43040216 -870.10968921 entropy T*S EENTRO = -0.05785198 eigenvalues EBANDS = -241.01713364 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -13.31579958 eV energy without entropy = -13.25794760 energy(sigma->0) = -13.29651558 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- eigenvalue-minimisations : 886 total energy-change (2. order) :-0.5773628E+01 (-0.6106867E+00) number of electron 26.9999963 magnetization 9.6650671 augmentation part 0.4025066 magnetization 1.5432938 Broyden mixing: rms(total) = 0.10303E+01 rms(broyden)= 0.10286E+01 rms(prec ) = 0.13203E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7032 1.9976 1.9976 1.6485 1.6485 1.0187 1.1378 1.1378 1.3050 1.3050 0.7745 0.7745 0.8438 0.8438 1.0425 0.8297 0.8297 0.3163 0.2720 0.2319 0.7960 0.7960 0.3374 0.3374 0.1358 0.1005 0.4881 0.4881 0.5977 0.5977 0.6111 0.6111 0.5432 0.5432 0.4445 0.4445 0.4859 0.0843 0.2543 0.2543 0.1894 0.2199 0.2199 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -628.67592620 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.25175607 PAW double counting = 862.35995254 -830.04387614 entropy T*S EENTRO = -0.00866141 eigenvalues EBANDS = -244.29200671 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -19.08942777 eV energy without entropy = -19.08076636 energy(sigma->0) = -19.08654063 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 53) --------------------------------------- eigenvalue-minimisations : 958 total energy-change (2. order) : 0.2564123E-01 (-0.5099285E-01) number of electron 26.9999966 magnetization 9.4517231 augmentation part 0.3961972 magnetization 1.2474385 Broyden mixing: rms(total) = 0.10386E+01 rms(broyden)= 0.10384E+01 rms(prec ) = 0.13260E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6994 1.9741 1.9741 1.6567 1.6567 1.0187 1.1788 1.1788 1.3053 1.3053 0.8142 0.8142 1.0431 0.8519 0.8519 0.8567 0.8567 0.8403 0.3262 0.2344 0.2344 0.2082 0.3682 0.3682 0.1281 0.1281 0.6711 0.6711 0.4834 0.4834 0.5920 0.5920 0.0923 0.2477 0.2477 0.6193 0.5407 0.5407 0.4939 0.4066 0.4066 0.4125 0.1974 0.2011 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -629.77506718 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 45.21228003 PAW double counting = 862.27406745 -829.98135242 entropy T*S EENTRO = -0.00182079 eigenvalues EBANDS = -243.11122771 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -19.06378654 eV energy without entropy = -19.06196575 energy(sigma->0) = -19.06317961 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 54) --------------------------------------- eigenvalue-minimisations : 972 total energy-change (2. order) :-0.1080106E+01 (-0.5492432E-01) number of electron 26.9999962 magnetization 9.3146303 augmentation part 0.3996874 magnetization 1.1743549 Broyden mixing: rms(total) = 0.98837E+00 rms(broyden)= 0.98827E+00 rms(prec ) = 0.12914E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6913 1.9846 1.9846 1.5947 1.5947 1.0175 1.1892 1.1892 1.3208 1.3208 0.7946 0.7946 1.0597 0.8605 0.8605 0.8520 0.8520 0.8767 0.3209 0.2629 0.2629 0.2519 0.3938 0.3938 0.1401 0.0989 0.6483 0.6483 0.4861 0.4861 0.3152 0.3152 0.6240 0.5746 0.5746 0.0967 0.5332 0.5332 0.4440 0.4440 0.4180 0.3709 0.2562 0.2061 0.1698 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -628.75783089 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.78595755 PAW double counting = 855.11488476 -822.80933555 entropy T*S EENTRO = -0.00967704 eigenvalues EBANDS = -244.78722559 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.14389266 eV energy without entropy = -20.13421562 energy(sigma->0) = -20.14066698 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 55) --------------------------------------- eigenvalue-minimisations : 882 total energy-change (2. order) :-0.3398447E+00 (-0.1871514E-01) number of electron 26.9999960 magnetization 8.0596842 augmentation part 0.4004871 magnetization -0.0927541 Broyden mixing: rms(total) = 0.97408E+00 rms(broyden)= 0.97404E+00 rms(prec ) = 0.12787E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6592 1.7462 1.7462 1.3482 1.1080 1.1080 1.0502 1.0502 1.1288 1.1288 0.5215 0.5215 0.9697 0.9697 0.8756 0.8756 0.9141 0.5295 0.5295 0.7758 0.7758 0.2522 0.1689 0.1689 0.7244 0.0550 0.4913 0.4913 0.2480 0.2480 0.0264 0.0814 0.5221 0.5221 0.5616 0.5388 0.4290 0.4290 0.2005 0.2005 0.3342 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -629.01911923 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.58913266 PAW double counting = 849.47049731 -817.15898928 entropy T*S EENTRO = -0.01031559 eigenvalues EBANDS = -244.67427727 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -20.48373732 eV energy without entropy = -20.47342173 energy(sigma->0) = -20.48029879 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 56) --------------------------------------- eigenvalue-minimisations : 850 total energy-change (2. order) : 0.1113744E+01 (-0.7586224E+00) number of electron 26.9999943 magnetization 8.3812035 augmentation part 0.4089246 magnetization 0.7074183 Broyden mixing: rms(total) = 0.12095E+01 rms(broyden)= 0.12052E+01 rms(prec ) = 0.14525E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6528 1.3377 1.7369 1.7369 1.1879 1.1879 1.0333 1.0333 0.9956 0.9956 1.1411 1.1411 0.5105 0.5105 0.7807 0.7807 0.9176 0.9176 0.5936 0.5936 0.2509 0.1761 0.1761 0.6914 0.6914 0.0677 0.1280 0.1280 0.5012 0.5012 0.0810 0.2152 0.2152 0.2079 0.2079 0.5358 0.5358 0.5501 0.5501 0.3470 0.4378 0.4378 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -624.99839663 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 44.61460262 PAW double counting = 797.34685681 -764.79690776 entropy T*S EENTRO = -0.04242339 eigenvalues EBANDS = -247.81305889 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -19.36999316 eV energy without entropy = -19.32756977 energy(sigma->0) = -19.35585203 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 57) --------------------------------------- eigenvalue-minimisations : 886 total energy-change (2. order) :-0.2280400E+01 (-0.2892144E+00) number of electron 26.9999949 magnetization 8.3965390 augmentation part 0.4681061 magnetization 0.4543363 Broyden mixing: rms(total) = 0.12540E+01 rms(broyden)= 0.12533E+01 rms(prec ) = 0.16521E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6437 1.3620 1.7730 1.7730 1.1601 1.1601 0.9680 0.9680 0.9957 0.9957 1.1099 1.1099 0.5781 0.5781 0.8205 0.8205 0.6631 0.6631 0.9124 0.9124 0.6798 0.6798 0.2520 0.1759 0.1759 0.1644 0.1644 0.5032 0.5032 0.5431 0.5431 0.5465 0.5465 0.0668 0.4379 0.4379 0.0631 0.0631 0.3520 0.1863 0.1863 0.2203 0.2203 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 2.18778779 Ewald energy TEWEN = -228.56720593 -Hartree energ DENC = -625.03773603 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 43.33590803 PAW double counting = 793.07133068 -760.57116280 entropy T*S EENTRO = -0.03129458 eigenvalues EBANDS = -248.73677289 atomic energy EATOM = 1002.69875223 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -21.65039351 eV energy without entropy = -21.61909893 energy(sigma->0) = -21.63996198 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 58) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 7 1 8 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------