vasp.6.2.1 16May21 (build Apr 11 2022 11:03:26) complex MD_VERSION_INFO: Compiled 2022-04-11T18:25:55-UTC in devlin.sd.materialsdesign. com:/home/medea2/data/build/vasp6.2.1/16685/x86_64/src/src/build/gpu from svn 1 6685 This VASP executable licensed from Materials Design, Inc. executed on Lin64 date 2024.08.18 15:29:58 running on 2 total cores distrk: each k-point on 2 cores, 1 groups distr: one band on NCORE= 1 cores, 2 groups -------------------------------------------------------------------------------------------------------- INCAR: SYSTEM = No title PREC = Normal ENCUT = 400.000 IBRION = -1 NSW = 0 ISIF = 2 NELMIN = 2 EDIFF = 1.0e-05 EDIFFG = -0.02 VOSKOWN = 1 NBLOCK = 1 NWRITE = 1 NELM = 95 LHFCALC = .TRUE. HFSCREEN = 0.2 PRECFOCK = Normal ALGO = Damped TIME = 0.4 LMAXFOCK = 4 NKREDX = 1 NKREDY = 1 NKREDZ = 1 ISPIN = 1 INIWAV = 1 ISTART = 0 ICHARG = 2 LWAVE = .TRUE. LCHARG = .FALSE. ADDGRID = .FALSE. ISMEAR = 1 SIGMA = 0.2 LREAL = Auto LSCALAPACK = .FALSE. RWIGS = 0.75 0.32 0.77 1.11 0.73 NPAR = 2 POTCAR: PAW_PBE N 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE C 08Apr2002 POTCAR: PAW_PBE Si 05Jan2001 POTCAR: PAW_PBE O 08Apr2002 POTCAR: PAW_PBE N 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE H 15Jun2001 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 POTCAR: PAW_PBE C 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE Si 05Jan2001 local pseudopotential read in partial core-charges read in partial kinetic energy density read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE O 08Apr2002 local pseudopotential read in partial core-charges read in partial kinetic energy density read in kinetic energy density of atom read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 Optimization of the real space projectors (new method) maximal supplied QI-value = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.65 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 10 10.053 79.467 0.76E-04 0.72E-04 0.56E-06 0 10 10.053 66.151 0.76E-04 0.72E-04 0.55E-06 1 10 10.053 8.350 0.25E-03 0.92E-03 0.41E-05 1 10 10.053 5.531 0.27E-03 0.10E-02 0.45E-05 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 9.92, 20.18] = [ 27.55,114.04] Ry Optimized for a Real-space Cutoff 1.26 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 9.919 19.460 0.50E-03 0.23E-03 0.29E-06 0 8 9.919 12.209 0.48E-03 0.23E-03 0.28E-06 1 7 9.919 4.655 0.17E-03 0.75E-03 0.30E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 25.13 optimisation between [QCUT,QGAM] = [ 10.05, 20.36] = [ 28.30,116.06] Ry Optimized for a Real-space Cutoff 1.30 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.053 115.676 0.49E-03 0.72E-03 0.18E-06 0 8 10.053 87.132 0.49E-03 0.71E-03 0.18E-06 1 7 10.053 4.429 0.32E-03 0.31E-03 0.18E-06 1 7 10.053 2.733 0.23E-03 0.19E-03 0.20E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 19.84 optimisation between [QCUT,QGAM] = [ 10.12, 20.44] = [ 28.68,116.96] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 7 10.119 159.560 0.56E-04 0.22E-03 0.45E-07 0 7 10.119 115.863 0.56E-04 0.21E-03 0.45E-07 1 7 10.119 88.339 0.34E-03 0.49E-03 0.11E-06 1 7 10.119 48.592 0.33E-03 0.48E-03 0.11E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 24.76 optimisation between [QCUT,QGAM] = [ 10.15, 20.30] = [ 28.85,115.39] Ry Optimized for a Real-space Cutoff 1.38 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.150 20.381 0.22E-03 0.32E-03 0.29E-06 0 8 10.150 15.268 0.23E-03 0.35E-03 0.30E-06 1 8 10.150 5.964 0.46E-03 0.53E-03 0.21E-06 1 8 10.150 5.382 0.38E-03 0.45E-03 0.19E-06 PAW_PBE N 08Apr2002 : energy of atom 1 EATOM= -264.5486 kinetic energy error for atom= 0.0736 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 2 EATOM= -12.4884 kinetic energy error for atom= 0.0098 (will be added to EATOM!!) PAW_PBE C 08Apr2002 : energy of atom 3 EATOM= -147.1560 kinetic energy error for atom= 0.0288 (will be added to EATOM!!) PAW_PBE Si 05Jan2001 : energy of atom 4 EATOM= -103.0669 kinetic energy error for atom= 0.0012 (will be added to EATOM!!) PAW_PBE O 08Apr2002 : energy of atom 5 EATOM= -432.3788 kinetic energy error for atom= 0.1156 (will be added to EATOM!!) POSCAR: No title positions in direct lattice No initial velocities read in exchange correlation table for LEXCH = 8 RHO(1)= 0.500 N(1) = 2000 RHO(2)= 100.500 N(2) = 4000 -------------------------------------------------------------------------------------------------------- ion position nearest neighbor table 1 0.682 0.462 0.595- 2 1.00 3 1.02 2 0.703 0.466 0.656- 1 1.00 3 0.720 0.459 0.550- 1 1.02 4 0.152 0.517 0.367- 25 1.10 5 0.209 0.609 0.367- 25 1.10 6 0.222 0.518 0.293- 25 1.10 7 0.112 0.548 0.658- 26 1.10 8 0.151 0.636 0.598- 26 1.10 9 0.108 0.550 0.539- 26 1.10 10 0.421 0.429 0.295- 27 1.10 11 0.421 0.339 0.373- 27 1.10 12 0.348 0.408 0.360- 27 1.10 13 0.529 0.269 0.573- 28 1.09 14 0.534 0.352 0.479- 28 1.09 15 0.592 0.420 0.592- 16 0.249 0.689 0.504- 29 1.10 17 0.306 0.666 0.591- 29 1.10 18 0.331 0.655 0.477- 29 1.10 19 0.442 0.618 0.558- 30 1.10 20 0.452 0.555 0.659- 30 1.10 21 0.515 0.551 0.575- 30 1.10 22 0.326 0.445 0.641- 24 1.10 23 0.311 0.377 0.547- 24 1.10 24 0.328 0.444 0.568- 22 1.10 23 1.10 32 1.86 31 1.87 25 0.205 0.536 0.360- 6 1.10 4 1.10 5 1.10 33 1.42 26 0.140 0.564 0.597- 7 1.10 9 1.10 8 1.10 34 1.42 27 0.403 0.407 0.360- 10 1.10 12 1.10 11 1.10 35 1.43 28 0.515 0.333 0.544- 13 1.09 14 1.09 36 1.36 29 0.290 0.644 0.524- 18 1.10 17 1.10 16 1.10 32 1.88 30 0.461 0.554 0.586- 19 1.10 20 1.10 21 1.10 31 1.86 31 0.417 0.458 0.534- 35 1.65 36 1.69 30 1.86 24 1.87 32 0.264 0.524 0.526- 34 1.66 33 1.66 24 1.86 29 1.88 33 0.246 0.490 0.423- 25 1.42 32 1.66 34 0.199 0.512 0.595- 26 1.42 32 1.66 35 0.428 0.468 0.425- 27 1.43 31 1.65 36 0.453 0.362 0.569- 28 1.36 31 1.69 LATTYP: Found a simple tetragonal cell. ALAT = 15.0000000000 C/A-ratio = 1.3333333333 Lattice vectors: A1 = ( 0.0000000000, 15.0000000000, 0.0000000000) A2 = ( 0.0000000000, 0.0000000000, 15.0000000000) A3 = ( 20.0000000000, 0.0000000000, 0.0000000000) Analysis of symmetry for initial positions (statically): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple tetragonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 16 trial point group operations. The static configuration has the point symmetry C_1 . Analysis of symmetry for dynamics (positions and initial velocities): ===================================================================== Subroutine PRICEL returns: Original cell was already a primitive cell. Routine SETGRP: Setting up the symmetry group for a simple tetragonal supercell. Subroutine GETGRP returns: Found 1 space group operations (whereof 1 operations were pure point group operations) out of a pool of 16 trial point group operations. The dynamic configuration has the point symmetry C_1 . Subroutine INISYM returns: Found 1 space group operations (whereof 1 operations are pure point group operations), and found 1 'primitive' translations ---------------------------------------------------------------------------------------- Primitive cell volume of cell : 4500.0000 direct lattice vectors reciprocal lattice vectors 20.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 length of vectors 20.000000000 15.000000000 15.000000000 0.050000000 0.066666667 0.066666667 position of ions in fractional coordinates (direct lattice) 0.682113240 0.461904050 0.595209280 0.702645180 0.465672610 0.655648600 0.719960870 0.458848870 0.549785760 0.152481310 0.516711620 0.367451420 0.209132770 0.609058630 0.367250660 0.222037270 0.517924570 0.293069800 0.111897270 0.547583860 0.657503740 0.151281780 0.636266240 0.597880500 0.107946680 0.550032000 0.538597070 0.420976950 0.428858120 0.294579540 0.420598210 0.338943460 0.372877130 0.347512220 0.407572790 0.359591300 0.529328630 0.269030140 0.573100810 0.533614270 0.351527300 0.478801640 0.592372690 0.420220530 0.591642470 0.248617630 0.688985330 0.504464050 0.306427920 0.665772770 0.590911320 0.331067030 0.654631890 0.477409880 0.442236220 0.617807220 0.558481360 0.451893130 0.555322730 0.658581670 0.514915470 0.550979280 0.575040320 0.326140040 0.444655110 0.641191470 0.311267780 0.376909120 0.547232760 0.327869640 0.444082830 0.567645540 0.205350590 0.536079260 0.360447040 0.140068170 0.564202120 0.596742520 0.402627240 0.407356110 0.360239770 0.514843960 0.332589250 0.544446770 0.289587230 0.644241830 0.524458480 0.460546770 0.554468310 0.586118080 0.417215530 0.457917570 0.534062670 0.263631140 0.524027740 0.526443160 0.245855950 0.490176260 0.423496740 0.199196600 0.511619120 0.595376570 0.427903770 0.468490520 0.425192010 0.453094890 0.361857740 0.568937750 ion indices of the primitive-cell ions primitive index ion index 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31 32 32 33 33 34 34 35 35 36 36 ---------------------------------------------------------------------------------------- KPOINTS: Automatic mesh Automatic generation of k-mesh. Grid dimensions read from file: generate k-points for: 1 1 1 Generating k-lattice: Cartesian coordinates Fractional coordinates (reciprocal lattice) 0.050000000 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 1.000000000 0.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 1.000000000 Length of vectors 0.050000000 0.066666667 0.066666667 Shift w.r.t. Gamma in fractional coordinates (k-lattice) 0.000000000 0.000000000 0.000000000 Subroutine IBZKPT returns following result: =========================================== Found 1 irreducible k-points: Following reciprocal coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 Following cartesian coordinates: Coordinates Weight 0.000000 0.000000 0.000000 1.000000 Subroutine IBZKPT_HF returns following result: ============================================== Found 1 k-points in 1st BZ the following 1 k-points will be used (e.g. in the exchange kernel) Following reciprocal coordinates: # in IRBZ 0.000000 0.000000 0.000000 1.00000000 1 t-inv F -------------------------------------------------------------------------------------------------------- Dimension of arrays: k-points NKPTS = 1 k-points in BZ NKDIM = 1 number of bands NBANDS= 62 number of dos NEDOS = 301 number of ions NIONS = 36 non local maximal LDIM = 4 non local SUM 2l+1 LMDIM = 8 total plane-waves NPLWV = 627200 max r-space proj IRMAX = 2665 max aug-charges IRDMAX= 4859 dimension x,y,z NGX = 98 NGY = 80 NGZ = 80 dimension x,y,z NGXF= 196 NGYF= 160 NGZF= 160 support grid NGXF= 196 NGYF= 160 NGZF= 160 ions per type = 1 22 7 2 4 NGX,Y,Z is equivalent to a cutoff of 8.15, 8.87, 8.87 a.u. NGXF,Y,Z is equivalent to a cutoff of 16.29, 17.73, 17.73 a.u. SYSTEM = No title POSCAR = No title Startparameter for this run: NWRITE = 1 write-flag & timer PREC = normal normal or accurate (medium, high low for compatibility) ISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ISPIN = 1 spin polarized calculation? LNONCOLLINEAR = F non collinear calculations LSORBIT = F spin-orbit coupling INIWAV = 1 electr: 0-lowe 1-rand 2-diag LASPH = F aspherical Exc in radial PAW Electronic Relaxation 1 ENCUT = 400.0 eV 29.40 Ry 5.42 a.u. 32.61 24.46 24.46*2*pi/ulx,y,z ENINI = 400.0 initial cutoff ENAUG = 644.9 eV augmentation charge cutoff NELM = 95; NELMIN= 2; NELMDL= -5 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM LREAL = T real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 1 Vosko Wilk Nusair interpolation ROPT = -0.00050 -0.00050 -0.00050 -0.00050 ROPT = -0.00050 Ionic relaxation EDIFFG = -.2E-01 stopping-criterion for IOM NSW = 0 number of steps for IOM NBLOCK = 1; KBLOCK = 1 inner block; outer block IBRION = -1 ionic relax: 0-MD 1-quasi-New 2-CG NFREE = 0 steps in history (QN), initial steepest desc. (CG) ISIF = 2 stress and relaxation IWAVPR = 10 prediction: 0-non 1-charg 2-wave 3-comb ISYM = 3 0-nonsym 1-usesym 2-fastsym LCORR = T Harris-Foulkes like correction to forces POTIM = 0.5000 time-step for ionic-motion TEIN = 0.0 initial temperature TEBEG = 0.0; TEEND = 0.0 temperature during run SMASS = -3.00 Nose mass-parameter (am) estimated Nose-frequenzy (Omega) = 0.10E-29 period in steps = 0.13E+47 mass= -0.914E-26a.u. SCALEE = 1.0000 scale energy and forces NPACO = 256; APACO = 16.0 distance and # of slots for P.C. PSTRESS= 0.0 pullay stress Mass of Ions in am POMASS = 14.00 1.00 12.01 28.09 16.00 Ionic Valenz ZVAL = 5.00 1.00 4.00 4.00 6.00 Atomic Wigner-Seitz radii RWIGS = 0.75 0.32 0.77 1.11 0.73 virtual crystal weights VCA = 1.00 1.00 1.00 1.00 1.00 NELECT = 87.0000 total number of electrons NUPDOWN= -1.0000 fix difference up-down DOS related values: EMIN = 10.00; EMAX =-10.00 energy-range for DOS EFERMI = 0.00 ISMEAR = 1; SIGMA = 0.20 broadening in eV -4-tet -1-fermi 0-gaus Electronic relaxation 2 (details) IALGO = 53 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.40E-07 absolut break condition DEPER = 0.30 relativ break condition TIME = 0.40 timestep for ELM volume/ion in A,a.u. = 125.00 843.54 Fermi-wavevector in a.u.,A,eV,Ry = 0.439383 0.830314 2.626712 0.193058 Thomas-Fermi vector in A = 1.413435 Write flags LWAVE = T 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 EXXOEP = 0 0=HF, 1=EXX-LHF (local Hartree Fock) 2=EXX OEP LHFCALC = T Hartree Fock is set to LSYMGRAD= F symmetrize gradient (conserves proper symmetry) PRECFOCK=normal Normal, Fast or Accurate (Low or Medium for compatibility) LRHFCALC= F long range Hartree Fock LRSCOR = F long range correlation only (use DFT for short range part) LTHOMAS = F Thomas Fermi screening in HF LMODELHF= F short range full HF, long range fraction AEXX ENCUT4O = -1.0 cutoff for four orbital integrals eV LMAXFOCK= 4 L truncation for augmentation on plane wave grid LMAXFOCKAE= -1 L truncation for all-electron charge restoration on plane wave grid NMAXFOCKAE= 1 number of basis functions for all-electron charge restoration LFOCKAEDFT= F apply the AE augmentation even for DFT NKREDX = 1 reduce k-point grid by NKREDY = 1 reduce k-point grid by NKREDZ = 1 reduce k-point grid by SHIFTRED= F shift reduced grid of Gamma HFKIDENT= F idential grid for each k-point ODDONLY = F use only odd q-grid points EVENONLY= F use only even q-grid points HFALPHA = -1.0000 decay constant for conv. correction MCALPHA = 0.0000 extent of test-charge in conv. correction in multipole expansion AEXX = 0.2500 exact exchange contribution HFSCREEN= 0.2000 screening length (either q_TF or 0.3 A-1) HFSCREENC= 0.2000 screening length for correlation (either q_TF or 0.3 A-1) HFRCUT = 0.0000 spherical cutoff for potential kernel ALDAX = 0.7500 LDA exchange part AGGAX = 0.7500 GGA exchange part ALDAC = 1.0000 LDA correlation AGGAC = 1.0000 GGA correlation NBANDSGWLOW= 1 first orbital included in HF term ENCUTFOCK= -1.0 apply spherical cutoff to Coloumb kernel 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 non-spin polarized calculation Conjugate gradient for all bands (Freysoldt, et al. PRB 79, 241103 (2009)) perform sub-space diagonalisation before 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 18 real space projection scheme for non local part use partial core corrections no Harris-corrections to forces use gradient corrections use of overlap-Matrix (Vanderbilt PP) Methfessel and Paxton Order N= 1 SIGMA = 0.20 -------------------------------------------------------------------------------------------------------- energy-cutoff : 400.00 volume of cell : 4500.00 direct lattice vectors reciprocal lattice vectors 20.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 length of vectors 20.000000000 15.000000000 15.000000000 0.050000000 0.066666667 0.066666667 k-points in units of 2pi/SCALE and weight: Automatic mesh 0.00000000 0.00000000 0.00000000 1.000 k-points in reciprocal lattice and weights: Automatic mesh 0.00000000 0.00000000 0.00000000 1.000 position of ions in fractional coordinates (direct lattice) 0.68211324 0.46190405 0.59520928 0.70264518 0.46567261 0.65564860 0.71996087 0.45884887 0.54978576 0.15248131 0.51671162 0.36745142 0.20913277 0.60905863 0.36725066 0.22203727 0.51792457 0.29306980 0.11189727 0.54758386 0.65750374 0.15128178 0.63626624 0.59788050 0.10794668 0.55003200 0.53859707 0.42097695 0.42885812 0.29457954 0.42059821 0.33894346 0.37287713 0.34751222 0.40757279 0.35959130 0.52932863 0.26903014 0.57310081 0.53361427 0.35152730 0.47880164 0.59237269 0.42022053 0.59164247 0.24861763 0.68898533 0.50446405 0.30642792 0.66577277 0.59091132 0.33106703 0.65463189 0.47740988 0.44223622 0.61780722 0.55848136 0.45189313 0.55532273 0.65858167 0.51491547 0.55097928 0.57504032 0.32614004 0.44465511 0.64119147 0.31126778 0.37690912 0.54723276 0.32786964 0.44408283 0.56764554 0.20535059 0.53607926 0.36044704 0.14006817 0.56420212 0.59674252 0.40262724 0.40735611 0.36023977 0.51484396 0.33258925 0.54444677 0.28958723 0.64424183 0.52445848 0.46054677 0.55446831 0.58611808 0.41721553 0.45791757 0.53406267 0.26363114 0.52402774 0.52644316 0.24585595 0.49017626 0.42349674 0.19919660 0.51161912 0.59537657 0.42790377 0.46849052 0.42519201 0.45309489 0.36185774 0.56893775 position of ions in cartesian coordinates (Angst): 13.64226480 6.92856075 8.92813920 14.05290360 6.98508915 9.83472900 14.39921740 6.88273305 8.24678640 3.04962620 7.75067430 5.51177130 4.18265540 9.13587945 5.50875990 4.44074540 7.76886855 4.39604700 2.23794540 8.21375790 9.86255610 3.02563560 9.54399360 8.96820750 2.15893360 8.25048000 8.07895605 8.41953900 6.43287180 4.41869310 8.41196420 5.08415190 5.59315695 6.95024440 6.11359185 5.39386950 10.58657260 4.03545210 8.59651215 10.67228540 5.27290950 7.18202460 11.84745380 6.30330795 8.87463705 4.97235260 10.33477995 7.56696075 6.12855840 9.98659155 8.86366980 6.62134060 9.81947835 7.16114820 8.84472440 9.26710830 8.37722040 9.03786260 8.32984095 9.87872505 10.29830940 8.26468920 8.62560480 6.52280080 6.66982665 9.61787205 6.22535560 5.65363680 8.20849140 6.55739280 6.66124245 8.51468310 4.10701180 8.04118890 5.40670560 2.80136340 8.46303180 8.95113780 8.05254480 6.11034165 5.40359655 10.29687920 4.98883875 8.16670155 5.79174460 9.66362745 7.86687720 9.21093540 8.31702465 8.79177120 8.34431060 6.86876355 8.01094005 5.27262280 7.86041610 7.89664740 4.91711900 7.35264390 6.35245110 3.98393200 7.67428680 8.93064855 8.55807540 7.02735780 6.37788015 9.06189780 5.42786610 8.53406625 -------------------------------------------------------------------------------------------------------- k-point 1 : 0.0000 0.0000 0.0000 plane waves: 81909 maximum and minimum number of plane-waves per node : 81909 81909 maximum number of plane-waves: 81909 maximum index in each direction: IXMAX= 32 IYMAX= 24 IZMAX= 24 IXMIN= -32 IYMIN= -24 IZMIN= -24 The following grids will avoid any aliasing or wrap around errors in the Hartre e energy - symmetry arguments have not been applied - exchange correlation energies might require even more grid points - we recommend to set PREC=Normal or Accurate and rely on VASP defaults WARNING: aliasing errors must be expected set NGX to 140 to avoid them WARNING: aliasing errors must be expected set NGY to 98 to avoid them WARNING: aliasing errors must be expected set NGZ to 98 to avoid them serial 3D FFT for wavefunctions parallel 3D FFT for charge: minimum data exchange during FFTs selected (reduces bandwidth) Radii for the augmentation spheres in the non-local exchange for species 1 augmentation radius 0.886 (default was 0.709) energy cutoff for augmentation 1600.0 for species 2 augmentation radius 0.650 (default was 0.520) energy cutoff for augmentation 1600.0 for species 3 augmentation radius 0.874 (default was 0.699) energy cutoff for augmentation 1600.0 for species 4 augmentation radius 1.237 (default was 0.989) energy cutoff for augmentation 1600.0 for species 5 augmentation radius 0.902 (default was 0.722) energy cutoff for augmentation 1600.0 SETUP_FOCK is finished total amount of memory used by VASP MPI-rank0 473576. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 5602. kBytes fftplans : 87706. kBytes grid : 227924. kBytes one-center: 110. kBytes HF : 23. kBytes wavefun : 122211. kBytes Broyden mixing: mesh for mixing (old mesh) NGX = 65 NGY = 49 NGZ = 49 (NGX =196 NGY =160 NGZ =160) gives a total of 156065 points initial charge density was supplied: charge density of overlapping atoms calculated number of electron 87.0000000 magnetization keeping initial charge density in first step -------------------------------------------------------------------------------------------------------- Maximum index for non-local projection operator 2525 Maximum index for augmentation-charges 2261 (set IRDMAX) -------------------------------------------------------------------------------------------------------- First call to EWALD: gamma= 0.107 Maximum number of real-space cells 2x 3x 3 Maximum number of reciprocal cells 3x 3x 3 --------------------------------------- Iteration 1( 1) --------------------------------------- eigenvalue-minimisations : 194 total energy-change (2. order) : 0.7287470E+03 (-0.1559388E+04) number of electron 87.0000000 magnetization augmentation part 87.0000000 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8312.24703371 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 254.03747802 PAW double counting = 2573.98116666 -3109.12938481 entropy T*S EENTRO = -0.00109136 eigenvalues EBANDS = 106.67394675 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = 728.74701473 eV energy without entropy = 728.74810609 energy(sigma->0) = 728.74737851 exchange ACFDT corr. = 0.00000000 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 2) --------------------------------------- eigenvalue-minimisations : 248 total energy-change (2. order) :-0.7456982E+03 (-0.6687710E+03) number of electron 87.0000000 magnetization augmentation part 87.0000000 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8312.24703371 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 254.03747802 PAW double counting = 2573.98116666 -3109.12938481 entropy T*S EENTRO = 0.01497359 eigenvalues EBANDS = -639.04027236 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -16.95113943 eV energy without entropy = -16.96611302 energy(sigma->0) = -16.95613063 exchange ACFDT corr. = 0.00000000 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 3) --------------------------------------- eigenvalue-minimisations : 248 total energy-change (2. order) :-0.5497211E+03 (-0.5246466E+03) number of electron 87.0000000 magnetization augmentation part 87.0000000 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8312.24703371 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 254.03747802 PAW double counting = 2573.98116666 -3109.12938481 entropy T*S EENTRO = -0.02672169 eigenvalues EBANDS = -1188.71969441 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -566.67225676 eV energy without entropy = -566.64553508 energy(sigma->0) = -566.66334954 exchange ACFDT corr. = 0.00000000 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 4) --------------------------------------- eigenvalue-minimisations : 248 total energy-change (2. order) :-0.2104695E+03 (-0.2085192E+03) number of electron 87.0000000 magnetization augmentation part 87.0000000 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8312.24703371 -exchange EXHF = 0.00000000 -V(xc)+E(xc) XCENC = 254.03747802 PAW double counting = 2573.98116666 -3109.12938481 entropy T*S EENTRO = -0.05641896 eigenvalues EBANDS = -1399.15952164 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -777.14178126 eV energy without entropy = -777.08536230 energy(sigma->0) = -777.12297494 exchange ACFDT corr. = 0.00000000 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 5) --------------------------------------- eigenvalue-minimisations : 240 total energy-change (2. order) : 0.5370967E+03 (-0.5056307E+02) number of electron 86.9999971 magnetization augmentation part 4.4808974 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8312.24703371 -exchange EXHF = 212.48597996 -V(xc)+E(xc) XCENC = 195.21676828 PAW double counting = 2623.18677848 -2602.98652763 entropy T*S EENTRO = -0.05641896 eigenvalues EBANDS = -1571.07656248 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -240.04508288 eV energy without entropy = -239.98866392 energy(sigma->0) = -240.02627656 exchange ACFDT corr. = 0.00000000 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 6) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) : 0.5624121E+02 (-0.7879031E+02) number of electron 86.9999978 magnetization augmentation part 3.5785590 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -9016.37816568 -exchange EXHF = 221.65443361 -V(xc)+E(xc) XCENC = 224.29739102 PAW double counting = 10552.99748346 -10537.68516851 entropy T*S EENTRO = -0.05641896 eigenvalues EBANDS = -844.06536552 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -183.80387739 eV energy without entropy = -183.74745843 energy(sigma->0) = -183.78507107 exchange ACFDT corr. = -0.35907589 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 7) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.2804106E+02 (-0.1201890E+02) number of electron 86.9999979 magnetization augmentation part 3.3920349 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8500.09096811 -exchange EXHF = 199.89568572 -V(xc)+E(xc) XCENC = 204.46056544 PAW double counting = 10543.18602670 -10525.17034634 entropy T*S EENTRO = -0.05641896 eigenvalues EBANDS = -1349.50141605 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -211.84493842 eV energy without entropy = -211.78851946 energy(sigma->0) = -211.82613210 exchange ACFDT corr. = -0.37596365 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 8) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.5501106E+01 (-0.8285690E+01) number of electron 86.9999980 magnetization augmentation part 3.4996308 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8567.06368223 -exchange EXHF = 201.28066595 -V(xc)+E(xc) XCENC = 205.43447535 PAW double counting = 10181.91471270 -10163.14098887 entropy T*S EENTRO = -0.05641896 eigenvalues EBANDS = -1291.14674200 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -217.34604488 eV energy without entropy = -217.28962592 energy(sigma->0) = -217.32723856 exchange ACFDT corr. = -0.37788790 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 9) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.3228288E+01 (-0.2168219E+01) number of electron 86.9999981 magnetization augmentation part 3.5526894 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8720.00274120 -exchange EXHF = 207.07433687 -V(xc)+E(xc) XCENC = 210.59819134 PAW double counting = 9739.06067487 -9720.44704447 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1152.23326469 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -220.57433305 eV energy without entropy = -220.51791411 energy(sigma->0) = -220.55552674 exchange ACFDT corr. = -0.37223790 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 10) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1132243E+01 (-0.1219490E+01) number of electron 86.9999982 magnetization augmentation part 3.5325624 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8738.75468089 -exchange EXHF = 208.01534697 -V(xc)+E(xc) XCENC = 211.45051683 PAW double counting = 9397.78476242 -9379.35663432 entropy T*S EENTRO = -0.05641807 eigenvalues EBANDS = -1136.22140119 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -221.70657594 eV energy without entropy = -221.65015787 energy(sigma->0) = -221.68776992 exchange ACFDT corr. = -0.36208168 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 11) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.5047292E+00 (-0.4526891E+00) number of electron 86.9999982 magnetization augmentation part 3.5316862 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.44946801 -exchange EXHF = 206.61779831 -V(xc)+E(xc) XCENC = 210.19956907 PAW double counting = 9173.75388836 -9155.28416493 entropy T*S EENTRO = -0.05641883 eigenvalues EBANDS = -1176.42444309 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.21130519 eV energy without entropy = -222.15488636 energy(sigma->0) = -222.19249892 exchange ACFDT corr. = -0.36638285 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 12) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.2081155E+00 (-0.1650507E+00) number of electron 86.9999982 magnetization augmentation part 3.5397157 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8684.77016699 -exchange EXHF = 206.38876972 -V(xc)+E(xc) XCENC = 209.99569057 PAW double counting = 9074.73978868 -9056.23824327 entropy T*S EENTRO = -0.05641889 eigenvalues EBANDS = -1187.91077373 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.41942066 eV energy without entropy = -222.36300177 energy(sigma->0) = -222.40061436 exchange ACFDT corr. = -0.36549455 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 13) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.8676873E-01 (-0.7381639E-01) number of electron 86.9999982 magnetization augmentation part 3.5384080 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.76681072 -exchange EXHF = 206.68194491 -V(xc)+E(xc) XCENC = 210.27650816 PAW double counting = 9046.96946236 -9028.48447809 entropy T*S EENTRO = -0.05641889 eigenvalues EBANDS = -1176.55833030 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.50618938 eV energy without entropy = -222.44977050 energy(sigma->0) = -222.48738309 exchange ACFDT corr. = -0.36666290 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 14) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.3502174E-01 (-0.3257504E-01) number of electron 86.9999982 magnetization augmentation part 3.5340179 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8701.73263031 -exchange EXHF = 206.71871759 -V(xc)+E(xc) XCENC = 210.32027476 PAW double counting = 9045.18590418 -9026.71025444 entropy T*S EENTRO = -0.05641892 eigenvalues EBANDS = -1171.69873720 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.54121112 eV energy without entropy = -222.48479220 energy(sigma->0) = -222.52240481 exchange ACFDT corr. = -0.36859277 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 15) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1887989E-01 (-0.1127797E-01) number of electron 86.9999982 magnetization augmentation part 3.5324281 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8697.71665115 -exchange EXHF = 206.62830391 -V(xc)+E(xc) XCENC = 210.23660057 PAW double counting = 9053.18152816 -9034.69965595 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1175.56573080 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.56009101 eV energy without entropy = -222.50367207 energy(sigma->0) = -222.54128470 exchange ACFDT corr. = -0.36740682 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 16) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.6023698E-02 (-0.5729011E-02) number of electron 86.9999982 magnetization augmentation part 3.5330997 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8695.07210171 -exchange EXHF = 206.58781750 -V(xc)+E(xc) XCENC = 210.19672063 PAW double counting = 9061.85026307 -9043.36041349 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1178.14391495 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.56611471 eV energy without entropy = -222.50969577 energy(sigma->0) = -222.54730840 exchange ACFDT corr. = -0.36692084 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 17) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.3443257E-02 (-0.2036750E-02) number of electron 86.9999982 magnetization augmentation part 3.5337413 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8695.82314546 -exchange EXHF = 206.60072651 -V(xc)+E(xc) XCENC = 210.20909709 PAW double counting = 9069.54324690 -9051.05266905 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1177.42232821 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.56955797 eV energy without entropy = -222.51313903 energy(sigma->0) = -222.55075165 exchange ACFDT corr. = -0.36781470 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 18) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1340470E-02 (-0.1115032E-02) number of electron 86.9999982 magnetization augmentation part 3.5335679 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.55453792 -exchange EXHF = 206.61474042 -V(xc)+E(xc) XCENC = 210.22312270 PAW double counting = 9076.14334761 -9057.65429551 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.71878997 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57089844 eV energy without entropy = -222.51447950 energy(sigma->0) = -222.55209212 exchange ACFDT corr. = -0.36832983 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 19) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.7240022E-03 (-0.4685334E-03) number of electron 86.9999982 magnetization augmentation part 3.5331973 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.38427272 -exchange EXHF = 206.60960389 -V(xc)+E(xc) XCENC = 210.21923134 PAW double counting = 9081.04429751 -9062.55521094 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.88078577 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57162244 eV energy without entropy = -222.51520350 energy(sigma->0) = -222.55281612 exchange ACFDT corr. = -0.36842115 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 20) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.3454722E-03 (-0.2487754E-03) number of electron 86.9999982 magnetization augmentation part 3.5331469 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.20895757 -exchange EXHF = 206.60166599 -V(xc)+E(xc) XCENC = 210.21260616 PAW double counting = 9083.77260415 -9065.28217376 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1177.04322712 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57196791 eV energy without entropy = -222.51554897 energy(sigma->0) = -222.55316160 exchange ACFDT corr. = -0.36869688 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 21) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1833790E-03 (-0.1245547E-03) number of electron 86.9999982 magnetization augmentation part 3.5332308 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.23409578 -exchange EXHF = 206.60070103 -V(xc)+E(xc) XCENC = 210.21207969 PAW double counting = 9085.13684689 -9066.64565351 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1177.01754386 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57215129 eV energy without entropy = -222.51573235 energy(sigma->0) = -222.55334498 exchange ACFDT corr. = -0.36898513 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 22) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1003595E-03 (-0.6796661E-04) number of electron 86.9999982 magnetization augmentation part 3.5332117 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.25986101 -exchange EXHF = 206.60312469 -V(xc)+E(xc) XCENC = 210.21439766 PAW double counting = 9085.88123203 -9067.38985818 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99680109 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57225165 eV energy without entropy = -222.51583271 energy(sigma->0) = -222.55344533 exchange ACFDT corr. = -0.36911161 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 23) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.5646421E-04 (-0.4132215E-04) number of electron 86.9999982 magnetization augmentation part 3.5331105 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.25994141 -exchange EXHF = 206.60360671 -V(xc)+E(xc) XCENC = 210.21496434 PAW double counting = 9086.24427860 -9067.75276232 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99796827 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57230811 eV energy without entropy = -222.51588917 energy(sigma->0) = -222.55350180 exchange ACFDT corr. = -0.36921004 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 24) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.3508174E-04 (-0.2434466E-04) number of electron 86.9999982 magnetization augmentation part 3.5330384 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.26000750 -exchange EXHF = 206.60252163 -V(xc)+E(xc) XCENC = 210.21414223 PAW double counting = 9086.35415193 -9067.86236704 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99629870 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57234319 eV energy without entropy = -222.51592425 energy(sigma->0) = -222.55353688 exchange ACFDT corr. = -0.36932669 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 25) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.2133523E-04 (-0.1516835E-04) number of electron 86.9999982 magnetization augmentation part 3.5330061 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.25904520 -exchange EXHF = 206.60206059 -V(xc)+E(xc) XCENC = 210.21381473 PAW double counting = 9086.34136745 -9067.84934058 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99673576 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57236453 eV energy without entropy = -222.51594559 energy(sigma->0) = -222.55355822 exchange ACFDT corr. = -0.36940797 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 26) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.1341030E-04 (-0.9472568E-05) number of electron 86.9999982 magnetization augmentation part 3.5329801 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.25833018 -exchange EXHF = 206.60220193 -V(xc)+E(xc) XCENC = 210.21398014 PAW double counting = 9086.28973956 -9067.79756019 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99792344 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57237794 eV energy without entropy = -222.51595900 energy(sigma->0) = -222.55357163 exchange ACFDT corr. = -0.36945885 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- --------------------------------------- Iteration 1( 27) --------------------------------------- eigenvalue-minimisations : 0 total energy-change (2. order) :-0.8524724E-05 (-0.6205726E-05) number of electron 86.9999982 magnetization augmentation part 3.5329505 magnetization Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 6.34847707 Ewald energy TEWEN = 5744.67698435 -Hartree energ DENC = -8696.26272132 -exchange EXHF = 206.60225821 -V(xc)+E(xc) XCENC = 210.21412352 PAW double counting = 9086.23751284 -9067.74522437 entropy T*S EENTRO = -0.05641894 eigenvalues EBANDS = -1176.99391005 atomic energy EATOM = 3464.40647177 --------------------------------------------------- free energy TOTEN = -222.57238646 eV energy without entropy = -222.51596752 energy(sigma->0) = -222.55358015 exchange ACFDT corr. = -0.36950510 see jH, gK, PRB 81, 115126 -------------------------------------------------------------------------------------------------------- average scaling for gradient 1.0194 average (electrostatic) potential at core the test charge radii are 0.7089 0.5201 0.6991 0.9892 0.7215 (the norm of the test charge is 1.0000) 1 -66.6559 2 -38.0336 3 -37.9452 4 -37.5531 5 -37.6647 6 -37.5647 7 -37.4684 8 -37.5455 9 -37.4975 10 -37.4631 11 -37.5254 12 -37.3455 13 -38.3034 14 -38.2207 15 -36.8919 16 -37.6530 17 -37.5971 18 -37.5722 19 -37.4954 20 -37.4756 21 -37.3281 22 -37.7609 23 -37.7471 24 -54.0268 25 -55.7426 26 -55.6746 27 -55.6039 28 -56.7905 29 -54.1408 30 -54.0200 31 -90.0450 32 -89.7835 33 -73.8818 34 -73.7891 35 -73.9612 36 -75.0915 E-fermi : -4.2851 XC(G=0): -0.8415 alpha+bet : -0.4834 k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -28.8912 2.00000 2 -27.7773 2.00000 3 -27.5871 2.00000 4 -27.3264 2.00000 5 -21.0790 2.00000 6 -19.7206 2.00000 7 -19.1276 2.00000 8 -18.7761 2.00000 9 -18.4689 2.00000 10 -18.4032 2.00000 11 -18.3478 2.00000 12 -18.0645 2.00000 13 -14.5556 2.00000 14 -14.0326 2.00000 15 -13.2312 2.00000 16 -12.9922 2.00000 17 -12.8876 2.00000 18 -12.2416 2.00000 19 -12.1845 2.00000 20 -12.0425 2.00000 21 -11.9097 2.00000 22 -11.8035 2.00000 23 -11.6227 2.00000 24 -11.3786 2.00000 25 -11.3044 2.00000 26 -11.2699 2.00000 27 -10.7135 2.00000 28 -10.6417 2.00000 29 -10.4342 2.00000 30 -10.3713 2.00000 31 -10.2502 2.00000 32 -9.8980 2.00000 33 -9.2424 2.00000 34 -8.9578 2.00000 35 -8.6510 2.00000 36 -8.3456 2.00000 37 -8.1735 2.00000 38 -8.1167 2.00000 39 -7.9711 2.00000 40 -7.6393 2.00000 41 -7.4569 2.00000 42 -7.3152 2.00000 43 -5.1491 2.00000 44 -4.2853 1.00000 45 -1.0409 -0.00000 46 -0.3251 -0.00000 47 0.0843 -0.00000 48 0.1953 -0.00000 49 0.2991 -0.00000 50 0.3940 -0.00000 51 0.4741 -0.00000 52 0.6186 -0.00000 53 0.6740 -0.00000 54 0.7783 -0.00000 55 0.8670 -0.00000 56 0.9485 -0.00000 57 0.9872 -0.00000 58 1.0228 -0.00000 59 1.0807 0.00000 60 1.1330 0.00000 61 1.2824 0.00000 62 1.3952 0.00000 -------------------------------------------------------------------------------------------------------- soft charge-density along one line, spin component 1 0 1 2 3 4 5 6 7 8 9 total charge-density along one line pseudopotential strength for first ion, spin component: 1 13.066 15.677 0.001 -0.002 -0.007 0.006 -0.006 -0.021 15.677 18.810 0.002 -0.003 -0.009 0.008 -0.007 -0.025 0.001 0.002 -6.503 -0.010 -0.031 -8.829 -0.015 -0.049 -0.002 -0.003 -0.010 -6.559 0.008 -0.015 -8.918 0.012 -0.007 -0.009 -0.031 0.008 -6.535 -0.049 0.012 -8.880 0.006 0.008 -8.829 -0.015 -0.049 -11.301 -0.023 -0.075 -0.006 -0.007 -0.015 -8.918 0.012 -0.023 -11.440 0.019 -0.021 -0.025 -0.049 0.012 -8.880 -0.075 0.019 -11.380 total augmentation occupancy for first ion, spin component: 1 7.029 -3.137 -0.070 0.294 1.166 0.002 -0.079 -0.331 -3.137 1.524 0.159 -0.240 -0.939 -0.009 0.049 0.202 -0.070 0.159 1.616 0.013 -0.229 -0.127 -0.023 0.004 0.294 -0.240 0.013 2.463 -0.046 -0.023 -0.493 0.035 1.166 -0.939 -0.229 -0.046 2.100 0.005 0.035 -0.343 0.002 -0.009 -0.127 -0.023 0.005 0.010 0.006 0.002 -0.079 0.049 -0.023 -0.493 0.035 0.006 0.099 -0.010 -0.331 0.202 0.004 0.035 -0.343 0.002 -0.010 0.062 ------------------------ aborting loop because EDIFF is reached ---------------------------------------- FORCE on cell =-STRESS in cart. coord. units (eV): Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 6.34848 6.34848 6.34848 Ewald 2946.00002 1035.95968 1762.71449 -1097.33958 186.52659 230.89377 Hartree 3789.84242 2219.88992 2686.53307 -879.33496 160.02882 222.38189 E(xc) -257.45601 -258.12475 -257.72330 -0.32231 -0.00534 0.03224 Local -7631.18169 -4183.19796 -5324.70895 1956.18412 -344.73746 -455.03382 n-local -103.78807 -105.90445 -106.59717 0.29088 2.53726 1.79409 augment 15.85278 17.52583 16.31431 0.99571 -0.55696 -0.47668 Kinetic 1321.13489 1359.84608 1312.14728 16.89594 -4.89881 0.07771 Fock -23.10266 -23.55735 -23.42685 -0.13973 0.37382 0.38595 ------------------------------------------------------------------------------------- Total 63.6501569 68.7854828 71.6013512 -2.7699374 -0.7320926 0.0551534 in kB 22.6619641 24.4903425 25.4929026 -0.9862069 -0.2606538 0.0196368 external PRESSURE = 24.2150697 kB Pullay stress = 0.0000000 kB VOLUME and BASIS-vectors are now : ----------------------------------------------------------------------------- energy-cutoff : 400.00 volume of cell : 4500.00 direct lattice vectors reciprocal lattice vectors 20.000000000 0.000000000 0.000000000 0.050000000 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 0.000000000 0.000000000 0.000000000 15.000000000 0.000000000 0.000000000 0.066666667 length of vectors 20.000000000 15.000000000 15.000000000 0.050000000 0.066666667 0.066666667 FORCES acting on ions: Electron-Ion Ewald-Force Non-Local-Force ----------------------------------------------------------------------------------------------- -.648E+02 -.307E+01 -.857E+01 0.419E+02 0.142E+02 0.222E+01 0.191E+02 -.910E+01 0.509E+01 -.427E+02 -.601E+01 -.693E+02 0.459E+02 0.609E+01 0.774E+02 -.296E+01 -.629E+00 -.678E+01 -.646E+02 0.119E+01 0.457E+02 0.696E+02 -.211E+01 -.508E+02 -.506E+01 0.641E-01 0.466E+01 0.857E+02 0.121E+02 0.257E+02 -.910E+02 -.136E+02 -.252E+02 0.524E+01 0.145E+01 -.479E+00 0.215E+02 -.760E+02 0.317E+02 -.211E+02 0.814E+02 -.313E+02 -.351E+00 -.540E+01 -.420E+00 0.476E+01 0.777E+01 0.872E+02 -.307E+01 -.915E+01 -.923E+02 -.167E+01 0.136E+01 0.509E+01 0.619E+02 0.541E+01 -.641E+02 -.647E+02 -.667E+01 0.686E+02 0.283E+01 0.125E+01 -.457E+01 0.260E+02 -.773E+02 -.175E+02 -.250E+02 0.826E+02 0.176E+02 -.104E+01 -.529E+01 -.124E+00 0.764E+02 0.204E+01 0.358E+02 -.796E+02 -.310E+01 -.400E+02 0.319E+01 0.105E+01 0.424E+01 -.269E+02 -.543E+01 0.884E+02 0.287E+02 0.704E+01 -.933E+02 -.183E+01 -.162E+01 0.495E+01 -.211E+02 0.806E+02 0.284E+02 0.229E+02 -.856E+02 -.275E+02 -.174E+01 0.505E+01 -.860E+00 0.515E+02 0.269E+02 0.433E+02 -.571E+02 -.268E+02 -.434E+02 0.554E+01 -.484E-01 0.738E-01 -.393E+02 0.751E+02 -.264E+02 0.411E+02 -.804E+02 0.290E+02 -.166E+01 0.515E+01 -.234E+01 -.572E+02 0.123E+02 0.564E+02 0.593E+02 -.107E+02 -.617E+02 -.214E+01 -.155E+01 0.514E+01 -.227E+02 0.240E+01 -.136E+02 0.239E+02 -.294E+01 0.133E+02 0.247E-01 -.242E+00 -.108E+00 0.412E+02 -.739E+02 0.145E+02 -.451E+02 0.771E+02 -.159E+02 0.384E+01 -.324E+01 0.141E+01 -.169E+02 -.536E+02 -.612E+02 0.185E+02 0.552E+02 0.659E+02 -.161E+01 -.159E+01 -.473E+01 -.425E+02 -.509E+02 0.389E+02 0.465E+02 0.517E+02 -.423E+02 -.396E+01 -.787E+00 0.337E+01 -.851E+01 -.772E+02 0.659E+01 0.682E+01 0.817E+02 -.858E+01 0.169E+01 -.459E+01 0.197E+01 -.106E+02 -.226E+02 -.795E+02 0.981E+01 0.227E+02 0.847E+02 0.759E+00 -.111E+00 -.518E+01 -.758E+02 -.258E+02 -.531E+01 0.810E+02 0.256E+02 0.453E+01 -.526E+01 0.206E+00 0.769E+00 0.747E+01 0.153E+02 -.874E+02 -.770E+01 -.153E+02 0.926E+02 0.224E+00 0.231E-01 -.516E+01 0.304E+02 0.823E+02 -.340E+01 -.320E+02 -.870E+02 0.202E+01 0.157E+01 0.472E+01 0.138E+01 0.310E+02 0.750E+02 -.103E+03 -.303E+02 -.726E+02 0.101E+03 -.556E+00 -.237E+01 0.195E+01 0.153E+03 -.665E+02 0.198E+03 -.157E+03 0.691E+02 -.202E+03 0.352E+01 -.265E+01 0.440E+01 0.224E+03 -.870E+02 -.712E+02 -.230E+03 0.902E+02 0.717E+02 0.530E+01 -.315E+01 -.459E+00 -.865E+00 0.121E+03 0.217E+03 -.192E+01 -.125E+03 -.222E+03 0.269E+01 0.378E+01 0.470E+01 -.219E+03 0.141E+03 -.304E+01 0.226E+03 -.143E+03 0.940E+00 -.723E+01 0.384E+01 0.253E+01 -.784E+01 -.178E+03 -.170E+02 0.722E+01 0.175E+03 0.171E+02 0.611E+00 0.272E+01 0.123E-02 -.100E+03 -.122E+03 -.855E+02 0.100E+03 0.121E+03 0.845E+02 0.399E+00 0.135E+01 0.107E+01 -.242E+02 -.338E+02 -.757E+02 0.228E+02 0.376E+02 0.760E+02 0.123E+01 -.336E+01 -.651E+00 0.711E+00 -.346E+02 -.469E+02 0.844E+00 0.364E+02 0.474E+02 -.186E+01 -.189E+01 -.525E+00 0.694E+02 0.136E+03 0.145E+03 -.535E+02 -.166E+03 -.149E+03 -.140E+02 0.265E+02 0.346E+01 0.125E+03 0.110E+03 -.204E+03 -.117E+03 -.137E+03 0.225E+03 -.693E+01 0.236E+02 -.182E+02 -.167E+03 -.130E+03 0.169E+03 0.186E+03 0.159E+03 -.176E+03 -.162E+02 -.259E+02 0.623E+01 -.225E+01 0.194E+03 -.178E+03 -.237E+02 -.207E+03 0.200E+03 0.238E+02 0.107E+02 -.199E+02 ----------------------------------------------------------------------------------------------- -.506E+01 -.231E+02 0.109E+02 -.853E-13 -.284E-13 -.853E-13 0.540E+01 0.192E+02 -.797E+01 POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 13.64226 6.92856 8.92814 -2.230800 1.118621 -0.812955 14.05290 6.98509 9.83473 0.135746 -0.569137 1.044184 14.39922 6.88273 8.24679 -0.348600 -0.847596 -0.227099 3.04963 7.75067 5.51177 0.216848 0.041943 0.004662 4.18266 9.13588 5.50876 0.009469 -0.201675 0.001529 4.44075 7.76887 4.39605 -0.043601 0.034979 0.195141 2.23795 8.21376 9.86256 0.152949 0.051225 -0.221758 3.02564 9.54399 8.96821 -0.028484 -0.255923 -0.002724 2.15893 8.25048 8.07896 0.168053 0.035546 0.204647 8.41954 6.43287 4.41869 -0.076080 -0.075001 0.246852 8.41196 5.08415 5.59316 -0.079983 0.245254 -0.043791 6.95024 6.11359 5.39387 0.271102 0.010683 0.024527 10.58657 4.03545 8.59651 0.047537 0.153529 0.117768 10.67229 5.27291 7.18202 -0.139768 -0.071178 0.110893 11.84745 6.30331 8.87464 1.261342 -0.792366 -0.456727 4.97235 10.33478 7.56696 0.162417 -0.161842 0.060166 6.12856 9.98659 8.86367 -0.071807 -0.074850 -0.198924 6.62134 9.81948 7.16115 -0.170526 -0.038349 0.134124 8.84472 9.26711 8.37722 0.071730 -0.226053 0.058325 9.03786 8.32984 9.87873 0.045172 -0.003256 -0.227635 10.29831 8.26469 8.62560 -0.240579 0.007685 0.014401 6.52280 6.66983 9.61787 0.002502 0.002296 -0.220615 6.22536 5.65364 8.20849 0.049423 0.212957 0.057213 6.55739 6.66124 8.51468 0.135738 -0.079449 0.028382 4.10701 8.04119 5.40671 0.113045 -0.121147 0.150405 2.80136 8.46303 8.95114 0.118153 -0.112313 -0.000458 8.05254 6.11034 5.40360 0.032158 0.136274 0.091484 10.29688 4.98884 8.16670 -0.121564 1.555380 0.563899 5.79174 9.66363 7.86688 0.021172 0.001878 0.007496 9.21094 8.31702 8.79177 0.052763 -0.079191 0.087520 8.34431 6.86876 8.01094 -0.145405 0.308864 -0.422391 5.27262 7.86042 7.89665 -0.406686 -0.110296 -0.120700 4.91712 7.35264 6.35245 -0.146977 0.418839 0.209000 3.98393 7.67429 8.93065 0.084583 0.324722 -0.364244 8.55808 7.02736 6.37788 -0.195193 -0.473855 0.320588 9.06190 5.42787 8.53407 1.294152 -0.367198 -0.413184 ----------------------------------------------------------------------------------- total drift: 0.006113 -0.020903 0.007734 -------------------------------------------------------------------------------------------------------- FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV) --------------------------------------------------- free energy TOTEN = -222.5723864647 eV energy without entropy= -222.5159675225 energy(sigma->0) = -222.55358015 -------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------- writing wavefunctions volume of typ 1: 0.0 % volume of typ 2: 0.1 % volume of typ 3: 0.3 % volume of typ 4: 0.3 % volume of typ 5: 0.1 % total charge # of ion s p d tot ------------------------------------------ 1 0.996 2.052 0.016 3.064 2 0.169 0.004 0.000 0.174 3 0.162 0.004 0.000 0.166 4 0.167 0.002 0.000 0.169 5 0.167 0.002 0.000 0.169 6 0.167 0.002 0.000 0.169 7 0.166 0.002 0.000 0.169 8 0.166 0.002 0.000 0.169 9 0.166 0.002 0.000 0.169 10 0.167 0.002 0.000 0.169 11 0.166 0.002 0.000 0.169 12 0.165 0.002 0.000 0.167 13 0.165 0.002 0.000 0.168 14 0.163 0.002 0.000 0.166 15 0.086 0.000 0.000 0.086 16 0.161 0.002 0.000 0.163 17 0.161 0.002 0.000 0.163 18 0.161 0.002 0.000 0.163 19 0.161 0.002 0.000 0.163 20 0.161 0.002 0.000 0.163 21 0.159 0.002 0.000 0.161 22 0.158 0.002 0.000 0.160 23 0.157 0.002 0.000 0.160 24 0.667 1.492 0.014 2.172 25 0.663 1.439 0.043 2.146 26 0.663 1.437 0.043 2.143 27 0.664 1.438 0.042 2.143 28 0.681 1.340 0.042 2.063 29 0.667 1.498 0.018 2.183 30 0.668 1.502 0.018 2.187 31 0.665 0.951 0.317 1.933 32 0.659 0.953 0.328 1.940 33 1.237 2.943 0.012 4.192 34 1.238 2.942 0.012 4.192 35 1.239 2.938 0.012 4.189 36 1.245 2.933 0.015 4.192 -------------------------------------------------- tot 15.47 25.91 0.93 42.31 total amount of memory used by VASP MPI-rank0 479840. kBytes ======================================================================= base : 30000. kBytes nonlr-proj: 11866. kBytes fftplans : 87706. kBytes grid : 227924. kBytes one-center: 110. kBytes HF : 23. kBytes wavefun : 122211. kBytes General timing and accounting informations for this job: ======================================================== Total CPU time used (sec): 291.613 User time (sec): 274.280 System time (sec): 17.333 Elapsed time (sec): 294.947 Maximum memory used (kb): 1535940. Average memory used (kb): N/A Minor page faults: 362126 Major page faults: 0 Voluntary context switches: 25483