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.03 07:52:13 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.040 0.539 0.566- 30 1.19 23 1.97 8 2.59 2 0.414 0.810 0.901- 3 0.720 0.930 0.201- 4 0.754 0.669 0.194- 5 0.045 0.036 0.800- 12 0.86 34 1.43 6 0.764 0.387 0.236- 7 0.837 0.513 0.766- 23 2.16 15 2.47 25 2.56 8 0.729 0.491 0.531- 30 2.39 1 2.59 9 0.478 0.713 0.394- 14 1.36 10 0.034 0.981 0.434- 11 0.420 0.418 0.425- 27 0.33 24 2.54 12 0.072 0.001 0.845- 34 0.76 5 0.86 13 0.773 0.193 0.116- 16 0.89 29 1.53 22 2.46 14 0.562 0.773 0.409- 9 1.36 26 1.96 15 0.903 0.424 0.920- 31 1.38 7 2.47 16 0.717 0.173 0.179- 13 0.89 29 2.04 17 0.088 0.653 0.948- 18 0.225 0.111 0.637- 19 0.701 0.069 0.753- 20 0.585 0.224 0.640- 24 2.49 21 0.425 0.074 0.318- 22 0.995 0.121 0.017- 29 0.98 13 2.46 23 0.929 0.605 0.679- 1 1.97 7 2.16 24 0.400 0.297 0.504- 27 2.45 20 2.49 11 2.54 25 0.538 0.456 0.788- 7 2.56 26 0.765 0.807 0.317- 14 1.96 27 0.380 0.414 0.431- 11 0.33 24 2.45 28 0.646 0.564 0.988- 29 0.887 0.143 0.044- 22 0.98 13 1.53 16 2.04 30 0.040 0.481 0.538- 1 1.19 8 2.39 31 0.981 0.378 0.999- 15 1.38 32 0.175 0.892 0.053- 33 0.295 0.500 0.182- 34 0.006 0.980 0.881- 12 0.76 5 1.43 35 0.925 0.869 0.665- ----------------------------------------------------------------------------- | | | 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.040084050 0.539101290 0.566442040 0.414100960 0.810354090 0.900721610 0.720185740 0.930404800 0.200723340 0.754178430 0.669000070 0.194067830 0.045226420 0.035589800 0.800431210 0.763875290 0.386933070 0.236363670 0.837251110 0.512787300 0.766339030 0.729343580 0.491322960 0.531035170 0.478278400 0.712690500 0.394064910 0.034425570 0.980550310 0.434354320 0.420481510 0.417879290 0.425097350 0.072280760 0.000685190 0.844506860 0.773075090 0.193270120 0.116252240 0.562050180 0.772965750 0.409190760 0.903160590 0.423617210 0.920218460 0.716933200 0.173387850 0.178518110 0.087870790 0.652577890 0.948092130 0.225144700 0.111255340 0.637257800 0.700932640 0.069357080 0.752769100 0.584514000 0.224354860 0.639596280 0.425361420 0.074349730 0.318023010 0.994860740 0.121042140 0.016552080 0.928960680 0.604922800 0.679040700 0.399567690 0.296650920 0.504246240 0.537532200 0.455929550 0.787674970 0.765102110 0.806556100 0.316896140 0.379739200 0.414080020 0.431084500 0.646498040 0.564230140 0.987555810 0.887305990 0.143314830 0.043941850 0.039659380 0.480778600 0.537730740 0.981281080 0.377966040 0.999224720 0.174621130 0.891900150 0.053411440 0.295445260 0.500124710 0.181681640 0.006275010 0.979879630 0.881409290 0.925275230 0.869485830 0.665480920 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.04008405 0.53910129 0.56644204 0.41410096 0.81035409 0.90072161 0.72018574 0.93040480 0.20072334 0.75417843 0.66900007 0.19406783 0.04522642 0.03558980 0.80043121 0.76387529 0.38693307 0.23636367 0.83725111 0.51278730 0.76633903 0.72934358 0.49132296 0.53103517 0.47827840 0.71269050 0.39406491 0.03442557 0.98055031 0.43435432 0.42048151 0.41787929 0.42509735 0.07228076 0.00068519 0.84450686 0.77307509 0.19327012 0.11625224 0.56205018 0.77296575 0.40919076 0.90316059 0.42361721 0.92021846 0.71693320 0.17338785 0.17851811 0.08787079 0.65257789 0.94809213 0.22514470 0.11125534 0.63725780 0.70093264 0.06935708 0.75276910 0.58451400 0.22435486 0.63959628 0.42536142 0.07434973 0.31802301 0.99486074 0.12104214 0.01655208 0.92896068 0.60492280 0.67904070 0.39956769 0.29665092 0.50424624 0.53753220 0.45592955 0.78767497 0.76510211 0.80655610 0.31689614 0.37973920 0.41408002 0.43108450 0.64649804 0.56423014 0.98755581 0.88730599 0.14331483 0.04394185 0.03965938 0.48077860 0.53773074 0.98128108 0.37796604 0.99922472 0.17462113 0.89190015 0.05341144 0.29544526 0.50012471 0.18168164 0.00627501 0.97987963 0.88140929 0.92527523 0.86948583 0.66548092 position of ions in cartesian coordinates (Angst): 0.30716808 10.60040258 6.13867966 3.17329707 15.93407351 9.76135428 5.51885534 18.29464262 2.17529102 5.77934473 13.15461528 2.10316353 0.34657458 0.69980580 8.67448113 5.85365273 7.60830365 2.56153455 6.41593898 10.08298796 8.30501532 5.58903279 9.66093249 5.75496621 3.66509521 14.01370457 4.27058389 0.26380659 19.28065880 4.70721070 3.22219186 8.21680227 4.60689051 0.55389469 0.01347296 9.15214041 5.92415172 3.80028970 1.25985575 4.30704673 15.19890284 4.43450666 6.92100992 8.32962756 9.97264671 5.49393080 3.40934263 1.93464717 0.67336265 12.83170431 10.27472092 1.72530635 2.18762488 6.90612847 5.37131691 1.36377520 8.15795446 4.47918923 4.41151205 6.93147118 3.25958710 1.46194618 3.44649805 7.62371734 2.38006370 0.17937919 7.11871859 11.89465751 7.35894062 3.06192717 5.83307671 5.46464761 4.11916300 8.96498833 8.53623845 5.86305398 15.85939325 3.43428586 2.90997946 8.14209684 4.67177481 4.95417913 11.09451367 10.70239908 6.79951453 2.81801383 0.47620925 0.30391379 9.45359769 5.82752784 7.51965504 7.43198404 10.82885809 1.33813918 17.53752184 0.57883366 2.26402657 9.83400219 1.96893117 0.04808603 19.26747115 9.55206164 7.09047662 17.09678682 7.21198976 -------------------------------------------------------------------------------------------------------- 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 186669. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 3164. 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 1350 Maximum index for augmentation-charges 1055 (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.1413368E+04 (-0.3841669E+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 = -729.09930942 -Hartree energ DENC = -1709.77674572 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -297.48019331 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.01709489 eigenvalues EBANDS = -191.51481576 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 1413.36823119 eV energy without entropy = 1413.38532608 energy(sigma->0) = 1413.37392949 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 1960 total energy-change (2. order) :-0.6310814E+03 (-0.6090385E+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 = -729.09930942 -Hartree energ DENC = -1709.77674572 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -297.48019331 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.00118476 eigenvalues EBANDS = -822.61213658 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 782.28682051 eV energy without entropy = 782.28800526 energy(sigma->0) = 782.28721543 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 2344 total energy-change (2. order) :-0.1495715E+03 (-0.1408349E+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 = -729.09930942 -Hartree energ DENC = -1709.77674572 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -297.48019331 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = -0.01340296 eigenvalues EBANDS = -972.17145174 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 632.71528713 eV energy without entropy = 632.72869009 energy(sigma->0) = 632.71975479 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 2272 total energy-change (2. order) :-0.2390537E+02 (-0.2326569E+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 = -729.09930942 -Hartree energ DENC = -1709.77674572 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -297.48019331 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = 0.02138923 eigenvalues EBANDS = -996.11161177 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 608.80991930 eV energy without entropy = 608.78853007 energy(sigma->0) = 608.80278956 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 2000 total energy-change (2. order) :-0.1182285E+01 (-0.1166125E+01) number of electron 135.9999951 magnetization 0.1011748 augmentation part -8.7608769 magnetization 0.0350274 Broyden mixing: rms(total) = 0.11970E+03 rms(broyden)= 0.11970E+03 rms(prec ) = 0.11979E+03 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1709.77674572 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -297.48019331 PAW double counting = 3170.77583059 -2609.60722799 entropy T*S EENTRO = 0.01702335 eigenvalues EBANDS = -997.28953077 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 607.62763441 eV energy without entropy = 607.61061106 energy(sigma->0) = 607.62195996 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 2688 total energy-change (2. order) :-0.3179880E+02 (-0.2128668E+03) number of electron 135.9999979 magnetization 0.1574952 augmentation part -8.7253008 magnetization -0.4671447 Broyden mixing: rms(total) = 0.42138E+02 rms(broyden)= 0.42138E+02 rms(prec ) = 0.42548E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7469 0.7469 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -2048.68214658 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -284.35051157 PAW double counting = 1639688.81070415 -1639134.94922529 entropy T*S EENTRO = 0.02103777 eigenvalues EBANDS = -696.00950092 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 575.82883584 eV energy without entropy = 575.80779807 energy(sigma->0) = 575.82182325 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 2464 total energy-change (2. order) : 0.7814956E+02 (-0.1064579E+03) number of electron 135.9999947 magnetization 0.1377005 augmentation part -7.8869472 magnetization 0.2501207 Broyden mixing: rms(total) = 0.14459E+02 rms(broyden)= 0.14459E+02 rms(prec ) = 0.15242E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7088 0.8935 0.5241 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1845.88203354 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -288.48277791 PAW double counting = 941099.81908662 -940546.22638630 entropy T*S EENTRO = -0.00291231 eigenvalues EBANDS = -816.23505454 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 653.97840030 eV energy without entropy = 653.98131260 energy(sigma->0) = 653.97937107 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 2168 total energy-change (2. order) : 0.4424803E+02 (-0.4295423E+02) number of electron 135.9999948 magnetization 0.1353062 augmentation part -7.4922524 magnetization 0.0240203 Broyden mixing: rms(total) = 0.28635E+02 rms(broyden)= 0.28635E+02 rms(prec ) = 0.28692E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6178 0.8509 0.6639 0.3386 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1888.39934742 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -287.51876096 PAW double counting = 978247.10973070 -977694.17646943 entropy T*S EENTRO = 0.01663361 eigenvalues EBANDS = -729.79383470 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 698.22643007 eV energy without entropy = 698.20979646 energy(sigma->0) = 698.22088553 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 2056 total energy-change (2. order) :-0.2487319E+02 (-0.1426877E+02) number of electron 135.9999951 magnetization 0.2100345 augmentation part -7.5951131 magnetization 0.2903495 Broyden mixing: rms(total) = 0.35721E+02 rms(broyden)= 0.35721E+02 rms(prec ) = 0.35871E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5619 0.8011 0.5483 0.4490 0.4490 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1892.65782309 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -287.77227697 PAW double counting = 955831.26919568 -955278.21850595 entropy T*S EENTRO = -0.05327203 eigenvalues EBANDS = -750.20256063 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 673.35323527 eV energy without entropy = 673.40650731 energy(sigma->0) = 673.37099262 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 2168 total energy-change (2. order) :-0.1813129E+02 (-0.1627768E+02) number of electron 135.9999962 magnetization 0.2420602 augmentation part -8.1003150 magnetization 0.2866653 Broyden mixing: rms(total) = 0.15102E+02 rms(broyden)= 0.15102E+02 rms(prec ) = 0.16033E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5982 0.7527 0.7527 0.6566 0.6566 0.1724 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1945.87102540 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -285.00311779 PAW double counting = 1074560.55320517 -1074008.08398917 entropy T*S EENTRO = -0.01291041 eigenvalues EBANDS = -717.34869865 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 655.22194201 eV energy without entropy = 655.23485242 energy(sigma->0) = 655.22624548 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 2184 total energy-change (2. order) : 0.1545109E+02 (-0.7388542E+01) number of electron 135.9999980 magnetization 0.2603692 augmentation part -9.2382078 magnetization 0.4878897 Broyden mixing: rms(total) = 0.16002E+02 rms(broyden)= 0.16002E+02 rms(prec ) = 0.16714E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6043 0.9293 0.9293 0.5963 0.5963 0.2872 0.2872 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1934.00232090 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -284.19379044 PAW double counting = 1044657.50054462 -1044105.38494803 entropy T*S EENTRO = -0.01416633 eigenvalues EBANDS = -714.22076585 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 670.67303135 eV energy without entropy = 670.68719768 energy(sigma->0) = 670.67775346 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 2128 total energy-change (2. order) :-0.1784531E+01 (-0.2962485E+01) number of electron 135.9999957 magnetization 0.3354639 augmentation part -7.8396721 magnetization 0.0319393 Broyden mixing: rms(total) = 0.18809E+02 rms(broyden)= 0.18809E+02 rms(prec ) = 0.19160E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6623 1.1606 1.1606 0.7317 0.5404 0.4189 0.4189 0.2053 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1898.16646836 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -285.22709759 PAW double counting = 1085723.49311002 -1085171.72459966 entropy T*S EENTRO = -0.03273626 eigenvalues EBANDS = -750.44218585 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 668.88850058 eV energy without entropy = 668.92123684 energy(sigma->0) = 668.89941267 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 1904 total energy-change (2. order) : 0.5664086E+01 (-0.2606950E+01) number of electron 135.9999961 magnetization 0.3756352 augmentation part -7.9953031 magnetization 0.2306351 Broyden mixing: rms(total) = 0.18676E+02 rms(broyden)= 0.18676E+02 rms(prec ) = 0.19135E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7654 2.1248 0.9026 0.8166 0.8166 0.4767 0.4767 0.2356 0.2740 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1930.40293227 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -282.31512892 PAW double counting = 1167030.07767179 -1166479.20432899 entropy T*S EENTRO = -0.01985794 eigenvalues EBANDS = -714.57131493 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 674.55258702 eV energy without entropy = 674.57244496 energy(sigma->0) = 674.55920634 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 2080 total energy-change (2. order) :-0.2143616E-01 (-0.6196707E+01) number of electron 135.9999979 magnetization 0.4858380 augmentation part -9.2054629 magnetization -1.0614532 Broyden mixing: rms(total) = 0.18215E+02 rms(broyden)= 0.18215E+02 rms(prec ) = 0.18628E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.7303 1.9390 1.3160 0.7663 0.7663 0.5028 0.5028 0.2162 0.2817 0.2817 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1950.41395556 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -280.24436332 PAW double counting = 1239765.90549958 -1239215.47698410 entropy T*S EENTRO = -0.00852278 eigenvalues EBANDS = -696.21900123 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 674.53115086 eV energy without entropy = 674.53967364 energy(sigma->0) = 674.53399179 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 1752 total energy-change (2. order) : 0.5784200E+01 (-0.4960339E+01) number of electron 135.9999976 magnetization 0.4274213 augmentation part -9.2943308 magnetization -1.4500905 Broyden mixing: rms(total) = 0.18854E+02 rms(broyden)= 0.18854E+02 rms(prec ) = 0.19222E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6692 2.1789 0.9465 0.7737 0.7737 0.4736 0.4736 0.2182 0.2886 0.2886 0.2763 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1949.80638113 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -278.14440998 PAW double counting = 1310301.27428453 -1309750.83830958 entropy T*S EENTRO = -0.02747824 eigenvalues EBANDS = -693.13083278 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 680.31535109 eV energy without entropy = 680.34282933 energy(sigma->0) = 680.32451051 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 1744 total energy-change (2. order) :-0.2512508E+01 (-0.6436635E+00) number of electron 135.9999979 magnetization 0.4371795 augmentation part -9.2868108 magnetization -0.8285402 Broyden mixing: rms(total) = 0.19596E+02 rms(broyden)= 0.19596E+02 rms(prec ) = 0.19957E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6198 2.2012 0.9436 0.7648 0.7648 0.4652 0.4652 0.2880 0.2880 0.2187 0.2094 0.2094 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1945.81053813 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.48072021 PAW double counting = 1276609.20312073 -1276058.84644981 entropy T*S EENTRO = -0.02470231 eigenvalues EBANDS = -698.22634553 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 677.80284300 eV energy without entropy = 677.82754531 energy(sigma->0) = 677.81107710 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 2192 total energy-change (2. order) :-0.1293255E+03 (-0.1025702E+03) number of electron 135.9999936 magnetization 0.4536390 augmentation part -9.1375336 magnetization -0.9317860 Broyden mixing: rms(total) = 0.21954E+02 rms(broyden)= 0.21926E+02 rms(prec ) = 0.22388E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5722 2.1936 0.9489 0.7660 0.7660 0.4650 0.4650 0.2881 0.2881 0.2188 0.2318 0.2318 0.0029 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1945.85572300 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.37170696 PAW double counting = 1278578.07961890 -1278027.72650643 entropy T*S EENTRO = -0.03856103 eigenvalues EBANDS = -827.59828831 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 548.47731144 eV energy without entropy = 548.51587247 energy(sigma->0) = 548.49016512 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 1632 total energy-change (2. order) : 0.1383414E+03 (-0.3049598E+03) number of electron 135.9999980 magnetization 0.2233417 augmentation part -9.2500017 magnetization -3.0349091 Broyden mixing: rms(total) = 0.20438E+02 rms(broyden)= 0.20418E+02 rms(prec ) = 0.20754E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5398 2.2032 0.9363 0.7646 0.7646 0.4704 0.4704 0.2933 0.2867 0.2867 0.2182 0.1577 0.1577 0.0080 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1946.10934468 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -275.76090979 PAW double counting = 1284482.95374918 -1283932.54521992 entropy T*S EENTRO = -0.02753466 eigenvalues EBANDS = -692.68048047 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.81873793 eV energy without entropy = 686.84627259 energy(sigma->0) = 686.82791615 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 1632 total energy-change (2. order) :-0.3070122E+01 (-0.1606574E+01) number of electron 135.9999977 magnetization 0.4837303 augmentation part -9.2998985 magnetization -1.7627273 Broyden mixing: rms(total) = 0.20764E+02 rms(broyden)= 0.20764E+02 rms(prec ) = 0.21083E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5083 2.2047 0.9357 0.7650 0.7650 0.4703 0.4703 0.3250 0.2866 0.2866 0.2183 0.1743 0.1743 0.0083 0.0313 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1944.57662122 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -278.38420265 PAW double counting = 1275097.34667365 -1274546.85116103 entropy T*S EENTRO = -0.01796290 eigenvalues EBANDS = -694.75658776 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 683.74861636 eV energy without entropy = 683.76657926 energy(sigma->0) = 683.75460399 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 1744 total energy-change (2. order) :-0.3958964E+01 (-0.9634934E+00) number of electron 135.9999980 magnetization 0.5040475 augmentation part -9.2931265 magnetization -1.1092657 Broyden mixing: rms(total) = 0.20199E+02 rms(broyden)= 0.20199E+02 rms(prec ) = 0.20535E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5445 2.2294 0.9364 0.7680 0.7680 0.5066 0.5066 0.4722 0.4722 0.2904 0.2904 0.2848 0.2848 0.2202 0.1288 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1944.23761835 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.43899339 PAW double counting = 1282414.34715778 -1281863.93138434 entropy T*S EENTRO = -0.02843904 eigenvalues EBANDS = -697.90954863 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 679.78965231 eV energy without entropy = 679.81809135 energy(sigma->0) = 679.79913199 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 2080 total energy-change (2. order) :-0.4392888E+00 (-0.2278169E+00) number of electron 135.9999980 magnetization 0.8284817 augmentation part -9.2919020 magnetization -0.7668945 Broyden mixing: rms(total) = 0.19913E+02 rms(broyden)= 0.19913E+02 rms(prec ) = 0.20253E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5686 2.2797 0.6841 0.6841 0.9423 0.7772 0.7772 0.4516 0.4516 0.4127 0.4127 0.3106 0.3106 0.2242 0.2242 0.1472 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1944.21556944 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -279.43875861 PAW double counting = 1283694.54054762 -1283144.12795990 entropy T*S EENTRO = -0.01733151 eigenvalues EBANDS = -698.37904288 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 679.35036356 eV energy without entropy = 679.36769507 energy(sigma->0) = 679.35614073 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 1528 total energy-change (2. order) : 0.7587218E+01 (-0.2077242E+01) number of electron 135.9999975 magnetization 1.3765030 augmentation part -9.3243021 magnetization -1.3893723 Broyden mixing: rms(total) = 0.19231E+02 rms(broyden)= 0.19231E+02 rms(prec ) = 0.19568E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5739 2.3436 0.8348 0.8348 0.9017 0.7622 0.7622 0.4542 0.4542 0.4627 0.4627 0.3198 0.3198 0.2191 0.2376 0.2376 0.1416 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1946.64629079 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -277.73563402 PAW double counting = 1315007.90697993 -1314457.62277333 entropy T*S EENTRO = -0.00166875 eigenvalues EBANDS = -689.95150982 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 686.93758148 eV energy without entropy = 686.93925024 energy(sigma->0) = 686.93813773 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 1512 total energy-change (2. order) : 0.1100060E+02 (-0.1457892E+01) number of electron 135.9999976 magnetization 1.9374573 augmentation part -9.3121147 magnetization -1.5882447 Broyden mixing: rms(total) = 0.18481E+02 rms(broyden)= 0.18481E+02 rms(prec ) = 0.18805E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5783 2.3412 0.9775 0.9775 0.9109 0.7526 0.7526 0.4581 0.4581 0.4749 0.4749 0.2966 0.2966 0.3091 0.3091 0.2375 0.2302 0.1440 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1947.90170591 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -272.13385413 PAW double counting = 1334238.40798107 -1333687.97451924 entropy T*S EENTRO = 0.01863249 eigenvalues EBANDS = -683.46683484 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 697.93817772 eV energy without entropy = 697.91954523 energy(sigma->0) = 697.93196689 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 1792 total energy-change (2. order) : 0.1027267E+02 (-0.1257317E+01) number of electron 135.9999974 magnetization 3.3872080 augmentation part -9.2261515 magnetization -0.6793885 Broyden mixing: rms(total) = 0.17340E+02 rms(broyden)= 0.17340E+02 rms(prec ) = 0.17652E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6427 2.4583 1.4100 1.4100 0.8912 0.7471 0.7471 0.5260 0.5260 0.5208 0.5208 0.4057 0.4057 0.3596 0.3596 0.2741 0.2741 0.2235 0.1435 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1948.34983199 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -267.21971694 PAW double counting = 1345735.27527773 -1345184.86093110 entropy T*S EENTRO = 0.01181158 eigenvalues EBANDS = -677.63424167 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 708.21084588 eV energy without entropy = 708.19903430 energy(sigma->0) = 708.20690869 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 1664 total energy-change (2. order) : 0.2418419E+02 (-0.3969721E+01) number of electron 135.9999962 magnetization 3.8683411 augmentation part -8.4267489 magnetization -1.1519827 Broyden mixing: rms(total) = 0.15668E+02 rms(broyden)= 0.15668E+02 rms(prec ) = 0.15851E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6298 2.4819 1.4808 1.4808 0.8976 0.7175 0.7175 0.5352 0.5352 0.5185 0.5185 0.4185 0.4185 0.3286 0.3286 0.2923 0.2923 0.2573 0.2245 0.1435 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1955.43328295 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -251.01656887 PAW double counting = 1405587.61088162 -1405037.10628264 entropy T*S EENTRO = -0.00911375 eigenvalues EBANDS = -662.63908036 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 732.39503132 eV energy without entropy = 732.40414507 energy(sigma->0) = 732.39806924 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 2096 total energy-change (2. order) : 0.4208042E+01 (-0.1089830E+01) number of electron 135.9999973 magnetization 4.4022043 augmentation part -9.0416989 magnetization -1.0119421 Broyden mixing: rms(total) = 0.15680E+02 rms(broyden)= 0.15680E+02 rms(prec ) = 0.15973E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.6138 2.4990 1.5202 1.5202 0.8955 0.7155 0.7155 0.5438 0.5438 0.5140 0.5140 0.4287 0.4287 0.3177 0.3177 0.2787 0.2787 0.2221 0.2423 0.2423 0.1435 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1961.89553619 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.82909020 PAW double counting = 1427437.24441589 -1426886.76245686 entropy T*S EENTRO = -0.00754835 eigenvalues EBANDS = -657.13518963 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 736.60307293 eV energy without entropy = 736.61062128 energy(sigma->0) = 736.60558904 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 1792 total energy-change (2. order) : 0.1328472E+02 (-0.9145819E+01) number of electron 135.9999960 magnetization 4.4483582 augmentation part -8.7464171 magnetization -1.1011312 Broyden mixing: rms(total) = 0.14982E+02 rms(broyden)= 0.14982E+02 rms(prec ) = 0.15197E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5997 2.5265 1.5250 1.5250 0.8892 0.7114 0.7114 0.5395 0.5395 0.4684 0.4684 0.4250 0.4250 0.2976 0.2976 0.3288 0.3288 0.2779 0.2779 0.2245 0.2541 0.1435 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1963.82839149 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -242.37400886 PAW double counting = 1456666.77278714 -1456116.20098609 entropy T*S EENTRO = -0.04642535 eigenvalues EBANDS = -645.42365604 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 749.88779757 eV energy without entropy = 749.93422292 energy(sigma->0) = 749.90327269 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 2096 total energy-change (2. order) :-0.7545818E+01 (-0.4663235E+01) number of electron 135.9999983 magnetization 4.3424203 augmentation part -8.7499338 magnetization -0.5424976 Broyden mixing: rms(total) = 0.15705E+02 rms(broyden)= 0.15705E+02 rms(prec ) = 0.15961E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5750 2.5301 1.5299 1.5299 0.8893 0.7106 0.7106 0.5410 0.5410 0.4766 0.4766 0.4271 0.4271 0.2854 0.2854 0.3274 0.3274 0.2755 0.2755 0.2245 0.2530 0.1435 0.0282 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1964.70389858 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -241.34343104 PAW double counting = 1461388.42989714 -1460838.09722821 entropy T*S EENTRO = -0.00574355 eigenvalues EBANDS = -652.92609458 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 742.34197945 eV energy without entropy = 742.34772300 energy(sigma->0) = 742.34389397 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 2248 total energy-change (2. order) :-0.8792537E+04 (-0.8015395E+04) number of electron 135.9999950 magnetization 4.3785797 augmentation part -7.8856500 magnetization 19.6843446 Broyden mixing: rms(total) = 0.17263E+02 rms(broyden)= 0.17222E+02 rms(prec ) = 0.17618E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5561 2.5324 1.5308 1.5308 0.8859 0.7130 0.7130 0.5439 0.5439 0.4829 0.4829 0.4365 0.4365 0.2995 0.2995 0.3269 0.3269 0.2739 0.2739 0.2585 0.2244 0.1435 0.0395 0.0395 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1964.79432031 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.70047448 PAW double counting = 1461127.46703546 -1460577.00420059 entropy T*S EENTRO = 0.01316726 eigenvalues EBANDS = -9443.16445122 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -8050.19476562 eV energy without entropy = -8050.20793287 energy(sigma->0) = -8050.19915470 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 1656 total energy-change (2. order) : 0.8819412E+04 (-0.3350962E+03) number of electron 135.9999977 magnetization 4.3668137 augmentation part -8.1948833 magnetization 16.5057352 Broyden mixing: rms(total) = 0.14899E+02 rms(broyden)= 0.14869E+02 rms(prec ) = 0.15085E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5342 2.5260 1.5250 1.5250 0.8883 0.7155 0.7155 0.5418 0.5418 0.4860 0.4860 0.4365 0.4365 0.2844 0.2844 0.3230 0.3230 0.2740 0.2740 0.2244 0.2611 0.1435 0.0601 0.0601 0.0082 0.0099 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1965.12514656 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -244.51203967 PAW double counting = 1466666.81130658 -1466116.38704884 entropy T*S EENTRO = 0.04240000 eigenvalues EBANDS = -622.60029045 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 769.21765934 eV energy without entropy = 769.17525934 energy(sigma->0) = 769.20352601 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 1888 total energy-change (2. order) :-0.1903142E+02 (-0.1674574E+02) number of electron 135.9999977 magnetization 4.3085895 augmentation part -8.5684800 magnetization 8.6585587 Broyden mixing: rms(total) = 0.15444E+02 rms(broyden)= 0.15442E+02 rms(prec ) = 0.15687E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5213 2.5543 1.5298 1.5298 0.8943 0.7026 0.7026 0.5438 0.5438 0.5034 0.5034 0.4448 0.4448 0.3181 0.3181 0.2949 0.2949 0.2654 0.2245 0.2292 0.2292 0.1435 0.1479 0.1479 0.0178 0.0178 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1965.17701437 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.62885549 PAW double counting = 1469357.53348677 -1468807.23032954 entropy T*S EENTRO = -0.02437558 eigenvalues EBANDS = -642.27515446 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 750.18623561 eV energy without entropy = 750.21061119 energy(sigma->0) = 750.19436080 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 1912 total energy-change (2. order) :-0.9144225E+01 (-0.6566575E+01) number of electron 135.9999965 magnetization 4.4159050 augmentation part -8.7913367 magnetization 4.0765517 Broyden mixing: rms(total) = 0.16112E+02 rms(broyden)= 0.16112E+02 rms(prec ) = 0.16325E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5110 2.5508 1.5424 1.5424 0.8887 0.7092 0.7092 0.5478 0.5478 0.4963 0.4963 0.4452 0.4452 0.3085 0.3085 0.3189 0.3189 0.2620 0.2246 0.2353 0.2353 0.2023 0.2023 0.1435 0.0665 0.0203 0.0203 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1964.77524714 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -244.57027434 PAW double counting = 1472183.69963996 -1471633.50494921 entropy T*S EENTRO = 0.00161219 eigenvalues EBANDS = -650.79724953 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 741.04201020 eV energy without entropy = 741.04039801 energy(sigma->0) = 741.04147281 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 1960 total energy-change (2. order) :-0.1649872E+01 (-0.1674594E+01) number of electron 135.9999965 magnetization 4.4142817 augmentation part -8.8109264 magnetization 1.0754384 Broyden mixing: rms(total) = 0.16170E+02 rms(broyden)= 0.16170E+02 rms(prec ) = 0.16380E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5025 2.5581 1.5329 1.5329 0.8887 0.7072 0.7072 0.5476 0.5476 0.5157 0.5157 0.4475 0.4475 0.2517 0.3114 0.3114 0.3188 0.3188 0.2603 0.2247 0.2362 0.2362 0.1435 0.1753 0.1753 0.1103 0.0201 0.0201 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1965.53026712 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.94459046 PAW double counting = 1480231.62947030 -1479681.33327427 entropy T*S EENTRO = -0.00864072 eigenvalues EBANDS = -652.40903807 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 739.39213792 eV energy without entropy = 739.40077865 energy(sigma->0) = 739.39501817 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 2344 total energy-change (2. order) :-0.3556644E+00 (-0.4583267E+00) number of electron 135.9999976 magnetization 4.0896464 augmentation part -8.8013891 magnetization 0.1504575 Broyden mixing: rms(total) = 0.15934E+02 rms(broyden)= 0.15934E+02 rms(prec ) = 0.16176E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4897 2.5292 1.4739 1.4739 0.8923 0.7008 0.7008 0.5434 0.5434 0.5722 0.5722 0.4478 0.4478 0.2606 0.2606 0.3538 0.3538 0.3219 0.3219 0.2647 0.2647 0.2233 0.1669 0.1434 0.1381 0.0202 0.0202 0.0082 0.0905 0.0905 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1965.52797357 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.87753603 PAW double counting = 1480029.91845790 -1479479.65954809 entropy T*S EENTRO = 0.00888811 eigenvalues EBANDS = -652.81429306 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 739.03647354 eV energy without entropy = 739.02758542 energy(sigma->0) = 739.03351083 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 1880 total energy-change (2. order) : 0.2547558E+01 (-0.8742625E+00) number of electron 135.9999973 magnetization 3.7229993 augmentation part -8.3931817 magnetization -0.7985509 Broyden mixing: rms(total) = 0.15609E+02 rms(broyden)= 0.15609E+02 rms(prec ) = 0.15763E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4868 2.4530 1.2635 1.2635 0.9695 0.7270 0.7270 0.6925 0.6925 0.5105 0.5105 0.4709 0.4709 0.3823 0.3823 0.3723 0.3723 0.3242 0.3242 0.2596 0.2596 0.2657 0.2243 0.1493 0.1435 0.1229 0.1115 0.1115 0.0202 0.0202 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1961.41795479 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -243.68621290 PAW double counting = 1491968.55708636 -1491417.35150770 entropy T*S EENTRO = -0.02769417 eigenvalues EBANDS = -655.47816316 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 741.58403192 eV energy without entropy = 741.61172608 energy(sigma->0) = 741.59326331 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 1832 total energy-change (2. order) :-0.3505892E+01 (-0.8780078E+00) number of electron 135.9999973 magnetization 3.1173532 augmentation part -8.9473686 magnetization -0.8986851 Broyden mixing: rms(total) = 0.14694E+02 rms(broyden)= 0.14694E+02 rms(prec ) = 0.14952E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.4927 2.4173 1.2302 1.2302 0.9619 0.8909 0.8909 0.6811 0.6811 0.4931 0.4931 0.4496 0.4496 0.4751 0.4751 0.3889 0.3889 0.3580 0.3580 0.2876 0.2876 0.2717 0.2242 0.1436 0.1827 0.1435 0.1435 0.1141 0.1141 0.0202 0.0202 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1962.67671573 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -245.72201813 PAW double counting = 1520965.38922918 -1520414.17039173 entropy T*S EENTRO = -0.02260991 eigenvalues EBANDS = -655.70783247 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 738.07813949 eV energy without entropy = 738.10074940 energy(sigma->0) = 738.08567613 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 1816 total energy-change (2. order) :-0.6378699E+01 (-0.1374152E+01) number of electron 135.9999982 magnetization 2.7593087 augmentation part -8.9780123 magnetization -0.9026009 Broyden mixing: rms(total) = 0.13471E+02 rms(broyden)= 0.13471E+02 rms(prec ) = 0.13756E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5109 2.4507 1.3557 1.3557 1.0012 1.0012 0.9491 0.6642 0.6642 0.4594 0.4594 0.4924 0.4924 0.4729 0.4729 0.4592 0.4592 0.3949 0.3949 0.3306 0.3306 0.2777 0.2777 0.2232 0.2124 0.1431 0.1435 0.1346 0.1130 0.1130 0.0202 0.0202 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1958.20393970 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -250.71144448 PAW double counting = 1542117.27245061 -1541562.54804241 entropy T*S EENTRO = 0.00011773 eigenvalues EBANDS = -665.09817942 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 731.69944060 eV energy without entropy = 731.69932286 energy(sigma->0) = 731.69940135 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 2088 total energy-change (2. order) :-0.2506971E+00 (-0.1162222E+01) number of electron 135.9999980 magnetization 2.1718758 augmentation part -9.0089315 magnetization -1.0569378 Broyden mixing: rms(total) = 0.11752E+02 rms(broyden)= 0.11752E+02 rms(prec ) = 0.12044E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5303 2.6188 1.4638 1.4638 1.0937 1.0937 0.9679 0.4513 0.4513 0.6279 0.6279 0.4946 0.4946 0.5608 0.5608 0.4937 0.4937 0.4009 0.4009 0.3698 0.3698 0.2985 0.2879 0.2879 0.2237 0.1429 0.2039 0.1435 0.1354 0.1132 0.1132 0.0202 0.0202 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1956.86099520 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -253.03166634 PAW double counting = 1594124.77286145 -1593573.50322610 entropy T*S EENTRO = -0.00132892 eigenvalues EBANDS = -660.91537970 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 731.44874346 eV energy without entropy = 731.45007238 energy(sigma->0) = 731.44918643 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 1976 total energy-change (2. order) : 0.4267454E+01 (-0.2727925E+00) number of electron 135.9999974 magnetization 1.4042462 augmentation part -8.8721705 magnetization -0.4878319 Broyden mixing: rms(total) = 0.10226E+02 rms(broyden)= 0.10226E+02 rms(prec ) = 0.10473E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.5315 2.6399 1.4300 1.4300 1.1529 1.1529 1.0025 0.6283 0.6283 0.6266 0.6266 0.4989 0.4989 0.4457 0.4457 0.4903 0.4903 0.3969 0.3969 0.3794 0.3794 0.3215 0.3215 0.2808 0.2808 0.2236 0.1429 0.2053 0.1435 0.1353 0.1132 0.1132 0.0202 0.0202 0.0082 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 24.76344702 Ewald energy TEWEN = -729.09930942 -Hartree energ DENC = -1952.96598709 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = -251.60821228 PAW double counting = 1647591.37253149 -1647039.03721255 entropy T*S EENTRO = -0.00333050 eigenvalues EBANDS = -663.03006964 atomic energy EATOM = 3755.32434067 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 735.71619770 eV energy without entropy = 735.71952820 energy(sigma->0) = 735.71730787 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- ----------------------------------------------------------------------------- | | | 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!!! <---- | | | -----------------------------------------------------------------------------