Program PWSCF v.5.4.0 starts on 20Mar2017 at 21:34:57 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 & OpenMP), running on 144 processor cores Number of MPI processes: 72 Threads/MPI process: 2 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 Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S renormalized Subspace diagonalization in iterative solution of the eigenvalue problem: one sub-group per band group 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 39 39 11 1596 1596 235 Max 41 41 12 1600 1600 240 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 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 Br read from file: /users/gautes/Pseudo/Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 5d639767687de88634f0d235e6526498 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 1213 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 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) Br 7.00 79.90400 Br( 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.33 Mb ( 410, 52) NL pseudopotentials 0.43 Mb ( 205, 138) Each V/rho on FFT grid 0.06 Mb ( 4096) Each G-vector array 0.01 Mb ( 1600) G-vector shells 0.00 Mb ( 470) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 1.30 Mb ( 410, 208) Each subspace H/S matrix 0.02 Mb ( 34, 34) Each matrix 0.22 Mb ( 138, 2, 52) Arrays for rho mixing 0.50 Mb ( 4096, 8) Initial potential from superposition of free atoms starting charge 43.96071, renormalised to 44.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 3.8 secs per-process dynamical memory: 8.9 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.6 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 4.39E-04, avg # of iterations = 3.4 total cpu time spent up to now is 9.8 secs total energy = -264.28894775 Ry Harris-Foulkes estimate = -264.42741777 Ry estimated scf accuracy < 0.19351380 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.40E-04, avg # of iterations = 5.1 total cpu time spent up to now is 13.5 secs total energy = -264.34252586 Ry Harris-Foulkes estimate = -264.37052589 Ry estimated scf accuracy < 0.04567846 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.04E-04, avg # of iterations = 5.6 total cpu time spent up to now is 17.2 secs total energy = -264.35427630 Ry Harris-Foulkes estimate = -264.35439099 Ry estimated scf accuracy < 0.00082857 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.88E-06, avg # of iterations = 7.2 total cpu time spent up to now is 21.8 secs total energy = -264.35453697 Ry Harris-Foulkes estimate = -264.35450492 Ry estimated scf accuracy < 0.00004155 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.44E-08, avg # of iterations = 3.8 total cpu time spent up to now is 25.3 secs total energy = -264.35454820 Ry Harris-Foulkes estimate = -264.35454636 Ry estimated scf accuracy < 0.00000118 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.69E-09, avg # of iterations = 3.9 total cpu time spent up to now is 29.5 secs total energy = -264.35454902 Ry Harris-Foulkes estimate = -264.35454895 Ry estimated scf accuracy < 0.00000004 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-10, avg # of iterations = 3.6 total cpu time spent up to now is 33.0 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 14363 PWs) bands (ev): -14.8528 -14.8528 -12.6489 -12.6489 -12.6489 -12.6489 -10.4195 -10.4195 -10.4195 -10.4195 -10.0783 -10.0783 -10.0047 -10.0047 -8.7566 -8.7566 -8.7566 -8.7566 -3.6217 -3.6217 -2.1994 -2.1994 -2.1994 -2.1994 -1.2101 -1.2101 0.0898 0.0898 0.4457 0.4457 0.4457 0.4457 2.6998 2.6998 2.7658 2.7658 2.7658 2.7658 2.9922 2.9922 3.2443 3.2443 3.2443 3.2443 8.4775 8.4775 9.0449 9.0449 9.0449 9.0449 9.8520 9.8520 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): -14.8505 -14.8505 -12.6489 -12.6489 -12.6324 -12.6324 -10.4193 -10.4193 -10.3669 -10.3669 -10.0058 -10.0058 -9.9660 -9.9660 -9.0077 -9.0077 -8.7553 -8.7553 -3.6083 -3.6083 -2.1743 -2.1743 -2.1732 -2.1732 -1.9713 -1.9713 0.2453 0.2453 0.4680 0.4680 1.5765 1.5765 2.6128 2.6128 2.6982 2.6982 2.7485 2.7485 2.9094 2.9094 3.1338 3.1338 3.2981 3.2981 9.1666 9.1666 9.5282 9.5282 9.5899 9.5900 9.8586 9.8587 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): -14.8482 -14.8482 -12.6490 -12.6490 -12.6162 -12.6162 -10.4191 -10.4191 -10.3329 -10.3329 -10.0068 -10.0068 -9.5755 -9.5755 -9.5182 -9.5182 -8.7539 -8.7539 -3.5960 -3.5960 -2.3721 -2.3721 -2.1486 -2.1486 -2.1427 -2.1427 0.3052 0.3052 0.4892 0.4892 2.0193 2.0193 2.5066 2.5066 2.6361 2.6361 2.7789 2.7789 2.8570 2.8570 2.9995 2.9995 4.2228 4.2228 8.0517 8.0517 9.9358 9.9358 10.2974 10.2974 10.4697 10.4697 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): -14.8485 -14.8485 -12.6409 -12.6409 -12.6237 -12.6237 -10.3979 -10.3979 -10.3073 -10.3073 -10.0068 -10.0068 -9.8891 -9.8891 -9.0985 -9.0985 -8.8919 -8.8919 -3.6063 -3.6063 -2.5648 -2.5648 -2.1819 -2.1819 -2.1502 -2.1502 0.4092 0.4092 1.2568 1.2568 1.7186 1.7186 2.4155 2.4155 2.6576 2.6576 2.7319 2.7319 2.9546 2.9546 3.2220 3.2220 3.4543 3.4543 9.2036 9.2037 9.6400 9.6400 9.8042 9.8042 10.0172 10.0172 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): -14.8465 -14.8465 -12.6389 -12.6389 -12.6093 -12.6093 -10.3937 -10.3937 -10.2421 -10.2421 -10.0078 -10.0078 -9.5631 -9.5631 -9.5203 -9.5203 -8.9357 -8.9357 -3.6050 -3.6050 -2.9288 -2.9288 -2.1762 -2.1762 -2.1473 -2.1473 0.4823 0.4823 1.5283 1.5283 1.8512 1.8512 2.4612 2.4612 2.6059 2.6059 2.6835 2.6835 2.9594 2.9594 3.2081 3.2081 4.4473 4.4473 7.7068 7.7068 9.7680 9.7680 10.4522 10.4528 10.5662 10.5675 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): -14.8448 -14.8448 -12.6335 -12.6335 -12.5978 -12.5978 -10.3864 -10.3864 -10.0942 -10.0942 -10.0088 -10.0088 -9.5327 -9.5327 -9.5137 -9.5137 -9.1850 -9.1850 -3.6133 -3.6133 -3.2938 -3.2938 -2.2073 -2.2073 -2.1470 -2.1470 0.5874 0.5874 1.6842 1.6842 2.3885 2.3885 2.5062 2.5062 2.5363 2.5363 2.6736 2.6736 2.9894 2.9894 3.4644 3.4644 5.2376 5.2376 6.9296 6.9296 8.2453 8.2453 10.4130 10.4130 10.7501 10.8320 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): -14.8468 -14.8468 -12.6238 -12.6238 -12.6238 -12.6238 -10.3080 -10.3080 -10.3059 -10.3059 -10.0079 -10.0079 -9.8296 -9.8296 -9.1001 -9.1001 -9.0914 -9.0914 -3.6154 -3.6154 -3.0809 -3.0809 -2.2077 -2.2077 -2.1554 -2.1554 1.2524 1.2524 1.3886 1.3886 1.7667 1.7667 2.3760 2.3760 2.4193 2.4193 2.9078 2.9078 3.1836 3.1836 3.3028 3.3028 3.5416 3.5416 8.6938 8.6938 9.5425 9.5425 9.7103 9.7103 10.0660 10.0660 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): -14.8450 -14.8450 -12.6239 -12.6239 -12.6069 -12.6069 -10.3066 -10.3066 -10.2149 -10.2149 -10.0089 -10.0089 -9.5537 -9.5537 -9.5189 -9.5189 -9.0930 -9.0930 -3.6254 -3.6254 -3.4318 -3.4318 -2.2145 -2.2145 -2.1810 -2.1810 1.4046 1.4046 1.7521 1.7521 1.8829 1.8829 2.2358 2.2358 2.5330 2.5330 2.9354 2.9354 3.1966 3.1966 3.2868 3.2868 4.5989 4.5989 7.3677 7.3677 9.1495 9.1495 9.5369 9.5369 10.2134 10.2134 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): -14.8435 -14.8435 -12.6157 -12.6157 -12.5980 -12.5980 -10.2697 -10.2697 -10.0944 -10.0944 -10.0100 -10.0100 -9.5288 -9.5288 -9.5093 -9.5093 -9.3194 -9.3194 -3.7916 -3.7916 -3.6432 -3.6432 -2.2472 -2.2472 -2.2189 -2.2189 1.6931 1.6931 1.8560 1.8560 2.1459 2.1459 2.2300 2.2300 2.9352 2.9352 2.9841 2.9841 3.2493 3.2493 3.4932 3.4932 5.3698 5.3698 6.6226 6.6226 7.8537 7.8537 9.0425 9.0425 10.5000 10.5001 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): -14.8422 -14.8422 -12.5981 -12.5981 -12.5981 -12.5981 -10.0947 -10.0947 -10.0947 -10.0947 -10.0112 -10.0112 -9.5339 -9.5339 -9.5048 -9.5048 -9.5048 -9.5048 -4.1608 -4.1608 -3.6727 -3.6727 -2.2864 -2.2864 -2.2864 -2.2864 2.0067 2.0067 2.0067 2.0067 2.1395 2.1395 2.9565 2.9565 2.9565 2.9565 3.2720 3.2720 3.5855 3.5855 3.5855 3.5855 5.8273 5.8273 6.2284 6.2284 7.4435 7.4435 7.4435 7.4435 9.9586 9.9586 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 6.0278 ev ! total energy = -264.35454904 Ry Harris-Foulkes estimate = -264.35454904 Ry estimated scf accuracy < 2.4E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -89.04190891 Ry hartree contribution = 64.57472590 Ry xc contribution = -64.94593138 Ry ewald contribution = -174.94143465 Ry smearing contrib. (-TS) = -0.00000000 Ry convergence has been achieved in 7 iterations Writing output data file CsPbBr3.save init_run : 3.29s CPU 1.84s WALL ( 1 calls) electrons : 53.17s CPU 29.27s WALL ( 1 calls) Called by init_run: wfcinit : 1.41s CPU 0.74s WALL ( 1 calls) potinit : 0.12s CPU 0.06s WALL ( 1 calls) Called by electrons: c_bands : 43.58s CPU 24.14s WALL ( 8 calls) sum_band : 5.50s CPU 2.98s WALL ( 8 calls) v_of_rho : 0.12s CPU 0.05s WALL ( 8 calls) v_h : 0.02s CPU 0.01s WALL ( 8 calls) v_xc : 0.10s CPU 0.05s WALL ( 8 calls) newd : 4.46s CPU 2.31s WALL ( 8 calls) mix_rho : 0.05s CPU 0.03s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.14s CPU 0.07s WALL ( 170 calls) cegterg : 41.98s CPU 23.33s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.39s CPU 0.20s WALL ( 80 calls) addusdens : 0.47s CPU 0.31s WALL ( 8 calls) Called by *egterg: h_psi : 24.52s CPU 13.12s WALL ( 462 calls) s_psi : 0.78s CPU 0.42s WALL ( 462 calls) g_psi : 0.05s CPU 0.02s WALL ( 372 calls) cdiaghg : 13.79s CPU 8.14s WALL ( 442 calls) cegterg:over : 1.55s CPU 0.88s WALL ( 372 calls) cegterg:upda : 1.22s CPU 0.69s WALL ( 372 calls) cegterg:last : 0.26s CPU 0.17s WALL ( 80 calls) cdiaghg:chol : 0.59s CPU 0.32s WALL ( 442 calls) cdiaghg:inve : 0.25s CPU 0.16s WALL ( 442 calls) cdiaghg:para : 0.84s CPU 0.52s WALL ( 884 calls) Called by h_psi: h_psi:vloc : 21.59s CPU 11.62s WALL ( 462 calls) h_psi:vnl : 2.91s CPU 1.48s WALL ( 462 calls) add_vuspsi : 1.78s CPU 0.88s WALL ( 462 calls) General routines calbec : 1.49s CPU 0.79s WALL ( 542 calls) fft : 0.32s CPU 0.18s WALL ( 154 calls) fftw : 24.58s CPU 13.25s WALL ( 56232 calls) Parallel routines fft_scatter : 20.58s CPU 11.11s WALL ( 56386 calls) PWSCF : 1m 0.41s CPU 0m36.64s WALL This run was terminated on: 21:35:33 20Mar2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=