vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.08.31 13:27:59 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.931 0.528 0.108- 19 2.31 14 2.39 3 2.42 2 0.850 0.233 0.428- 4 2.34 13 2.38 20 2.44 3 0.137 0.448 0.195- 8 2.31 17 2.35 7 2.37 1 2.42 4 0.093 0.282 0.333- 18 2.29 2 2.34 7 2.35 8 2.66 5 0.818 0.333 0.026- 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.40 11 2.42 9 2.61 7 0.081 0.333 0.137- 21 2.34 5 2.35 4 2.35 3 2.37 8 0.198 0.407 0.390- 3 2.31 6 2.33 22 2.38 18 2.60 4 2.66 9 0.831 0.545 0.521- 11 2.36 27 2.36 6 2.61 10 0.771 0.244 0.883- 5 2.37 28 2.39 12 2.40 31 2.46 11 0.061 0.510 0.653- 16 2.36 9 2.36 6 2.42 25 2.43 12 0.013 0.205 0.761- 26 2.34 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.444 0.955- 16 2.33 5 2.34 32 2.35 1 2.39 15 0.138 0.323 0.737- 29 2.33 16 2.38 13 2.40 26 2.51 12 2.52 16 0.098 0.434 0.819- 14 2.33 30 2.34 11 2.36 15 2.38 17 0.424 0.460 0.121- 19 2.30 30 2.32 3 2.35 23 2.48 18 0.363 0.292 0.421- 4 2.29 29 2.33 20 2.35 8 2.60 19 0.690 0.484 0.209- 17 2.30 1 2.31 24 2.34 20 0.643 0.313 0.341- 18 2.35 23 2.35 24 2.38 2 2.44 21 0.324 0.316 0.009- 7 2.34 26 2.34 23 2.38 30 2.38 22 0.435 0.438 0.521- 24 2.38 8 2.38 29 2.44 27 2.47 25 2.49 23 0.561 0.346 0.141- 5 2.34 20 2.35 21 2.38 17 2.48 24 0.693 0.429 0.400- 19 2.34 22 2.38 20 2.38 6 2.40 25 0.319 0.555 0.551- 35 1.80 11 2.43 22 2.49 27 2.55 26 0.273 0.231 0.864- 21 2.34 12 2.34 28 2.40 15 2.51 27 0.600 0.517 0.656- 35 1.67 32 2.36 9 2.36 22 2.47 25 2.55 28 0.538 0.209 0.751- 10 2.39 26 2.40 31 2.52 29 0.396 0.329 0.624- 18 2.33 15 2.33 31 2.38 22 2.44 30 0.368 0.426 0.920- 32 2.32 17 2.32 16 2.34 21 2.38 31 0.670 0.326 0.727- 13 2.32 29 2.38 32 2.38 10 2.46 28 2.52 32 0.631 0.436 0.815- 30 2.32 14 2.35 27 2.36 31 2.38 33 0.654 0.705 0.622- 34 0.72 34 0.671 0.717 0.561- 33 0.72 35 0.459 0.580 0.677- 27 1.67 25 1.80 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.930880880 0.527973460 0.108085240 0.849679860 0.232506810 0.427890010 0.136890400 0.447556120 0.194880590 0.092573260 0.282090830 0.332529210 0.817969700 0.332898830 0.026361630 0.956583240 0.421118100 0.518185960 0.081331080 0.333153270 0.136830110 0.198451900 0.406825800 0.390163080 0.831374990 0.544730400 0.521414410 0.771308040 0.243660840 0.882657700 0.061083380 0.510443190 0.652930850 0.012946930 0.205110760 0.760956220 0.902110570 0.309795080 0.592636740 0.864308600 0.443844800 0.955429590 0.137572490 0.323214710 0.736823100 0.098272960 0.434374430 0.819307590 0.424021670 0.459500650 0.121432190 0.362601340 0.292448920 0.421102590 0.690397190 0.484334400 0.209328280 0.642613110 0.313456040 0.341409350 0.323923980 0.316497440 0.009488670 0.435081280 0.438310990 0.521100270 0.561049100 0.345937490 0.140772700 0.692504610 0.428584550 0.400281860 0.318911070 0.555358010 0.551209090 0.273423920 0.230913370 0.863956940 0.600034640 0.516529540 0.656357880 0.537884650 0.208965920 0.751386060 0.396378460 0.328744720 0.624049620 0.368036830 0.425626020 0.919841810 0.669879730 0.325718960 0.727170860 0.630757280 0.435756690 0.815482770 0.654008740 0.705113800 0.622391660 0.671434980 0.717486530 0.561097900 0.459179480 0.580090690 0.676949750 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.93088088 0.52797346 0.10808524 0.84967986 0.23250681 0.42789001 0.13689040 0.44755612 0.19488059 0.09257326 0.28209083 0.33252921 0.81796970 0.33289883 0.02636163 0.95658324 0.42111810 0.51818596 0.08133108 0.33315327 0.13683011 0.19845190 0.40682580 0.39016308 0.83137499 0.54473040 0.52141441 0.77130804 0.24366084 0.88265770 0.06108338 0.51044319 0.65293085 0.01294693 0.20511076 0.76095622 0.90211057 0.30979508 0.59263674 0.86430860 0.44384480 0.95542959 0.13757249 0.32321471 0.73682310 0.09827296 0.43437443 0.81930759 0.42402167 0.45950065 0.12143219 0.36260134 0.29244892 0.42110259 0.69039719 0.48433440 0.20932828 0.64261311 0.31345604 0.34140935 0.32392398 0.31649744 0.00948867 0.43508128 0.43831099 0.52110027 0.56104910 0.34593749 0.14077270 0.69250461 0.42858455 0.40028186 0.31891107 0.55535801 0.55120909 0.27342392 0.23091337 0.86395694 0.60003464 0.51652954 0.65635788 0.53788465 0.20896592 0.75138606 0.39637846 0.32874472 0.62404962 0.36803683 0.42562602 0.91984181 0.66987973 0.32571896 0.72717086 0.63075728 0.43575669 0.81548277 0.65400874 0.70511380 0.62239166 0.67143498 0.71748653 0.56109790 0.45917948 0.58009069 0.67694975 position of ions in cartesian coordinates (Angst): 7.13343327 10.38159494 1.17134785 6.51118174 4.57180466 4.63715529 1.04900482 8.80034074 2.11197162 0.70939815 5.54678020 3.60370551 6.26818361 6.54582298 0.28568784 7.33039303 8.28048731 5.61571598 0.62324820 6.55082606 1.48286348 1.52075675 7.99945639 4.22829874 6.37090969 10.71108833 5.65070353 5.91061064 4.79112746 9.56559099 0.46808805 10.03689549 7.07598138 0.09921362 4.03311338 8.24668040 6.91296351 6.09153164 6.42255844 6.62328323 8.72736469 10.35423888 1.05423175 6.35540316 7.98514351 0.75307552 8.54114785 8.87904937 3.24932046 9.03520723 1.31599222 2.77865033 5.75045236 4.56359825 5.29058271 9.52351574 2.26854500 4.92440852 6.16351746 3.69994189 2.48226185 6.22332081 0.10283118 3.33407136 8.61855283 5.64729911 4.29937536 6.80220346 1.52559035 5.30673208 8.42730087 4.33795859 2.44384742 10.92006009 5.97359622 2.09527484 4.54047269 9.36292599 4.59812545 10.15657200 7.11312100 4.12186386 4.10891778 8.14296609 3.03748778 6.46414030 6.76298798 2.82030303 8.36912699 9.96856485 5.13335536 6.40464448 7.88053967 4.83355611 8.56832737 8.83759880 5.01173438 13.86472316 6.74502024 5.14527340 14.10800939 6.08076383 3.51873827 11.40638125 7.33628045 -------------------------------------------------------------------------------------------------------- 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 254501. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3169. 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 1361 Maximum index for augmentation-charges 1468 (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.5983562E+03 (-0.3890998E+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 = 8971.05836900 -Hartree energ DENC = -12252.94656813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.20546154 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00309395 eigenvalues EBANDS = -197.00515932 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 598.35615427 eV energy without entropy = 598.35924822 energy(sigma->0) = 598.35718559 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 2265 total energy-change (2. order) :-0.6874379E+03 (-0.6531641E+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 = 8971.05836900 -Hartree energ DENC = -12252.94656813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.20546154 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.00083274 eigenvalues EBANDS = -884.44536769 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -89.08179289 eV energy without entropy = -89.08096015 energy(sigma->0) = -89.08151531 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2118 total energy-change (2. order) :-0.7324869E+02 (-0.7147841E+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 = 8971.05836900 -Hartree energ DENC = -12252.94656813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.20546154 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02199052 eigenvalues EBANDS = -957.67289506 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -162.33047804 eV energy without entropy = -162.30848752 energy(sigma->0) = -162.32314787 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 2208 total energy-change (2. order) :-0.2986220E+01 (-0.2967862E+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 = 8971.05836900 -Hartree energ DENC = -12252.94656813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.20546154 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02392750 eigenvalues EBANDS = -960.65717843 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.31669840 eV energy without entropy = -165.29277090 energy(sigma->0) = -165.30872256 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 2424 total energy-change (2. order) :-0.1494322E+00 (-0.1493438E+00) number of electron 136.0000018 magnetization 30.2848637 augmentation part -6.9918095 magnetization 27.1456437 Broyden mixing: rms(total) = 0.25334E+01 rms(broyden)= 0.25333E+01 rms(prec ) = 0.26722E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12252.94656813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.20546154 PAW double counting = 3167.08814863 -2605.71786812 entropy T*S EENTRO = -0.02394904 eigenvalues EBANDS = -960.80658910 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.46613060 eV energy without entropy = -165.44218156 energy(sigma->0) = -165.45814759 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 2409 total energy-change (2. order) :-0.7100282E+04 (-0.6715548E+04) number of electron 136.0000031 magnetization 30.1893960 augmentation part -6.6841460 magnetization 26.8906455 Broyden mixing: rms(total) = 0.12139E+02 rms(broyden)= 0.12041E+02 rms(prec ) = 0.12610E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0130 0.0130 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12361.71718651 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.93324442 PAW double counting = 6197.23937244 -5637.90929279 entropy T*S EENTRO = 0.00742628 eigenvalues EBANDS = -7968.58134457 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -7265.74811289 eV energy without entropy = -7265.75553916 energy(sigma->0) = -7265.75058831 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 1983 total energy-change (2. order) : 0.7145852E+04 (-0.7016661E+03) number of electron 136.0000012 magnetization 27.3027060 augmentation part -6.9774002 magnetization 22.2000457 Broyden mixing: rms(total) = 0.25276E+01 rms(broyden)= 0.22002E+01 rms(prec ) = 0.22355E+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.2730 0.5208 0.0252 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12358.59246459 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -237.57657398 PAW double counting = 6338.80800530 -5779.15741550 entropy T*S EENTRO = -0.00297793 eigenvalues EBANDS = -834.52127120 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -119.89654121 eV energy without entropy = -119.89356328 energy(sigma->0) = -119.89554857 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 2310 total energy-change (2. order) :-0.1840305E+02 (-0.9832439E+01) number of electron 136.0000012 magnetization 22.5224820 augmentation part -7.0584493 magnetization 19.7087530 Broyden mixing: rms(total) = 0.15074E+01 rms(broyden)= 0.14870E+01 rms(prec ) = 0.15198E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5851 1.3076 0.0246 0.4232 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12397.32075761 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.20820272 PAW double counting = 9256.84985910 -8697.05498748 entropy T*S EENTRO = -0.01821820 eigenvalues EBANDS = -808.69343631 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -138.29958653 eV energy without entropy = -138.28136833 energy(sigma->0) = -138.29351380 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 2247 total energy-change (2. order) :-0.4127749E+01 (-0.1411305E+01) number of electron 136.0000012 magnetization 18.6480654 augmentation part -6.9937502 magnetization 16.4804779 Broyden mixing: rms(total) = 0.93657E+00 rms(broyden)= 0.93511E+00 rms(prec ) = 0.96462E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7694 1.8153 0.0246 0.7108 0.5268 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12447.37132284 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -240.86096663 PAW double counting = 12597.49344887 -12038.81189337 entropy T*S EENTRO = -0.03320972 eigenvalues EBANDS = -763.98954890 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -142.42733590 eV energy without entropy = -142.39412618 energy(sigma->0) = -142.41626599 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 2208 total energy-change (2. order) :-0.3091702E+01 (-0.2732230E+00) number of electron 136.0000012 magnetization 14.5610432 augmentation part -6.9701113 magnetization 12.7860764 Broyden mixing: rms(total) = 0.78132E+00 rms(broyden)= 0.78121E+00 rms(prec ) = 0.80470E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8260 1.9902 0.0246 0.4947 0.8337 0.7868 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12481.75613616 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -240.76391453 PAW double counting = 13653.05303735 -13095.71382412 entropy T*S EENTRO = -0.02577642 eigenvalues EBANDS = -731.45858114 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -145.51903833 eV energy without entropy = -145.49326191 energy(sigma->0) = -145.51044619 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 1935 total energy-change (2. order) : 0.1449020E+00 (-0.8303093E+02) number of electron 136.0000012 magnetization 13.7142936 augmentation part -6.9710605 magnetization 12.7292765 Broyden mixing: rms(total) = 0.86686E+00 rms(broyden)= 0.86678E+00 rms(prec ) = 0.88114E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7285 2.0134 0.0246 0.8123 0.8123 0.4897 0.2184 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12508.04468128 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -242.87066021 PAW double counting = 13662.80513463 -13106.26505749 entropy T*S EENTRO = -0.00015220 eigenvalues EBANDS = -702.14487643 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -145.37413629 eV energy without entropy = -145.37398409 energy(sigma->0) = -145.37408556 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 2451 total energy-change (2. order) :-0.6206774E+03 (-0.4836807E+03) number of electron 136.0000034 magnetization 13.5918627 augmentation part -6.8258727 magnetization 13.1442760 Broyden mixing: rms(total) = 0.66143E+01 rms(broyden)= 0.65329E+01 rms(prec ) = 0.68169E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6345 2.0184 0.8265 0.8265 0.0246 0.4908 0.2355 0.0194 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12510.12082752 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.54546530 PAW double counting = 13591.18974659 -13034.97485863 entropy T*S EENTRO = -0.00237305 eigenvalues EBANDS = -1319.74387757 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -766.05149879 eV energy without entropy = -766.04912574 energy(sigma->0) = -766.05070778 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 2631 total energy-change (2. order) :-0.2420163E+04 (-0.3291549E+04) number of electron 136.0000027 magnetization 13.5469327 augmentation part -6.9145439 magnetization 11.4174850 Broyden mixing: rms(total) = 0.10226E+02 rms(broyden)= 0.99700E+01 rms(prec ) = 0.10759E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5555 2.0189 0.8264 0.8264 0.4910 0.0246 0.2344 0.0206 0.0019 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12511.01570597 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -239.43733824 PAW double counting = 13659.36443425 -13095.03730128 entropy T*S EENTRO = 0.02512181 eigenvalues EBANDS = -3751.26027709 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -3186.21490983 eV energy without entropy = -3186.24003164 energy(sigma->0) = -3186.22328376 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 2061 total energy-change (2. order) : 0.3046467E+04 (-0.1110521E+04) number of electron 136.0000012 magnetization 9.6587476 augmentation part -6.9400305 magnetization 9.0764722 Broyden mixing: rms(total) = 0.17843E+01 rms(broyden)= 0.62503E+00 rms(prec ) = 0.63657E+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.6779 2.6219 1.2638 0.8721 0.5480 0.5359 0.0246 0.2140 0.0158 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12511.95707863 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -236.22869825 PAW double counting = 13630.98432297 -13073.92907393 entropy T*S EENTRO = -0.00037200 eigenvalues EBANDS = -699.76293940 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -139.74768254 eV energy without entropy = -139.74731054 energy(sigma->0) = -139.74755854 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 2151 total energy-change (2. order) :-0.8302819E+01 (-0.1349305E+02) number of electron 136.0000013 magnetization 9.4507509 augmentation part -6.9591694 magnetization 9.1530283 Broyden mixing: rms(total) = 0.76460E+00 rms(broyden)= 0.69161E+00 rms(prec ) = 0.73002E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6410 2.5767 1.3381 0.8504 0.5702 0.5287 0.0246 0.2508 0.2508 0.0157 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12521.40186771 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.28891148 PAW double counting = 13846.43949767 -13296.70005305 entropy T*S EENTRO = -0.00275995 eigenvalues EBANDS = -682.24256337 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -148.05050121 eV energy without entropy = -148.04774125 energy(sigma->0) = -148.04958122 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 3081 total energy-change (2. order) :-0.8155522E+04 (-0.7511017E+04) number of electron 136.0000063 magnetization 9.2977004 augmentation part -6.7022979 magnetization 12.5154132 Broyden mixing: rms(total) = 0.95647E+01 rms(broyden)= 0.93615E+01 rms(prec ) = 0.10164E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5878 2.6355 1.3521 0.8433 0.5681 0.5285 0.0246 0.2435 0.2435 0.0157 0.0044 0.0061 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12517.48953067 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -246.84799735 PAW double counting = 13739.59035413 -13168.89035139 entropy T*S EENTRO = 0.01079977 eigenvalues EBANDS = -8861.09242946 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -8303.57299828 eV energy without entropy = -8303.58379805 energy(sigma->0) = -8303.57659820 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 1947 total energy-change (2. order) : 0.8159257E+04 (-0.8029973E+03) number of electron 136.0000012 magnetization 7.5392385 augmentation part -6.7698301 magnetization 7.3440118 Broyden mixing: rms(total) = 0.17386E+01 rms(broyden)= 0.75087E+00 rms(prec ) = 0.76920E+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.5800 3.0497 1.1227 0.7516 0.5298 0.5298 0.5238 0.2608 0.0246 0.1381 0.0157 0.0044 0.0090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12518.90527532 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -240.86784366 PAW double counting = 13753.96881151 -13196.68876961 entropy T*S EENTRO = -0.00360729 eigenvalues EBANDS = -692.96586225 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -144.31638993 eV energy without entropy = -144.31278264 energy(sigma->0) = -144.31518750 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 2184 total energy-change (2. order) :-0.6248796E+01 (-0.4885820E+02) number of electron 136.0000015 magnetization 7.2063364 augmentation part -6.9137265 magnetization 7.3069375 Broyden mixing: rms(total) = 0.83466E+00 rms(broyden)= 0.77481E+00 rms(prec ) = 0.79454E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5553 3.0274 1.0630 0.7355 0.7355 0.5366 0.5366 0.2518 0.0246 0.1394 0.1394 0.0157 0.0044 0.0090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12521.74502341 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -247.35774262 PAW double counting = 14773.93963375 -14216.68913918 entropy T*S EENTRO = 0.01952948 eigenvalues EBANDS = -689.87860083 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -150.56518611 eV energy without entropy = -150.58471559 energy(sigma->0) = -150.57169594 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 2283 total energy-change (2. order) :-0.8223028E+02 (-0.3838502E+02) number of electron 136.0000027 magnetization 6.5149040 augmentation part -6.8592572 magnetization 6.5415865 Broyden mixing: rms(total) = 0.52246E+01 rms(broyden)= 0.51697E+01 rms(prec ) = 0.54314E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5491 3.1198 1.2505 0.6963 0.6963 0.7236 0.5066 0.2338 0.0246 0.1568 0.1568 0.0938 0.0157 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12517.18833601 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -249.03289292 PAW double counting = 14780.74842982 -14222.52451213 entropy T*S EENTRO = 0.00183548 eigenvalues EBANDS = -775.94614676 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -232.79546583 eV energy without entropy = -232.79730131 energy(sigma->0) = -232.79607765 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 2331 total energy-change (2. order) : 0.8172568E+02 (-0.2456037E+03) number of electron 136.0000011 magnetization 5.4954869 augmentation part -6.8947549 magnetization 5.5355256 Broyden mixing: rms(total) = 0.99613E+00 rms(broyden)= 0.74466E+00 rms(prec ) = 0.75959E+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.5812 3.3251 1.1159 1.1159 1.0577 0.5562 0.5562 0.2648 0.2648 0.0246 0.1606 0.1491 0.0981 0.0157 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12514.77680466 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -247.77310432 PAW double counting = 14533.98954101 -13974.89718156 entropy T*S EENTRO = 0.01662478 eigenvalues EBANDS = -698.77501494 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -151.06978299 eV energy without entropy = -151.08640778 energy(sigma->0) = -151.07532459 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 2229 total energy-change (2. order) :-0.1466076E+02 (-0.9329467E+01) number of electron 136.0000013 magnetization 3.2359032 augmentation part -6.9111264 magnetization 3.3500022 Broyden mixing: rms(total) = 0.50336E+00 rms(broyden)= 0.48980E+00 rms(prec ) = 0.50970E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6833 4.5592 1.9539 0.8974 0.8974 0.5578 0.5578 0.4801 0.2886 0.2886 0.0246 0.1501 0.1501 0.0982 0.0157 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12507.08723514 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -251.68127085 PAW double counting = 14526.31830922 -13968.45183808 entropy T*S EENTRO = 0.02480859 eigenvalues EBANDS = -715.99947300 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.73054256 eV energy without entropy = -165.75535115 energy(sigma->0) = -165.73881209 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 2118 total energy-change (2. order) :-0.3699963E+01 (-0.8987635E+00) number of electron 136.0000012 magnetization 1.9715349 augmentation part -6.8984948 magnetization 1.9936012 Broyden mixing: rms(total) = 0.40208E+00 rms(broyden)= 0.40105E+00 rms(prec ) = 0.42539E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7509 5.3332 2.3310 0.7852 0.7852 0.7399 0.7399 0.5006 0.5006 0.2989 0.2989 0.0246 0.1502 0.1502 0.0970 0.0157 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12502.89359491 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.21589378 PAW double counting = 13920.43316027 -13361.99864508 entropy T*S EENTRO = 0.02640976 eigenvalues EBANDS = -721.92809808 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.43050513 eV energy without entropy = -169.45691489 energy(sigma->0) = -169.43930838 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 2046 total energy-change (2. order) :-0.3773693E+01 (-0.3236631E+00) number of electron 136.0000011 magnetization 1.5021822 augmentation part -6.8902662 magnetization 1.4585497 Broyden mixing: rms(total) = 0.39300E+00 rms(broyden)= 0.39278E+00 rms(prec ) = 0.42013E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7859 5.4098 2.3691 1.0521 1.0521 0.7511 0.7511 0.6147 0.6147 0.4967 0.2912 0.2912 0.0246 0.1505 0.1505 0.0970 0.0157 0.0044 0.0090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12502.36107645 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.68742816 PAW double counting = 13298.20465407 -12739.76678953 entropy T*S EENTRO = 0.02640522 eigenvalues EBANDS = -724.76611972 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.20419787 eV energy without entropy = -173.23060309 energy(sigma->0) = -173.21299961 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 2055 total energy-change (2. order) : 0.5774024E+00 (-0.2843367E+00) number of electron 136.0000012 magnetization 1.1922675 augmentation part -6.8951301 magnetization 1.1472216 Broyden mixing: rms(total) = 0.30921E+00 rms(broyden)= 0.30916E+00 rms(prec ) = 0.33630E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7985 5.5592 2.5468 1.1636 1.1636 0.7368 0.7368 0.6066 0.6066 0.5083 0.5083 0.2916 0.2916 0.0246 0.1505 0.1505 0.0970 0.0157 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12501.42288654 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.44012949 PAW double counting = 12846.24687594 -12287.83859748 entropy T*S EENTRO = 0.02700943 eigenvalues EBANDS = -724.34522400 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.62679546 eV energy without entropy = -172.65380488 energy(sigma->0) = -172.63579860 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 2424 total energy-change (2. order) :-0.4346009E+02 (-0.2360787E+02) number of electron 136.0000025 magnetization 1.2197963 augmentation part -6.8843696 magnetization 3.2769576 Broyden mixing: rms(total) = 0.49853E+01 rms(broyden)= 0.49026E+01 rms(prec ) = 0.51899E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7637 5.5681 2.5545 1.1679 1.1679 0.7341 0.7341 0.6121 0.6121 0.5436 0.5268 0.2914 0.2914 0.0246 0.1505 0.1505 0.0970 0.0157 0.0189 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12496.79369782 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.89374651 PAW double counting = 12536.08641341 -11977.79661531 entropy T*S EENTRO = -0.02441661 eigenvalues EBANDS = -771.81098342 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -216.08688958 eV energy without entropy = -216.06247296 energy(sigma->0) = -216.07875071 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 1911 total energy-change (2. order) : 0.4997036E+02 (-0.2995289E+03) number of electron 136.0000013 magnetization 0.9603629 augmentation part -6.8882424 magnetization 0.9074656 Broyden mixing: rms(total) = 0.77074E+00 rms(broyden)= 0.27444E+00 rms(prec ) = 0.29530E+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.7783 5.8020 2.6530 1.1087 1.1087 0.9741 0.9741 0.6009 0.6009 0.5354 0.5354 0.3950 0.2904 0.2904 0.0246 0.1504 0.1504 0.0970 0.0157 0.0235 0.0044 0.0090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12496.58668268 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.78718506 PAW double counting = 12510.93197844 -11952.72366844 entropy T*S EENTRO = 0.02721576 eigenvalues EBANDS = -724.12433953 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -166.11652481 eV energy without entropy = -166.14374057 energy(sigma->0) = -166.12559673 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 1905 total energy-change (2. order) :-0.2475286E+01 (-0.2059680E+01) number of electron 136.0000014 magnetization 0.9652900 augmentation part -6.9072007 magnetization 0.9255326 Broyden mixing: rms(total) = 0.45194E+00 rms(broyden)= 0.42939E+00 rms(prec ) = 0.45919E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7543 5.7940 2.7910 1.0467 1.0467 1.0218 1.0218 0.5452 0.5452 0.5644 0.5415 0.2898 0.2898 0.3101 0.3101 0.0246 0.1504 0.1504 0.0970 0.0157 0.0238 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12490.44167339 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.12665714 PAW double counting = 12403.68886320 -11846.64877506 entropy T*S EENTRO = 0.02728515 eigenvalues EBANDS = -729.23700985 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -168.59181039 eV energy without entropy = -168.61909554 energy(sigma->0) = -168.60090544 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 1941 total energy-change (2. order) :-0.1595431E+00 (-0.6539114E+01) number of electron 136.0000012 magnetization 0.8689898 augmentation part -6.9038213 magnetization 0.8348359 Broyden mixing: rms(total) = 0.37698E+00 rms(broyden)= 0.37469E+00 rms(prec ) = 0.39156E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7729 6.1948 2.8996 1.3499 1.0564 1.0564 0.9493 0.5937 0.5937 0.6468 0.4816 0.4816 0.3421 0.2923 0.2923 0.1504 0.1504 0.0246 0.0970 0.0709 0.0157 0.0239 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12489.92573661 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.27439464 PAW double counting = 12391.54051092 -11833.25150367 entropy T*S EENTRO = 0.02727305 eigenvalues EBANDS = -731.01365923 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -168.75135349 eV energy without entropy = -168.77862654 energy(sigma->0) = -168.76044451 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 2019 total energy-change (2. order) :-0.5648267E+01 (-0.2334527E+01) number of electron 136.0000012 magnetization 0.6220954 augmentation part -6.8995220 magnetization 0.5879526 Broyden mixing: rms(total) = 0.29381E+00 rms(broyden)= 0.29364E+00 rms(prec ) = 0.31602E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7836 6.5357 2.7664 1.6804 1.0358 1.0358 0.8208 0.8208 0.5422 0.5422 0.5864 0.5006 0.4028 0.4028 0.2924 0.2924 0.1504 0.1504 0.0246 0.0970 0.0727 0.0157 0.0239 0.0090 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12485.45613498 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.43848635 PAW double counting = 12410.42529856 -11852.14063979 entropy T*S EENTRO = 0.02703555 eigenvalues EBANDS = -740.96285022 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.39962053 eV energy without entropy = -174.42665609 energy(sigma->0) = -174.40863239 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 2190 total energy-change (2. order) :-0.1057123E+02 (-0.3772582E+01) number of electron 136.0000008 magnetization 0.5933305 augmentation part -6.8810693 magnetization 0.5485669 Broyden mixing: rms(total) = 0.81885E+00 rms(broyden)= 0.81733E+00 rms(prec ) = 0.84630E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7527 6.5439 2.7862 1.6745 1.0278 1.0278 0.8191 0.8191 0.5180 0.5180 0.5689 0.5033 0.4220 0.4220 0.2923 0.2923 0.1504 0.1504 0.0246 0.0970 0.0157 0.0239 0.0358 0.0044 0.0090 0.0723 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12484.66762803 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.51927383 PAW double counting = 12411.57696201 -11853.46528329 entropy T*S EENTRO = 0.02697240 eigenvalues EBANDS = -752.06875331 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -184.97084736 eV energy without entropy = -184.99781976 energy(sigma->0) = -184.97983816 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 1836 total energy-change (2. order) :-0.2235499E+03 (-0.2208630E+03) number of electron 136.0000007 magnetization 0.5901566 augmentation part -6.7974613 magnetization 1.6738205 Broyden mixing: rms(total) = 0.62312E+01 rms(broyden)= 0.61143E+01 rms(prec ) = 0.64086E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7242 6.5478 2.7812 1.6724 1.0274 1.0274 0.8289 0.8289 0.5158 0.5158 0.5600 0.5100 0.4215 0.4215 0.2924 0.2924 0.1504 0.1504 0.0970 0.0246 0.0722 0.0373 0.0239 0.0157 0.0090 0.0044 0.0020 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12484.48448849 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.51999828 PAW double counting = 12392.72529870 -11835.57143341 entropy T*S EENTRO = -0.02539918 eigenvalues EBANDS = -974.79090493 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -408.52076889 eV energy without entropy = -408.49536972 energy(sigma->0) = -408.51230250 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 1926 total energy-change (2. order) : 0.2510157E+03 (-0.3761436E+03) number of electron 136.0000014 magnetization 0.6576444 augmentation part -6.8411120 magnetization 0.0013579 Broyden mixing: rms(total) = 0.10899E+01 rms(broyden)= 0.51463E+00 rms(prec ) = 0.52827E+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.7001 6.5313 2.7983 1.6812 1.0247 1.0247 0.8243 0.8243 0.6137 0.5147 0.5147 0.4907 0.4212 0.4212 0.2923 0.2923 0.1504 0.1504 0.0970 0.0246 0.0239 0.0157 0.0702 0.0416 0.0416 0.0090 0.0044 0.0045 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12484.69614265 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.75668229 PAW double counting = 12392.35099086 -11834.73860506 entropy T*S EENTRO = -0.01198633 eigenvalues EBANDS = -727.79880268 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -157.50507144 eV energy without entropy = -157.49308511 energy(sigma->0) = -157.50107600 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 2070 total energy-change (2. order) :-0.1494034E+02 (-0.3741754E+01) number of electron 136.0000012 magnetization 0.6550024 augmentation part -6.8854795 magnetization 0.3273993 Broyden mixing: rms(total) = 0.39036E+00 rms(broyden)= 0.34037E+00 rms(prec ) = 0.35651E+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.6825 6.5187 2.8045 1.7005 1.0205 1.0205 0.8223 0.8223 0.6385 0.4962 0.4962 0.4872 0.4294 0.4294 0.2924 0.2924 0.1504 0.1504 0.1228 0.1228 0.0246 0.0970 0.0741 0.0401 0.0157 0.0239 0.0090 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12482.98357481 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.23237874 PAW double counting = 12414.55115059 -11856.33972444 entropy T*S EENTRO = 0.00388303 eigenvalues EBANDS = -741.59092337 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.44541105 eV energy without entropy = -172.44929408 energy(sigma->0) = -172.44670539 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 2082 total energy-change (2. order) :-0.3211115E+00 (-0.6161285E+00) number of electron 136.0000012 magnetization 0.5341240 augmentation part -6.9011939 magnetization 0.5088733 Broyden mixing: rms(total) = 0.26577E+00 rms(broyden)= 0.26288E+00 rms(prec ) = 0.28111E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7038 6.4943 2.7345 2.0355 1.0136 1.0136 0.8562 0.8562 0.5612 0.5612 0.6253 0.6253 0.5304 0.5304 0.2921 0.2921 0.3053 0.3053 0.1504 0.1504 0.1803 0.0246 0.0970 0.0760 0.0417 0.0157 0.0239 0.0090 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12482.66220966 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.61580131 PAW double counting = 12431.55581520 -11873.41851485 entropy T*S EENTRO = 0.02657914 eigenvalues EBANDS = -741.79854779 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.76652258 eV energy without entropy = -172.79310172 energy(sigma->0) = -172.77538230 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 1911 total energy-change (2. order) : 0.4150902E+01 (-0.4881542E+01) number of electron 136.0000013 magnetization 0.4971577 augmentation part -6.9123535 magnetization 0.4873983 Broyden mixing: rms(total) = 0.35982E+00 rms(broyden)= 0.35973E+00 rms(prec ) = 0.37684E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6961 6.4915 2.7885 2.1048 1.0505 1.0505 0.7822 0.7822 0.6874 0.6874 0.5558 0.5558 0.4783 0.4783 0.4147 0.4147 0.2913 0.2913 0.1504 0.1504 0.1904 0.1904 0.0246 0.0970 0.0157 0.0239 0.0417 0.0759 0.0090 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12480.66959358 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.70806245 PAW double counting = 12485.77262729 -11927.63545649 entropy T*S EENTRO = 0.02669299 eigenvalues EBANDS = -739.54798529 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -168.61562084 eV energy without entropy = -168.64231383 energy(sigma->0) = -168.62451850 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 2571 total energy-change (2. order) :-0.1511321E+05 (-0.1424768E+05) number of electron 136.0000011 magnetization 0.5761418 augmentation part -6.6700644 magnetization -14.1216983 Broyden mixing: rms(total) = 0.76711E+01 rms(broyden)= 0.75538E+01 rms(prec ) = 0.78552E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6731 6.4799 2.7798 2.0959 1.0524 1.0524 0.7587 0.7587 0.7100 0.7100 0.5566 0.5566 0.4796 0.4796 0.4158 0.4158 0.2913 0.2913 0.1504 0.1504 0.1917 0.1917 0.0246 0.0970 0.0759 0.0417 0.0157 0.0239 0.0090 0.0019 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12478.94975483 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.75213482 PAW double counting = 12515.23064383 -11956.59203189 entropy T*S EENTRO = -0.00742685 eigenvalues EBANDS = -15854.89826670 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -15281.82281456 eV energy without entropy = -15281.81538771 energy(sigma->0) = -15281.82033895 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 1896 total energy-change (2. order) : 0.1514177E+05 (-0.4870384E+03) number of electron 136.0000012 magnetization 0.8336509 augmentation part -6.7465441 magnetization -8.5871337 Broyden mixing: rms(total) = 0.16047E+01 rms(broyden)= 0.11908E+01 rms(prec ) = 0.12162E+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.6513 6.4295 2.8202 2.0929 1.0493 1.0493 0.8666 0.7271 0.7271 0.5563 0.5563 0.5442 0.5442 0.4521 0.4233 0.4233 0.2913 0.2913 0.1504 0.1504 0.1910 0.1910 0.0970 0.0246 0.0759 0.0417 0.0157 0.0239 0.0089 0.0089 0.0090 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12479.05203191 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.30519642 PAW double counting = 12531.55826207 -11973.45272941 entropy T*S EENTRO = -0.00083056 eigenvalues EBANDS = -716.94643160 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -140.05280114 eV energy without entropy = -140.05197057 energy(sigma->0) = -140.05252428 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 2160 total energy-change (2. order) :-0.2918211E+02 (-0.1501433E+02) number of electron 136.0000012 magnetization 0.1476643 augmentation part -6.8602591 magnetization -0.2596520 Broyden mixing: rms(total) = 0.40315E+00 rms(broyden)= 0.34172E+00 rms(prec ) = 0.36179E+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.6391 6.5037 2.6600 2.2663 1.0611 1.0611 0.7987 0.7673 0.7673 0.5520 0.5520 0.5985 0.5195 0.4504 0.4504 0.3802 0.2911 0.2911 0.1976 0.1976 0.1504 0.1504 0.0246 0.0970 0.0611 0.0611 0.0759 0.0417 0.0157 0.0239 0.0090 0.0064 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12478.38507940 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.27949376 PAW double counting = 12576.66081002 -12018.64499829 entropy T*S EENTRO = -0.00575173 eigenvalues EBANDS = -742.72655193 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.23490839 eV energy without entropy = -169.22915666 energy(sigma->0) = -169.23299115 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 2460 total energy-change (2. order) :-0.4727652E+03 (-0.3982147E+03) number of electron 136.0000051 magnetization 0.1577505 augmentation part -6.8884136 magnetization 1.9893197 Broyden mixing: rms(total) = 0.77097E+01 rms(broyden)= 0.75673E+01 rms(prec ) = 0.82524E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6207 6.4983 2.6652 2.2704 1.0569 1.0569 0.8190 0.7621 0.7621 0.5588 0.5588 0.5804 0.5205 0.4343 0.4343 0.3887 0.2910 0.2910 0.1960 0.1960 0.1504 0.1504 0.0970 0.0671 0.0671 0.0246 0.0759 0.0157 0.0239 0.0417 0.0237 0.0090 0.0065 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12476.30426896 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.74866051 PAW double counting = 12695.37875178 -12134.80370544 entropy T*S EENTRO = -0.01690840 eigenvalues EBANDS = -1219.65145206 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -642.00008690 eV energy without entropy = -641.98317850 energy(sigma->0) = -641.99445077 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 2940 total energy-change (2. order) :-0.1014564E+06 (-0.4570928E+05) number of electron 136.0000066 magnetization 0.1656805 augmentation part -6.9523653 magnetization -4.4606361 Broyden mixing: rms(total) = 0.12512E+02 rms(broyden)= 0.12400E+02 rms(prec ) = 0.13041E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6035 6.4928 2.7287 2.2372 1.0540 1.0540 0.8370 0.7466 0.7466 0.5997 0.5652 0.5652 0.5092 0.4437 0.4437 0.3600 0.2909 0.2909 0.2035 0.2035 0.1504 0.1504 0.0970 0.0246 0.0624 0.0624 0.0759 0.0417 0.0157 0.0239 0.0192 0.0018 0.0090 0.0065 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12477.09289712 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.46779406 PAW double counting = 12758.46989177 -12132.08419155 entropy T*S EENTRO = -0.01049230 eigenvalues EBANDS = -102744.31146732 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -102098.35079389 eV energy without entropy = -102098.34030159 energy(sigma->0) = -102098.34729646 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 2034 total energy-change (2. order) : 0.8336749E+05 (-0.1876697E+05) number of electron 136.0000054 magnetization 0.1485810 augmentation part -6.9575607 magnetization -0.0006799 Broyden mixing: rms(total) = 0.11679E+02 rms(broyden)= 0.11437E+02 rms(prec ) = 0.12054E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5866 6.4919 2.7332 2.2330 1.0535 1.0535 0.8329 0.7439 0.7439 0.6050 0.5676 0.5676 0.5090 0.4431 0.4431 0.3595 0.2909 0.2909 0.2031 0.2031 0.1504 0.1504 0.0970 0.0621 0.0621 0.0246 0.0759 0.0417 0.0239 0.0157 0.0187 0.0026 0.0002 0.0090 0.0065 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12478.66475280 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -248.97386341 PAW double counting = 12699.73029837 -12141.06071252 entropy T*S EENTRO = -0.01082168 eigenvalues EBANDS = -19313.02372293 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -18730.85741827 eV energy without entropy = -18730.84659659 energy(sigma->0) = -18730.85381104 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 2862 total energy-change (2. order) : 0.1642388E+05 (-0.2036139E+05) number of electron 136.0000089 magnetization 0.3935382 augmentation part -6.9479610 magnetization -4.4451447 Broyden mixing: rms(total) = 0.99430E+01 rms(broyden)= 0.97855E+01 rms(prec ) = 0.10511E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5711 6.4853 2.7407 2.2340 1.0541 1.0541 0.8408 0.7506 0.7506 0.5618 0.5618 0.5867 0.5124 0.4417 0.4417 0.3653 0.2909 0.2909 0.2029 0.2029 0.1504 0.1504 0.0970 0.0628 0.0628 0.0759 0.0246 0.0417 0.0239 0.0157 0.0191 0.0055 0.0055 0.0008 0.0090 0.0065 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12478.41553807 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.63195691 PAW double counting = 12682.85149365 -12117.25341504 entropy T*S EENTRO = -0.00305551 eigenvalues EBANDS = -2899.67602293 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -2306.98233811 eV energy without entropy = -2306.97928260 energy(sigma->0) = -2306.98131960 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 2106 total energy-change (2. order) : 0.2150580E+04 (-0.1115624E+04) number of electron 136.0000016 magnetization 0.6775313 augmentation part -6.8603169 magnetization -0.7459384 Broyden mixing: rms(total) = 0.21920E+01 rms(broyden)= 0.47905E+00 rms(prec ) = 0.49018E+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.5607 6.4821 2.6987 2.2816 1.0548 1.0548 0.8865 0.7562 0.7562 0.5696 0.5696 0.5277 0.5277 0.4238 0.4238 0.3712 0.2909 0.2909 0.1989 0.1989 0.1504 0.1504 0.1255 0.1255 0.0970 0.0759 0.0246 0.0447 0.0416 0.0239 0.0157 0.0161 0.0161 0.0009 0.0107 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12479.13422448 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -249.20480170 PAW double counting = 12694.67202095 -12136.92490492 entropy T*S EENTRO = -0.02471155 eigenvalues EBANDS = -736.93191715 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -156.40238216 eV energy without entropy = -156.37767062 energy(sigma->0) = -156.39414498 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 2109 total energy-change (2. order) :-0.1563311E+02 (-0.6338511E+01) number of electron 136.0000012 magnetization 0.7128558 augmentation part -6.9025523 magnetization 0.3992212 Broyden mixing: rms(total) = 0.53806E+00 rms(broyden)= 0.32971E+00 rms(prec ) = 0.34676E+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.5744 6.4633 2.5611 2.5611 1.0965 1.0965 0.9185 0.6804 0.6310 0.6310 0.6429 0.6429 0.4733 0.4733 0.4873 0.4089 0.4089 0.3321 0.2911 0.2911 0.2662 0.1504 0.1504 0.1588 0.1588 0.0970 0.0246 0.0803 0.0755 0.0417 0.0239 0.0157 0.0144 0.0144 0.0122 0.0009 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12475.29062206 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.07333720 PAW double counting = 12993.86341789 -12436.70493321 entropy T*S EENTRO = 0.00760350 eigenvalues EBANDS = -747.98378035 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -172.03549474 eV energy without entropy = -172.04309824 energy(sigma->0) = -172.03802924 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 2127 total energy-change (2. order) : 0.4354575E+01 (-0.7251735E+01) number of electron 136.0000014 magnetization 0.7054978 augmentation part -6.9214258 magnetization 0.6934509 Broyden mixing: rms(total) = 0.53847E+00 rms(broyden)= 0.53194E+00 rms(prec ) = 0.54454E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5724 6.4727 2.5505 2.5505 1.1290 1.1290 0.8484 0.8484 0.5458 0.5458 0.5956 0.5956 0.5955 0.5955 0.5011 0.4401 0.4401 0.4113 0.2908 0.2908 0.2375 0.2375 0.1504 0.1504 0.1615 0.1615 0.0970 0.0246 0.0781 0.0752 0.0239 0.0157 0.0417 0.0144 0.0144 0.0122 0.0009 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12472.50540637 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.95440484 PAW double counting = 12942.18465397 -12384.59859685 entropy T*S EENTRO = 0.02616063 eigenvalues EBANDS = -745.97948306 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -167.68091984 eV energy without entropy = -167.70708047 energy(sigma->0) = -167.68964005 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 2073 total energy-change (2. order) :-0.6478765E+01 (-0.3893993E+01) number of electron 136.0000012 magnetization 0.6735468 augmentation part -6.9096250 magnetization 0.6467516 Broyden mixing: rms(total) = 0.30357E+00 rms(broyden)= 0.30240E+00 rms(prec ) = 0.32269E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5911 6.3948 2.8121 2.8121 1.2030 1.2030 0.8464 0.8464 0.6442 0.6442 0.6925 0.6925 0.6027 0.6027 0.5146 0.4346 0.4346 0.3715 0.3715 0.2911 0.2911 0.2407 0.2407 0.1504 0.1504 0.1622 0.1622 0.0970 0.0246 0.0785 0.0753 0.0239 0.0157 0.0417 0.0144 0.0144 0.0122 0.0009 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12471.55648762 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.02286582 PAW double counting = 12963.50671966 -12405.24737794 entropy T*S EENTRO = 0.02574088 eigenvalues EBANDS = -754.01157108 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.15968523 eV energy without entropy = -174.18542611 energy(sigma->0) = -174.16826552 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 2040 total energy-change (2. order) : 0.4217250E+01 (-0.7786895E+00) number of electron 136.0000013 magnetization 0.6838320 augmentation part -6.9182648 magnetization 0.6663201 Broyden mixing: rms(total) = 0.21968E+00 rms(broyden)= 0.21957E+00 rms(prec ) = 0.22901E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6010 6.2398 3.0688 2.7923 1.2848 1.2848 1.0206 0.7453 0.7453 0.7395 0.7395 0.5803 0.5803 0.5120 0.5120 0.5337 0.4977 0.4977 0.3755 0.3755 0.2906 0.2906 0.2434 0.2434 0.1504 0.1504 0.1621 0.1621 0.0970 0.0785 0.0753 0.0246 0.0239 0.0157 0.0417 0.0144 0.0144 0.0122 0.0009 0.0063 0.0043 0.0044 0.0090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12471.19969056 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.06156721 PAW double counting = 12853.34463172 -12294.98509078 entropy T*S EENTRO = 0.02615250 eigenvalues EBANDS = -750.21302725 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.94243489 eV energy without entropy = -169.96858739 energy(sigma->0) = -169.95115239 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 2037 total energy-change (2. order) :-0.7788840E+00 (-0.1178801E+01) number of electron 136.0000013 magnetization 0.6667051 augmentation part -6.9167733 magnetization 0.6484634 Broyden mixing: rms(total) = 0.43719E+00 rms(broyden)= 0.43707E+00 rms(prec ) = 0.45048E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6003 6.2082 3.3113 2.7683 1.2358 1.2358 0.9633 0.7911 0.7911 0.8561 0.7186 0.7186 0.5671 0.5671 0.5078 0.5078 0.4995 0.4995 0.3641 0.3641 0.2909 0.2909 0.2791 0.2143 0.2143 0.1504 0.1504 0.1628 0.1628 0.0970 0.0246 0.0785 0.0753 0.0417 0.0239 0.0157 0.0144 0.0144 0.0122 0.0009 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12470.99327475 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.06466366 PAW double counting = 12798.35662640 -12240.01014781 entropy T*S EENTRO = 0.02636499 eigenvalues EBANDS = -751.18238076 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -170.72131889 eV energy without entropy = -170.74768388 energy(sigma->0) = -170.73010722 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 1986 total energy-change (2. order) :-0.3242096E+01 (-0.1216032E+01) number of electron 136.0000012 magnetization 0.6473090 augmentation part -6.9099008 magnetization 0.6238279 Broyden mixing: rms(total) = 0.28710E+00 rms(broyden)= 0.28691E+00 rms(prec ) = 0.30599E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5999 6.1817 3.4293 2.7147 1.2415 1.2415 1.0876 0.8049 0.8049 0.7949 0.7662 0.7662 0.5480 0.5480 0.5073 0.5073 0.4708 0.4708 0.4179 0.4179 0.3292 0.2910 0.2910 0.2877 0.2144 0.2144 0.1504 0.1504 0.1624 0.1624 0.0970 0.0246 0.0785 0.0753 0.0417 0.0239 0.0157 0.0144 0.0144 0.0122 0.0009 0.0090 0.0063 0.0043 0.0044 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12470.89780377 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.06945846 PAW double counting = 12812.32400812 -12253.95159110 entropy T*S EENTRO = 0.02642820 eigenvalues EBANDS = -754.54115485 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -173.96341517 eV energy without entropy = -173.98984337 energy(sigma->0) = -173.97222457 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 2082 total energy-change (2. order) :-0.4829140E+00 (-0.1980280E-01) number of electron 136.0000012 magnetization 0.3786337 augmentation part -6.9086482 magnetization 0.3547408 Broyden mixing: rms(total) = 0.28770E+00 rms(broyden)= 0.28767E+00 rms(prec ) = 0.30751E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5608 6.0540 3.3713 1.8212 1.0006 1.0006 0.8460 0.8460 0.6583 0.6583 0.6722 0.5376 0.5376 0.6044 0.4095 0.4095 0.3695 0.2759 0.3045 0.3045 0.2716 0.2716 0.2053 0.2053 0.1748 0.1748 0.0985 0.0985 0.0576 0.0576 0.0242 0.0380 0.0100 0.0140 0.0140 0.0142 0.0083 0.0007 0.0018 0.0044 0.0036 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12470.94862588 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.08131249 PAW double counting = 12762.94078107 -12204.65672312 entropy T*S EENTRO = 0.02654862 eigenvalues EBANDS = -754.87315407 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -174.44632918 eV energy without entropy = -174.47287780 energy(sigma->0) = -174.45517872 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 2583 total energy-change (2. order) :-0.1597333E+05 (-0.1488043E+05) number of electron 136.0033119 magnetization 0.3724748 augmentation part -6.6981262 magnetization -10.8426415 Broyden mixing: rms(total) = 0.81676E+01 rms(broyden)= 0.80570E+01 rms(prec ) = 0.83598E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5477 6.0671 3.3646 1.8323 1.0011 1.0011 0.8539 0.8539 0.6577 0.6577 0.6304 0.6304 0.5340 0.5340 0.4097 0.4097 0.3709 0.3047 0.3047 0.2813 0.2698 0.2698 0.2057 0.2057 0.1744 0.1744 0.0972 0.0972 0.0576 0.0576 0.0238 0.0394 0.0101 0.0199 0.0130 0.0130 0.0092 0.0008 0.0021 0.0049 0.0049 0.0042 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 8971.05836900 -Hartree energ DENC = -12468.82422981 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.11959354 PAW double counting = 12613.79270921 -12056.12490405 entropy T*S EENTRO = -0.00796510 eigenvalues EBANDS = -16729.63575405 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -16147.77358066 eV energy without entropy = -16147.76561555 energy(sigma->0) = -16147.77092562 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 8 1 9 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------