vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.09.02 08:32:07 running on 4 total cores distrk: each k-point on 4 cores, 1 groups distr: one band on NCORE= 1 cores, 4 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 200 ALGO = Normal (blocked Davidson) ISPIN = 2 MAGMOM = 34*0 0.107 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 1.11 0.73 0.32 NPAR = 4 POTCAR: PAW_PBE Si 05Jan2001 POTCAR: PAW_PBE O 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE Si 05Jan2001 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE 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 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 Optimization of the real space projectors (new method) maximal supplied QI-value = 19.84 optimisation between [QCUT,QGAM] = [ 10.12, 20.44] = [ 28.68,116.96] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 7 10.119 159.560 0.56E-04 0.22E-03 0.45E-07 0 7 10.119 115.863 0.56E-04 0.21E-03 0.45E-07 1 7 10.119 88.339 0.34E-03 0.49E-03 0.11E-06 1 7 10.119 48.592 0.33E-03 0.48E-03 0.11E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 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 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 PAW_PBE Si 05Jan2001 : energy of atom 1 EATOM= -103.0669 kinetic energy error for atom= 0.0012 (will be added to EATOM!!) PAW_PBE O 08Apr2002 : energy of atom 2 EATOM= -432.3788 kinetic energy error for atom= 0.1156 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 3 EATOM= -12.4884 kinetic energy error for atom= 0.0098 (will be added to EATOM!!) POSCAR: No title positions in direct lattice No initial velocities read in exchange correlation table for LEXCH = 8 RHO(1)= 0.500 N(1) = 2000 RHO(2)= 100.500 N(2) = 4000 -------------------------------------------------------------------------------------------------------- ion position nearest neighbor table 1 0.842 0.974 0.160- 24 1.66 2 0.013 0.422 0.520- 30 2.63 3 0.349 0.986 0.774- 22 0.20 20 0.77 4 0.144 0.357 0.155- 8 0.78 11 2.21 5 0.241 0.231 0.895- 6 0.303 0.341 0.521- 7 0.947 0.211 0.723- 32 2.50 23 2.57 8 0.144 0.328 0.204- 4 0.78 11 1.95 9 0.371 0.642 0.704- 25 2.18 21 2.53 10 0.696 0.200 0.499- 32 0.84 23 1.74 11 0.375 0.344 0.275- 8 1.95 4 2.21 12 0.822 0.085 0.971- 26 2.33 13 0.152 0.793 0.053- 14 0.594 0.162 0.152- 15 0.682 0.824 0.347- 16 0.405 0.864 0.554- 35 0.75 17 0.631 0.687 0.534- 33 0.55 21 1.99 18 0.177 0.175 0.398- 29 1.54 19 0.485 0.619 0.961- 27 2.15 20 0.286 1.000 0.823- 22 0.75 3 0.77 21 0.390 0.724 0.526- 17 1.99 33 2.00 9 2.53 22 0.358 0.981 0.790- 3 0.20 20 0.75 23 0.772 0.277 0.561- 10 1.74 32 1.98 7 2.57 24 0.880 0.931 0.288- 1 1.66 25 0.133 0.602 0.788- 34 0.61 9 2.18 30 2.64 26 0.766 0.199 0.935- 12 2.33 27 0.418 0.529 0.859- 19 2.15 28 0.669 0.396 0.076- 29 0.346 0.195 0.332- 18 1.54 30 0.069 0.551 0.568- 2 2.63 25 2.64 31 0.054 0.607 0.284- 32 0.736 0.177 0.558- 10 0.84 23 1.98 7 2.50 33 0.598 0.665 0.556- 17 0.55 21 2.00 34 0.183 0.593 0.748- 25 0.61 35 0.317 0.880 0.545- 16 0.75 ----------------------------------------------------------------------------- | | | W W AA RRRRR N N II N N GGGG !!! | | W W A A R R NN N II NN N G G !!! | | W W A A R R N N N II N N N G !!! | | W WW W AAAAAA RRRRR N N N II N N N G GGG ! | | WW WW A A R R N NN II N NN G G | | W W A A R R N N II N N GGGG !!! | | | | The distance between some ions is very small. Please check the | | nearest-neighbor list in the OUTCAR file. | | I HOPE YOU KNOW WHAT YOU ARE DOING! | | | ----------------------------------------------------------------------------- LATTYP: Found a simple orthorhombic cell. ALAT = 7.6631000000 B/A-ratio = 1.4142135689 C/A-ratio = 2.5659458966 Lattice vectors: A1 = ( -7.6631000000, 0.0000000000, 0.0000000000) A2 = ( 0.0000000000, 0.0000000000, 10.8372600000) A3 = ( 0.0000000000, 19.6631000000, 0.0000000000) Analysis of symmetry for initial positions (statically): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The static configuration has the point symmetry C_1 . Analysis of symmetry for dynamics (positions and initial velocities): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The dynamic configuration has the point symmetry C_1 . Analysis of structural, dynamic, and magnetic symmetry: ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple orthorhombic supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 8 trial point group operations. The overall configuration has the point symmetry C_1 . Subroutine INISYM returns: Found 1 space group operations (whereof 1 operations are pure point group operations), and found 1 'primitive' translations ---------------------------------------------------------------------------------------- Primitive cell volume of cell : 1632.9616 direct lattice vectors reciprocal lattice vectors 7.663100000 0.000000000 0.000000000 0.130495491 0.000000000 0.000000000 0.000000000 19.663100000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 0.000000000 10.837260000 0.000000000 0.000000000 0.092274246 length of vectors 7.663100000 19.663100000 10.837260000 0.130495491 0.050856681 0.092274246 position of ions in fractional coordinates (direct lattice) 0.841825130 0.974111680 0.159570970 0.013427760 0.421750780 0.519904320 0.349091350 0.985620810 0.774058930 0.143954380 0.357226140 0.154565710 0.241302320 0.231226110 0.895063410 0.302832430 0.341128360 0.520580960 0.947283590 0.210805090 0.722920360 0.144009080 0.328270890 0.204224410 0.370955720 0.641579330 0.703943070 0.695590830 0.200405580 0.499212460 0.374626220 0.343990600 0.275283650 0.822386460 0.084612720 0.971482380 0.151652130 0.793155080 0.052868040 0.594301650 0.162222600 0.151552940 0.682135440 0.823533660 0.347156080 0.404858500 0.864124940 0.554370040 0.630768370 0.686670620 0.533730620 0.176635390 0.175025900 0.398341280 0.484924640 0.619441790 0.960503850 0.285652420 0.999725960 0.822769630 0.390480980 0.724457720 0.526280540 0.357559300 0.981433480 0.790011710 0.771700840 0.276766240 0.560543850 0.879805640 0.930587550 0.287580260 0.133164060 0.602167970 0.788076140 0.766243550 0.199264970 0.935447850 0.418339890 0.528982780 0.859491030 0.669234540 0.396158180 0.075842850 0.346478890 0.194913020 0.332027590 0.069070070 0.550873960 0.568014030 0.053589530 0.606718150 0.283548220 0.736057120 0.176985110 0.557973400 0.598109930 0.664917690 0.556149690 0.183422110 0.592905470 0.748176010 0.316506690 0.880054060 0.545197590 ion indices of the primitive-cell ions primitive index ion index 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31 32 32 33 33 34 34 35 35 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 2 1 2 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.065247746 0.000000000 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.046137123 0.000000000 0.000000000 0.500000000 Length of vectors 0.065247746 0.050856681 0.046137123 Shift w.r.t. Gamma in fractional coordinates (k-lattice) 0.000000000 0.000000000 0.000000000 Subroutine IBZKPT returns following result: =========================================== Found 4 irreducible k-points: Following reciprocal coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.000000 0.500000 1.000000 0.500000 0.000000 0.500000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.065248 0.000000 0.000000 1.000000 0.000000 0.000000 0.046137 1.000000 0.065248 0.000000 0.046137 1.000000 -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 4 k-points in BZ NKDIM = 4 number of bands NBANDS= 88 number of dos NEDOS = 301 number of ions NIONS = 35 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 207360 max r-space proj IRMAX = 1449 max aug-charges IRDMAX= 4445 dimension x,y,z NGX = 40 NGY = 96 NGZ = 54 dimension x,y,z NGXF= 80 NGYF= 192 NGZF= 108 support grid NGXF= 80 NGYF= 192 NGZF= 108 ions per type = 32 1 2 NGX,Y,Z is equivalent to a cutoff of 8.68, 8.12, 8.28 a.u. NGXF,Y,Z is equivalent to a cutoff of 17.36, 16.23, 16.57 a.u. SYSTEM = No title POSCAR = No title Startparameter for this run: NWRITE = 1 write-flag & timer PREC = normal normal or accurate (medium, high low for compatibility) ISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ISPIN = 2 spin polarized calculation? LNONCOLLINEAR = F non collinear calculations LSORBIT = F spin-orbit coupling INIWAV = 1 electr: 0-lowe 1-rand 2-diag LASPH = F aspherical Exc in radial PAW Electronic Relaxation 1 ENCUT = 400.0 eV 29.40 Ry 5.42 a.u. 12.50 32.07 17.67*2*pi/ulx,y,z ENINI = 400.0 initial cutoff ENAUG = 605.4 eV augmentation charge cutoff NELM = 200; NELMIN= 2; NELMDL= -5 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM LREAL = T real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 1 Vosko Wilk Nusair interpolation ROPT = -0.00050 -0.00050 -0.00050 Ionic relaxation EDIFFG = -.2E-01 stopping-criterion for IOM NSW = 0 number of steps for IOM NBLOCK = 1; KBLOCK = 1 inner block; outer block IBRION = -1 ionic relax: 0-MD 1-quasi-New 2-CG NFREE = 0 steps in history (QN), initial steepest desc. (CG) ISIF = 2 stress and relaxation IWAVPR = 10 prediction: 0-non 1-charg 2-wave 3-comb ISYM = 2 0-nonsym 1-usesym 2-fastsym LCORR = T Harris-Foulkes like correction to forces POTIM = 0.5000 time-step for ionic-motion TEIN = 0.0 initial temperature TEBEG = 0.0; TEEND = 0.0 temperature during run SMASS = -3.00 Nose mass-parameter (am) estimated Nose-frequenzy (Omega) = 0.10E-29 period in steps = 0.13E+47 mass= -0.134E-26a.u. SCALEE = 1.0000 scale energy and forces NPACO = 256; APACO = 16.0 distance and # of slots for P.C. PSTRESS= 0.0 pullay stress Mass of Ions in am POMASS = 28.09 16.00 1.00 Ionic Valenz ZVAL = 4.00 6.00 1.00 Atomic Wigner-Seitz radii RWIGS = 1.11 0.73 0.32 virtual crystal weights VCA = 1.00 1.00 1.00 NELECT = 136.0000 total number of electrons NUPDOWN= -1.0000 fix difference up-down DOS related values: EMIN = 10.00; EMAX =-10.00 energy-range for DOS EFERMI = 0.00 ISMEAR = 1; SIGMA = 0.20 broadening in eV -4-tet -1-fermi 0-gaus Electronic relaxation 2 (details) IALGO = 38 algorithm LDIAG = T sub-space diagonalisation (order eigenvalues) LSUBROT= F optimize rotation matrix (better conditioning) TURBO = 0 0=normal 1=particle mesh IRESTART = 0 0=no restart 2=restart with 2 vectors NREBOOT = 0 no. of reboots NMIN = 0 reboot dimension EREF = 0.00 reference energy to select bands IMIX = 4 mixing-type and parameters AMIX = 0.40; BMIX = 1.00 AMIX_MAG = 1.60; BMIX_MAG = 1.00 AMIN = 0.10 WC = 100.; INIMIX= 1; MIXPRE= 1; MAXMIX= -45 Intra band minimization: WEIMIN = 0.0000 energy-eigenvalue tresh-hold EBREAK = 0.28E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 46.66 314.85 Fermi-wavevector in a.u.,A,eV,Ry = 0.714928 1.351019 6.954243 0.511122 Thomas-Fermi vector in A = 1.802956 Write flags LWAVE = F write WAVECAR LDOWNSAMPLE = F k-point downsampling of WAVECAR LCHARG = F write CHGCAR LVTOT = F write LOCPOT, total local potential LVHAR = F write LOCPOT, Hartree potential only LELF = F write electronic localiz. function (ELF) LORBIT = 0 0 simple, 1 ext, 2 COOP (PROOUT), +10 PAW based schemes Dipole corrections LMONO = F monopole corrections only (constant potential shift) LDIPOL = F correct potential (dipole corrections) IDIPOL = 0 1-x, 2-y, 3-z, 4-all directions EPSILON= 1.0000000 bulk dielectric constant Exchange correlation treatment: GGA = -- GGA type LEXCH = 8 internal setting for exchange type VOSKOWN= 1 Vosko Wilk Nusair interpolation LHFCALC = F Hartree Fock is set to LHFONE = F Hartree Fock one center treatment AEXX = 0.0000 exact exchange contribution Linear response parameters LEPSILON= F determine dielectric tensor LRPA = F only Hartree local field effects (RPA) LNABLA = F use nabla operator in PAW spheres LVEL = F velocity operator in full k-point grid LINTERFAST= F fast interpolation KINTER = 0 interpolate to denser k-point grid CSHIFT =0.1000 complex shift for real part using Kramers Kronig OMEGAMAX= -1.0 maximum frequency DEG_THRESHOLD= 0.2000000E-02 threshold for treating states as degnerate RTIME = -0.100 relaxation time in fs (WPLASMAI= 0.000 imaginary part of plasma frequency in eV, 0.658/RTIME) DFIELD = 0.0000000 0.0000000 0.0000000 field for delta impulse in time Orbital magnetization related: ORBITALMAG= F switch on orbital magnetization LCHIMAG = F perturbation theory with respect to B field DQ = 0.001000 dq finite difference perturbation B field LLRAUG = F two centre corrections for induced B field -------------------------------------------------------------------------------------------------------- Static calculation charge density and potential will be updated during run spin polarized calculation Variant of blocked Davidson Davidson routine will perform the subspace rotation perform sub-space diagonalisation after iterative eigenvector-optimisation modified Broyden-mixing scheme, WC = 100.0 initial mixing is a Kerker type mixing with AMIX = 0.4000 and BMIX = 1.0000 Hartree-type preconditioning will be used using additional bands 20 real space projection scheme for non local part use partial core corrections calculate Harris-corrections to forces (improved forces if not selfconsistent) use gradient corrections use of overlap-Matrix (Vanderbilt PP) Methfessel and Paxton Order N= 1 SIGMA = 0.20 -------------------------------------------------------------------------------------------------------- energy-cutoff : 400.00 volume of cell : 1632.96 direct lattice vectors reciprocal lattice vectors 7.663100000 0.000000000 0.000000000 0.130495491 0.000000000 0.000000000 0.000000000 19.663100000 0.000000000 0.000000000 0.050856681 0.000000000 0.000000000 0.000000000 10.837260000 0.000000000 0.000000000 0.092274246 length of vectors 7.663100000 19.663100000 10.837260000 0.130495491 0.050856681 0.092274246 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 0.250 0.06524775 0.00000000 0.00000000 0.250 0.00000000 0.00000000 0.04613712 0.250 0.06524775 0.00000000 0.04613712 0.250 k-points in reciprocal lattice and weights: Automatic mesh 0.00000000 0.00000000 0.00000000 0.250 0.50000000 0.00000000 0.00000000 0.250 0.00000000 0.00000000 0.50000000 0.250 0.50000000 0.00000000 0.50000000 0.250 position of ions in fractional coordinates (direct lattice) 0.84182513 0.97411168 0.15957097 0.01342776 0.42175078 0.51990432 0.34909135 0.98562081 0.77405893 0.14395438 0.35722614 0.15456571 0.24130232 0.23122611 0.89506341 0.30283243 0.34112836 0.52058096 0.94728359 0.21080509 0.72292036 0.14400908 0.32827089 0.20422441 0.37095572 0.64157933 0.70394307 0.69559083 0.20040558 0.49921246 0.37462622 0.34399060 0.27528365 0.82238646 0.08461272 0.97148238 0.15165213 0.79315508 0.05286804 0.59430165 0.16222260 0.15155294 0.68213544 0.82353366 0.34715608 0.40485850 0.86412494 0.55437004 0.63076837 0.68667062 0.53373062 0.17663539 0.17502590 0.39834128 0.48492464 0.61944179 0.96050385 0.28565242 0.99972596 0.82276963 0.39048098 0.72445772 0.52628054 0.35755930 0.98143348 0.79001171 0.77170084 0.27676624 0.56054385 0.87980564 0.93058755 0.28758026 0.13316406 0.60216797 0.78807614 0.76624355 0.19926497 0.93544785 0.41833989 0.52898278 0.85949103 0.66923454 0.39615818 0.07584285 0.34647889 0.19491302 0.33202759 0.06907007 0.55087396 0.56801403 0.05358953 0.60671815 0.28354822 0.73605712 0.17698511 0.55797340 0.59810993 0.66491769 0.55614969 0.18342211 0.59290547 0.74817601 0.31650669 0.88005406 0.54519759 position of ions in cartesian coordinates (Angst): 6.45099015 19.15405538 1.72931209 0.10289827 8.29292776 5.63433829 2.67512192 19.38036055 8.38867788 1.10313681 7.02417331 1.67506879 1.84912381 4.54662212 9.70003489 2.32063519 6.70764106 5.64167121 7.25912888 4.14508157 7.83447590 1.10355598 6.45482334 2.21323303 2.84267078 12.61543852 7.62881407 5.33038209 3.94059496 5.41009522 2.87079819 6.76392157 2.98332049 6.30202968 1.66374837 10.52820714 1.16212544 15.59588765 0.57294470 4.55419297 3.18979921 1.64241861 5.22727209 16.19322471 3.76222070 3.10247117 16.99137511 6.00785226 4.83364110 13.50207307 5.78417750 1.35357466 3.44155177 4.31692802 3.71602601 12.18014586 10.40922995 2.18898306 19.65771152 8.91656840 2.99229480 14.24508459 5.70343904 2.74001267 19.29802466 8.56156230 5.91362071 5.44208225 6.07475944 6.74203860 18.29823605 3.11658205 1.02044951 11.84048901 8.54058603 5.87180095 3.91816703 10.13769157 3.20578041 10.40144130 9.31452776 5.12841120 7.78969791 0.82192868 2.65510238 3.83259420 3.59826932 0.52929085 10.83188976 6.15571573 0.41066193 11.92995966 3.07288578 5.64047932 3.48007592 6.04690281 4.58337620 13.07434303 6.02713879 1.40558197 11.65835955 8.10817795 2.42542242 17.30459099 5.90844803 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 29735 k-point 2 : 0.5000 0.0000 0.0000 plane waves: 29568 k-point 3 : 0.0000 0.0000 0.5000 plane waves: 29692 k-point 4 : 0.5000 0.0000 0.5000 plane waves: 29632 maximum and minimum number of plane-waves per node : 29735 29568 maximum number of plane-waves: 29735 maximum index in each direction: IXMAX= 12 IYMAX= 32 IZMAX= 17 IXMIN= -12 IYMIN= -32 IZMIN= -18 The following grids will avoid any aliasing or wrap around errors in the Hartre e energy - symmetry arguments have not been applied - exchange correlation energies might require even more grid points - we recommend to set PREC=Normal or Accurate and rely on VASP defaults WARNING: aliasing errors must be expected set NGX to 50 to avoid them WARNING: aliasing errors must be expected set NGY to 140 to avoid them WARNING: aliasing errors must be expected set NGZ to 72 to avoid them serial 3D FFT for wavefunctions parallel 3D FFT for charge: minimum data exchange during FFTs selected (reduces bandwidth) total amount of memory used by VASP MPI-rank0 186670. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3165. kBytes fftplans : 15017. kBytes grid : 53760. kBytes one-center: 215. kBytes wavefun : 84513. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 25 NGY = 65 NGZ = 35 (NGX = 80 NGY =192 NGZ =108) gives a total of 56875 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 136.0000000 magnetization 0.1070000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 1353 Maximum index for augmentation-charges 1052 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.151 Maximum number of real-space cells 4x 2x 3 Maximum number of reciprocal cells 2x 4x 3 ----------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 1408 total energy-change (2. order) : 0.2381211E+04 (-0.3867803E+04) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -1955.38954846 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -299.15520249 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.01319442 eigenvalues EBANDS = -216.86228312 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 2381.21111426 eV energy without entropy = 2381.22430868 energy(sigma->0) = 2381.21551240 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 1992 total energy-change (2. order) :-0.7917324E+03 (-0.7647788E+03) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -1955.38954846 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -299.15520249 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.01682404 eigenvalues EBANDS = -1008.59108049 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1589.47868728 eV energy without entropy = 1589.49551132 energy(sigma->0) = 1589.48429529 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2080 total energy-change (2. order) :-0.2890415E+03 (-0.2836334E+03) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -1955.38954846 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -299.15520249 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = 0.00934575 eigenvalues EBANDS = -1297.65879741 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1300.43714014 eV energy without entropy = 1300.42779439 energy(sigma->0) = 1300.43402489 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 1712 total energy-change (2. order) :-0.7028003E+02 (-0.6903493E+02) number of electron 136.0000000 magnetization 0.1070000 augmentation part 136.0000000 magnetization 0.1070000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -1955.38954846 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -299.15520249 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.02580985 eigenvalues EBANDS = -1367.90367136 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1230.15711059 eV energy without entropy = 1230.18292044 energy(sigma->0) = 1230.16571387 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 1688 total energy-change (2. order) :-0.1658484E+02 (-0.1627892E+02) number of electron 136.0000110 magnetization -0.0168414 augmentation part -8.1573801 magnetization 0.7331182 Broyden mixing: rms(total) = 0.41468E+03 rms(broyden)= 0.41468E+03 rms(prec ) = 0.41470E+03 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -1955.38954846 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -299.15520249 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.00376770 eigenvalues EBANDS = -1384.51055037 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1213.57227373 eV energy without entropy = 1213.57604143 energy(sigma->0) = 1213.57352963 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 1888 total energy-change (2. order) : 0.2102764E+03 (-0.4922837E+02) number of electron 136.0000101 magnetization 0.1668904 augmentation part -7.3851929 magnetization -1.0355076 Broyden mixing: rms(total) = 0.13554E+03 rms(broyden)= 0.13554E+03 rms(prec ) = 0.13558E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0114 1.0114 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2308.64449003 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -284.95906823 PAW double counting = 16532999.71170914-16532450.01959492 entropy T*S EENTRO = -0.01187473 eigenvalues EBANDS = -823.69079653 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1423.84862484 eV energy without entropy = 1423.86049958 energy(sigma->0) = 1423.85258309 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 1736 total energy-change (2. order) :-0.4418670E+02 (-0.3372462E+02) number of electron 136.0000095 magnetization 0.2783272 augmentation part -7.4074857 magnetization -0.5293322 Broyden mixing: rms(total) = 0.98826E+02 rms(broyden)= 0.98826E+02 rms(prec ) = 0.98923E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 1.0239 1.4804 0.5674 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2303.39669449 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -281.83723984 PAW double counting = 20133464.12381430-20132918.53674290 entropy T*S EENTRO = -0.00751951 eigenvalues EBANDS = -872.14642851 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1379.66192919 eV energy without entropy = 1379.66944870 energy(sigma->0) = 1379.66443569 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 2056 total energy-change (2. order) : 0.3638123E+02 (-0.2774153E+02) number of electron 136.0000092 magnetization 0.2567567 augmentation part -7.1939488 magnetization 0.5296379 Broyden mixing: rms(total) = 0.76391E+02 rms(broyden)= 0.76391E+02 rms(prec ) = 0.76500E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7069 1.4760 0.5773 0.0673 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2373.50098069 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.65808369 PAW double counting = 26966802.38278396-26966259.31923831 entropy T*S EENTRO = 0.00079253 eigenvalues EBANDS = -767.32485967 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1416.04315428 eV energy without entropy = 1416.04236175 energy(sigma->0) = 1416.04289010 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 1664 total energy-change (2. order) : 0.1780021E+02 (-0.5800642E+01) number of electron 136.0000094 magnetization 0.2519985 augmentation part -7.4138541 magnetization 0.4113423 Broyden mixing: rms(total) = 0.51881E+02 rms(broyden)= 0.51881E+02 rms(prec ) = 0.51898E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7082 1.4880 0.5509 0.3970 0.3970 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2345.82456096 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -278.40575662 PAW double counting = 26149252.28415061-26148709.25252195 entropy T*S EENTRO = -0.01863389 eigenvalues EBANDS = -776.40204810 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1433.84336923 eV energy without entropy = 1433.86200312 energy(sigma->0) = 1433.84958053 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 1888 total energy-change (2. order) :-0.2392078E+02 (-0.3638949E+01) number of electron 136.0000091 magnetization 0.2390828 augmentation part -7.2936819 magnetization 0.4377846 Broyden mixing: rms(total) = 0.78147E+02 rms(broyden)= 0.78147E+02 rms(prec ) = 0.78290E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6493 1.4723 0.5584 0.5418 0.5418 0.1322 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2383.87348099 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.20534057 PAW double counting = 26485060.49845576-26484517.70780880 entropy T*S EENTRO = -0.00619028 eigenvalues EBANDS = -763.24578131 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1409.92259396 eV energy without entropy = 1409.92878424 energy(sigma->0) = 1409.92465739 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 1832 total energy-change (2. order) : 0.1857671E+02 (-0.4258956E+01) number of electron 136.0000094 magnetization 0.2452820 augmentation part -7.4736301 magnetization 0.2260620 Broyden mixing: rms(total) = 0.68481E+02 rms(broyden)= 0.68481E+02 rms(prec ) = 0.68555E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7007 1.4993 0.8011 0.8011 0.5919 0.2554 0.2554 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2372.73305734 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.53866106 PAW double counting = 25871687.02868491-25871144.18576537 entropy T*S EENTRO = 0.00937422 eigenvalues EBANDS = -755.54401330 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1428.49930222 eV energy without entropy = 1428.48992800 energy(sigma->0) = 1428.49617748 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 1784 total energy-change (2. order) :-0.1006462E+02 (-0.3753790E+01) number of electron 136.0000093 magnetization 0.3389101 augmentation part -7.4128500 magnetization -0.1249968 Broyden mixing: rms(total) = 0.65305E+02 rms(broyden)= 0.65305E+02 rms(prec ) = 0.65418E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8883 2.4350 1.4721 0.5804 0.5804 0.4668 0.4668 0.2167 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2379.63978921 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.79301403 PAW double counting = 25293870.95972940-25293329.06561509 entropy T*S EENTRO = -0.00136987 eigenvalues EBANDS = -758.48800149 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1418.43467986 eV energy without entropy = 1418.43604973 energy(sigma->0) = 1418.43513648 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 2088 total energy-change (2. order) :-0.1694430E+02 (-0.1331346E+02) number of electron 136.0000095 magnetization 0.3763828 augmentation part -7.2994621 magnetization 0.0489765 Broyden mixing: rms(total) = 0.27331E+02 rms(broyden)= 0.27331E+02 rms(prec ) = 0.27686E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8806 3.0478 1.2515 0.6864 0.6864 0.4152 0.4152 0.3269 0.2155 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2374.84486607 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.73033813 PAW double counting = 22943585.11909617-22943044.16983756 entropy T*S EENTRO = 0.02738802 eigenvalues EBANDS = -779.37380208 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1401.49038050 eV energy without entropy = 1401.46299248 energy(sigma->0) = 1401.48125116 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 2168 total energy-change (2. order) : 0.3210828E+02 (-0.1427393E+02) number of electron 136.0000094 magnetization 0.4005043 augmentation part -7.3619918 magnetization -0.6261008 Broyden mixing: rms(total) = 0.19787E+02 rms(broyden)= 0.19787E+02 rms(prec ) = 0.19882E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.8117 2.3259 1.4764 0.7378 0.7378 0.5157 0.5157 0.3892 0.3892 0.2176 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2351.15538220 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.03548749 PAW double counting = 23987639.96542582-23987099.56438215 entropy T*S EENTRO = -0.03583110 eigenvalues EBANDS = -771.03842644 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1433.59865660 eV energy without entropy = 1433.63448769 energy(sigma->0) = 1433.61060029 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 2000 total energy-change (2. order) :-0.6072389E+01 (-0.5784953E+01) number of electron 136.0000092 magnetization 0.3850286 augmentation part -7.2812128 magnetization -1.7278724 Broyden mixing: rms(total) = 0.28826E+02 rms(broyden)= 0.28826E+02 rms(prec ) = 0.28969E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7404 2.3702 1.5243 0.7114 0.7114 0.5095 0.5095 0.3869 0.3869 0.2176 0.0763 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2358.82472527 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.90304940 PAW double counting = 23696198.39815871-23695657.27583111 entropy T*S EENTRO = -0.01272447 eigenvalues EBANDS = -769.31830073 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1427.52626789 eV energy without entropy = 1427.53899235 energy(sigma->0) = 1427.53050938 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 1984 total energy-change (2. order) :-0.2191023E+01 (-0.1137762E+01) number of electron 136.0000093 magnetization 0.4127688 augmentation part -7.3547431 magnetization -1.2840807 Broyden mixing: rms(total) = 0.34793E+02 rms(broyden)= 0.34793E+02 rms(prec ) = 0.34926E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6788 2.2846 1.5887 0.7043 0.7043 0.5113 0.5113 0.3859 0.3859 0.2176 0.0865 0.0865 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2358.99309112 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.89086086 PAW double counting = 23600656.56059153-23600115.44913286 entropy T*S EENTRO = -0.02869163 eigenvalues EBANDS = -771.32630989 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1425.33524532 eV energy without entropy = 1425.36393694 energy(sigma->0) = 1425.34480919 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 2096 total energy-change (2. order) : 0.2541617E+01 (-0.4006406E+00) number of electron 136.0000092 magnetization 0.4149188 augmentation part -7.3596140 magnetization -0.7981680 Broyden mixing: rms(total) = 0.39617E+02 rms(broyden)= 0.39617E+02 rms(prec ) = 0.39710E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6436 2.3269 1.5354 0.7172 0.7172 0.5022 0.5022 0.3826 0.3826 0.2177 0.1766 0.1766 0.0859 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2351.23035269 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.13873992 PAW double counting = 23250778.52429308-23250237.22768312 entropy T*S EENTRO = -0.01413749 eigenvalues EBANDS = -776.49925766 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1427.87686236 eV energy without entropy = 1427.89099984 energy(sigma->0) = 1427.88157485 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 1992 total energy-change (2. order) : 0.7632402E+00 (-0.8323187E-01) number of electron 136.0000092 magnetization 0.3921616 augmentation part -7.3629024 magnetization -0.9450988 Broyden mixing: rms(total) = 0.41329E+02 rms(broyden)= 0.41329E+02 rms(prec ) = 0.41414E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6314 2.2667 1.5389 0.7420 0.7420 0.5051 0.5051 0.3151 0.3768 0.3768 0.2175 0.2342 0.2342 0.1539 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2349.81343140 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.18101660 PAW double counting = 23190877.49366003-23190336.15418147 entropy T*S EENTRO = -0.00657219 eigenvalues EBANDS = -777.16109600 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1428.64010252 eV energy without entropy = 1428.64667471 energy(sigma->0) = 1428.64229325 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 1896 total energy-change (2. order) : 0.1636339E+01 (-0.2130802E+00) number of electron 136.0000092 magnetization 0.4470798 augmentation part -7.3616696 magnetization -1.6180010 Broyden mixing: rms(total) = 0.37699E+02 rms(broyden)= 0.37699E+02 rms(prec ) = 0.37785E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6399 2.4348 1.4647 0.8008 0.4265 0.6503 0.4462 0.4462 0.4787 0.4787 0.3820 0.3820 0.2174 0.1750 0.1750 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2353.75872519 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.00853081 PAW double counting = 23352478.97240304-23351937.67951985 entropy T*S EENTRO = 0.00980333 eigenvalues EBANDS = -771.72172938 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1430.27644128 eV energy without entropy = 1430.26663796 energy(sigma->0) = 1430.27317351 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 2016 total energy-change (2. order) :-0.1544510E+01 (-0.2921193E+00) number of electron 136.0000093 magnetization 0.4799578 augmentation part -7.3710132 magnetization -1.3687943 Broyden mixing: rms(total) = 0.34302E+02 rms(broyden)= 0.34302E+02 rms(prec ) = 0.34424E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6417 2.5933 1.3436 0.8877 0.4815 0.5992 0.5992 0.6358 0.4589 0.4589 0.3633 0.3633 0.2150 0.2150 0.2172 0.1937 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2358.69715549 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.80178842 PAW double counting = 23553499.71616391-23552958.47953717 entropy T*S EENTRO = 0.00024942 eigenvalues EBANDS = -768.46874127 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1428.73193113 eV energy without entropy = 1428.73168171 energy(sigma->0) = 1428.73184799 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 1832 total energy-change (2. order) : 0.8648344E+01 (-0.1329509E+01) number of electron 136.0000093 magnetization 0.4696878 augmentation part -7.4045064 magnetization -2.2172994 Broyden mixing: rms(total) = 0.26015E+02 rms(broyden)= 0.26015E+02 rms(prec ) = 0.26054E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6121 2.6413 1.3199 0.9052 0.4787 0.6105 0.6105 0.6307 0.4594 0.4594 0.3601 0.3601 0.2269 0.2269 0.2169 0.2057 0.0815 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2361.20297727 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.24661323 PAW double counting = 23792306.89132150-23791765.72490803 entropy T*S EENTRO = -0.00311723 eigenvalues EBANDS = -757.79617114 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1437.38027474 eV energy without entropy = 1437.38339197 energy(sigma->0) = 1437.38131382 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 1904 total energy-change (2. order) :-0.6476969E-01 (-0.9986298E-01) number of electron 136.0000093 magnetization 0.4142931 augmentation part -7.4131990 magnetization -1.9622413 Broyden mixing: rms(total) = 0.24762E+02 rms(broyden)= 0.24762E+02 rms(prec ) = 0.24799E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5888 2.6483 1.3721 0.8385 0.6096 0.6096 0.4729 0.6592 0.4614 0.4614 0.3578 0.3578 0.2323 0.2323 0.2171 0.2021 0.1875 0.0891 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2360.39211245 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.28177636 PAW double counting = 23784255.56506182-23783714.40845809 entropy T*S EENTRO = 0.00279843 eigenvalues EBANDS = -758.63274845 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1437.31550505 eV energy without entropy = 1437.31270661 energy(sigma->0) = 1437.31457224 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 1960 total energy-change (2. order) : 0.4319453E+00 (-0.9337484E-01) number of electron 136.0000093 magnetization 0.5172081 augmentation part -7.4068505 magnetization -1.1477934 Broyden mixing: rms(total) = 0.25559E+02 rms(broyden)= 0.25559E+02 rms(prec ) = 0.25590E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5929 2.6286 1.3121 0.9277 0.4880 0.5695 0.5695 0.6362 0.4018 0.4018 0.4572 0.4572 0.3514 0.3514 0.2535 0.2535 0.2175 0.1972 0.1972 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2359.85103439 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.28565520 PAW double counting = 23554938.25144410-23554397.12036925 entropy T*S EENTRO = -0.03284702 eigenvalues EBANDS = -758.67682809 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1437.74745030 eV energy without entropy = 1437.78029732 energy(sigma->0) = 1437.75839931 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 1600 total energy-change (2. order) :-0.3600105E+00 (-0.3506048E+00) number of electron 136.0000093 magnetization 0.2699098 augmentation part -7.4049467 magnetization -2.4918674 Broyden mixing: rms(total) = 0.25008E+02 rms(broyden)= 0.25008E+02 rms(prec ) = 0.25049E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5705 2.6501 1.2996 0.9380 0.4788 0.6351 0.5640 0.5640 0.4552 0.4552 0.3562 0.3562 0.3467 0.3467 0.2590 0.2590 0.2177 0.2070 0.2070 0.2446 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2361.34313095 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.25280210 PAW double counting = 23686682.66025286-23686141.53653443 entropy T*S EENTRO = -0.06273080 eigenvalues EBANDS = -757.54035493 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1437.38743979 eV energy without entropy = 1437.45017059 energy(sigma->0) = 1437.40835006 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 1656 total energy-change (2. order) :-0.1422271E+01 (-0.8591817E+00) number of electron 136.0000093 magnetization 0.0040438 augmentation part -7.4039801 magnetization -1.1096515 Broyden mixing: rms(total) = 0.24267E+02 rms(broyden)= 0.24267E+02 rms(prec ) = 0.24327E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5927 2.5678 1.2402 1.0700 0.5895 0.5460 0.5759 0.5759 0.6213 0.5348 0.5348 0.4473 0.4473 0.3507 0.3507 0.2821 0.2821 0.2069 0.2069 0.2172 0.2056 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2363.35950188 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -276.12155578 PAW double counting = 23681327.77587884-23680786.70220335 entropy T*S EENTRO = -0.01172144 eigenvalues EBANDS = -757.07846824 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1435.96516830 eV energy without entropy = 1435.97688974 energy(sigma->0) = 1435.96907544 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 1464 total energy-change (2. order) : 0.2076933E+01 (-0.3532670E+01) number of electron 136.0000093 magnetization 0.1106358 augmentation part -7.4117501 magnetization 0.7584408 Broyden mixing: rms(total) = 0.24957E+02 rms(broyden)= 0.24957E+02 rms(prec ) = 0.24996E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5693 2.5685 1.2364 1.0741 0.6271 0.5305 0.5879 0.5879 0.5384 0.5384 0.6211 0.4485 0.4485 0.3524 0.3524 0.2791 0.2791 0.2172 0.2051 0.2048 0.2048 0.0531 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2363.29268177 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -275.34902565 PAW double counting = 23683409.60805951-23682868.56691530 entropy T*S EENTRO = 0.00881459 eigenvalues EBANDS = -755.82889012 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1438.04210141 eV energy without entropy = 1438.03328682 energy(sigma->0) = 1438.03916321 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 1448 total energy-change (2. order) :-0.7079345E+00 (-0.1726665E+00) number of electron 136.0000093 magnetization 1.0357835 augmentation part -7.4094147 magnetization 1.3198978 Broyden mixing: rms(total) = 0.25264E+02 rms(broyden)= 0.25264E+02 rms(prec ) = 0.25302E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6225 2.5795 1.2117 1.2117 0.8457 0.8457 0.7321 0.7321 0.5037 0.5480 0.5480 0.5961 0.4366 0.4366 0.3599 0.3599 0.3114 0.3114 0.2611 0.2200 0.2200 0.2176 0.2060 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2361.71906067 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -275.79448939 PAW double counting = 23663127.47286566-23662586.41147107 entropy T*S EENTRO = -0.00232255 eigenvalues EBANDS = -757.67409526 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1437.33416687 eV energy without entropy = 1437.33648943 energy(sigma->0) = 1437.33494106 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 1512 total energy-change (2. order) : 0.5864695E+01 (-0.1120371E+02) number of electron 136.0000094 magnetization 1.2552511 augmentation part -7.4462361 magnetization -1.1686761 Broyden mixing: rms(total) = 0.22714E+02 rms(broyden)= 0.22714E+02 rms(prec ) = 0.22730E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6451 2.6691 1.3028 1.3028 0.9460 0.9460 0.8639 0.8639 0.5052 0.5761 0.5761 0.5681 0.4356 0.4356 0.3789 0.3789 0.3251 0.3251 0.2874 0.2874 0.2207 0.2207 0.2173 0.2053 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2356.56507508 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -273.90988420 PAW double counting = 23661049.36723836-23660508.16860741 entropy T*S EENTRO = -0.03972309 eigenvalues EBANDS = -758.94782669 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1443.19886204 eV energy without entropy = 1443.23858513 energy(sigma->0) = 1443.21210307 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 1520 total energy-change (2. order) : 0.1157998E+02 (-0.2234148E+01) number of electron 136.0000094 magnetization 1.1326471 augmentation part -7.4250453 magnetization -1.2884671 Broyden mixing: rms(total) = 0.17495E+02 rms(broyden)= 0.17495E+02 rms(prec ) = 0.17518E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6202 2.6457 1.2988 1.2988 0.9562 0.9562 0.8663 0.8663 0.5050 0.5832 0.5832 0.5657 0.4348 0.4348 0.3796 0.3796 0.3226 0.3226 0.2775 0.2775 0.2208 0.2208 0.2173 0.2056 0.0658 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2365.60033861 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -265.19238517 PAW double counting = 23848676.56318804-23848135.44836010 entropy T*S EENTRO = 0.00833216 eigenvalues EBANDS = -747.01433713 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1454.77883934 eV energy without entropy = 1454.77050718 energy(sigma->0) = 1454.77606195 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 1704 total energy-change (2. order) : 0.1466604E+01 (-0.3215599E+00) number of electron 136.0000094 magnetization 1.2563050 augmentation part -7.4341047 magnetization -0.9871098 Broyden mixing: rms(total) = 0.17784E+02 rms(broyden)= 0.17784E+02 rms(prec ) = 0.17802E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6130 2.6590 1.2848 1.2848 0.9891 0.9891 0.8803 0.8803 0.5049 0.5883 0.5883 0.5675 0.4361 0.4361 0.3828 0.3828 0.3258 0.3258 0.2861 0.2861 0.2206 0.2206 0.2174 0.2052 0.1923 0.1923 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2362.88883468 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -264.45326557 PAW double counting = 23800158.58395031-23799617.42852033 entropy T*S EENTRO = -0.04563520 eigenvalues EBANDS = -748.98499153 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1456.24544315 eV energy without entropy = 1456.29107835 energy(sigma->0) = 1456.26065489 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 1832 total energy-change (2. order) : 0.1611608E+01 (-0.1412149E+00) number of electron 136.0000094 magnetization 1.1910150 augmentation part -7.4217330 magnetization -0.8917902 Broyden mixing: rms(total) = 0.18171E+02 rms(broyden)= 0.18171E+02 rms(prec ) = 0.18193E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6141 2.6652 1.2775 1.2775 1.0023 1.0023 0.8942 0.8942 0.5049 0.5727 0.5727 0.5691 0.3595 0.3595 0.4341 0.4341 0.3830 0.3830 0.3297 0.3297 0.3013 0.3013 0.2200 0.2200 0.2561 0.2176 0.2045 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2361.75213012 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -262.76929178 PAW double counting = 23714541.16277242-23713999.99037715 entropy T*S EENTRO = -0.03325571 eigenvalues EBANDS = -750.22340695 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1457.85705086 eV energy without entropy = 1457.89030657 energy(sigma->0) = 1457.86813610 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 1544 total energy-change (2. order) :-0.7702772E+01 (-0.5701886E+00) number of electron 136.0000094 magnetization 1.7341549 augmentation part -7.4349992 magnetization -0.4656896 Broyden mixing: rms(total) = 0.17284E+02 rms(broyden)= 0.17284E+02 rms(prec ) = 0.17300E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6536 2.7503 1.1080 1.1080 1.2212 1.2212 1.1250 1.1250 0.5049 0.7098 0.7098 0.5565 0.5565 0.5724 0.4197 0.4197 0.3934 0.3934 0.3198 0.3198 0.3117 0.3117 0.3128 0.3128 0.2202 0.2202 0.2175 0.2048 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2363.16344084 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -269.10850927 PAW double counting = 23824730.43120521-23824189.25670666 entropy T*S EENTRO = -0.02555212 eigenvalues EBANDS = -750.18545731 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1450.15427917 eV energy without entropy = 1450.17983128 energy(sigma->0) = 1450.16279654 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 1520 total energy-change (2. order) : 0.1960126E+01 (-0.1665630E+01) number of electron 136.0000094 magnetization 2.3582889 augmentation part -7.4324959 magnetization -0.3934488 Broyden mixing: rms(total) = 0.15078E+02 rms(broyden)= 0.15078E+02 rms(prec ) = 0.15094E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6326 2.7445 1.1114 1.1114 1.1425 1.1425 1.2200 1.2200 0.5049 0.7104 0.7104 0.5584 0.5584 0.5737 0.4204 0.4204 0.3926 0.3926 0.3176 0.3176 0.3110 0.3110 0.3117 0.3117 0.2202 0.2202 0.2175 0.2048 0.0347 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2365.21360043 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -268.29502714 PAW double counting = 23967122.74071884-23966581.56433105 entropy T*S EENTRO = -0.00983594 eigenvalues EBANDS = -747.00625897 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1452.11440545 eV energy without entropy = 1452.12424139 energy(sigma->0) = 1452.11768410 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 1784 total energy-change (2. order) : 0.9236700E+01 (-0.2664531E+01) number of electron 136.0000095 magnetization 3.0028203 augmentation part -7.3350796 magnetization -1.1687729 Broyden mixing: rms(total) = 0.14157E+02 rms(broyden)= 0.14157E+02 rms(prec ) = 0.14167E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6206 2.7480 1.1139 1.1139 1.2201 1.2201 1.1307 1.1307 0.5049 0.7139 0.7139 0.5547 0.5547 0.5737 0.4196 0.4196 0.3936 0.3936 0.3230 0.3230 0.3099 0.3099 0.3127 0.3127 0.2202 0.2202 0.2175 0.2048 0.1616 0.1616 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2364.64193928 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -259.81585051 PAW double counting = 24010527.18400338-24009985.89675249 entropy T*S EENTRO = -0.01202972 eigenvalues EBANDS = -746.92906566 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1461.35110588 eV energy without entropy = 1461.36313560 energy(sigma->0) = 1461.35511579 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 1712 total energy-change (2. order) : 0.9420770E+01 (-0.3686427E+00) number of electron 136.0000095 magnetization 3.9501197 augmentation part -7.3412528 magnetization -0.1291554 Broyden mixing: rms(total) = 0.14093E+02 rms(broyden)= 0.14093E+02 rms(prec ) = 0.14103E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6092 2.7493 1.1153 1.1153 1.2199 1.2199 1.1240 1.1240 0.7145 0.7145 0.5049 0.5528 0.5528 0.5735 0.4187 0.4187 0.3938 0.3938 0.3236 0.3236 0.3117 0.3117 0.3126 0.3126 0.2202 0.2202 0.2175 0.2048 0.2364 0.2364 0.1396 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2363.00332000 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -249.18378377 PAW double counting = 24005996.58226701-24005456.16011847 entropy T*S EENTRO = -0.02375219 eigenvalues EBANDS = -748.90215703 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1470.77187571 eV energy without entropy = 1470.79562790 energy(sigma->0) = 1470.77979311 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 1672 total energy-change (2. order) : 0.1439435E+02 (-0.2608576E+00) number of electron 136.0000095 magnetization 4.3612795 augmentation part -7.3575524 magnetization 0.0592727 Broyden mixing: rms(total) = 0.13696E+02 rms(broyden)= 0.13696E+02 rms(prec ) = 0.13704E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6429 2.7253 1.1233 1.1233 1.0978 1.0978 1.2004 1.2004 0.8184 0.8184 0.5049 0.6763 0.6763 0.5707 0.5707 0.5794 0.4283 0.4283 0.3232 0.3232 0.3917 0.3917 0.3482 0.3482 0.3307 0.3307 0.3307 0.3097 0.2202 0.2202 0.2175 0.2048 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2360.86237714 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -234.62559198 PAW double counting = 24110969.94525092-24110429.93374831 entropy T*S EENTRO = -0.02509423 eigenvalues EBANDS = -750.79495103 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1485.16622837 eV energy without entropy = 1485.19132261 energy(sigma->0) = 1485.17459312 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 2088 total energy-change (2. order) : 0.5954360E+01 (-0.5896193E-01) number of electron 136.0000095 magnetization 4.5528430 augmentation part -7.3676549 magnetization 0.1804241 Broyden mixing: rms(total) = 0.13698E+02 rms(broyden)= 0.13698E+02 rms(prec ) = 0.13707E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6228 2.7257 1.1233 1.1233 1.2005 1.2005 1.0973 1.0973 0.8163 0.8163 0.5049 0.6763 0.6763 0.5705 0.5705 0.5794 0.4283 0.4283 0.3233 0.3233 0.3916 0.3916 0.3481 0.3481 0.3309 0.3309 0.3309 0.3096 0.2202 0.2202 0.2175 0.2048 0.0048 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2360.12665367 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -229.29430044 PAW double counting = 24158525.59739518-24157985.09305359 entropy T*S EENTRO = -0.02168275 eigenvalues EBANDS = -751.40385661 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1491.12058826 eV energy without entropy = 1491.14227101 energy(sigma->0) = 1491.12781584 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 2000 total energy-change (2. order) : 0.1130483E+01 (-0.1566251E-01) number of electron 136.0000095 magnetization 5.6128224 augmentation part -7.3711571 magnetization 1.1172885 Broyden mixing: rms(total) = 0.13612E+02 rms(broyden)= 0.13612E+02 rms(prec ) = 0.13621E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6330 2.7530 1.1151 1.1151 1.1931 1.1931 1.0834 1.0834 0.9518 0.9518 0.5049 0.6878 0.6878 0.5596 0.5596 0.5802 0.3673 0.3673 0.4221 0.4221 0.3262 0.3262 0.3821 0.3821 0.3345 0.3345 0.3486 0.3486 0.3396 0.3043 0.2202 0.2202 0.2175 0.2048 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2359.98522014 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -228.08707894 PAW double counting = 24215505.16213197-24214964.69966977 entropy T*S EENTRO = -0.01529742 eigenvalues EBANDS = -751.58653472 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1492.25107113 eV energy without entropy = 1492.26636855 energy(sigma->0) = 1492.25617027 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 1632 total energy-change (2. order) : 0.5418596E+01 (-0.1909291E+00) number of electron 136.0000096 magnetization 6.5894687 augmentation part -7.3941886 magnetization 1.7352446 Broyden mixing: rms(total) = 0.13072E+02 rms(broyden)= 0.13072E+02 rms(prec ) = 0.13078E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6698 2.7260 1.4154 1.4154 1.0627 1.0627 1.1155 1.1155 1.2125 1.2125 0.5049 0.7680 0.7680 0.6272 0.6272 0.5285 0.5285 0.5814 0.4172 0.4172 0.3305 0.3305 0.3950 0.3950 0.3509 0.3509 0.3235 0.3235 0.3549 0.3549 0.2961 0.2202 0.2202 0.2175 0.2048 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2359.26914830 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -221.95788370 PAW double counting = 24451100.39359318-24450560.50864451 entropy T*S EENTRO = 0.01045693 eigenvalues EBANDS = -752.46144635 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1497.66966741 eV energy without entropy = 1497.65921048 energy(sigma->0) = 1497.66618176 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 1808 total energy-change (2. order) : 0.2711978E+01 (-0.1178867E+00) number of electron 136.0000096 magnetization 7.7078696 augmentation part -7.4110555 magnetization 2.6217869 Broyden mixing: rms(total) = 0.12570E+02 rms(broyden)= 0.12570E+02 rms(prec ) = 0.12576E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6871 2.7873 1.5880 1.5880 1.0603 1.0603 1.1808 1.1808 1.1520 1.1520 0.8972 0.8972 0.5049 0.6371 0.6371 0.5625 0.5625 0.6144 0.4145 0.4145 0.3323 0.3323 0.3762 0.3762 0.4065 0.4065 0.3254 0.3254 0.3763 0.3763 0.3615 0.2989 0.2202 0.2202 0.2048 0.2175 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = 511.37495245 -Hartree energ DENC = -2357.74655001 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -218.44587716 PAW double counting = 24623115.88262162-24622576.18532139 entropy T*S EENTRO = 0.00980686 eigenvalues EBANDS = -754.59577438 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1500.38164568 eV energy without entropy = 1500.37183883 energy(sigma->0) = 1500.37837673 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- ----------------------------------------------------------------------------- | | | EEEEEEE RRRRRR RRRRRR OOOOOOO RRRRRR ### ### ### | | E R R R R O O R R ### ### ### | | E R R R R O O R R ### ### ### | | EEEEE RRRRRR RRRRRR O O RRRRRR # # # | | E R R R R O O R R | | E R R R R O O R R ### ### ### | | EEEEEEE R R R R OOOOOOO R R ### ### ### | | | | Error EDDDAV: Call to ZHEGV failed. Returncode = 7 1 8 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------