Program PWSCF v.5.3.0 (svn rev. 11974) starts on 31Jan2017 at 14:32:24 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 72 processors R & G space division: proc/nbgrp/npool/nimage = 72 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 Sn.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 5S 4D 4D 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: 6* 6 procs) Parallelization info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Min 48 48 13 2162 2162 313 Max 49 49 14 2169 2169 320 Sum 3513 3513 973 155955 155955 22743 bravais-lattice index = 14 lattice parameter (alat) = 11.7271 a.u. unit-cell volume = 1612.7566 (a.u.)^3 number of atoms/cell = 5 number of atomic types = 3 number of electrons = 44.00 number of Kohn-Sham states= 52 kinetic-energy cutoff = 80.0000 Ry charge density cutoff = 320.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)= 11.727073 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 Sn read from file: /users/gautes/Pseudo/Sn.rel-pbe-dn-rrkjus_psl.1.0.0.UPF MD5 check sum: 9499e19b89e678f7f9f46ad7de03aa82 Pseudo is Ultrasoft + core correction, Zval = 14.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1243 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 PseudoPot. # 3 for Cs read from file: /users/gautes/Pseudo/Cs.rel-pbe-oncvpsp.UPF MD5 check sum: b26d75608c6aa7e713949ad9f7a49707 Pseudo is Norm-conserving, Zval = 9.0 Generated using ONCVPSP code by D. R. Hamann Using radial grid of 2414 points, 14 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 l(11) = 3 l(12) = 3 l(13) = 3 l(14) = 3 atomic species valence mass pseudopotential Sn 14.00 118.71000 Sn( 1.00) I 7.00 126.90450 I( 1.00) Cs 9.00 132.90550 Cs( 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= 10 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0156250 k( 2) = ( 0.0000000 0.0000000 0.2500000), wk = 0.0937500 k( 3) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0468750 k( 4) = ( 0.0000000 0.2500000 0.2500000), wk = 0.1875000 k( 5) = ( 0.0000000 0.2500000 -0.5000000), wk = 0.1875000 k( 6) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0468750 k( 7) = ( 0.2500000 0.2500000 0.2500000), wk = 0.1250000 k( 8) = ( 0.2500000 0.2500000 -0.5000000), wk = 0.1875000 k( 9) = ( 0.2500000 -0.5000000 -0.5000000), wk = 0.0937500 k( 10) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0156250 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0156250 k( 2) = ( 0.0000000 0.0000000 0.2500000), wk = 0.0937500 k( 3) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0468750 k( 4) = ( 0.0000000 0.2500000 0.2500000), wk = 0.1875000 k( 5) = ( 0.0000000 0.2500000 -0.5000000), wk = 0.1875000 k( 6) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0468750 k( 7) = ( 0.2500000 0.2500000 0.2500000), wk = 0.1250000 k( 8) = ( 0.2500000 0.2500000 -0.5000000), wk = 0.1875000 k( 9) = ( 0.2500000 -0.5000000 -0.5000000), wk = 0.0937500 k( 10) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0156250 Dense grid: 155955 G-vectors FFT dimensions: ( 72, 72, 72) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.45 Mb ( 564, 52) NL pseudopotentials 0.85 Mb ( 282, 198) Each V/rho on FFT grid 0.08 Mb ( 5184) Each G-vector array 0.02 Mb ( 2169) G-vector shells 0.00 Mb ( 589) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 1.79 Mb ( 564, 208) Each subspace H/S matrix 0.02 Mb ( 34, 34) Each matrix 0.31 Mb ( 198, 2, 52) Arrays for rho mixing 0.63 Mb ( 5184, 8) Initial potential from superposition of free atoms starting charge 43.96048, renormalised to 44.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 5.7 secs per-process dynamical memory: 100.5 Mb Self-consistent Calculation iteration # 1 ecut= 80.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 = 2.94E-04, avg # of iterations = 5.3 total cpu time spent up to now is 11.2 secs total energy = -302.82468963 Ry Harris-Foulkes estimate = -302.91536737 Ry estimated scf accuracy < 0.13132026 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.98E-04, avg # of iterations = 2.7 total cpu time spent up to now is 14.0 secs total energy = -302.85849994 Ry Harris-Foulkes estimate = -302.87717364 Ry estimated scf accuracy < 0.03076536 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.99E-05, avg # of iterations = 3.3 total cpu time spent up to now is 17.1 secs total energy = -302.86697388 Ry Harris-Foulkes estimate = -302.86732867 Ry estimated scf accuracy < 0.00114568 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.60E-06, avg # of iterations = 8.2 total cpu time spent up to now is 21.3 secs total energy = -302.86729501 Ry Harris-Foulkes estimate = -302.86726002 Ry estimated scf accuracy < 0.00005571 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.27E-07, avg # of iterations = 3.0 total cpu time spent up to now is 24.1 secs total energy = -302.86730750 Ry Harris-Foulkes estimate = -302.86730593 Ry estimated scf accuracy < 0.00000212 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.81E-09, avg # of iterations = 4.1 total cpu time spent up to now is 27.6 secs total energy = -302.86730886 Ry Harris-Foulkes estimate = -302.86730877 Ry estimated scf accuracy < 0.00000005 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.07E-10, avg # of iterations = 4.0 total cpu time spent up to now is 30.9 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 19549 PWs) bands (ev): -18.3504 -18.3504 -18.3504 -18.3504 -17.3177 -17.3177 -17.3177 -17.3177 -17.3051 -17.3051 -16.7141 -16.7141 -8.9235 -8.9235 -8.1990 -8.1990 -8.1990 -8.1990 -5.4487 -5.4487 -3.9912 -3.9912 -3.9912 -3.9912 -1.8708 -1.8708 -0.5045 -0.5045 -0.5045 -0.5045 -0.3194 -0.3194 1.4818 1.4818 1.8506 1.8506 1.8506 1.8506 2.2973 2.2973 2.4867 2.4867 2.4867 2.4867 6.5914 6.5914 7.3480 7.3481 7.3482 7.3482 7.6988 7.7042 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 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 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.2500 ( 19529 PWs) bands (ev): -18.3504 -18.3504 -18.3502 -18.3502 -17.3177 -17.3177 -17.3171 -17.3171 -17.3052 -17.3052 -16.7138 -16.7138 -8.8464 -8.8464 -8.2925 -8.2925 -8.1974 -8.1974 -5.4438 -5.4438 -3.9853 -3.9853 -3.9810 -3.9810 -2.4452 -2.4452 -0.4727 -0.4727 -0.3220 -0.3220 0.5594 0.5594 1.4080 1.4080 1.8110 1.8110 1.8133 1.8133 2.1851 2.1851 2.3792 2.3792 2.4358 2.4358 7.1538 7.1538 7.2348 7.2348 7.7130 7.7130 7.8244 7.8247 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 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 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 ( 19506 PWs) bands (ev): -18.3504 -18.3504 -18.3499 -18.3499 -17.3177 -17.3177 -17.3166 -17.3166 -17.3052 -17.3052 -16.7135 -16.7135 -8.7267 -8.7267 -8.4288 -8.4288 -8.1958 -8.1958 -5.4393 -5.4393 -3.9794 -3.9794 -3.9722 -3.9722 -2.7711 -2.7711 -0.4396 -0.4396 -0.2403 -0.2403 1.0784 1.0784 1.3114 1.3114 1.7893 1.7893 1.9661 1.9661 2.0433 2.0433 2.2964 2.2964 2.8722 2.8722 5.4362 5.4362 8.0595 8.0595 8.2937 8.2937 8.3943 8.3943 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 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 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.2500 0.2500 ( 19496 PWs) bands (ev): -18.3503 -18.3503 -18.3501 -18.3501 -17.3174 -17.3174 -17.3169 -17.3169 -17.3052 -17.3052 -16.7135 -16.7135 -8.7708 -8.7708 -8.3188 -8.3188 -8.2512 -8.2512 -5.4424 -5.4424 -3.9859 -3.9859 -3.9770 -3.9770 -2.8962 -2.8962 -0.3492 -0.3492 0.2687 0.2687 0.7680 0.7680 1.3799 1.3799 1.6312 1.6312 1.8260 1.8260 2.1664 2.1664 2.3775 2.3775 2.5179 2.5179 6.9219 6.9219 7.2152 7.2152 7.6135 7.6135 7.9771 7.9771 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 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 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.2500-0.5000 ( 19462 PWs) bands (ev): -18.3503 -18.3503 -18.3499 -18.3499 -17.3173 -17.3173 -17.3164 -17.3164 -17.3053 -17.3053 -16.7132 -16.7132 -8.6479 -8.6479 -8.4296 -8.4296 -8.2682 -8.2682 -5.4411 -5.4411 -3.9832 -3.9832 -3.9773 -3.9773 -3.1955 -3.1955 -0.2839 -0.2839 0.4827 0.4827 0.9632 0.9632 1.4055 1.4055 1.6644 1.6644 1.8452 1.8452 2.1523 2.1523 2.4077 2.4077 3.0422 3.0422 5.2903 5.2903 6.9774 6.9774 8.2983 8.2984 8.4346 8.4347 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 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 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 ( 19496 PWs) bands (ev): -18.3502 -18.3502 -18.3497 -18.3497 -17.3172 -17.3172 -17.3161 -17.3161 -17.3054 -17.3054 -16.7129 -16.7129 -8.4797 -8.4797 -8.4320 -8.4320 -8.4287 -8.4287 -5.4427 -5.4427 -3.9907 -3.9907 -3.9772 -3.9772 -3.4956 -3.4956 -0.1975 -0.1975 0.6535 0.6535 1.3783 1.3783 1.4899 1.4899 1.5988 1.5988 1.8811 1.8811 2.0443 2.0443 2.6135 2.6135 3.6653 3.6653 4.9243 4.9243 5.2997 5.2997 8.3051 8.3051 8.5477 8.5478 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 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 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.2500 0.2500 0.2500 ( 19468 PWs) bands (ev): -18.3501 -18.3501 -18.3501 -18.3501 -17.3169 -17.3169 -17.3169 -17.3169 -17.3053 -17.3053 -16.7132 -16.7132 -8.7012 -8.7012 -8.3170 -8.3170 -8.3144 -8.3144 -5.4449 -5.4449 -3.9949 -3.9949 -3.9780 -3.9780 -3.2935 -3.2935 0.2292 0.2292 0.3946 0.3946 0.8698 0.8698 1.4812 1.4812 1.4890 1.4890 1.7126 1.7126 2.3702 2.3702 2.4410 2.4410 2.5402 2.5402 6.6844 6.6844 6.9301 6.9302 7.1141 7.1141 7.7782 7.7783 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 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 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.2500 0.2500-0.5000 ( 19514 PWs) bands (ev): -18.3501 -18.3501 -18.3498 -18.3498 -17.3169 -17.3169 -17.3163 -17.3163 -17.3054 -17.3054 -16.7129 -16.7129 -8.5842 -8.5842 -8.4278 -8.4278 -8.3126 -8.3126 -5.4475 -5.4475 -3.9976 -3.9976 -3.9861 -3.9861 -3.5806 -3.5806 0.3415 0.3415 0.6338 0.6338 1.0101 1.0101 1.3281 1.3281 1.5990 1.5990 1.8642 1.8642 2.4286 2.4286 2.5032 2.5032 3.1322 3.1322 5.1152 5.1152 6.6843 6.6843 6.9459 6.9459 7.9975 7.9975 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 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 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.2500-0.5000-0.5000 ( 19496 PWs) bands (ev): -18.3500 -18.3500 -18.3497 -18.3497 -17.3166 -17.3166 -17.3161 -17.3161 -17.3054 -17.3054 -16.7127 -16.7127 -8.4515 -8.4515 -8.4278 -8.4278 -8.4241 -8.4241 -5.4536 -5.4536 -4.0059 -4.0059 -3.9993 -3.9993 -3.8731 -3.8731 0.5317 0.5317 0.7078 0.7078 1.2237 1.2237 1.3707 1.3707 1.8836 1.8836 1.8990 1.8990 2.4786 2.4786 2.6851 2.6851 3.7596 3.7596 4.7377 4.7377 5.0829 5.0829 6.7138 6.7138 7.8911 7.8911 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 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 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 ( 19400 PWs) bands (ev): -18.3497 -18.3497 -18.3497 -18.3497 -17.3161 -17.3161 -17.3161 -17.3161 -17.3055 -17.3055 -16.7124 -16.7124 -8.4261 -8.4261 -8.4194 -8.4194 -8.4194 -8.4194 -5.4645 -5.4645 -4.1696 -4.1696 -4.0210 -4.0210 -4.0210 -4.0210 0.7465 0.7465 0.7465 0.7465 1.2746 1.2746 1.8991 1.8991 1.8991 1.8991 2.1179 2.1179 2.8301 2.8301 2.8301 2.8301 4.3431 4.3431 4.4400 4.4400 4.8143 4.8143 4.8143 4.8143 7.4228 7.4228 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 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 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9723 0.9723 0.0277 0.0277 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 4.3916 ev ! total energy = -302.86730890 Ry Harris-Foulkes estimate = -302.86730889 Ry estimated scf accuracy < 2.1E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -121.85125596 Ry hartree contribution = 78.68187057 Ry xc contribution = -101.69048031 Ry ewald contribution = -158.00743529 Ry smearing contrib. (-TS) = -0.00000791 Ry convergence has been achieved in 7 iterations Writing output data file CsSnI3.save init_run : 1.15s CPU 1.29s WALL ( 1 calls) electrons : 24.91s CPU 25.22s WALL ( 1 calls) Called by init_run: wfcinit : 0.65s CPU 0.69s WALL ( 1 calls) potinit : 0.03s CPU 0.03s WALL ( 1 calls) Called by electrons: c_bands : 20.54s CPU 20.73s WALL ( 8 calls) sum_band : 3.12s CPU 3.18s WALL ( 8 calls) v_of_rho : 0.04s CPU 0.04s WALL ( 8 calls) v_h : 0.00s CPU 0.00s WALL ( 8 calls) v_xc : 0.03s CPU 0.04s WALL ( 8 calls) newd : 1.32s CPU 1.35s WALL ( 8 calls) mix_rho : 0.03s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.10s CPU 0.11s WALL ( 170 calls) cegterg : 19.00s CPU 19.15s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.37s CPU 0.36s WALL ( 80 calls) addusdens : 0.44s CPU 0.45s WALL ( 8 calls) Called by *egterg: h_psi : 12.35s CPU 12.47s WALL ( 427 calls) s_psi : 0.54s CPU 0.59s WALL ( 427 calls) g_psi : 0.02s CPU 0.03s WALL ( 337 calls) cdiaghg : 4.60s CPU 4.56s WALL ( 407 calls) cegterg:over : 0.75s CPU 0.75s WALL ( 337 calls) cegterg:upda : 0.60s CPU 0.65s WALL ( 337 calls) cegterg:last : 0.18s CPU 0.18s WALL ( 81 calls) cdiaghg:chol : 0.21s CPU 0.21s WALL ( 407 calls) cdiaghg:inve : 0.11s CPU 0.11s WALL ( 407 calls) cdiaghg:para : 0.28s CPU 0.26s WALL ( 814 calls) Called by h_psi: h_psi:vloc : 10.73s CPU 10.80s WALL ( 427 calls) h_psi:vnl : 1.57s CPU 1.63s WALL ( 427 calls) add_vuspsi : 0.90s CPU 0.88s WALL ( 427 calls) General routines calbec : 0.90s CPU 0.99s WALL ( 507 calls) fft : 0.07s CPU 0.07s WALL ( 154 calls) fftw : 12.24s CPU 12.34s WALL ( 57404 calls) Parallel routines fft_scatter : 7.52s CPU 7.59s WALL ( 57558 calls) PWSCF : 29.95s CPU 34.02s WALL This run was terminated on: 14:32:58 31Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=