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 2023.12.29 09:39:58 running on 3 total cores distrk: each k-point on 3 cores, 1 groups distr: one band on NCORE= 1 cores, 3 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 = 60 ALGO = Normal (blocked Davidson) ISPIN = 2 MAGMOM = 1.656 1.706 1.667 1.666 1.799 1.141 1.806 1.71 2*1.667 1.799 1.304 1.809 1.714 1.668 1.701 1.801 1.622 1.672 1.707 1.662 1.656 1.799 0.9 1.69 1.715 1.669 1.672 1.801 1.85 1.806 1.71 1.668 1.666 1.799 1.288 1.666 1.702 1.665 1.605 1.798 0.084 1.671 1.706 1.662 1.656 1.798 0.897 1.653 1.707 1.666 1.665 1.799 1.139 0.101 2*0.001 0.059 -0.008 0.06 -0.007 0.014 2*0.001 0 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 1.16 0.73 0.75 0.32 NPAR = 3 POTCAR: PAW_PBE Co 02Aug2007 POTCAR: PAW_PBE O 08Apr2002 POTCAR: PAW_PBE N 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE Co 02Aug2007 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= 2 read in real space projection operators read in non local Contribution for L= 2 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= 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 = 6 number of lm-projection operators is LMMAX = 18 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 N 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 Optimization of the real space projectors (new method) maximal supplied QI-value = 16.25 optimisation between [QCUT,QGAM] = [ 10.24, 20.48] = [ 29.35,117.42] Ry Optimized for a Real-space Cutoff 1.34 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 2 7 10.239 65.824 0.41E-03 0.10E-02 0.11E-06 2 7 10.239 61.429 0.41E-03 0.10E-02 0.11E-06 0 8 10.239 56.174 0.42E-03 0.17E-03 0.96E-07 0 8 10.239 31.263 0.38E-03 0.16E-03 0.87E-07 1 8 10.239 19.828 0.20E-03 0.56E-03 0.13E-06 1 8 10.239 15.139 0.15E-03 0.46E-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 = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.65 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 10 10.053 79.467 0.76E-04 0.72E-04 0.56E-06 0 10 10.053 66.151 0.76E-04 0.72E-04 0.55E-06 1 10 10.053 8.350 0.25E-03 0.92E-03 0.41E-05 1 10 10.053 5.531 0.27E-03 0.10E-02 0.45E-05 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 PAW_PBE Co 02Aug2007 : energy of atom 1 EATOM= -813.3670 kinetic energy error for atom= 0.0093 (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 N 08Apr2002 : energy of atom 3 EATOM= -264.5486 kinetic energy error for atom= 0.0736 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 4 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.119 0.229 0.645- 25 2.37 7 2.42 42 2.43 19 2.47 4 2.48 6 2.49 13 2.53 37 2.54 10 2.54 49 2.58 40 2.65 2 0.226 0.112 0.363- 3 2.47 21 2.47 5 2.47 23 2.47 15 2.47 17 2.47 38 2.47 14 2.47 8 2.47 50 2.47 26 2.47 20 2.47 3 0.115 0.223 0.454- 40 2.45 38 2.47 2 2.47 8 2.47 39 2.47 21 2.47 15 2.47 51 2.47 27 2.47 9 2.47 10 2.59 4 2.62 4 0.230 0.113 0.553- 22 2.37 28 2.44 16 2.46 1 2.48 13 2.49 40 2.51 10 2.53 52 2.54 21 2.60 15 2.62 3 2.62 19 2.77 5 0.115 0.223 0.272- 2 2.47 38 2.47 8 2.47 23 2.47 41 2.47 17 2.47 53 2.47 29 2.47 11 2.47 6 0.230 0.114 0.736- 64 1.59 62 1.76 58 1.81 19 2.32 24 2.33 13 2.47 42 2.48 1 2.49 12 2.51 18 2.58 54 2.61 7 0.117 0.554 0.646- 48 2.29 31 2.41 1 2.42 13 2.47 37 2.50 46 2.52 10 2.52 25 2.57 16 2.66 8 0.226 0.445 0.363- 27 2.47 29 2.47 3 2.47 5 2.47 9 2.47 11 2.47 44 2.47 38 2.47 2 2.47 26 2.47 14 2.47 32 2.47 9 0.115 0.557 0.454- 16 2.46 44 2.47 14 2.47 8 2.47 27 2.47 45 2.47 15 2.47 33 2.47 39 2.47 3 2.47 10 2.53 46 2.59 10 0.231 0.455 0.550- 34 2.33 16 2.44 28 2.47 46 2.50 7 2.52 40 2.53 4 2.53 9 2.53 27 2.53 1 2.54 3 2.59 25 2.69 11 0.115 0.557 0.272- 44 2.47 14 2.47 8 2.47 29 2.47 47 2.47 17 2.47 35 2.47 41 2.47 5 2.47 12 0.261 0.460 0.760- 61 2.11 55 2.20 36 2.21 18 2.43 6 2.51 25 2.66 48 2.77 13 0.122 0.890 0.645- 54 2.36 31 2.40 6 2.47 7 2.47 4 2.49 1 2.53 49 2.57 52 2.58 16 2.61 19 2.69 14 0.226 0.779 0.363- 33 2.47 15 2.47 35 2.47 17 2.47 9 2.47 11 2.47 50 2.47 2 2.47 44 2.47 20 2.47 32 2.47 8 2.47 15 0.115 0.890 0.454- 16 2.44 50 2.47 14 2.47 2 2.47 33 2.47 51 2.47 3 2.47 21 2.47 9 2.47 45 2.47 52 2.57 4 2.62 16 0.224 0.780 0.545- 34 2.43 10 2.44 15 2.44 22 2.44 46 2.45 52 2.46 4 2.46 9 2.46 33 2.46 31 2.47 13 2.61 7 2.66 17 0.115 0.890 0.272- 50 2.47 14 2.47 2 2.47 35 2.47 53 2.47 5 2.47 23 2.47 11 2.47 47 2.47 18 0.225 0.769 0.757- 64 1.04 55 1.32 60 1.72 58 1.75 24 2.37 12 2.43 48 2.50 54 2.55 36 2.57 6 2.58 19 0.466 0.265 0.657- 43 1.97 25 2.05 24 2.26 6 2.32 1 2.47 37 2.50 28 2.63 13 2.69 4 2.77 20 0.559 0.112 0.363- 39 2.47 41 2.47 33 2.47 35 2.47 21 2.47 23 2.47 38 2.47 32 2.47 26 2.47 44 2.47 14 2.47 2 2.47 21 0.448 0.223 0.454- 2 2.47 26 2.47 20 2.47 3 2.47 33 2.47 15 2.47 39 2.47 45 2.47 27 2.47 22 2.53 28 2.54 4 2.60 22 0.549 0.114 0.551- 4 2.37 37 2.37 28 2.42 31 2.43 16 2.44 21 2.53 46 2.54 34 2.54 40 2.55 33 2.57 39 2.61 23 0.448 0.223 0.272- 2 2.47 26 2.47 20 2.47 5 2.47 35 2.47 17 2.47 41 2.47 47 2.47 29 2.47 24 0.528 0.084 0.729- 63 1.22 59 1.71 58 1.77 55 1.93 19 2.26 6 2.33 18 2.37 31 2.41 37 2.60 36 2.61 42 2.63 48 2.71 25 0.451 0.533 0.653- 19 2.05 43 2.06 36 2.23 1 2.37 49 2.50 7 2.57 28 2.60 31 2.64 12 2.66 10 2.69 34 2.73 26 0.559 0.445 0.363- 45 2.47 47 2.47 21 2.47 23 2.47 27 2.47 29 2.47 44 2.47 20 2.47 2 2.47 50 2.47 8 2.47 32 2.47 27 0.448 0.557 0.454- 8 2.47 32 2.47 26 2.47 9 2.47 33 2.47 51 2.47 3 2.47 45 2.47 21 2.47 34 2.53 10 2.53 28 2.55 28 0.557 0.444 0.550- 22 2.42 52 2.44 4 2.44 46 2.46 34 2.46 10 2.47 21 2.54 27 2.55 43 2.55 45 2.56 25 2.60 19 2.63 29 0.448 0.557 0.272- 8 2.47 32 2.47 26 2.47 11 2.47 35 2.47 53 2.47 5 2.47 47 2.47 23 2.47 30 0.586 0.423 0.918- 31 0.439 0.880 0.638- 63 2.08 55 2.13 34 2.38 13 2.40 24 2.41 7 2.41 22 2.43 49 2.47 16 2.47 37 2.52 36 2.63 25 2.64 32 0.559 0.779 0.363- 51 2.47 53 2.47 27 2.47 29 2.47 33 2.47 35 2.47 50 2.47 38 2.47 20 2.47 14 2.47 8 2.47 26 2.47 33 0.448 0.890 0.454- 16 2.46 14 2.47 20 2.47 32 2.47 15 2.47 21 2.47 27 2.47 9 2.47 39 2.47 51 2.47 34 2.55 22 2.57 34 0.546 0.769 0.550- 10 2.33 31 2.38 16 2.43 49 2.45 28 2.46 27 2.53 22 2.54 52 2.54 33 2.55 40 2.56 51 2.61 25 2.73 35 0.448 0.890 0.272- 14 2.47 20 2.47 32 2.47 17 2.47 23 2.47 29 2.47 11 2.47 41 2.47 53 2.47 36 0.551 0.744 0.733- 55 1.36 59 1.69 12 2.21 25 2.23 54 2.45 49 2.50 42 2.52 18 2.57 24 2.61 31 2.63 37 0.775 0.224 0.636- 63 1.49 46 2.37 22 2.37 48 2.47 40 2.50 19 2.50 7 2.50 43 2.51 31 2.52 1 2.54 49 2.54 24 2.60 42 2.61 38 0.893 0.112 0.363- 3 2.47 39 2.47 5 2.47 41 2.47 51 2.47 53 2.47 2 2.47 20 2.47 44 2.47 50 2.47 8 2.47 32 2.47 39 0.782 0.223 0.454- 40 2.44 38 2.47 20 2.47 44 2.47 3 2.47 51 2.47 33 2.47 9 2.47 21 2.47 45 2.47 22 2.61 46 2.62 40 0.892 0.113 0.544- 39 2.44 51 2.45 3 2.45 37 2.50 4 2.51 46 2.53 10 2.53 52 2.53 22 2.55 34 2.56 49 2.59 1 2.65 41 0.782 0.223 0.272- 20 2.47 38 2.47 44 2.47 5 2.47 53 2.47 35 2.47 11 2.47 23 2.47 47 2.47 42 0.885 0.089 0.730- 61 1.64 62 1.75 59 1.80 63 1.95 49 2.40 1 2.43 6 2.48 54 2.48 36 2.52 37 2.61 48 2.63 24 2.63 43 0.729 0.533 0.654- 19 1.97 25 2.06 54 2.18 48 2.29 49 2.45 37 2.51 28 2.55 46 2.76 52 2.77 44 0.893 0.445 0.363- 9 2.47 11 2.47 39 2.47 41 2.47 45 2.47 47 2.47 8 2.47 26 2.47 20 2.47 38 2.47 14 2.47 50 2.47 45 0.782 0.557 0.454- 26 2.47 50 2.47 44 2.47 9 2.47 15 2.47 51 2.47 27 2.47 21 2.47 39 2.47 52 2.52 28 2.56 46 2.56 46 0.893 0.453 0.552- 37 2.37 52 2.37 16 2.45 28 2.46 10 2.50 7 2.52 40 2.53 22 2.54 45 2.56 9 2.59 39 2.62 43 2.76 47 0.782 0.557 0.272- 26 2.47 50 2.47 44 2.47 11 2.47 17 2.47 53 2.47 29 2.47 23 2.47 41 2.47 48 0.895 0.448 0.727- 61 1.56 60 1.84 57 1.90 7 2.29 43 2.29 54 2.32 37 2.47 18 2.50 42 2.63 24 2.71 12 2.77 49 0.772 0.880 0.640- 42 2.40 43 2.45 34 2.45 54 2.46 52 2.47 31 2.47 36 2.50 25 2.50 37 2.54 13 2.57 1 2.58 40 2.59 50 0.893 0.779 0.363- 15 2.47 51 2.47 17 2.47 53 2.47 45 2.47 47 2.47 14 2.47 32 2.47 2 2.47 38 2.47 26 2.47 44 2.47 51 0.782 0.890 0.454- 40 2.45 32 2.47 50 2.47 38 2.47 15 2.47 39 2.47 3 2.47 45 2.47 27 2.47 33 2.47 52 2.57 34 2.61 52 0.889 0.770 0.550- 46 2.37 28 2.44 16 2.46 49 2.47 45 2.52 40 2.53 4 2.54 34 2.54 15 2.57 51 2.57 13 2.58 43 2.77 53 0.782 0.890 0.272- 32 2.47 50 2.47 38 2.47 17 2.47 41 2.47 5 2.47 47 2.47 29 2.47 35 2.47 54 0.887 0.756 0.727- 60 1.80 62 1.90 43 2.18 48 2.32 13 2.36 36 2.45 49 2.46 42 2.48 18 2.55 6 2.61 55 0.392 0.785 0.730- 18 1.32 36 1.36 24 1.93 31 2.13 12 2.20 56 0.678 0.659 0.046- 57 0.708 0.388 0.792- 48 1.90 58 0.348 0.018 0.789- 18 1.75 24 1.77 6 1.81 59 0.641 0.972 0.770- 36 1.69 24 1.71 42 1.80 60 0.976 0.649 0.786- 18 1.72 54 1.80 48 1.84 61 0.934 0.300 0.768- 48 1.56 42 1.64 12 2.11 62 0.001 0.991 0.779- 42 1.75 6 1.76 54 1.90 63 0.662 0.114 0.692- 24 1.22 37 1.49 42 1.95 31 2.08 64 0.189 0.887 0.751- 18 1.04 6 1.59 65 0.312 0.375 0.917- LATTYP: Found a hexagonal cell. ALAT = 7.4168100000 C/A-ratio = 2.9755541803 Lattice vectors: A1 = ( 7.4168100000, 0.0000000000, 0.0000000000) A2 = ( -3.7084050000, 6.4231458800, 0.0000000000) A3 = ( 0.0000000000, 0.0000000000, 22.0691200000) 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 hexagonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 24 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 hexagonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 24 trial point group operations. The dynamic configuration has the point symmetry C_1 . Analysis of constrained symmetry for selective dynamics: ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a hexagonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 24 trial point group operations. The constrained 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 hexagonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 24 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 : 1051.3564 direct lattice vectors reciprocal lattice vectors 7.416810000 0.000000000 0.000000000 0.134828855 0.077843476 0.000000000 -3.708405000 6.423145880 0.000000000 0.000000000 0.155686951 0.000000000 0.000000000 0.000000000 22.069120000 0.000000000 0.000000000 0.045312183 length of vectors 7.416810000 7.416810004 22.069120000 0.155686951 0.155686951 0.045312183 position of ions in fractional coordinates (direct lattice) 0.118848800 0.229224780 0.644504390 0.226119500 0.112066890 0.363033770 0.115008390 0.223178000 0.454180130 0.229552420 0.113181990 0.553277060 0.115008390 0.223178000 0.271887410 0.229866100 0.114019460 0.736021000 0.117145950 0.554421070 0.645825970 0.226119500 0.445400220 0.363033770 0.115008390 0.556511340 0.454180130 0.231006280 0.455299350 0.549953370 0.115008390 0.556511340 0.271887410 0.260963650 0.460065830 0.759824430 0.121964850 0.889709000 0.645344560 0.226119500 0.778733560 0.363033770 0.115008390 0.889844670 0.454180130 0.224404940 0.780209380 0.544696400 0.115008390 0.889844670 0.271887410 0.225078000 0.768822830 0.756702200 0.466384630 0.265253570 0.657396500 0.559452840 0.112066890 0.363033770 0.448341720 0.223178000 0.454180130 0.548726080 0.113584360 0.550905220 0.448341720 0.223178000 0.271887410 0.527760990 0.084048690 0.729101680 0.451381120 0.533023080 0.652857070 0.559452840 0.445400220 0.363033770 0.448341720 0.556511340 0.454180130 0.557455840 0.443677230 0.549934920 0.448341720 0.556511340 0.271887410 0.586343870 0.422867820 0.918247100 0.439273570 0.880255110 0.637658390 0.559452840 0.778733560 0.363033770 0.448341720 0.889844670 0.454180130 0.545696430 0.769381730 0.550499030 0.448341720 0.889844670 0.271887410 0.550649350 0.744136290 0.732774460 0.775196630 0.223915240 0.635962930 0.892786170 0.112066890 0.363033770 0.781675060 0.223178000 0.454180130 0.891797700 0.112951240 0.544275170 0.781675060 0.223178000 0.271887410 0.884592730 0.089110700 0.730447080 0.728947770 0.532997280 0.654238140 0.892786170 0.445400220 0.363033770 0.781675060 0.556511340 0.454180130 0.893234880 0.453357110 0.552149740 0.781675060 0.556511340 0.271887410 0.895140040 0.448084040 0.726793240 0.771645400 0.880201940 0.640183500 0.892786170 0.778733560 0.363033770 0.781675060 0.889844670 0.454180130 0.888995370 0.770397590 0.550116170 0.781675060 0.889844670 0.271887410 0.886609380 0.755962500 0.727277200 0.391665750 0.784681120 0.730149230 0.678059490 0.658707920 0.045749470 0.707922320 0.387995940 0.792312670 0.347906550 0.017989900 0.789313960 0.641075220 0.972147490 0.770432880 0.975996610 0.649237180 0.785623630 0.933740550 0.299754600 0.767922620 0.000837950 0.991101180 0.779414980 0.661605070 0.114300910 0.691999050 0.188632600 0.886974300 0.751396610 0.312022920 0.375224120 0.916597860 ion indices of the primitive-cell ions primitive index ion index 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31 32 32 33 33 34 34 35 35 36 36 37 37 38 38 39 39 40 40 41 41 42 42 43 43 44 44 45 45 46 46 47 47 48 48 49 49 50 50 51 51 52 52 53 53 54 54 55 55 56 56 57 57 58 58 59 59 60 60 61 61 62 62 63 63 64 64 65 65 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 2 2 1 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.067414427 0.038921738 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.077843476 0.000000000 0.000000000 0.500000000 0.000000000 0.000000000 0.000000000 0.045312183 0.000000000 0.000000000 1.000000000 Length of vectors 0.077843476 0.077843476 0.045312183 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.500000 0.000000 1.000000 0.500000 0.500000 0.000000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 0.067414 0.038922 0.000000 1.000000 0.000000 0.077843 0.000000 1.000000 0.067414 0.116765 0.000000 1.000000 -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 4 k-points in BZ NKDIM = 4 number of bands NBANDS= 369 number of dos NEDOS = 301 number of ions NIONS = 65 non local maximal LDIM = 6 non local SUM 2l+1 LMDIM = 18 total plane-waves NPLWV = 139968 max r-space proj IRMAX = 2547 max aug-charges IRDMAX= 5449 dimension x,y,z NGX = 36 NGY = 36 NGZ = 108 dimension x,y,z NGXF= 72 NGYF= 72 NGZF= 216 support grid NGXF= 72 NGYF= 72 NGZF= 216 ions per type = 54 8 1 2 NGX,Y,Z is equivalent to a cutoff of 8.07, 8.07, 8.14 a.u. NGXF,Y,Z is equivalent to a cutoff of 16.14, 16.14, 16.27 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.09 12.09 35.99*2*pi/ulx,y,z ENINI = 400.0 initial cutoff ENAUG = 627.1 eV augmentation charge cutoff NELM = 60; NELMIN= 2; NELMDL= -5 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM LREAL = T real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 1 Vosko Wilk Nusair interpolation ROPT = -0.00050 -0.00050 -0.00050 -0.00050 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.126E-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 = 58.93 16.00 14.00 1.00 Ionic Valenz ZVAL = 9.00 6.00 5.00 1.00 Atomic Wigner-Seitz radii RWIGS = 1.16 0.73 0.75 0.32 virtual crystal weights VCA = 1.00 1.00 1.00 1.00 NELECT = 541.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.68E-08 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 16.17 109.15 Fermi-wavevector in a.u.,A,eV,Ry = 1.311869 2.479074 23.415644 1.721001 Thomas-Fermi vector in A = 2.442301 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 98 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 : 1051.36 direct lattice vectors reciprocal lattice vectors 7.416810000 0.000000000 0.000000000 0.134828855 0.077843476 0.000000000 -3.708405000 6.423145880 0.000000000 0.000000000 0.155686951 0.000000000 0.000000000 0.000000000 22.069120000 0.000000000 0.000000000 0.045312183 length of vectors 7.416810000 7.416810004 22.069120000 0.155686951 0.155686951 0.045312183 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 0.250 0.06741443 0.03892174 0.00000000 0.250 0.00000000 0.07784348 0.00000000 0.250 0.06741443 0.11676521 0.00000000 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.50000000 0.00000000 0.250 0.50000000 0.50000000 0.00000000 0.250 position of ions in fractional coordinates (direct lattice) 0.11884880 0.22922478 0.64450439 0.22611950 0.11206689 0.36303377 0.11500839 0.22317800 0.45418013 0.22955242 0.11318199 0.55327706 0.11500839 0.22317800 0.27188741 0.22986610 0.11401946 0.73602100 0.11714595 0.55442107 0.64582597 0.22611950 0.44540022 0.36303377 0.11500839 0.55651134 0.45418013 0.23100628 0.45529935 0.54995337 0.11500839 0.55651134 0.27188741 0.26096365 0.46006583 0.75982443 0.12196485 0.88970900 0.64534456 0.22611950 0.77873356 0.36303377 0.11500839 0.88984467 0.45418013 0.22440494 0.78020938 0.54469640 0.11500839 0.88984467 0.27188741 0.22507800 0.76882283 0.75670220 0.46638463 0.26525357 0.65739650 0.55945284 0.11206689 0.36303377 0.44834172 0.22317800 0.45418013 0.54872608 0.11358436 0.55090522 0.44834172 0.22317800 0.27188741 0.52776099 0.08404869 0.72910168 0.45138112 0.53302308 0.65285707 0.55945284 0.44540022 0.36303377 0.44834172 0.55651134 0.45418013 0.55745584 0.44367723 0.54993492 0.44834172 0.55651134 0.27188741 0.58634387 0.42286782 0.91824710 0.43927357 0.88025511 0.63765839 0.55945284 0.77873356 0.36303377 0.44834172 0.88984467 0.45418013 0.54569643 0.76938173 0.55049903 0.44834172 0.88984467 0.27188741 0.55064935 0.74413629 0.73277446 0.77519663 0.22391524 0.63596293 0.89278617 0.11206689 0.36303377 0.78167506 0.22317800 0.45418013 0.89179770 0.11295124 0.54427517 0.78167506 0.22317800 0.27188741 0.88459273 0.08911070 0.73044708 0.72894777 0.53299728 0.65423814 0.89278617 0.44540022 0.36303377 0.78167506 0.55651134 0.45418013 0.89323488 0.45335711 0.55214974 0.78167506 0.55651134 0.27188741 0.89514004 0.44808404 0.72679324 0.77164540 0.88020194 0.64018350 0.89278617 0.77873356 0.36303377 0.78167506 0.88984467 0.45418013 0.88899537 0.77039759 0.55011617 0.78167506 0.88984467 0.27188741 0.88660938 0.75596250 0.72727720 0.39166575 0.78468112 0.73014923 0.67805949 0.65870792 0.04574947 0.70792232 0.38799594 0.79231267 0.34790655 0.01798990 0.78931396 0.64107522 0.97214749 0.77043288 0.97599661 0.64923718 0.78562363 0.93374055 0.29975460 0.76792262 0.00083795 0.99110118 0.77941498 0.66160507 0.11430091 0.69199905 0.18863260 0.88697430 0.75139661 0.31202292 0.37522412 0.91659786 position of ions in cartesian coordinates (Angst): 0.03142065 1.47234420 14.22364472 1.26149595 0.71982198 8.01183583 0.02536097 1.43350485 10.02335579 1.28282203 0.72698443 12.21033783 0.02536097 1.43350485 6.00031588 1.28204285 0.73236362 16.24333577 -1.18716861 3.56112741 14.25281083 0.02536097 2.86087059 8.01183583 -1.21077406 3.57455352 10.02335579 0.02489530 2.92445414 12.13698692 -1.21077406 3.57455352 6.00031588 0.22940738 2.95506994 16.76865652 -2.39481119 5.71473070 14.24218654 -1.21077406 5.00191926 8.01183583 -2.44690905 5.71560213 10.02335579 -1.22896356 5.01139866 12.02097022 -2.44690905 5.71560213 6.00031588 -1.18174567 4.93826119 16.69975166 2.47541852 1.70376238 14.50816225 3.73376600 0.71982198 8.01183583 2.49763094 1.43350485 10.02335579 3.64858027 0.72956891 12.15799341 2.49763094 1.43350485 6.00031588 3.60261641 0.53985700 16.09063247 1.37114255 3.42368500 14.40798102 2.49763102 2.86087059 8.01183583 1.26149592 3.57455352 10.02335579 2.48920919 2.84980357 12.13657974 1.26149592 3.57455352 6.00031588 2.78063594 2.71614170 20.26490544 -0.00633384 5.65400698 14.07255953 1.26149599 5.00191926 8.01183583 0.02536093 5.71560213 10.02335579 1.19414768 4.94185109 12.14902915 0.02536093 5.71560213 6.00031588 1.32450287 4.77969595 16.17168749 4.91911772 1.43824025 14.03514222 6.20603598 0.71982198 8.01183583 4.96990099 1.43350485 10.02335579 6.19542516 0.72550229 12.01167404 4.96990099 1.43350485 6.00031588 6.23039764 0.57237103 16.12032426 3.42989733 3.42351928 14.43846002 4.96990099 2.86087059 8.01183583 3.73376597 3.57455352 10.02335579 4.94372162 2.91197885 12.18545887 3.73376597 3.57455352 6.00031588 4.97740651 2.87810916 16.03968723 2.45900204 5.65366546 14.12828648 3.73376597 5.00191926 8.01183583 2.49763098 5.71560213 10.02335579 3.73656348 4.94837611 12.14057977 2.49763098 5.71560213 6.00031588 3.77239820 4.85565742 16.05036780 -0.00500494 5.04012130 16.11375097 2.58628266 4.23097706 1.00965054 3.81167926 2.49215452 17.48564339 2.51364294 0.11555175 17.41946450 1.14961649 6.24424515 17.00277568 4.83114701 4.17014512 17.33802217 5.81376479 1.92536752 16.94737645 -3.66918966 6.36598746 17.20100272 4.48312503 0.73417142 15.27181007 -1.89020777 5.69716532 16.58266195 0.92273171 2.41011926 20.22850816 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 18983 k-point 2 : 0.5000-0.0000 0.0000 plane waves: 19132 k-point 3 : 0.0000 0.5000 0.0000 plane waves: 19132 k-point 4 : 0.5000 0.5000 0.0000 plane waves: 19132 maximum and minimum number of plane-waves per node : 19132 18983 maximum number of plane-waves: 19132 maximum index in each direction: IXMAX= 12 IYMAX= 12 IZMAX= 35 IXMIN= -12 IYMIN= -12 IZMIN= -35 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 50 to avoid them WARNING: aliasing errors must be expected set NGZ to 144 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 425456. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 15029. kBytes fftplans : 13188. kBytes grid : 47347. kBytes one-center: 2021. kBytes wavefun : 317871. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 25 NGY = 25 NGZ = 71 (NGX = 72 NGY = 72 NGZ =216) gives a total of 44375 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 541.0000000 magnetization 87.3950000 keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 2419 Maximum index for augmentation-charges 1716 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.174 Maximum number of real-space cells 4x 4x 2 Maximum number of reciprocal cells 2x 2x 5 ----------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 5913 total energy-change (2. order) : 0.4760465E+04 (-0.1985385E+05) number of electron 541.0000000 magnetization 87.3950000 augmentation part 541.0000000 magnetization 87.3950000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154972.11165768 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1166.35087493 PAW double counting = 41899.85718576 -44777.76369754 entropy T*S EENTRO = -0.01735667 eigenvalues EBANDS = 5945.17777469 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 4760.46474155 eV energy without entropy = 4760.48209822 energy(sigma->0) = 4760.47052711 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 5904 total energy-change (2. order) :-0.4387877E+04 (-0.4122881E+04) number of electron 541.0000000 magnetization 87.3950000 augmentation part 541.0000000 magnetization 87.3950000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154972.11165768 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1166.35087493 PAW double counting = 41899.85718576 -44777.76369754 entropy T*S EENTRO = -0.01366074 eigenvalues EBANDS = 1557.29661640 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 372.58727920 eV energy without entropy = 372.60093994 energy(sigma->0) = 372.59183278 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 6696 total energy-change (2. order) :-0.7250089E+03 (-0.6013191E+03) number of electron 541.0000000 magnetization 87.3950000 augmentation part 541.0000000 magnetization 87.3950000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154972.11165768 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1166.35087493 PAW double counting = 41899.85718576 -44777.76369754 entropy T*S EENTRO = 0.11228632 eigenvalues EBANDS = 832.16179157 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -352.42159858 eV energy without entropy = -352.53388490 energy(sigma->0) = -352.45902735 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 7632 total energy-change (2. order) :-0.1148412E+03 (-0.1063069E+03) number of electron 541.0000000 magnetization 87.3950000 augmentation part 541.0000000 magnetization 87.3950000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154972.11165768 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1166.35087493 PAW double counting = 41899.85718576 -44777.76369754 entropy T*S EENTRO = 0.20939315 eigenvalues EBANDS = 717.22346694 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -467.26281638 eV energy without entropy = -467.47220952 energy(sigma->0) = -467.33261409 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 6417 total energy-change (2. order) :-0.4296452E+01 (-0.4256403E+01) number of electron 540.9999690 magnetization 85.9182136 augmentation part 253.0895079 magnetization 69.5397112 Broyden mixing: rms(total) = 0.13522E+02 rms(broyden)= 0.13521E+02 rms(prec ) = 0.14194E+02 weight for this iteration 100.00 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154972.11165768 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1166.35087493 PAW double counting = 41899.85718576 -44777.76369754 entropy T*S EENTRO = 0.20843757 eigenvalues EBANDS = 712.92797101 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -471.55926789 eV energy without entropy = -471.76770546 energy(sigma->0) = -471.62874708 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 7920 total energy-change (2. order) :-0.4069428E+03 (-0.8418783E+02) number of electron 540.9999824 magnetization 85.7495656 augmentation part 270.5943008 magnetization 74.3064561 Broyden mixing: rms(total) = 0.71618E+02 rms(broyden)= 0.71617E+02 rms(prec ) = 0.73520E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0311 0.0311 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155020.45010970 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1131.42697441 PAW double counting = 47585.37811348 -50181.79952900 entropy T*S EENTRO = 0.09468596 eigenvalues EBANDS = 107.87617098 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -878.50207580 eV energy without entropy = -878.59676176 energy(sigma->0) = -878.53363778 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 7992 total energy-change (2. order) : 0.1027501E+03 (-0.5202463E+02) number of electron 540.9999697 magnetization 85.8018292 augmentation part 285.1928520 magnetization 73.7283597 Broyden mixing: rms(total) = 0.84747E+02 rms(broyden)= 0.84746E+02 rms(prec ) = 0.85313E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0321 0.0321 0.0321 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154856.21093117 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1130.80134479 PAW double counting = 48849.61342909 -51448.26163710 entropy T*S EENTRO = -0.20049885 eigenvalues EBANDS = 49.53468222 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -775.75199296 eV energy without entropy = -775.55149411 energy(sigma->0) = -775.68516001 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 7074 total energy-change (2. order) :-0.3441735E+02 (-0.5026600E+01) number of electron 540.9999780 magnetization 85.8362540 augmentation part 282.5097441 magnetization 73.9009160 Broyden mixing: rms(total) = 0.85312E+02 rms(broyden)= 0.85312E+02 rms(prec ) = 0.86362E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0741 0.0630 0.0797 0.0797 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154902.17956535 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1131.87723366 PAW double counting = 49061.28140134 -51663.29499090 entropy T*S EENTRO = 0.03549829 eigenvalues EBANDS = 63.13946189 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -810.16934300 eV energy without entropy = -810.20484129 energy(sigma->0) = -810.18117577 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 8649 total energy-change (2. order) :-0.2851680E+03 (-0.7086541E+01) number of electron 541.0000038 magnetization 86.3524335 augmentation part 278.9567767 magnetization 79.6174579 Broyden mixing: rms(total) = 0.87736E+02 rms(broyden)= 0.87736E+02 rms(prec ) = 0.91541E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0861 0.0831 0.0831 0.0892 0.0892 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154800.81064069 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1131.77044580 PAW double counting = 48757.02255727 -51354.54770840 entropy T*S EENTRO = -0.09517732 eigenvalues EBANDS = -327.64840986 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1095.33731513 eV energy without entropy = -1095.24213781 energy(sigma->0) = -1095.30558936 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 9414 total energy-change (2. order) :-0.1018060E+03 (-0.5069653E+01) number of electron 541.0000120 magnetization 86.1162811 augmentation part 281.5705899 magnetization 78.7360964 Broyden mixing: rms(total) = 0.87693E+02 rms(broyden)= 0.87693E+02 rms(prec ) = 0.91938E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.1103 0.1227 0.1283 0.1283 0.0998 0.0723 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154749.85528874 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1134.62616342 PAW double counting = 48721.47265808 -51315.37614820 entropy T*S EENTRO = -0.09190687 eigenvalues EBANDS = -486.89038510 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1197.14328933 eV energy without entropy = -1197.05138246 energy(sigma->0) = -1197.11265371 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 8361 total energy-change (2. order) :-0.3484578E+02 (-0.3085983E+01) number of electron 541.0000255 magnetization 85.8395006 augmentation part 282.0509086 magnetization 76.8012651 Broyden mixing: rms(total) = 0.87769E+02 rms(broyden)= 0.87769E+02 rms(prec ) = 0.92576E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.0982 0.1309 0.1354 0.1354 0.1008 0.0734 0.0135 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154669.57928012 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1133.21578997 PAW double counting = 48723.88031768 -51317.42605316 entropy T*S EENTRO = 0.00028723 eigenvalues EBANDS = -601.05175150 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1231.98907183 eV energy without entropy = -1231.98935906 energy(sigma->0) = -1231.98916757 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 6984 total energy-change (2. order) : 0.9294029E+01 (-0.1397167E+01) number of electron 541.0000218 magnetization 82.4866476 augmentation part 283.6372866 magnetization 73.7172305 Broyden mixing: rms(total) = 0.87561E+02 rms(broyden)= 0.87561E+02 rms(prec ) = 0.92027E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.1657 0.2148 0.2511 0.2511 0.1622 0.1622 0.0591 0.0591 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154695.31024086 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1131.78415950 PAW double counting = 48718.18468741 -51310.88404011 entropy T*S EENTRO = -0.09656636 eigenvalues EBANDS = -565.34466039 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1222.69504274 eV energy without entropy = -1222.59847638 energy(sigma->0) = -1222.66285395 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 8784 total energy-change (2. order) :-0.6577101E+02 (-0.1520316E+02) number of electron 541.0000060 magnetization 82.0222800 augmentation part 267.5139748 magnetization 69.3548738 Broyden mixing: rms(total) = 0.90479E+02 rms(broyden)= 0.90478E+02 rms(prec ) = 0.96033E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.1627 0.2348 0.2348 0.1874 0.1874 0.1704 0.1704 0.0580 0.0580 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154739.47601262 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1136.01775976 PAW double counting = 48145.53618627 -50772.99291226 entropy T*S EENTRO = -0.01963238 eigenvalues EBANDS = -556.50305709 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1288.46605024 eV energy without entropy = -1288.44641786 energy(sigma->0) = -1288.45950612 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 7587 total energy-change (2. order) : 0.2423969E+02 (-0.5410951E+01) number of electron 541.0000126 magnetization 80.1488545 augmentation part 266.5108544 magnetization 66.1476089 Broyden mixing: rms(total) = 0.90199E+02 rms(broyden)= 0.90199E+02 rms(prec ) = 0.95690E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2062 0.3639 0.3639 0.2572 0.2176 0.2176 0.1602 0.1602 0.0577 0.0577 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154732.39013326 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1136.63911744 PAW double counting = 48120.18676819 -50748.93127987 entropy T*S EENTRO = -0.02043529 eigenvalues EBANDS = -538.68201637 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1264.22636109 eV energy without entropy = -1264.20592580 energy(sigma->0) = -1264.21954933 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 9072 total energy-change (2. order) :-0.7647923E+02 (-0.8163061E+01) number of electron 541.0000015 magnetization 77.7384286 augmentation part 274.5438432 magnetization 65.6955607 Broyden mixing: rms(total) = 0.90102E+02 rms(broyden)= 0.90102E+02 rms(prec ) = 0.95376E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2199 0.4484 0.4484 0.2686 0.2482 0.2482 0.1715 0.1715 0.0577 0.0577 0.0790 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154559.60543461 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1136.84875299 PAW double counting = 48003.70943876 -50606.76322436 entropy T*S EENTRO = -0.00093057 eigenvalues EBANDS = -813.86581327 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1340.70559298 eV energy without entropy = -1340.70466242 energy(sigma->0) = -1340.70528280 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 8964 total energy-change (2. order) : 0.1758246E+02 (-0.1008232E+02) number of electron 541.0000015 magnetization 77.5536567 augmentation part 274.5485493 magnetization 64.5429504 Broyden mixing: rms(total) = 0.90214E+02 rms(broyden)= 0.90214E+02 rms(prec ) = 0.95200E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2000 0.4417 0.4417 0.2678 0.2466 0.2466 0.1712 0.1712 0.0576 0.0576 0.0813 0.0165 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154478.77462983 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1138.46871427 PAW double counting = 48144.38445193 -50736.13510688 entropy T*S EENTRO = -0.01213210 eigenvalues EBANDS = -890.02604370 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1323.12312825 eV energy without entropy = -1323.11099615 energy(sigma->0) = -1323.11908421 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 7047 total energy-change (2. order) :-0.1158532E+01 (-0.1707370E+01) number of electron 541.0000019 magnetization 74.9205233 augmentation part 274.2277944 magnetization 61.7339660 Broyden mixing: rms(total) = 0.90333E+02 rms(broyden)= 0.90333E+02 rms(prec ) = 0.95339E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2143 0.4935 0.4935 0.2736 0.2501 0.2501 0.1678 0.1678 0.1354 0.1354 0.0576 0.0576 0.0894 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154472.89664654 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1138.82542739 PAW double counting = 48173.82076537 -50765.04044937 entropy T*S EENTRO = 0.00479612 eigenvalues EBANDS = -897.96717079 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1324.28165976 eV energy without entropy = -1324.28645588 energy(sigma->0) = -1324.28325846 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 8127 total energy-change (2. order) :-0.7769928E+02 (-0.4622802E+01) number of electron 540.9999985 magnetization 74.5935140 augmentation part 274.3009216 magnetization 62.1357388 Broyden mixing: rms(total) = 0.91696E+02 rms(broyden)= 0.91696E+02 rms(prec ) = 0.96974E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.1979 0.4919 0.4919 0.2734 0.2506 0.2506 0.1679 0.1679 0.1328 0.1328 0.0576 0.0576 0.0896 0.0080 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154396.59996505 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1146.45067105 PAW double counting = 48692.06951016 -51278.31096588 entropy T*S EENTRO = 0.01308634 eigenvalues EBANDS = -1064.57488983 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1401.98093514 eV energy without entropy = -1401.99402148 energy(sigma->0) = -1401.98529725 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 7362 total energy-change (2. order) : 0.1069126E+01 (-0.1719558E+01) number of electron 540.9999992 magnetization 73.9432059 augmentation part 273.7362698 magnetization 60.8766230 Broyden mixing: rms(total) = 0.91735E+02 rms(broyden)= 0.91735E+02 rms(prec ) = 0.97053E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2272 0.5470 0.5470 0.2771 0.2892 0.2892 0.2335 0.2335 0.1756 0.1756 0.0576 0.0576 0.1025 0.1025 0.0929 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154385.54116039 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1147.34786808 PAW double counting = 48796.01932491 -51381.99509760 entropy T*S EENTRO = -0.00258816 eigenvalues EBANDS = -1075.71177434 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1400.91180943 eV energy without entropy = -1400.90922127 energy(sigma->0) = -1400.91094671 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 8190 total energy-change (2. order) :-0.5107274E+02 (-0.5595020E+01) number of electron 540.9999960 magnetization 72.1376465 augmentation part 278.3188197 magnetization 59.7841346 Broyden mixing: rms(total) = 0.92182E+02 rms(broyden)= 0.92182E+02 rms(prec ) = 0.97422E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2462 0.6554 0.6554 0.2784 0.3594 0.3594 0.2502 0.2502 0.1786 0.1786 0.1198 0.1198 0.0576 0.0576 0.0861 0.0861 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154349.69715775 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1147.69563493 PAW double counting = 49041.67925769 -51624.13510568 entropy T*S EENTRO = 0.01548926 eigenvalues EBANDS = -1166.51428258 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1451.98454606 eV energy without entropy = -1452.00003531 energy(sigma->0) = -1451.98970914 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 8289 total energy-change (2. order) : 0.1094152E+02 (-0.7710472E+01) number of electron 540.9999975 magnetization 71.0893353 augmentation part 276.7076361 magnetization 57.4835802 Broyden mixing: rms(total) = 0.93147E+02 rms(broyden)= 0.93147E+02 rms(prec ) = 0.98219E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2341 0.6659 0.6659 0.2785 0.3636 0.3636 0.2511 0.2511 0.1777 0.1777 0.1209 0.1209 0.0576 0.0576 0.0848 0.0848 0.0234 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154345.65736294 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.66565564 PAW double counting = 50128.48771320 -52717.65887123 entropy T*S EENTRO = -0.00510906 eigenvalues EBANDS = -1155.84666683 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1441.04302314 eV energy without entropy = -1441.03791408 energy(sigma->0) = -1441.04132012 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 8415 total energy-change (2. order) : 0.1760832E+01 (-0.2380517E+01) number of electron 540.9999958 magnetization 71.6114316 augmentation part 276.0722066 magnetization 57.7795611 Broyden mixing: rms(total) = 0.93626E+02 rms(broyden)= 0.93626E+02 rms(prec ) = 0.98678E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2387 0.6125 0.6125 0.3550 0.3550 0.3321 0.2819 0.2530 0.2530 0.1785 0.1785 0.1315 0.1315 0.0576 0.0576 0.0884 0.0894 0.0894 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154304.53346167 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1154.14419672 PAW double counting = 51053.73821400 -53642.34414621 entropy T*S EENTRO = -0.03322204 eigenvalues EBANDS = -1199.22538952 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1439.28219066 eV energy without entropy = -1439.24896862 energy(sigma->0) = -1439.27111665 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 8262 total energy-change (2. order) :-0.2808450E+01 (-0.1325611E+01) number of electron 540.9999966 magnetization 72.6481080 augmentation part 276.9569889 magnetization 58.4117516 Broyden mixing: rms(total) = 0.94483E+02 rms(broyden)= 0.94483E+02 rms(prec ) = 0.99451E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2471 0.5940 0.6206 0.6206 0.2794 0.3627 0.3627 0.2561 0.2561 0.1799 0.1799 0.1376 0.1376 0.0576 0.0576 0.0915 0.0915 0.0846 0.0771 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154327.49266349 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1152.39168576 PAW double counting = 50529.74183835 -53119.92854877 entropy T*S EENTRO = -0.01510045 eigenvalues EBANDS = -1175.75946989 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1442.09064042 eV energy without entropy = -1442.07553996 energy(sigma->0) = -1442.08560693 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 7848 total energy-change (2. order) : 0.5862737E+02 (-0.2110321E+01) number of electron 540.9999958 magnetization 73.5380634 augmentation part 276.8392175 magnetization 60.1422817 Broyden mixing: rms(total) = 0.91985E+02 rms(broyden)= 0.91985E+02 rms(prec ) = 0.96805E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2539 0.7502 0.6490 0.6490 0.2793 0.3757 0.3757 0.2595 0.2595 0.1748 0.1748 0.1454 0.1454 0.1091 0.1091 0.0576 0.0576 0.0866 0.0866 0.0794 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154387.53752223 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.81248641 PAW double counting = 49668.67138419 -52265.41755585 entropy T*S EENTRO = 0.01177079 eigenvalues EBANDS = -1048.97544886 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1383.46326749 eV energy without entropy = -1383.47503827 energy(sigma->0) = -1383.46719108 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 8289 total energy-change (2. order) : 0.4782887E+02 (-0.2943042E+01) number of electron 541.0000020 magnetization 74.3328531 augmentation part 276.1124249 magnetization 60.6007816 Broyden mixing: rms(total) = 0.92052E+02 rms(broyden)= 0.92052E+02 rms(prec ) = 0.96835E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2650 0.9130 0.6879 0.6879 0.2792 0.3931 0.3931 0.2626 0.2626 0.1830 0.1830 0.1744 0.1744 0.1224 0.1224 0.0576 0.0576 0.0929 0.0929 0.0887 0.0713 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154467.56792172 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.18794002 PAW double counting = 48970.54003894 -51575.23542240 entropy T*S EENTRO = 0.01471815 eigenvalues EBANDS = -912.54536652 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1335.63439545 eV energy without entropy = -1335.64911360 energy(sigma->0) = -1335.63930150 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 8334 total energy-change (2. order) : 0.5831829E+02 (-0.3739299E+01) number of electron 541.0000046 magnetization 75.2026222 augmentation part 277.2300650 magnetization 61.2813192 Broyden mixing: rms(total) = 0.90634E+02 rms(broyden)= 0.90634E+02 rms(prec ) = 0.95233E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2758 1.0680 0.7284 0.7284 0.4118 0.4118 0.2792 0.2653 0.2653 0.2218 0.2218 0.1756 0.1756 0.1364 0.1364 0.0576 0.0576 0.1009 0.1009 0.0883 0.0883 0.0722 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154550.31668787 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.06600548 PAW double counting = 48401.04675763 -51012.57330470 entropy T*S EENTRO = 0.02839926 eigenvalues EBANDS = -764.53888944 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1277.31610156 eV energy without entropy = -1277.34450082 energy(sigma->0) = -1277.32556798 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 8442 total energy-change (2. order) : 0.8226232E+02 (-0.2822443E+01) number of electron 541.0000013 magnetization 75.4928689 augmentation part 275.9294549 magnetization 61.2467895 Broyden mixing: rms(total) = 0.89745E+02 rms(broyden)= 0.89745E+02 rms(prec ) = 0.94027E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2907 1.2500 0.7830 0.7830 0.4391 0.4391 0.2792 0.2919 0.2919 0.2483 0.2483 0.1861 0.1861 0.1480 0.1480 0.1090 0.1090 0.0576 0.0576 0.0903 0.0903 0.0882 0.0717 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154603.96739870 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1149.86161051 PAW double counting = 48155.93073856 -50772.44943232 entropy T*S EENTRO = 0.01282814 eigenvalues EBANDS = -623.41374980 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1195.05378554 eV energy without entropy = -1195.06661367 energy(sigma->0) = -1195.05806158 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 28) --------------------------------------- eigenvalue-minimisations : 7686 total energy-change (2. order) : 0.1882484E+02 (-0.1281355E+01) number of electron 541.0000027 magnetization 75.6911519 augmentation part 274.7861835 magnetization 61.1163134 Broyden mixing: rms(total) = 0.89387E+02 rms(broyden)= 0.89387E+02 rms(prec ) = 0.93675E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3037 1.4673 0.8422 0.8422 0.4696 0.4696 0.2792 0.3170 0.3170 0.2585 0.2585 0.1954 0.1954 0.1525 0.1525 0.1144 0.1144 0.0576 0.0576 0.0925 0.0925 0.0843 0.0835 0.0714 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154628.67934485 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.09086280 PAW double counting = 48181.42313641 -50798.51844606 entropy T*S EENTRO = 0.02912628 eigenvalues EBANDS = -579.54590072 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1176.22894807 eV energy without entropy = -1176.25807435 energy(sigma->0) = -1176.23865683 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 29) --------------------------------------- eigenvalue-minimisations : 7740 total energy-change (2. order) : 0.1059930E+02 (-0.2592390E+00) number of electron 541.0000026 magnetization 75.1126697 augmentation part 274.5128684 magnetization 60.5136492 Broyden mixing: rms(total) = 0.89168E+02 rms(broyden)= 0.89168E+02 rms(prec ) = 0.93424E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.2897 1.3922 0.8207 0.8207 0.4623 0.4623 0.2792 0.3128 0.3128 0.2577 0.2577 0.1978 0.1978 0.1064 0.1511 0.1511 0.1146 0.1146 0.0576 0.0576 0.0934 0.0934 0.0885 0.0796 0.0716 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154647.31226425 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1150.36467334 PAW double counting = 48235.29433060 -50851.53487886 entropy T*S EENTRO = 0.02218773 eigenvalues EBANDS = -551.43531326 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1165.62964662 eV energy without entropy = -1165.65183435 energy(sigma->0) = -1165.63704253 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 30) --------------------------------------- eigenvalue-minimisations : 8631 total energy-change (2. order) :-0.3935635E+02 (-0.1185794E+01) number of electron 541.0000012 magnetization 73.0646681 augmentation part 274.9865264 magnetization 58.8527521 Broyden mixing: rms(total) = 0.90120E+02 rms(broyden)= 0.90120E+02 rms(prec ) = 0.94396E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3074 1.0650 1.0650 0.7726 0.7726 0.4728 0.4728 0.2792 0.3264 0.3264 0.2677 0.2677 0.2146 0.2146 0.1513 0.1513 0.1188 0.1188 0.0576 0.0576 0.0975 0.0975 0.0887 0.0887 0.0717 0.0686 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154586.76643908 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1149.84849610 PAW double counting = 48208.71785329 -50829.74538515 entropy T*S EENTRO = -0.00734874 eigenvalues EBANDS = -646.00478708 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1204.98599258 eV energy without entropy = -1204.97864384 energy(sigma->0) = -1204.98354300 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 31) --------------------------------------- eigenvalue-minimisations : 8550 total energy-change (2. order) :-0.7616344E+02 (-0.2519368E+01) number of electron 541.0000017 magnetization 71.0549750 augmentation part 275.0891756 magnetization 56.2984473 Broyden mixing: rms(total) = 0.91238E+02 rms(broyden)= 0.91238E+02 rms(prec ) = 0.95708E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3303 1.3644 1.3644 0.8014 0.8014 0.4946 0.4946 0.2792 0.3418 0.3418 0.2734 0.2734 0.2217 0.2217 0.1612 0.1612 0.1255 0.1255 0.0576 0.0576 0.1048 0.1048 0.0917 0.0917 0.0885 0.0716 0.0716 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154452.88972909 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1154.72115449 PAW double counting = 49670.60241592 -52305.07375696 entropy T*S EENTRO = -0.01150113 eigenvalues EBANDS = -847.46962957 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1281.14942827 eV energy without entropy = -1281.13792714 energy(sigma->0) = -1281.14559456 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 32) --------------------------------------- eigenvalue-minimisations : 8550 total energy-change (2. order) :-0.5971983E+02 (-0.2918890E+01) number of electron 541.0000026 magnetization 72.3429444 augmentation part 273.3962573 magnetization 57.5321127 Broyden mixing: rms(total) = 0.92656E+02 rms(broyden)= 0.92656E+02 rms(prec ) = 0.97257E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3128 1.2806 1.2806 0.7933 0.7933 0.4900 0.4900 0.2792 0.3386 0.3386 0.2730 0.2730 0.2219 0.2219 0.1592 0.1592 0.0633 0.1249 0.1249 0.1051 0.1051 0.0576 0.0576 0.0919 0.0919 0.0884 0.0715 0.0715 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154362.20003097 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1161.91479179 PAW double counting = 52684.28088945 -55330.91729725 entropy T*S EENTRO = -0.00765818 eigenvalues EBANDS = -992.91157534 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1340.86926244 eV energy without entropy = -1340.86160426 energy(sigma->0) = -1340.86670972 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 33) --------------------------------------- eigenvalue-minimisations : 8262 total energy-change (2. order) : 0.4867181E+02 (-0.1387104E+01) number of electron 541.0000001 magnetization 72.8769341 augmentation part 273.8963922 magnetization 57.8253375 Broyden mixing: rms(total) = 0.91505E+02 rms(broyden)= 0.91505E+02 rms(prec ) = 0.95990E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3043 1.2946 1.2946 0.7953 0.7953 0.4904 0.4904 0.2792 0.3384 0.3384 0.2728 0.2728 0.2216 0.2216 0.1587 0.1587 0.0542 0.0542 0.1247 0.1247 0.1050 0.1050 0.0576 0.0576 0.0919 0.0919 0.0883 0.0715 0.0715 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154445.53923139 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1158.06472054 PAW double counting = 50119.77872447 -52759.62314458 entropy T*S EENTRO = -0.00897938 eigenvalues EBANDS = -863.84116424 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1292.19745651 eV energy without entropy = -1292.18847713 energy(sigma->0) = -1292.19446338 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 34) --------------------------------------- eigenvalue-minimisations : 8397 total energy-change (2. order) : 0.4184497E+02 (-0.1360338E+01) number of electron 541.0000039 magnetization 72.8750824 augmentation part 274.2649685 magnetization 58.8654512 Broyden mixing: rms(total) = 0.90680E+02 rms(broyden)= 0.90680E+02 rms(prec ) = 0.95107E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3161 1.3804 1.3804 0.8076 0.8076 0.4862 0.4862 0.2792 0.3293 0.3293 0.2569 0.2569 0.2716 0.2716 0.2233 0.2233 0.1571 0.1571 0.1247 0.1247 0.1068 0.1068 0.0576 0.0576 0.0922 0.0922 0.0887 0.0715 0.0703 0.0703 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154524.77154392 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1158.41757405 PAW double counting = 48922.27854382 -51557.67910314 entropy T*S EENTRO = -0.00042964 eigenvalues EBANDS = -747.56914277 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1250.35248352 eV energy without entropy = -1250.35205388 energy(sigma->0) = -1250.35234031 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 35) --------------------------------------- eigenvalue-minimisations : 7569 total energy-change (2. order) : 0.7027877E+01 (-0.3055191E+00) number of electron 541.0000043 magnetization 72.8466425 augmentation part 273.9821432 magnetization 58.8733195 Broyden mixing: rms(total) = 0.90537E+02 rms(broyden)= 0.90537E+02 rms(prec ) = 0.94954E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3074 1.2867 1.2867 0.7942 0.7942 0.4922 0.4922 0.2792 0.3467 0.3358 0.3358 0.2685 0.2685 0.2093 0.2093 0.2202 0.2202 0.1558 0.1558 0.1240 0.1240 0.1073 0.1073 0.0576 0.0576 0.0924 0.0924 0.0891 0.0719 0.0719 0.0745 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154539.56564427 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1159.05768297 PAW double counting = 48766.61702014 -51402.03883715 entropy T*S EENTRO = 0.02462304 eigenvalues EBANDS = -726.39106948 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1243.32460667 eV energy without entropy = -1243.34922972 energy(sigma->0) = -1243.33281436 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 36) --------------------------------------- eigenvalue-minimisations : 7659 total energy-change (2. order) : 0.8669765E+01 (-0.1350951E+00) number of electron 541.0000040 magnetization 72.7997347 augmentation part 274.2104505 magnetization 58.9199422 Broyden mixing: rms(total) = 0.90335E+02 rms(broyden)= 0.90335E+02 rms(prec ) = 0.94735E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3514 1.7364 1.2252 1.2252 0.7836 0.7836 0.5190 0.5190 0.2792 0.3773 0.3773 0.2917 0.2917 0.2735 0.2735 0.2253 0.2253 0.1549 0.1549 0.1249 0.1249 0.1112 0.1112 0.0576 0.0576 0.0936 0.0936 0.0873 0.0873 0.0831 0.0717 0.0717 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154549.88731401 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1159.51562117 PAW double counting = 48571.28320121 -51205.81319916 entropy T*S EENTRO = 0.01069280 eigenvalues EBANDS = -708.73546222 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1234.65484215 eV energy without entropy = -1234.66553494 energy(sigma->0) = -1234.65840641 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 37) --------------------------------------- eigenvalue-minimisations : 8145 total energy-change (2. order) : 0.3561826E+02 (-0.1574710E+01) number of electron 541.0000007 magnetization 72.4506570 augmentation part 275.7970787 magnetization 58.8858265 Broyden mixing: rms(total) = 0.89259E+02 rms(broyden)= 0.89259E+02 rms(prec ) = 0.93500E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3779 2.2194 1.3925 1.3925 0.8143 0.8143 0.5426 0.5426 0.2792 0.4061 0.4061 0.2999 0.2999 0.2951 0.2951 0.2368 0.2368 0.1584 0.1584 0.1303 0.1303 0.1160 0.1160 0.0576 0.0576 0.0993 0.0993 0.0917 0.0917 0.0898 0.0786 0.0717 0.0717 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154605.80924555 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1161.54491923 PAW double counting = 48065.70766871 -50695.32769622 entropy T*S EENTRO = 0.03067197 eigenvalues EBANDS = -624.15452150 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1199.03658530 eV energy without entropy = -1199.06725726 energy(sigma->0) = -1199.04680928 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 38) --------------------------------------- eigenvalue-minimisations : 8235 total energy-change (2. order) : 0.4016546E+02 (-0.1956307E+01) number of electron 540.9999990 magnetization 72.5062250 augmentation part 275.3150167 magnetization 58.4569249 Broyden mixing: rms(total) = 0.88119E+02 rms(broyden)= 0.88119E+02 rms(prec ) = 0.92338E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3715 2.2950 1.4097 1.4097 0.8181 0.8181 0.5451 0.5451 0.2792 0.4090 0.4090 0.3011 0.3011 0.2960 0.2960 0.2371 0.2371 0.1588 0.1588 0.1307 0.1307 0.1159 0.1159 0.0576 0.0576 0.0992 0.0992 0.0913 0.0913 0.0896 0.0782 0.0717 0.0717 0.0339 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154707.77430794 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1169.99479217 PAW double counting = 48659.41184960 -51289.69841405 entropy T*S EENTRO = 0.00045450 eigenvalues EBANDS = -489.77712186 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1158.87112951 eV energy without entropy = -1158.87158401 energy(sigma->0) = -1158.87128101 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 39) --------------------------------------- eigenvalue-minimisations : 7191 total energy-change (2. order) : 0.8893650E+01 (-0.6049268E+00) number of electron 540.9999995 magnetization 72.3172782 augmentation part 274.5810014 magnetization 57.8938245 Broyden mixing: rms(total) = 0.87935E+02 rms(broyden)= 0.87935E+02 rms(prec ) = 0.92151E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3610 2.2782 1.4057 1.4057 0.8172 0.8172 0.5448 0.5448 0.2792 0.4087 0.4087 0.3007 0.3007 0.2964 0.2964 0.2372 0.2372 0.1588 0.1588 0.0269 0.1308 0.1308 0.1160 0.1160 0.0576 0.0576 0.0993 0.0993 0.0912 0.0912 0.0895 0.0780 0.0717 0.0717 0.0486 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154726.15204730 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1171.12612383 PAW double counting = 48928.79304257 -51559.42196547 entropy T*S EENTRO = -0.00487741 eigenvalues EBANDS = -463.28937426 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1149.97747997 eV energy without entropy = -1149.97260256 energy(sigma->0) = -1149.97585416 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 40) --------------------------------------- eigenvalue-minimisations : 8379 total energy-change (2. order) :-0.2909792E+02 (-0.2072220E+00) number of electron 540.9999995 magnetization 73.0889788 augmentation part 273.8052243 magnetization 58.7718379 Broyden mixing: rms(total) = 0.88727E+02 rms(broyden)= 0.88727E+02 rms(prec ) = 0.93025E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3948 2.8429 1.5079 1.5079 0.8427 0.8427 0.5450 0.5450 0.2792 0.4120 0.4120 0.3690 0.3690 0.2894 0.2894 0.2885 0.2885 0.2366 0.2366 0.1553 0.1553 0.1306 0.1306 0.1190 0.1190 0.0576 0.0576 0.0988 0.0988 0.0927 0.0927 0.0913 0.0913 0.0717 0.0717 0.0797 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154680.53652093 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1168.64908087 PAW double counting = 48298.37731146 -50932.39735338 entropy T*S EENTRO = -0.00518883 eigenvalues EBANDS = -532.13434267 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1179.07539541 eV energy without entropy = -1179.07020657 energy(sigma->0) = -1179.07366580 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 41) --------------------------------------- eigenvalue-minimisations : 10035 total energy-change (2. order) :-0.3978794E+03 (-0.2788310E+03) number of electron 540.9999974 magnetization 73.1168187 augmentation part 270.3370849 magnetization 53.2578357 Broyden mixing: rms(total) = 0.17849E+03 rms(broyden)= 0.17849E+03 rms(prec ) = 0.18083E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3835 2.8332 1.5007 1.5007 0.8430 0.8430 0.5420 0.5420 0.2792 0.4176 0.4176 0.3720 0.3720 0.2893 0.2893 0.2883 0.2883 0.2361 0.2361 0.1560 0.1560 0.0029 0.1302 0.1302 0.1189 0.1189 0.0576 0.0576 0.0989 0.0989 0.0926 0.0926 0.0911 0.0911 0.0717 0.0717 0.0797 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154804.60910315 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1171.98520920 PAW double counting = 50589.26673365 -53216.36496593 entropy T*S EENTRO = 0.02528373 eigenvalues EBANDS = -816.22957483 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1576.95479927 eV energy without entropy = -1576.98008299 energy(sigma->0) = -1576.96322718 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 42) --------------------------------------- eigenvalue-minimisations : 7362 total energy-change (2. order) : 0.4626581E+03 (-0.1675112E+03) number of electron 540.9999959 magnetization 75.2029759 augmentation part 274.2985866 magnetization 63.2181245 Broyden mixing: rms(total) = 0.88192E+02 rms(broyden)= 0.88188E+02 rms(prec ) = 0.92290E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3873 3.0602 1.5261 1.5261 0.8541 0.8541 0.5348 0.5348 0.4749 0.4749 0.2792 0.3839 0.3839 0.2904 0.2904 0.2963 0.2963 0.2386 0.2386 0.1555 0.1555 0.1317 0.1317 0.1188 0.1188 0.0027 0.0576 0.0576 0.1021 0.1021 0.0923 0.0923 0.0889 0.0889 0.0802 0.0717 0.0717 0.0709 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -154827.78999348 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1174.87838479 PAW double counting = 50761.24028310 -53390.22973463 entropy T*S EENTRO = 0.01718515 eigenvalues EBANDS = -331.38440506 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1114.29666206 eV energy without entropy = -1114.31384720 energy(sigma->0) = -1114.30239044 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 43) --------------------------------------- eigenvalue-minimisations : 8694 total energy-change (2. order) : 0.1116302E+03 (-0.1400711E+02) number of electron 540.9999997 magnetization 73.1489212 augmentation part 267.7400947 magnetization 59.2275239 Broyden mixing: rms(total) = 0.84001E+02 rms(broyden)= 0.84001E+02 rms(prec ) = 0.87881E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3692 1.8567 1.3821 1.3821 1.1505 0.8384 0.8384 0.5703 0.5703 0.2792 0.4661 0.4661 0.3750 0.3750 0.2891 0.2891 0.2992 0.2992 0.2411 0.2411 0.1521 0.1521 0.0027 0.1324 0.1324 0.1207 0.1207 0.0576 0.0576 0.1112 0.0944 0.0944 0.0970 0.0949 0.0949 0.0717 0.0717 0.0818 0.0792 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155242.56391480 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1187.17754131 PAW double counting = 64961.24957215 -67599.73722528 entropy T*S EENTRO = -0.02708190 eigenvalues EBANDS = 192.26307428 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1002.66641616 eV energy without entropy = -1002.63933426 energy(sigma->0) = -1002.65738886 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 44) --------------------------------------- eigenvalue-minimisations : 8622 total energy-change (2. order) :-0.1296946E+03 (-0.5753077E+01) number of electron 540.9999989 magnetization 72.1832957 augmentation part 266.4574694 magnetization 60.0355291 Broyden mixing: rms(total) = 0.87442E+02 rms(broyden)= 0.87442E+02 rms(prec ) = 0.91559E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3651 1.5633 1.5633 1.3832 1.3832 0.8441 0.8441 0.5779 0.5779 0.4725 0.4725 0.2792 0.3776 0.3776 0.2893 0.2893 0.3006 0.3006 0.2413 0.2413 0.1531 0.1531 0.0027 0.1329 0.1329 0.1205 0.1205 0.0576 0.0576 0.1047 0.1047 0.0942 0.0942 0.0933 0.0933 0.0717 0.0717 0.0805 0.0789 0.0435 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155176.73904693 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1176.52721258 PAW double counting = 53812.54620731 -56448.77830499 entropy T*S EENTRO = -0.00425098 eigenvalues EBANDS = 5.11554538 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1132.36101956 eV energy without entropy = -1132.35676858 energy(sigma->0) = -1132.35960257 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 45) --------------------------------------- eigenvalue-minimisations : 8361 total energy-change (2. order) :-0.1119239E+02 (-0.1892582E+01) number of electron 540.9999972 magnetization 71.2296259 augmentation part 266.5168595 magnetization 58.0029072 Broyden mixing: rms(total) = 0.89060E+02 rms(broyden)= 0.89060E+02 rms(prec ) = 0.93368E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3625 1.6019 1.6019 1.4012 1.4012 0.8499 0.8499 0.5828 0.5828 0.4764 0.4764 0.2792 0.3794 0.3794 0.2894 0.2894 0.3007 0.3007 0.2407 0.2407 0.1533 0.1533 0.1327 0.1327 0.1204 0.1204 0.0027 0.0783 0.0783 0.0576 0.0576 0.1052 0.1038 0.0942 0.0942 0.0937 0.0937 0.0717 0.0717 0.0793 0.0818 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155164.73099155 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1177.84526922 PAW double counting = 51814.91287753 -54450.11297389 entropy T*S EENTRO = -0.03244165 eigenvalues EBANDS = -20.40676876 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1143.55341104 eV energy without entropy = -1143.52096939 energy(sigma->0) = -1143.54259716 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 46) --------------------------------------- eigenvalue-minimisations : 7947 total energy-change (2. order) :-0.1154545E+02 (-0.2207565E+01) number of electron 541.0000022 magnetization 71.4199069 augmentation part 266.1352486 magnetization 58.7154826 Broyden mixing: rms(total) = 0.88885E+02 rms(broyden)= 0.88885E+02 rms(prec ) = 0.93111E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3545 1.7018 1.3883 1.3883 1.3906 0.8437 0.8437 0.5795 0.5795 0.4738 0.4738 0.2792 0.3786 0.3786 0.2895 0.2895 0.3014 0.3014 0.2410 0.2410 0.1638 0.1539 0.1539 0.0027 0.1331 0.1331 0.1202 0.1202 0.0576 0.0576 0.0931 0.0931 0.1047 0.1047 0.0940 0.0940 0.0933 0.0933 0.0717 0.0717 0.0793 0.0821 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155174.66354316 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1181.07702445 PAW double counting = 50620.64950261 -53249.03020815 entropy T*S EENTRO = -0.04553284 eigenvalues EBANDS = -32.05772214 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1155.09886117 eV energy without entropy = -1155.05332833 energy(sigma->0) = -1155.08368356 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 47) --------------------------------------- eigenvalue-minimisations : 7785 total energy-change (2. order) :-0.5715785E+01 (-0.5220710E+00) number of electron 541.0000027 magnetization 71.5719395 augmentation part 265.8996213 magnetization 58.3640218 Broyden mixing: rms(total) = 0.89250E+02 rms(broyden)= 0.89250E+02 rms(prec ) = 0.93503E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3504 1.9350 1.3822 1.3822 0.8156 0.8380 0.8380 0.6042 0.5781 0.5781 0.4728 0.4728 0.2792 0.3792 0.3792 0.2896 0.2896 0.3021 0.3021 0.2409 0.2409 0.1541 0.1541 0.0027 0.1149 0.1149 0.1325 0.1325 0.1201 0.1201 0.0576 0.0576 0.1045 0.1045 0.0942 0.0942 0.0937 0.0937 0.0817 0.0793 0.0717 0.0717 0.0673 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155176.57771393 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1179.14932129 PAW double counting = 50657.38805463 -53287.99621304 entropy T*S EENTRO = -0.04121236 eigenvalues EBANDS = -31.70850119 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1160.81464654 eV energy without entropy = -1160.77343418 energy(sigma->0) = -1160.80090908 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 48) --------------------------------------- eigenvalue-minimisations : 7821 total energy-change (2. order) : 0.1107920E+02 (-0.1484271E+01) number of electron 541.0000016 magnetization 70.8011502 augmentation part 266.8699079 magnetization 57.8422712 Broyden mixing: rms(total) = 0.88277E+02 rms(broyden)= 0.88277E+02 rms(prec ) = 0.92473E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3786 2.4405 1.4185 1.4185 1.0132 1.0132 0.8387 0.8387 0.5897 0.5897 0.4846 0.4846 0.2792 0.3866 0.3866 0.2901 0.2901 0.3009 0.3009 0.2387 0.2387 0.2073 0.2073 0.1509 0.1509 0.0027 0.1343 0.1343 0.1212 0.1212 0.0576 0.0576 0.1065 0.1065 0.0998 0.0998 0.0932 0.0932 0.0929 0.0929 0.0717 0.0717 0.0793 0.0831 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155061.11529984 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1175.91141053 PAW double counting = 50387.64103836 -53025.08510099 entropy T*S EENTRO = -0.03601264 eigenvalues EBANDS = -126.02310137 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1149.73544788 eV energy without entropy = -1149.69943524 energy(sigma->0) = -1149.72344367 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 49) --------------------------------------- eigenvalue-minimisations : 6966 total energy-change (2. order) : 0.2200215E+02 (-0.7448154E+02) number of electron 541.0000016 magnetization 68.4743304 augmentation part 272.5630753 magnetization 59.5169726 Broyden mixing: rms(total) = 0.65900E+02 rms(broyden)= 0.65899E+02 rms(prec ) = 0.71059E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3848 2.6551 1.4734 1.4734 1.0774 1.0774 0.8443 0.8443 0.5994 0.5994 0.4953 0.4953 0.2792 0.3939 0.3939 0.2910 0.2910 0.3041 0.3041 0.2369 0.2369 0.2279 0.2279 0.1510 0.1510 0.0027 0.1354 0.1354 0.1208 0.1208 0.0576 0.0576 0.1053 0.1053 0.1033 0.1033 0.0939 0.0939 0.0933 0.0933 0.0717 0.0717 0.0822 0.0793 0.0802 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155002.17586756 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1183.73701899 PAW double counting = 51171.14436817 -53808.96601781 entropy T*S EENTRO = -0.00942004 eigenvalues EBANDS = -170.43500061 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1127.73330080 eV energy without entropy = -1127.72388076 energy(sigma->0) = -1127.73016078 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 50) --------------------------------------- eigenvalue-minimisations : 8037 total energy-change (2. order) :-0.1087482E+02 (-0.1524194E+02) number of electron 541.0000044 magnetization 44.9746200 augmentation part 266.2756329 magnetization 32.2437713 Broyden mixing: rms(total) = 0.80761E+02 rms(broyden)= 0.80761E+02 rms(prec ) = 0.85433E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3404 1.4149 1.4149 0.9117 0.9117 0.9089 0.7268 0.7268 0.5370 0.5370 0.5442 0.5442 0.4208 0.3506 0.3506 0.3031 0.3031 0.2975 0.1813 0.1813 0.0023 0.1610 0.1357 0.1357 0.1068 0.1068 0.1119 0.1119 0.1168 0.1130 0.1130 0.0787 0.0787 0.0715 0.0715 0.0907 0.0907 0.0956 0.0956 0.0803 0.0803 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -155023.60083521 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 1196.94578655 PAW double counting = 49943.69453592 -52603.70984418 entropy T*S EENTRO = 0.00523780 eigenvalues EBANDS = -150.91461970 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = -1138.60812075 eV energy without entropy = -1138.61335855 energy(sigma->0) = -1138.60986668 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 51) --------------------------------------- eigenvalue-minimisations : 8127 total energy-change (2. order) : 0.4949620E+04 (-0.2643910E+04) number of electron 540.9999173 magnetization 67.2930153 augmentation part 279.0127924 magnetization 17.2467186 Broyden mixing: rms(total) = 0.13321E+03 rms(broyden)= 0.13321E+03 rms(prec ) = 0.13401E+03 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3347 1.2606 1.2606 0.9941 0.9941 0.9423 0.6728 0.6728 0.7163 0.5380 0.5380 0.4533 0.4533 0.3056 0.3056 0.3272 0.3272 0.3084 0.2499 0.1831 0.1831 0.0022 0.1609 0.1376 0.1376 0.1139 0.1139 0.0734 0.0734 0.1151 0.1151 0.1153 0.1000 0.1000 0.0716 0.0716 0.0930 0.0930 0.0932 0.0932 0.0815 0.0804 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -156382.99109767 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 4197.12677541 PAW double counting = 2029273.17108105 -2029450.21086103 entropy T*S EENTRO = -0.00561469 eigenvalues EBANDS = 674.94983928 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 3811.01174042 eV energy without entropy = 3811.01735511 energy(sigma->0) = 3811.01361198 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 52) --------------------------------------- eigenvalue-minimisations : 8946 total energy-change (2. order) :-0.3291533E+04 (-0.5810447E+03) number of electron 540.9999563 magnetization 78.8599642 augmentation part 273.7984954 magnetization 37.1374470 Broyden mixing: rms(total) = 0.97316E+02 rms(broyden)= 0.97316E+02 rms(prec ) = 0.98215E+02 weight for this iteration 100.00 eigenvalues of (default mixing * dielectric matrix) average eigenvalue GAMMA= 0.3319 1.2504 1.2504 1.0323 1.0323 0.8082 0.8082 0.6724 0.6724 0.5215 0.5215 0.4711 0.4711 0.3058 0.3058 0.3078 0.3078 0.2963 0.2483 0.2483 0.0022 0.1891 0.1891 0.1609 0.1292 0.1292 0.0963 0.0963 0.1204 0.1204 0.0790 0.0790 0.1212 0.1107 0.1107 0.0716 0.0716 0.0918 0.0918 0.0943 0.0943 0.0801 0.0806 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 1416.70056035 Ewald energy TEWEN = 106413.41558026 -Hartree energ DENC = -156357.90131641 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 2807.76008393 PAW double counting = 1494831.42055863 -1497298.54978158 entropy T*S EENTRO = 0.00385099 eigenvalues EBANDS = 1037.77411466 atomic energy EATOM = 47668.85547746 Solvation Ediel_sol = 0.00000000 --------------------------------------------------- free energy TOTEN = 519.47912829 eV energy without entropy = 519.47527731 energy(sigma->0) = 519.47784463 -------------------------------------------------------------------------------------------------------- ----------------------------------------- Iteration 1( 53) --------------------------------------- ----------------------------------------------------------------------------- | | | 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 = 8 1 9 | | | | ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- | | | -----------------------------------------------------------------------------