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.17 11:12:47 running on 2 total cores distrk: each k-point on 2 cores, 1 groups distr: one band on NCORE= 1 cores, 2 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 95 ALGO = Normal (blocked Davidson) ISPIN = 2 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 0.75 0.32 0.77 1.11 0.73 NPAR = 2 POTCAR: PAW_PBE N 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE C 08Apr2002 POTCAR: PAW_PBE Si 05Jan2001 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 N 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 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 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 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 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 = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.65 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 10 10.053 79.467 0.76E-04 0.72E-04 0.56E-06 0 10 10.053 66.151 0.76E-04 0.72E-04 0.55E-06 1 10 10.053 8.350 0.25E-03 0.92E-03 0.41E-05 1 10 10.053 5.531 0.27E-03 0.10E-02 0.45E-05 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.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 = 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 = 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 N 08Apr2002 : energy of atom 1 EATOM= -264.5486 kinetic energy error for atom= 0.0736 (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 C 08Apr2002 : energy of atom 3 EATOM= -147.1560 kinetic energy error for atom= 0.0288 (will be added to EATOM!!) PAW_PBE Si 05Jan2001 : energy of atom 4 EATOM= -103.0669 kinetic energy error for atom= 0.0012 (will be added to EATOM!!) PAW_PBE O 08Apr2002 : energy of atom 5 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.812 0.391 0.590- 3 1.03 2 1.04 2 0.829 0.421 0.648- 1 1.04 3 0.851 0.402 0.546- 1 1.03 4 0.143 0.520 0.364- 25 1.10 5 0.195 0.617 0.372- 25 1.10 6 0.214 0.532 0.293- 25 1.10 7 0.101 0.557 0.657- 26 1.10 8 0.143 0.645 0.600- 26 1.10 9 0.097 0.562 0.538- 26 1.10 10 0.412 0.438 0.294- 27 1.10 11 0.409 0.347 0.371- 27 1.10 12 0.337 0.417 0.356- 27 1.10 13 0.516 0.273 0.573- 28 1.10 14 0.524 0.350 0.483- 28 1.10 15 0.543 0.385 0.594- 28 1.12 16 0.239 0.696 0.505- 29 1.10 17 0.295 0.671 0.594- 29 1.10 18 0.322 0.661 0.481- 29 1.10 19 0.433 0.621 0.559- 30 1.10 20 0.441 0.558 0.659- 30 1.10 21 0.504 0.553 0.576- 30 1.10 22 0.315 0.451 0.641- 24 1.10 23 0.299 0.382 0.548- 24 1.10 24 0.316 0.450 0.568- 23 1.10 22 1.10 32 1.86 31 1.88 25 0.195 0.544 0.361- 6 1.10 4 1.10 5 1.10 33 1.42 26 0.130 0.573 0.597- 7 1.10 8 1.10 9 1.10 34 1.42 27 0.392 0.416 0.358- 10 1.10 12 1.10 11 1.10 35 1.43 28 0.509 0.343 0.553- 13 1.10 14 1.10 15 1.12 36 1.42 29 0.279 0.650 0.526- 18 1.10 16 1.10 17 1.10 32 1.88 30 0.450 0.557 0.587- 19 1.10 20 1.10 21 1.10 31 1.87 31 0.406 0.460 0.534- 36 1.66 35 1.66 30 1.87 24 1.88 32 0.253 0.530 0.526- 33 1.66 34 1.66 24 1.86 29 1.88 33 0.235 0.498 0.423- 25 1.42 32 1.66 34 0.187 0.518 0.594- 26 1.42 32 1.66 35 0.416 0.476 0.425- 27 1.43 31 1.66 36 0.441 0.365 0.566- 28 1.42 31 1.66 LATTYP: Found a simple tetragonal cell. ALAT = 15.0000000000 C/A-ratio = 1.3333333333 Lattice vectors: A1 = ( 0.0000000000, 15.0000000000, 0.0000000000) A2 = ( 0.0000000000, 0.0000000000, 15.0000000000) A3 = ( 20.0000000000, 0.0000000000, 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 tetragonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 16 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 tetragonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 16 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 tetragonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 16 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 : 4500.0000 direct lattice vectors reciprocal lattice vectors 20.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 length of vectors 20.000000000 15.000000000 15.000000000 0.050000000 0.066666667 0.066666667 position of ions in fractional coordinates (direct lattice) 0.812103130 0.391084680 0.589759290 0.828671710 0.421094070 0.647859800 0.851359510 0.402071890 0.546228850 0.142576110 0.520074830 0.363647310 0.194596440 0.616977670 0.372451530 0.214042430 0.532225890 0.293136500 0.101056740 0.556524170 0.656904670 0.142948470 0.644791270 0.599854510 0.097405060 0.562428030 0.538072860 0.412057300 0.438384900 0.293781590 0.409346630 0.347445790 0.370915550 0.337066810 0.416825470 0.355662110 0.516377870 0.273210010 0.573237390 0.524270250 0.350406240 0.482536020 0.542572200 0.385345350 0.594144380 0.239057160 0.695703560 0.505082030 0.294802500 0.671175140 0.593540250 0.321786910 0.660517460 0.480913350 0.432953970 0.620566400 0.558691160 0.441159280 0.558249370 0.659276820 0.504389370 0.552545710 0.576009030 0.315004330 0.451193460 0.641350510 0.299171560 0.382451800 0.548465630 0.316375810 0.449650660 0.567790780 0.194604980 0.544314330 0.360712800 0.129819480 0.573343730 0.596685840 0.392196320 0.416280230 0.358492800 0.508679640 0.342921900 0.552663380 0.279257570 0.650275490 0.526301090 0.450073310 0.556898080 0.586895040 0.406314670 0.460353090 0.534170890 0.252730650 0.530255080 0.526370470 0.235349340 0.497596160 0.422897160 0.187448980 0.517878840 0.594066380 0.416023160 0.476347760 0.425278490 0.440606380 0.364918410 0.566063370 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 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 1 1 1 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.050000000 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 1.000000000 Length of vectors 0.050000000 0.066666667 0.066666667 Shift w.r.t. Gamma in fractional coordinates (k-lattice) 0.000000000 0.000000000 0.000000000 Subroutine IBZKPT returns following result: =========================================== Found 1 irreducible k-points: Following reciprocal coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 1 k-points in BZ NKDIM = 1 number of bands NBANDS= 80 number of dos NEDOS = 301 number of ions NIONS = 36 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 627200 max r-space proj IRMAX = 2665 max aug-charges IRDMAX= 4859 dimension x,y,z NGX = 98 NGY = 80 NGZ = 80 dimension x,y,z NGXF= 196 NGYF= 160 NGZF= 160 support grid NGXF= 196 NGYF= 160 NGZF= 160 ions per type = 1 22 7 2 4 NGX,Y,Z is equivalent to a cutoff of 8.15, 8.87, 8.87 a.u. NGXF,Y,Z is equivalent to a cutoff of 16.29, 17.73, 17.73 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. 32.61 24.46 24.46*2*pi/ulx,y,z ENINI = 400.0 initial cutoff ENAUG = 644.9 eV augmentation charge cutoff NELM = 95; NELMIN= 2; NELMDL= -5 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM LREAL = T real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 1 Vosko Wilk Nusair interpolation ROPT = -0.00050 -0.00050 -0.00050 -0.00050 ROPT = -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.914E-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 = 14.00 1.00 12.01 28.09 16.00 Ionic Valenz ZVAL = 5.00 1.00 4.00 4.00 6.00 Atomic Wigner-Seitz radii RWIGS = 0.75 0.32 0.77 1.11 0.73 virtual crystal weights VCA = 1.00 1.00 1.00 1.00 1.00 NELECT = 87.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.31E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 125.00 843.54 Fermi-wavevector in a.u.,A,eV,Ry = 0.439383 0.830314 2.626712 0.193058 Thomas-Fermi vector in A = 1.413435 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 36 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 : 4500.00 direct lattice vectors reciprocal lattice vectors 20.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 length of vectors 20.000000000 15.000000000 15.000000000 0.050000000 0.066666667 0.066666667 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 1.000 k-points in reciprocal lattice and weights: Automatic mesh 0.00000000 0.00000000 0.00000000 1.000 position of ions in fractional coordinates (direct lattice) 0.81210313 0.39108468 0.58975929 0.82867171 0.42109407 0.64785980 0.85135951 0.40207189 0.54622885 0.14257611 0.52007483 0.36364731 0.19459644 0.61697767 0.37245153 0.21404243 0.53222589 0.29313650 0.10105674 0.55652417 0.65690467 0.14294847 0.64479127 0.59985451 0.09740506 0.56242803 0.53807286 0.41205730 0.43838490 0.29378159 0.40934663 0.34744579 0.37091555 0.33706681 0.41682547 0.35566211 0.51637787 0.27321001 0.57323739 0.52427025 0.35040624 0.48253602 0.54257220 0.38534535 0.59414438 0.23905716 0.69570356 0.50508203 0.29480250 0.67117514 0.59354025 0.32178691 0.66051746 0.48091335 0.43295397 0.62056640 0.55869116 0.44115928 0.55824937 0.65927682 0.50438937 0.55254571 0.57600903 0.31500433 0.45119346 0.64135051 0.29917156 0.38245180 0.54846563 0.31637581 0.44965066 0.56779078 0.19460498 0.54431433 0.36071280 0.12981948 0.57334373 0.59668584 0.39219632 0.41628023 0.35849280 0.50867964 0.34292190 0.55266338 0.27925757 0.65027549 0.52630109 0.45007331 0.55689808 0.58689504 0.40631467 0.46035309 0.53417089 0.25273065 0.53025508 0.52637047 0.23534934 0.49759616 0.42289716 0.18744898 0.51787884 0.59406638 0.41602316 0.47634776 0.42527849 0.44060638 0.36491841 0.56606337 position of ions in cartesian coordinates (Angst): 16.24206260 5.86627020 8.84638935 16.57343420 6.31641105 9.71789700 17.02719020 6.03107835 8.19343275 2.85152220 7.80112245 5.45470965 3.89192880 9.25466505 5.58677295 4.28084860 7.98338835 4.39704750 2.02113480 8.34786255 9.85357005 2.85896940 9.67186905 8.99781765 1.94810120 8.43642045 8.07109290 8.24114600 6.57577350 4.40672385 8.18693260 5.21168685 5.56373325 6.74133620 6.25238205 5.33493165 10.32755740 4.09815015 8.59856085 10.48540500 5.25609360 7.23804030 10.85144400 5.78018025 8.91216570 4.78114320 10.43555340 7.57623045 5.89605000 10.06762710 8.90310375 6.43573820 9.90776190 7.21370025 8.65907940 9.30849600 8.38036740 8.82318560 8.37374055 9.88915230 10.08778740 8.28818565 8.64013545 6.30008660 6.76790190 9.62025765 5.98343120 5.73677700 8.22698445 6.32751620 6.74475990 8.51686170 3.89209960 8.16471495 5.41069200 2.59638960 8.60015595 8.95028760 7.84392640 6.24420345 5.37739200 10.17359280 5.14382850 8.28995070 5.58515140 9.75413235 7.89451635 9.00146620 8.35347120 8.80342560 8.12629340 6.90529635 8.01256335 5.05461300 7.95382620 7.89555705 4.70698680 7.46394240 6.34345740 3.74897960 7.76818260 8.91099570 8.32046320 7.14521640 6.37917735 8.81212760 5.47377615 8.49095055 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 81909 maximum and minimum number of plane-waves per node : 81909 81909 maximum number of plane-waves: 81909 maximum index in each direction: IXMAX= 32 IYMAX= 24 IZMAX= 24 IXMIN= -32 IYMIN= -24 IZMIN= -24 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 140 to avoid them WARNING: aliasing errors must be expected set NGY to 98 to avoid them WARNING: aliasing errors must be expected set NGZ to 98 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 541232. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 4076. kBytes fftplans : 87706. kBytes grid : 314101. kBytes one-center: 221. kBytes wavefun : 105128. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 65 NGY = 49 NGZ = 49 (NGX =196 NGY =160 NGZ =160) gives a total of 156065 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 87.0000000 magnetization 36.0000000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 2540 Maximum index for augmentation-charges 2257 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.107 Maximum number of real-space cells 2x 3x 3 Maximum number of reciprocal cells 3x 3x 3 ----------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) : 0.6218386E+03 (-0.2139969E+04) number of electron 87.0000000 magnetization 36.0000000 augmentation part 87.0000000 magnetization 36.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8128.13327426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.08260407 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = -0.00209303 eigenvalues EBANDS = -611.04773734 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 621.83860152 eV energy without entropy = 621.84069455 energy(sigma->0) = 621.83929920 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 536 total energy-change (2. order) :-0.5409294E+03 (-0.4858286E+03) number of electron 87.0000000 magnetization 36.0000000 augmentation part 87.0000000 magnetization 36.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8128.13327426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.08260407 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.01524961 eigenvalues EBANDS = -1151.99447853 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 80.90920297 eV energy without entropy = 80.89395336 energy(sigma->0) = 80.90411977 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 432 total energy-change (2. order) :-0.2286351E+03 (-0.2269427E+03) number of electron 87.0000000 magnetization 36.0000000 augmentation part 87.0000000 magnetization 36.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8128.13327426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.08260407 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = -0.00244250 eigenvalues EBANDS = -1380.61186352 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -147.72587413 eV energy without entropy = -147.72343163 energy(sigma->0) = -147.72505996 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 432 total energy-change (2. order) :-0.2155583E+02 (-0.2139941E+02) number of electron 87.0000000 magnetization 36.0000000 augmentation part 87.0000000 magnetization 36.0000000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8128.13327426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.08260407 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.01677114 eigenvalues EBANDS = -1402.18691065 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.28170761 eV energy without entropy = -169.29847875 energy(sigma->0) = -169.28729799 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.4671432E+00 (-0.4658968E+00) number of electron 86.9999976 magnetization 30.9695109 augmentation part 4.2675890 magnetization 30.1118629 Broyden mixing: rms(total) = 0.41844E+01 rms(broyden)= 0.41819E+01 rms(prec ) = 0.43574E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8128.13327426 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.08260407 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.01724222 eigenvalues EBANDS = -1402.65452489 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.74885077 eV energy without entropy = -169.76609299 energy(sigma->0) = -169.75459818 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 368 total energy-change (2. order) : 0.9842731E+02 (-0.2405903E+02) number of electron 86.9999981 magnetization 26.4600822 augmentation part 3.8668578 magnetization 25.1560879 Broyden mixing: rms(total) = 0.20308E+01 rms(broyden)= 0.20297E+01 rms(prec ) = 0.21149E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9133 0.9133 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8330.80271071 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 360.65466345 PAW double counting = 4241.97855853 -4221.50874832 entropy T*S EENTRO = 0.01726111 eigenvalues EBANDS = -1195.86930539 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -71.32154003 eV energy without entropy = -71.33880114 energy(sigma->0) = -71.32729374 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 416 total energy-change (2. order) :-0.5019087E+02 (-0.4026228E+01) number of electron 86.9999983 magnetization 22.4639704 augmentation part 3.5910981 magnetization 21.2296125 Broyden mixing: rms(total) = 0.12850E+01 rms(broyden)= 0.12847E+01 rms(prec ) = 0.13234E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9143 1.0465 0.7821 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8435.63974144 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 324.73079954 PAW double counting = 5756.37312804 -5736.64814688 entropy T*S EENTRO = 0.01705428 eigenvalues EBANDS = -1104.55424146 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -121.51240663 eV energy without entropy = -121.52946091 energy(sigma->0) = -121.51809139 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 392 total energy-change (2. order) :-0.1492576E+02 (-0.6408839E+00) number of electron 86.9999983 magnetization 16.9184598 augmentation part 3.5750819 magnetization 15.7592463 Broyden mixing: rms(total) = 0.82367E+00 rms(broyden)= 0.82360E+00 rms(prec ) = 0.84745E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1620 1.7124 1.1174 0.6562 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8475.23637703 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 315.87523097 PAW double counting = 6667.97942121 -6648.13270216 entropy T*S EENTRO = 0.01910985 eigenvalues EBANDS = -1071.15159273 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -136.43816859 eV energy without entropy = -136.45727844 energy(sigma->0) = -136.44453854 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) :-0.2862621E+02 (-0.1233255E+01) number of electron 86.9999983 magnetization 13.7583681 augmentation part 3.5897925 magnetization 12.6607253 Broyden mixing: rms(total) = 0.40822E+00 rms(broyden)= 0.40805E+00 rms(prec ) = 0.41829E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2400 2.2078 1.3096 0.7633 0.6794 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8511.74745977 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 295.65352680 PAW double counting = 7493.50495048 -7473.49031064 entropy T*S EENTRO = 0.02626140 eigenvalues EBANDS = -1043.22009271 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -165.06438315 eV energy without entropy = -165.09064455 energy(sigma->0) = -165.07313695 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 408 total energy-change (2. order) :-0.1160499E+02 (-0.5634871E+00) number of electron 86.9999983 magnetization 10.4381662 augmentation part 3.5620237 magnetization 9.3590599 Broyden mixing: rms(total) = 0.31059E+00 rms(broyden)= 0.31047E+00 rms(prec ) = 0.31916E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2120 1.8941 1.8941 0.8274 0.8274 0.6169 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8526.21735340 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 288.05661218 PAW double counting = 7652.22430727 -7632.16247008 entropy T*S EENTRO = 0.03328220 eigenvalues EBANDS = -1032.81248886 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -176.66936941 eV energy without entropy = -176.70265161 energy(sigma->0) = -176.68046348 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.6846845E+01 (-0.3289738E+00) number of electron 86.9999983 magnetization 6.5548507 augmentation part 3.5381327 magnetization 5.5260167 Broyden mixing: rms(total) = 0.21270E+00 rms(broyden)= 0.21261E+00 rms(prec ) = 0.21847E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4133 2.5623 2.5623 1.1490 0.8875 0.7103 0.6082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8527.96024886 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 282.80586411 PAW double counting = 7560.81864373 -7540.68437891 entropy T*S EENTRO = 0.02985468 eigenvalues EBANDS = -1032.73469085 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -183.51621481 eV energy without entropy = -183.54606949 energy(sigma->0) = -183.52616637 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.4694766E+01 (-0.2342249E+00) number of electron 86.9999983 magnetization 5.0723215 augmentation part 3.5318236 magnetization 4.1140230 Broyden mixing: rms(total) = 0.14224E+00 rms(broyden)= 0.14209E+00 rms(prec ) = 0.14690E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3747 3.0253 2.5029 1.2189 0.7916 0.6415 0.7215 0.7215 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8529.15515868 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 278.89075275 PAW double counting = 7476.87949661 -7456.65597170 entropy T*S EENTRO = 0.02904098 eigenvalues EBANDS = -1032.40788194 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -188.21098069 eV energy without entropy = -188.24002167 energy(sigma->0) = -188.22066101 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 424 total energy-change (2. order) :-0.1103045E+01 (-0.2551668E-01) number of electron 86.9999983 magnetization 3.4757366 augmentation part 3.5292830 magnetization 2.5048897 Broyden mixing: rms(total) = 0.93160E-01 rms(broyden)= 0.93132E-01 rms(prec ) = 0.97048E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4670 4.2176 2.2985 1.3316 0.8374 0.8374 0.8421 0.6188 0.7530 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8530.38522744 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.92959575 PAW double counting = 7450.06688003 -7429.84079818 entropy T*S EENTRO = 0.03512360 eigenvalues EBANDS = -1031.32834033 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -189.31402528 eV energy without entropy = -189.34914888 energy(sigma->0) = -189.32573315 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.1040361E+01 (-0.2356447E-01) number of electron 86.9999983 magnetization 2.7722780 augmentation part 3.5330180 magnetization 1.8243024 Broyden mixing: rms(total) = 0.73146E-01 rms(broyden)= 0.73006E-01 rms(prec ) = 0.77935E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4597 4.6354 2.1266 1.7445 0.9604 0.9604 0.7023 0.7023 0.6893 0.6159 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8530.12752482 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.87405329 PAW double counting = 7410.42642532 -7390.19592029 entropy T*S EENTRO = 0.03311891 eigenvalues EBANDS = -1031.57328015 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -190.35438644 eV energy without entropy = -190.38750536 energy(sigma->0) = -190.36542608 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 392 total energy-change (2. order) :-0.3907845E+00 (-0.5158773E-02) number of electron 86.9999983 magnetization 1.9382000 augmentation part 3.5327338 magnetization 1.0143619 Broyden mixing: rms(total) = 0.44069E-01 rms(broyden)= 0.44012E-01 rms(prec ) = 0.46434E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5087 5.3487 2.0030 2.0030 1.1404 0.9327 0.9327 0.7132 0.7132 0.6798 0.6199 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8530.62746138 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.50330700 PAW double counting = 7404.73296062 -7384.51662640 entropy T*S EENTRO = 0.03551400 eigenvalues EBANDS = -1031.08160605 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -190.74517091 eV energy without entropy = -190.78068491 energy(sigma->0) = -190.75700891 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.5965728E+00 (-0.4151560E-02) number of electron 86.9999983 magnetization 1.4458335 augmentation part 3.5305070 magnetization 0.5348026 Broyden mixing: rms(total) = 0.29418E-01 rms(broyden)= 0.29381E-01 rms(prec ) = 0.31285E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5702 5.8499 2.3314 2.3314 1.2967 1.1423 0.6958 0.6958 0.8745 0.7706 0.6240 0.6603 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8530.08240599 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.85231541 PAW double counting = 7385.58295143 -7365.37130368 entropy T*S EENTRO = 0.03570415 eigenvalues EBANDS = -1031.56774636 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.34174376 eV energy without entropy = -191.37744791 energy(sigma->0) = -191.35364515 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.3713424E+00 (-0.1468102E-02) number of electron 86.9999983 magnetization 1.2079893 augmentation part 3.5278983 magnetization 0.3093667 Broyden mixing: rms(total) = 0.15847E-01 rms(broyden)= 0.15811E-01 rms(prec ) = 0.17084E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5912 6.3048 2.5476 2.5476 1.3664 0.9924 0.9924 0.6895 0.6895 0.9387 0.7607 0.6325 0.6325 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8529.99092816 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.49798392 PAW double counting = 7394.70762560 -7374.50150555 entropy T*S EENTRO = 0.03576293 eigenvalues EBANDS = -1031.67076623 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.71308620 eV energy without entropy = -191.74884913 energy(sigma->0) = -191.72500718 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.2030564E+00 (-0.3867676E-03) number of electron 86.9999983 magnetization 1.0605790 augmentation part 3.5291531 magnetization 0.1693803 Broyden mixing: rms(total) = 0.91941E-02 rms(broyden)= 0.91804E-02 rms(prec ) = 0.10084E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6211 6.6478 2.9100 2.3677 1.4209 1.3083 1.3083 0.6885 0.6885 0.8712 0.8712 0.7243 0.6201 0.6471 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8529.40144365 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.26934693 PAW double counting = 7396.65186093 -7376.44289153 entropy T*S EENTRO = 0.03587011 eigenvalues EBANDS = -1032.23762668 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.91614259 eV energy without entropy = -191.95201270 energy(sigma->0) = -191.92809930 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.1396296E+00 (-0.3011927E-03) number of electron 86.9999983 magnetization 1.0190114 augmentation part 3.5305851 magnetization 0.1305486 Broyden mixing: rms(total) = 0.58654E-02 rms(broyden)= 0.58562E-02 rms(prec ) = 0.65557E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6348 6.8471 3.4654 2.1121 2.1121 1.4219 0.9486 0.9486 0.6878 0.6878 0.8600 0.8600 0.6204 0.6578 0.6578 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.71978613 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.10191440 PAW double counting = 7400.42875982 -7380.21471114 entropy T*S EENTRO = 0.03587734 eigenvalues EBANDS = -1032.89656776 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.05577217 eV energy without entropy = -192.09164951 energy(sigma->0) = -192.06773128 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.6061386E-01 (-0.8227795E-04) number of electron 86.9999983 magnetization 1.0060083 augmentation part 3.5298210 magnetization 0.1176813 Broyden mixing: rms(total) = 0.28515E-02 rms(broyden)= 0.28490E-02 rms(prec ) = 0.34235E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6344 6.9463 3.7068 2.2280 2.2280 1.4354 1.0381 1.0381 0.6904 0.6904 0.9148 0.9148 0.7382 0.7004 0.6228 0.6228 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.62459985 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.04714898 PAW double counting = 7405.27804130 -7385.06424386 entropy T*S EENTRO = 0.03587038 eigenvalues EBANDS = -1032.99734428 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.11638603 eV energy without entropy = -192.15225641 energy(sigma->0) = -192.12834282 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) :-0.1109290E-01 (-0.3247690E-04) number of electron 86.9999983 magnetization 1.0011312 augmentation part 3.5290546 magnetization 0.1120828 Broyden mixing: rms(total) = 0.16633E-02 rms(broyden)= 0.16603E-02 rms(prec ) = 0.21293E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6133 6.9980 3.8466 2.3337 2.3337 1.3090 1.0275 1.0275 1.0930 0.6908 0.6908 0.9045 0.8185 0.8185 0.6778 0.6217 0.6217 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.54584725 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.03685147 PAW double counting = 7405.95693943 -7385.74370431 entropy T*S EENTRO = 0.03587894 eigenvalues EBANDS = -1033.07633850 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.12747892 eV energy without entropy = -192.16335787 energy(sigma->0) = -192.13943857 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.2789606E-02 (-0.7188826E-05) number of electron 86.9999983 magnetization 0.9992482 augmentation part 3.5291484 magnetization 0.1094322 Broyden mixing: rms(total) = 0.10542E-02 rms(broyden)= 0.10533E-02 rms(prec ) = 0.14549E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6225 7.0244 3.9988 2.4004 2.4004 1.5243 1.5243 0.9953 0.9953 0.6904 0.6904 0.9279 0.8837 0.8837 0.7443 0.6583 0.6202 0.6202 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.43198039 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.03025405 PAW double counting = 7405.18494863 -7384.97149384 entropy T*S EENTRO = 0.03589611 eigenvalues EBANDS = -1033.18663439 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13026853 eV energy without entropy = -192.16616464 energy(sigma->0) = -192.14223390 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.1077085E-02 (-0.6270398E-05) number of electron 86.9999983 magnetization 0.9987631 augmentation part 3.5293006 magnetization 0.1083063 Broyden mixing: rms(total) = 0.54604E-03 rms(broyden)= 0.54451E-03 rms(prec ) = 0.87066E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6130 7.0378 4.1459 2.4825 2.4825 1.6194 1.6194 1.0256 1.0256 1.0707 0.6906 0.6906 0.8441 0.8441 0.8417 0.7153 0.6593 0.6198 0.6198 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.33638570 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02687013 PAW double counting = 7404.64682038 -7384.43352197 entropy T*S EENTRO = 0.03590237 eigenvalues EBANDS = -1033.27977212 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13134562 eV energy without entropy = -192.16724798 energy(sigma->0) = -192.14331307 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 328 total energy-change (2. order) :-0.2328675E-03 (-0.1601958E-05) number of electron 86.9999983 magnetization 0.9994340 augmentation part 3.5293149 magnetization 0.1086326 Broyden mixing: rms(total) = 0.38243E-03 rms(broyden)= 0.38153E-03 rms(prec ) = 0.63837E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6114 7.0357 4.2505 2.5463 2.5463 1.7605 1.7605 1.0454 1.0454 1.1137 0.6905 0.6905 0.9777 0.9062 0.9062 0.7700 0.6856 0.6515 0.6173 0.6173 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.28945145 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02548129 PAW double counting = 7404.36301181 -7384.14985656 entropy T*S EENTRO = 0.03590134 eigenvalues EBANDS = -1033.32540621 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13157848 eV energy without entropy = -192.16747983 energy(sigma->0) = -192.14354560 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 288 total energy-change (2. order) :-0.3847417E-03 (-0.1962271E-05) number of electron 86.9999983 magnetization 0.9999385 augmentation part 3.5293225 magnetization 0.1088301 Broyden mixing: rms(total) = 0.26218E-03 rms(broyden)= 0.26117E-03 rms(prec ) = 0.40685E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6226 7.0365 4.3871 2.8526 2.3746 2.3746 1.5500 1.2575 1.2575 1.0302 1.0302 0.6906 0.6906 0.8904 0.8904 0.8222 0.7593 0.6811 0.6416 0.6171 0.6171 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.23931477 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02395353 PAW double counting = 7404.05301701 -7383.84005103 entropy T*S EENTRO = 0.03590267 eigenvalues EBANDS = -1033.37421194 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13196322 eV energy without entropy = -192.16786590 energy(sigma->0) = -192.14393078 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 304 total energy-change (2. order) :-0.2927978E-03 (-0.8031407E-06) number of electron 86.9999983 magnetization 0.9998906 augmentation part 3.5293299 magnetization 0.1087544 Broyden mixing: rms(total) = 0.19677E-03 rms(broyden)= 0.19627E-03 rms(prec ) = 0.27929E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6189 7.0389 4.5096 3.0909 2.4478 2.4478 1.4559 1.4559 1.0335 1.0335 0.6906 0.6906 1.0562 1.0562 1.0258 0.8463 0.8463 0.7317 0.6717 0.6265 0.6265 0.6138 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.21772375 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02326506 PAW double counting = 7404.03374754 -7383.82085292 entropy T*S EENTRO = 0.03590470 eigenvalues EBANDS = -1033.39533794 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13225602 eV energy without entropy = -192.16816072 energy(sigma->0) = -192.14422426 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 296 total energy-change (2. order) :-0.1884850E-03 (-0.2803065E-06) number of electron 86.9999983 magnetization 0.9998817 augmentation part 3.5293373 magnetization 0.1088627 Broyden mixing: rms(total) = 0.12446E-03 rms(broyden)= 0.12422E-03 rms(prec ) = 0.18050E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6642 7.0396 4.8174 3.5863 2.5101 2.5101 1.7417 1.7417 1.2410 1.2410 1.0243 1.0243 0.6906 0.6906 0.8941 0.8941 0.9017 0.8196 0.7134 0.6727 0.6208 0.6208 0.6155 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.20901288 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02260759 PAW double counting = 7404.10147042 -7383.88856022 entropy T*S EENTRO = 0.03590584 eigenvalues EBANDS = -1033.40359656 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13244451 eV energy without entropy = -192.16835035 energy(sigma->0) = -192.14441312 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 296 total energy-change (2. order) : 0.2821908E-05 (-0.1335016E-03) number of electron 86.9999983 magnetization 0.9999330 augmentation part 3.5293328 magnetization 0.1090225 Broyden mixing: rms(total) = 0.11200E-03 rms(broyden)= 0.11182E-03 rms(prec ) = 0.13825E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6513 7.0383 5.1078 3.8362 2.5508 2.5508 1.7876 1.7876 1.2332 1.2332 1.0204 1.0204 0.6906 0.6906 0.9272 0.9272 0.8654 0.8654 0.7085 0.7085 0.6415 0.6203 0.6039 0.5653 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.20537410 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02211497 PAW double counting = 7404.18507391 -7383.97214089 entropy T*S EENTRO = 0.03590698 eigenvalues EBANDS = -1033.40676384 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13244169 eV energy without entropy = -192.16834866 energy(sigma->0) = -192.14441068 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.1606773E-03 (-0.2699711E-03) number of electron 86.9999983 magnetization 0.9999518 augmentation part 3.5293480 magnetization 0.1090771 Broyden mixing: rms(total) = 0.15503E-02 rms(broyden)= 0.15501E-02 rms(prec ) = 0.15523E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5954 7.0369 5.2414 3.9012 2.5499 2.5499 1.7706 1.7706 1.2398 1.2398 1.0203 1.0203 0.6906 0.6906 0.9518 0.9518 0.8623 0.8623 0.7303 0.6921 0.6513 0.6098 0.6098 0.6276 0.0194 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.20793392 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02215896 PAW double counting = 7404.21926804 -7384.00633394 entropy T*S EENTRO = 0.03590710 eigenvalues EBANDS = -1033.40440989 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13260236 eV energy without entropy = -192.16850946 energy(sigma->0) = -192.14457140 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 368 total energy-change (2. order) : 0.1011653E-05 (-0.1205429E-03) number of electron 86.9999983 magnetization 0.9999596 augmentation part 3.5293263 magnetization 0.1090600 Broyden mixing: rms(total) = 0.42649E-03 rms(broyden)= 0.42616E-03 rms(prec ) = 0.43096E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5350 7.0368 5.2568 3.9085 2.5465 2.5465 1.7701 1.7701 1.2458 1.2458 1.0208 1.0208 0.6906 0.6906 0.9451 0.9451 0.8653 0.8653 0.7250 0.6972 0.6538 0.6229 0.6097 0.6097 0.0642 0.0233 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.20780892 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02214786 PAW double counting = 7404.19850894 -7383.98557739 entropy T*S EENTRO = 0.03590710 eigenvalues EBANDS = -1033.40452023 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13260135 eV energy without entropy = -192.16850845 energy(sigma->0) = -192.14457038 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 472 total energy-change (2. order) :-0.1205782E+04 (-0.1198838E+04) number of electron 87.0059636 magnetization 0.9936051 augmentation part 3.0723092 magnetization -2.0731697 Broyden mixing: rms(total) = 0.45533E+02 rms(broyden)= 0.45533E+02 rms(prec ) = 0.45535E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4760 7.0369 5.2551 3.9085 2.5460 2.5460 1.7712 1.7712 1.2465 1.2465 1.0208 1.0208 0.6906 0.6906 0.9446 0.9446 0.8657 0.8657 0.7244 0.6998 0.6533 0.6222 0.6098 0.6098 0.0606 0.0000 0.0236 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.20839674 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02214889 PAW double counting = 7404.19278453 -7383.97985728 entropy T*S EENTRO = -0.04080015 eigenvalues EBANDS = -2239.10876947 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1397.91414893 eV energy without entropy = -1397.87334878 energy(sigma->0) = -1397.90054888 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 504 total energy-change (2. order) :-0.1767713E+06 (-0.1638696E+06) number of electron 87.4433049 magnetization 0.9884267 augmentation part 0.6488036 magnetization -38.3103045 Broyden mixing: rms(total) = 0.57693E+02 rms(broyden)= 0.57692E+02 rms(prec ) = 0.57821E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4260 7.0370 5.2795 3.9197 2.5441 2.5441 1.7830 1.7830 1.2539 1.2539 1.0222 1.0222 0.6906 0.6906 0.9488 0.9488 0.8670 0.8670 0.7393 0.7015 0.6497 0.6217 0.6133 0.6133 0.0864 0.0210 0.0002 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.16856014 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.14205667 PAW double counting = 7412.28043344 -7392.06398358 entropy T*S EENTRO = -0.03335076 eigenvalues EBANDS = -179010.60038814 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -178169.23505122 eV energy without entropy = -178169.20170046 energy(sigma->0) = -178169.22393430 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 456 total energy-change (2. order) :-0.5441207E+06 (-0.7223133E+06) number of electron 87.7601617 magnetization 0.9968569 augmentation part 0.9613740 magnetization -23.4709616 Broyden mixing: rms(total) = 0.74647E+02 rms(broyden)= 0.74647E+02 rms(prec ) = 0.74719E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3771 7.0373 5.2858 3.9263 2.5485 2.5485 1.7832 1.7832 1.2577 1.2577 1.0226 1.0226 0.6906 0.6906 0.9511 0.9511 0.8644 0.8644 0.7343 0.7045 0.6433 0.6217 0.6110 0.6110 0.1243 0.0217 0.0002 0.0002 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.26805438 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.14750740 PAW double counting = 7408.83172085 -7388.61632758 entropy T*S EENTRO = -0.03322073 eigenvalues EBANDS = -723132.24817135 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -722289.97780450 eV energy without entropy = -722289.94458378 energy(sigma->0) = -722289.96673093 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 456 total energy-change (2. order) : 0.6634818E+06 (-0.5590640E+05) number of electron 87.9784794 magnetization 0.9940746 augmentation part 1.3376039 magnetization -17.6586904 Broyden mixing: rms(total) = 0.67061E+02 rms(broyden)= 0.67061E+02 rms(prec ) = 0.67118E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3295 7.0373 5.2859 3.9263 2.5486 2.5486 1.7831 1.7831 1.2577 1.2577 1.0225 1.0225 0.6906 0.6906 0.9513 0.9513 0.8643 0.8643 0.7342 0.7047 0.6434 0.6218 0.6110 0.6110 0.1222 0.0217 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.06238368 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.20134047 PAW double counting = 7406.22256191 -7386.00680837 entropy T*S EENTRO = -0.00748527 eigenvalues EBANDS = -59650.77161765 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -58808.21565131 eV energy without entropy = -58808.20816603 energy(sigma->0) = -58808.21315621 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 448 total energy-change (2. order) :-0.5352663E+05 (-0.9387321E+05) number of electron 87.6282204 magnetization 0.9907899 augmentation part 0.8414985 magnetization -24.0368544 Broyden mixing: rms(total) = 0.85408E+02 rms(broyden)= 0.85408E+02 rms(prec ) = 0.85476E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2853 7.0371 5.2892 3.9267 2.5449 2.5449 1.7845 1.7845 1.2573 1.2573 1.0227 1.0227 0.6906 0.6906 0.9490 0.9490 0.8667 0.8667 0.7331 0.7058 0.6448 0.6177 0.6127 0.6127 0.1258 0.0215 0.0003 0.0003 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8527.98463939 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.18895799 PAW double counting = 7404.13695483 -7383.92086457 entropy T*S EENTRO = -0.04465310 eigenvalues EBANDS = -113177.42696168 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -112334.84246462 eV energy without entropy = -112334.79781152 energy(sigma->0) = -112334.82758025 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 440 total energy-change (2. order) :-0.2966926E+05 (-0.1235830E+06) number of electron 87.9258550 magnetization 0.9930579 augmentation part 0.7490411 magnetization -28.3266542 Broyden mixing: rms(total) = 0.11852E+03 rms(broyden)= 0.11852E+03 rms(prec ) = 0.11858E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2442 7.0370 5.2936 3.9279 2.5449 2.5449 1.7842 1.7842 1.2574 1.2574 1.0225 1.0225 0.6906 0.6906 0.9500 0.9500 0.8662 0.8662 0.7361 0.7037 0.6449 0.6188 0.6134 0.6134 0.1264 0.0216 0.0003 0.0003 0.0002 0.0002 0.0001 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.03984917 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.95722731 PAW double counting = 7405.15175348 -7384.93557152 entropy T*S EENTRO = -0.03005958 eigenvalues EBANDS = -142846.41174278 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -142004.09950097 eV energy without entropy = -142004.06944138 energy(sigma->0) = -142004.08948111 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) : 0.7762858E+05 (-0.6139464E+05) number of electron 88.4247109 magnetization 0.9871414 augmentation part 1.1522062 magnetization -21.4693803 Broyden mixing: rms(total) = 0.97984E+02 rms(broyden)= 0.97984E+02 rms(prec ) = 0.98023E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2066 7.0370 5.3045 3.9342 2.5437 2.5437 1.7854 1.7854 1.2579 1.2579 1.0223 1.0223 0.6906 0.6906 0.9480 0.9480 0.8682 0.8682 0.7305 0.7119 0.6439 0.6145 0.6145 0.6145 0.1528 0.0206 0.0003 0.0003 0.0003 0.0003 0.0003 0.0002 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8527.96140022 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.14258298 PAW double counting = 7405.22376489 -7385.00728323 entropy T*S EENTRO = 0.00841920 eigenvalues EBANDS = -65218.13669655 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -64375.52187163 eV energy without entropy = -64375.53029084 energy(sigma->0) = -64375.52467804 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 472 total energy-change (2. order) :-0.4910001E+05 (-0.8906969E+05) number of electron 87.4818510 magnetization 0.9710204 augmentation part 0.6483454 magnetization -34.8910084 Broyden mixing: rms(total) = 0.13704E+03 rms(broyden)= 0.13704E+03 rms(prec ) = 0.13710E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1705 7.0368 5.3094 3.9346 2.5430 2.5430 1.7857 1.7857 1.2577 1.2577 1.0223 1.0223 0.6906 0.6906 0.9419 0.9419 0.8732 0.8732 0.7348 0.7101 0.6448 0.6135 0.6142 0.6142 0.1608 0.0206 0.0004 0.0004 0.0004 0.0003 0.0003 0.0000 0.0003 0.0003 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.00467707 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.34185538 PAW double counting = 7414.77594765 -7394.55603006 entropy T*S EENTRO = -0.02155978 eigenvalues EBANDS = -114318.27797712 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -113475.53369971 eV energy without entropy = -113475.51213993 energy(sigma->0) = -113475.52651312 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.9565548E+05 (-0.1960691E+06) number of electron 87.6559469 magnetization 0.9694954 augmentation part 0.2732369 magnetization -44.2063120 Broyden mixing: rms(total) = 0.21663E+03 rms(broyden)= 0.21663E+03 rms(prec ) = 0.21668E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1359 7.0368 5.3096 3.9359 2.5414 2.5414 1.7875 1.7875 1.2571 1.2571 1.0224 1.0224 0.6906 0.6906 0.9389 0.9389 0.8760 0.8760 0.7322 0.7128 0.6451 0.6118 0.6147 0.6147 0.1541 0.0206 0.0005 0.0005 0.0005 0.0003 0.0003 0.0002 0.0002 0.0000 0.0003 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.08770035 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.91807959 PAW double counting = 7415.99790561 -7395.77781866 entropy T*S EENTRO = 0.00379926 eigenvalues EBANDS = -209973.27467913 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -209131.01167240 eV energy without entropy = -209131.01547166 energy(sigma->0) = -209131.01293882 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 488 total energy-change (2. order) : 0.1343696E+06 (-0.7382722E+05) number of electron 87.9995319 magnetization 1.0027219 augmentation part 0.5629769 magnetization -38.7444879 Broyden mixing: rms(total) = 0.20177E+03 rms(broyden)= 0.20177E+03 rms(prec ) = 0.20181E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1037 7.0370 5.3084 3.9378 2.5405 2.5405 1.7891 1.7891 1.2584 1.2584 1.0229 1.0229 0.6906 0.6906 0.9450 0.9450 0.8716 0.8716 0.7276 0.7170 0.6448 0.6112 0.6147 0.6147 0.1558 0.0208 0.0005 0.0005 0.0005 0.0004 0.0004 0.0004 0.0004 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.49401833 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.19113921 PAW double counting = 7416.19435722 -7395.97590002 entropy T*S EENTRO = 0.00646990 eigenvalues EBANDS = -75603.58244007 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -74761.45165081 eV energy without entropy = -74761.45812070 energy(sigma->0) = -74761.45380744 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 416 total energy-change (2. order) :-0.3032727E+05 (-0.1070157E+06) number of electron 88.4905692 magnetization 0.9333769 augmentation part 0.9811238 magnetization -22.7817612 Broyden mixing: rms(total) = 0.17394E+03 rms(broyden)= 0.17394E+03 rms(prec ) = 0.17397E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0731 7.0369 5.3079 3.9378 2.5395 2.5395 1.7897 1.7897 1.2583 1.2583 1.0229 1.0229 0.6906 0.6906 0.9454 0.9454 0.8717 0.8717 0.7244 0.7199 0.6450 0.6107 0.6151 0.6151 0.1585 0.0207 0.0006 0.0006 0.0004 0.0004 0.0004 0.0004 0.0005 0.0004 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8528.24738163 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.37696773 PAW double counting = 7413.40982306 -7393.19094607 entropy T*S EENTRO = -0.03110718 eigenvalues EBANDS = -105931.24821188 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -105088.72211467 eV energy without entropy = -105088.69100749 energy(sigma->0) = -105088.71174561 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 344 total energy-change (2. order) : 0.1048934E+06 (-0.1588344E+04) number of electron 88.3745718 magnetization -0.2682131 augmentation part 3.4671408 magnetization -0.7996294 Broyden mixing: rms(total) = 0.86003E+01 rms(broyden)= 0.85978E+01 rms(prec ) = 0.86186E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0445 7.0368 5.3231 3.9432 2.5385 2.5385 1.7889 1.7889 1.2515 1.2515 1.0236 1.0236 0.6906 0.6906 0.9400 0.9400 0.8775 0.8775 0.7212 0.7212 0.6470 0.6122 0.6164 0.6164 0.1574 0.0197 0.0026 0.0026 0.0010 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8529.84562857 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.63151939 PAW double counting = 7418.22708975 -7398.01283870 entropy T*S EENTRO = 0.02119425 eigenvalues EBANDS = -1036.59674814 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -195.36667072 eV energy without entropy = -195.38786497 energy(sigma->0) = -195.37373547 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.3028906E+04 (-0.2212901E+04) number of electron 88.4212518 magnetization -0.6109569 augmentation part 1.8918586 magnetization -18.3734079 Broyden mixing: rms(total) = 0.56513E+02 rms(broyden)= 0.56513E+02 rms(prec ) = 0.56554E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0173 7.0368 5.3224 3.9429 2.5385 2.5385 1.7887 1.7887 1.2519 1.2519 1.0236 1.0236 0.6906 0.6906 0.9405 0.9405 0.8769 0.8769 0.7212 0.7212 0.6469 0.6123 0.6165 0.6165 0.1587 0.0192 0.0084 0.0084 0.0010 0.0012 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8504.08047504 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.72336177 PAW double counting = 7979.55014952 -7959.28725100 entropy T*S EENTRO = -0.01088073 eigenvalues EBANDS = -4092.37638996 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -3224.27274415 eV energy without entropy = -3224.26186342 energy(sigma->0) = -3224.26911724 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 456 total energy-change (2. order) :-0.3128674E+04 (-0.6503826E+04) number of electron 88.5217671 magnetization -0.3329969 augmentation part 2.0267040 magnetization -17.3342008 Broyden mixing: rms(total) = 0.67166E+02 rms(broyden)= 0.67166E+02 rms(prec ) = 0.67199E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9984 7.0374 5.3554 3.9626 2.5481 2.5481 1.7817 1.7817 1.2602 1.2602 1.0236 1.0236 0.6906 0.6906 0.9508 0.9508 0.8673 0.8673 0.7240 0.7240 0.6462 0.6111 0.6193 0.6193 0.2387 0.0460 0.0460 0.0413 0.0160 0.0013 0.0010 0.0008 0.0008 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 0.0007 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8510.55549080 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 278.04800225 PAW double counting = 7996.82660752 -7976.59709155 entropy T*S EENTRO = -0.04020672 eigenvalues EBANDS = -7214.83695565 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -6352.94639364 eV energy without entropy = -6352.90618693 energy(sigma->0) = -6352.93299140 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 424 total energy-change (2. order) : 0.3151091E+04 (-0.3335980E+04) number of electron 88.1640173 magnetization -0.0222515 augmentation part 2.5101747 magnetization -11.7127967 Broyden mixing: rms(total) = 0.41845E+02 rms(broyden)= 0.41845E+02 rms(prec ) = 0.41874E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9751 7.0372 5.3668 3.9671 2.5491 2.5491 1.7845 1.7845 1.2627 1.2627 1.0240 1.0240 0.6906 0.6906 0.9343 0.9343 0.8831 0.8831 0.7290 0.7290 0.6440 0.6226 0.6226 0.6059 0.2503 0.0398 0.0398 0.0347 0.0347 0.0163 0.0013 0.0010 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8503.52155058 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.59813325 PAW double counting = 7979.23137054 -7958.98366635 entropy T*S EENTRO = -0.00992564 eigenvalues EBANDS = -4070.37851628 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -3201.85541376 eV energy without entropy = -3201.84548812 energy(sigma->0) = -3201.85210521 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.3440598E+05 (-0.5679702E+04) number of electron 88.2160528 magnetization -0.4624576 augmentation part 1.4108991 magnetization -36.3772906 Broyden mixing: rms(total) = 0.95446E+02 rms(broyden)= 0.95446E+02 rms(prec ) = 0.95512E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9521 7.0373 5.3664 3.9685 2.5507 2.5507 1.7852 1.7852 1.2633 1.2633 1.0238 1.0238 0.6906 0.6906 0.9433 0.9433 0.8742 0.8742 0.7289 0.7289 0.6440 0.6224 0.6224 0.6059 0.2532 0.0397 0.0371 0.0371 0.0363 0.0248 0.0161 0.0013 0.0010 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8512.33679261 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.58197790 PAW double counting = 7892.34662951 -7872.16784225 entropy T*S EENTRO = -0.00904520 eigenvalues EBANDS = -38467.45699308 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -37607.83332444 eV energy without entropy = -37607.82427923 energy(sigma->0) = -37607.83030937 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 504 total energy-change (2. order) :-0.3608170E+05 (-0.7217261E+05) number of electron 88.2012718 magnetization -0.4865476 augmentation part 0.6982073 magnetization -44.2292872 Broyden mixing: rms(total) = 0.15133E+03 rms(broyden)= 0.15133E+03 rms(prec ) = 0.15138E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9313 7.0373 5.3739 3.9707 2.5515 2.5515 1.7842 1.7842 1.2647 1.2647 1.0239 1.0239 0.6906 0.6906 0.9442 0.9442 0.8736 0.8736 0.7278 0.7278 0.6432 0.6215 0.6215 0.6060 0.2638 0.0472 0.0472 0.0439 0.0381 0.0381 0.0162 0.0183 0.0013 0.0010 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8515.71409199 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 278.08703317 PAW double counting = 7938.57401123 -7918.38799412 entropy T*S EENTRO = -0.01050938 eigenvalues EBANDS = -74546.29309630 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -73689.53590610 eV energy without entropy = -73689.52539672 energy(sigma->0) = -73689.53240298 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.8997760E+04 (-0.8233479E+05) number of electron 87.9003470 magnetization 0.1182358 augmentation part 0.5488023 magnetization -44.9326612 Broyden mixing: rms(total) = 0.14325E+03 rms(broyden)= 0.14325E+03 rms(prec ) = 0.14331E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9097 7.0373 5.3738 3.9706 2.5515 2.5515 1.7841 1.7841 1.2648 1.2648 1.0239 1.0239 0.6906 0.6906 0.9441 0.9441 0.8737 0.8737 0.7278 0.7278 0.6432 0.6215 0.6215 0.6060 0.2640 0.0476 0.0476 0.0442 0.0381 0.0381 0.0184 0.0162 0.0013 0.0010 0.0011 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0003 0.0003 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8513.73702648 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.75200538 PAW double counting = 7886.25458733 -7866.05028039 entropy T*S EENTRO = -0.03091637 eigenvalues EBANDS = -83545.69347520 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -82687.29636444 eV energy without entropy = -82687.26544807 energy(sigma->0) = -82687.28605898 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 504 total energy-change (2. order) :-0.1270123E+06 (-0.1985833E+06) number of electron 87.1522900 magnetization 0.0987121 augmentation part -0.4364013 magnetization -57.5514111 Broyden mixing: rms(total) = 0.22807E+03 rms(broyden)= 0.22807E+03 rms(prec ) = 0.22814E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8913 7.0375 5.3798 3.9735 2.5559 2.5559 1.7851 1.7851 1.2643 1.2643 1.0241 1.0241 0.6906 0.6906 0.9446 0.9446 0.8733 0.8733 0.7293 0.7293 0.6447 0.6189 0.6189 0.6081 0.2653 0.0585 0.0480 0.0480 0.0496 0.0496 0.0222 0.0222 0.0162 0.0163 0.0013 0.0010 0.0008 0.0008 0.0007 0.0005 0.0005 0.0005 0.0000 0.0003 0.0003 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8513.85928681 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.45004915 PAW double counting = 7785.12695272 -7765.03951306 entropy T*S EENTRO = -0.04954335 eigenvalues EBANDS = -210557.47496307 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -209699.63756313 eV energy without entropy = -209699.58801978 energy(sigma->0) = -209699.62104868 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 488 total energy-change (2. order) :-0.4321197E+06 (-0.6391163E+06) number of electron 86.9459119 magnetization -1.5978594 augmentation part -1.6995297 magnetization -68.8261140 Broyden mixing: rms(total) = 0.60942E+03 rms(broyden)= 0.60942E+03 rms(prec ) = 0.60945E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8264 6.9847 5.3081 3.9224 2.6517 1.9979 1.3739 1.2515 1.2515 1.0302 1.0302 0.9338 0.6933 0.6933 0.7559 0.7028 0.6453 0.6172 0.6094 0.2569 0.0617 0.0509 0.0509 0.0452 0.0452 0.0247 0.0247 0.0194 0.0151 0.0012 0.0009 0.0009 0.0007 0.0007 0.0005 0.0005 0.0003 0.0003 0.0000 0.0003 0.0002 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5563.78049948 -Hartree energ DENC = -8514.51143762 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.08345562 PAW double counting = 7797.38349224 -7777.29629475 entropy T*S EENTRO = -0.02447285 eigenvalues EBANDS = -642676.14105798 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -641819.29757406 eV energy without entropy = -641819.27310121 energy(sigma->0) = -641819.28941644 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- ----------------------------------------------------------------------------- | | | 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 = 12 2 16 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------