Program PWSCF v.5.3.0 (svn rev. 11974) starts on 24Dec2016 at 18:31:38 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 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 52 52 15 1706 1706 261 Max 53 53 16 1709 1709 264 Sum 1885 1885 545 61445 61445 9435 bravais-lattice index = 14 lattice parameter (alat) = 7.5778 a.u. unit-cell volume = 435.1406 (a.u.)^3 number of atoms/cell = 4 number of atomic types = 2 number of electrons = 12.00 number of Kohn-Sham states= 20 kinetic-energy cutoff = 103.0000 Ry charge density cutoff = 412.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.577801 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 Li read from file: /users/gautes/Pseudo/Li.rel-pbe-s-rrkjus_psl.1.0.0.UPF MD5 check sum: 9da6379fed0e9079eb53e17c51ce5ca0 Pseudo is Ultrasoft, Zval = 3.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1017 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 Q(r) pseudized with 0 coefficients PseudoPot. # 2 for Al read from file: /users/gautes/Pseudo/Al.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: d344aa71f13a435d235890da883ec65e Pseudo is Ultrasoft + core correction, Zval = 3.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1135 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 Li 3.00 6.94100 Li( 1.00) Al 3.00 26.98150 Al( 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: 61445 G-vectors FFT dimensions: ( 50, 50, 50) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.13 Mb ( 440, 20) NL pseudopotentials 0.39 Mb ( 220, 116) Each V/rho on FFT grid 0.08 Mb ( 5000) Each G-vector array 0.01 Mb ( 1709) G-vector shells 0.00 Mb ( 372) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.54 Mb ( 440, 80) Each subspace H/S matrix 0.01 Mb ( 20, 20) Each matrix 0.07 Mb ( 116, 2, 20) Arrays for rho mixing 0.61 Mb ( 5000, 8) Initial potential from superposition of free atoms starting charge 11.98860, renormalised to 12.00000 Starting wfc are 34 randomized atomic wfcs total cpu time spent up to now is 2.2 secs per-process dynamical memory: 23.3 Mb Self-consistent Calculation iteration # 1 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 3.6 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.46E-04, avg # of iterations = 3.6 total cpu time spent up to now is 5.0 secs total energy = -29.99859248 Ry Harris-Foulkes estimate = -30.00279279 Ry estimated scf accuracy < 0.01805302 Ry iteration # 2 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.50E-04, avg # of iterations = 1.1 total cpu time spent up to now is 6.1 secs total energy = -29.99898670 Ry Harris-Foulkes estimate = -29.99920975 Ry estimated scf accuracy < 0.00127023 Ry iteration # 3 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.06E-05, avg # of iterations = 2.4 total cpu time spent up to now is 7.3 secs total energy = -29.99909062 Ry Harris-Foulkes estimate = -29.99910050 Ry estimated scf accuracy < 0.00002979 Ry iteration # 4 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.48E-07, avg # of iterations = 3.1 total cpu time spent up to now is 8.7 secs total energy = -29.99910606 Ry Harris-Foulkes estimate = -29.99910812 Ry estimated scf accuracy < 0.00000438 Ry iteration # 5 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.65E-08, avg # of iterations = 2.1 total cpu time spent up to now is 9.9 secs total energy = -29.99910674 Ry Harris-Foulkes estimate = -29.99910715 Ry estimated scf accuracy < 0.00000077 Ry iteration # 6 ecut= 103.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.43E-09, avg # of iterations = 2.0 total cpu time spent up to now is 11.1 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 7689 PWs) bands (ev): -39.5732 -39.5732 -3.5010 -3.5010 4.4952 4.4952 4.4952 4.4952 4.5044 4.5044 5.7625 5.7625 5.7625 5.7625 7.7794 7.7794 12.7663 12.7663 12.7664 12.7665 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.1667 ( 7705 PWs) bands (ev): -39.5731 -39.5731 -3.2696 -3.2696 2.6074 2.6074 4.7266 4.7266 4.7330 4.7330 6.0121 6.0121 7.0204 7.0204 8.9574 8.9574 10.4947 10.4947 10.4948 10.4948 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.3333 ( 7712 PWs) bands (ev): -39.5730 -39.5730 -2.6056 -2.6056 0.5917 0.5917 5.4082 5.4082 5.4163 5.4163 6.7528 6.7528 7.4502 7.4502 8.3975 8.3975 8.3989 8.3989 9.6201 9.6201 occupation numbers 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 0.0000 0.0000 k = 0.0000 0.0000-0.5000 ( 7672 PWs) bands (ev): -39.5729 -39.5729 -1.9243 -1.9243 -0.5495 -0.5495 6.2182 6.2182 6.2363 6.2363 7.1033 7.1033 7.1147 7.1147 7.6035 7.6035 7.8284 7.8284 8.1239 8.1239 occupation numbers 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 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 7699 PWs) bands (ev): -39.5731 -39.5731 -3.0373 -3.0373 2.5984 2.5984 3.1124 3.1124 4.9612 4.9612 6.6716 6.6716 8.5423 8.5423 8.6621 8.6621 9.4668 9.4668 10.5742 10.5742 occupation numbers 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 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 7700 PWs) bands (ev): -39.5729 -39.5729 -2.3688 -2.3688 0.7946 0.7946 3.5529 3.5529 5.6430 5.6430 6.3507 6.3507 7.7048 7.7048 8.6268 8.6268 9.2799 9.2799 10.8071 10.8071 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9947 0.9947 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667-0.5000 ( 7686 PWs) bands (ev): -39.5728 -39.5728 -1.6735 -1.6735 -0.3532 -0.3532 4.2551 4.2551 5.3326 5.3326 6.4559 6.4559 7.3418 7.3418 8.0961 8.0961 9.2549 9.2549 10.2717 10.2717 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0754 0.0754 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 7635 PWs) bands (ev): -39.5728 -39.5728 -1.6819 -1.6819 1.1276 1.1276 1.7856 1.7856 4.5820 4.5820 6.3229 6.3229 8.2425 8.2425 9.1804 9.1804 9.4272 9.4272 10.6114 10.6114 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9993 0.9993 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 7680 PWs) bands (ev): -39.5727 -39.5727 -0.9287 -0.9287 0.1820 0.1820 2.1042 2.1042 3.5710 3.5710 7.1328 7.1328 8.0138 8.0138 8.8341 8.8341 10.0068 10.0068 10.0283 10.0284 occupation numbers 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 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 7692 PWs) bands (ev): -39.5726 -39.5726 0.0689 0.0689 0.6275 0.6275 0.6308 0.6308 2.7038 2.7038 7.9782 7.9782 8.7157 8.7157 8.7158 8.7158 9.8692 9.8692 10.2180 10.2239 occupation numbers 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 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 7679 PWs) bands (ev): -39.5730 -39.5730 -2.8041 -2.8041 2.8407 2.8407 2.8444 2.8444 3.6219 3.6219 8.4279 8.4279 8.4281 8.4281 8.9975 8.9975 9.0189 9.0189 9.0248 9.0248 occupation numbers 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 0.0000 0.0000 k = 0.1667 0.1667 0.3333 ( 7668 PWs) bands (ev): -39.5728 -39.5728 -2.1316 -2.1316 1.0048 1.0048 3.5687 3.5687 3.9785 3.9785 6.4963 6.4963 7.1741 7.1741 9.5572 9.5572 9.5678 9.5678 10.0579 10.0579 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0042 0.0042 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 7678 PWs) bands (ev): -39.5728 -39.5728 -1.4231 -1.4231 -0.1517 -0.1517 4.4255 4.4255 4.5683 4.5683 5.0381 5.0381 6.3280 6.3280 10.1420 10.1420 10.4130 10.4130 10.6284 10.6284 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 0.9990 0.9990 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333 0.3333 ( 7670 PWs) bands (ev): -39.5727 -39.5727 -1.4427 -1.4427 1.3751 1.3751 1.9669 1.9669 4.2443 4.2443 5.0345 5.0345 7.2098 7.2098 7.9411 7.9411 10.3508 10.3508 10.9851 10.9851 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 7688 PWs) bands (ev): -39.5726 -39.5726 -0.6800 -0.6800 0.4017 0.4017 2.3477 2.3477 3.6504 3.6504 5.2076 5.2076 6.4206 6.4206 8.1270 8.1270 9.9021 9.9021 11.1885 11.1885 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 0.5210 0.5210 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667-0.5000-0.5000 ( 7672 PWs) bands (ev): -39.5726 -39.5726 0.3140 0.3140 0.8759 0.8759 0.8791 0.8791 2.7933 2.7933 6.1346 6.1346 6.7379 6.7379 6.7405 6.7405 9.1426 9.1426 12.1359 12.1360 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 7679 PWs) bands (ev): -39.5726 -39.5726 -0.7443 -0.7443 2.1096 2.1096 2.1126 2.1126 2.7527 2.7527 5.0836 5.0836 5.0868 5.0868 6.3802 6.3802 9.2328 9.2328 13.6122 13.6142 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 0.9551 0.9551 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 7708 PWs) bands (ev): -39.5725 -39.5725 0.0580 0.0580 1.0214 1.0214 3.0200 3.0200 3.1256 3.1256 3.5996 3.5996 4.8568 4.8568 6.0327 6.0327 8.5330 8.5330 14.5294 14.5295 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 k = 0.3333-0.5000-0.5000 ( 7688 PWs) bands (ev): -39.5724 -39.5724 1.0396 1.0396 1.6145 1.6145 1.6173 1.6173 2.9982 2.9982 4.0357 4.0357 4.6028 4.6028 4.6062 4.6062 7.9622 7.9622 15.4521 15.4523 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 k =-0.5000-0.5000-0.5000 ( 7664 PWs) bands (ev): -39.5724 -39.5724 1.9781 1.9781 2.6287 2.6287 2.6287 2.6287 2.6287 2.6287 3.1151 3.1151 3.1151 3.1151 3.1247 3.1247 7.4870 7.4870 18.2920 18.2921 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 the Fermi energy is 6.4218 ev ! total energy = -29.99910687 Ry Harris-Foulkes estimate = -29.99910687 Ry estimated scf accuracy < 3.3E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -0.85306113 Ry hartree contribution = 4.00790807 Ry xc contribution = -11.37235297 Ry ewald contribution = -21.78138882 Ry smearing contrib. (-TS) = -0.00021203 Ry convergence has been achieved in 6 iterations Writing output data file LiAl3.save init_run : 0.76s CPU 0.84s WALL ( 1 calls) electrons : 8.77s CPU 8.95s WALL ( 1 calls) Called by init_run: wfcinit : 0.51s CPU 0.54s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 6.91s CPU 7.02s WALL ( 7 calls) sum_band : 1.53s CPU 1.56s WALL ( 7 calls) v_of_rho : 0.02s CPU 0.03s WALL ( 7 calls) v_h : 0.01s CPU 0.00s WALL ( 7 calls) v_xc : 0.01s CPU 0.02s WALL ( 7 calls) newd : 0.32s CPU 0.35s WALL ( 7 calls) mix_rho : 0.01s CPU 0.01s WALL ( 7 calls) Called by c_bands: init_us_2 : 0.04s CPU 0.04s WALL ( 300 calls) cegterg : 5.98s CPU 6.03s WALL ( 140 calls) Called by sum_band: sum_band:bec : 0.42s CPU 0.42s WALL ( 140 calls) addusdens : 0.16s CPU 0.17s WALL ( 7 calls) Called by *egterg: h_psi : 4.95s CPU 5.00s WALL ( 518 calls) s_psi : 0.19s CPU 0.19s WALL ( 518 calls) g_psi : 0.01s CPU 0.01s WALL ( 358 calls) cdiaghg : 0.84s CPU 0.88s WALL ( 478 calls) cegterg:over : 0.10s CPU 0.10s WALL ( 358 calls) cegterg:upda : 0.11s CPU 0.13s WALL ( 358 calls) cegterg:last : 0.06s CPU 0.05s WALL ( 140 calls) cdiaghg:chol : 0.04s CPU 0.07s WALL ( 478 calls) cdiaghg:inve : 0.00s CPU 0.01s WALL ( 478 calls) cdiaghg:para : 0.06s CPU 0.07s WALL ( 956 calls) Called by h_psi: h_psi:vloc : 4.47s CPU 4.50s WALL ( 518 calls) h_psi:vnl : 0.46s CPU 0.48s WALL ( 518 calls) add_vuspsi : 0.26s CPU 0.24s WALL ( 518 calls) General routines calbec : 0.28s CPU 0.31s WALL ( 658 calls) fft : 0.05s CPU 0.04s WALL ( 135 calls) fftw : 5.06s CPU 5.08s WALL ( 36596 calls) Parallel routines fft_scatter : 1.81s CPU 1.83s WALL ( 36731 calls) PWSCF : 11.20s CPU 12.14s WALL This run was terminated on: 18:31:50 24Dec2016 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=