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 14:36:42 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.687 0.463 0.594- 2 1.04 3 1.05 2 0.710 0.458 0.656- 1 1.04 3 0.726 0.449 0.550- 1 1.05 4 0.152 0.517 0.368- 25 1.10 5 0.208 0.610 0.368- 25 1.10 6 0.222 0.519 0.293- 25 1.10 7 0.112 0.548 0.658- 26 1.10 8 0.151 0.637 0.598- 26 1.11 9 0.108 0.551 0.539- 26 1.11 10 0.420 0.430 0.295- 27 1.10 11 0.420 0.340 0.373- 27 1.11 12 0.347 0.409 0.360- 27 1.11 13 0.536 0.278 0.573- 28 1.10 14 0.535 0.349 0.478- 28 1.10 15 0.562 0.392 0.582- 28 1.10 16 0.248 0.689 0.505- 29 1.10 17 0.306 0.666 0.591- 29 1.10 18 0.331 0.655 0.478- 29 1.10 19 0.442 0.621 0.559- 30 1.10 20 0.451 0.557 0.658- 30 1.10 21 0.514 0.554 0.575- 30 1.10 22 0.326 0.444 0.642- 24 1.11 23 0.311 0.377 0.547- 24 1.11 24 0.329 0.444 0.568- 22 1.11 23 1.11 32 1.87 31 1.88 25 0.205 0.537 0.361- 6 1.10 4 1.10 5 1.10 33 1.42 26 0.140 0.565 0.597- 7 1.10 9 1.11 8 1.11 34 1.42 27 0.402 0.409 0.361- 10 1.10 11 1.11 12 1.11 35 1.42 28 0.526 0.347 0.550- 14 1.10 15 1.10 13 1.10 36 1.43 29 0.289 0.645 0.525- 18 1.10 16 1.10 17 1.10 32 1.88 30 0.460 0.557 0.586- 19 1.10 21 1.10 20 1.10 31 1.86 31 0.419 0.458 0.535- 35 1.66 36 1.66 30 1.86 24 1.88 32 0.264 0.524 0.527- 33 1.67 34 1.67 24 1.87 29 1.88 33 0.245 0.490 0.423- 25 1.42 32 1.67 34 0.199 0.512 0.595- 26 1.42 32 1.67 35 0.427 0.469 0.426- 27 1.42 31 1.66 36 0.458 0.366 0.570- 28 1.43 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.686613100 0.463358430 0.593958480 0.710088430 0.458245080 0.655632020 0.725763080 0.449150570 0.549804740 0.151923900 0.516767120 0.367725210 0.208334250 0.609629990 0.367795840 0.221539190 0.518696220 0.293110410 0.111585600 0.547929020 0.657507560 0.150964930 0.636685430 0.597993040 0.107577330 0.550507190 0.538546740 0.420342410 0.429752640 0.294851770 0.419582510 0.339632820 0.372939680 0.346845370 0.408642310 0.359818630 0.535824640 0.277915490 0.572704120 0.534812770 0.349039380 0.478151200 0.562062080 0.391803460 0.582069300 0.247851870 0.688941800 0.504519100 0.305721140 0.666220290 0.591052450 0.330511290 0.655289320 0.477558640 0.441992060 0.620871170 0.559011580 0.451317050 0.557424350 0.658271430 0.513974080 0.553584890 0.574687850 0.325882290 0.444500000 0.641570420 0.311402070 0.377019670 0.547359660 0.328759950 0.443968510 0.567915250 0.204804420 0.536522320 0.360689760 0.139772570 0.564547780 0.596744870 0.402167900 0.408502490 0.360865800 0.526035880 0.346541280 0.550225190 0.289033350 0.644605820 0.524566950 0.459567270 0.556850990 0.585621080 0.418665130 0.457646640 0.535224900 0.263823350 0.524049040 0.526655920 0.245480380 0.490472170 0.423425300 0.198834340 0.511742400 0.595376270 0.427082000 0.469094240 0.425640030 0.457718050 0.366176570 0.570318440 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.68661310 0.46335843 0.59395848 0.71008843 0.45824508 0.65563202 0.72576308 0.44915057 0.54980474 0.15192390 0.51676712 0.36772521 0.20833425 0.60962999 0.36779584 0.22153919 0.51869622 0.29311041 0.11158560 0.54792902 0.65750756 0.15096493 0.63668543 0.59799304 0.10757733 0.55050719 0.53854674 0.42034241 0.42975264 0.29485177 0.41958251 0.33963282 0.37293968 0.34684537 0.40864231 0.35981863 0.53582464 0.27791549 0.57270412 0.53481277 0.34903938 0.47815120 0.56206208 0.39180346 0.58206930 0.24785187 0.68894180 0.50451910 0.30572114 0.66622029 0.59105245 0.33051129 0.65528932 0.47755864 0.44199206 0.62087117 0.55901158 0.45131705 0.55742435 0.65827143 0.51397408 0.55358489 0.57468785 0.32588229 0.44450000 0.64157042 0.31140207 0.37701967 0.54735966 0.32875995 0.44396851 0.56791525 0.20480442 0.53652232 0.36068976 0.13977257 0.56454778 0.59674487 0.40216790 0.40850249 0.36086580 0.52603588 0.34654128 0.55022519 0.28903335 0.64460582 0.52456695 0.45956727 0.55685099 0.58562108 0.41866513 0.45764664 0.53522490 0.26382335 0.52404904 0.52665592 0.24548038 0.49047217 0.42342530 0.19883434 0.51174240 0.59537627 0.42708200 0.46909424 0.42564003 0.45771805 0.36617657 0.57031844 position of ions in cartesian coordinates (Angst): 13.73226200 6.95037645 8.90937720 14.20176860 6.87367620 9.83448030 14.51526160 6.73725855 8.24707110 3.03847800 7.75150680 5.51587815 4.16668500 9.14444985 5.51693760 4.43078380 7.78044330 4.39665615 2.23171200 8.21893530 9.86261340 3.01929860 9.55028145 8.96989560 2.15154660 8.25760785 8.07820110 8.40684820 6.44628960 4.42277655 8.39165020 5.09449230 5.59409520 6.93690740 6.12963465 5.39727945 10.71649280 4.16873235 8.59056180 10.69625540 5.23559070 7.17226800 11.24124160 5.87705190 8.73103950 4.95703740 10.33412700 7.56778650 6.11442280 9.99330435 8.86578675 6.61022580 9.82933980 7.16337960 8.83984120 9.31306755 8.38517370 9.02634100 8.36136525 9.87407145 10.27948160 8.30377335 8.62031775 6.51764580 6.66750000 9.62355630 6.22804140 5.65529505 8.21039490 6.57519900 6.65952765 8.51872875 4.09608840 8.04783480 5.41034640 2.79545140 8.46821670 8.95117305 8.04335800 6.12753735 5.41298700 10.52071760 5.19811920 8.25337785 5.78066700 9.66908730 7.86850425 9.19134540 8.35276485 8.78431620 8.37330260 6.86469960 8.02837350 5.27646700 7.86073560 7.89983880 4.90960760 7.35708255 6.35137950 3.97668680 7.67613600 8.93064405 8.54164000 7.03641360 6.38460045 9.15436100 5.49264855 8.55477660 -------------------------------------------------------------------------------------------------------- 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 541229. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 4073. 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 2537 Maximum index for augmentation-charges 2259 (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 : 384 total energy-change (2. order) : 0.6884333E+03 (-0.2068665E+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 = 5741.54159136 -Hartree energ DENC = -8307.17406959 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.06252480 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.02426899 eigenvalues EBANDS = -543.17958993 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 688.43332823 eV energy without entropy = 688.40905923 energy(sigma->0) = 688.42523856 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 536 total energy-change (2. order) :-0.5832926E+03 (-0.5253589E+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 = 5741.54159136 -Hartree energ DENC = -8307.17406959 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.06252480 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = -0.00388628 eigenvalues EBANDS = -1126.44402118 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 105.14074170 eV energy without entropy = 105.14462798 energy(sigma->0) = 105.14203713 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 440 total energy-change (2. order) :-0.2514706E+03 (-0.2494164E+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 = 5741.54159136 -Hartree energ DENC = -8307.17406959 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.06252480 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.02012361 eigenvalues EBANDS = -1377.93864783 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -146.32987506 eV energy without entropy = -146.34999867 energy(sigma->0) = -146.33658293 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 432 total energy-change (2. order) :-0.2288330E+02 (-0.2276294E+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 = 5741.54159136 -Hartree energ DENC = -8307.17406959 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.06252480 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.01395088 eigenvalues EBANDS = -1400.81577948 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.21317944 eV energy without entropy = -169.22713032 energy(sigma->0) = -169.21782973 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 472 total energy-change (2. order) :-0.4948457E+00 (-0.4936388E+00) number of electron 87.0000014 magnetization 30.9695109 augmentation part 4.2625937 magnetization 30.1096750 Broyden mixing: rms(total) = 0.41782E+01 rms(broyden)= 0.41755E+01 rms(prec ) = 0.43504E+01 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8307.17406959 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 264.06252480 PAW double counting = 2564.72054280 -2541.99018201 entropy T*S EENTRO = 0.01446685 eigenvalues EBANDS = -1401.31114113 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -169.70802511 eV energy without entropy = -169.72249197 energy(sigma->0) = -169.71284740 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) : 0.9968865E+02 (-0.2438936E+02) number of electron 87.0000018 magnetization 26.4491117 augmentation part 3.8654943 magnetization 25.1414310 Broyden mixing: rms(total) = 0.20341E+01 rms(broyden)= 0.20330E+01 rms(prec ) = 0.21186E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9146 0.9146 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8509.86595341 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 361.19211105 PAW double counting = 4231.58447536 -4211.10727011 entropy T*S EENTRO = 0.01906496 eigenvalues EBANDS = -1193.81164071 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -70.01937969 eV energy without entropy = -70.03844465 energy(sigma->0) = -70.02573467 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 416 total energy-change (2. order) :-0.4981096E+02 (-0.3992076E+01) number of electron 87.0000018 magnetization 22.4380434 augmentation part 3.5907278 magnetization 21.1964133 Broyden mixing: rms(total) = 0.12880E+01 rms(broyden)= 0.12877E+01 rms(prec ) = 0.13266E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.9165 1.0527 0.7804 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8615.93181150 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 325.63525987 PAW double counting = 5738.41124504 -5718.68369680 entropy T*S EENTRO = 0.01712110 eigenvalues EBANDS = -1101.24829265 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -119.83034177 eV energy without entropy = -119.84746286 energy(sigma->0) = -119.83604880 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) :-0.1513999E+02 (-0.6684331E+00) number of electron 87.0000018 magnetization 16.8766330 augmentation part 3.5716807 magnetization 15.7101348 Broyden mixing: rms(total) = 0.82569E+00 rms(broyden)= 0.82562E+00 rms(prec ) = 0.84974E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1668 1.7209 1.1229 0.6567 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8656.77444855 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 316.82056590 PAW double counting = 6649.56118750 -6629.71752152 entropy T*S EENTRO = 0.02005848 eigenvalues EBANDS = -1066.85000751 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -134.97033253 eV energy without entropy = -134.99039100 energy(sigma->0) = -134.97701868 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.2958113E+02 (-0.1303625E+01) number of electron 87.0000018 magnetization 13.7426444 augmentation part 3.5868416 magnetization 12.6356437 Broyden mixing: rms(total) = 0.40997E+00 rms(broyden)= 0.40980E+00 rms(prec ) = 0.42055E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2427 2.2049 1.3263 0.7606 0.6791 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8693.90862033 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 295.97119890 PAW double counting = 7466.29380635 -7446.27424607 entropy T*S EENTRO = 0.02885769 eigenvalues EBANDS = -1038.63229558 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -164.55146587 eV energy without entropy = -164.58032356 energy(sigma->0) = -164.56108510 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 392 total energy-change (2. order) :-0.1194704E+02 (-0.5837318E+00) number of electron 87.0000018 magnetization 10.2145061 augmentation part 3.5598743 magnetization 9.1251038 Broyden mixing: rms(total) = 0.31222E+00 rms(broyden)= 0.31210E+00 rms(prec ) = 0.32151E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2293 1.9234 1.9234 0.8396 0.8396 0.6205 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8708.31706197 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 288.11347585 PAW double counting = 7614.63253582 -7594.56233243 entropy T*S EENTRO = 0.03518142 eigenvalues EBANDS = -1028.37013852 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -176.49850666 eV energy without entropy = -176.53368808 energy(sigma->0) = -176.51023380 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.7447470E+01 (-0.3715822E+00) number of electron 87.0000018 magnetization 6.3431302 augmentation part 3.5359432 magnetization 5.3085742 Broyden mixing: rms(total) = 0.20920E+00 rms(broyden)= 0.20911E+00 rms(prec ) = 0.21518E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4286 2.6020 2.6020 1.1793 0.8670 0.7148 0.6068 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8710.31818273 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 282.32758144 PAW double counting = 7509.04652675 -7488.89977382 entropy T*S EENTRO = 0.03149444 eigenvalues EBANDS = -1028.10345587 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -183.94597661 eV energy without entropy = -183.97747106 energy(sigma->0) = -183.95647476 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.4409621E+01 (-0.2355838E+00) number of electron 87.0000018 magnetization 5.0295460 augmentation part 3.5290087 magnetization 4.0624035 Broyden mixing: rms(total) = 0.14012E+00 rms(broyden)= 0.13997E+00 rms(prec ) = 0.14520E+00 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3797 3.0668 2.4887 1.2444 0.7759 0.6446 0.7186 0.7186 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8712.57049810 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 278.68256574 PAW double counting = 7428.48269624 -7408.25120008 entropy T*S EENTRO = 0.03000883 eigenvalues EBANDS = -1026.69900311 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -188.35559731 eV energy without entropy = -188.38560614 energy(sigma->0) = -188.36560025 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 416 total energy-change (2. order) :-0.9647231E+00 (-0.2161758E-01) number of electron 87.0000018 magnetization 3.4608893 augmentation part 3.5269822 magnetization 2.4811816 Broyden mixing: rms(total) = 0.92805E-01 rms(broyden)= 0.92774E-01 rms(prec ) = 0.96919E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4706 4.2319 2.3024 1.3286 0.8543 0.8543 0.8421 0.6191 0.7322 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.12711171 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 277.84633692 PAW double counting = 7407.23751576 -7387.00745490 entropy T*S EENTRO = 0.03673006 eigenvalues EBANDS = -1025.27616974 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -189.32032044 eV energy without entropy = -189.35705049 energy(sigma->0) = -189.33256379 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 368 total energy-change (2. order) :-0.1019778E+01 (-0.2256579E-01) number of electron 87.0000018 magnetization 2.7962745 augmentation part 3.5313106 magnetization 1.8489824 Broyden mixing: rms(total) = 0.75887E-01 rms(broyden)= 0.75738E-01 rms(prec ) = 0.81162E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4515 4.6075 2.1162 1.7426 0.9657 0.9657 0.6984 0.6177 0.6749 0.6749 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.21914585 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.79919672 PAW double counting = 7367.45611003 -7347.22045151 entropy T*S EENTRO = 0.03499453 eigenvalues EBANDS = -1025.16063563 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -190.34009853 eV energy without entropy = -190.37509307 energy(sigma->0) = -190.35176338 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 392 total energy-change (2. order) :-0.3693101E+00 (-0.5650654E-02) number of electron 87.0000018 magnetization 1.9955695 augmentation part 3.5309753 magnetization 1.0639133 Broyden mixing: rms(total) = 0.44851E-01 rms(broyden)= 0.44768E-01 rms(prec ) = 0.47378E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5001 5.3087 2.0013 2.0013 1.1456 0.9379 0.9379 0.6877 0.6196 0.6806 0.6806 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8715.06657977 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.45885915 PAW double counting = 7363.65836181 -7343.43627432 entropy T*S EENTRO = 0.03714643 eigenvalues EBANDS = -1024.33075509 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -190.70940863 eV energy without entropy = -190.74655506 energy(sigma->0) = -190.72179077 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 352 total energy-change (2. order) :-0.5776746E+00 (-0.3912704E-02) number of electron 87.0000018 magnetization 1.4943190 augmentation part 3.5283483 magnetization 0.5710118 Broyden mixing: rms(total) = 0.29798E-01 rms(broyden)= 0.29777E-01 rms(prec ) = 0.31748E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5571 5.7792 2.2982 2.2982 1.3006 1.1593 0.6691 0.6691 0.8318 0.8318 0.6693 0.6217 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.85033146 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.83169398 PAW double counting = 7345.26286663 -7325.04474175 entropy T*S EENTRO = 0.03702468 eigenvalues EBANDS = -1024.49342851 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.28708325 eV energy without entropy = -191.32410793 energy(sigma->0) = -191.29942481 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.3897379E+00 (-0.1604829E-02) number of electron 87.0000018 magnetization 1.2399945 augmentation part 3.5256246 magnetization 0.3297884 Broyden mixing: rms(total) = 0.17407E-01 rms(broyden)= 0.17355E-01 rms(prec ) = 0.18912E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5762 6.2488 2.5087 2.5087 1.3670 1.0188 1.0188 0.6490 0.6490 0.8784 0.7979 0.6371 0.6317 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.99392254 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.46025723 PAW double counting = 7353.69635602 -7333.48243417 entropy T*S EENTRO = 0.03707226 eigenvalues EBANDS = -1024.36398307 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.67682111 eV energy without entropy = -191.71389337 energy(sigma->0) = -191.68917853 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.2225436E+00 (-0.4689603E-03) number of electron 87.0000018 magnetization 1.0774815 augmentation part 3.5271033 magnetization 0.1760520 Broyden mixing: rms(total) = 0.98609E-02 rms(broyden)= 0.98417E-02 rms(prec ) = 0.10791E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6181 6.6388 2.9384 2.3678 1.3689 1.3689 1.3124 0.6545 0.6545 0.8683 0.8683 0.7272 0.6480 0.6194 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.53377645 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.20979273 PAW double counting = 7354.35474518 -7334.13749256 entropy T*S EENTRO = 0.03720447 eigenvalues EBANDS = -1024.79967123 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.89936470 eV energy without entropy = -191.93656917 energy(sigma->0) = -191.91176619 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.1560780E+00 (-0.2864080E-03) number of electron 87.0000018 magnetization 1.0288788 augmentation part 3.5284777 magnetization 0.1337717 Broyden mixing: rms(total) = 0.64302E-02 rms(broyden)= 0.64140E-02 rms(prec ) = 0.71965E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6216 6.8345 3.4067 2.2160 1.8205 1.5682 0.9625 0.9625 0.6493 0.6493 0.8325 0.8325 0.6895 0.6590 0.6188 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.02479081 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.02941801 PAW double counting = 7358.07104127 -7337.84940500 entropy T*S EENTRO = 0.03723372 eigenvalues EBANDS = -1025.28877307 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.05544272 eV energy without entropy = -192.09267644 energy(sigma->0) = -192.06785396 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.6907190E-01 (-0.7158119E-04) number of electron 87.0000018 magnetization 1.0114346 augmentation part 3.5275810 magnetization 0.1189719 Broyden mixing: rms(total) = 0.31962E-02 rms(broyden)= 0.31841E-02 rms(prec ) = 0.38145E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6233 6.9564 3.6069 2.1466 2.1466 1.5433 1.0515 1.0515 0.6511 0.6511 0.9098 0.9098 0.7734 0.7037 0.6244 0.6244 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.00944185 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.96730543 PAW double counting = 7363.28233235 -7343.06106206 entropy T*S EENTRO = 0.03720997 eigenvalues EBANDS = -1025.31069162 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.12451462 eV energy without entropy = -192.16172459 energy(sigma->0) = -192.13691794 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.1449298E-01 (-0.2732341E-04) number of electron 87.0000018 magnetization 1.0029109 augmentation part 3.5267442 magnetization 0.1118625 Broyden mixing: rms(total) = 0.20569E-02 rms(broyden)= 0.20470E-02 rms(prec ) = 0.25568E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6033 7.0311 3.7413 2.2469 2.2469 1.5029 1.0228 1.0228 1.0063 1.0063 0.6525 0.6525 0.7840 0.7840 0.7025 0.6247 0.6247 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.97484557 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.95442614 PAW double counting = 7364.03406176 -7343.81331689 entropy T*S EENTRO = 0.03718523 eigenvalues EBANDS = -1025.34635143 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13900760 eV energy without entropy = -192.17619283 energy(sigma->0) = -192.15140267 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.4029031E-02 ( 0.3114787E-06) number of electron 87.0000018 magnetization 0.9969254 augmentation part 3.5268528 magnetization 0.1071248 Broyden mixing: rms(total) = 0.13171E-02 rms(broyden)= 0.13111E-02 rms(prec ) = 0.17444E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6139 7.0790 3.8733 2.3280 2.3280 1.7293 1.1918 1.1918 1.0191 1.0191 0.6519 0.6519 0.8522 0.8522 0.7552 0.6740 0.6194 0.6194 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.88191817 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.94695120 PAW double counting = 7363.39845170 -7343.17740142 entropy T*S EENTRO = 0.03719989 eigenvalues EBANDS = -1025.43615300 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14303663 eV energy without entropy = -192.18023652 energy(sigma->0) = -192.15543659 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 360 total energy-change (2. order) :-0.1748961E-02 ( 0.8475948E-06) number of electron 87.0000018 magnetization 0.9953085 augmentation part 3.5270978 magnetization 0.1066891 Broyden mixing: rms(total) = 0.65431E-03 rms(broyden)= 0.64776E-03 rms(prec ) = 0.10147E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6005 7.1110 3.9947 2.3995 2.3995 1.8107 1.2940 1.2940 1.0280 1.0280 0.6517 0.6517 0.8622 0.8622 0.8095 0.7193 0.6501 0.6162 0.6268 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.79477597 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.94232030 PAW double counting = 7362.96681795 -7342.74566622 entropy T*S EENTRO = 0.03721603 eigenvalues EBANDS = -1025.52053085 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14478559 eV energy without entropy = -192.18200162 energy(sigma->0) = -192.15719093 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 320 total energy-change (2. order) :-0.3034105E-03 ( 0.8369209E-05) number of electron 87.0000018 magnetization 0.9942566 augmentation part 3.5271021 magnetization 0.1063037 Broyden mixing: rms(total) = 0.48235E-03 rms(broyden)= 0.47957E-03 rms(prec ) = 0.77593E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6002 7.1359 4.0816 2.7268 2.1311 2.0102 1.4111 1.4111 1.0427 1.0427 0.9406 0.9406 0.6519 0.6519 0.8591 0.7363 0.7363 0.6625 0.6153 0.6153 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.76434658 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.94127430 PAW double counting = 7362.69151615 -7342.47058184 entropy T*S EENTRO = 0.03721316 eigenvalues EBANDS = -1025.54999736 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14508900 eV energy without entropy = -192.18230216 energy(sigma->0) = -192.15749339 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 312 total energy-change (2. order) :-0.4738988E-03 ( 0.7681346E-05) number of electron 87.0000018 magnetization 0.9939321 augmentation part 3.5270968 magnetization 0.1069698 Broyden mixing: rms(total) = 0.31382E-03 rms(broyden)= 0.30883E-03 rms(prec ) = 0.49031E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6368 7.1641 4.1944 3.3367 2.3092 2.3092 1.4414 1.4414 1.0351 1.0351 1.1266 0.6519 0.6519 0.9817 0.8545 0.8545 0.7915 0.6880 0.6349 0.6170 0.6170 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.72395607 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93961739 PAW double counting = 7362.25308129 -7342.03238102 entropy T*S EENTRO = 0.03721335 eigenvalues EBANDS = -1025.58897101 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14556290 eV energy without entropy = -192.18277624 energy(sigma->0) = -192.15796735 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 312 total energy-change (2. order) :-0.4454035E-03 ( 0.9040035E-05) number of electron 87.0000018 magnetization 0.9936744 augmentation part 3.5271008 magnetization 0.1075160 Broyden mixing: rms(total) = 0.23963E-03 rms(broyden)= 0.23456E-03 rms(prec ) = 0.33092E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6488 7.1964 4.4329 3.7549 2.3789 2.3789 1.5091 1.5091 1.0385 1.0385 0.6518 0.6518 1.0921 1.0921 0.8610 0.8610 0.9132 0.6938 0.6938 0.6402 0.6130 0.6236 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.70079056 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93850664 PAW double counting = 7362.13452441 -7341.91389255 entropy T*S EENTRO = 0.03721950 eigenvalues EBANDS = -1025.61140891 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14600830 eV energy without entropy = -192.18322780 energy(sigma->0) = -192.15841480 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 280 total energy-change (2. order) :-0.2651648E-03 ( 0.9997579E-05) number of electron 87.0000018 magnetization 0.9935871 augmentation part 3.5271085 magnetization 0.1079222 Broyden mixing: rms(total) = 0.16700E-03 rms(broyden)= 0.16357E-03 rms(prec ) = 0.22545E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.7112 7.2865 5.4107 3.9824 2.4767 2.4767 1.6901 1.6901 1.2495 1.2495 1.0266 1.0266 0.6519 0.6519 0.9077 0.9077 0.8299 0.8299 0.7699 0.6759 0.6227 0.6227 0.6100 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68898937 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93753636 PAW double counting = 7362.15486358 -7341.93417972 entropy T*S EENTRO = 0.03722087 eigenvalues EBANDS = -1025.62255836 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14627347 eV energy without entropy = -192.18349434 energy(sigma->0) = -192.15868042 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 272 total energy-change (2. order) :-0.2166975E-03 ( 0.1006922E-04) number of electron 87.0000018 magnetization 0.9934970 augmentation part 3.5271027 magnetization 0.1084033 Broyden mixing: rms(total) = 0.11775E-03 rms(broyden)= 0.11248E-03 rms(prec ) = 0.14328E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.7697 7.7336 6.5778 4.0726 2.5454 2.5454 1.8012 1.8012 1.3033 1.3033 1.0281 1.0281 0.6519 0.6519 0.9773 0.9773 0.8481 0.8481 0.7808 0.7161 0.6516 0.6100 0.6250 0.6250 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68505008 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93693556 PAW double counting = 7362.31439921 -7342.09367089 entropy T*S EENTRO = 0.03722097 eigenvalues EBANDS = -1025.62615810 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14649016 eV energy without entropy = -192.18371113 energy(sigma->0) = -192.15889715 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 296 total energy-change (2. order) :-0.2931150E-04 ( 0.1036219E-04) number of electron 87.0000018 magnetization 0.9933572 augmentation part 3.5271007 magnetization 0.1084474 Broyden mixing: rms(total) = 0.71195E-04 rms(broyden)= 0.68868E-04 rms(prec ) = 0.93762E-04 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.7893 8.2918 6.8229 4.1263 2.5758 2.5758 1.8543 1.8543 1.4130 1.4130 1.0320 1.0320 0.6519 0.6519 1.0410 1.0410 0.8597 0.8597 0.8791 0.7283 0.7283 0.6617 0.6212 0.6212 0.6062 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68715444 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93687511 PAW double counting = 7362.36148643 -7342.14074788 entropy T*S EENTRO = 0.03722144 eigenvalues EBANDS = -1025.62403331 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14651948 eV energy without entropy = -192.18374091 energy(sigma->0) = -192.15892662 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 304 total energy-change (2. order) :-0.1131975E-04 ( 0.1041875E-04) number of electron 87.0000018 magnetization 0.9932753 augmentation part 3.5270970 magnetization 0.1085231 Broyden mixing: rms(total) = 0.49990E-04 rms(broyden)= 0.48304E-04 rms(prec ) = 0.65273E-04 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.7898 8.5503 6.8806 4.2030 2.6101 2.6101 1.9734 1.9734 1.4304 1.4304 1.2798 1.0308 1.0308 1.1047 0.6519 0.6519 0.8771 0.8771 0.8335 0.7874 0.7874 0.6753 0.6400 0.6232 0.6232 0.6077 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69258998 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93692283 PAW double counting = 7362.36557772 -7342.14484709 entropy T*S EENTRO = 0.03722144 eigenvalues EBANDS = -1025.61864890 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14653080 eV energy without entropy = -192.18375224 energy(sigma->0) = -192.15893794 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 280 total energy-change (2. order) :-0.6224837E-05 ( 0.1047095E-04) number of electron 87.0000018 magnetization 0.9932004 augmentation part 3.5270968 magnetization 0.1085520 Broyden mixing: rms(total) = 0.47753E-04 rms(broyden)= 0.47055E-04 rms(prec ) = 0.57552E-04 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.7548 8.7791 6.9109 4.2065 2.6502 2.6502 1.9977 1.9977 1.3647 1.3647 1.0301 1.0301 1.2231 1.2231 0.6519 0.6519 0.8936 0.8936 0.7740 0.7740 0.7978 0.6540 0.6540 0.6507 0.6079 0.6079 0.5850 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69645656 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93696411 PAW double counting = 7362.36683997 -7342.14611239 entropy T*S EENTRO = 0.03722137 eigenvalues EBANDS = -1025.61482669 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14653702 eV energy without entropy = -192.18375839 energy(sigma->0) = -192.15894414 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 336 total energy-change (2. order) :-0.1254227E-03 (-0.1090784E-03) number of electron 87.0000018 magnetization 0.9931993 augmentation part 3.5270622 magnetization 0.1086472 Broyden mixing: rms(total) = 0.99727E-03 rms(broyden)= 0.99650E-03 rms(prec ) = 0.10005E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6897 8.7785 6.9108 4.2066 2.6505 2.6505 1.9979 1.9979 1.3651 1.3651 1.0301 1.0301 1.2223 1.2223 0.6519 0.6519 0.8936 0.8936 0.7740 0.7740 0.7977 0.6540 0.6540 0.6505 0.6079 0.6079 0.5843 0.0001 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69863339 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93699979 PAW double counting = 7362.37573713 -7342.15501113 entropy T*S EENTRO = 0.03722157 eigenvalues EBANDS = -1025.61280958 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14666244 eV energy without entropy = -192.18388401 energy(sigma->0) = -192.15906963 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 376 total energy-change (2. order) :-0.8432627E-04 (-0.7629025E-04) number of electron 87.0000018 magnetization 0.9931992 augmentation part 3.5270991 magnetization 0.1086095 Broyden mixing: rms(total) = 0.41399E-03 rms(broyden)= 0.41189E-03 rms(prec ) = 0.41685E-03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.6295 8.7787 6.9109 4.2066 2.6506 2.6506 1.9978 1.9978 1.3653 1.3653 1.0301 1.0301 1.2221 1.2221 0.6519 0.6519 0.8935 0.8935 0.7739 0.7739 0.7978 0.6541 0.6541 0.6504 0.6079 0.6079 0.5842 0.0019 0.0009 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69862997 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93699883 PAW double counting = 7362.37463861 -7342.15391266 entropy T*S EENTRO = 0.03722157 eigenvalues EBANDS = -1025.61289631 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14674677 eV energy without entropy = -192.18396834 energy(sigma->0) = -192.15915396 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 320 total energy-change (2. order) :-0.5024037E-03 (-0.3732143E-03) number of electron 87.0000018 magnetization 0.9931990 augmentation part 3.5271553 magnetization 0.1085529 Broyden mixing: rms(total) = 0.30880E-02 rms(broyden)= 0.30870E-02 rms(prec ) = 0.30891E-02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5734 8.7786 6.9109 4.2065 2.6513 2.6513 1.9984 1.9984 1.3665 1.3665 1.0301 1.0301 1.2208 1.2208 0.6519 0.6519 0.8930 0.8930 0.7738 0.7738 0.7983 0.6541 0.6541 0.6509 0.6068 0.6068 0.5860 0.0013 0.0013 0.0016 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69862510 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93699898 PAW double counting = 7362.37481424 -7342.15408807 entropy T*S EENTRO = 0.03722156 eigenvalues EBANDS = -1025.61340395 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.14724917 eV energy without entropy = -192.18447074 energy(sigma->0) = -192.15965636 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) :-0.1614545E+03 (-0.1613881E+03) number of electron 86.9999784 magnetization 0.9937567 augmentation part 3.4350029 magnetization 0.2012661 Broyden mixing: rms(total) = 0.41981E+01 rms(broyden)= 0.41977E+01 rms(prec ) = 0.42030E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.5207 8.7730 6.9102 4.2068 2.6545 2.6545 2.0003 2.0003 1.3695 1.3695 1.0301 1.0301 1.2126 1.2126 0.6519 0.6519 0.8935 0.8935 0.7748 0.7748 0.7984 0.6546 0.6546 0.6501 0.6082 0.6082 0.5794 0.0016 0.0000 0.0003 0.0018 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.69866848 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93700077 PAW double counting = 7362.37541671 -7342.15469034 entropy T*S EENTRO = 0.03721976 eigenvalues EBANDS = -1187.06786278 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -353.60175118 eV energy without entropy = -353.63897094 energy(sigma->0) = -353.61415777 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 400 total energy-change (2. order) : 0.1618303E+03 (-0.4654001E+02) number of electron 86.9999988 magnetization 0.9690160 augmentation part 3.5063433 magnetization 0.1053795 Broyden mixing: rms(total) = 0.87130E-01 rms(broyden)= 0.78335E-01 rms(prec ) = 0.79823E-01 weight for this iteration 100.00 WARNING: grid for Broyden might be to small eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4720 8.7709 6.9097 4.2074 2.6562 2.6562 1.9986 1.9986 1.3710 1.3710 1.0301 1.0301 1.2069 1.2069 0.6519 0.6519 0.8936 0.8936 0.7743 0.7743 0.8016 0.6569 0.6569 0.6496 0.6082 0.6082 0.5814 0.0041 0.0041 0.0046 0.0000 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68425407 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.98460394 PAW double counting = 7363.58653906 -7343.36578269 entropy T*S EENTRO = 0.03721507 eigenvalues EBANDS = -1025.29962811 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -191.77147362 eV energy without entropy = -191.80868869 energy(sigma->0) = -191.78387864 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 336 total energy-change (2. order) :-0.3634844E+00 (-0.2943369E+00) number of electron 87.0000017 magnetization 0.9803836 augmentation part 3.5311432 magnetization 0.1073768 Broyden mixing: rms(total) = 0.41611E-01 rms(broyden)= 0.39722E-01 rms(prec ) = 0.40127E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4289 8.7729 6.9096 4.2085 2.6572 2.6572 1.9912 1.9912 1.3689 1.3689 1.0301 1.0301 1.2096 1.2096 0.6519 0.6519 0.8956 0.8956 0.7747 0.7747 0.8026 0.6527 0.6527 0.6504 0.6083 0.6083 0.5913 0.0374 0.0374 0.0229 0.0089 0.0000 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.12627530 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.90118156 PAW double counting = 7371.16071628 -7350.91864580 entropy T*S EENTRO = 0.03731644 eigenvalues EBANDS = -1026.15908434 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.13495799 eV energy without entropy = -192.17227443 energy(sigma->0) = -192.14739681 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 408 total energy-change (2. order) :-0.1986123E-01 (-0.1481045E-02) number of electron 87.0000018 magnetization 0.9869213 augmentation part 3.5306335 magnetization 0.1087267 Broyden mixing: rms(total) = 0.20779E-01 rms(broyden)= 0.20587E-01 rms(prec ) = 0.21236E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.4044 8.7409 6.9053 4.2137 2.6720 2.6720 1.9908 1.9908 1.3582 1.3582 1.2347 1.2347 1.0294 1.0294 0.8959 0.8959 0.6519 0.6519 0.7780 0.7780 0.7932 0.6612 0.6612 0.6558 0.6147 0.5970 0.5854 0.3086 0.3086 0.0277 0.0399 0.0091 0.0000 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.35833869 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.90936224 PAW double counting = 7367.34460749 -7347.10732314 entropy T*S EENTRO = 0.03730418 eigenvalues EBANDS = -1025.95026447 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.15481922 eV energy without entropy = -192.19212340 energy(sigma->0) = -192.16725395 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 384 total energy-change (2. order) :-0.1667958E-04 ( 0.6153707E-02) number of electron 87.0000018 magnetization 0.9911866 augmentation part 3.5289396 magnetization 0.1099392 Broyden mixing: rms(total) = 0.15118E-01 rms(broyden)= 0.15100E-01 rms(prec ) = 0.15470E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3990 8.6253 6.8888 4.2221 2.7038 2.7038 2.0143 2.0143 0.7528 0.7528 1.3727 1.3727 1.2019 1.2019 1.0295 1.0295 0.6519 0.6519 0.8934 0.8934 0.7844 0.7844 0.7854 0.6628 0.6628 0.6457 0.6099 0.6023 0.6023 0.3709 0.0278 0.0405 0.0091 0.0000 0.0015 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.64694942 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.92616737 PAW double counting = 7368.20524461 -7347.97206001 entropy T*S EENTRO = 0.03727480 eigenvalues EBANDS = -1025.67434643 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.15483590 eV energy without entropy = -192.19211070 energy(sigma->0) = -192.16726084 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 472 total energy-change (2. order) :-0.2127618E-01 (-0.1037690E-01) number of electron 87.0000018 magnetization 0.9909936 augmentation part 3.5291040 magnetization 0.1065530 Broyden mixing: rms(total) = 0.15192E-01 rms(broyden)= 0.15184E-01 rms(prec ) = 0.15385E-01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3596 8.6254 6.8888 4.2222 2.7030 2.7030 2.0138 2.0138 1.3724 1.3724 0.7489 0.7489 1.2033 1.2033 1.0294 1.0294 0.8927 0.8927 0.6519 0.6519 0.7839 0.7839 0.7861 0.6635 0.6635 0.6456 0.6102 0.6027 0.6027 0.3755 0.0404 0.0279 0.0229 0.0091 0.0015 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.67870770 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93037944 PAW double counting = 7369.22457249 -7348.99337343 entropy T*S EENTRO = 0.03726278 eigenvalues EBANDS = -1025.66607884 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -192.17611208 eV energy without entropy = -192.21337487 energy(sigma->0) = -192.18853301 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 448 total energy-change (2. order) :-0.3864106E+04 (-0.3864244E+04) number of electron 87.0036542 magnetization 0.9907170 augmentation part 2.9401522 magnetization 1.5440576 Broyden mixing: rms(total) = 0.56202E+01 rms(broyden)= 0.56193E+01 rms(prec ) = 0.56572E+01 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.3221 8.6243 6.8883 4.2230 2.7069 2.7069 2.0147 2.0147 1.3731 1.3731 0.7433 0.7433 1.1992 1.1992 1.0294 1.0294 0.8939 0.8939 0.6519 0.6519 0.7839 0.7839 0.7865 0.6694 0.6694 0.6423 0.6058 0.6015 0.6015 0.3779 0.0404 0.0332 0.0277 0.0091 0.0015 0.0002 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68626888 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 274.93079224 PAW double counting = 7369.45843740 -7349.22747557 entropy T*S EENTRO = -0.02331947 eigenvalues EBANDS = -4889.70365020 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -4056.28165132 eV energy without entropy = -4056.25833185 energy(sigma->0) = -4056.27387816 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 464 total energy-change (2. order) :-0.1709643E+06 (-0.1727702E+06) number of electron 86.7833307 magnetization 0.9821187 augmentation part 1.2807382 magnetization 16.6048493 Broyden mixing: rms(total) = 0.25099E+02 rms(broyden)= 0.25098E+02 rms(prec ) = 0.25246E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2863 8.6209 6.8880 4.2229 2.7039 2.7039 2.0128 2.0128 1.3723 1.3723 0.7506 0.7506 1.2041 1.2041 1.0293 1.0293 0.8928 0.8928 0.6519 0.6519 0.7825 0.7825 0.7869 0.6674 0.6674 0.6423 0.6097 0.6011 0.6011 0.3799 0.0404 0.0279 0.0270 0.0091 0.0015 0.0003 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.68143524 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.06964603 PAW double counting = 7369.74306861 -7349.51202412 entropy T*S EENTRO = -0.00748828 eigenvalues EBANDS = -175854.13579061 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -175020.55419044 eV energy without entropy = -175020.54670216 energy(sigma->0) = -175020.55169435 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 424 total energy-change (2. order) : 0.1456781E+06 (-0.9420046E+04) number of electron 87.4832624 magnetization 0.9821335 augmentation part 1.1361413 magnetization 20.7254041 Broyden mixing: rms(total) = 0.29368E+02 rms(broyden)= 0.29368E+02 rms(prec ) = 0.29500E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2535 8.6213 6.8880 4.2241 2.7041 2.7041 2.0155 2.0155 1.3721 1.3721 0.7692 0.7692 1.2012 1.2012 1.0293 1.0293 0.6519 0.6519 0.8911 0.8911 0.7821 0.7821 0.7903 0.6699 0.6699 0.6456 0.6096 0.6019 0.6019 0.3841 0.0404 0.0278 0.0147 0.0091 0.0015 0.0013 0.0002 0.0001 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.63903669 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.43941857 PAW double counting = 7373.22850136 -7352.99572085 entropy T*S EENTRO = -0.01869277 eigenvalues EBANDS = -30176.39475404 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -29342.41045125 eV energy without entropy = -29342.39175848 energy(sigma->0) = -29342.40422033 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 488 total energy-change (2. order) :-0.1172297E+06 (-0.1378586E+06) number of electron 87.9426667 magnetization 0.9821919 augmentation part 0.6866639 magnetization 39.1682203 Broyden mixing: rms(total) = 0.58636E+02 rms(broyden)= 0.58636E+02 rms(prec ) = 0.58766E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.2215 8.6203 6.8878 4.2241 2.7046 2.7046 2.0153 2.0153 0.7683 0.7683 1.3724 1.3724 1.2008 1.2008 1.0292 1.0292 0.6519 0.6519 0.8912 0.8912 0.7821 0.7821 0.7897 0.6699 0.6699 0.6455 0.6100 0.6016 0.6016 0.3852 0.0404 0.0279 0.0195 0.0091 0.0015 0.0010 0.0002 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.63676353 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.81867279 PAW double counting = 7373.22017569 -7352.98741198 entropy T*S EENTRO = -0.03935965 eigenvalues EBANDS = -147406.42834734 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -146572.08320085 eV energy without entropy = -146572.04384120 energy(sigma->0) = -146572.07008097 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 504 total energy-change (2. order) : 0.7057535E+05 (-0.7676665E+05) number of electron 86.8595051 magnetization 0.9961592 augmentation part 0.0159575 magnetization 44.9693961 Broyden mixing: rms(total) = 0.10523E+03 rms(broyden)= 0.10523E+03 rms(prec ) = 0.10532E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1914 8.6177 6.8877 4.2241 2.7051 2.7051 2.0166 2.0166 1.3725 1.3725 0.7663 0.7663 1.2005 1.2005 1.0291 1.0291 0.8914 0.8914 0.6519 0.6519 0.7819 0.7819 0.7885 0.6694 0.6694 0.6454 0.6054 0.6063 0.6063 0.3876 0.0404 0.0368 0.0278 0.0091 0.0015 0.0005 0.0002 0.0001 0.0000 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8713.83188252 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.22582282 PAW double counting = 7376.08595685 -7355.85325957 entropy T*S EENTRO = -0.01479596 eigenvalues EBANDS = -76831.31068076 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -75996.72900597 eV energy without entropy = -75996.71421002 energy(sigma->0) = -75996.72407399 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 488 total energy-change (2. order) : 0.5260737E+05 (-0.2367499E+05) number of electron 88.0472656 magnetization 1.0160605 augmentation part 0.4328842 magnetization 40.9174052 Broyden mixing: rms(total) = 0.85037E+02 rms(broyden)= 0.85037E+02 rms(prec ) = 0.85125E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1625 8.6212 6.8881 4.2242 2.7043 2.7043 2.0153 2.0153 1.3714 1.3714 0.7662 0.7662 1.2029 1.2029 1.0291 1.0291 0.8910 0.8910 0.6519 0.6519 0.7816 0.7816 0.7891 0.6707 0.6707 0.6453 0.6058 0.6064 0.6064 0.3874 0.0404 0.0390 0.0278 0.0091 0.0015 0.0010 0.0002 0.0001 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8714.79140981 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 275.82588540 PAW double counting = 7373.61022907 -7353.37917941 entropy T*S EENTRO = 0.02241291 eigenvalues EBANDS = -24222.62139820 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -23389.36362686 eV energy without entropy = -23389.38603978 energy(sigma->0) = -23389.37109783 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 448 total energy-change (2. order) :-0.3950953E+05 (-0.5294870E+05) number of electron 88.0978541 magnetization 1.0344334 augmentation part 0.4580736 magnetization 28.9375866 Broyden mixing: rms(total) = 0.11162E+03 rms(broyden)= 0.11162E+03 rms(prec ) = 0.11168E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1349 8.6224 6.8881 4.2244 2.7051 2.7051 2.0176 2.0176 1.3726 1.3726 0.7649 0.7649 1.1996 1.1996 1.0291 1.0291 0.8917 0.8917 0.6519 0.6519 0.7822 0.7822 0.7887 0.6707 0.6707 0.6452 0.6055 0.6064 0.6064 0.3873 0.0398 0.0404 0.0278 0.0091 0.0015 0.0014 0.0004 0.0004 0.0001 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8716.26563035 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.44062079 PAW double counting = 7381.28502032 -7361.06462294 entropy T*S EENTRO = -0.03969946 eigenvalues EBANDS = -63731.22210420 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -62898.89658267 eV energy without entropy = -62898.85688322 energy(sigma->0) = -62898.88334952 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 480 total energy-change (2. order) :-0.4083019E+06 (-0.4056850E+06) number of electron 87.6503927 magnetization 1.0343570 augmentation part 0.1629173 magnetization 44.1292183 Broyden mixing: rms(total) = 0.17400E+03 rms(broyden)= 0.17400E+03 rms(prec ) = 0.17405E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.1094 8.6248 6.8886 4.2250 2.7078 2.7078 2.0179 2.0179 1.3708 1.3708 0.7882 0.7882 1.2010 1.2010 1.0294 1.0294 0.8915 0.8915 0.6519 0.6519 0.7842 0.7842 0.7851 0.6730 0.6730 0.6476 0.6091 0.6091 0.6020 0.3860 0.0404 0.0278 0.0094 0.0091 0.0022 0.0017 0.0017 0.0015 0.0003 0.0001 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8716.52239533 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.40531393 PAW double counting = 7382.40494839 -7362.18349389 entropy T*S EENTRO = -0.01526380 eigenvalues EBANDS = -472032.86770109 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -471200.80875863 eV energy without entropy = -471200.79349482 energy(sigma->0) = -471200.80367069 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 464 total energy-change (2. order) :-0.2304218E+06 (-0.6848004E+06) number of electron 87.5062142 magnetization 1.0452208 augmentation part 0.3226153 magnetization 37.2963541 Broyden mixing: rms(total) = 0.21346E+03 rms(broyden)= 0.21346E+03 rms(prec ) = 0.21350E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0840 8.6211 6.8882 4.2250 2.7081 2.7081 2.0179 2.0179 1.3703 1.3703 0.7846 0.7846 1.2013 1.2013 1.0294 1.0294 0.8915 0.8915 0.6519 0.6519 0.7849 0.7849 0.7849 0.6732 0.6732 0.6479 0.6089 0.6089 0.6017 0.3845 0.0404 0.0278 0.0127 0.0091 0.0020 0.0016 0.0016 0.0015 0.0003 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8716.42522896 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.31574521 PAW double counting = 7380.77623936 -7360.55449980 entropy T*S EENTRO = -0.03229142 eigenvalues EBANDS = -702454.66054345 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -701622.61074588 eV energy without entropy = -701622.57845447 energy(sigma->0) = -701622.59998208 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 496 total energy-change (2. order) : 0.5746308E+06 (-0.1217963E+06) number of electron 87.9159461 magnetization 1.0510031 augmentation part 0.5551103 magnetization 42.5057271 Broyden mixing: rms(total) = 0.19526E+03 rms(broyden)= 0.19526E+03 rms(prec ) = 0.19530E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0432 8.7159 6.8713 4.1519 2.8607 2.4868 0.8031 0.8031 1.5027 1.3694 1.3694 1.2886 1.0361 1.0361 0.9311 0.6561 0.6561 0.7727 0.7727 0.3865 0.6742 0.6742 0.6210 0.6131 0.5962 0.0347 0.0126 0.0110 0.0054 0.0054 0.0054 0.0015 0.0008 0.0003 0.0001 0.0001 0.0003 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8716.69938991 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.27845366 PAW double counting = 7380.86706619 -7360.64556795 entropy T*S EENTRO = 0.00936630 eigenvalues EBANDS = -127823.59897779 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -126991.81921632 eV energy without entropy = -126991.82858262 energy(sigma->0) = -126991.82233842 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 496 total energy-change (2. order) :-0.6051996E+05 (-0.1719922E+06) number of electron 86.4084993 magnetization 1.3882697 augmentation part -1.6354901 magnetization 60.6869374 Broyden mixing: rms(total) = 0.55416E+03 rms(broyden)= 0.55416E+03 rms(prec ) = 0.55418E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0176 8.7044 6.8685 4.1521 2.8585 2.4877 0.8047 0.8047 1.5082 1.3687 1.3687 1.2925 1.0362 1.0362 0.9291 0.6561 0.6561 0.7718 0.7718 0.3874 0.6733 0.6733 0.6161 0.6161 0.5975 0.0337 0.0149 0.0130 0.0050 0.0050 0.0026 0.0026 0.0015 0.0006 0.0004 0.0004 0.0003 0.0001 0.0001 0.0001 0.0000 0.0000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5741.54159136 -Hartree energ DENC = -8717.03881338 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 276.62471716 PAW double counting = 7376.03309890 -7355.81501390 entropy T*S EENTRO = -0.04552007 eigenvalues EBANDS = -188343.50413839 atomic energy EATOM = 3504.07976474 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -187511.77583652 eV energy without entropy = -187511.73031645 energy(sigma->0) = -187511.76066316 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 9 2 16 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------