vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.09.02 02:11:32 running on 1 total cores distrk: each k-point on 1 cores, 1 groups distr: one band on NCORE= 1 cores, 1 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 200 ALGO = Normal (blocked Davidson) ISPIN = 2 MAGMOM = 36*0 0.107 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.77 0.32 POTCAR: PAW_PBE Si 05Jan2001 POTCAR: PAW_PBE C 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 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 C 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 Optimization of the real space projectors (new method) maximal supplied QI-value = 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 = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.30 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.053 115.676 0.49E-03 0.72E-03 0.18E-06 0 8 10.053 87.132 0.49E-03 0.71E-03 0.18E-06 1 7 10.053 4.429 0.32E-03 0.31E-03 0.18E-06 1 7 10.053 2.733 0.23E-03 0.19E-03 0.20E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 PAW_PBE Si 05Jan2001 : energy of atom 1 EATOM= -103.0669 kinetic energy error for atom= 0.0012 (will be added to EATOM!!) PAW_PBE C 08Apr2002 : energy of atom 2 EATOM= -147.1560 kinetic energy error for atom= 0.0288 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 3 EATOM= -12.4884 kinetic energy error for atom= 0.0098 (will be added to EATOM!!) 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.899 0.490 0.111- 3 1.71 2 0.744 0.138 0.369- 3 0.004 0.543 0.212- 1 1.71 17 2.61 4 0.980 0.161 0.184- 5 0.914 0.372 0.332- 36 1.19 6 1.50 16 2.41 15 2.61 6 0.936 0.441 0.391- 36 0.50 5 1.50 16 1.82 15 1.85 7 0.376 0.406 0.304- 23 0.98 20 2.23 29 2.55 8 0.546 0.521 0.635- 24 1.64 32 1.72 28 2.28 9 0.375 0.243 0.885- 26 1.27 12 2.18 21 2.50 10 0.965 0.774 0.831- 11 0.661 0.684 0.429- 12 0.319 0.257 0.689- 9 2.18 26 2.59 13 0.618 0.287 0.502- 25 1.50 20 2.05 14 0.772 0.424 0.833- 28 2.19 15 0.065 0.420 0.530- 16 0.32 29 1.74 6 1.85 5 2.61 16 0.076 0.408 0.512- 15 0.32 36 1.63 29 1.77 6 1.82 5 2.41 17 0.287 0.494 0.111- 19 2.19 3 2.61 18 0.129 0.219 0.479- 19 0.550 0.536 0.120- 17 2.19 20 0.394 0.308 0.406- 25 0.98 33 1.97 13 2.05 7 2.23 23 2.24 21 0.487 0.329 0.036- 35 1.48 9 2.50 22 0.401 0.019 0.132- 27 1.75 23 0.493 0.409 0.343- 7 0.98 20 2.24 24 0.349 0.513 0.578- 32 1.10 29 1.47 8 1.64 31 2.31 25 0.464 0.267 0.426- 20 0.98 13 1.50 26 0.502 0.203 0.865- 9 1.27 12 2.59 27 0.292 0.992 0.266- 22 1.75 28 0.521 0.420 0.738- 14 2.19 8 2.28 29 0.258 0.463 0.499- 24 1.47 15 1.74 16 1.77 32 2.14 7 2.55 30 0.188 0.389 0.914- 31 0.306 0.623 0.643- 32 2.07 24 2.31 32 0.411 0.553 0.522- 24 1.10 8 1.72 31 2.07 29 2.14 33 0.399 0.362 0.559- 20 1.97 34 0.912 0.274 0.403- 35 0.361 0.384 0.063- 21 1.48 36 0.973 0.421 0.382- 6 0.50 5 1.19 16 1.63 37 0.811 0.978 0.727- ----------------------------------------------------------------------------- | | | 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 !!! | | | | The distance between some ions is very small. Please check the | | nearest-neighbor list in the OUTCAR file. | | I HOPE YOU KNOW WHAT YOU ARE DOING! | | | ----------------------------------------------------------------------------- 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.898697230 0.489521300 0.110994430 0.743825880 0.137923800 0.369380460 0.003663460 0.542561330 0.211609630 0.979620230 0.161406380 0.184416190 0.913532600 0.372407380 0.332381750 0.935699720 0.441272890 0.390839800 0.376362510 0.405943840 0.304453550 0.545700130 0.520525190 0.634892900 0.375280060 0.243430400 0.884603590 0.964567820 0.773652730 0.831153880 0.661088710 0.683881780 0.429111720 0.319010570 0.256518360 0.688933070 0.618406960 0.287242310 0.502073670 0.771992410 0.424489670 0.833474740 0.064801160 0.420067910 0.529929340 0.075952870 0.407947970 0.512115080 0.286505870 0.493665580 0.110747460 0.128552040 0.218783440 0.478963760 0.550218870 0.536307020 0.120374540 0.394399280 0.307677310 0.405744460 0.486512330 0.328852320 0.035724680 0.400527680 0.019327000 0.132331480 0.492568700 0.409301660 0.342797550 0.348760590 0.512731740 0.577618470 0.463567060 0.267457180 0.425751490 0.501578520 0.202785350 0.864714680 0.292432110 0.992456550 0.266284710 0.520859110 0.419789230 0.738038730 0.258051270 0.463087170 0.498510460 0.187734170 0.389100240 0.913574840 0.306057260 0.623298130 0.643266730 0.410925310 0.552997480 0.522165820 0.399093260 0.361519870 0.558639330 0.911815570 0.273795140 0.402910120 0.361356720 0.383889500 0.062767960 0.972950840 0.421103100 0.382247100 0.811173840 0.978487820 0.727249030 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 36 36 37 37 ---------------------------------------------------------------------------------------- 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= 87 number of dos NEDOS = 301 number of ions NIONS = 37 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 207360 max r-space proj IRMAX = 1222 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 1 4 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 = 644.9 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 12.01 1.00 Ionic Valenz ZVAL = 4.00 4.00 1.00 Atomic Wigner-Seitz radii RWIGS = 1.11 0.77 0.32 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.29E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 44.13 297.83 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 19 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.89869723 0.48952130 0.11099443 0.74382588 0.13792380 0.36938046 0.00366346 0.54256133 0.21160963 0.97962023 0.16140638 0.18441619 0.91353260 0.37240738 0.33238175 0.93569972 0.44127289 0.39083980 0.37636251 0.40594384 0.30445355 0.54570013 0.52052519 0.63489290 0.37528006 0.24343040 0.88460359 0.96456782 0.77365273 0.83115388 0.66108871 0.68388178 0.42911172 0.31901057 0.25651836 0.68893307 0.61840696 0.28724231 0.50207367 0.77199241 0.42448967 0.83347474 0.06480116 0.42006791 0.52992934 0.07595287 0.40794797 0.51211508 0.28650587 0.49366558 0.11074746 0.12855204 0.21878344 0.47896376 0.55021887 0.53630702 0.12037454 0.39439928 0.30767731 0.40574446 0.48651233 0.32885232 0.03572468 0.40052768 0.01932700 0.13233148 0.49256870 0.40930166 0.34279755 0.34876059 0.51273174 0.57761847 0.46356706 0.26745718 0.42575149 0.50157852 0.20278535 0.86471468 0.29243211 0.99245655 0.26628471 0.52085911 0.41978923 0.73803873 0.25805127 0.46308717 0.49851046 0.18773417 0.38910024 0.91357484 0.30605726 0.62329813 0.64326673 0.41092531 0.55299748 0.52216582 0.39909326 0.36151987 0.55863933 0.91181557 0.27379514 0.40291012 0.36135672 0.38388950 0.06276796 0.97295084 0.42110310 0.38224710 0.81117384 0.97848782 0.72724903 position of ions in cartesian coordinates (Angst): 6.88680674 9.62550627 1.20287550 5.70001210 2.71200947 4.00307208 0.02807346 10.66843769 2.29326858 7.50692778 3.17374979 1.99856620 7.00049167 7.32268355 3.60210744 7.17036052 8.67679296 4.23563253 2.88410355 7.98211432 3.29944228 4.18175467 10.23513886 6.88049943 2.87580863 4.78659630 9.58667910 7.39157966 15.21241100 9.00743070 5.06598889 13.44723583 4.65039528 2.44460990 5.04394616 7.46614680 4.73891438 5.64807427 5.44110290 5.91585504 8.34678283 9.03258246 0.49657777 8.25983732 5.74298204 0.58203444 8.02152173 5.54992427 2.19552313 9.70699567 1.20019902 0.98510714 4.30196066 5.19065480 4.21638222 10.54545856 1.30453019 3.02232112 6.04988971 4.39715821 3.72819264 6.46625605 0.38715765 3.06928366 0.38002873 1.43411065 3.77460320 8.04813947 3.71498618 2.67258728 10.08189548 6.25980154 3.55236074 5.25903728 4.61397959 3.84364636 3.98738862 9.37113781 2.24093650 19.51477239 2.88579664 3.99139545 8.25435761 7.99831761 1.97747269 9.10572933 5.40248747 1.43862572 7.65091693 9.90064807 2.34534739 12.25597346 6.97124880 3.14896174 10.87364475 5.65884675 3.05829156 7.10860136 6.05411967 6.98733389 5.38366122 4.36644173 2.76911268 7.54845763 0.68023270 7.45581958 8.28019237 4.14251121 6.21610625 19.24010385 7.88138682 -------------------------------------------------------------------------------------------------------- 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 590071. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3145. kBytes fftplans : 57504. kBytes grid : 164904. kBytes one-center: 227. kBytes wavefun : 334291. 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 0.1070000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 1135 Maximum index for augmentation-charges 4133 (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 : 1392 total energy-change (2. order) : 0.1551480E+04 (-0.3923999E+04) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7622.48113964 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -290.08167409 PAW double counting = 2937.01645563 -2374.63634374 entropy T*S EENTRO = 0.02374079 eigenvalues EBANDS = -214.46943000 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1551.48019517 eV energy without entropy = 1551.45645437 energy(sigma->0) = 1551.47228157 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 1952 total energy-change (2. order) :-0.7856059E+03 (-0.7516931E+03) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7622.48113964 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -290.08167409 PAW double counting = 2937.01645563 -2374.63634374 entropy T*S EENTRO = -0.00822114 eigenvalues EBANDS = -1000.04334131 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 765.87432192 eV energy without entropy = 765.88254306 energy(sigma->0) = 765.87706230 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2656 total energy-change (2. order) :-0.1902366E+03 (-0.1843931E+03) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7622.48113964 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -290.08167409 PAW double counting = 2937.01645563 -2374.63634374 entropy T*S EENTRO = 0.01787891 eigenvalues EBANDS = -1190.30599724 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 575.63776605 eV energy without entropy = 575.61988714 energy(sigma->0) = 575.63180641 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 2072 total energy-change (2. order) :-0.2032266E+02 (-0.1976858E+02) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7622.48113964 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -290.08167409 PAW double counting = 2937.01645563 -2374.63634374 entropy T*S EENTRO = 0.01884322 eigenvalues EBANDS = -1210.62962165 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 555.31510595 eV energy without entropy = 555.29626273 energy(sigma->0) = 555.30882488 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 2112 total energy-change (2. order) :-0.1093700E+01 (-0.1078279E+01) number of electron 136.0000007 magnetization 0.1044053 augmentation part -9.7362304 magnetization 0.0158369 Broyden mixing: rms(total) = 0.12439E+03 rms(broyden)= 0.12439E+03 rms(prec ) = 0.12448E+03 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7622.48113964 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -290.08167409 PAW double counting = 2937.01645563 -2374.63634374 entropy T*S EENTRO = 0.01972998 eigenvalues EBANDS = -1211.72420817 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 554.22140618 eV energy without entropy = 554.20167620 energy(sigma->0) = 554.21482952 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 2302 total energy-change (2. order) : 0.5480788E+02 (-0.1144442E+03) number of electron 136.0000001 magnetization 0.1130529 augmentation part -8.4601637 magnetization -0.0573169 Broyden mixing: rms(total) = 0.23018E+02 rms(broyden)= 0.23018E+02 rms(prec ) = 0.23448E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8559 0.8559 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8143.29366108 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.77457032 PAW double counting = 1782792.55391401 -1782237.48983356 entropy T*S EENTRO = 0.00369272 eigenvalues EBANDS = -639.07884423 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 609.02928375 eV energy without entropy = 609.02559103 energy(sigma->0) = 609.02805284 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 2129 total energy-change (2. order) : 0.4460636E+02 (-0.3563018E+02) number of electron 136.0000008 magnetization 0.1398717 augmentation part -8.7316774 magnetization -0.2235660 Broyden mixing: rms(total) = 0.34389E+02 rms(broyden)= 0.34389E+02 rms(prec ) = 0.34598E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8287 1.1892 0.4681 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7800.32739384 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.08274289 PAW double counting = 1225230.38118695 -1224676.43193640 entropy T*S EENTRO = -0.04879066 eigenvalues EBANDS = -936.96326801 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 653.63564137 eV energy without entropy = 653.68443203 energy(sigma->0) = 653.65190492 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 2152 total energy-change (2. order) : 0.2022255E+02 (-0.2582321E+02) number of electron 136.0000002 magnetization 0.1414456 augmentation part -8.5431903 magnetization -0.1234065 Broyden mixing: rms(total) = 0.22743E+02 rms(broyden)= 0.22743E+02 rms(prec ) = 0.22840E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5849 1.2205 0.4693 0.0650 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7976.15731928 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.92829342 PAW double counting = 1678518.42793696 -1677965.99479018 entropy T*S EENTRO = 0.00919229 eigenvalues EBANDS = -741.60712095 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 673.85819163 eV energy without entropy = 673.84899935 energy(sigma->0) = 673.85512754 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 2056 total energy-change (2. order) :-0.6191290E+01 (-0.7827317E+01) number of electron 136.0000002 magnetization 0.0723349 augmentation part -8.7807324 magnetization 0.0476839 Broyden mixing: rms(total) = 0.18110E+02 rms(broyden)= 0.18110E+02 rms(prec ) = 0.18254E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6038 1.2254 0.5872 0.3013 0.3013 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7970.92064998 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.77223556 PAW double counting = 1702473.56603938 -1701921.21640687 entropy T*S EENTRO = 0.00334081 eigenvalues EBANDS = -753.10177204 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 667.66690194 eV energy without entropy = 667.66356113 energy(sigma->0) = 667.66578834 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 1784 total energy-change (2. order) : 0.3271588E+01 (-0.1670798E+01) number of electron 136.0000003 magnetization 0.0055461 augmentation part -8.7085165 magnetization 0.2013321 Broyden mixing: rms(total) = 0.22319E+02 rms(broyden)= 0.22319E+02 rms(prec ) = 0.22400E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6364 1.2412 0.5632 0.5632 0.5158 0.2984 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7978.93863281 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -275.81472682 PAW double counting = 1556525.69453239 -1555974.11320902 entropy T*S EENTRO = 0.05250857 eigenvalues EBANDS = -742.05056833 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 670.93849018 eV energy without entropy = 670.88598161 energy(sigma->0) = 670.92098732 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 1720 total energy-change (2. order) :-0.9073063E+01 (-0.1867663E+01) number of electron 136.0000001 magnetization -0.1501630 augmentation part -8.6320355 magnetization 0.4163614 Broyden mixing: rms(total) = 0.26675E+02 rms(broyden)= 0.26675E+02 rms(prec ) = 0.26758E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5697 1.2341 0.6433 0.6433 0.5111 0.2221 0.1643 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8010.06743689 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -275.38987366 PAW double counting = 1497888.14883121 -1497336.94991808 entropy T*S EENTRO = -0.03195276 eigenvalues EBANDS = -719.95280911 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 661.86542692 eV energy without entropy = 661.89737968 energy(sigma->0) = 661.87607784 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 2208 total energy-change (2. order) :-0.4125017E+01 (-0.6633657E+00) number of electron 136.0000002 magnetization -1.0569838 augmentation part -8.5823542 magnetization -0.1255947 Broyden mixing: rms(total) = 0.30146E+02 rms(broyden)= 0.30146E+02 rms(prec ) = 0.30237E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8874 1.4872 1.4872 1.4723 0.6553 0.3660 0.3660 0.3775 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8019.35855098 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -274.76445993 PAW double counting = 1398480.11005240 -1397929.27668634 entropy T*S EENTRO = 0.03095936 eigenvalues EBANDS = -715.10949074 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 657.74040998 eV energy without entropy = 657.70945062 energy(sigma->0) = 657.73009019 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 1801 total energy-change (2. order) :-0.1921249E+02 (-0.1651223E+02) number of electron 136.0000002 magnetization -1.4723867 augmentation part -8.4051138 magnetization 0.0545427 Broyden mixing: rms(total) = 0.41986E+02 rms(broyden)= 0.41986E+02 rms(prec ) = 0.42150E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8981 1.6690 1.6690 1.6452 0.6428 0.3744 0.3744 0.4049 0.4049 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8032.81704624 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -265.23304300 PAW double counting = 1215391.62441135 -1214843.01400246 entropy T*S EENTRO = 0.02109970 eigenvalues EBANDS = -728.16208527 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 638.52792028 eV energy without entropy = 638.50682058 energy(sigma->0) = 638.52088705 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 1985 total energy-change (2. order) : 0.1175568E+01 (-0.1097350E+02) number of electron 136.0000004 magnetization -1.5087813 augmentation part -8.4119196 magnetization 0.1067494 Broyden mixing: rms(total) = 0.51893E+02 rms(broyden)= 0.51893E+02 rms(prec ) = 0.51978E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8167 1.7351 1.5744 1.5744 0.5380 0.4358 0.4358 0.4166 0.3203 0.3203 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8016.83734771 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -261.17593961 PAW double counting = 1236713.09642242 -1236165.26995429 entropy T*S EENTRO = -0.00510146 eigenvalues EBANDS = -746.21317723 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 639.70348833 eV energy without entropy = 639.70858979 energy(sigma->0) = 639.70518882 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 2008 total energy-change (2. order) :-0.5025389E+01 (-0.2832409E+01) number of electron 136.0000005 magnetization -1.3730464 augmentation part -8.4375241 magnetization 0.2430860 Broyden mixing: rms(total) = 0.55997E+02 rms(broyden)= 0.55997E+02 rms(prec ) = 0.56106E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7748 1.7449 1.5312 1.5312 0.5201 0.5201 0.4203 0.4203 0.3824 0.3824 0.2953 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8012.75739204 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -261.03021332 PAW double counting = 1336952.43173842 -1336404.74946363 entropy T*S EENTRO = -0.00428397 eigenvalues EBANDS = -755.32087210 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 634.67809957 eV energy without entropy = 634.68238354 energy(sigma->0) = 634.67952756 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 1872 total energy-change (2. order) :-0.3556597E+00 (-0.1151015E+01) number of electron 136.0000007 magnetization -1.2012781 augmentation part -8.4993480 magnetization 0.1241388 Broyden mixing: rms(total) = 0.52520E+02 rms(broyden)= 0.52520E+02 rms(prec ) = 0.52636E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7595 1.7412 1.4688 1.4688 0.6309 0.6309 0.5319 0.5319 0.3734 0.3734 0.3016 0.3016 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7983.05908042 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -263.98259259 PAW double counting = 1490036.63568443 -1489488.71437600 entropy T*S EENTRO = 0.00388278 eigenvalues EBANDS = -782.66966452 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 634.32243989 eV energy without entropy = 634.31855711 energy(sigma->0) = 634.32114563 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 1831 total energy-change (2. order) : 0.1267298E+02 (-0.1064250E+01) number of electron 136.0000006 magnetization -1.1792717 augmentation part -8.4947971 magnetization -0.4055284 Broyden mixing: rms(total) = 0.40531E+02 rms(broyden)= 0.40531E+02 rms(prec ) = 0.40637E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7345 1.7628 1.4671 1.4671 0.7573 0.7573 0.5161 0.5161 0.3767 0.3767 0.3155 0.3155 0.1855 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7992.36631970 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -265.99153525 PAW double counting = 1737653.29915595 -1737105.01828925 entropy T*S EENTRO = -0.00119991 eigenvalues EBANDS = -759.03497483 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 646.99542322 eV energy without entropy = 646.99662313 energy(sigma->0) = 646.99582319 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 2024 total energy-change (2. order) : 0.7928683E+01 (-0.4080341E+00) number of electron 136.0000006 magnetization -1.7111756 augmentation part -8.5348267 magnetization -0.9276461 Broyden mixing: rms(total) = 0.33401E+02 rms(broyden)= 0.33401E+02 rms(prec ) = 0.33513E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7564 2.1939 1.1922 1.1922 0.7874 0.7874 0.6433 0.6433 0.5877 0.3621 0.3621 0.3878 0.3878 0.3061 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7996.86301438 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -263.92674068 PAW double counting = 1851381.26307011 -1850832.80393886 entropy T*S EENTRO = 0.00043222 eigenvalues EBANDS = -748.85428829 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 654.92410633 eV energy without entropy = 654.92367412 energy(sigma->0) = 654.92396226 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 1494 total energy-change (2. order) : 0.1755059E+02 (-0.3384492E+01) number of electron 136.0000005 magnetization -2.0071469 augmentation part -8.5344371 magnetization 0.5822890 Broyden mixing: rms(total) = 0.33840E+02 rms(broyden)= 0.33840E+02 rms(prec ) = 0.33929E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7853 2.7732 1.0640 1.0640 1.0707 1.0707 0.6012 0.6012 0.5669 0.3674 0.3674 0.4111 0.4111 0.3366 0.2887 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8042.69437492 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -247.75570165 PAW double counting = 1919034.08773508 -1918486.38253591 entropy T*S EENTRO = 0.01730098 eigenvalues EBANDS = -700.90631490 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 672.47469490 eV energy without entropy = 672.45739392 energy(sigma->0) = 672.46892790 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 1880 total energy-change (2. order) :-0.3591048E+01 (-0.1337143E+01) number of electron 136.0000006 magnetization -1.9946001 augmentation part -8.4228031 magnetization 0.8058395 Broyden mixing: rms(total) = 0.36031E+02 rms(broyden)= 0.36031E+02 rms(prec ) = 0.36157E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7763 2.8918 1.1619 1.1619 1.1870 0.9936 0.5717 0.5717 0.5608 0.3997 0.3997 0.3498 0.3498 0.3685 0.3685 0.3083 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8041.20034409 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.43926058 PAW double counting = 1976946.04548993 -1976398.70203243 entropy T*S EENTRO = 0.02605479 eigenvalues EBANDS = -707.95484703 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 668.88364680 eV energy without entropy = 668.85759202 energy(sigma->0) = 668.87496188 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 1936 total energy-change (2. order) : 0.7385735E+01 (-0.5661413E+00) number of electron 136.0000005 magnetization -2.0113972 augmentation part -8.5165225 magnetization 0.5428773 Broyden mixing: rms(total) = 0.30837E+02 rms(broyden)= 0.30837E+02 rms(prec ) = 0.30924E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7847 2.4473 1.3019 1.3019 1.5436 1.0595 0.6110 0.6110 0.5108 0.5108 0.3613 0.3613 0.5079 0.4043 0.4043 0.3094 0.3094 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8041.72095029 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -247.37513242 PAW double counting = 2106423.36232575 -2105875.75622289 entropy T*S EENTRO = 0.00884396 eigenvalues EBANDS = -698.35806869 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 676.26938162 eV energy without entropy = 676.26053766 energy(sigma->0) = 676.26643363 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 1912 total energy-change (2. order) :-0.3851336E+00 (-0.4603300E+00) number of electron 136.0000005 magnetization -2.0607897 augmentation part -8.5917937 magnetization 0.3945584 Broyden mixing: rms(total) = 0.23934E+02 rms(broyden)= 0.23934E+02 rms(prec ) = 0.24022E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7533 2.5303 1.5247 1.3221 1.3221 0.9579 0.6170 0.6170 0.5348 0.5348 0.5243 0.3635 0.3635 0.3819 0.3819 0.3039 0.2634 0.2634 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8021.92664397 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.78679799 PAW double counting = 2202127.15542489 -2201579.22196786 entropy T*S EENTRO = -0.01618975 eigenvalues EBANDS = -711.42816352 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 675.88424802 eV energy without entropy = 675.90043777 energy(sigma->0) = 675.88964460 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 2040 total energy-change (2. order) : 0.8299801E+00 (-0.1684688E+00) number of electron 136.0000005 magnetization -2.0364499 augmentation part -8.6009534 magnetization 0.4362132 Broyden mixing: rms(total) = 0.21288E+02 rms(broyden)= 0.21288E+02 rms(prec ) = 0.21371E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7192 2.3889 1.5594 1.3216 1.3216 0.9849 0.6226 0.6226 0.5327 0.5327 0.5254 0.3623 0.3623 0.3874 0.3874 0.3022 0.3022 0.2145 0.2145 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8019.23158474 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.88555604 PAW double counting = 2210579.33782871 -2210031.33334417 entropy T*S EENTRO = -0.00463866 eigenvalues EBANDS = -712.27706316 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 676.71422815 eV energy without entropy = 676.71886680 energy(sigma->0) = 676.71577437 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 1872 total energy-change (2. order) :-0.6637792E+00 (-0.2207303E-01) number of electron 136.0000005 magnetization -2.0016826 augmentation part -8.6104202 magnetization 0.4496542 Broyden mixing: rms(total) = 0.21211E+02 rms(broyden)= 0.21211E+02 rms(prec ) = 0.21294E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7165 2.4336 1.3323 1.3323 1.4758 1.0031 0.6161 0.6161 0.3835 0.5461 0.5461 0.5279 0.3673 0.3673 0.3929 0.3929 0.3354 0.3354 0.3050 0.3050 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8015.16077125 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.00782493 PAW double counting = 2193104.02338525 -2192555.93219027 entropy T*S EENTRO = -0.00564399 eigenvalues EBANDS = -715.97509208 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 676.05044893 eV energy without entropy = 676.05609292 energy(sigma->0) = 676.05233026 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 2216 total energy-change (2. order) :-0.3200907E+00 (-0.1282946E-01) number of electron 136.0000005 magnetization -2.1455984 augmentation part -8.6134609 magnetization 0.3622042 Broyden mixing: rms(total) = 0.20394E+02 rms(broyden)= 0.20394E+02 rms(prec ) = 0.20479E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7478 2.3045 1.7981 1.3309 1.3309 0.7397 0.7397 0.9098 0.6048 0.6048 0.5887 0.5887 0.4991 0.3635 0.3635 0.3916 0.3916 0.3957 0.3957 0.3069 0.3069 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8015.97253203 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.41905695 PAW double counting = 2218038.83846101 -2217490.75404555 entropy T*S EENTRO = -0.00659145 eigenvalues EBANDS = -715.06446300 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 675.73035822 eV energy without entropy = 675.73694968 energy(sigma->0) = 675.73255538 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 1648 total energy-change (2. order) : 0.4413612E+01 (-0.2034319E+00) number of electron 136.0000004 magnetization -2.1995346 augmentation part -8.6235115 magnetization 0.7870727 Broyden mixing: rms(total) = 0.19744E+02 rms(broyden)= 0.19744E+02 rms(prec ) = 0.19811E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7397 2.3676 1.7805 1.3092 1.3092 0.9283 0.9283 0.8979 0.5993 0.5993 0.6010 0.6010 0.4978 0.3619 0.3619 0.3882 0.3882 0.3909 0.3909 0.3063 0.3063 0.2207 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8005.67655095 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.76156146 PAW double counting = 2175590.11507014 -2175041.83375353 entropy T*S EENTRO = -0.02895765 eigenvalues EBANDS = -722.77886223 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 680.14397053 eV energy without entropy = 680.17292818 energy(sigma->0) = 680.15362308 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 1928 total energy-change (2. order) : 0.3603176E+00 (-0.3944848E-01) number of electron 136.0000005 magnetization -2.3608222 augmentation part -8.6290079 magnetization 0.4955998 Broyden mixing: rms(total) = 0.19832E+02 rms(broyden)= 0.19832E+02 rms(prec ) = 0.19904E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7752 2.2663 1.8593 1.4661 1.4661 1.2076 1.2076 0.8758 0.5948 0.5948 0.6673 0.6673 0.5111 0.3971 0.3971 0.3621 0.3621 0.3973 0.3973 0.3696 0.3696 0.3090 0.3090 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -8003.00698301 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.16358758 PAW double counting = 2152748.68534019 -2152200.39404597 entropy T*S EENTRO = -0.02246886 eigenvalues EBANDS = -725.70255289 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 680.50428809 eV energy without entropy = 680.52675695 energy(sigma->0) = 680.51177771 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 1544 total energy-change (2. order) : 0.1175833E+01 (-0.2084569E+00) number of electron 136.0000004 magnetization -2.4522317 augmentation part -8.6420029 magnetization 0.7347336 Broyden mixing: rms(total) = 0.21082E+02 rms(broyden)= 0.21082E+02 rms(prec ) = 0.21145E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7907 2.2921 1.7777 1.7777 1.7768 1.2374 1.2374 0.9296 0.5913 0.5913 0.6971 0.6971 0.4451 0.4451 0.3620 0.3620 0.4902 0.4102 0.4102 0.3945 0.3945 0.3051 0.3051 0.2560 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7993.89930937 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -254.80311305 PAW double counting = 2049572.09800982 -2049023.68739292 entropy T*S EENTRO = -0.00553844 eigenvalues EBANDS = -734.13112075 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 681.68012149 eV energy without entropy = 681.68565994 energy(sigma->0) = 681.68196764 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 1848 total energy-change (2. order) : 0.7849797E+00 (-0.8289566E-01) number of electron 136.0000004 magnetization -2.5868242 augmentation part -8.6527668 magnetization 0.5248694 Broyden mixing: rms(total) = 0.21668E+02 rms(broyden)= 0.21668E+02 rms(prec ) = 0.21728E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8266 2.1486 2.1486 1.9853 1.9853 1.2996 1.2996 1.0217 0.5886 0.5886 0.7300 0.7300 0.5084 0.5084 0.4786 0.4786 0.5119 0.3612 0.3612 0.3904 0.3904 0.3139 0.3139 0.3478 0.3478 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7989.26720156 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.09539236 PAW double counting = 2002934.19423555 -2002385.71242817 entropy T*S EENTRO = -0.00987250 eigenvalues EBANDS = -737.75282593 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 682.46510124 eV energy without entropy = 682.47497373 energy(sigma->0) = 682.46839207 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 1624 total energy-change (2. order) : 0.3575637E+01 (-0.4439068E+00) number of electron 136.0000004 magnetization -2.6396892 augmentation part -8.6686345 magnetization 1.0842546 Broyden mixing: rms(total) = 0.20837E+02 rms(broyden)= 0.20837E+02 rms(prec ) = 0.20885E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8118 2.2388 2.2388 1.9754 1.9754 1.2886 1.2886 1.0201 0.5883 0.5883 0.6621 0.6621 0.6277 0.6277 0.4205 0.4205 0.3612 0.3612 0.4947 0.3904 0.3904 0.3938 0.3938 0.3059 0.3059 0.2755 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7982.16575838 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.41520886 PAW double counting = 1959975.88450903 -1959427.16670521 entropy T*S EENTRO = 0.00542014 eigenvalues EBANDS = -741.21010496 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.04073797 eV energy without entropy = 686.03531783 energy(sigma->0) = 686.03893126 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 1883 total energy-change (2. order) : 0.2982298E+00 (-0.6121294E-01) number of electron 136.0000003 magnetization -2.6515460 augmentation part -8.6735742 magnetization 1.0610477 Broyden mixing: rms(total) = 0.20246E+02 rms(broyden)= 0.20246E+02 rms(prec ) = 0.20294E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7956 2.2947 2.2947 1.9793 1.9793 1.2807 1.2807 1.0037 0.5882 0.5882 0.6273 0.6273 0.6661 0.6661 0.4111 0.4111 0.4917 0.3619 0.3619 0.3968 0.3968 0.3890 0.3890 0.3092 0.3092 0.2912 0.2912 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7981.26001587 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.31999225 PAW double counting = 1945075.58355551 -1944526.84782167 entropy T*S EENTRO = 0.00184111 eigenvalues EBANDS = -741.92718524 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.33896780 eV energy without entropy = 686.33712669 energy(sigma->0) = 686.33835410 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 1880 total energy-change (2. order) :-0.3391912E+00 (-0.7443675E-02) number of electron 136.0000003 magnetization -2.7012432 augmentation part -8.6741602 magnetization 1.0284095 Broyden mixing: rms(total) = 0.20330E+02 rms(broyden)= 0.20330E+02 rms(prec ) = 0.20377E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8180 2.4683 2.4683 1.9938 1.9938 1.2603 1.2603 1.0066 0.7527 0.7527 0.5920 0.5920 0.6931 0.6931 0.4826 0.4826 0.4597 0.4597 0.3613 0.3613 0.4945 0.3840 0.3840 0.3862 0.3862 0.3074 0.3074 0.3011 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7981.35272655 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.64712164 PAW double counting = 1940163.97713146 -1939615.23697826 entropy T*S EENTRO = 0.00118769 eigenvalues EBANDS = -741.85030227 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 685.99977664 eV energy without entropy = 685.99858895 energy(sigma->0) = 685.99938074 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 1633 total energy-change (2. order) : 0.6153036E-01 (-0.2536454E-01) number of electron 136.0000003 magnetization -3.0292107 augmentation part -8.6755624 magnetization 0.7513618 Broyden mixing: rms(total) = 0.20418E+02 rms(broyden)= 0.20418E+02 rms(prec ) = 0.20465E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8917 3.1601 3.1601 1.7931 1.7931 1.3045 1.3045 1.0967 1.0967 1.2148 0.5901 0.5901 0.7592 0.7592 0.5907 0.5907 0.4579 0.4579 0.3614 0.3614 0.5122 0.4525 0.4525 0.3952 0.3952 0.3742 0.3082 0.3176 0.3176 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7981.33068596 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -255.95956655 PAW double counting = 1930834.11328927 -1930285.34951414 entropy T*S EENTRO = 0.00392508 eigenvalues EBANDS = -741.52472691 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.06130700 eV energy without entropy = 686.05738192 energy(sigma->0) = 686.05999864 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 1532 total energy-change (2. order) : 0.3776616E+01 (-0.6706336E+00) number of electron 136.0000003 magnetization -3.0335240 augmentation part -8.6975841 magnetization 0.3612170 Broyden mixing: rms(total) = 0.18463E+02 rms(broyden)= 0.18463E+02 rms(prec ) = 0.18501E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8821 3.3728 3.3728 1.6448 1.6448 1.3500 1.3500 1.3458 1.1092 1.1092 0.5900 0.5900 0.7369 0.7369 0.6263 0.6263 0.4561 0.4561 0.3614 0.3614 0.5063 0.4492 0.4492 0.3937 0.3937 0.3515 0.3099 0.3230 0.3230 0.2414 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7977.00873643 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -256.35094622 PAW double counting = 1941200.07314051 -1940651.06001626 entropy T*S EENTRO = -0.01688048 eigenvalues EBANDS = -741.90722448 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 689.83792285 eV energy without entropy = 689.85480333 energy(sigma->0) = 689.84354968 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 2024 total energy-change (2. order) :-0.4959189E+00 (-0.1401086E+00) number of electron 136.0000003 magnetization -3.0202472 augmentation part -8.7142877 magnetization 0.1952083 Broyden mixing: rms(total) = 0.16645E+02 rms(broyden)= 0.16645E+02 rms(prec ) = 0.16683E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8640 3.4006 3.4006 1.6486 1.6486 1.3599 1.3599 1.3323 1.1128 1.1128 0.5900 0.5900 0.7339 0.7339 0.6292 0.6292 0.4558 0.4558 0.3614 0.3614 0.4970 0.4495 0.4495 0.3957 0.3957 0.3673 0.3083 0.3177 0.3177 0.2524 0.2524 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7977.22673103 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.44948138 PAW double counting = 1957968.38090693 -1957419.32637425 entropy T*S EENTRO = -0.00716795 eigenvalues EBANDS = -741.13773458 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 689.34200396 eV energy without entropy = 689.34917191 energy(sigma->0) = 689.34439328 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 2168 total energy-change (2. order) :-0.1556133E+00 (-0.7553642E-02) number of electron 136.0000003 magnetization -3.0193025 augmentation part -8.7160328 magnetization 0.1722492 Broyden mixing: rms(total) = 0.16484E+02 rms(broyden)= 0.16484E+02 rms(prec ) = 0.16521E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8376 3.4071 3.4071 1.6432 1.6432 1.3596 1.3596 1.3379 1.1137 1.1137 0.5900 0.5900 0.7327 0.7327 0.6309 0.6309 0.4559 0.4559 0.3614 0.3614 0.4970 0.4502 0.4502 0.3957 0.3957 0.3666 0.3083 0.3179 0.3179 0.0678 0.2359 0.2359 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7977.45781876 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.65161733 PAW double counting = 1960965.94802733 -1960416.89180769 entropy T*S EENTRO = -0.00482856 eigenvalues EBANDS = -740.86415052 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 689.18639067 eV energy without entropy = 689.19121923 energy(sigma->0) = 689.18800019 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 1816 total energy-change (2. order) :-0.5602061E-03 (-0.5334672E-02) number of electron 136.0000003 magnetization -2.9596522 augmentation part -8.7182170 magnetization 0.1853591 Broyden mixing: rms(total) = 0.16014E+02 rms(broyden)= 0.16014E+02 rms(prec ) = 0.16051E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8720 3.4725 3.4725 1.7781 1.7781 1.3387 1.3387 1.0838 1.0838 1.1906 0.8169 0.8169 0.5902 0.5902 0.7178 0.7178 0.5990 0.5990 0.4568 0.4568 0.3614 0.3614 0.4983 0.4400 0.4400 0.4001 0.4001 0.3933 0.3933 0.3772 0.3078 0.3163 0.3163 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7977.77886169 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -257.94442054 PAW double counting = 1976537.86998496 -1975988.80182253 entropy T*S EENTRO = -0.00100117 eigenvalues EBANDS = -740.26663478 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 689.18583046 eV energy without entropy = 689.18683163 energy(sigma->0) = 689.18616418 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 1554 total energy-change (2. order) :-0.6545188E+00 (-0.1396298E-01) number of electron 136.0000003 magnetization -2.7171685 augmentation part -8.7174811 magnetization 0.4987070 Broyden mixing: rms(total) = 0.15643E+02 rms(broyden)= 0.15643E+02 rms(prec ) = 0.15681E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8989 3.6149 3.6149 1.8475 1.8475 1.5860 1.5860 1.2484 1.2484 1.0522 0.8810 0.8810 0.5902 0.5902 0.7331 0.7331 0.5753 0.5753 0.4540 0.4540 0.4618 0.4618 0.5297 0.3614 0.3614 0.4404 0.4404 0.3959 0.3959 0.3818 0.3818 0.3075 0.3158 0.3158 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7979.46523056 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -258.52242930 PAW double counting = 1996710.37044535 -1996161.33841885 entropy T*S EENTRO = -0.01088766 eigenvalues EBANDS = -738.61075358 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 688.53131162 eV energy without entropy = 688.54219928 energy(sigma->0) = 688.53494084 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 1664 total energy-change (2. order) :-0.1328233E+01 (-0.1818247E+00) number of electron 136.0000003 magnetization -2.6609682 augmentation part -8.7311740 magnetization 0.3536001 Broyden mixing: rms(total) = 0.13025E+02 rms(broyden)= 0.13025E+02 rms(prec ) = 0.13062E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8725 3.5706 3.5706 1.8599 1.8599 1.5292 1.5292 1.2490 1.2490 0.8658 0.8658 1.0315 0.5902 0.5902 0.7263 0.7263 0.2063 0.5732 0.5732 0.4561 0.4561 0.4724 0.4724 0.5365 0.3614 0.3614 0.4368 0.4368 0.3954 0.3954 0.3894 0.3894 0.3076 0.3158 0.3158 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7979.58359485 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -261.26629784 PAW double counting = 2098118.92866632 -2097569.78465415 entropy T*S EENTRO = -0.01224283 eigenvalues EBANDS = -737.18738420 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 687.20307865 eV energy without entropy = 687.21532148 energy(sigma->0) = 687.20715959 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 2064 total energy-change (2. order) :-0.2590225E-01 (-0.4724359E-01) number of electron 136.0000003 magnetization -2.6128888 augmentation part -8.7479594 magnetization 0.3576487 Broyden mixing: rms(total) = 0.11594E+02 rms(broyden)= 0.11594E+02 rms(prec ) = 0.11634E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8664 3.5977 3.5977 1.8680 1.8680 1.5899 1.5899 1.2637 1.2637 0.8687 0.8687 1.0077 0.5902 0.5902 0.7212 0.7212 0.3366 0.3366 0.5735 0.5735 0.4560 0.4560 0.5395 0.4745 0.4745 0.3614 0.3614 0.4314 0.4314 0.3931 0.3931 0.3921 0.3921 0.3076 0.3159 0.3159 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7979.92945187 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -261.65209894 PAW double counting = 2117484.08709237 -2116934.92075661 entropy T*S EENTRO = -0.01523007 eigenvalues EBANDS = -736.50096470 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 687.17717640 eV energy without entropy = 687.19240647 energy(sigma->0) = 687.18225309 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 1888 total energy-change (2. order) :-0.4172065E+00 (-0.6701849E-02) number of electron 136.0000003 magnetization -2.6260077 augmentation part -8.7478931 magnetization 0.3882242 Broyden mixing: rms(total) = 0.11509E+02 rms(broyden)= 0.11509E+02 rms(prec ) = 0.11550E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8680 3.5795 3.5795 1.8754 1.8754 1.5674 1.5674 1.2631 1.2631 0.9084 0.9084 0.9895 0.7352 0.7352 0.5903 0.5903 0.6580 0.6580 0.5800 0.5800 0.4570 0.4570 0.5378 0.4877 0.4877 0.3614 0.3614 0.4248 0.4248 0.3939 0.3939 0.3919 0.3919 0.3076 0.3158 0.3158 0.2351 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7980.59878916 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -262.03718357 PAW double counting = 2125809.23073365 -2125260.07324501 entropy T*S EENTRO = -0.02274734 eigenvalues EBANDS = -735.84738490 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.75996989 eV energy without entropy = 686.78271722 energy(sigma->0) = 686.76755233 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 1744 total energy-change (2. order) : 0.2773423E+00 (-0.2810171E-02) number of electron 136.0000003 magnetization -2.5206263 augmentation part -8.7497377 magnetization 0.4801794 Broyden mixing: rms(total) = 0.11487E+02 rms(broyden)= 0.11487E+02 rms(prec ) = 0.11527E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9424 3.9992 3.9992 2.1381 2.1381 1.9986 1.9986 1.3136 1.3136 1.2016 1.2016 0.9256 0.7853 0.7853 0.5902 0.5902 0.6515 0.6515 0.5990 0.5990 0.4587 0.4587 0.5565 0.5136 0.5136 0.3614 0.3614 0.4626 0.4626 0.3921 0.3921 0.3951 0.3951 0.3076 0.3160 0.3160 0.3638 0.3638 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.95083145 Ewald energy TEWEN = 5596.01212093 -Hartree energ DENC = -7980.78018056 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -261.89866933 PAW double counting = 2124981.46914431 -2124432.30416395 entropy T*S EENTRO = -0.02276010 eigenvalues EBANDS = -735.53464435 atomic energy EATOM = 3495.14563384 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 687.03731223 eV energy without entropy = 687.06007233 energy(sigma->0) = 687.04489893 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 7 1 8 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------