Program PWSCF v.5.3.0 (svn rev. 11974) starts on 30Jan2017 at 11:17:55 This program is part of the open-source Quantum ESPRESSO suite for quantum simulation of materials; please cite "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009); URL http://www.quantum-espresso.org", in publications or presentations arising from this work. More details at http://www.quantum-espresso.org/quote Parallel version (MPI), running on 36 processors R & G space division: proc/nbgrp/npool/nimage = 36 Waiting for input... Reading input from standard input Current dimensions of program PWSCF are: Max number of different atomic species (ntypx) = 10 Max number of k-points (npk) = 40000 Max angular momentum in pseudopotentials (lmaxx) = 3 file K.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 3S renormalized file I.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 5S renormalized Subspace diagonalization in iterative solution of the eigenvalue problem: one sub-group per k-point group (pool) will be used scalapack distributed-memory algorithm (size of sub-group: 4* 4 procs) Parallelization info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Min 36 21 6 1004 447 76 Max 37 22 7 1011 464 79 Sum 1321 777 241 36257 16375 2801 bravais-lattice index = 14 lattice parameter (alat) = 7.7346 a.u. unit-cell volume = 462.7237 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 16.00 number of Kohn-Sham states= 24 kinetic-energy cutoff = 41.0000 Ry charge density cutoff = 278.0000 Ry convergence threshold = 1.0E-08 mixing beta = 0.7000 number of iterations used = 8 plain mixing Exchange-correlation = PBE ( 1 4 3 4 0 0) Non magnetic calculation with spin-orbit celldm(1)= 7.734648 celldm(2)= 1.000000 celldm(3)= 1.000000 celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000 crystal axes: (cart. coord. in units of alat) a(1) = ( 1.000000 0.000000 0.000000 ) a(2) = ( 0.000000 1.000000 0.000000 ) a(3) = ( 0.000000 0.000000 1.000000 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 1.000000 -0.000000 -0.000000 ) b(2) = ( 0.000000 1.000000 -0.000000 ) b(3) = ( 0.000000 0.000000 1.000000 ) PseudoPot. # 1 for K read from file: /users/gautes/Pseudo/K.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 3dd1b7807415f540a701f2a402ca1bb0 Pseudo is Ultrasoft + core correction, Zval = 9.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1165 points, 10 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 l(7) = 2 l(8) = 2 l(9) = 2 l(10) = 2 Q(r) pseudized with 0 coefficients PseudoPot. # 2 for I read from file: /users/gautes/Pseudo/I.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 7c65668ffc3f6431347fa26e98076fd5 Pseudo is Ultrasoft + core correction, Zval = 7.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1247 points, 10 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 l(7) = 2 l(8) = 2 l(9) = 2 l(10) = 2 Q(r) pseudized with 0 coefficients atomic species valence mass pseudopotential K 9.00 39.09830 K( 1.00) I 7.00 126.90450 I( 1.00) 48 Sym. Ops., with inversion, found s frac. trans. isym = 1 identity cryst. s( 1) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 2 180 deg rotation - cart. axis [0,0,1] cryst. s( 2) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 3 180 deg rotation - cart. axis [0,1,0] cryst. s( 3) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s( 3) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 4 180 deg rotation - cart. axis [1,0,0] cryst. s( 4) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s( 4) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 5 180 deg rotation - cart. axis [1,1,0] cryst. s( 5) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s( 5) = ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 6 180 deg rotation - cart. axis [1,-1,0] cryst. s( 6) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s( 6) = ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 7 90 deg rotation - cart. axis [0,0,-1] cryst. s( 7) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s( 7) = ( -0.0000000 1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 8 90 deg rotation - cart. axis [0,0,1] cryst. s( 8) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s( 8) = ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 9 180 deg rotation - cart. axis [1,0,1] cryst. s( 9) = ( 0 0 1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s( 9) = ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 10 180 deg rotation - cart. axis [-1,0,1] cryst. s(10) = ( 0 0 -1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(10) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 11 90 deg rotation - cart. axis [0,1,0] cryst. s(11) = ( 0 0 -1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(11) = ( -0.0000000 0.0000000 1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 12 90 deg rotation - cart. axis [0,-1,0] cryst. s(12) = ( 0 0 1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(12) = ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 13 180 deg rotation - cart. axis [0,1,1] cryst. s(13) = ( -1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(13) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 14 180 deg rotation - cart. axis [0,1,-1] cryst. s(14) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(14) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 15 90 deg rotation - cart. axis [-1,0,0] cryst. s(15) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(15) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 16 90 deg rotation - cart. axis [1,0,0] cryst. s(16) = ( 1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(16) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 17 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(17) = ( 0 0 1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(17) = ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 18 120 deg rotation - cart. axis [-1,1,1] cryst. s(18) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(18) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 19 120 deg rotation - cart. axis [1,1,-1] cryst. s(19) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(19) = ( -0.0000000 1.0000000 -0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 20 120 deg rotation - cart. axis [1,-1,1] cryst. s(20) = ( 0 0 1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(20) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 21 120 deg rotation - cart. axis [1,1,1] cryst. s(21) = ( 0 1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(21) = ( 0.0000000 0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 22 120 deg rotation - cart. axis [-1,1,-1] cryst. s(22) = ( 0 -1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(22) = ( -0.0000000 -0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 23 120 deg rotation - cart. axis [1,-1,-1] cryst. s(23) = ( 0 -1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(23) = ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 24 120 deg rotation - cart. axis [-1,-1,1] cryst. s(24) = ( 0 1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(24) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 25 inversion cryst. s(25) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s(25) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 26 inv. 180 deg rotation - cart. axis [0,0,1] cryst. s(26) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s(26) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 27 inv. 180 deg rotation - cart. axis [0,1,0] cryst. s(27) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s(27) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 28 inv. 180 deg rotation - cart. axis [1,0,0] cryst. s(28) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s(28) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 29 inv. 180 deg rotation - cart. axis [1,1,0] cryst. s(29) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s(29) = ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 30 inv. 180 deg rotation - cart. axis [1,-1,0] cryst. s(30) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s(30) = ( 0.0000000 1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 31 inv. 90 deg rotation - cart. axis [0,0,-1] cryst. s(31) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s(31) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 32 inv. 90 deg rotation - cart. axis [0,0,1] cryst. s(32) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s(32) = ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 33 inv. 180 deg rotation - cart. axis [1,0,1] cryst. s(33) = ( 0 0 -1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(33) = ( -0.0000000 0.0000000 -1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 34 inv. 180 deg rotation - cart. axis [-1,0,1] cryst. s(34) = ( 0 0 1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(34) = ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 35 inv. 90 deg rotation - cart. axis [0,1,0] cryst. s(35) = ( 0 0 1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(35) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 36 inv. 90 deg rotation - cart. axis [0,-1,0] cryst. s(36) = ( 0 0 -1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s(36) = ( -0.0000000 -0.0000000 1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 37 inv. 180 deg rotation - cart. axis [0,1,1] cryst. s(37) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(37) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 38 inv. 180 deg rotation - cart. axis [0,1,-1] cryst. s(38) = ( 1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(38) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 39 inv. 90 deg rotation - cart. axis [-1,0,0] cryst. s(39) = ( -1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(39) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 40 inv. 90 deg rotation - cart. axis [1,0,0] cryst. s(40) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(40) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 41 inv. 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(41) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(41) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 42 inv. 120 deg rotation - cart. axis [-1,1,1] cryst. s(42) = ( 0 0 1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(42) = ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 43 inv. 120 deg rotation - cart. axis [1,1,-1] cryst. s(43) = ( 0 0 1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(43) = ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 44 inv. 120 deg rotation - cart. axis [1,-1,1] cryst. s(44) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(44) = ( -0.0000000 1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 45 inv. 120 deg rotation - cart. axis [1,1,1] cryst. s(45) = ( 0 -1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(45) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 46 inv. 120 deg rotation - cart. axis [-1,1,-1] cryst. s(46) = ( 0 1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(46) = ( 0.0000000 0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 47 inv. 120 deg rotation - cart. axis [1,-1,-1] cryst. s(47) = ( 0 1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(47) = ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 48 inv. 120 deg rotation - cart. axis [-1,-1,1] cryst. s(48) = ( 0 -1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(48) = ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) double point group O_h (m-3m) there are 16 classes and 6 irreducible representations the character table: E -E 8C3 -8C3 3C2 6C4 -6C4 6C2' i -i 8S6 -8S6 -3C2 -6C2' G_6+ 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 2.00 -2.00 1.00 -1.00 G_7+ 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 2.00 -2.00 1.00 -1.00 G_8+ 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 4.00 -4.00 -1.00 1.00 G_6- 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 -2.00 2.00 -1.00 1.00 G_7- 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 -2.00 2.00 -1.00 1.00 G_8- 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 -4.00 4.00 1.00 -1.00 3s_h 6S4 -6S4 6s_d -3s_h -6s_d G_6+ 0.00 1.41 -1.41 0.00 G_7+ 0.00 -1.41 1.41 0.00 G_8+ 0.00 0.00 0.00 0.00 G_6- 0.00 -1.41 1.41 0.00 G_7- 0.00 1.41 -1.41 0.00 G_8- 0.00 0.00 0.00 0.00 the symmetry operations in each class and the name of the first element: E 1 identity -E -1 identity E 8C3 17 19 20 18 24 21 22 23 120 deg rotation - cart. axis [-1,-1,-1] -8C3 -17 -19 -20 -18 -24 -21 -22 -23 120 deg rotation - cart. axis [-1,-1,-1] E 3C2 -3C2 2 -2 4 -4 3 -3 180 deg rotation - cart. axis [0,0,1] 6C4 7 8 15 16 12 11 90 deg rotation - cart. axis [0,0,-1] -6C4 -7 -8 -15 -16 -12 -11 90 deg rotation - cart. axis [0,0,-1] E 6C2'-6C2' 5 -5 6 -6 14 -13 -14 13 -10 -9 10 9 180 deg rotation - cart. axis [1,1,0] i 25 inversion -i -25 inversion E 8S6 41 43 44 42 48 45 46 47 inv. 120 deg rotation - cart. axis [-1,-1,-1] -8S6 -41 -43 -44 -42 -48 -45 -46 -47 inv. 120 deg rotation - cart. axis [-1,-1,-1] E 3s_h-3s_h 26 -26 28 -28 27 -27 inv. 180 deg rotation - cart. axis [0,0,1] 6S4 31 32 39 40 36 35 inv. 90 deg rotation - cart. axis [0,0,-1] -6S4 -31 -32 -39 -40 -36 -35 inv. 90 deg rotation - cart. axis [0,0,-1] E 6s_d-6s_d 29 -29 30 -30 38 -37 -38 37 -34 -33 34 33 inv. 180 deg rotation - cart. axis [1,1,0] Cartesian axes number of k points= 20 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 Dense grid: 36257 G-vectors FFT dimensions: ( 45, 45, 45) Smooth grid: 16375 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.05 Mb ( 124, 24) NL pseudopotentials 0.06 Mb ( 62, 68) Each V/rho on FFT grid 0.06 Mb ( 4050) Each G-vector array 0.01 Mb ( 1006) G-vector shells 0.00 Mb ( 253) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.18 Mb ( 124, 96) Each subspace H/S matrix 0.01 Mb ( 24, 24) Each matrix 0.05 Mb ( 68, 2, 24) Arrays for rho mixing 0.49 Mb ( 4050, 8) Check: negative/imaginary core charge= -0.000002 0.000000 Initial potential from superposition of free atoms starting charge 15.98109, renormalised to 16.00000 Starting wfc are 24 randomized atomic wfcs total cpu time spent up to now is 3.2 secs per-process dynamical memory: 28.7 Mb Self-consistent Calculation iteration # 1 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 6.5 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.29E-04, avg # of iterations = 4.8 total cpu time spent up to now is 4.7 secs total energy = -92.35570894 Ry Harris-Foulkes estimate = -92.36631166 Ry estimated scf accuracy < 0.01960508 Ry iteration # 2 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.23E-04, avg # of iterations = 3.6 total cpu time spent up to now is 5.4 secs total energy = -92.36039095 Ry Harris-Foulkes estimate = -92.36151046 Ry estimated scf accuracy < 0.00245910 Ry iteration # 3 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.54E-05, avg # of iterations = 3.1 total cpu time spent up to now is 6.0 secs total energy = -92.36096048 Ry Harris-Foulkes estimate = -92.36090141 Ry estimated scf accuracy < 0.00010494 Ry iteration # 4 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.56E-07, avg # of iterations = 6.0 total cpu time spent up to now is 6.9 secs total energy = -92.36098135 Ry Harris-Foulkes estimate = -92.36098108 Ry estimated scf accuracy < 0.00000275 Ry iteration # 5 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.72E-08, avg # of iterations = 2.3 total cpu time spent up to now is 7.4 secs total energy = -92.36098220 Ry Harris-Foulkes estimate = -92.36098182 Ry estimated scf accuracy < 0.00000029 Ry iteration # 6 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.80E-09, avg # of iterations = 2.4 total cpu time spent up to now is 8.0 secs total energy = -92.36098233 Ry Harris-Foulkes estimate = -92.36098235 Ry estimated scf accuracy < 0.00000003 Ry iteration # 7 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.04E-10, avg # of iterations = 2.0 total cpu time spent up to now is 8.5 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2103 PWs) bands (ev): -25.1046 -25.1046 -8.7417 -8.7417 -8.4749 -8.4749 -8.4749 -8.4749 -7.3761 -7.3761 3.0140 3.0140 4.0704 4.0704 4.0704 4.0704 7.9296 7.9296 8.6055 8.6055 8.6055 8.6055 12.6805 12.6954 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.1667 ( 2032 PWs) bands (ev): -25.1019 -25.1019 -8.7772 -8.7772 -8.5339 -8.5339 -8.4695 -8.4695 -7.2297 -7.2297 2.6359 2.6359 3.6434 3.6434 4.0836 4.0836 8.4897 8.4897 8.8272 8.8272 9.5425 9.5425 12.5924 12.5925 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.3333 ( 2040 PWs) bands (ev): -25.0973 -25.0973 -8.8805 -8.8805 -8.5821 -8.5821 -8.4590 -8.4590 -6.9634 -6.9634 1.6814 1.6814 3.4881 3.4881 4.1089 4.1089 9.4192 9.4192 9.4740 9.4740 11.2865 11.2865 12.2331 12.2332 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000-0.5000 ( 2040 PWs) bands (ev): -25.0950 -25.0950 -8.9374 -8.9374 -8.5896 -8.5896 -8.4537 -8.4537 -6.8349 -6.8349 1.2478 1.2478 3.4757 3.4757 4.1210 4.1210 10.0662 10.0662 10.2175 10.2175 11.1052 11.1052 11.9249 11.9249 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 2037 PWs) bands (ev): -25.0997 -25.0997 -8.7955 -8.7955 -8.5584 -8.5584 -8.4829 -8.4829 -7.1292 -7.1292 2.5033 2.5033 3.2971 3.2971 3.9333 3.9333 8.9150 8.9150 9.4806 9.4806 9.8444 9.8444 12.0984 12.0986 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 2045 PWs) bands (ev): -25.0951 -25.0951 -8.8608 -8.8608 -8.5945 -8.5945 -8.4760 -8.4760 -6.9369 -6.9369 1.7871 1.7871 3.1264 3.1264 3.9748 3.9748 9.6094 9.6094 9.7526 9.7526 11.0794 11.0794 11.8413 11.8413 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667-0.5000 ( 2052 PWs) bands (ev): -25.0929 -25.0929 -8.9029 -8.9029 -8.6002 -8.6002 -8.4693 -8.4693 -6.8419 -6.8419 1.3978 1.3978 3.1559 3.1559 4.0115 4.0115 9.7631 9.7631 10.1877 10.1877 11.4287 11.4287 11.9572 11.9587 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 2051 PWs) bands (ev): -25.0907 -25.0907 -8.8428 -8.8428 -8.6179 -8.6179 -8.4845 -8.4845 -6.8847 -6.8847 1.8580 1.8580 2.5104 2.5104 4.0599 4.0599 8.9459 8.9459 9.9092 9.9092 11.2378 11.2378 11.9328 11.9328 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 2058 PWs) bands (ev): -25.0886 -25.0886 -8.8495 -8.8495 -8.6195 -8.6195 -8.4808 -8.4808 -6.8585 -6.8585 1.6803 1.6803 2.4832 2.4832 4.1347 4.1347 8.4995 8.4995 10.0137 10.0137 11.1711 11.1711 11.3395 11.3396 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 2048 PWs) bands (ev): -25.0863 -25.0863 -8.8385 -8.8385 -8.6180 -8.6180 -8.4786 -8.4786 -6.8684 -6.8684 1.7210 1.7210 2.3137 2.3137 4.2361 4.2361 7.9292 7.9292 9.9050 9.9050 10.8899 10.8900 10.9003 10.9003 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 2064 PWs) bands (ev): -25.0975 -25.0975 -8.8045 -8.8045 -8.5660 -8.5660 -8.5041 -8.5041 -7.0546 -7.0546 2.4400 2.4400 3.2714 3.2714 3.4655 3.4655 9.6722 9.6722 9.8997 9.8997 9.9363 9.9363 11.9412 11.9603 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.3333 ( 2040 PWs) bands (ev): -25.0929 -25.0929 -8.8462 -8.8462 -8.5945 -8.5945 -8.4993 -8.4993 -6.9073 -6.9073 1.8722 1.8722 3.0469 3.0469 3.5311 3.5311 9.9798 9.9798 10.1616 10.1616 10.7259 10.7259 11.7487 11.7488 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 2054 PWs) bands (ev): -25.0907 -25.0907 -8.8775 -8.8775 -8.6015 -8.6015 -8.4901 -8.4901 -6.8332 -6.8332 1.5115 1.5115 3.0565 3.0565 3.6212 3.6212 10.0010 10.0010 10.0640 10.0640 10.7750 10.7750 11.9610 11.9611 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333 0.3333 ( 2054 PWs) bands (ev): -25.0886 -25.0886 -8.8353 -8.8353 -8.6045 -8.6045 -8.5141 -8.5141 -6.8467 -6.8467 1.8897 1.8897 2.5508 2.5508 3.4571 3.4571 9.2346 9.2346 10.7041 10.7041 10.8484 10.8484 10.9240 10.9240 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 2062 PWs) bands (ev): -25.0865 -25.0865 -8.8420 -8.8420 -8.6082 -8.6082 -8.5088 -8.5088 -6.8154 -6.8154 1.7016 1.7016 2.4937 2.4937 3.5245 3.5245 8.7880 8.7880 10.4620 10.4620 10.7186 10.7186 11.2899 11.2899 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667-0.5000-0.5000 ( 2052 PWs) bands (ev): -25.0843 -25.0843 -8.8364 -8.8364 -8.6085 -8.6085 -8.5072 -8.5072 -6.8061 -6.8061 1.6962 1.6962 2.3100 2.3100 3.5692 3.5692 8.2982 8.2982 10.5706 10.5706 10.5737 10.5737 10.9827 10.9827 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 2064 PWs) bands (ev): -25.0844 -25.0844 -8.8307 -8.8307 -8.5821 -8.5821 -8.5567 -8.5567 -6.7664 -6.7664 1.8376 1.8376 2.6440 2.6440 2.6737 2.6737 9.6023 9.6023 10.0859 10.0859 10.0948 10.0948 12.2728 12.2729 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 2046 PWs) bands (ev): -25.0822 -25.0822 -8.8351 -8.8351 -8.5886 -8.5886 -8.5528 -8.5528 -6.7235 -6.7235 1.6871 1.6871 2.4538 2.4538 2.7316 2.7316 9.3083 9.3083 9.8604 9.8604 10.2986 10.2986 12.7048 12.7049 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333-0.5000-0.5000 ( 2048 PWs) bands (ev): -25.0802 -25.0802 -8.8371 -8.8371 -8.5894 -8.5894 -8.5582 -8.5582 -6.6746 -6.6746 1.6040 1.6040 2.3013 2.3013 2.6262 2.6262 9.2240 9.2240 10.0752 10.0752 10.0844 10.0844 12.7410 12.7410 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k =-0.5000-0.5000-0.5000 ( 2008 PWs) bands (ev): -25.0779 -25.0779 -8.8409 -8.8409 -8.5797 -8.5797 -8.5797 -8.5797 -6.6047 -6.6047 1.5053 1.5053 2.2962 2.2962 2.2962 2.2962 9.8667 9.8667 9.8886 9.8886 9.8886 9.8886 13.2284 13.2284 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 6.0580 ev ! total energy = -92.36098234 Ry Harris-Foulkes estimate = -92.36098234 Ry estimated scf accuracy < 2.1E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -20.22198578 Ry hartree contribution = 15.78644594 Ry xc contribution = -27.17379926 Ry ewald contribution = -60.75164323 Ry smearing contrib. (-TS) = 0.00000000 Ry convergence has been achieved in 7 iterations Writing output data file KI.save init_run : 0.70s CPU 0.75s WALL ( 1 calls) electrons : 5.19s CPU 5.37s WALL ( 1 calls) Called by init_run: wfcinit : 0.13s CPU 0.14s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 3.99s CPU 4.13s WALL ( 8 calls) sum_band : 0.88s CPU 0.88s WALL ( 8 calls) v_of_rho : 0.02s CPU 0.02s WALL ( 8 calls) v_h : 0.00s CPU 0.00s WALL ( 8 calls) v_xc : 0.02s CPU 0.02s WALL ( 8 calls) newd : 0.30s CPU 0.32s WALL ( 8 calls) mix_rho : 0.01s CPU 0.01s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.00s CPU 0.01s WALL ( 340 calls) cegterg : 3.80s CPU 3.89s WALL ( 160 calls) Called by sum_band: sum_band:bec : 0.36s CPU 0.32s WALL ( 160 calls) addusdens : 0.19s CPU 0.20s WALL ( 8 calls) Called by *egterg: h_psi : 1.73s CPU 1.83s WALL ( 794 calls) s_psi : 0.14s CPU 0.11s WALL ( 794 calls) g_psi : 0.02s CPU 0.01s WALL ( 614 calls) cdiaghg : 1.71s CPU 1.75s WALL ( 754 calls) cegterg:over : 0.10s CPU 0.10s WALL ( 614 calls) cegterg:upda : 0.08s CPU 0.08s WALL ( 614 calls) cegterg:last : 0.04s CPU 0.03s WALL ( 161 calls) cdiaghg:chol : 0.07s CPU 0.09s WALL ( 754 calls) cdiaghg:inve : 0.02s CPU 0.02s WALL ( 754 calls) cdiaghg:para : 0.12s CPU 0.13s WALL ( 1508 calls) Called by h_psi: h_psi:vloc : 1.49s CPU 1.58s WALL ( 794 calls) h_psi:vnl : 0.23s CPU 0.24s WALL ( 794 calls) add_vuspsi : 0.12s CPU 0.14s WALL ( 794 calls) General routines calbec : 0.15s CPU 0.13s WALL ( 954 calls) fft : 0.06s CPU 0.05s WALL ( 242 calls) ffts : 0.01s CPU 0.00s WALL ( 64 calls) fftw : 1.58s CPU 1.70s WALL ( 51728 calls) interpolate : 0.03s CPU 0.02s WALL ( 64 calls) Parallel routines fft_scatter : 1.01s CPU 1.09s WALL ( 52034 calls) PWSCF : 7.93s CPU 9.39s WALL This run was terminated on: 11:18: 5 30Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=