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 19:42:27 running on 3 total cores distrk: each k-point on 3 cores, 1 groups distr: one band on NCORE= 1 cores, 3 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 200 ALGO = Normal (blocked Davidson) ISPIN = 2 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 1.11 0.32 0.73 NPAR = 3 POTCAR: PAW_PBE Si 05Jan2001 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE O 08Apr2002 ----------------------------------------------------------------------------- | | | W W AA RRRRR N N II N N GGGG !!! | | W W A A R R NN N II NN N G G !!! | | W W A A R R N N N II N N N G !!! | | W WW W AAAAAA RRRRR N N N II N N N G GGG ! | | WW WW A A R R N NN II N NN G G | | W W A A R R N N II N N GGGG !!! | | | | You use a magnetic or noncollinear calculation, but did not specify | | the initial magnetic moment with the MAGMOM tag. Note that a | | default of 1 will be used for all atoms. This ferromagnetic setup | | may break the symmetry of the crystal, in particular it may rule | | out finding an antiferromagnetic solution. Thence, we recommend | | setting the initial magnetic moment manually or verifying carefully | | that this magnetic setup is desired. | | | ----------------------------------------------------------------------------- POTCAR: PAW_PBE Si 05Jan2001 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 POTCAR: PAW_PBE O 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in kinetic energy density of atom read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 Optimization of the real space projectors (new method) maximal supplied QI-value = 19.84 optimisation between [QCUT,QGAM] = [ 10.12, 20.44] = [ 28.68,116.96] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 7 10.119 159.560 0.56E-04 0.22E-03 0.45E-07 0 7 10.119 115.863 0.56E-04 0.21E-03 0.45E-07 1 7 10.119 88.339 0.34E-03 0.49E-03 0.11E-06 1 7 10.119 48.592 0.33E-03 0.48E-03 0.11E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 24.76 optimisation between [QCUT,QGAM] = [ 10.15, 20.30] = [ 28.85,115.39] Ry Optimized for a Real-space Cutoff 1.38 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.150 20.381 0.22E-03 0.32E-03 0.29E-06 0 8 10.150 15.268 0.23E-03 0.35E-03 0.30E-06 1 8 10.150 5.964 0.46E-03 0.53E-03 0.21E-06 1 8 10.150 5.382 0.38E-03 0.45E-03 0.19E-06 PAW_PBE Si 05Jan2001 : energy of atom 1 EATOM= -103.0669 kinetic energy error for atom= 0.0012 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 2 EATOM= -12.4884 kinetic energy error for atom= 0.0098 (will be added to EATOM!!) PAW_PBE O 08Apr2002 : energy of atom 3 EATOM= -432.3788 kinetic energy error for atom= 0.1156 (will be added to EATOM!!) POSCAR: No title positions in direct lattice No initial velocities read in exchange correlation table for LEXCH = 8 RHO(1)= 0.500 N(1) = 2000 RHO(2)= 100.500 N(2) = 4000 -------------------------------------------------------------------------------------------------------- ion position nearest neighbor table 1 0.930 0.528 0.109- 19 2.31 14 2.39 3 2.42 2 0.851 0.233 0.427- 4 2.34 13 2.38 20 2.43 3 0.137 0.448 0.195- 8 2.31 17 2.35 7 2.38 1 2.42 4 0.095 0.282 0.333- 18 2.29 2 2.34 7 2.35 8 2.66 5 0.817 0.334 0.026- 23 2.34 14 2.34 7 2.35 10 2.37 6 0.957 0.421 0.518- 8 2.33 13 2.37 24 2.40 11 2.41 9 2.63 7 0.082 0.334 0.137- 21 2.34 4 2.35 5 2.35 3 2.38 8 0.200 0.407 0.390- 3 2.31 6 2.33 22 2.38 18 2.61 4 2.66 9 0.831 0.546 0.521- 11 2.36 27 2.36 6 2.63 10 0.772 0.244 0.883- 5 2.37 28 2.39 12 2.40 31 2.45 11 0.060 0.511 0.652- 16 2.36 9 2.36 25 2.41 6 2.41 12 0.013 0.206 0.762- 26 2.35 10 2.40 15 2.53 13 0.902 0.310 0.592- 31 2.32 6 2.37 2 2.38 15 2.40 14 0.865 0.445 0.956- 16 2.33 5 2.34 32 2.35 1 2.39 15 0.138 0.324 0.737- 29 2.34 16 2.38 13 2.40 26 2.50 12 2.53 16 0.097 0.435 0.818- 14 2.33 30 2.35 11 2.36 15 2.38 17 0.423 0.461 0.121- 19 2.30 30 2.33 3 2.35 23 2.50 18 0.364 0.293 0.421- 4 2.29 29 2.33 20 2.35 8 2.61 19 0.689 0.485 0.209- 17 2.30 1 2.31 24 2.34 20 0.643 0.313 0.340- 23 2.34 18 2.35 24 2.42 2 2.43 21 0.326 0.318 0.010- 26 2.34 7 2.34 23 2.37 30 2.38 22 0.435 0.438 0.522- 24 2.38 8 2.38 29 2.45 25 2.48 27 2.48 23 0.561 0.346 0.141- 5 2.34 20 2.34 21 2.37 17 2.50 24 0.692 0.430 0.400- 19 2.34 22 2.38 6 2.40 20 2.42 25 0.318 0.555 0.551- 35 1.78 11 2.41 22 2.48 27 2.56 26 0.274 0.232 0.864- 21 2.34 12 2.35 28 2.41 15 2.50 27 0.601 0.517 0.657- 35 1.67 32 2.35 9 2.36 22 2.48 25 2.56 28 0.538 0.210 0.751- 10 2.39 26 2.41 31 2.52 29 0.398 0.328 0.625- 18 2.33 15 2.34 31 2.38 22 2.45 30 0.368 0.427 0.919- 32 2.32 17 2.33 16 2.35 21 2.38 31 0.672 0.326 0.728- 13 2.32 29 2.38 32 2.38 10 2.45 28 2.52 32 0.632 0.437 0.815- 30 2.32 14 2.35 27 2.35 31 2.38 33 0.656 0.697 0.624- 34 0.72 34 0.664 0.710 0.562- 33 0.72 35 0.455 0.580 0.676- 27 1.67 25 1.78 LATTYP: Found a simple orthorhombic cell. ALAT = 7.6631000000 B/A-ratio = 1.4142135689 C/A-ratio = 2.5659458966 Lattice vectors: A1 = ( -7.6631000000, 0.0000000000, 0.0000000000) A2 = ( 0.0000000000, 0.0000000000, 10.8372600000) A3 = ( 0.0000000000, 19.6631000000, 0.0000000000) Analysis of symmetry for initial positions (statically): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The static configuration has the point symmetry C_1 . Analysis of symmetry for dynamics (positions and initial velocities): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The dynamic configuration has the point symmetry C_1 . Analysis of structural, dynamic, and magnetic symmetry: ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The overall configuration has the point symmetry C_1 . Subroutine INISYM returns: Found 1 space group operations (whereof 1 operations are pure point group operations), and found 1 'primitive' translations ---------------------------------------------------------------------------------------- Primitive cell volume of cell : 1632.9616 direct lattice vectors reciprocal lattice vectors 7.663100000 0.000000000 0.000000000 0.130495491 0.000000000 0.000000000 0.000000000 19.663100000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 0.000000000 10.837260000 0.000000000 0.000000000 0.092274246 length of vectors 7.663100000 19.663100000 10.837260000 0.130495491 0.050856681 0.092274246 position of ions in fractional coordinates (direct lattice) 0.930273680 0.528472380 0.108600600 0.851080040 0.233005280 0.427363270 0.136661500 0.448255330 0.195492910 0.094539390 0.282363280 0.332534300 0.817445110 0.333761340 0.026019990 0.956912860 0.421497320 0.517673160 0.081617350 0.333600100 0.136688110 0.199611300 0.407448980 0.390138540 0.830548460 0.545798590 0.520696650 0.771507560 0.244397730 0.882937460 0.060020920 0.510800300 0.651762660 0.013386240 0.205548650 0.761962260 0.902427990 0.310245310 0.592383940 0.864581110 0.444663970 0.955624190 0.137646260 0.323713840 0.736697570 0.097296190 0.435195240 0.818403980 0.423014390 0.460648700 0.121456720 0.364182950 0.292608490 0.421130370 0.689474280 0.484557350 0.208916490 0.643301950 0.312683830 0.340176050 0.325784310 0.317890090 0.009657850 0.435198130 0.438198720 0.522123960 0.561489560 0.346275930 0.141243590 0.692070810 0.429536060 0.400453220 0.317631120 0.554523220 0.551468980 0.274315110 0.231907380 0.864485220 0.600555990 0.517333300 0.656614700 0.538423170 0.209894200 0.751154120 0.397689950 0.327909950 0.624552190 0.368420620 0.427144510 0.918973560 0.671679940 0.326489620 0.727596970 0.632210700 0.436576870 0.815201770 0.655633280 0.696605040 0.623554100 0.663646170 0.709507430 0.562101090 0.455181910 0.579613800 0.676051730 ion indices of the primitive-cell ions primitive index ion index 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31 32 32 33 33 34 34 35 35 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 2 1 2 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.065247746 0.000000000 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.046137123 0.000000000 0.000000000 0.500000000 Length of vectors 0.065247746 0.050856681 0.046137123 Shift w.r.t. Gamma in fractional coordinates (k-lattice) 0.000000000 0.000000000 0.000000000 Subroutine IBZKPT returns following result: =========================================== Found 4 irreducible k-points: Following reciprocal coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.000000 0.500000 1.000000 0.500000 0.000000 0.500000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.065248 0.000000 0.000000 1.000000 0.000000 0.000000 0.046137 1.000000 0.065248 0.000000 0.046137 1.000000 -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 4 k-points in BZ NKDIM = 4 number of bands NBANDS= 105 number of dos NEDOS = 301 number of ions NIONS = 35 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 207360 max r-space proj IRMAX = 1449 max aug-charges IRDMAX= 4445 dimension x,y,z NGX = 40 NGY = 96 NGZ = 54 dimension x,y,z NGXF= 80 NGYF= 192 NGZF= 108 support grid NGXF= 80 NGYF= 192 NGZF= 108 ions per type = 32 2 1 NGX,Y,Z is equivalent to a cutoff of 8.68, 8.12, 8.28 a.u. NGXF,Y,Z is equivalent to a cutoff of 17.36, 16.23, 16.57 a.u. SYSTEM = No title POSCAR = No title Startparameter for this run: NWRITE = 1 write-flag & timer PREC = normal normal or accurate (medium, high low for compatibility) ISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ISPIN = 2 spin polarized calculation? LNONCOLLINEAR = F non collinear calculations LSORBIT = F spin-orbit coupling INIWAV = 1 electr: 0-lowe 1-rand 2-diag LASPH = F aspherical Exc in radial PAW Electronic Relaxation 1 ENCUT = 400.0 eV 29.40 Ry 5.42 a.u. 12.50 32.07 17.67*2*pi/ulx,y,z ENINI = 400.0 initial cutoff ENAUG = 605.4 eV augmentation charge cutoff NELM = 200; NELMIN= 2; NELMDL= -5 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM LREAL = T real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 1 Vosko Wilk Nusair interpolation ROPT = -0.00050 -0.00050 -0.00050 Ionic relaxation EDIFFG = -.2E-01 stopping-criterion for IOM NSW = 0 number of steps for IOM NBLOCK = 1; KBLOCK = 1 inner block; outer block IBRION = -1 ionic relax: 0-MD 1-quasi-New 2-CG NFREE = 0 steps in history (QN), initial steepest desc. (CG) ISIF = 2 stress and relaxation IWAVPR = 10 prediction: 0-non 1-charg 2-wave 3-comb ISYM = 2 0-nonsym 1-usesym 2-fastsym LCORR = T Harris-Foulkes like correction to forces POTIM = 0.5000 time-step for ionic-motion TEIN = 0.0 initial temperature TEBEG = 0.0; TEEND = 0.0 temperature during run SMASS = -3.00 Nose mass-parameter (am) estimated Nose-frequenzy (Omega) = 0.10E-29 period in steps = 0.13E+47 mass= -0.134E-26a.u. SCALEE = 1.0000 scale energy and forces NPACO = 256; APACO = 16.0 distance and # of slots for P.C. PSTRESS= 0.0 pullay stress Mass of Ions in am POMASS = 28.09 1.00 16.00 Ionic Valenz ZVAL = 4.00 1.00 6.00 Atomic Wigner-Seitz radii RWIGS = 1.11 0.32 0.73 virtual crystal weights VCA = 1.00 1.00 1.00 NELECT = 136.0000 total number of electrons NUPDOWN= -1.0000 fix difference up-down DOS related values: EMIN = 10.00; EMAX =-10.00 energy-range for DOS EFERMI = 0.00 ISMEAR = 1; SIGMA = 0.20 broadening in eV -4-tet -1-fermi 0-gaus Electronic relaxation 2 (details) IALGO = 38 algorithm LDIAG = T sub-space diagonalisation (order eigenvalues) LSUBROT= F optimize rotation matrix (better conditioning) TURBO = 0 0=normal 1=particle mesh IRESTART = 0 0=no restart 2=restart with 2 vectors NREBOOT = 0 no. of reboots NMIN = 0 reboot dimension EREF = 0.00 reference energy to select bands IMIX = 4 mixing-type and parameters AMIX = 0.40; BMIX = 1.00 AMIX_MAG = 1.60; BMIX_MAG = 1.00 AMIN = 0.10 WC = 100.; INIMIX= 1; MIXPRE= 1; MAXMIX= -45 Intra band minimization: WEIMIN = 0.0000 energy-eigenvalue tresh-hold EBREAK = 0.24E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 46.66 314.85 Fermi-wavevector in a.u.,A,eV,Ry = 0.714928 1.351019 6.954243 0.511122 Thomas-Fermi vector in A = 1.802956 Write flags LWAVE = F write WAVECAR LDOWNSAMPLE = F k-point downsampling of WAVECAR LCHARG = F write CHGCAR LVTOT = F write LOCPOT, total local potential LVHAR = F write LOCPOT, Hartree potential only LELF = F write electronic localiz. function (ELF) LORBIT = 0 0 simple, 1 ext, 2 COOP (PROOUT), +10 PAW based schemes Dipole corrections LMONO = F monopole corrections only (constant potential shift) LDIPOL = F correct potential (dipole corrections) IDIPOL = 0 1-x, 2-y, 3-z, 4-all directions EPSILON= 1.0000000 bulk dielectric constant Exchange correlation treatment: GGA = -- GGA type LEXCH = 8 internal setting for exchange type VOSKOWN= 1 Vosko Wilk Nusair interpolation LHFCALC = F Hartree Fock is set to LHFONE = F Hartree Fock one center treatment AEXX = 0.0000 exact exchange contribution Linear response parameters LEPSILON= F determine dielectric tensor LRPA = F only Hartree local field effects (RPA) LNABLA = F use nabla operator in PAW spheres LVEL = F velocity operator in full k-point grid LINTERFAST= F fast interpolation KINTER = 0 interpolate to denser k-point grid CSHIFT =0.1000 complex shift for real part using Kramers Kronig OMEGAMAX= -1.0 maximum frequency DEG_THRESHOLD= 0.2000000E-02 threshold for treating states as degnerate RTIME = -0.100 relaxation time in fs (WPLASMAI= 0.000 imaginary part of plasma frequency in eV, 0.658/RTIME) DFIELD = 0.0000000 0.0000000 0.0000000 field for delta impulse in time Orbital magnetization related: ORBITALMAG= F switch on orbital magnetization LCHIMAG = F perturbation theory with respect to B field DQ = 0.001000 dq finite difference perturbation B field LLRAUG = F two centre corrections for induced B field -------------------------------------------------------------------------------------------------------- Static calculation charge density and potential will be updated during run spin polarized calculation Variant of blocked Davidson Davidson routine will perform the subspace rotation perform sub-space diagonalisation after iterative eigenvector-optimisation modified Broyden-mixing scheme, WC = 100.0 initial mixing is a Kerker type mixing with AMIX = 0.4000 and BMIX = 1.0000 Hartree-type preconditioning will be used using additional bands 37 real space projection scheme for non local part use partial core corrections calculate Harris-corrections to forces (improved forces if not selfconsistent) use gradient corrections use of overlap-Matrix (Vanderbilt PP) Methfessel and Paxton Order N= 1 SIGMA = 0.20 -------------------------------------------------------------------------------------------------------- energy-cutoff : 400.00 volume of cell : 1632.96 direct lattice vectors reciprocal lattice vectors 7.663100000 0.000000000 0.000000000 0.130495491 0.000000000 0.000000000 0.000000000 19.663100000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 0.000000000 10.837260000 0.000000000 0.000000000 0.092274246 length of vectors 7.663100000 19.663100000 10.837260000 0.130495491 0.050856681 0.092274246 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 0.250 0.06524775 0.00000000 0.00000000 0.250 0.00000000 0.00000000 0.04613712 0.250 0.06524775 0.00000000 0.04613712 0.250 k-points in reciprocal lattice and weights: Automatic mesh 0.00000000 0.00000000 0.00000000 0.250 0.50000000 0.00000000 0.00000000 0.250 0.00000000 0.00000000 0.50000000 0.250 0.50000000 0.00000000 0.50000000 0.250 position of ions in fractional coordinates (direct lattice) 0.93027368 0.52847238 0.10860060 0.85108004 0.23300528 0.42736327 0.13666150 0.44825533 0.19549291 0.09453939 0.28236328 0.33253430 0.81744511 0.33376134 0.02601999 0.95691286 0.42149732 0.51767316 0.08161735 0.33360010 0.13668811 0.19961130 0.40744898 0.39013854 0.83054846 0.54579859 0.52069665 0.77150756 0.24439773 0.88293746 0.06002092 0.51080030 0.65176266 0.01338624 0.20554865 0.76196226 0.90242799 0.31024531 0.59238394 0.86458111 0.44466397 0.95562419 0.13764626 0.32371384 0.73669757 0.09729619 0.43519524 0.81840398 0.42301439 0.46064870 0.12145672 0.36418295 0.29260849 0.42113037 0.68947428 0.48455735 0.20891649 0.64330195 0.31268383 0.34017605 0.32578431 0.31789009 0.00965785 0.43519813 0.43819872 0.52212396 0.56148956 0.34627593 0.14124359 0.69207081 0.42953606 0.40045322 0.31763112 0.55452322 0.55146898 0.27431511 0.23190738 0.86448522 0.60055599 0.51733330 0.65661470 0.53842317 0.20989420 0.75115412 0.39768995 0.32790995 0.62455219 0.36842062 0.42714451 0.91897356 0.67167994 0.32648962 0.72759697 0.63221070 0.43657687 0.81520177 0.65563328 0.69660504 0.62355410 0.66364617 0.70950743 0.56210109 0.45518191 0.57961380 0.67605173 position of ions in cartesian coordinates (Angst): 7.12878024 10.39140526 1.17693294 6.52191145 4.58160612 4.63144687 1.04725074 8.81408938 2.11860749 0.72446480 5.55213741 3.60376067 6.26416362 6.56278260 0.28198540 7.33291894 8.28794395 5.61015863 0.62544191 6.55961213 1.48132459 1.52964135 8.01171004 4.22803279 6.36457590 10.73209226 5.64292498 5.91213958 4.80561700 9.56862282 0.45994631 10.04391738 7.06332140 0.10258010 4.04172366 8.25758312 6.91539593 6.10038456 6.41981878 6.62537150 8.74347211 10.35634781 1.05479706 6.36521761 7.98378311 0.74559043 8.55728752 8.86925672 3.24160157 9.05778145 1.31625805 2.79077036 5.75359000 4.56389931 5.28351036 9.52789963 2.26408232 4.92968717 6.14833342 3.68657630 2.49651775 6.25070463 0.10466463 3.33496679 8.61634525 5.65839311 4.30275065 6.80885824 1.53069351 5.30340782 8.44601050 4.33981566 2.43403904 10.90364553 5.97641272 2.10210412 4.56001800 9.36865110 4.60212061 10.17237641 7.11590422 4.12599059 4.12717064 8.14045250 3.04753786 6.44772614 6.76843447 2.82324405 8.39898521 9.95915540 5.14715055 6.41979805 7.88515754 4.84469382 8.58445465 8.83455353 5.02418339 13.69741456 6.75761791 5.08558697 13.95111555 6.09163566 3.48810449 11.39700411 7.32654837 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 29735 k-point 2 : 0.5000 0.0000 0.0000 plane waves: 29568 k-point 3 : 0.0000 0.0000 0.5000 plane waves: 29692 k-point 4 : 0.5000 0.0000 0.5000 plane waves: 29632 maximum and minimum number of plane-waves per node : 29735 29568 maximum number of plane-waves: 29735 maximum index in each direction: IXMAX= 12 IYMAX= 32 IZMAX= 17 IXMIN= -12 IYMIN= -32 IZMIN= -18 The following grids will avoid any aliasing or wrap around errors in the Hartre e energy - symmetry arguments have not been applied - exchange correlation energies might require even more grid points - we recommend to set PREC=Normal or Accurate and rely on VASP defaults WARNING: aliasing errors must be expected set NGX to 50 to avoid them WARNING: aliasing errors must be expected set NGY to 140 to avoid them WARNING: aliasing errors must be expected set NGZ to 72 to avoid them serial 3D FFT for wavefunctions parallel 3D FFT for charge: minimum data exchange during FFTs selected (reduces bandwidth) total amount of memory used by VASP MPI-rank0 254498. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3166. kBytes fftplans : 19938. kBytes grid : 66726. kBytes one-center: 215. kBytes wavefun : 134453. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 25 NGY = 65 NGZ = 35 (NGX = 80 NGY =192 NGZ =108) gives a total of 56875 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 136.0000000 magnetization 35.0000000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 1356 Maximum index for augmentation-charges 1467 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.151 Maximum number of real-space cells 4x 2x 3 Maximum number of reciprocal cells 2x 4x 3 ----------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 1680 total energy-change (2. order) : 0.5988269E+03 (-0.3891358E+04) number of electron 136.0000000 magnetization 35.0000000 augmentation part 136.0000000 magnetization 35.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12279.45473997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.16209756 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00764930 eigenvalues EBANDS = -197.01325871 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 598.82694486 eV energy without entropy = 598.83459415 energy(sigma->0) = 598.82949462 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 2253 total energy-change (2. order) :-0.6875908E+03 (-0.6532013E+03) number of electron 136.0000000 magnetization 35.0000000 augmentation part 136.0000000 magnetization 35.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12279.45473997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.16209756 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00202307 eigenvalues EBANDS = -884.60964537 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -88.76381557 eV energy without entropy = -88.76179250 energy(sigma->0) = -88.76314122 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2109 total energy-change (2. order) :-0.7343822E+02 (-0.7164734E+02) number of electron 136.0000000 magnetization 35.0000000 augmentation part 136.0000000 magnetization 35.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12279.45473997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.16209756 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02579104 eigenvalues EBANDS = -958.02410174 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -162.20203991 eV energy without entropy = -162.17624888 energy(sigma->0) = -162.19344290 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 2199 total energy-change (2. order) :-0.3039521E+01 (-0.3020558E+01) number of electron 136.0000000 magnetization 35.0000000 augmentation part 136.0000000 magnetization 35.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12279.45473997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.16209756 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02688956 eigenvalues EBANDS = -961.06252415 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.24156085 eV energy without entropy = -165.21467128 energy(sigma->0) = -165.23259766 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 2406 total energy-change (2. order) :-0.1544509E+00 (-0.1543657E+00) number of electron 136.0000031 magnetization 30.2909572 augmentation part -6.9878317 magnetization 27.1084142 Broyden mixing: rms(total) = 0.25340E+01 rms(broyden)= 0.25338E+01 rms(prec ) = 0.26712E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12279.45473997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.16209756 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02690776 eigenvalues EBANDS = -961.21695688 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.39601178 eV energy without entropy = -165.36910401 energy(sigma->0) = -165.38704252 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 1890 total energy-change (2. order) : 0.4477052E+01 (-0.1568910E+02) number of electron 136.0000028 magnetization 27.4364136 augmentation part -6.9915020 magnetization 23.9880181 Broyden mixing: rms(total) = 0.21351E+01 rms(broyden)= 0.21327E+01 rms(prec ) = 0.21984E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5206 0.5206 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12386.97183411 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.92024739 PAW double counting = 6200.89536062 -5641.45338544 entropy T*S EENTRO = -0.01037488 eigenvalues EBANDS = -865.55288852 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -160.91895983 eV energy without entropy = -160.90858495 energy(sigma->0) = -160.91550154 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 2346 total energy-change (2. order) :-0.1582920E+04 (-0.1470229E+04) number of electron 136.0000045 magnetization 27.3222029 augmentation part -6.9824377 magnetization 24.3908260 Broyden mixing: rms(total) = 0.74446E+01 rms(broyden)= 0.73402E+01 rms(prec ) = 0.76305E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3726 0.6526 0.0926 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12417.55554229 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -239.76681798 PAW double counting = 8015.08131221 -7453.98991609 entropy T*S EENTRO = -0.00643745 eigenvalues EBANDS = -2425.69552752 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1743.83851923 eV energy without entropy = -1743.83208178 energy(sigma->0) = -1743.83637341 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 2205 total energy-change (2. order) : 0.1612362E+04 (-0.4193052E+03) number of electron 136.0000023 magnetization 22.5954836 augmentation part -7.0467403 magnetization 20.8557423 Broyden mixing: rms(total) = 0.17669E+01 rms(broyden)= 0.14663E+01 rms(prec ) = 0.14979E+01 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5386 1.1148 0.4690 0.0320 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12416.43900242 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -236.62490398 PAW double counting = 8116.53900001 -7557.51625997 entropy T*S EENTRO = 0.00810876 eigenvalues EBANDS = -815.53801088 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -131.47665859 eV energy without entropy = -131.48476736 energy(sigma->0) = -131.47936151 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 2145 total energy-change (2. order) :-0.3955854E+01 (-0.2333608E+01) number of electron 136.0000022 magnetization 18.8605798 augmentation part -7.0369991 magnetization 16.4318987 Broyden mixing: rms(total) = 0.98771E+00 rms(broyden)= 0.96772E+00 rms(prec ) = 0.99031E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6996 1.5781 0.5937 0.5937 0.0330 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12458.06779251 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -236.21440736 PAW double counting = 11039.45526242 -10480.93844490 entropy T*S EENTRO = -0.02352558 eigenvalues EBANDS = -777.73801477 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -135.43251282 eV energy without entropy = -135.40898724 energy(sigma->0) = -135.42467096 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 2091 total energy-change (2. order) :-0.8065270E+01 (-0.4906650E+00) number of electron 136.0000022 magnetization 14.1714866 augmentation part -6.9860708 magnetization 12.2667881 Broyden mixing: rms(total) = 0.74680E+00 rms(broyden)= 0.74577E+00 rms(prec ) = 0.75662E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8627 1.9719 0.9028 0.9028 0.5028 0.0330 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.85836072 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -240.04696792 PAW double counting = 12435.18174260 -11877.29386797 entropy T*S EENTRO = -0.02669510 eigenvalues EBANDS = -749.54804391 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -143.49778314 eV energy without entropy = -143.47108804 energy(sigma->0) = -143.48888477 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 2019 total energy-change (2. order) :-0.8033362E+01 (-0.4304504E+00) number of electron 136.0000023 magnetization 10.5842761 augmentation part -6.9564361 magnetization 9.4807557 Broyden mixing: rms(total) = 0.55419E+00 rms(broyden)= 0.55410E+00 rms(prec ) = 0.55940E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9626 2.5775 1.1858 0.0330 0.8518 0.5636 0.5636 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12519.16078884 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -242.79560070 PAW double counting = 13213.51702880 -12656.48394146 entropy T*S EENTRO = -0.00190841 eigenvalues EBANDS = -724.70034401 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -151.53114473 eV energy without entropy = -151.52923632 energy(sigma->0) = -151.53050859 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 1836 total energy-change (2. order) : 0.2576141E+01 (-0.1626506E+02) number of electron 136.0000022 magnetization 10.0769722 augmentation part -6.9682412 magnetization 9.2301518 Broyden mixing: rms(total) = 0.78580E+00 rms(broyden)= 0.78558E+00 rms(prec ) = 0.84227E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8719 2.5802 1.2535 0.9380 0.5363 0.5363 0.0330 0.2262 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12532.04013575 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.69000457 PAW double counting = 13399.84539737 -12847.20063573 entropy T*S EENTRO = -0.01681991 eigenvalues EBANDS = -701.94721540 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -148.95500410 eV energy without entropy = -148.93818420 energy(sigma->0) = -148.94939747 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 2658 total energy-change (2. order) :-0.4370431E+04 (-0.4104210E+04) number of electron 136.0000036 magnetization 10.0409691 augmentation part -6.8497048 magnetization 10.3284831 Broyden mixing: rms(total) = 0.70010E+01 rms(broyden)= 0.68864E+01 rms(prec ) = 0.71655E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7659 2.5940 1.2340 0.9526 0.5333 0.5333 0.0330 0.2318 0.0150 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12532.27369782 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -246.46405927 PAW double counting = 13372.54787004 -12814.98321829 entropy T*S EENTRO = 0.01373468 eigenvalues EBANDS = -5076.32079293 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -4519.38575370 eV energy without entropy = -4519.39948838 energy(sigma->0) = -4519.39033193 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 1908 total energy-change (2. order) : 0.4377562E+04 (-0.1915552E+04) number of electron 136.0000023 magnetization 9.7178129 augmentation part -6.8881251 magnetization 9.5626631 Broyden mixing: rms(total) = 0.12581E+01 rms(broyden)= 0.75269E+00 rms(prec ) = 0.77568E+00 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6834 2.6391 1.0751 1.0410 0.5397 0.5397 0.0330 0.1989 0.0692 0.0152 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12532.37172353 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -242.22466232 PAW double counting = 13393.93628500 -12836.12475998 entropy T*S EENTRO = 0.00623517 eigenvalues EBANDS = -703.13994177 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -141.82415753 eV energy without entropy = -141.83039270 energy(sigma->0) = -141.82623592 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 2667 total energy-change (2. order) :-0.6880791E+04 (-0.6446142E+04) number of electron 136.0000034 magnetization 9.5494322 augmentation part -6.7464751 magnetization 2.4511019 Broyden mixing: rms(total) = 0.66128E+01 rms(broyden)= 0.64457E+01 rms(prec ) = 0.67444E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6367 2.7455 1.0267 0.9166 0.5947 0.5947 0.3126 0.0330 0.0641 0.0641 0.0153 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12532.38680604 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.99056991 PAW double counting = 13727.43083169 -13169.64836343 entropy T*S EENTRO = 0.01171517 eigenvalues EBANDS = -7580.12631620 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -7022.61509883 eV energy without entropy = -7022.62681400 energy(sigma->0) = -7022.61900388 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 2130 total energy-change (2. order) : 0.6880453E+04 (-0.4431727E+03) number of electron 136.0000023 magnetization 9.0256000 augmentation part -6.8180690 magnetization 4.9444544 Broyden mixing: rms(total) = 0.13728E+01 rms(broyden)= 0.72194E+00 rms(prec ) = 0.74031E+00 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6118 2.9643 1.1614 0.6566 0.6566 0.5250 0.3034 0.3034 0.0330 0.0644 0.0467 0.0148 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12530.73423777 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -242.21414161 PAW double counting = 13778.03034948 -13220.17958452 entropy T*S EENTRO = 0.02110968 eigenvalues EBANDS = -705.18037655 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -142.16247141 eV energy without entropy = -142.18358109 energy(sigma->0) = -142.16950797 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 2301 total energy-change (2. order) :-0.1561034E+02 (-0.6742278E+01) number of electron 136.0000023 magnetization 5.9501551 augmentation part -6.9084511 magnetization 5.5526782 Broyden mixing: rms(total) = 0.46236E+00 rms(broyden)= 0.39731E+00 rms(prec ) = 0.40458E+00 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6978 3.5568 1.0845 0.8930 0.7245 0.7245 0.4696 0.3817 0.3817 0.0330 0.0550 0.0550 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12530.51230117 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -247.51911255 PAW double counting = 14275.67274076 -13718.17254897 entropy T*S EENTRO = 0.00312676 eigenvalues EBANDS = -715.33912577 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -157.77281104 eV energy without entropy = -157.77593780 energy(sigma->0) = -157.77385330 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 2022 total energy-change (2. order) :-0.1452876E+01 (-0.2997204E+02) number of electron 136.0000023 magnetization 5.1019064 augmentation part -6.9176769 magnetization 5.1180299 Broyden mixing: rms(total) = 0.49660E+00 rms(broyden)= 0.49403E+00 rms(prec ) = 0.50897E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6955 3.5666 1.3449 0.8118 0.8118 0.6496 0.6496 0.5143 0.2681 0.2681 0.0330 0.0545 0.0545 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12525.56617275 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -251.06437519 PAW double counting = 14519.34151095 -13961.38714713 entropy T*S EENTRO = 0.01361169 eigenvalues EBANDS = -718.65752491 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -159.22568745 eV energy without entropy = -159.23929914 energy(sigma->0) = -159.23022468 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 2082 total energy-change (2. order) :-0.7838686E+01 (-0.4551670E+01) number of electron 136.0000024 magnetization 3.7385457 augmentation part -6.9113901 magnetization 3.8010211 Broyden mixing: rms(total) = 0.44666E+00 rms(broyden)= 0.44642E+00 rms(prec ) = 0.45281E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8335 4.6039 2.3338 0.9876 0.9876 0.5545 0.5545 0.4653 0.4653 0.2794 0.2794 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12518.84077785 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -252.78314056 PAW double counting = 14352.97814709 -13794.75417782 entropy T*S EENTRO = 0.02190948 eigenvalues EBANDS = -731.78074397 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -167.06437375 eV energy without entropy = -167.08628322 energy(sigma->0) = -167.07167690 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 2046 total energy-change (2. order) :-0.3508497E+01 (-0.3532732E+00) number of electron 136.0000024 magnetization 2.5489758 augmentation part -6.9105100 magnetization 2.5813744 Broyden mixing: rms(total) = 0.45939E+00 rms(broyden)= 0.45936E+00 rms(prec ) = 0.46845E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8584 5.0190 2.4010 0.9315 0.5643 0.5643 0.7465 0.7465 0.6907 0.5065 0.2741 0.2741 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12508.16806827 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.02045711 PAW double counting = 14139.84833159 -13581.34783436 entropy T*S EENTRO = 0.02951877 eigenvalues EBANDS = -744.00877086 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -170.57287035 eV energy without entropy = -170.60238912 energy(sigma->0) = -170.58270994 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 2001 total energy-change (2. order) :-0.1168704E+01 (-0.1019608E+00) number of electron 136.0000024 magnetization 1.9789613 augmentation part -6.9115849 magnetization 1.9642047 Broyden mixing: rms(total) = 0.37441E+00 rms(broyden)= 0.37440E+00 rms(prec ) = 0.38530E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9060 5.5008 2.4730 0.9653 0.9653 0.8828 0.8828 0.5571 0.5571 0.5013 0.5013 0.2765 0.2765 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12506.92033147 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.13841989 PAW double counting = 13916.33905141 -13357.74937234 entropy T*S EENTRO = 0.02332005 eigenvalues EBANDS = -745.39023244 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.74157479 eV energy without entropy = -171.76489484 energy(sigma->0) = -171.74934814 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 2046 total energy-change (2. order) :-0.1388446E+01 (-0.7007136E-01) number of electron 136.0000024 magnetization 1.3184985 augmentation part -6.9096752 magnetization 1.2847765 Broyden mixing: rms(total) = 0.36178E+00 rms(broyden)= 0.36178E+00 rms(prec ) = 0.37621E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9742 5.9247 2.5817 1.2863 1.2863 1.0464 1.0464 0.5083 0.5083 0.5199 0.5722 0.5722 0.2756 0.2756 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12503.56184476 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.83530090 PAW double counting = 13627.08593398 -13068.52375384 entropy T*S EENTRO = 0.02398914 eigenvalues EBANDS = -749.41345414 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.13002064 eV energy without entropy = -173.15400978 energy(sigma->0) = -173.13801702 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 2004 total energy-change (2. order) :-0.1842808E+00 (-0.2390025E-01) number of electron 136.0000024 magnetization 0.4802291 augmentation part -6.9067932 magnetization 0.4413943 Broyden mixing: rms(total) = 0.28642E+00 rms(broyden)= 0.28642E+00 rms(prec ) = 0.30346E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0312 6.6537 2.5113 1.5304 1.5304 1.3723 0.5326 0.5326 0.7569 0.6402 0.6402 0.6057 0.5475 0.2758 0.2758 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12497.76426273 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.46292572 PAW double counting = 13182.64510309 -12624.06219570 entropy T*S EENTRO = 0.02599953 eigenvalues EBANDS = -754.79042976 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.31430140 eV energy without entropy = -173.34030093 energy(sigma->0) = -173.32296791 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 1938 total energy-change (2. order) : 0.5815150E+01 (-0.1984289E+02) number of electron 136.0000023 magnetization 0.4879639 augmentation part -6.9150712 magnetization 0.4599986 Broyden mixing: rms(total) = 0.40690E+00 rms(broyden)= 0.40670E+00 rms(prec ) = 0.41501E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0008 6.7266 2.6039 1.4735 1.4735 1.4064 0.8808 0.5389 0.5389 0.5634 0.5634 0.5682 0.5682 0.3999 0.2759 0.2759 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12495.98768211 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.77220599 PAW double counting = 12907.03159571 -12348.61016301 entropy T*S EENTRO = 0.02618229 eigenvalues EBANDS = -750.28128793 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -167.49915114 eV energy without entropy = -167.52533343 energy(sigma->0) = -167.50787857 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 1899 total energy-change (2. order) :-0.3613661E+01 (-0.2424858E+01) number of electron 136.0000023 magnetization 0.4351888 augmentation part -6.9122155 magnetization 0.4164038 Broyden mixing: rms(total) = 0.12553E+00 rms(broyden)= 0.12544E+00 rms(prec ) = 0.13331E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9911 6.8943 2.5177 1.5376 1.5376 1.0797 0.7777 0.7777 0.8079 0.6530 0.6530 0.5023 0.5023 0.5192 0.2760 0.2760 0.3540 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12494.90169169 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.76445315 PAW double counting = 12917.76514676 -12359.37181929 entropy T*S EENTRO = 0.02558033 eigenvalues EBANDS = -754.95998534 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.11281250 eV energy without entropy = -171.13839283 energy(sigma->0) = -171.12133927 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 1959 total energy-change (2. order) :-0.3003508E+01 (-0.8137852E+00) number of electron 136.0000024 magnetization 0.1634542 augmentation part -6.9092221 magnetization 0.1535210 Broyden mixing: rms(total) = 0.27997E+00 rms(broyden)= 0.27981E+00 rms(prec ) = 0.29857E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9942 7.2250 2.3012 1.7116 1.7116 1.3991 0.8477 0.8477 0.7378 0.6181 0.5399 0.4829 0.4829 0.4796 0.4796 0.2760 0.2760 0.3041 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12494.92228573 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.78028446 PAW double counting = 12909.17040878 -12350.82795814 entropy T*S EENTRO = 0.02548528 eigenvalues EBANDS = -757.87609580 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.11632017 eV energy without entropy = -174.14180545 energy(sigma->0) = -174.12481527 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 2043 total energy-change (2. order) :-0.1659521E+00 (-0.1379417E-01) number of electron 136.0000023 magnetization -0.0922222 augmentation part -6.9085842 magnetization -0.0960608 Broyden mixing: rms(total) = 0.27040E+00 rms(broyden)= 0.27038E+00 rms(prec ) = 0.28982E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0115 7.5116 2.0631 2.0631 1.9868 0.9617 0.9617 0.9802 0.9802 0.7092 0.5133 0.5133 0.4565 0.4565 0.5228 0.5228 0.2760 0.2760 0.3403 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12494.42445379 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.86025583 PAW double counting = 12787.17813942 -12228.88275528 entropy T*S EENTRO = 0.02598043 eigenvalues EBANDS = -758.41333713 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.28227228 eV energy without entropy = -174.30825271 energy(sigma->0) = -174.29093242 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 2139 total energy-change (2. order) :-0.4430267E+00 (-0.4920161E-02) number of electron 136.0000023 magnetization -0.2291027 augmentation part -6.9074434 magnetization -0.2294885 Broyden mixing: rms(total) = 0.26679E+00 rms(broyden)= 0.26679E+00 rms(prec ) = 0.28790E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0647 7.7968 2.4716 2.4716 2.2245 1.1491 1.1491 0.8136 0.8136 0.7385 0.5353 0.5353 0.5320 0.5320 0.5776 0.5776 0.5290 0.2760 0.2760 0.3316 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12493.57316951 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.95768693 PAW double counting = 12620.81079675 -12062.67296765 entropy T*S EENTRO = 0.02650073 eigenvalues EBANDS = -759.45318223 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.72529894 eV energy without entropy = -174.75179968 energy(sigma->0) = -174.73413252 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 2121 total energy-change (2. order) : 0.6036995E+00 (-0.8644140E-02) number of electron 136.0000023 magnetization -0.2691679 augmentation part -6.9062630 magnetization -0.2688079 Broyden mixing: rms(total) = 0.23147E+00 rms(broyden)= 0.23147E+00 rms(prec ) = 0.25212E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0739 7.8720 2.6837 2.6837 2.2245 1.2405 1.2405 0.8274 0.8274 0.6768 0.6768 0.5578 0.5578 0.5080 0.5080 0.6206 0.5211 0.5065 0.2760 0.2760 0.3310 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12493.20917270 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00014988 PAW double counting = 12497.79487362 -11939.51237546 entropy T*S EENTRO = 0.02675616 eigenvalues EBANDS = -759.31594107 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.12159944 eV energy without entropy = -174.14835560 energy(sigma->0) = -174.13051816 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 1977 total energy-change (2. order) : 0.1815957E+01 (-0.2785242E+00) number of electron 136.0000023 magnetization -0.2652619 augmentation part -6.9084722 magnetization -0.2700283 Broyden mixing: rms(total) = 0.14009E+00 rms(broyden)= 0.13996E+00 rms(prec ) = 0.15343E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0567 7.8731 2.7281 2.7281 2.3200 1.2866 1.2866 0.8231 0.8231 0.6586 0.5859 0.5859 0.5259 0.5259 0.5663 0.5663 0.5166 0.4854 0.4854 0.2760 0.2760 0.3372 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.99600576 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.01101813 PAW double counting = 12441.36136277 -11883.08594657 entropy T*S EENTRO = 0.02698050 eigenvalues EBANDS = -757.69542552 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.30564281 eV energy without entropy = -172.33262331 energy(sigma->0) = -172.31463631 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 1938 total energy-change (2. order) :-0.1528247E+01 (-0.2115627E+00) number of electron 136.0000023 magnetization -0.2826657 augmentation part -6.9050104 magnetization -0.2847712 Broyden mixing: rms(total) = 0.20661E+00 rms(broyden)= 0.20657E+00 rms(prec ) = 0.22597E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0600 7.9316 2.8671 2.8671 2.3453 1.3364 1.3364 0.8279 0.8279 0.7464 0.7464 0.6972 0.5360 0.5360 0.5342 0.5342 0.5486 0.5486 0.5007 0.2760 0.2760 0.3242 0.2591 0.0330 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.90567128 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00767921 PAW double counting = 12435.18025847 -11876.94257369 entropy T*S EENTRO = 0.02698679 eigenvalues EBANDS = -759.27962029 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.83388933 eV energy without entropy = -173.86087611 energy(sigma->0) = -173.84288492 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 1929 total energy-change (2. order) : 0.9732388E+00 (-0.5088260E-01) number of electron 136.0000023 magnetization -0.2805401 augmentation part -6.9066632 magnetization -0.2814314 Broyden mixing: rms(total) = 0.16397E+00 rms(broyden)= 0.16396E+00 rms(prec ) = 0.17870E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0307 7.9494 2.8801 2.8801 2.3877 1.3438 1.3438 0.8338 0.8338 0.7666 0.7666 0.6827 0.5302 0.5302 0.5367 0.5367 0.5563 0.5563 0.5026 0.2760 0.2760 0.3261 0.1228 0.0330 0.2525 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.80634043 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00362311 PAW double counting = 12424.28303430 -11866.22202544 entropy T*S EENTRO = 0.02697156 eigenvalues EBANDS = -758.23307730 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.86065054 eV energy without entropy = -172.88762210 energy(sigma->0) = -172.86964106 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 1923 total energy-change (2. order) :-0.3164520E-01 (-0.1032531E-02) number of electron 136.0000023 magnetization -0.2742732 augmentation part -6.9070147 magnetization -0.2753212 Broyden mixing: rms(total) = 0.16202E+00 rms(broyden)= 0.16202E+00 rms(prec ) = 0.17702E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0097 7.9224 2.8967 2.8967 2.4436 1.3670 1.3670 0.8149 0.8149 0.7538 0.7538 0.7124 0.5343 0.5343 0.5348 0.5348 0.5511 0.5511 0.5054 0.2760 0.2760 0.3259 0.2444 0.2444 0.0330 0.2582 0.0546 0.0546 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.81682117 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00348328 PAW double counting = 12426.92641231 -11868.85399727 entropy T*S EENTRO = 0.02697144 eigenvalues EBANDS = -758.26578766 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.89229574 eV energy without entropy = -172.91926718 energy(sigma->0) = -172.90128622 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 1680 total energy-change (2. order) :-0.2411050E+02 (-0.4805191E+03) number of electron 136.0000027 magnetization -0.2439667 augmentation part -7.0158891 magnetization -0.1761022 Broyden mixing: rms(total) = 0.10135E+01 rms(broyden)= 0.10131E+01 rms(prec ) = 0.12876E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9698 7.9237 2.8192 2.8192 2.5119 1.3324 1.3324 0.8382 0.8382 0.7543 0.7543 0.7202 0.5348 0.5348 0.5373 0.5373 0.5586 0.5586 0.5031 0.2760 0.2760 0.3288 0.2201 0.2201 0.2314 0.0330 0.0546 0.0546 0.0149 0.0076 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.88046323 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00270967 PAW double counting = 12429.96602329 -11910.45548867 entropy T*S EENTRO = -0.04052284 eigenvalues EBANDS = -743.68404239 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -197.00279363 eV energy without entropy = -196.96227079 energy(sigma->0) = -196.98928602 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 2217 total energy-change (2. order) : 0.2487091E+02 (-0.4049328E+01) number of electron 136.0000025 magnetization -0.2336845 augmentation part -6.9010907 magnetization -0.2354836 Broyden mixing: rms(total) = 0.44983E+00 rms(broyden)= 0.44898E+00 rms(prec ) = 0.48788E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9429 7.9141 2.7921 2.7921 2.5462 1.3396 1.3396 0.8800 0.8800 0.7216 0.7216 0.5677 0.5677 0.5246 0.5246 0.6295 0.5921 0.5921 0.5155 0.3319 0.2760 0.2760 0.2816 0.2816 0.2033 0.0330 0.0546 0.0546 0.0192 0.0192 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.69625896 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.99296108 PAW double counting = 12451.33186332 -11893.07523986 entropy T*S EENTRO = 0.02681488 eigenvalues EBANDS = -757.82050740 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.13187921 eV energy without entropy = -172.15869409 energy(sigma->0) = -172.14081751 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 1992 total energy-change (2. order) :-0.2829952E+01 (-0.1402438E+01) number of electron 136.0000024 magnetization -0.2150816 augmentation part -6.9038026 magnetization -0.2160013 Broyden mixing: rms(total) = 0.27063E+00 rms(broyden)= 0.27052E+00 rms(prec ) = 0.28952E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9426 7.9063 2.6685 2.6685 2.6481 1.3779 1.3779 1.0729 1.0729 0.7634 0.7634 0.7093 0.5505 0.5505 0.5383 0.5383 0.5917 0.5917 0.5230 0.4632 0.4632 0.2760 0.2760 0.3317 0.2164 0.0330 0.0546 0.0546 0.0272 0.0272 0.0149 0.0713 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.41874588 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.01164505 PAW double counting = 12441.31675733 -11883.02744698 entropy T*S EENTRO = 0.02700321 eigenvalues EBANDS = -760.94216350 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.96183099 eV energy without entropy = -174.98883420 energy(sigma->0) = -174.97083206 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 1956 total energy-change (2. order) : 0.1051399E+01 (-0.1139568E+00) number of electron 136.0000023 magnetization -0.1704441 augmentation part -6.9059086 magnetization -0.1738808 Broyden mixing: rms(total) = 0.20146E+00 rms(broyden)= 0.20140E+00 rms(prec ) = 0.22049E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9625 7.7858 2.8730 2.8875 2.8875 1.5204 1.5204 1.2333 1.2333 0.8267 0.7295 0.7295 0.5384 0.5384 0.5474 0.5474 0.5605 0.5605 0.4947 0.4947 0.4518 0.4518 0.2760 0.2760 0.3279 0.2149 0.0330 0.0546 0.0546 0.0275 0.0275 0.0149 0.0809 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12491.85607395 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.02209666 PAW double counting = 12467.56293067 -11909.27913080 entropy T*S EENTRO = 0.02697445 eigenvalues EBANDS = -760.43744609 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.91043248 eV energy without entropy = -173.93740693 energy(sigma->0) = -173.91942397 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 2322 total energy-change (2. order) :-0.3077698E+02 (-0.1873433E+02) number of electron 136.0000050 magnetization -0.1695167 augmentation part -6.8235719 magnetization -0.7886665 Broyden mixing: rms(total) = 0.28846E+01 rms(broyden)= 0.28707E+01 rms(prec ) = 0.29500E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9332 7.7870 2.8712 2.8849 2.8849 1.5207 1.5207 1.2343 1.2343 0.8270 0.7299 0.7299 0.5381 0.5381 0.5474 0.5474 0.5612 0.5612 0.4938 0.4938 0.4511 0.4511 0.2760 0.2760 0.3279 0.2149 0.0330 0.0546 0.0546 0.0809 0.0275 0.0275 0.0149 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12491.91775318 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.00320468 PAW double counting = 12503.04847766 -11946.42056728 entropy T*S EENTRO = -0.01066886 eigenvalues EBANDS = -789.47810936 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -204.68741582 eV energy without entropy = -204.67674696 energy(sigma->0) = -204.68385953 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 2544 total energy-change (2. order) :-0.4530689E+04 (-0.2334175E+04) number of electron 135.9999944 magnetization -0.1682888 augmentation part -6.8444767 magnetization 1.7195547 Broyden mixing: rms(total) = 0.98868E+01 rms(broyden)= 0.97606E+01 rms(prec ) = 0.99667E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9056 7.7882 2.8694 2.8752 2.8752 1.5283 1.5283 1.2364 1.2364 0.8243 0.7317 0.7317 0.5376 0.5376 0.5472 0.5472 0.5630 0.5630 0.4940 0.4940 0.4472 0.4472 0.2760 0.2760 0.3279 0.2147 0.0330 0.0546 0.0546 0.0810 0.0275 0.0275 0.0149 0.0005 0.0003 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12491.93729345 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.53443639 PAW double counting = 12502.99123443 -11947.13121264 entropy T*S EENTRO = -0.02014680 eigenvalues EBANDS = -5319.83904905 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -4735.37649401 eV energy without entropy = -4735.35634721 energy(sigma->0) = -4735.36977841 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 1932 total energy-change (2. order) : 0.4578537E+04 (-0.9192350E+03) number of electron 136.0000020 magnetization -0.3802470 augmentation part -6.8660909 magnetization 1.1345543 Broyden mixing: rms(total) = 0.13604E+01 rms(broyden)= 0.47923E+00 rms(prec ) = 0.49764E+00 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8794 7.7885 2.8563 2.8741 2.8741 1.5256 1.5256 1.2376 1.2376 0.8268 0.7318 0.7318 0.5377 0.5377 0.5472 0.5472 0.5620 0.5620 0.4933 0.4933 0.4488 0.4488 0.2760 0.2760 0.3279 0.2145 0.0330 0.0546 0.0546 0.0810 0.0275 0.0275 0.0149 0.0010 0.0010 0.0009 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12492.27265989 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -252.51379804 PAW double counting = 12504.36127215 -11947.27936421 entropy T*S EENTRO = -0.01625948 eigenvalues EBANDS = -747.21277647 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -156.83917605 eV energy without entropy = -156.82291657 energy(sigma->0) = -156.83375622 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 2100 total energy-change (2. order) :-0.1721297E+02 (-0.5675210E+01) number of electron 136.0000029 magnetization -0.2603348 augmentation part -6.9055148 magnetization -0.2415249 Broyden mixing: rms(total) = 0.34724E+00 rms(broyden)= 0.23345E+00 rms(prec ) = 0.24910E+00 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8564 7.7883 2.8477 2.8757 2.8757 1.5280 1.5280 1.2374 1.2374 0.8272 0.7315 0.7315 0.5376 0.5376 0.5472 0.5472 0.5630 0.5630 0.4936 0.4936 0.4483 0.4483 0.3279 0.2760 0.2760 0.2148 0.0330 0.0809 0.0546 0.0546 0.0275 0.0275 0.0267 0.0267 0.0008 0.0008 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.75011142 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.36976094 PAW double counting = 12788.23379096 -12230.11888526 entropy T*S EENTRO = 0.02544182 eigenvalues EBANDS = -763.16703346 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.05214841 eV energy without entropy = -174.07759023 energy(sigma->0) = -174.06062902 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 2109 total energy-change (2. order) :-0.3551032E+00 (-0.6575198E+00) number of electron 136.0000024 magnetization -0.2635469 augmentation part -6.9094101 magnetization -0.2614451 Broyden mixing: rms(total) = 0.23552E+00 rms(broyden)= 0.22972E+00 rms(prec ) = 0.24655E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8384 7.7890 2.8456 2.8749 2.8749 1.5292 1.5292 1.2371 1.2371 0.8272 0.7317 0.7317 0.5373 0.5373 0.5474 0.5474 0.5624 0.5624 0.4934 0.4934 0.4479 0.4479 0.3279 0.2760 0.2760 0.2143 0.0899 0.0899 0.0820 0.0330 0.0655 0.0546 0.0546 0.0274 0.0274 0.0008 0.0008 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12490.49245275 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93382944 PAW double counting = 12842.78122770 -12284.82798693 entropy T*S EENTRO = 0.02624684 eigenvalues EBANDS = -762.05486696 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.40725166 eV energy without entropy = -174.43349850 energy(sigma->0) = -174.41600061 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 2073 total energy-change (2. order) : 0.2531294E+00 (-0.7681237E-01) number of electron 136.0000023 magnetization -0.2556143 augmentation part -6.9120792 magnetization -0.2515885 Broyden mixing: rms(total) = 0.19699E+00 rms(broyden)= 0.19671E+00 rms(prec ) = 0.20888E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8687 7.8214 2.7211 2.7211 2.7134 1.7707 1.7707 1.1945 1.1945 0.8237 0.8237 0.8090 0.6994 0.6994 0.5716 0.5716 0.5280 0.5280 0.5415 0.5415 0.5697 0.5697 0.4943 0.4943 0.3265 0.2760 0.2760 0.2533 0.2066 0.2066 0.0330 0.0798 0.0546 0.0546 0.0274 0.0274 0.0149 0.0008 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.66112855 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.02160887 PAW double counting = 12765.76929190 -12208.45528793 entropy T*S EENTRO = 0.02639709 eigenvalues EBANDS = -761.90619580 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.15412227 eV energy without entropy = -174.18051936 energy(sigma->0) = -174.16292130 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 2145 total energy-change (2. order) :-0.5528401E+00 (-0.2342777E-01) number of electron 136.0000023 magnetization -0.2483835 augmentation part -6.9096695 magnetization -0.2453545 Broyden mixing: rms(total) = 0.21166E+00 rms(broyden)= 0.21165E+00 rms(prec ) = 0.22773E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8971 7.9204 2.7022 2.5372 2.5372 2.1993 2.1993 1.4281 1.4281 0.9657 0.9657 0.8651 0.6349 0.6349 0.6992 0.6992 0.5547 0.5547 0.5117 0.5117 0.5503 0.5503 0.4999 0.4999 0.4786 0.2760 0.2760 0.3267 0.2636 0.2116 0.2116 0.0330 0.0798 0.0546 0.0546 0.0274 0.0274 0.0149 0.0008 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.77174097 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.03914128 PAW double counting = 12732.44420130 -12174.87590062 entropy T*S EENTRO = 0.02650524 eigenvalues EBANDS = -762.58529589 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.70696233 eV energy without entropy = -174.73346756 energy(sigma->0) = -174.71579740 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 2100 total energy-change (2. order) :-0.8239294E+00 (-0.5059524E-01) number of electron 136.0000024 magnetization -0.2071128 augmentation part -6.9072554 magnetization -0.2036085 Broyden mixing: rms(total) = 0.27142E+00 rms(broyden)= 0.27139E+00 rms(prec ) = 0.29082E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8949 7.8912 2.6641 2.5171 2.5171 2.4428 2.4428 1.4415 1.4415 0.9890 0.9890 0.8524 0.7351 0.7351 0.6252 0.6252 0.5517 0.5517 0.5112 0.5112 0.5697 0.5697 0.5166 0.4590 0.4590 0.3414 0.2760 0.2760 0.3051 0.2727 0.2117 0.2117 0.0330 0.0798 0.0546 0.0546 0.0274 0.0274 0.0008 0.0008 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.85169137 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.04205284 PAW double counting = 12682.69387750 -12124.91462787 entropy T*S EENTRO = 0.02655142 eigenvalues EBANDS = -763.53735843 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -175.53089168 eV energy without entropy = -175.55744311 energy(sigma->0) = -175.53974216 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 2028 total energy-change (2. order) : 0.1170360E+00 (-0.7286618E-02) number of electron 136.0000024 magnetization -0.1568867 augmentation part -6.9073795 magnetization -0.1539092 Broyden mixing: rms(total) = 0.27531E+00 rms(broyden)= 0.27530E+00 rms(prec ) = 0.29355E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9087 7.7552 2.9132 2.6749 2.6749 2.2282 2.2282 1.5280 1.5280 1.1535 1.1535 0.8611 0.7513 0.7513 0.6359 0.6359 0.5949 0.5949 0.5617 0.5617 0.5164 0.5164 0.5611 0.5611 0.4969 0.4969 0.4639 0.2760 0.2760 0.3256 0.2646 0.2116 0.2116 0.0330 0.0798 0.0546 0.0546 0.0274 0.0274 0.0008 0.0008 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.77770727 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.03872490 PAW double counting = 12676.48898099 -12118.71452869 entropy T*S EENTRO = 0.02654532 eigenvalues EBANDS = -763.49283103 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -175.41385569 eV energy without entropy = -175.44040101 energy(sigma->0) = -175.42270413 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 1707 total energy-change (2. order) : 0.6990274E+01 (-0.5667803E+01) number of electron 136.0000021 magnetization -0.1423990 augmentation part -6.9248637 magnetization -0.1401781 Broyden mixing: rms(total) = 0.41970E+00 rms(broyden)= 0.41928E+00 rms(prec ) = 0.44675E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8845 7.6728 3.1542 2.8940 2.8940 1.8289 1.8289 1.5386 1.5386 1.1900 1.1900 0.8759 0.7286 0.7286 0.6095 0.6095 0.5841 0.5841 0.5482 0.5482 0.5228 0.5228 0.5585 0.5585 0.5034 0.4892 0.4892 0.3256 0.2760 0.2760 0.2646 0.2116 0.2116 0.0330 0.0981 0.0798 0.0546 0.0546 0.0274 0.0274 0.0149 0.0008 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.60554426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.03993896 PAW double counting = 12674.00895514 -12118.92659417 entropy T*S EENTRO = 0.02658800 eigenvalues EBANDS = -753.98145777 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -168.42358214 eV energy without entropy = -168.45017014 energy(sigma->0) = -168.43244481 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 1806 total energy-change (2. order) :-0.3032583E+01 (-0.2493852E+01) number of electron 136.0000023 magnetization -0.1442680 augmentation part -6.9160709 magnetization -0.1501956 Broyden mixing: rms(total) = 0.20959E+00 rms(broyden)= 0.20931E+00 rms(prec ) = 0.21342E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8646 7.6688 3.2099 2.8404 2.8404 1.8132 1.8132 1.5762 1.5762 1.1971 1.1971 0.8713 0.7307 0.7307 0.6128 0.6128 0.5781 0.5781 0.5453 0.5453 0.5213 0.5213 0.5639 0.5639 0.4977 0.4918 0.4918 0.3255 0.2760 0.2760 0.2648 0.2116 0.2116 0.0330 0.0274 0.0274 0.0798 0.0790 0.0546 0.0546 0.0500 0.0008 0.0008 0.0149 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.46255777 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.04211339 PAW double counting = 12674.25383682 -12116.18731121 entropy T*S EENTRO = 0.02656630 eigenvalues EBANDS = -760.13899531 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.45616467 eV energy without entropy = -171.48273097 energy(sigma->0) = -171.46502011 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 1734 total energy-change (2. order) :-0.2122296E+01 (-0.7465237E+00) number of electron 136.0000024 magnetization -0.1080263 augmentation part -6.9109994 magnetization -0.1089877 Broyden mixing: rms(total) = 0.17649E+00 rms(broyden)= 0.17625E+00 rms(prec ) = 0.18913E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8783 7.5809 3.2736 3.2736 3.2037 2.2431 1.6161 1.4318 1.4318 1.2512 1.2512 0.8408 0.7240 0.7240 0.6076 0.6076 0.6888 0.6888 0.5510 0.5510 0.5213 0.5213 0.5439 0.5439 0.5263 0.5263 0.4765 0.3268 0.2760 0.2760 0.2618 0.2376 0.2376 0.2115 0.2115 0.1137 0.0330 0.0798 0.0546 0.0546 0.0274 0.0274 0.0149 0.0008 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.41054168 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.04151566 PAW double counting = 12668.55239662 -12110.71235799 entropy T*S EENTRO = 0.02656773 eigenvalues EBANDS = -762.08741913 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.57846021 eV energy without entropy = -173.60502794 energy(sigma->0) = -173.58731612 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 1953 total energy-change (2. order) : 0.2558506E+01 (-0.1063886E+01) number of electron 136.0000022 magnetization -0.0683694 augmentation part -6.9166846 magnetization -0.0752704 Broyden mixing: rms(total) = 0.21954E+00 rms(broyden)= 0.21937E+00 rms(prec ) = 0.22427E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8583 7.5078 2.9718 2.9298 2.9298 2.0038 1.7458 1.7458 1.1014 1.1014 0.6561 0.6561 0.7303 0.7303 0.5956 0.5956 0.6237 0.6237 0.5064 0.5064 0.5483 0.4962 0.4962 0.2019 0.2019 0.4236 0.3775 0.1986 0.1986 0.2340 0.1874 0.1874 0.1391 0.0528 0.0528 0.0379 0.0151 0.0151 0.0012 0.0012 0.0017 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.48626089 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.04105393 PAW double counting = 12665.03399755 -12107.14334768 entropy T*S EENTRO = 0.02654173 eigenvalues EBANDS = -759.50424131 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.01995464 eV energy without entropy = -171.04649638 energy(sigma->0) = -171.02880189 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 1911 total energy-change (2. order) :-0.3132667E+01 (-0.8170598E+00) number of electron 136.0000024 magnetization -0.0580712 augmentation part -6.9115508 magnetization -0.0599851 Broyden mixing: rms(total) = 0.18142E+00 rms(broyden)= 0.18122E+00 rms(prec ) = 0.19435E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8618 7.4835 3.1024 3.0059 3.0059 1.9743 1.8248 1.8248 1.0999 1.0999 0.6632 0.6632 0.7009 0.7009 0.5870 0.5870 0.6582 0.6582 0.6113 0.5349 0.5349 0.4916 0.4916 0.2557 0.2557 0.4390 0.2954 0.2954 0.3267 0.3267 0.1576 0.1576 0.1605 0.1605 0.0621 0.0621 0.0361 0.0149 0.0149 0.0080 0.0006 0.0006 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8998.00662219 -Hartree energ DENC = -12489.62920887 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.02459861 PAW double counting = 12701.77839605 -12144.34245236 entropy T*S EENTRO = 0.02648586 eigenvalues EBANDS = -762.05565380 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.15262184 eV energy without entropy = -174.17910770 energy(sigma->0) = -174.16145046 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 8 1 9 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------