vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.09.02 06:45:34 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.929 0.527 0.107- 19 2.31 14 2.39 3 2.43 2 0.850 0.233 0.428- 4 2.34 13 2.38 20 2.42 3 0.136 0.448 0.196- 8 2.31 17 2.35 7 2.38 1 2.43 4 0.092 0.283 0.334- 18 2.29 2 2.34 7 2.35 8 2.63 5 0.818 0.332 0.027- 14 2.34 23 2.34 7 2.35 10 2.37 6 0.957 0.421 0.518- 8 2.33 13 2.37 24 2.39 11 2.42 9 2.62 7 0.082 0.333 0.137- 21 2.34 4 2.35 5 2.35 3 2.38 8 0.199 0.406 0.390- 3 2.31 6 2.33 22 2.38 4 2.63 18 2.64 9 0.832 0.545 0.524- 11 2.36 27 2.36 6 2.62 10 0.772 0.243 0.882- 5 2.37 28 2.38 12 2.40 31 2.46 11 0.063 0.510 0.653- 16 2.35 9 2.36 25 2.42 6 2.42 12 0.014 0.205 0.760- 26 2.35 10 2.40 15 2.52 13 0.902 0.310 0.593- 31 2.32 6 2.37 2 2.38 15 2.40 14 0.864 0.443 0.955- 16 2.33 5 2.34 32 2.35 1 2.39 15 0.139 0.323 0.737- 29 2.33 16 2.38 13 2.40 26 2.50 12 2.52 16 0.098 0.434 0.819- 14 2.33 30 2.35 11 2.35 15 2.38 17 0.422 0.461 0.122- 19 2.30 30 2.33 3 2.35 23 2.50 18 0.363 0.290 0.422- 4 2.29 29 2.32 20 2.35 8 2.64 19 0.691 0.483 0.208- 17 2.30 1 2.31 24 2.34 20 0.642 0.313 0.341- 23 2.34 18 2.35 24 2.42 2 2.42 21 0.326 0.317 0.011- 7 2.34 26 2.35 30 2.38 23 2.38 22 0.434 0.439 0.521- 24 2.37 8 2.38 29 2.44 27 2.48 25 2.48 23 0.562 0.347 0.142- 20 2.34 5 2.34 21 2.38 17 2.50 24 0.692 0.429 0.401- 19 2.34 22 2.37 6 2.39 20 2.42 25 0.317 0.555 0.549- 35 1.79 11 2.42 22 2.48 27 2.57 26 0.274 0.231 0.865- 12 2.35 21 2.35 28 2.40 15 2.50 27 0.599 0.517 0.657- 35 1.67 9 2.36 32 2.36 22 2.48 25 2.57 28 0.538 0.208 0.753- 10 2.38 26 2.40 31 2.52 29 0.396 0.329 0.623- 18 2.32 15 2.33 31 2.37 22 2.44 30 0.368 0.426 0.920- 32 2.32 17 2.33 16 2.35 21 2.38 31 0.668 0.325 0.727- 13 2.32 29 2.37 32 2.38 10 2.46 28 2.52 32 0.630 0.435 0.815- 30 2.32 14 2.35 27 2.36 31 2.38 33 0.655 0.707 0.621- 34 0.71 34 0.674 0.719 0.562- 33 0.71 35 0.457 0.581 0.672- 27 1.67 25 1.79 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.929406980 0.527387170 0.106886350 0.849783440 0.233030560 0.427849660 0.136154000 0.447842870 0.195880950 0.092245840 0.283251420 0.333987250 0.818017580 0.331903940 0.026918890 0.956558470 0.421053510 0.517990960 0.081965280 0.333162280 0.137219590 0.198842230 0.406332160 0.390213410 0.831909820 0.545119510 0.524080000 0.771921900 0.243266240 0.881815510 0.062908570 0.509854010 0.652884890 0.014070110 0.205285450 0.759857240 0.901994970 0.309923740 0.593388480 0.863587970 0.442721760 0.955365840 0.138936910 0.323266660 0.736871930 0.098094180 0.434234840 0.819276300 0.422185930 0.460867770 0.121857810 0.363495500 0.289955220 0.421521580 0.690735080 0.482574400 0.208252370 0.641778380 0.312616680 0.341209770 0.326137240 0.316679540 0.010860400 0.433803500 0.438532180 0.520898380 0.562383700 0.346786890 0.142224680 0.692166720 0.429325360 0.401433560 0.316636390 0.555038110 0.549098790 0.273594930 0.231104860 0.864634900 0.598745770 0.516847910 0.656702450 0.537909390 0.208401410 0.752556330 0.396241120 0.328953300 0.622714620 0.367853070 0.425656000 0.920424590 0.668398150 0.325368570 0.726611590 0.630150820 0.435392240 0.815204460 0.654977680 0.706602660 0.621078020 0.674453910 0.719373190 0.561774960 0.457414800 0.580959710 0.672345670 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.92940698 0.52738717 0.10688635 0.84978344 0.23303056 0.42784966 0.13615400 0.44784287 0.19588095 0.09224584 0.28325142 0.33398725 0.81801758 0.33190394 0.02691889 0.95655847 0.42105351 0.51799096 0.08196528 0.33316228 0.13721959 0.19884223 0.40633216 0.39021341 0.83190982 0.54511951 0.52408000 0.77192190 0.24326624 0.88181551 0.06290857 0.50985401 0.65288489 0.01407011 0.20528545 0.75985724 0.90199497 0.30992374 0.59338848 0.86358797 0.44272176 0.95536584 0.13893691 0.32326666 0.73687193 0.09809418 0.43423484 0.81927630 0.42218593 0.46086777 0.12185781 0.36349550 0.28995522 0.42152158 0.69073508 0.48257440 0.20825237 0.64177838 0.31261668 0.34120977 0.32613724 0.31667954 0.01086040 0.43380350 0.43853218 0.52089838 0.56238370 0.34678689 0.14222468 0.69216672 0.42932536 0.40143356 0.31663639 0.55503811 0.54909879 0.27359493 0.23110486 0.86463490 0.59874577 0.51684791 0.65670245 0.53790939 0.20840141 0.75255633 0.39624112 0.32895330 0.62271462 0.36785307 0.42565600 0.92042459 0.66839815 0.32536857 0.72661159 0.63015082 0.43539224 0.81520446 0.65497768 0.70660266 0.62107802 0.67445391 0.71937319 0.56177496 0.45741480 0.58095971 0.67234567 position of ions in cartesian coordinates (Angst): 7.12213863 10.37006666 1.15835517 6.51197548 4.58210320 4.63671801 1.04336172 8.80597914 2.12281278 0.70688910 5.56960100 3.61950666 6.26855052 6.52626036 0.29172701 7.33020321 8.27921727 5.61360271 0.62810814 6.55100323 1.48708437 1.52374789 7.98974990 4.22884418 6.37500814 10.71873944 5.67959122 5.91531471 4.78336840 9.55646395 0.48207466 10.02531038 7.07548330 0.10782066 4.03654833 8.23477047 6.91207765 6.09406149 6.43070524 6.61776097 8.70528224 10.35354800 1.06468744 6.35642466 7.98567269 0.75170551 8.53840308 8.87871027 3.23525300 9.06208905 1.32060477 2.78550237 5.70141849 4.56813896 5.29317199 9.48890868 2.25688508 4.91801190 6.14701304 3.69777899 2.49922228 6.22690146 0.11769698 3.32427960 8.62290211 5.64511118 4.30960253 6.81890530 1.54132584 5.30414279 8.44186749 4.35043986 2.42641632 10.91376986 5.95072635 2.09658531 4.54423797 9.37027322 4.58824871 10.16283214 7.11685519 4.12205345 4.09781776 8.15564861 3.03643533 6.46824163 6.74852024 2.81889486 8.36971649 9.97488059 5.12200186 6.39775473 7.87447872 4.82890875 8.56116115 8.83458269 5.01915946 13.89399876 6.73078398 5.16840776 14.14510697 6.08810130 3.50521535 11.42346887 7.28638484 -------------------------------------------------------------------------------------------------------- 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 254499. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3167. 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 1358 Maximum index for augmentation-charges 1470 (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.5987356E+03 (-0.3891245E+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 = 8963.22719592 -Hartree energ DENC = -12244.49188456 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.17710369 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00451378 eigenvalues EBANDS = -197.27618299 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 598.73557911 eV energy without entropy = 598.74009289 energy(sigma->0) = 598.73708370 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 2271 total energy-change (2. order) :-0.6877514E+03 (-0.6533468E+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 = 8963.22719592 -Hartree energ DENC = -12244.49188456 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.17710369 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00014127 eigenvalues EBANDS = -885.03194126 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -89.01580665 eV energy without entropy = -89.01566537 energy(sigma->0) = -89.01575956 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2100 total energy-change (2. order) :-0.7313717E+02 (-0.7136096E+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 = 8963.22719592 -Hartree energ DENC = -12244.49188456 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.17710369 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02495916 eigenvalues EBANDS = -958.14429131 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -162.15297459 eV energy without entropy = -162.12801543 energy(sigma->0) = -162.14465487 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 2226 total energy-change (2. order) :-0.3028839E+01 (-0.3011272E+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 = 8963.22719592 -Hartree energ DENC = -12244.49188456 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.17710369 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02713018 eigenvalues EBANDS = -961.17095946 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.18181376 eV energy without entropy = -165.15468358 energy(sigma->0) = -165.17277037 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 2397 total energy-change (2. order) :-0.1481265E+00 (-0.1480491E+00) number of electron 135.9999985 magnetization 30.2871384 augmentation part -6.9899398 magnetization 27.1324386 Broyden mixing: rms(total) = 0.25347E+01 rms(broyden)= 0.25345E+01 rms(prec ) = 0.26729E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12244.49188456 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.17710369 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02718023 eigenvalues EBANDS = -961.31903590 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.32994025 eV energy without entropy = -165.30276002 energy(sigma->0) = -165.32088017 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 1809 total energy-change (2. order) : 0.2755493E+02 (-0.5772139E+01) number of electron 135.9999986 magnetization 24.6835586 augmentation part -7.0365733 magnetization 21.3427954 Broyden mixing: rms(total) = 0.15933E+01 rms(broyden)= 0.15932E+01 rms(prec ) = 0.16807E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0254 1.0254 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12353.22406331 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.88620005 PAW double counting = 6201.98535049 -5642.48767300 entropy T*S EENTRO = -0.00642519 eigenvalues EBANDS = -841.47098256 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -137.77501000 eV energy without entropy = -137.76858480 energy(sigma->0) = -137.77286827 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 2235 total energy-change (2. order) :-0.2838734E+03 (-0.1991659E+03) number of electron 135.9999980 magnetization 24.3783163 augmentation part -7.1117369 magnetization 23.3301853 Broyden mixing: rms(total) = 0.52206E+01 rms(broyden)= 0.51149E+01 rms(prec ) = 0.54239E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5458 1.0410 0.0506 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12450.23196935 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -240.35701699 PAW double counting = 10050.00437181 -9491.82849135 entropy T*S EENTRO = -0.03972612 eigenvalues EBANDS = -1032.51054129 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -421.64838967 eV energy without entropy = -421.60866355 energy(sigma->0) = -421.63514763 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 1998 total energy-change (2. order) : 0.2921021E+03 (-0.2648378E+03) number of electron 135.9999988 magnetization 21.0698157 augmentation part -7.0852278 magnetization 18.5893141 Broyden mixing: rms(total) = 0.14318E+01 rms(broyden)= 0.11516E+01 rms(prec ) = 0.11914E+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.7133 1.4175 0.6467 0.0758 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12447.17517543 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -235.45892785 PAW double counting = 10172.05222617 -9613.49204320 entropy T*S EENTRO = -0.03598808 eigenvalues EBANDS = -748.75132773 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -129.54625249 eV energy without entropy = -129.51026441 energy(sigma->0) = -129.53425646 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 2382 total energy-change (2. order) :-0.1172616E+02 (-0.1518062E+01) number of electron 135.9999988 magnetization 17.2375172 augmentation part -7.0282149 magnetization 15.1972464 Broyden mixing: rms(total) = 0.88196E+00 rms(broyden)= 0.86509E+00 rms(prec ) = 0.89881E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8478 1.8525 0.8579 0.6055 0.0752 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12466.18280441 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -239.17795591 PAW double counting = 11665.72175240 -11107.07920066 entropy T*S EENTRO = -0.02931895 eigenvalues EBANDS = -737.83986524 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -141.27240915 eV energy without entropy = -141.24309019 energy(sigma->0) = -141.26263616 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 2847 total energy-change (2. order) :-0.8328421E+03 (-0.7303142E+03) number of electron 136.0000014 magnetization 17.2379919 augmentation part -6.8972186 magnetization 12.9355156 Broyden mixing: rms(total) = 0.83395E+01 rms(broyden)= 0.81624E+01 rms(prec ) = 0.89507E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6793 1.8534 0.8582 0.6054 0.0752 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12486.57232274 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -241.11030895 PAW double counting = 12754.68536873 -12190.44851975 entropy T*S EENTRO = 0.02125196 eigenvalues EBANDS = -1554.00499050 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -974.11453764 eV energy without entropy = -974.13578959 energy(sigma->0) = -974.12162162 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 2163 total energy-change (2. order) : 0.8338200E+03 (-0.6575499E+03) number of electron 135.9999987 magnetization 12.8291333 augmentation part -6.9190924 magnetization 11.6767973 Broyden mixing: rms(total) = 0.15518E+01 rms(broyden)= 0.69512E+00 rms(prec ) = 0.71975E+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.8450 2.6162 0.8821 0.8821 0.6113 0.0751 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12487.72573792 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -235.86783374 PAW double counting = 12727.76731934 -12170.67610217 entropy T*S EENTRO = -0.00333148 eigenvalues EBANDS = -717.10385482 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -140.29455717 eV energy without entropy = -140.29122569 energy(sigma->0) = -140.29344667 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 3027 total energy-change (2. order) :-0.3153374E+04 (-0.2842174E+04) number of electron 135.9999995 magnetization 12.7422825 augmentation part -6.7587138 magnetization 10.4828897 Broyden mixing: rms(total) = 0.96488E+01 rms(broyden)= 0.95440E+01 rms(prec ) = 0.10313E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7256 2.6245 0.8790 0.8790 0.6129 0.0751 0.0031 0.0058 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12513.03253310 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -241.90008961 PAW double counting = 14198.52623870 -13623.32299040 entropy T*S EENTRO = -0.01298454 eigenvalues EBANDS = -3857.24117803 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -3293.66855336 eV energy without entropy = -3293.65556882 energy(sigma->0) = -3293.66422518 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 2355 total energy-change (2. order) : 0.3150889E+04 (-0.7239612E+03) number of electron 135.9999988 magnetization 10.6659653 augmentation part -6.8718154 magnetization 9.1390448 Broyden mixing: rms(total) = 0.15333E+01 rms(broyden)= 0.59738E+00 rms(prec ) = 0.62147E+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.7397 2.9178 0.9852 0.7044 0.7044 0.5123 0.0751 0.0156 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12514.41361788 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -236.83761683 PAW double counting = 14228.42156097 -13671.64450567 entropy T*S EENTRO = -0.00020500 eigenvalues EBANDS = -691.61970155 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -142.77910234 eV energy without entropy = -142.77889734 energy(sigma->0) = -142.77903400 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 2382 total energy-change (2. order) :-0.1206123E+02 (-0.4159479E+01) number of electron 135.9999988 magnetization 8.2241268 augmentation part -6.9244103 magnetization 7.7716469 Broyden mixing: rms(total) = 0.55760E+00 rms(broyden)= 0.48429E+00 rms(prec ) = 0.50520E+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.7774 2.9111 1.2730 0.8230 0.8230 0.5366 0.5366 0.0751 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12514.78181180 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -244.91052490 PAW double counting = 14791.80924939 -14234.09235894 entropy T*S EENTRO = 0.00214510 eigenvalues EBANDS = -696.18201925 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -154.84033678 eV energy without entropy = -154.84248189 energy(sigma->0) = -154.84105182 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 2280 total energy-change (2. order) :-0.1581428E+02 (-0.3219553E+01) number of electron 135.9999988 magnetization 7.4889901 augmentation part -6.8890120 magnetization 7.4335459 Broyden mixing: rms(total) = 0.80370E+00 rms(broyden)= 0.79896E+00 rms(prec ) = 0.84633E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7624 2.9028 1.6322 0.8477 0.7592 0.5360 0.5360 0.0751 0.3172 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12511.13709028 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -249.05709337 PAW double counting = 14441.20689255 -13882.28182098 entropy T*S EENTRO = 0.00972316 eigenvalues EBANDS = -712.71021019 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -170.65461550 eV energy without entropy = -170.66433866 energy(sigma->0) = -170.65785655 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 2082 total energy-change (2. order) : 0.6873696E+01 (-0.2254581E+01) number of electron 135.9999988 magnetization 4.2694282 augmentation part -6.9135683 magnetization 4.2798188 Broyden mixing: rms(total) = 0.42542E+00 rms(broyden)= 0.42495E+00 rms(prec ) = 0.45121E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8370 2.7278 2.7278 0.8603 0.8603 0.6574 0.5237 0.5237 0.0751 0.2325 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12512.26298723 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -249.98391029 PAW double counting = 14304.97960703 -13746.91619886 entropy T*S EENTRO = 0.01476678 eigenvalues EBANDS = -702.92718088 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -163.78091985 eV energy without entropy = -163.79568662 energy(sigma->0) = -163.78584211 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 1992 total energy-change (2. order) :-0.7178140E+01 (-0.5790999E+00) number of electron 135.9999988 magnetization 1.3100821 augmentation part -6.8945899 magnetization 1.3865323 Broyden mixing: rms(total) = 0.40683E+00 rms(broyden)= 0.40680E+00 rms(prec ) = 0.43663E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0382 5.0328 2.5890 0.9790 0.9790 0.8478 0.6569 0.5220 0.5220 0.0751 0.2357 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12506.15303236 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.74003771 PAW double counting = 13691.20165138 -13132.68661945 entropy T*S EENTRO = 0.02705460 eigenvalues EBANDS = -712.92305975 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -170.95905967 eV energy without entropy = -170.98611427 energy(sigma->0) = -170.96807787 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 2010 total energy-change (2. order) :-0.1501385E+01 (-0.4483425E+00) number of electron 135.9999988 magnetization 1.1192223 augmentation part -6.8961558 magnetization 1.1263228 Broyden mixing: rms(total) = 0.27089E+00 rms(broyden)= 0.27080E+00 rms(prec ) = 0.29827E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0765 5.6678 2.5532 1.0919 1.0919 0.8231 0.8231 0.5734 0.5206 0.5206 0.0751 0.2354 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12502.68117099 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.91454111 PAW double counting = 13389.55789103 -12831.19493442 entropy T*S EENTRO = 0.02529571 eigenvalues EBANDS = -715.56796901 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.46044517 eV energy without entropy = -172.48574088 energy(sigma->0) = -172.46887707 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 2166 total energy-change (2. order) :-0.5244646E+01 (-0.2646386E+01) number of electron 135.9999988 magnetization 1.0513866 augmentation part -6.8848693 magnetization 1.0241183 Broyden mixing: rms(total) = 0.71235E+00 rms(broyden)= 0.71123E+00 rms(prec ) = 0.73130E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0229 5.7754 2.5967 1.1920 1.1920 0.7549 0.7549 0.5750 0.5114 0.5114 0.0751 0.2342 0.1296 0.0153 0.0031 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12497.93627052 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.50467681 PAW double counting = 13241.12706536 -12683.14492519 entropy T*S EENTRO = 0.02694840 eigenvalues EBANDS = -724.58821604 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -177.70509118 eV energy without entropy = -177.73203957 energy(sigma->0) = -177.71407398 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 1824 total energy-change (2. order) : 0.1161327E+02 (-0.1432196E+02) number of electron 135.9999987 magnetization 1.0477444 augmentation part -6.9094137 magnetization 1.0443137 Broyden mixing: rms(total) = 0.51454E+00 rms(broyden)= 0.51251E+00 rms(prec ) = 0.54089E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9742 5.7093 2.6753 1.2271 1.2271 0.7473 0.7473 0.5722 0.5267 0.5267 0.0751 0.2479 0.2479 0.0153 0.0031 0.0649 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12496.85021651 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.61618029 PAW double counting = 13186.30815702 -12628.79066421 entropy T*S EENTRO = 0.02714072 eigenvalues EBANDS = -713.48504524 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -166.09182489 eV energy without entropy = -166.11896561 energy(sigma->0) = -166.10087180 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 1869 total energy-change (2. order) :-0.6765330E+01 (-0.2834702E+01) number of electron 135.9999988 magnetization 0.6685750 augmentation part -6.8946680 magnetization 0.6596631 Broyden mixing: rms(total) = 0.22475E+00 rms(broyden)= 0.22446E+00 rms(prec ) = 0.25248E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0252 6.1488 3.0012 1.5699 0.9823 0.8595 0.8595 0.6565 0.6565 0.5244 0.5244 0.0751 0.2295 0.2295 0.0153 0.0031 0.0674 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12494.19856231 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.77167161 PAW double counting = 13078.52604166 -12520.26674834 entropy T*S EENTRO = 0.02718738 eigenvalues EBANDS = -723.48838562 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.85715522 eV energy without entropy = -172.88434260 energy(sigma->0) = -172.86621768 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 1866 total energy-change (2. order) : 0.9334455E+00 (-0.1909822E+00) number of electron 135.9999988 magnetization 0.4502415 augmentation part -6.9018913 magnetization 0.4465308 Broyden mixing: rms(total) = 0.21206E+00 rms(broyden)= 0.21204E+00 rms(prec ) = 0.23477E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0216 6.4036 2.8308 1.8228 1.0777 0.8810 0.8810 0.6769 0.6769 0.5346 0.5346 0.3823 0.2520 0.2520 0.0751 0.0153 0.0031 0.0671 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12486.31221090 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.24541172 PAW double counting = 12862.50505544 -12304.37062725 entropy T*S EENTRO = 0.02731872 eigenvalues EBANDS = -729.84281762 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.92370971 eV energy without entropy = -171.95102843 energy(sigma->0) = -171.93281595 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 3084 total energy-change (2. order) :-0.6704392E+02 (-0.4089340E+02) number of electron 135.9999985 magnetization 0.4447005 augmentation part -6.8614694 magnetization -0.0693505 Broyden mixing: rms(total) = 0.56391E+01 rms(broyden)= 0.55456E+01 rms(prec ) = 0.58886E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9650 6.4059 2.8300 1.8278 1.0728 0.8805 0.8805 0.6779 0.6779 0.5346 0.5346 0.3804 0.2523 0.2523 0.0751 0.0153 0.0031 0.0671 0.0011 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12482.33796160 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.41891945 PAW double counting = 12889.42996307 -12329.38987072 entropy T*S EENTRO = -0.00892739 eigenvalues EBANDS = -802.55690118 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -238.96763365 eV energy without entropy = -238.95870626 energy(sigma->0) = -238.96465785 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 2316 total energy-change (2. order) :-0.4557782E+03 (-0.1796517E+03) number of electron 135.9999978 magnetization 0.4439313 augmentation part -6.9972909 magnetization -0.9815643 Broyden mixing: rms(total) = 0.95635E+01 rms(broyden)= 0.94530E+01 rms(prec ) = 0.10238E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9138 6.4444 2.7610 1.8504 1.0024 0.8766 0.8766 0.7065 0.7065 0.5312 0.5312 0.4073 0.2489 0.2489 0.0751 0.0153 0.0671 0.0031 0.0029 0.0058 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12482.44347140 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.78512205 PAW double counting = 12889.62060997 -12328.49988941 entropy T*S EENTRO = -0.01898024 eigenvalues EBANDS = -1261.93399943 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -694.74586895 eV energy without entropy = -694.72688871 energy(sigma->0) = -694.73954220 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 2391 total energy-change (2. order) :-0.1058731E+04 (-0.1621577E+03) number of electron 135.9999973 magnetization 0.4404811 augmentation part -7.0443108 magnetization -0.2343289 Broyden mixing: rms(total) = 0.11139E+02 rms(broyden)= 0.11090E+02 rms(prec ) = 0.11838E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8679 6.4488 2.7465 1.8587 1.0196 0.8776 0.8776 0.6949 0.6949 0.5322 0.5322 0.4069 0.2494 0.2494 0.0751 0.0671 0.0153 0.0031 0.0060 0.0029 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12483.60627976 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -251.50987073 PAW double counting = 12892.60911547 -12330.10044465 entropy T*S EENTRO = -0.01944803 eigenvalues EBANDS = -2324.16537637 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1753.47732046 eV energy without entropy = -1753.45787243 energy(sigma->0) = -1753.47083778 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 2151 total energy-change (2. order) : 0.1584447E+04 (-0.9603065E+03) number of electron 135.9999988 magnetization 0.1345556 augmentation part -6.8871000 magnetization 0.1408601 Broyden mixing: rms(total) = 0.24036E+01 rms(broyden)= 0.52832E+00 rms(prec ) = 0.57019E+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.8639 6.5820 2.6603 2.2066 0.9801 0.9184 0.9184 0.7298 0.7298 0.5099 0.5099 0.5311 0.0751 0.2468 0.2264 0.2264 0.0671 0.0153 0.0031 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12484.65717780 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -248.43637876 PAW double counting = 12882.88224118 -12324.11762474 entropy T*S EENTRO = 0.02519659 eigenvalues EBANDS = -738.04126111 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.03002103 eV energy without entropy = -169.05521761 energy(sigma->0) = -169.03841989 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 1950 total energy-change (2. order) :-0.2679965E+01 (-0.2528617E+01) number of electron 135.9999988 magnetization 0.1710955 augmentation part -6.9054618 magnetization 0.1765110 Broyden mixing: rms(total) = 0.50093E+00 rms(broyden)= 0.21034E+00 rms(prec ) = 0.23574E+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.8741 6.6197 2.7274 2.3574 0.9753 0.9753 0.9084 0.9084 0.6177 0.5906 0.5906 0.5110 0.5110 0.2764 0.2467 0.2467 0.0751 0.0671 0.0153 0.0031 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12474.89352819 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.83044564 PAW double counting = 12778.65452550 -12220.37497695 entropy T*S EENTRO = 0.02793570 eigenvalues EBANDS = -741.60847963 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -171.70998559 eV energy without entropy = -171.73792130 energy(sigma->0) = -171.71929749 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 3081 total energy-change (2. order) :-0.3989886E+04 (-0.3660774E+04) number of electron 135.9999999 magnetization 0.1274282 augmentation part -6.7038894 magnetization 2.7347931 Broyden mixing: rms(total) = 0.10338E+02 rms(broyden)= 0.10208E+02 rms(prec ) = 0.10963E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8358 6.6618 2.5143 2.5143 0.9956 0.9956 0.9032 0.9032 0.6077 0.5805 0.5805 0.5092 0.5092 0.2826 0.2457 0.2457 0.0751 0.0671 0.0153 0.0031 0.0074 0.0038 0.0015 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12474.00203121 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.66519341 PAW double counting = 12798.32956866 -12236.19441825 entropy T*S EENTRO = -0.02231581 eigenvalues EBANDS = -4735.35657885 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -4161.59598525 eV energy without entropy = -4161.57366944 energy(sigma->0) = -4161.58854665 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 2067 total energy-change (2. order) : 0.3997347E+04 (-0.8505005E+03) number of electron 135.9999988 magnetization 0.0231305 augmentation part -6.8268787 magnetization -0.2241368 Broyden mixing: rms(total) = 0.14531E+01 rms(broyden)= 0.40442E+00 rms(prec ) = 0.41805E+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.8009 6.6642 2.5037 2.5037 0.9934 0.9934 0.9084 0.9084 0.5948 0.5855 0.5855 0.5061 0.5061 0.2818 0.2462 0.2462 0.0751 0.0671 0.0031 0.0153 0.0136 0.0136 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12474.99288021 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.37379330 PAW double counting = 12809.94259754 -12251.67544784 entropy T*S EENTRO = -0.00610750 eigenvalues EBANDS = -737.45810608 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -164.24875378 eV energy without entropy = -164.24264628 energy(sigma->0) = -164.24671795 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 2226 total energy-change (2. order) :-0.1174810E+02 (-0.3578508E+01) number of electron 135.9999988 magnetization 0.0896399 augmentation part -6.8967287 magnetization 0.1067939 Broyden mixing: rms(total) = 0.38543E+00 rms(broyden)= 0.26270E+00 rms(prec ) = 0.28762E+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.8209 6.7404 2.6175 2.6175 1.1164 1.1164 0.8799 0.8799 0.6572 0.5030 0.5030 0.5075 0.5075 0.4486 0.4486 0.3224 0.2407 0.2407 0.0751 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12473.74801207 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.34657239 PAW double counting = 13150.83194885 -12592.23146568 entropy T*S EENTRO = 0.02777671 eigenvalues EBANDS = -746.84550820 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -175.99684916 eV energy without entropy = -176.02462587 energy(sigma->0) = -176.00610806 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 2217 total energy-change (2. order) : 0.2777107E+01 (-0.3672735E+00) number of electron 135.9999988 magnetization 0.1754174 augmentation part -6.8990430 magnetization 0.1805003 Broyden mixing: rms(total) = 0.22608E+00 rms(broyden)= 0.21939E+00 rms(prec ) = 0.24494E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8586 6.7276 2.9087 2.9087 1.2416 1.2416 0.8176 0.8176 0.8197 0.6812 0.6812 0.6631 0.5473 0.5015 0.5015 0.3057 0.3057 0.2385 0.2385 0.0751 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12471.04537793 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.81175803 PAW double counting = 13162.68055924 -12604.38919328 entropy T*S EENTRO = 0.02775760 eigenvalues EBANDS = -745.99671379 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.21974257 eV energy without entropy = -173.24750017 energy(sigma->0) = -173.22899510 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 2253 total energy-change (2. order) :-0.7335465E+00 (-0.7222794E-01) number of electron 135.9999988 magnetization 0.1598448 augmentation part -6.9048303 magnetization 0.1608605 Broyden mixing: rms(total) = 0.21402E+00 rms(broyden)= 0.21374E+00 rms(prec ) = 0.23924E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8601 6.7438 3.0191 3.0191 1.2076 1.2076 0.9296 0.9296 0.8467 0.6877 0.6877 0.6472 0.5302 0.5302 0.4702 0.4702 0.3212 0.3212 0.2387 0.2387 0.0751 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12469.86655077 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.89374406 PAW double counting = 13093.77543907 -12535.47519278 entropy T*S EENTRO = 0.02765808 eigenvalues EBANDS = -747.83588224 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.95328909 eV energy without entropy = -173.98094717 energy(sigma->0) = -173.96250845 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 2163 total energy-change (2. order) :-0.5234552E+00 (-0.2258901E-01) number of electron 135.9999988 magnetization 0.1367982 augmentation part -6.9043561 magnetization 0.1373341 Broyden mixing: rms(total) = 0.22235E+00 rms(broyden)= 0.22231E+00 rms(prec ) = 0.24944E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8658 6.7474 3.0842 3.0842 1.0321 1.0321 1.0554 1.0554 0.8737 0.8737 0.7452 0.7452 0.5677 0.5677 0.5340 0.4875 0.4875 0.3072 0.3072 0.2386 0.2386 0.0751 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12469.99659845 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.88557774 PAW double counting = 13086.72925720 -12528.43110829 entropy T*S EENTRO = 0.02762668 eigenvalues EBANDS = -748.23532726 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.47674424 eV energy without entropy = -174.50437092 energy(sigma->0) = -174.48595314 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 2145 total energy-change (2. order) : 0.9107768E+00 (-0.6724863E-01) number of electron 135.9999988 magnetization 0.1122938 augmentation part -6.9065002 magnetization 0.1162640 Broyden mixing: rms(total) = 0.20185E+00 rms(broyden)= 0.20182E+00 rms(prec ) = 0.22568E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8729 6.7706 3.1032 3.1032 1.3950 1.3950 1.1172 1.1172 0.7828 0.7828 0.6559 0.6559 0.6480 0.6480 0.5095 0.5095 0.5351 0.0751 0.3135 0.3135 0.3059 0.2386 0.2386 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12469.65608225 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.88893590 PAW double counting = 13058.86984379 -12500.55977069 entropy T*S EENTRO = 0.02760233 eigenvalues EBANDS = -747.67360838 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.56596747 eV energy without entropy = -173.59356981 energy(sigma->0) = -173.57516825 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 2046 total energy-change (2. order) :-0.4674273E+00 (-0.1176357E-01) number of electron 135.9999988 magnetization 0.0497395 augmentation part -6.9056771 magnetization 0.0538532 Broyden mixing: rms(total) = 0.20050E+00 rms(broyden)= 0.20050E+00 rms(prec ) = 0.22567E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8987 6.8089 3.0860 3.0860 1.7159 1.7159 1.1913 1.1913 0.8298 0.8298 0.6715 0.6715 0.6844 0.6844 0.4836 0.4836 0.5278 0.5278 0.4955 0.3109 0.3109 0.2386 0.2386 0.0751 0.0671 0.0153 0.0031 0.0097 0.0015 0.0008 0.0038 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12469.22089990 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.88445571 PAW double counting = 13027.07197157 -12468.76741169 entropy T*S EENTRO = 0.02757824 eigenvalues EBANDS = -748.57516094 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.03339482 eV energy without entropy = -174.06097306 energy(sigma->0) = -174.04258756 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 2574 total energy-change (2. order) :-0.1710798E+02 (-0.5683635E+01) number of electron 135.9999989 magnetization 0.0440333 augmentation part -6.8744005 magnetization 0.0407368 Broyden mixing: rms(total) = 0.96356E+00 rms(broyden)= 0.96192E+00 rms(prec ) = 0.10198E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8700 6.8084 3.1047 3.1047 1.6971 1.6971 1.1711 1.1711 0.8431 0.8431 0.6647 0.6647 0.7067 0.7067 0.4836 0.4836 0.5476 0.4870 0.4870 0.3105 0.3105 0.2386 0.2386 0.0751 0.0671 0.0222 0.0153 0.0097 0.0031 0.0038 0.0008 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12468.15991438 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.88852314 PAW double counting = 12962.08807720 -12403.64990378 entropy T*S EENTRO = 0.02752323 eigenvalues EBANDS = -766.87362169 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.14137894 eV energy without entropy = -191.16890217 energy(sigma->0) = -191.15055335 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 2427 total energy-change (2. order) :-0.4199442E+02 (-0.2008576E+02) number of electron 135.9999983 magnetization 0.0414136 augmentation part -6.9045184 magnetization 0.8760075 Broyden mixing: rms(total) = 0.58538E+01 rms(broyden)= 0.57407E+01 rms(prec ) = 0.60880E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8431 6.8068 3.1075 3.1075 1.7007 1.7007 1.1699 1.1699 0.8444 0.8444 0.6644 0.6644 0.7068 0.7068 0.4834 0.4834 0.5463 0.4883 0.4883 0.3105 0.3105 0.2386 0.2386 0.0751 0.0671 0.0153 0.0203 0.0097 0.0038 0.0031 0.0008 0.0018 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12468.07520887 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.90467400 PAW double counting = 12955.50933220 -12397.02530802 entropy T*S EENTRO = -0.01188620 eigenvalues EBANDS = -808.94304004 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -233.13580131 eV energy without entropy = -233.12391511 energy(sigma->0) = -233.13183925 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 2496 total energy-change (2. order) : 0.3944244E+02 (-0.1351158E+03) number of electron 136.0000031 magnetization 0.0338711 augmentation part -6.8448638 magnetization -1.4767141 Broyden mixing: rms(total) = 0.26363E+01 rms(broyden)= 0.24466E+01 rms(prec ) = 0.25868E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8177 6.8060 3.1037 3.1037 1.7038 1.7038 1.1687 1.1687 0.8417 0.8417 0.6653 0.6653 0.7092 0.7092 0.4834 0.4834 0.5466 0.4868 0.4868 0.3105 0.3105 0.2386 0.2386 0.0751 0.0671 0.0239 0.0153 0.0097 0.0031 0.0038 0.0032 0.0032 0.0008 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12468.17100610 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.74148044 PAW double counting = 12956.17885022 -12397.69134176 entropy T*S EENTRO = -0.00151205 eigenvalues EBANDS = -772.58185027 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -193.69335678 eV energy without entropy = -193.69184473 energy(sigma->0) = -193.69285276 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 2217 total energy-change (2. order) :-0.9884043E+01 (-0.2752764E+02) number of electron 135.9999981 magnetization -0.0088901 augmentation part -6.8853891 magnetization 1.3742982 Broyden mixing: rms(total) = 0.48373E+01 rms(broyden)= 0.47785E+01 rms(prec ) = 0.51106E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7941 6.8056 3.1065 3.1065 1.7046 1.7046 1.1584 1.1584 0.8435 0.8435 0.6645 0.6645 0.7124 0.7124 0.4833 0.4833 0.5551 0.4856 0.4856 0.3105 0.3105 0.2386 0.2386 0.0751 0.0671 0.0295 0.0153 0.0078 0.0078 0.0031 0.0097 0.0015 0.0008 0.0038 0.0030 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12467.97148923 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.18901505 PAW double counting = 12950.27509654 -12391.50381493 entropy T*S EENTRO = -0.03056930 eigenvalues EBANDS = -780.47259152 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -203.57739987 eV energy without entropy = -203.54683057 energy(sigma->0) = -203.56721010 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 1848 total energy-change (2. order) : 0.3503876E+02 (-0.2401633E+03) number of electron 135.9999987 magnetization 0.0332621 augmentation part -6.9007398 magnetization -0.2309717 Broyden mixing: rms(total) = 0.74041E+00 rms(broyden)= 0.22005E+00 rms(prec ) = 0.24056E+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.7728 6.8121 3.0935 3.0935 1.6761 1.6761 1.1886 1.0899 0.8389 0.8389 0.6627 0.6627 0.7492 0.7492 0.5605 0.4836 0.4836 0.4828 0.4828 0.3105 0.3105 0.2386 0.2386 0.0751 0.0549 0.0549 0.0671 0.0153 0.0171 0.0171 0.0031 0.0097 0.0044 0.0038 0.0015 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12467.66959945 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.69333697 PAW double counting = 12941.10673613 -12382.80885837 entropy T*S EENTRO = 0.00351286 eigenvalues EBANDS = -746.79207434 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -168.53863653 eV energy without entropy = -168.54214939 energy(sigma->0) = -168.53980748 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 1935 total energy-change (2. order) : 0.2362711E+01 (-0.6174085E+02) number of electron 135.9999987 magnetization -0.0285580 augmentation part -6.9204360 magnetization -0.0356009 Broyden mixing: rms(total) = 0.56792E+00 rms(broyden)= 0.54888E+00 rms(prec ) = 0.56784E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7556 6.8118 3.1032 3.1032 1.6538 1.6538 1.1214 1.1214 0.8711 0.8711 0.6617 0.6617 0.7276 0.7276 0.5511 0.4838 0.4838 0.4824 0.4824 0.3105 0.3105 0.2386 0.2386 0.0751 0.1262 0.1262 0.0671 0.0534 0.0216 0.0216 0.0153 0.0031 0.0097 0.0015 0.0044 0.0038 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12466.51939480 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.72461132 PAW double counting = 13000.46710558 -12442.10020536 entropy T*S EENTRO = 0.03042466 eigenvalues EBANDS = -743.64422811 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -166.17592574 eV energy without entropy = -166.20635040 energy(sigma->0) = -166.18606729 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 1980 total energy-change (2. order) :-0.6580774E+01 (-0.5866071E+01) number of electron 135.9999988 magnetization 0.0322635 augmentation part -6.9083297 magnetization 0.0479111 Broyden mixing: rms(total) = 0.23789E+00 rms(broyden)= 0.23607E+00 rms(prec ) = 0.25164E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7463 6.8123 3.0678 3.0678 1.4877 1.4877 1.5002 1.0098 0.8436 0.8436 0.7325 0.7325 0.6679 0.6679 0.5857 0.4872 0.4872 0.4688 0.4688 0.3445 0.3445 0.3100 0.3100 0.2386 0.2386 0.0751 0.0928 0.0928 0.0671 0.0215 0.0215 0.0153 0.0097 0.0031 0.0044 0.0038 0.0008 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12465.84013232 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.90194875 PAW double counting = 12946.77813718 -12388.46549755 entropy T*S EENTRO = 0.02696094 eigenvalues EBANDS = -750.66920293 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.75669981 eV energy without entropy = -172.78366075 energy(sigma->0) = -172.76568679 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 2010 total energy-change (2. order) : 0.2858843E+01 (-0.1038153E+01) number of electron 135.9999988 magnetization -0.0009684 augmentation part -6.9142325 magnetization 0.0110729 Broyden mixing: rms(total) = 0.31299E+00 rms(broyden)= 0.31289E+00 rms(prec ) = 0.32786E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7386 6.8114 3.0750 3.0750 1.4853 1.4853 1.4844 1.0752 0.8553 0.8553 0.6703 0.6703 0.6914 0.6914 0.4845 0.4845 0.5736 0.4863 0.4863 0.4787 0.4787 0.3105 0.3105 0.2386 0.2386 0.0751 0.1297 0.1087 0.1087 0.0671 0.0215 0.0215 0.0153 0.0031 0.0097 0.0044 0.0038 0.0015 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12465.83814380 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.95363649 PAW double counting = 12976.62030321 -12418.36096547 entropy T*S EENTRO = 0.02725718 eigenvalues EBANDS = -747.70765541 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.89785716 eV energy without entropy = -169.92511434 energy(sigma->0) = -169.90694288 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 2001 total energy-change (2. order) :-0.4072442E+01 (-0.7431926E+00) number of electron 135.9999988 magnetization -0.0092275 augmentation part -6.9073110 magnetization -0.0000685 Broyden mixing: rms(total) = 0.19981E+00 rms(broyden)= 0.19957E+00 rms(prec ) = 0.22036E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7488 6.7986 3.0839 3.0839 1.5746 1.5746 1.4826 1.0044 0.7198 0.7198 0.8508 0.8508 0.7333 0.7333 0.6605 0.6605 0.5786 0.4871 0.4871 0.4457 0.4457 0.3462 0.3462 0.3081 0.3081 0.2386 0.2386 0.0751 0.1085 0.1085 0.0671 0.0215 0.0215 0.0153 0.0097 0.0031 0.0044 0.0038 0.0008 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.74310293 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94245260 PAW double counting = 12919.57170089 -12361.28111439 entropy T*S EENTRO = 0.02729903 eigenvalues EBANDS = -752.91761282 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.97029921 eV energy without entropy = -173.99759824 energy(sigma->0) = -173.97939889 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 2190 total energy-change (2. order) :-0.8923122E+01 (-0.2256274E+01) number of electron 135.9999988 magnetization -0.0105029 augmentation part -6.8894979 magnetization -0.0092622 Broyden mixing: rms(total) = 0.64140E+00 rms(broyden)= 0.64049E+00 rms(prec ) = 0.68507E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7591 6.7630 3.1425 3.1425 1.5435 1.5435 1.3996 1.0140 1.0140 0.9726 0.9726 0.7198 0.7198 0.6939 0.6939 0.7211 0.6295 0.4883 0.4883 0.5335 0.4940 0.4940 0.3626 0.3137 0.3137 0.2718 0.2387 0.2387 0.0751 0.1084 0.1084 0.0671 0.0215 0.0215 0.0153 0.0031 0.0097 0.0044 0.0038 0.0008 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.22240143 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94077476 PAW double counting = 12887.99243990 -12329.15242984 entropy T*S EENTRO = 0.02733461 eigenvalues EBANDS = -762.91257362 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -182.89342152 eV energy without entropy = -182.92075613 energy(sigma->0) = -182.90253306 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 1797 total energy-change (2. order) : 0.1337272E+02 (-0.1672835E+03) number of electron 135.9999988 magnetization -0.0154130 augmentation part -6.9151396 magnetization -0.0555093 Broyden mixing: rms(total) = 0.63003E+00 rms(broyden)= 0.62964E+00 rms(prec ) = 0.63574E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7417 6.7586 3.1263 3.1263 1.5530 1.5530 1.2972 1.0477 1.0477 1.0216 1.0216 0.7138 0.7138 0.7213 0.6985 0.6985 0.6263 0.5308 0.4951 0.4951 0.4725 0.4725 0.3547 0.3134 0.3134 0.2680 0.2387 0.2387 0.0751 0.1085 0.1085 0.0671 0.0502 0.0215 0.0215 0.0153 0.0031 0.0097 0.0015 0.0044 0.0038 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.39311933 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.96096542 PAW double counting = 12894.15768868 -12335.50535127 entropy T*S EENTRO = 0.03568576 eigenvalues EBANDS = -749.16962610 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.52070407 eV energy without entropy = -169.55638983 energy(sigma->0) = -169.53259932 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 2358 total energy-change (2. order) :-0.1431817E+02 (-0.8504718E+01) number of electron 135.9999989 magnetization -0.0269522 augmentation part -6.8814342 magnetization -0.0425908 Broyden mixing: rms(total) = 0.78640E+00 rms(broyden)= 0.78573E+00 rms(prec ) = 0.83619E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7277 6.7644 3.1065 3.1065 1.5461 1.5461 1.3211 1.0796 1.0796 1.0183 1.0183 0.7082 0.7082 0.7138 0.7138 0.6570 0.6570 0.5477 0.4961 0.4961 0.4580 0.4580 0.3527 0.3132 0.3132 0.2633 0.2387 0.2387 0.0751 0.1021 0.1021 0.1084 0.1084 0.0671 0.0215 0.0215 0.0153 0.0031 0.0097 0.0015 0.0044 0.0038 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.26202552 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.95622632 PAW double counting = 12875.18695834 -12316.21507679 entropy T*S EENTRO = 0.03163694 eigenvalues EBANDS = -763.93912197 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -183.83887169 eV energy without entropy = -183.87050863 energy(sigma->0) = -183.84941734 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 2490 total energy-change (2. order) :-0.4593559E+03 (-0.4006101E+03) number of electron 136.0000017 magnetization -0.0061021 augmentation part -6.8691858 magnetization -1.2199750 Broyden mixing: rms(total) = 0.88606E+01 rms(broyden)= 0.86860E+01 rms(prec ) = 0.94806E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7108 6.7647 3.1061 3.1061 1.5474 1.5474 1.3149 1.0795 1.0795 1.0174 1.0174 0.7080 0.7080 0.7138 0.7138 0.6600 0.6600 0.5464 0.4961 0.4961 0.4577 0.4577 0.3519 0.3132 0.3132 0.2634 0.2387 0.2387 0.0751 0.1024 0.1024 0.1084 0.1084 0.0671 0.0215 0.0215 0.0153 0.0097 0.0031 0.0044 0.0038 0.0015 0.0008 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.40202815 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94985704 PAW double counting = 12874.88900022 -12314.78998540 entropy T*S EENTRO = -0.00923363 eigenvalues EBANDS = -1224.24760728 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -643.19472766 eV energy without entropy = -643.18549403 energy(sigma->0) = -643.19164979 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 2448 total energy-change (2. order) :-0.1784346E+04 (-0.6246830E+03) number of electron 135.9999973 magnetization -0.0887754 augmentation part -7.0401500 magnetization -0.0281822 Broyden mixing: rms(total) = 0.11017E+02 rms(broyden)= 0.10930E+02 rms(prec ) = 0.11673E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6949 6.7699 3.1042 3.1042 1.5328 1.5328 1.2011 1.2011 1.0900 1.0900 0.8251 0.8251 0.7174 0.7174 0.6989 0.6989 0.5795 0.5795 0.4956 0.4956 0.4623 0.4623 0.3559 0.3134 0.3134 0.2646 0.2387 0.2387 0.0751 0.1084 0.1084 0.0976 0.0976 0.0671 0.0215 0.0215 0.0153 0.0165 0.0165 0.0031 0.0097 0.0015 0.0044 0.0038 0.0008 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12465.54464251 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.63716795 PAW double counting = 12885.52091462 -12324.24021083 entropy T*S EENTRO = -0.01224922 eigenvalues EBANDS = -3012.94199256 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -2427.54036484 eV energy without entropy = -2427.52811561 energy(sigma->0) = -2427.53628176 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 2241 total energy-change (2. order) : 0.1834453E+04 (-0.9598856E+02) number of electron 135.9999996 magnetization -0.0797968 augmentation part -7.0060523 magnetization 1.2175093 Broyden mixing: rms(total) = 0.84588E+01 rms(broyden)= 0.83874E+01 rms(prec ) = 0.90095E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6448 6.5754 3.3432 2.2256 1.6311 1.6311 0.9437 0.9437 0.8452 0.8452 0.6502 0.6502 0.5790 0.5790 0.5907 0.5907 0.3846 0.3846 0.3686 0.3686 0.3619 0.2530 0.2530 0.1521 0.1521 0.1125 0.1125 0.1104 0.0599 0.0124 0.0165 0.0165 0.0162 0.0083 0.0083 0.0026 0.0014 0.0051 0.0035 0.0006 0.0016 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12467.55809100 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -248.84733154 PAW double counting = 12891.83552455 -12331.57454005 entropy T*S EENTRO = -0.00900930 eigenvalues EBANDS = -1180.24911445 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -593.08757817 eV energy without entropy = -593.07856887 energy(sigma->0) = -593.08457507 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 1884 total energy-change (2. order) : 0.4200513E+03 (-0.3943206E+03) number of electron 135.9999988 magnetization -0.0501695 augmentation part -6.8811274 magnetization -0.0429576 Broyden mixing: rms(total) = 0.21636E+01 rms(broyden)= 0.67974E+00 rms(prec ) = 0.72526E+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.6293 6.5773 3.3677 2.2153 1.6320 1.6320 0.9648 0.9648 0.7885 0.7885 0.6164 0.6164 0.5959 0.5959 0.5900 0.5900 0.3906 0.3906 0.3795 0.3795 0.3538 0.2466 0.2466 0.0967 0.0967 0.1467 0.1467 0.1266 0.0621 0.0621 0.0244 0.0244 0.0250 0.0250 0.0159 0.0095 0.0030 0.0012 0.0047 0.0039 0.0010 0.0014 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12466.62143196 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -250.24357760 PAW double counting = 12870.01097747 -12311.02537800 entropy T*S EENTRO = 0.02537440 eigenvalues EBANDS = -758.49718971 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.03624179 eV energy without entropy = -173.06161619 energy(sigma->0) = -173.04469992 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- eigenvalue-minimisations : 1839 total energy-change (2. order) : 0.8842145E+01 (-0.2045608E+02) number of electron 135.9999987 magnetization -0.0983709 augmentation part -6.9252505 magnetization -0.0843640 Broyden mixing: rms(total) = 0.63151E+00 rms(broyden)= 0.50522E+00 rms(prec ) = 0.52953E+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.6320 6.5686 3.4280 2.0631 1.6274 1.6274 0.9708 0.9708 0.7530 0.7530 0.3873 0.3873 0.6805 0.6805 0.5703 0.5703 0.6051 0.6051 0.4034 0.4034 0.3981 0.3981 0.3620 0.2545 0.2545 0.1632 0.1632 0.1077 0.0846 0.0846 0.0617 0.0444 0.0444 0.0145 0.0145 0.0161 0.0087 0.0019 0.0019 0.0041 0.0032 0.0009 0.0014 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12465.80171682 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.15561550 PAW double counting = 12975.68913620 -12418.39218834 entropy T*S EENTRO = 0.02750024 eigenvalues EBANDS = -741.87619612 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -164.19409674 eV energy without entropy = -164.22159697 energy(sigma->0) = -164.20326348 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 53) --------------------------------------- eigenvalue-minimisations : 2082 total energy-change (2. order) :-0.9483890E+01 (-0.3210157E+01) number of electron 135.9999988 magnetization -0.1385520 augmentation part -6.9113919 magnetization -0.1238380 Broyden mixing: rms(total) = 0.22019E+00 rms(broyden)= 0.20254E+00 rms(prec ) = 0.21683E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6441 6.5635 3.4151 2.2413 1.7060 1.7060 0.8024 1.0053 1.0053 0.7670 0.7670 0.6843 0.6843 0.6203 0.6203 0.4918 0.4918 0.3313 0.3313 0.4080 0.4080 0.4745 0.4745 0.3215 0.2537 0.2537 0.2677 0.1323 0.1323 0.0986 0.0389 0.0620 0.0254 0.0254 0.0217 0.0217 0.0162 0.0034 0.0091 0.0012 0.0051 0.0040 0.0014 0.0009 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.50339204 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.87490445 PAW double counting = 12906.81139918 -12348.57898801 entropy T*S EENTRO = 0.02755426 eigenvalues EBANDS = -753.87463952 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.67798697 eV energy without entropy = -173.70554123 energy(sigma->0) = -173.68717172 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 54) --------------------------------------- eigenvalue-minimisations : 2085 total energy-change (2. order) :-0.1632826E+01 (-0.4003679E+00) number of electron 135.9999988 magnetization -0.1187083 augmentation part -6.9065337 magnetization -0.1036001 Broyden mixing: rms(total) = 0.29365E+00 rms(broyden)= 0.29282E+00 rms(prec ) = 0.31352E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6438 6.5952 3.1848 2.5720 1.6280 1.6280 1.0743 1.0743 0.6155 0.6155 0.8154 0.8154 0.7413 0.4588 0.4588 0.5128 0.5128 0.5585 0.5585 0.4570 0.4570 0.4720 0.4720 0.3171 0.3171 0.2619 0.2619 0.2044 0.1573 0.1573 0.0884 0.0576 0.0576 0.0583 0.0412 0.0146 0.0146 0.0170 0.0079 0.0023 0.0018 0.0042 0.0038 0.0007 0.0014 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.87777517 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.96673761 PAW double counting = 12893.89203616 -12335.65293376 entropy T*S EENTRO = 0.02747496 eigenvalues EBANDS = -756.04786133 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -175.31081313 eV energy without entropy = -175.33828808 energy(sigma->0) = -175.31997144 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 55) --------------------------------------- eigenvalue-minimisations : 2235 total energy-change (2. order) : 0.9430203E+00 (-0.4558096E-01) number of electron 135.9999988 magnetization -0.0919762 augmentation part -6.9077929 magnetization -0.0793821 Broyden mixing: rms(total) = 0.23869E+00 rms(broyden)= 0.23866E+00 rms(prec ) = 0.25809E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6163 5.3961 2.6261 2.6261 1.3293 1.3293 0.8779 0.8779 0.6305 0.6305 0.8862 0.7293 0.7293 0.5369 0.5369 0.4289 0.4289 0.6173 0.5640 0.4919 0.4513 0.3293 0.3293 0.3438 0.2001 0.2001 0.1969 0.0580 0.0580 0.0577 0.0487 0.0331 0.0331 0.0138 0.0138 0.0069 0.0001 0.0001 0.0014 0.0021 0.0021 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.46216539 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.96240571 PAW double counting = 12895.23105119 -12337.02425314 entropy T*S EENTRO = 0.02749903 eigenvalues EBANDS = -754.49250244 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.36779285 eV energy without entropy = -174.39529188 energy(sigma->0) = -174.37695920 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 56) --------------------------------------- eigenvalue-minimisations : 2175 total energy-change (2. order) :-0.1809667E+01 (-0.1713448E+00) number of electron 135.9999988 magnetization -0.0759608 augmentation part -6.9022875 magnetization -0.0680339 Broyden mixing: rms(total) = 0.33900E+00 rms(broyden)= 0.33890E+00 rms(prec ) = 0.35967E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6274 5.4550 2.6557 2.6557 1.3995 1.3995 0.6373 0.6373 0.8831 0.8831 0.9598 0.7235 0.7235 0.5394 0.5394 0.7140 0.4236 0.4236 0.6219 0.4936 0.4894 0.4894 0.3695 0.3695 0.2585 0.2585 0.1901 0.1901 0.0526 0.0526 0.0567 0.0567 0.0405 0.0405 0.0141 0.0141 0.0065 0.0040 0.0011 0.0011 0.0001 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.57447771 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94436392 PAW double counting = 12782.89187087 -12224.67076190 entropy T*S EENTRO = 0.02743422 eigenvalues EBANDS = -757.22214455 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -176.17745938 eV energy without entropy = -176.20489359 energy(sigma->0) = -176.18660412 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 57) --------------------------------------- eigenvalue-minimisations : 1914 total energy-change (2. order) : 0.2883258E+01 (-0.4857646E+00) number of electron 135.9999988 magnetization -0.0577214 augmentation part -6.9085932 magnetization -0.0546685 Broyden mixing: rms(total) = 0.17185E+00 rms(broyden)= 0.17151E+00 rms(prec ) = 0.18894E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6404 5.4733 2.6508 2.6508 1.4762 1.3303 1.3303 0.9003 0.9003 0.6326 0.6326 0.8563 0.7251 0.7251 0.5467 0.5467 0.4219 0.4219 0.6440 0.5774 0.5346 0.4817 0.4817 0.3600 0.3600 0.2476 0.2476 0.2083 0.1922 0.0579 0.0579 0.0557 0.0557 0.0374 0.0374 0.0140 0.0140 0.0001 0.0000 0.0011 0.0011 0.0040 0.0065 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.84777373 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94643423 PAW double counting = 12761.42291892 -12203.28025459 entropy T*S EENTRO = 0.02742997 eigenvalues EBANDS = -753.98507089 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.29420093 eV energy without entropy = -173.32163091 energy(sigma->0) = -173.30334426 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 58) --------------------------------------- eigenvalue-minimisations : 2037 total energy-change (2. order) :-0.1097154E-01 (-0.3222860E-02) number of electron 135.9999988 magnetization -0.0588543 augmentation part -6.9088600 magnetization -0.0576354 Broyden mixing: rms(total) = 0.17305E+00 rms(broyden)= 0.17304E+00 rms(prec ) = 0.19000E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6538 5.4859 2.6482 2.6482 1.9730 1.3053 1.3053 0.8921 0.8921 0.9632 0.6176 0.6176 0.7505 0.7505 0.4429 0.4429 0.5041 0.5041 0.6745 0.6081 0.6081 0.5245 0.4696 0.4696 0.3827 0.3827 0.2472 0.2472 0.2164 0.2015 0.0539 0.0539 0.0557 0.0557 0.0387 0.0387 0.0143 0.0143 0.0001 0.0000 0.0011 0.0011 0.0065 0.0040 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.87358951 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94228951 PAW double counting = 12738.69544477 -12180.55074486 entropy T*S EENTRO = 0.02741334 eigenvalues EBANDS = -753.97639032 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.30517247 eV energy without entropy = -173.33258581 energy(sigma->0) = -173.31431025 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 59) --------------------------------------- eigenvalue-minimisations : 2235 total energy-change (2. order) : 0.1617927E+00 (-0.1557770E-02) number of electron 135.9999988 magnetization -0.0692746 augmentation part -6.9094435 magnetization -0.0678768 Broyden mixing: rms(total) = 0.17528E+00 rms(broyden)= 0.17528E+00 rms(prec ) = 0.19137E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6701 5.4873 2.5785 2.5785 2.4897 1.2816 1.2816 1.0253 1.0253 0.8677 0.8677 0.6268 0.6268 0.7506 0.7506 0.4623 0.4623 0.4777 0.4777 0.6445 0.6445 0.5797 0.5521 0.4800 0.4800 0.3898 0.3898 0.2338 0.2338 0.1976 0.1976 0.0574 0.0574 0.0602 0.0540 0.0369 0.0369 0.0142 0.0142 0.0064 0.0043 0.0011 0.0011 0.0000 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.04159597 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93826074 PAW double counting = 12725.67925339 -12167.55622375 entropy T*S EENTRO = 0.02738305 eigenvalues EBANDS = -753.62891944 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.14337982 eV energy without entropy = -173.17076287 energy(sigma->0) = -173.15250751 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 60) --------------------------------------- eigenvalue-minimisations : 1758 total energy-change (2. order) : 0.6800247E+01 (-0.4613110E+03) number of electron 135.9999987 magnetization -0.0713715 augmentation part -6.9309826 magnetization 0.0529189 Broyden mixing: rms(total) = 0.54982E+00 rms(broyden)= 0.54961E+00 rms(prec ) = 0.57163E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5143 3.2025 2.3976 1.9982 0.9278 0.9278 0.9647 0.9647 0.6306 0.6306 0.8905 0.6727 0.6727 0.6304 0.6304 0.5567 0.5567 0.5805 0.3668 0.3668 0.4314 0.4314 0.3269 0.1525 0.1525 0.1200 0.1200 0.0777 0.0485 0.0485 0.0294 0.0294 0.0110 0.0110 0.0067 0.0039 0.0010 0.0010 0.0012 0.0012 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.05292836 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93587334 PAW double counting = 12720.87232500 -12161.94845808 entropy T*S EENTRO = 0.00314688 eigenvalues EBANDS = -747.59632815 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -166.34313243 eV energy without entropy = -166.34627931 energy(sigma->0) = -166.34418139 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 61) --------------------------------------- eigenvalue-minimisations : 1755 total energy-change (2. order) :-0.6746452E+01 (-0.6631484E+01) number of electron 135.9999988 magnetization -0.0687463 augmentation part -6.9082899 magnetization -0.0662103 Broyden mixing: rms(total) = 0.21152E+00 rms(broyden)= 0.21108E+00 rms(prec ) = 0.22710E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5090 3.2235 2.4240 1.9975 0.9596 0.9596 0.9751 0.9751 0.6153 0.6153 0.8908 0.6690 0.6690 0.6416 0.6416 0.5535 0.5535 0.5756 0.4651 0.4651 0.3735 0.3735 0.2848 0.2029 0.2029 0.1361 0.1361 0.0921 0.0517 0.0269 0.0260 0.0221 0.0221 0.0156 0.0156 0.0060 0.0060 0.0012 0.0011 0.0011 0.0000 0.0002 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.97383201 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93422100 PAW double counting = 12712.81934356 -12154.67667641 entropy T*S EENTRO = 0.02738065 eigenvalues EBANDS = -753.66656321 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.08958480 eV energy without entropy = -173.11696545 energy(sigma->0) = -173.09871169 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 62) --------------------------------------- eigenvalue-minimisations : 1944 total energy-change (2. order) : 0.2046924E+00 (-0.2440820E+00) number of electron 135.9999988 magnetization -0.0732060 augmentation part -6.9099695 magnetization -0.0712180 Broyden mixing: rms(total) = 0.17746E+00 rms(broyden)= 0.17740E+00 rms(prec ) = 0.19233E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5253 3.3735 2.4931 1.9848 1.0841 0.9842 0.9842 0.9562 0.9562 0.6077 0.6077 0.7084 0.7084 0.5217 0.5217 0.6255 0.6255 0.5288 0.5288 0.5042 0.4820 0.4820 0.3400 0.3400 0.2692 0.1936 0.1936 0.0850 0.0850 0.0798 0.0461 0.0461 0.0364 0.0364 0.0121 0.0121 0.0056 0.0068 0.0030 0.0012 0.0012 0.0001 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.13133779 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93463611 PAW double counting = 12715.10923384 -12157.02109344 entropy T*S EENTRO = 0.02740313 eigenvalues EBANDS = -753.24944566 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.88489240 eV energy without entropy = -172.91229553 energy(sigma->0) = -172.89402678 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 63) --------------------------------------- eigenvalue-minimisations : 1989 total energy-change (2. order) :-0.2254420E+00 (-0.6526781E-02) number of electron 135.9999988 magnetization -0.0843945 augmentation part -6.9094267 magnetization -0.0826186 Broyden mixing: rms(total) = 0.17318E+00 rms(broyden)= 0.17318E+00 rms(prec ) = 0.18914E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5431 3.3968 2.6253 1.9870 1.2565 1.0831 1.0831 1.0538 1.0538 0.6121 0.6121 0.7162 0.7162 0.5584 0.5584 0.6322 0.6322 0.5553 0.5553 0.5317 0.5317 0.4192 0.4192 0.3336 0.3336 0.2721 0.1913 0.1913 0.0747 0.0747 0.0795 0.0467 0.0467 0.0390 0.0390 0.0114 0.0114 0.0058 0.0067 0.0034 0.0012 0.0012 0.0002 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12463.06712450 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.92785562 PAW double counting = 12711.64770844 -12153.52390438 entropy T*S EENTRO = 0.02739815 eigenvalues EBANDS = -753.58154005 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.11033435 eV energy without entropy = -173.13773250 energy(sigma->0) = -173.11946707 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 64) --------------------------------------- eigenvalue-minimisations : 1983 total energy-change (2. order) :-0.1981740E+00 (-0.3805518E-02) number of electron 135.9999988 magnetization -0.0913287 augmentation part -6.9090519 magnetization -0.0896915 Broyden mixing: rms(total) = 0.16901E+00 rms(broyden)= 0.16900E+00 rms(prec ) = 0.18593E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5621 3.4337 2.6127 1.9906 1.8430 0.7435 0.7435 1.1165 0.9595 0.9595 1.0098 0.6974 0.6974 0.6900 0.6900 0.4049 0.4049 0.4749 0.4749 0.6426 0.6426 0.5552 0.5552 0.4441 0.4099 0.4099 0.2552 0.2552 0.1537 0.1537 0.0542 0.0542 0.0523 0.0364 0.0364 0.0306 0.0129 0.0129 0.0067 0.0070 0.0016 0.0016 0.0001 0.0001 0.0007 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.87138726 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.92907576 PAW double counting = 12706.02630208 -12147.86858926 entropy T*S EENTRO = 0.02739171 eigenvalues EBANDS = -754.00813346 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.30850833 eV energy without entropy = -173.33590004 energy(sigma->0) = -173.31763890 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 65) --------------------------------------- eigenvalue-minimisations : 1986 total energy-change (2. order) :-0.1468492E+00 (-0.1520704E-02) number of electron 135.9999988 magnetization -0.0115922 augmentation part -6.9088270 magnetization -0.0099092 Broyden mixing: rms(total) = 0.16911E+00 rms(broyden)= 0.16911E+00 rms(prec ) = 0.18660E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4876 1.8820 1.8820 1.5018 0.6090 0.6090 0.9912 0.9912 0.9815 0.7066 0.7066 0.7753 0.7753 0.5245 0.5245 0.6979 0.6348 0.6348 0.5639 0.5639 0.5815 0.5406 0.3734 0.3097 0.3097 0.2060 0.2060 0.1363 0.1285 0.0479 0.0235 0.0235 0.0140 0.0140 0.0118 0.0118 0.0035 0.0035 0.0002 0.0014 0.0014 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12462.76438453 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.93024342 PAW double counting = 12705.96282908 -12147.77905358 entropy T*S EENTRO = 0.02738643 eigenvalues EBANDS = -754.28687513 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.45535754 eV energy without entropy = -173.48274397 energy(sigma->0) = -173.46448635 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 66) --------------------------------------- eigenvalue-minimisations : 2040 total energy-change (2. order) :-0.1230662E+01 (-0.9785158E-01) number of electron 135.9999988 magnetization -0.0837047 augmentation part -6.9060241 magnetization -0.0816619 Broyden mixing: rms(total) = 0.23795E+00 rms(broyden)= 0.23790E+00 rms(prec ) = 0.25847E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5077 2.0939 1.8143 1.8143 1.2740 0.6067 0.6067 0.7134 0.7134 0.8950 0.8520 0.8520 0.8027 0.8027 0.6989 0.6117 0.6117 0.4706 0.4706 0.5711 0.5711 0.5982 0.5415 0.3683 0.3237 0.3237 0.1400 0.2021 0.2021 0.1120 0.0455 0.0236 0.0236 0.0134 0.0134 0.0126 0.0126 0.0043 0.0043 0.0004 0.0013 0.0013 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8963.22719592 -Hartree energ DENC = -12464.14492647 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.95067929 PAW double counting = 12770.19859784 -12211.95571392 entropy T*S EENTRO = 0.02755648 eigenvalues EBANDS = -754.17583729 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.68601904 eV energy without entropy = -174.71357551 energy(sigma->0) = -174.69520453 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 67) --------------------------------------- ----------------------------------------------------------------------------- | | | 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!!! <---- | | | -----------------------------------------------------------------------------