Program PWSCF v.5.3.0 (svn rev. 11974) starts on 31Jan2017 at 14:32:25 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 Pb.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 6S 5D 5D 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 2134 2134 309 Max 50 50 15 2137 2137 314 Sum 3489 3489 965 153699 153699 22335 bravais-lattice index = 14 lattice parameter (alat) = 11.6785 a.u. unit-cell volume = 1592.8025 (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.678507 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 Pb read from file: /users/gautes/Pseudo/Pb.rel-pbe-dn-rrkjus_psl.1.0.0.UPF MD5 check sum: 75f9fe62f948a5acf869e43db2c29526 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 1281 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 Pb 14.00 207.20000 Pb( 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: 153699 G-vectors FFT dimensions: ( 72, 72, 72) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.44 Mb ( 552, 52) NL pseudopotentials 0.83 Mb ( 276, 198) Each V/rho on FFT grid 0.08 Mb ( 5184) Each G-vector array 0.02 Mb ( 2136) G-vector shells 0.00 Mb ( 589) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 1.75 Mb ( 552, 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.96054, renormalised to 44.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 5.2 secs per-process dynamical memory: 100.4 Mb Self-consistent Calculation iteration # 1 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 4.3 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 3.85E-04, avg # of iterations = 5.2 total cpu time spent up to now is 10.6 secs total energy = -274.90703625 Ry Harris-Foulkes estimate = -275.03141140 Ry estimated scf accuracy < 0.17290802 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.93E-04, avg # of iterations = 3.8 total cpu time spent up to now is 13.8 secs total energy = -274.95341277 Ry Harris-Foulkes estimate = -274.98111902 Ry estimated scf accuracy < 0.04453602 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.01E-04, avg # of iterations = 3.1 total cpu time spent up to now is 16.8 secs total energy = -274.96549813 Ry Harris-Foulkes estimate = -274.96589121 Ry estimated scf accuracy < 0.00116793 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.65E-06, avg # of iterations = 9.0 total cpu time spent up to now is 21.2 secs total energy = -274.96585792 Ry Harris-Foulkes estimate = -274.96582996 Ry estimated scf accuracy < 0.00005744 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.31E-07, avg # of iterations = 3.1 total cpu time spent up to now is 24.1 secs total energy = -274.96587147 Ry Harris-Foulkes estimate = -274.96586879 Ry estimated scf accuracy < 0.00000213 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.83E-09, avg # of iterations = 4.7 total cpu time spent up to now is 27.7 secs total energy = -274.96587293 Ry Harris-Foulkes estimate = -274.96587282 Ry estimated scf accuracy < 0.00000005 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.06E-10, avg # of iterations = 3.7 total cpu time spent up to now is 31.0 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 19309 PWs) bands (ev): -16.3919 -16.3919 -14.1163 -14.1163 -14.1163 -14.1163 -11.6041 -11.6041 -11.6041 -11.6041 -11.5263 -11.5263 -8.7686 -8.7686 -7.8046 -7.8046 -7.8046 -7.8046 -5.1289 -5.1289 -3.6728 -3.6728 -3.6728 -3.6728 -2.8245 -2.8245 -0.2749 -0.2749 -0.1399 -0.1399 -0.1399 -0.1399 1.7725 1.7725 2.1538 2.1538 2.1538 2.1538 2.5906 2.5906 2.7995 2.7995 2.7995 2.7995 6.9021 6.9021 7.6355 7.6355 7.6355 7.6356 8.0868 8.0910 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 ( 19277 PWs) bands (ev): -16.3916 -16.3916 -14.1164 -14.1164 -14.1143 -14.1143 -11.6041 -11.6041 -11.5941 -11.5941 -11.5267 -11.5267 -8.6749 -8.6749 -7.9416 -7.9416 -7.8029 -7.8029 -5.1238 -5.1238 -3.6665 -3.6665 -3.6631 -3.6631 -3.2453 -3.2453 -0.1661 -0.1661 -0.1075 -0.1075 0.7349 0.7349 1.6972 1.6972 2.1084 2.1084 2.1144 2.1144 2.4873 2.4873 2.6927 2.6927 2.7308 2.7308 7.4698 7.4698 7.5801 7.5801 8.0107 8.0108 8.1264 8.1267 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 ( 19242 PWs) bands (ev): -16.3912 -16.3912 -14.1164 -14.1164 -14.1123 -14.1123 -11.6042 -11.6042 -11.5841 -11.5841 -11.5270 -11.5270 -8.5199 -8.5199 -8.1411 -8.1411 -7.8013 -7.8013 -5.1196 -5.1196 -3.6715 -3.6715 -3.6600 -3.6600 -3.5144 -3.5144 -0.0904 -0.0904 -0.0740 -0.0740 1.3561 1.3561 1.6089 1.6089 2.0912 2.0912 2.2462 2.2462 2.3640 2.3640 2.6015 2.6015 3.0768 3.0768 5.8932 5.8932 8.4096 8.4096 8.6062 8.6062 8.6941 8.6942 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 ( 19280 PWs) bands (ev): -16.3912 -16.3912 -14.1154 -14.1154 -14.1133 -14.1133 -11.5993 -11.5993 -11.5889 -11.5889 -11.5270 -11.5270 -8.5833 -8.5833 -7.9801 -7.9801 -7.8824 -7.8824 -5.1223 -5.1223 -3.6782 -3.6782 -3.6654 -3.6654 -3.6137 -3.6137 -0.0755 -0.0755 0.5951 0.5951 0.9038 0.9038 1.6356 1.6356 1.9428 1.9428 2.1279 2.1279 2.4769 2.4769 2.6509 2.6509 2.8196 2.8196 7.0739 7.0739 7.9033 7.9033 8.1038 8.1038 8.3524 8.3526 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 ( 19254 PWs) bands (ev): -16.3909 -16.3909 -14.1152 -14.1152 -14.1115 -14.1115 -11.5980 -11.5980 -11.5799 -11.5799 -11.5274 -11.5274 -8.4242 -8.4242 -8.1420 -8.1420 -7.9077 -7.9077 -5.1213 -5.1213 -3.9216 -3.9216 -3.6618 -3.6618 -3.6522 -3.6522 -0.0080 -0.0080 0.7787 0.7787 1.2181 1.2181 1.6947 1.6947 1.9847 1.9847 2.1474 2.1474 2.4866 2.4866 2.7051 2.7051 3.1902 3.1902 5.6930 5.6930 7.6296 7.6296 8.7154 8.7155 8.8096 8.8098 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 ( 19152 PWs) bands (ev): -16.3906 -16.3906 -14.1145 -14.1145 -14.1102 -14.1102 -11.5947 -11.5947 -11.5729 -11.5729 -11.5277 -11.5277 -8.2020 -8.2020 -8.1508 -8.1508 -8.1348 -8.1348 -5.1223 -5.1223 -4.2040 -4.2040 -3.6720 -3.6720 -3.6577 -3.6577 0.0754 0.0754 0.9809 0.9809 1.6220 1.6220 1.7828 1.7828 1.9540 1.9540 2.1917 2.1917 2.4859 2.4859 2.9025 2.9025 3.6772 3.6772 5.0458 5.0458 6.2376 6.2376 8.8636 8.8638 8.8791 8.8792 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 ( 19261 PWs) bands (ev): -16.3909 -16.3909 -14.1133 -14.1133 -14.1133 -14.1133 -11.5890 -11.5890 -11.5889 -11.5889 -11.5274 -11.5274 -8.4986 -8.4986 -7.9814 -7.9814 -7.9719 -7.9719 -5.1251 -5.1251 -4.0075 -4.0075 -3.6764 -3.6764 -3.6554 -3.6554 0.6346 0.6346 0.6797 0.6797 1.0198 1.0198 1.7102 1.7102 1.8118 1.8118 1.9806 1.9806 2.6332 2.6332 2.7195 2.7195 2.8643 2.8643 6.7066 6.7066 7.7427 7.7428 7.8771 7.8771 8.1756 8.1756 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 ( 19238 PWs) bands (ev): -16.3906 -16.3906 -14.1134 -14.1134 -14.1113 -14.1113 -11.5890 -11.5890 -11.5785 -11.5785 -11.5277 -11.5277 -8.3464 -8.3464 -8.1402 -8.1402 -7.9731 -7.9731 -5.1279 -5.1279 -4.2925 -4.2925 -3.6775 -3.6775 -3.6658 -3.6658 0.7155 0.7155 0.9413 0.9413 1.2285 1.2285 1.6330 1.6330 1.9132 1.9132 2.1841 2.1841 2.7460 2.7460 2.7984 2.7984 3.2466 3.2466 5.4757 5.4757 7.2235 7.2235 7.6700 7.6700 8.3110 8.3110 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 ( 19244 PWs) bands (ev): -16.3903 -16.3903 -14.1124 -14.1124 -14.1103 -14.1103 -11.5839 -11.5839 -11.5730 -11.5730 -11.5281 -11.5281 -8.1721 -8.1721 -8.1438 -8.1438 -8.1299 -8.1299 -5.1339 -5.1339 -4.5909 -4.5909 -3.6882 -3.6882 -3.6812 -3.6812 0.8890 0.8890 1.0422 1.0422 1.4885 1.4885 1.6895 1.6895 2.1745 2.1745 2.3362 2.3362 2.8181 2.8181 2.9731 2.9731 3.7309 3.7309 4.8305 4.8305 5.9906 5.9906 7.2914 7.2914 8.1959 8.1959 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 ( 19208 PWs) bands (ev): -16.3900 -16.3900 -14.1103 -14.1103 -14.1103 -14.1103 -11.5731 -11.5731 -11.5731 -11.5731 -11.5285 -11.5285 -8.1506 -8.1506 -8.1250 -8.1250 -8.1250 -8.1250 -5.1455 -5.1455 -4.9119 -4.9119 -3.7044 -3.7044 -3.7044 -3.7044 1.0868 1.0868 1.0868 1.0868 1.5934 1.5934 2.3644 2.3644 2.3644 2.3644 2.3960 2.3960 3.1557 3.1557 3.1557 3.1557 4.1214 4.1214 4.2383 4.2383 5.6933 5.6933 5.6933 5.6933 7.7188 7.7188 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.9866 0.9866 0.0134 0.0134 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 4.1798 ev ! total energy = -274.96587297 Ry Harris-Foulkes estimate = -274.96587297 Ry estimated scf accuracy < 1.5E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -96.42668251 Ry hartree contribution = 66.14163590 Ry xc contribution = -86.01630068 Ry ewald contribution = -158.66452123 Ry smearing contrib. (-TS) = -0.00000445 Ry convergence has been achieved in 7 iterations Writing output data file CsPbI3.save init_run : 1.16s CPU 1.27s WALL ( 1 calls) electrons : 25.53s CPU 25.85s WALL ( 1 calls) Called by init_run: wfcinit : 0.67s CPU 0.69s WALL ( 1 calls) potinit : 0.02s CPU 0.03s WALL ( 1 calls) Called by electrons: c_bands : 21.15s CPU 21.39s WALL ( 8 calls) sum_band : 3.14s CPU 3.17s 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.04s CPU 0.04s WALL ( 8 calls) newd : 1.32s CPU 1.35s WALL ( 8 calls) mix_rho : 0.02s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.11s CPU 0.10s WALL ( 170 calls) cegterg : 19.63s CPU 19.83s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.36s CPU 0.36s WALL ( 80 calls) addusdens : 0.44s CPU 0.44s WALL ( 8 calls) Called by *egterg: h_psi : 12.76s CPU 12.89s WALL ( 459 calls) s_psi : 0.52s CPU 0.59s WALL ( 459 calls) g_psi : 0.04s CPU 0.03s WALL ( 369 calls) cdiaghg : 4.81s CPU 4.84s WALL ( 439 calls) cegterg:over : 0.74s CPU 0.74s WALL ( 369 calls) cegterg:upda : 0.67s CPU 0.68s WALL ( 369 calls) cegterg:last : 0.17s CPU 0.16s WALL ( 80 calls) cdiaghg:chol : 0.20s CPU 0.22s WALL ( 439 calls) cdiaghg:inve : 0.08s CPU 0.11s WALL ( 439 calls) cdiaghg:para : 0.34s CPU 0.27s WALL ( 878 calls) Called by h_psi: h_psi:vloc : 11.02s CPU 11.19s WALL ( 459 calls) h_psi:vnl : 1.69s CPU 1.65s WALL ( 459 calls) add_vuspsi : 0.92s CPU 0.89s WALL ( 459 calls) General routines calbec : 1.00s CPU 0.99s WALL ( 539 calls) fft : 0.05s CPU 0.07s WALL ( 154 calls) fftw : 12.49s CPU 12.75s WALL ( 57756 calls) Parallel routines fft_scatter : 7.86s CPU 7.95s WALL ( 57910 calls) PWSCF : 30.52s CPU 33.80s WALL This run was terminated on: 14:32:59 31Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=