Program PWSCF v.5.3.0 (svn rev. 11974) starts on 24Dec2016 at 19:39:45 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 Sc.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 3S 3P 3P 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 18 5 990 354 57 Max 37 19 6 995 366 62 Sum 1313 665 197 35729 12965 2109 bravais-lattice index = 14 lattice parameter (alat) = 6.4024 a.u. unit-cell volume = 262.4380 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 14.00 number of Kohn-Sham states= 22 kinetic-energy cutoff = 51.0000 Ry charge density cutoff = 402.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)= 6.402392 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 Sc read from file: /users/gautes/Pseudo/Sc.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: e27f768e72ba10e20459589e24f2b6a7 Pseudo is Ultrasoft + core correction, Zval = 11.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1173 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 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 Sc 11.00 44.95590 Sc( 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.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 Dense grid: 35729 G-vectors FFT dimensions: ( 45, 45, 45) Smooth grid: 12965 G-vectors FFT dimensions: ( 30, 30, 30) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.03 Mb ( 98, 22) NL pseudopotentials 0.05 Mb ( 49, 68) Each V/rho on FFT grid 0.06 Mb ( 4050) Each G-vector array 0.01 Mb ( 992) G-vector shells 0.00 Mb ( 254) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.13 Mb ( 98, 88) Each subspace H/S matrix 0.01 Mb ( 22, 22) Each matrix 0.05 Mb ( 68, 2, 22) Arrays for rho mixing 0.49 Mb ( 4050, 8) Check: negative/imaginary core charge= -0.000007 0.000000 Initial potential from superposition of free atoms starting charge 13.99517, renormalised to 14.00000 Starting wfc are 28 randomized atomic wfcs total cpu time spent up to now is 2.4 secs per-process dynamical memory: 28.3 Mb Self-consistent Calculation iteration # 1 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 3.1 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 6.62E-05, avg # of iterations = 4.3 total cpu time spent up to now is 3.7 secs total energy = -102.98786687 Ry Harris-Foulkes estimate = -102.99679278 Ry estimated scf accuracy < 0.01838162 Ry iteration # 2 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.31E-04, avg # of iterations = 3.8 total cpu time spent up to now is 4.3 secs total energy = -102.98969494 Ry Harris-Foulkes estimate = -102.99414880 Ry estimated scf accuracy < 0.00737186 Ry iteration # 3 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.27E-05, avg # of iterations = 3.6 total cpu time spent up to now is 4.9 secs total energy = -102.99189343 Ry Harris-Foulkes estimate = -102.99217336 Ry estimated scf accuracy < 0.00069806 Ry iteration # 4 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.99E-06, avg # of iterations = 4.6 total cpu time spent up to now is 5.6 secs total energy = -102.99204653 Ry Harris-Foulkes estimate = -102.99204614 Ry estimated scf accuracy < 0.00000188 Ry iteration # 5 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.35E-08, avg # of iterations = 5.7 total cpu time spent up to now is 6.4 secs total energy = -102.99206117 Ry Harris-Foulkes estimate = -102.99206323 Ry estimated scf accuracy < 0.00000478 Ry iteration # 6 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.35E-08, avg # of iterations = 3.7 total cpu time spent up to now is 7.0 secs total energy = -102.99206163 Ry Harris-Foulkes estimate = -102.99206175 Ry estimated scf accuracy < 0.00000023 Ry iteration # 7 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.64E-09, avg # of iterations = 4.2 total cpu time spent up to now is 7.7 secs total energy = -102.99206173 Ry Harris-Foulkes estimate = -102.99206175 Ry estimated scf accuracy < 0.00000004 Ry iteration # 8 ecut= 51.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.97E-10, avg # of iterations = 3.8 total cpu time spent up to now is 8.2 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1575 PWs) bands (ev): -40.0687 -40.0687 -18.9866 -18.9866 -18.4789 -18.4789 -18.4789 -18.4789 1.0203 1.0203 8.8100 8.8100 8.8100 8.8100 12.1013 12.1013 12.1013 12.1014 12.1221 12.1221 13.0313 13.0323 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 0.0000 0.0000 k = 0.0000 0.0000 0.1429 ( 1615 PWs) bands (ev): -40.0661 -40.0661 -18.9957 -18.9957 -18.4995 -18.4995 -18.4762 -18.4762 1.3048 1.3048 8.4565 8.4565 8.9523 8.9523 11.4139 11.4139 11.4317 11.4317 12.1487 12.1487 12.4833 12.4833 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 0.0000 0.0000 k = 0.0000 0.0000 0.2857 ( 1630 PWs) bands (ev): -40.0604 -40.0604 -19.0198 -19.0198 -18.5406 -18.5406 -18.4698 -18.4698 2.1244 2.1244 7.1194 7.1194 9.3161 9.3161 10.6700 10.6700 10.6855 10.6855 12.2419 12.2419 12.4554 12.4554 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.9989 0.9989 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4286 ( 1622 PWs) bands (ev): -40.0557 -40.0557 -19.0430 -19.0430 -18.5687 -18.5687 -18.4645 -18.4645 3.3343 3.3343 5.3969 5.3969 9.6714 9.6714 10.2629 10.2629 10.2766 10.2766 12.3193 12.3193 13.0415 13.0415 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 0.0000 0.0000 k = 0.0000 0.1429 0.1429 ( 1623 PWs) bands (ev): -40.0636 -40.0636 -19.0039 -19.0039 -18.5087 -18.5087 -18.4846 -18.4846 1.5811 1.5811 8.5090 8.5090 8.7436 8.7436 10.7788 10.7788 11.2818 11.2818 11.8766 11.8766 12.9706 12.9706 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 0.0000 0.0000 k = 0.0000 0.1429 0.2857 ( 1604 PWs) bands (ev): -40.0579 -40.0579 -19.0259 -19.0259 -18.5463 -18.5463 -18.4817 -18.4817 2.3738 2.3738 7.3031 7.3031 9.0193 9.0193 10.3373 10.3373 10.7564 10.7564 11.8776 11.8776 12.7192 12.7192 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 0.0000 0.0000 k = 0.0000 0.1429 0.4286 ( 1616 PWs) bands (ev): -40.0533 -40.0533 -19.0476 -19.0476 -18.5746 -18.5746 -18.4768 -18.4768 3.5211 3.5211 5.6467 5.6467 9.4890 9.4890 9.9621 9.9621 10.3855 10.3855 11.9809 11.9809 13.0513 13.0513 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.0028 0.0028 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.2857 ( 1610 PWs) bands (ev): -40.0522 -40.0522 -19.0434 -19.0434 -18.5713 -18.5713 -18.4927 -18.4927 3.0758 3.0758 7.6214 7.6214 7.9603 7.9603 10.1694 10.1694 10.8526 10.8526 11.3854 11.3854 12.9132 12.9132 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 0.0000 0.0000 k = 0.0000 0.2857 0.4286 ( 1613 PWs) bands (ev): -40.0477 -40.0477 -19.0611 -19.0611 -18.5962 -18.5962 -18.4923 -18.4923 4.0147 4.0147 6.3521 6.3521 8.4333 8.4333 9.8289 9.8289 10.6750 10.6750 11.3828 11.3828 13.0138 13.0138 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 0.0000 0.0000 k = 0.0000 0.4286 0.4286 ( 1624 PWs) bands (ev): -40.0432 -40.0432 -19.0755 -19.0755 -18.6185 -18.6185 -18.4959 -18.4959 4.5637 4.5637 7.1963 7.1963 7.5991 7.5991 9.1790 9.1790 10.9272 10.9272 11.0185 11.0185 13.0208 13.0208 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 k = 0.1429 0.1429 0.1429 ( 1629 PWs) bands (ev): -40.0610 -40.0610 -19.0115 -19.0115 -18.5071 -18.5071 -18.5037 -18.5037 1.8494 1.8494 8.6164 8.6164 8.6248 8.6248 10.3520 10.3520 11.2397 11.2397 11.2457 11.2457 13.7202 13.7243 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 0.0000 0.0000 k = 0.1429 0.1429 0.2857 ( 1613 PWs) bands (ev): -40.0554 -40.0554 -19.0318 -19.0318 -18.5478 -18.5478 -18.4980 -18.4980 2.6164 2.6164 7.4782 7.4782 8.9050 8.9050 10.2253 10.2253 10.4809 10.4809 11.4310 11.4310 13.0813 13.0813 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 0.0000 0.0000 k = 0.1429 0.1429 0.4286 ( 1616 PWs) bands (ev): -40.0508 -40.0508 -19.0521 -19.0521 -18.5772 -18.5772 -18.4922 -18.4922 3.7072 3.7072 5.8886 5.8886 9.2019 9.2019 10.0045 10.0045 10.4239 10.4239 11.5620 11.5620 13.1236 13.1236 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 0.0000 0.0000 k = 0.1429 0.2857 0.2857 ( 1611 PWs) bands (ev): -40.0498 -40.0498 -19.0482 -19.0482 -18.5676 -18.5676 -18.5141 -18.5141 3.2991 3.2991 7.8389 7.8389 7.9529 7.9529 10.2134 10.2134 10.3401 10.3401 11.1000 11.1000 13.3524 13.3524 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 0.0000 0.0000 k = 0.1429 0.2857 0.4286 ( 1607 PWs) bands (ev): -40.0453 -40.0453 -19.0649 -19.0649 -18.5933 -18.5933 -18.5131 -18.5131 4.2074 4.2074 6.5455 6.5455 8.2940 8.2940 9.9902 9.9902 10.6589 10.6589 11.0981 11.0981 13.1708 13.1708 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 0.0000 0.0000 k = 0.1429 0.4286 0.4286 ( 1608 PWs) bands (ev): -40.0408 -40.0408 -19.0787 -19.0787 -18.6144 -18.6144 -18.5177 -18.5177 4.7832 4.7832 7.2562 7.2562 7.4843 7.4843 9.3943 9.3943 11.0149 11.0149 11.1547 11.1547 13.0969 13.0969 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.7475 0.7475 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.2857 0.2857 ( 1604 PWs) bands (ev): -40.0442 -40.0442 -19.0619 -19.0619 -18.5612 -18.5612 -18.5567 -18.5567 3.9242 3.9242 7.9079 7.9079 7.9194 7.9194 9.5451 9.5451 10.7429 10.7429 10.7600 10.7600 13.3468 13.3468 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 0.0000 0.0000 k = 0.2857 0.2857 0.4286 ( 1619 PWs) bands (ev): -40.0398 -40.0398 -19.0763 -19.0763 -18.5912 -18.5912 -18.5521 -18.5521 4.7406 4.7406 6.9624 6.9624 7.8590 7.8590 10.2442 10.2442 10.6109 10.6109 11.0691 11.0691 12.5954 12.5954 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 0.0000 0.0000 k = 0.2857 0.4286 0.4286 ( 1609 PWs) bands (ev): -40.0354 -40.0354 -19.0883 -19.0883 -18.6052 -18.6052 -18.5640 -18.5640 5.4034 5.4034 7.2399 7.2399 7.2457 7.2457 10.2001 10.2001 11.0401 11.0401 11.4708 11.4708 11.6582 11.6582 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 0.0000 0.0000 k = 0.4286 0.4286 0.4286 ( 1610 PWs) bands (ev): -40.0310 -40.0310 -19.0987 -19.0987 -18.5982 -18.5982 -18.5974 -18.5974 6.2185 6.2185 7.0578 7.0578 7.0637 7.0637 9.7499 9.7499 11.5511 11.5511 11.7246 11.7246 11.7567 11.7567 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 0.0000 0.0000 the Fermi energy is 9.4090 ev ! total energy = -102.99206174 Ry Harris-Foulkes estimate = -102.99206174 Ry estimated scf accuracy < 4.4E-11 Ry The total energy is the sum of the following terms: one-electron contribution = -41.47534675 Ry hartree contribution = 24.46196406 Ry xc contribution = -20.10064551 Ry ewald contribution = -65.87795145 Ry smearing contrib. (-TS) = -0.00008208 Ry convergence has been achieved in 8 iterations Writing output data file ScAl.save init_run : 0.30s CPU 0.35s WALL ( 1 calls) electrons : 5.66s CPU 5.88s WALL ( 1 calls) Called by init_run: wfcinit : 0.14s CPU 0.16s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 4.32s CPU 4.51s WALL ( 9 calls) sum_band : 0.96s CPU 0.97s WALL ( 9 calls) v_of_rho : 0.02s CPU 0.03s WALL ( 9 calls) v_h : 0.00s CPU 0.00s WALL ( 9 calls) v_xc : 0.02s CPU 0.02s WALL ( 9 calls) newd : 0.34s CPU 0.36s WALL ( 9 calls) mix_rho : 0.02s CPU 0.01s WALL ( 9 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.01s WALL ( 380 calls) cegterg : 4.12s CPU 4.26s WALL ( 180 calls) Called by sum_band: sum_band:bec : 0.37s CPU 0.35s WALL ( 180 calls) addusdens : 0.22s CPU 0.22s WALL ( 9 calls) Called by *egterg: h_psi : 1.97s CPU 2.10s WALL ( 937 calls) s_psi : 0.10s CPU 0.11s WALL ( 937 calls) g_psi : 0.00s CPU 0.00s WALL ( 737 calls) cdiaghg : 1.85s CPU 1.86s WALL ( 897 calls) cegterg:over : 0.08s CPU 0.10s WALL ( 737 calls) cegterg:upda : 0.14s CPU 0.09s WALL ( 737 calls) cegterg:last : 0.04s CPU 0.03s WALL ( 187 calls) cdiaghg:chol : 0.11s CPU 0.11s WALL ( 897 calls) cdiaghg:inve : 0.03s CPU 0.02s WALL ( 897 calls) cdiaghg:para : 0.13s CPU 0.13s WALL ( 1794 calls) Called by h_psi: h_psi:vloc : 1.75s CPU 1.84s WALL ( 937 calls) h_psi:vnl : 0.22s CPU 0.25s WALL ( 937 calls) add_vuspsi : 0.15s CPU 0.14s WALL ( 937 calls) General routines calbec : 0.10s CPU 0.14s WALL ( 1117 calls) fft : 0.06s CPU 0.05s WALL ( 273 calls) ffts : 0.00s CPU 0.00s WALL ( 72 calls) fftw : 1.86s CPU 1.98s WALL ( 54532 calls) interpolate : 0.02s CPU 0.02s WALL ( 72 calls) Parallel routines fft_scatter : 1.06s CPU 1.10s WALL ( 54877 calls) PWSCF : 7.95s CPU 9.23s WALL This run was terminated on: 19:39:54 24Dec2016 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=