Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 21:16: 0 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 Te.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 5S 5P 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 28 17 5 688 317 53 Max 29 18 6 691 336 56 Sum 1033 621 185 24837 11753 1935 bravais-lattice index = 14 lattice parameter (alat) = 7.0071 a.u. unit-cell volume = 344.0453 (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 = 40.0000 Ry charge density cutoff = 263.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.007104 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 Te read from file: /users/gautes/Pseudo/Te.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: d7fe42880ce6991a805b4693d1f2335c Pseudo is Ultrasoft + core correction, Zval = 6.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1245 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 Sr read from file: /users/gautes/Pseudo/Sr.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 971034ad0a3914f6282d12395b0824f8 Pseudo is Ultrasoft + core correction, Zval = 10.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1221 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 Te 6.00 127.60000 Te( 1.00) Sr 10.00 87.62000 Sr( 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: 24837 G-vectors FFT dimensions: ( 40, 40, 40) Smooth grid: 11753 G-vectors FFT dimensions: ( 30, 30, 30) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.03 Mb ( 82, 24) NL pseudopotentials 0.04 Mb ( 41, 68) Each V/rho on FFT grid 0.05 Mb ( 3200) Each G-vector array 0.01 Mb ( 691) G-vector shells 0.00 Mb ( 185) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.12 Mb ( 82, 96) Each subspace H/S matrix 0.01 Mb ( 24, 24) Each matrix 0.05 Mb ( 68, 2, 24) Arrays for rho mixing 0.39 Mb ( 3200, 8) Initial potential from superposition of free atoms starting charge 15.99218, renormalised to 16.00000 Starting wfc are 24 randomized atomic wfcs total cpu time spent up to now is 2.4 secs per-process dynamical memory: 24.5 Mb Self-consistent Calculation iteration # 1 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 5.9 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.17E-04, avg # of iterations = 4.5 total cpu time spent up to now is 4.0 secs total energy = -98.91094010 Ry Harris-Foulkes estimate = -98.91620430 Ry estimated scf accuracy < 0.01816591 Ry iteration # 2 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.14E-04, avg # of iterations = 3.4 total cpu time spent up to now is 4.6 secs total energy = -98.91393430 Ry Harris-Foulkes estimate = -98.91557091 Ry estimated scf accuracy < 0.00306815 Ry iteration # 3 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.92E-05, avg # of iterations = 2.7 total cpu time spent up to now is 5.1 secs total energy = -98.91492405 Ry Harris-Foulkes estimate = -98.91489413 Ry estimated scf accuracy < 0.00037593 Ry iteration # 4 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.35E-06, avg # of iterations = 4.3 total cpu time spent up to now is 5.8 secs total energy = -98.91497015 Ry Harris-Foulkes estimate = -98.91497081 Ry estimated scf accuracy < 0.00000392 Ry iteration # 5 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.45E-08, avg # of iterations = 3.8 total cpu time spent up to now is 6.5 secs total energy = -98.91497268 Ry Harris-Foulkes estimate = -98.91497286 Ry estimated scf accuracy < 0.00000041 Ry iteration # 6 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.56E-09, avg # of iterations = 2.4 total cpu time spent up to now is 7.0 secs total energy = -98.91497275 Ry Harris-Foulkes estimate = -98.91497275 Ry estimated scf accuracy < 0.00000001 Ry iteration # 7 ecut= 40.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.56E-11, avg # of iterations = 3.6 total cpu time spent up to now is 7.6 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1419 PWs) bands (ev): -24.4518 -24.4518 -7.1870 -7.1870 -6.0436 -6.0436 -6.0436 -6.0436 -1.3097 -1.3097 9.7489 9.7489 10.7120 10.7120 10.7120 10.7120 11.1061 11.1061 11.1061 11.1061 14.8570 14.8570 16.2511 16.2511 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.3351 0.3351 0.3351 0.3351 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.1429 ( 1450 PWs) bands (ev): -24.4450 -24.4450 -7.2212 -7.2212 -6.1335 -6.1335 -6.0306 -6.0306 -1.0667 -1.0667 8.6695 8.6695 10.0341 10.0341 10.5643 10.5643 11.2744 11.2744 12.3127 12.3127 15.8366 15.8366 16.3571 16.3574 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.2857 ( 1454 PWs) bands (ev): -24.4299 -24.4299 -7.3207 -7.3207 -6.2868 -6.2868 -6.0010 -6.0010 -0.4782 -0.4782 6.8456 6.8456 9.7218 9.7218 10.2426 10.2426 11.7391 11.7391 13.9700 13.9700 16.5554 16.5554 16.7749 16.7749 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.4286 ( 1458 PWs) bands (ev): -24.4177 -24.4177 -7.4214 -7.4214 -6.3686 -6.3686 -5.9771 -5.9771 0.0835 0.0835 5.6134 5.6134 9.5099 9.5099 9.9955 9.9955 12.2817 12.2817 15.8373 15.8373 16.0608 16.0608 16.7410 16.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.0000 0.1429 0.1429 ( 1440 PWs) bands (ev): -24.4384 -24.4384 -7.2482 -7.2482 -6.1728 -6.1728 -6.0615 -6.0615 -0.8603 -0.8603 8.4827 8.4827 9.1007 9.1007 10.2721 10.2721 11.8689 11.8689 13.0512 13.0512 15.3558 15.3558 15.6727 15.6727 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.1429 0.2857 ( 1465 PWs) bands (ev): -24.4234 -24.4234 -7.3298 -7.3298 -6.3049 -6.3049 -6.0447 -6.0447 -0.3651 -0.3651 6.9828 6.9828 8.9917 8.9917 10.1255 10.1255 12.2664 12.2664 14.4197 14.4197 15.4174 15.4174 15.4740 15.4740 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.1429 0.4286 ( 1464 PWs) bands (ev): -24.4113 -24.4113 -7.4160 -7.4160 -6.3846 -6.3846 -6.0197 -6.0197 0.0940 0.0940 5.8762 5.8762 9.0301 9.0301 9.9976 9.9976 12.4707 12.4707 15.6902 15.6902 15.7277 15.7277 16.4263 16.4263 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.2857 0.2857 ( 1474 PWs) bands (ev): -24.4088 -24.4088 -7.3685 -7.3685 -6.3780 -6.3780 -6.0780 -6.0780 -0.1068 -0.1068 7.0627 7.0627 7.9766 7.9766 10.3119 10.3119 12.5942 12.5942 14.6451 14.6451 14.6918 14.6918 15.5255 15.5255 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.2857 0.4286 ( 1477 PWs) bands (ev): -24.3970 -24.3970 -7.4198 -7.4198 -6.4395 -6.4395 -6.0689 -6.0689 0.1109 0.1109 6.5355 6.5355 8.0249 8.0249 10.4715 10.4715 12.2584 12.2584 14.6192 14.6192 14.9987 14.9987 15.8360 15.8360 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.4286 0.4286 ( 1480 PWs) bands (ev): -24.3855 -24.3855 -7.4420 -7.4420 -6.4874 -6.4874 -6.0740 -6.0740 0.1177 0.1177 7.1050 7.1050 7.5110 7.5110 11.0583 11.0583 11.4608 11.4608 14.0923 14.0923 14.5224 14.5224 14.6165 14.6165 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.9442 0.9442 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.1429 ( 1460 PWs) bands (ev): -24.4318 -24.4318 -7.2700 -7.2700 -6.1819 -6.1819 -6.1166 -6.1166 -0.6833 -0.6833 8.4509 8.4509 9.0442 9.0442 9.2090 9.2090 12.8283 12.8283 12.9543 12.9543 15.9461 15.9481 16.0314 16.0316 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.1429 0.1429 0.2857 ( 1471 PWs) bands (ev): -24.4170 -24.4170 -7.3382 -7.3382 -6.3056 -6.3056 -6.1005 -6.1005 -0.2602 -0.2602 7.1440 7.1440 8.8067 8.8067 9.4035 9.4035 12.7178 12.7178 14.2405 14.2405 15.2426 15.2426 16.6196 16.6196 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.1429 0.1429 0.4286 ( 1475 PWs) bands (ev): -24.4051 -24.4051 -7.4126 -7.4126 -6.3875 -6.3875 -6.0704 -6.0704 0.1260 0.1260 6.1003 6.1003 8.7093 8.7093 9.7424 9.7424 12.7186 12.7186 14.5048 14.5048 15.4891 15.4891 17.0340 17.0340 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.1429 0.2857 0.2857 ( 1459 PWs) bands (ev): -24.4026 -24.4026 -7.3730 -7.3730 -6.3566 -6.3566 -6.1477 -6.1477 -0.0153 -0.0153 7.2480 7.2480 7.8822 7.8822 9.3429 9.3429 12.9503 12.9503 14.1656 14.1656 15.3585 15.3585 16.1241 16.1241 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.1429 0.2857 0.4286 ( 1478 PWs) bands (ev): -24.3910 -24.3910 -7.4189 -7.4189 -6.4201 -6.4201 -6.1336 -6.1336 0.1968 0.1968 6.6222 6.6222 7.8426 7.8426 9.6957 9.6957 12.6016 12.6016 14.0690 14.0690 14.6157 14.6157 16.8568 16.8568 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.1429 0.4286 0.4286 ( 1469 PWs) bands (ev): -24.3796 -24.3796 -7.4408 -7.4408 -6.4638 -6.4638 -6.1399 -6.1399 0.2527 0.2527 6.9665 6.9665 7.3502 7.3502 9.8014 9.8014 11.8851 11.8851 14.1981 14.1981 14.2416 14.2416 16.3866 16.3866 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.2857 0.2857 0.2857 ( 1469 PWs) bands (ev): -24.3887 -24.3887 -7.3952 -7.3952 -6.3283 -6.3283 -6.2625 -6.2625 0.2069 0.2069 7.2907 7.2907 7.7738 7.7738 8.2237 8.2237 13.7439 13.7439 13.8120 13.8120 14.1597 14.1597 17.4301 17.4302 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.2857 0.2857 0.4286 ( 1476 PWs) bands (ev): -24.3774 -24.3774 -7.4272 -7.4272 -6.3893 -6.3893 -6.2477 -6.2477 0.4082 0.4082 6.7242 6.7242 7.3928 7.3928 8.5572 8.5572 13.2436 13.2436 13.6799 13.6799 14.3843 14.3843 17.9045 17.9045 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.2857 0.4286 0.4286 ( 1477 PWs) bands (ev): -24.3665 -24.3665 -7.4462 -7.4462 -6.4101 -6.4101 -6.2733 -6.2733 0.6014 0.6014 6.6142 6.6142 7.0135 7.0135 8.0289 8.0289 13.0610 13.0610 14.3877 14.3877 14.5286 14.5286 18.2850 18.2851 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.4286 0.4286 0.4286 ( 1478 PWs) bands (ev): -24.3559 -24.3559 -7.4604 -7.4604 -6.3697 -6.3697 -6.3607 -6.3607 0.9482 0.9482 6.2750 6.2750 6.7358 6.7358 6.9733 6.9733 14.3370 14.3370 14.7421 14.7421 14.8818 14.8818 18.7642 18.7643 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 11.0968 ev ! total energy = -98.91497275 Ry Harris-Foulkes estimate = -98.91497276 Ry estimated scf accuracy < 1.5E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -4.92244360 Ry hartree contribution = 10.18056562 Ry xc contribution = -35.37078592 Ry ewald contribution = -68.80228635 Ry smearing contrib. (-TS) = -0.00002251 Ry convergence has been achieved in 7 iterations Writing output data file SrTe.save init_run : 0.26s CPU 0.31s WALL ( 1 calls) electrons : 5.02s CPU 5.22s 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.96s CPU 4.09s WALL ( 8 calls) sum_band : 0.80s CPU 0.82s WALL ( 8 calls) v_of_rho : 0.00s CPU 0.02s WALL ( 8 calls) v_h : 0.00s CPU 0.00s WALL ( 8 calls) v_xc : 0.00s CPU 0.02s WALL ( 8 calls) newd : 0.25s CPU 0.27s WALL ( 8 calls) mix_rho : 0.02s CPU 0.01s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.02s CPU 0.01s WALL ( 340 calls) cegterg : 3.86s CPU 3.92s WALL ( 160 calls) Called by sum_band: sum_band:bec : 0.29s CPU 0.32s WALL ( 160 calls) addusdens : 0.14s CPU 0.14s WALL ( 8 calls) Called by *egterg: h_psi : 1.89s CPU 1.93s WALL ( 790 calls) s_psi : 0.10s CPU 0.10s WALL ( 790 calls) g_psi : 0.00s CPU 0.00s WALL ( 610 calls) cdiaghg : 1.75s CPU 1.72s WALL ( 750 calls) cegterg:over : 0.06s CPU 0.09s WALL ( 610 calls) cegterg:upda : 0.07s CPU 0.07s WALL ( 610 calls) cegterg:last : 0.02s CPU 0.03s WALL ( 162 calls) cdiaghg:chol : 0.05s CPU 0.09s WALL ( 750 calls) cdiaghg:inve : 0.01s CPU 0.02s WALL ( 750 calls) cdiaghg:para : 0.11s CPU 0.12s WALL ( 1500 calls) Called by h_psi: h_psi:vloc : 1.63s CPU 1.70s WALL ( 790 calls) h_psi:vnl : 0.26s CPU 0.23s WALL ( 790 calls) add_vuspsi : 0.15s CPU 0.13s WALL ( 790 calls) General routines calbec : 0.13s CPU 0.13s WALL ( 950 calls) fft : 0.06s CPU 0.04s WALL ( 242 calls) ffts : 0.00s CPU 0.00s WALL ( 64 calls) fftw : 1.78s CPU 1.86s WALL ( 51568 calls) interpolate : 0.02s CPU 0.01s WALL ( 64 calls) Parallel routines fft_scatter : 0.94s CPU 1.04s WALL ( 51874 calls) PWSCF : 7.36s CPU 8.97s WALL This run was terminated on: 21:16: 9 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=