Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 16:58:45 This program is part of the open-source Quantum ESPRESSO suite for quantum simulation of materials; please cite "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009); URL http://www.quantum-espresso.org", in publications or presentations arising from this work. More details at http://www.quantum-espresso.org/quote Parallel version (MPI), running on 36 processors R & G space division: proc/nbgrp/npool/nimage = 36 Waiting for input... Reading input from standard input Current dimensions of program PWSCF are: Max number of different atomic species (ntypx) = 10 Max number of k-points (npk) = 40000 Max angular momentum in pseudopotentials (lmaxx) = 3 file Pb.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 6S 5D 5D 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 79 79 22 3195 3195 474 Max 80 80 23 3200 3200 479 Sum 2869 2869 793 115121 115121 17071 bravais-lattice index = 14 lattice parameter (alat) = 10.5919 a.u. unit-cell volume = 1188.2925 (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)= 10.591914 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 Cl read from file: /users/gautes/Pseudo/Cl.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 689bc0dc0caaad1d61892f37b0005eca 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 1157 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 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. # 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 Cl 7.00 35.45300 Cl( 1.00) Pb 14.00 207.20000 Pb( 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: 115121 G-vectors FFT dimensions: ( 64, 64, 64) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.65 Mb ( 816, 52) NL pseudopotentials 1.23 Mb ( 408, 198) Each V/rho on FFT grid 0.12 Mb ( 8192) Each G-vector array 0.02 Mb ( 3197) G-vector shells 0.00 Mb ( 624) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 2.59 Mb ( 816, 208) Each subspace H/S matrix 0.04 Mb ( 52, 52) Each matrix 0.31 Mb ( 198, 2, 52) Arrays for rho mixing 1.00 Mb ( 8192, 8) Initial potential from superposition of free atoms starting charge 43.96076, renormalised to 44.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 3.5 secs per-process dynamical memory: 101.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 = 3.7 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 7.34E-04, avg # of iterations = 2.1 total cpu time spent up to now is 8.5 secs total energy = -278.68946002 Ry Harris-Foulkes estimate = -278.92423453 Ry estimated scf accuracy < 0.31919921 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.25E-04, avg # of iterations = 3.5 total cpu time spent up to now is 11.6 secs total energy = -278.77807240 Ry Harris-Foulkes estimate = -278.84783859 Ry estimated scf accuracy < 0.12371420 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.81E-04, avg # of iterations = 2.4 total cpu time spent up to now is 14.2 secs total energy = -278.80771145 Ry Harris-Foulkes estimate = -278.80784083 Ry estimated scf accuracy < 0.00152809 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 1 eigenvalues not converged ethr = 3.47E-06, avg # of iterations = 11.1 total cpu time spent up to now is 18.7 secs total energy = -278.80820491 Ry Harris-Foulkes estimate = -278.80816452 Ry estimated scf accuracy < 0.00005714 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.30E-07, avg # of iterations = 3.7 total cpu time spent up to now is 21.8 secs total energy = -278.80821988 Ry Harris-Foulkes estimate = -278.80821775 Ry estimated scf accuracy < 0.00000243 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.52E-09, avg # of iterations = 3.4 total cpu time spent up to now is 25.0 secs total energy = -278.80822107 Ry Harris-Foulkes estimate = -278.80822101 Ry estimated scf accuracy < 0.00000017 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.87E-10, avg # of iterations = 3.1 total cpu time spent up to now is 28.1 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 14363 PWs) bands (ev): -15.8565 -15.8565 -13.7839 -13.7839 -13.7839 -13.7839 -11.4716 -11.4716 -11.4716 -11.4716 -11.1623 -11.1623 -10.9303 -10.9303 -9.9421 -9.9421 -9.9421 -9.9421 -4.5959 -4.5959 -3.1525 -3.1525 -3.1525 -3.1525 -2.2971 -2.2971 -0.7111 -0.7111 -0.3223 -0.3223 -0.3223 -0.3223 1.5725 1.5725 1.5725 1.5725 1.6098 1.6098 1.6301 1.6301 1.7650 1.7650 1.7650 1.7650 7.6880 7.6880 8.4050 8.4050 8.4050 8.4050 9.0597 9.0597 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 ( 14330 PWs) bands (ev): -15.8547 -15.8547 -13.7840 -13.7840 -13.7725 -13.7725 -11.4715 -11.4715 -11.4175 -11.4175 -11.1630 -11.1630 -10.8486 -10.8486 -10.1334 -10.1334 -9.9412 -9.9412 -4.5881 -4.5881 -3.1441 -3.1441 -3.1349 -3.1349 -2.9076 -2.9076 -0.5656 -0.5656 -0.2957 -0.2957 0.6150 0.6150 1.5134 1.5134 1.5220 1.5220 1.5812 1.5812 1.5927 1.5927 1.7006 1.7006 1.8312 1.8312 8.0468 8.0468 8.3934 8.3934 8.8052 8.8052 8.8639 8.8639 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 ( 14454 PWs) bands (ev): -15.8529 -15.8529 -13.7841 -13.7841 -13.7611 -13.7611 -11.4714 -11.4714 -11.3674 -11.3674 -11.1638 -11.1638 -10.5752 -10.5752 -10.5125 -10.5125 -9.9403 -9.9403 -4.5805 -4.5805 -3.2824 -3.2824 -3.1171 -3.1171 -3.1170 -3.1170 -0.5177 -0.5177 -0.2691 -0.2691 1.0320 1.0320 1.4542 1.4542 1.4722 1.4722 1.5511 1.5511 1.5612 1.5612 1.6201 1.6201 2.5744 2.5744 6.8871 6.8871 8.8624 8.8624 9.3006 9.3006 9.4060 9.4060 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 ( 14359 PWs) bands (ev): -15.8531 -15.8531 -13.7784 -13.7784 -13.7665 -13.7665 -11.4478 -11.4478 -11.3739 -11.3739 -11.1637 -11.1637 -10.7873 -10.7873 -10.2039 -10.2039 -10.0444 -10.0444 -4.5883 -4.5883 -3.4069 -3.4069 -3.1465 -3.1465 -3.1194 -3.1194 -0.4010 -0.4010 0.3099 0.3099 0.7757 0.7757 1.3154 1.3154 1.4942 1.4942 1.5620 1.5620 1.6272 1.6272 1.8112 1.8112 2.0005 2.0005 7.6246 7.6246 8.6009 8.6009 8.9666 8.9667 9.1722 9.1722 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 ( 14374 PWs) bands (ev): -15.8515 -15.8515 -13.7769 -13.7769 -13.7565 -13.7565 -11.4422 -11.4422 -11.3102 -11.3102 -11.1644 -11.1644 -10.5591 -10.5591 -10.5141 -10.5141 -10.0775 -10.0775 -4.5886 -4.5886 -3.7303 -3.7303 -3.1480 -3.1480 -3.1194 -3.1194 -0.3396 -0.3396 0.5412 0.5412 0.9132 0.9132 1.3543 1.3543 1.4535 1.4535 1.4946 1.4946 1.5960 1.5960 1.8327 1.8327 2.8079 2.8079 6.5513 6.5513 8.3708 8.3708 9.3794 9.3794 9.5165 9.5165 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 ( 14332 PWs) bands (ev): -15.8501 -15.8501 -13.7730 -13.7730 -13.7488 -13.7488 -11.4302 -11.4302 -11.2161 -11.2161 -11.1651 -11.1651 -10.5232 -10.5232 -10.5092 -10.5092 -10.2617 -10.2617 -4.5967 -4.5967 -4.0476 -4.0476 -3.1739 -3.1739 -3.1199 -3.1199 -0.2578 -0.2578 0.8030 0.8030 1.3485 1.3485 1.3987 1.3987 1.4009 1.4009 1.4216 1.4216 1.4960 1.4960 2.0075 2.0075 3.4554 3.4554 5.9032 5.9032 7.2796 7.2796 9.2462 9.2462 9.4852 9.4853 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 ( 14360 PWs) bands (ev): -15.8517 -15.8517 -13.7666 -13.7666 -13.7666 -13.7666 -11.3743 -11.3743 -11.3732 -11.3732 -11.1645 -11.1645 -10.7401 -10.7401 -10.2054 -10.2054 -10.1987 -10.1987 -4.5961 -4.5961 -3.8308 -3.8308 -3.1662 -3.1662 -3.1285 -3.1285 0.2468 0.2468 0.4445 0.4445 0.8380 0.8380 1.2770 1.2770 1.2843 1.2843 1.7378 1.7378 1.8030 1.8030 1.8219 1.8219 2.1526 2.1526 7.1728 7.1728 8.4827 8.4827 8.5022 8.5022 9.1816 9.1816 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 ( 14332 PWs) bands (ev): -15.8503 -15.8503 -13.7667 -13.7667 -13.7549 -13.7549 -11.3736 -11.3736 -11.2915 -11.2915 -11.1652 -11.1652 -10.5483 -10.5483 -10.5134 -10.5134 -10.2002 -10.2002 -4.6038 -4.6038 -4.1373 -4.1373 -3.1740 -3.1740 -3.1480 -3.1480 0.3910 0.3910 0.7417 0.7417 0.9424 0.9424 1.1591 1.1591 1.3327 1.3327 1.7227 1.7227 1.7815 1.7815 1.8752 1.8752 2.9970 2.9970 6.2193 6.2193 7.8664 7.8664 8.3380 8.3380 9.4644 9.4644 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 ( 14360 PWs) bands (ev): -15.8490 -15.8490 -13.7610 -13.7610 -13.7489 -13.7489 -11.3387 -11.3387 -11.2162 -11.2162 -11.1659 -11.1659 -10.5207 -10.5207 -10.5063 -10.5063 -10.3700 -10.3700 -4.6188 -4.6188 -4.4410 -4.4410 -3.2005 -3.2005 -3.1746 -3.1746 0.6258 0.6258 0.9311 0.9311 1.1319 1.1319 1.1378 1.1378 1.4894 1.4894 1.7617 1.7617 1.8187 1.8187 2.0150 2.0150 3.6389 3.6389 5.5857 5.5857 6.9356 6.9356 7.8536 7.8536 9.8421 9.8421 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 ( 14400 PWs) bands (ev): -15.8479 -15.8479 -13.7490 -13.7490 -13.7490 -13.7490 -11.2164 -11.2164 -11.2164 -11.2164 -11.1667 -11.1667 -10.5249 -10.5249 -10.5033 -10.5033 -10.5033 -10.5033 -4.7442 -4.7442 -4.6409 -4.6409 -3.2267 -3.2267 -3.2267 -3.2267 1.0448 1.0448 1.0643 1.0643 1.0643 1.0643 1.4925 1.4925 1.4925 1.4925 1.9734 1.9734 2.0402 2.0402 2.0402 2.0402 4.2787 4.2787 4.9673 4.9673 6.6137 6.6137 6.6137 6.6137 9.4897 9.4898 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 the Fermi energy is 4.6382 ev ! total energy = -278.80822115 Ry Harris-Foulkes estimate = -278.80822114 Ry estimated scf accuracy < 3.7E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -103.54538687 Ry hartree contribution = 70.86273906 Ry xc contribution = -71.18413870 Ry ewald contribution = -174.94143465 Ry smearing contrib. (-TS) = -0.00000000 Ry convergence has been achieved in 7 iterations Writing output data file CsPbCl3.save init_run : 1.20s CPU 1.31s WALL ( 1 calls) electrons : 24.25s CPU 24.62s WALL ( 1 calls) Called by init_run: wfcinit : 0.65s CPU 0.68s WALL ( 1 calls) potinit : 0.02s CPU 0.03s WALL ( 1 calls) Called by electrons: c_bands : 19.51s CPU 19.75s WALL ( 8 calls) sum_band : 3.32s CPU 3.34s WALL ( 8 calls) v_of_rho : 0.05s CPU 0.05s 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.51s CPU 1.54s WALL ( 8 calls) mix_rho : 0.03s CPU 0.03s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.15s CPU 0.16s WALL ( 170 calls) cegterg : 17.24s CPU 17.38s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.35s CPU 0.36s WALL ( 80 calls) addusdens : 0.66s CPU 0.66s WALL ( 8 calls) Called by *egterg: h_psi : 12.08s CPU 12.19s WALL ( 420 calls) s_psi : 0.74s CPU 0.75s WALL ( 420 calls) g_psi : 0.05s CPU 0.04s WALL ( 330 calls) cdiaghg : 2.50s CPU 2.59s WALL ( 400 calls) cegterg:over : 0.76s CPU 0.73s WALL ( 330 calls) cegterg:upda : 0.83s CPU 0.80s WALL ( 330 calls) cegterg:last : 0.19s CPU 0.20s WALL ( 80 calls) cdiaghg:chol : 0.17s CPU 0.16s WALL ( 400 calls) cdiaghg:inve : 0.08s CPU 0.10s WALL ( 400 calls) cdiaghg:para : 0.15s CPU 0.16s WALL ( 800 calls) Called by h_psi: h_psi:vloc : 10.14s CPU 10.29s WALL ( 420 calls) h_psi:vnl : 1.88s CPU 1.84s WALL ( 420 calls) add_vuspsi : 1.07s CPU 1.01s WALL ( 420 calls) General routines calbec : 1.08s CPU 1.10s WALL ( 500 calls) fft : 0.07s CPU 0.06s WALL ( 154 calls) fftw : 11.28s CPU 11.46s WALL ( 55468 calls) Parallel routines fft_scatter : 3.88s CPU 3.78s WALL ( 55622 calls) PWSCF : 28.52s CPU 30.11s WALL This run was terminated on: 16:59:15 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=