Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 17:32:20 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 F.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 2S renormalized 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 57 57 16 1971 1971 288 Max 58 58 17 1976 1976 290 Sum 2077 2077 577 71015 71015 10395 bravais-lattice index = 14 lattice parameter (alat) = 9.0231 a.u. unit-cell volume = 734.6188 (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)= 9.023064 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 F read from file: /users/gautes/Pseudo/F.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: ce20c52fa94def3fe434d52ed9f7ea44 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 1105 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. # 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 F 7.00 18.99840 F( 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.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 Dense grid: 71015 G-vectors FFT dimensions: ( 54, 54, 54) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.41 Mb ( 518, 52) NL pseudopotentials 0.55 Mb ( 259, 138) Each V/rho on FFT grid 0.09 Mb ( 5832) Each G-vector array 0.02 Mb ( 1974) G-vector shells 0.00 Mb ( 432) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 1.64 Mb ( 518, 208) Each subspace H/S matrix 0.04 Mb ( 52, 52) Each matrix 0.22 Mb ( 138, 2, 52) Arrays for rho mixing 0.71 Mb ( 5832, 8) Check: negative/imaginary core charge= -0.000002 0.000000 Initial potential from superposition of free atoms starting charge 43.96078, renormalised to 44.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 3.0 secs per-process dynamical memory: 64.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.7 total cpu time spent up to now is 5.1 secs total energy = -320.58027817 Ry Harris-Foulkes estimate = -321.27097285 Ry estimated scf accuracy < 0.87320568 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.98E-03, avg # of iterations = 3.8 total cpu time spent up to now is 7.2 secs total energy = -320.65999324 Ry Harris-Foulkes estimate = -321.37131542 Ry estimated scf accuracy < 1.61981331 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.98E-03, avg # of iterations = 2.0 total cpu time spent up to now is 9.0 secs total energy = -321.00285465 Ry Harris-Foulkes estimate = -321.00891497 Ry estimated scf accuracy < 0.01617389 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.68E-05, avg # of iterations = 6.0 total cpu time spent up to now is 11.3 secs total energy = -321.00686117 Ry Harris-Foulkes estimate = -321.00681739 Ry estimated scf accuracy < 0.00014570 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.31E-07, avg # of iterations = 5.3 total cpu time spent up to now is 14.0 secs total energy = -321.00694430 Ry Harris-Foulkes estimate = -321.00694723 Ry estimated scf accuracy < 0.00001474 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.35E-08, avg # of iterations = 2.1 total cpu time spent up to now is 15.8 secs total energy = -321.00694867 Ry Harris-Foulkes estimate = -321.00694787 Ry estimated scf accuracy < 0.00000120 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.72E-09, avg # of iterations = 2.8 total cpu time spent up to now is 17.9 secs total energy = -321.00694899 Ry Harris-Foulkes estimate = -321.00694898 Ry estimated scf accuracy < 0.00000002 Ry iteration # 8 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.18E-11, avg # of iterations = 4.0 total cpu time spent up to now is 20.5 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 8829 PWs) bands (ev): -16.5947 -16.5947 -16.4205 -16.4205 -16.4205 -16.4205 -13.7121 -13.7121 -11.1242 -11.1242 -11.1242 -11.1242 -8.6718 -8.6718 -8.4995 -8.4995 -8.4995 -8.4995 -2.5206 -2.5206 -1.1518 -1.1518 -1.1518 -1.1518 0.3927 0.3927 0.4660 0.4660 0.8926 0.8926 0.8926 0.8926 2.3234 2.3234 2.3234 2.3234 2.3517 2.3517 2.4197 2.4197 2.5498 2.5498 2.5498 2.5498 10.6324 10.6324 10.6324 10.6324 10.9176 10.9176 11.3489 11.3489 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.2000 ( 8878 PWs) bands (ev): -16.5684 -16.5684 -16.4200 -16.4200 -16.4093 -16.4093 -13.7163 -13.7163 -11.1574 -11.1574 -11.1247 -11.1247 -8.6750 -8.6750 -8.5695 -8.5695 -8.5000 -8.5000 -2.5164 -2.5164 -1.2298 -1.2298 -1.0928 -1.0928 -0.4327 -0.4327 0.5920 0.5920 0.8676 0.8676 1.7787 1.7787 2.2922 2.2922 2.3049 2.3049 2.3310 2.3310 2.4269 2.4269 2.5148 2.5148 2.7743 2.7743 10.2877 10.2877 10.8454 10.8454 11.2493 11.2493 11.4795 11.4795 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.4000 ( 8866 PWs) bands (ev): -16.5291 -16.5291 -16.4194 -16.4194 -16.3866 -16.3866 -13.7231 -13.7231 -11.2146 -11.2146 -11.1255 -11.1255 -8.6803 -8.6803 -8.6791 -8.6791 -8.5008 -8.5008 -2.5169 -2.5169 -1.4185 -1.4185 -0.9885 -0.9885 -0.7423 -0.7423 0.5712 0.5712 0.8111 0.8111 1.9220 1.9220 2.2383 2.2383 2.2836 2.2836 2.3108 2.3108 2.3777 2.3777 2.4580 2.4580 3.9411 3.9411 9.3681 9.3681 11.1111 11.1111 11.2276 11.2276 11.7939 11.7939 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.2000 0.2000 ( 8874 PWs) bands (ev): -16.5424 -16.5424 -16.4135 -16.4135 -16.4044 -16.4044 -13.7190 -13.7190 -11.1737 -11.1737 -11.1416 -11.1416 -8.6780 -8.6780 -8.6014 -8.6014 -8.5363 -8.5363 -2.5218 -2.5218 -1.3646 -1.3646 -1.1138 -1.1138 -0.8175 -0.8175 0.7047 0.7047 1.2479 1.2479 1.8352 1.8352 2.2269 2.2269 2.2967 2.2967 2.3720 2.3720 2.4214 2.4214 2.5480 2.5480 3.1242 3.1242 10.2404 10.2404 10.9126 10.9126 11.3323 11.3323 11.6444 11.6444 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.2000 0.4000 ( 8869 PWs) bands (ev): -16.5050 -16.5050 -16.4096 -16.4096 -16.3837 -16.3837 -13.7232 -13.7232 -11.2248 -11.2248 -11.1479 -11.1479 -8.6964 -8.6964 -8.6827 -8.6827 -8.5489 -8.5489 -2.5367 -2.5367 -1.6812 -1.6812 -1.0868 -1.0868 -0.9753 -0.9753 0.7162 0.7162 1.4291 1.4291 1.8971 1.8971 2.2311 2.2311 2.2751 2.2751 2.3513 2.3513 2.3784 2.3784 2.5127 2.5127 4.2031 4.2031 9.4672 9.4672 10.8937 10.8937 11.5341 11.5341 11.7713 11.7713 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.4000 0.4000 ( 8890 PWs) bands (ev): -16.4758 -16.4758 -16.3857 -16.3857 -16.3757 -16.3757 -13.7235 -13.7235 -11.2569 -11.2569 -11.1719 -11.1719 -8.7494 -8.7494 -8.6870 -8.6870 -8.5988 -8.5988 -2.5648 -2.5648 -2.0716 -2.0716 -1.2530 -1.2530 -1.0028 -1.0028 0.7903 0.7903 1.8460 1.8460 2.1341 2.1341 2.2383 2.2383 2.2616 2.2616 2.3037 2.3037 2.4379 2.4379 2.5999 2.5999 4.9462 4.9462 9.2739 9.2739 10.5511 10.5511 11.7130 11.7130 11.9524 11.9525 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.2000 0.2000 0.2000 ( 8874 PWs) bands (ev): -16.5154 -16.5154 -16.4056 -16.4056 -16.4023 -16.4023 -13.7198 -13.7198 -11.1743 -11.1743 -11.1741 -11.1741 -8.6811 -8.6811 -8.6021 -8.6021 -8.6018 -8.6018 -2.5360 -2.5360 -1.6439 -1.6439 -1.1698 -1.1698 -0.9563 -0.9563 1.1117 1.1117 1.3238 1.3238 1.8412 1.8412 2.2398 2.2398 2.2456 2.2456 2.4346 2.4346 2.4735 2.4735 2.5433 2.5433 3.4699 3.4699 10.0970 10.0970 11.0443 11.0443 11.0572 11.0572 11.8263 11.8263 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.2000 0.2000 0.4000 ( 8884 PWs) bands (ev): -16.4772 -16.4772 -16.4034 -16.4034 -16.3815 -16.3815 -13.7212 -13.7212 -11.2304 -11.2304 -11.1750 -11.1750 -8.7054 -8.7054 -8.6861 -8.6861 -8.6027 -8.6027 -2.5648 -2.5648 -2.0144 -2.0144 -1.1786 -1.1786 -1.0621 -1.0621 1.1924 1.1924 1.6034 1.6034 1.9291 1.9291 2.2361 2.2361 2.2681 2.2681 2.4109 2.4109 2.4661 2.4661 2.5114 2.5114 4.4449 4.4449 9.3537 9.3537 10.5872 10.5872 11.4204 11.4204 11.9516 11.9516 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.2000 0.4000 0.4000 ( 8876 PWs) bands (ev): -16.4432 -16.4432 -16.3842 -16.3842 -16.3747 -16.3747 -13.7180 -13.7180 -11.2577 -11.2577 -11.2039 -11.2039 -8.7503 -8.7503 -8.6911 -8.6911 -8.6578 -8.6578 -2.6071 -2.6071 -2.4072 -2.4072 -1.2916 -1.2916 -1.1395 -1.1395 1.3982 1.3982 1.9876 1.9876 2.1558 2.1558 2.2508 2.2508 2.2861 2.2861 2.4314 2.4314 2.4918 2.4918 2.5995 2.5995 5.1447 5.1447 8.9130 8.9130 10.0994 10.0994 10.8082 10.8082 12.2409 12.2409 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.4000 0.4000 0.4000 ( 8889 PWs) bands (ev): -16.3944 -16.3944 -16.3777 -16.3777 -16.3716 -16.3716 -13.7091 -13.7091 -11.2591 -11.2591 -11.2590 -11.2590 -8.7520 -8.7520 -8.7516 -8.7516 -8.6977 -8.6977 -2.8078 -2.8078 -2.6538 -2.6538 -1.3494 -1.3494 -1.3214 -1.3214 2.1661 2.1661 2.1752 2.1752 2.2109 2.2109 2.2980 2.2980 2.3406 2.3406 2.5643 2.5643 2.5985 2.5985 2.6078 2.6078 5.7423 5.7423 8.2708 8.2708 9.5107 9.5107 9.6994 9.6994 13.2198 13.2198 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 7.6466 ev ! total energy = -321.00694901 Ry Harris-Foulkes estimate = -321.00694901 Ry estimated scf accuracy < 3.9E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -149.36767661 Ry hartree contribution = 101.25993747 Ry xc contribution = -67.54050559 Ry ewald contribution = -205.35870428 Ry smearing contrib. (-TS) = 0.00000000 Ry convergence has been achieved in 8 iterations Writing output data file CsPbF3.save init_run : 0.90s CPU 1.27s WALL ( 1 calls) electrons : 16.86s CPU 17.74s WALL ( 1 calls) Called by init_run: wfcinit : 0.47s CPU 0.63s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 13.82s CPU 14.31s WALL ( 8 calls) sum_band : 1.97s CPU 1.99s WALL ( 8 calls) v_of_rho : 0.04s CPU 0.04s WALL ( 9 calls) v_h : 0.00s CPU 0.00s WALL ( 9 calls) v_xc : 0.04s CPU 0.04s WALL ( 9 calls) newd : 1.10s CPU 1.11s WALL ( 9 calls) mix_rho : 0.01s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.06s CPU 0.07s WALL ( 170 calls) cegterg : 12.81s CPU 13.02s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.12s CPU 0.12s WALL ( 80 calls) addusdens : 0.22s CPU 0.23s WALL ( 8 calls) Called by *egterg: h_psi : 9.11s CPU 9.13s WALL ( 397 calls) s_psi : 0.29s CPU 0.33s WALL ( 397 calls) g_psi : 0.03s CPU 0.03s WALL ( 307 calls) cdiaghg : 2.42s CPU 2.42s WALL ( 387 calls) cegterg:over : 0.50s CPU 0.49s WALL ( 307 calls) cegterg:upda : 0.42s CPU 0.48s WALL ( 307 calls) cegterg:last : 0.10s CPU 0.13s WALL ( 81 calls) cdiaghg:chol : 0.12s CPU 0.15s WALL ( 387 calls) cdiaghg:inve : 0.10s CPU 0.09s WALL ( 387 calls) cdiaghg:para : 0.12s CPU 0.15s WALL ( 774 calls) Called by h_psi: h_psi:vloc : 7.96s CPU 8.03s WALL ( 397 calls) h_psi:vnl : 1.11s CPU 1.07s WALL ( 397 calls) add_vuspsi : 0.64s CPU 0.58s WALL ( 397 calls) General routines calbec : 0.61s CPU 0.63s WALL ( 477 calls) fft : 0.04s CPU 0.06s WALL ( 167 calls) fftw : 8.86s CPU 8.93s WALL ( 58104 calls) Parallel routines fft_scatter : 3.24s CPU 3.09s WALL ( 58271 calls) PWSCF : 20.02s CPU 23.10s WALL This run was terminated on: 17:32:43 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=