Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 17:48: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 file Rb.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S 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 36 10 1961 976 150 Max 59 37 11 1964 995 153 Sum 2077 1305 373 70655 35513 5449 bravais-lattice index = 14 lattice parameter (alat) = 9.0518 a.u. unit-cell volume = 741.6569 (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 = 50.0000 Ry charge density cutoff = 317.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.051787 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 Rb read from file: /users/gautes/Pseudo/Rb.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 948b8f0a9070089a35782939cf30e963 Pseudo is Ultrasoft + core correction, Zval = 9.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1219 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 atomic species valence mass pseudopotential F 7.00 18.99840 F( 1.00) Pb 14.00 207.20000 Pb( 1.00) Rb 9.00 85.46780 Rb( 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: 70655 G-vectors FFT dimensions: ( 54, 54, 54) Smooth grid: 35513 G-vectors FFT dimensions: ( 45, 45, 45) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.21 Mb ( 264, 52) NL pseudopotentials 0.22 Mb ( 132, 110) Each V/rho on FFT grid 0.09 Mb ( 5832) Each G-vector array 0.01 Mb ( 1962) G-vector shells 0.00 Mb ( 423) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.84 Mb ( 264, 208) Each subspace H/S matrix 0.04 Mb ( 52, 52) Each matrix 0.17 Mb ( 110, 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.97612, renormalised to 44.00000 Starting wfc are 58 randomized atomic wfcs total cpu time spent up to now is 3.2 secs per-process dynamical memory: 48.3 Mb Self-consistent Calculation iteration # 1 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 3.3 total cpu time spent up to now is 4.4 secs total energy = -340.32081850 Ry Harris-Foulkes estimate = -340.97669484 Ry estimated scf accuracy < 0.83633997 Ry iteration # 2 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.90E-03, avg # of iterations = 3.6 total cpu time spent up to now is 5.8 secs total energy = -340.40070459 Ry Harris-Foulkes estimate = -341.10113797 Ry estimated scf accuracy < 1.61467268 Ry iteration # 3 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.90E-03, avg # of iterations = 2.0 total cpu time spent up to now is 6.9 secs total energy = -340.73289202 Ry Harris-Foulkes estimate = -340.73844236 Ry estimated scf accuracy < 0.01320801 Ry iteration # 4 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.00E-05, avg # of iterations = 5.4 total cpu time spent up to now is 8.5 secs total energy = -340.73660450 Ry Harris-Foulkes estimate = -340.73694563 Ry estimated scf accuracy < 0.00078107 Ry iteration # 5 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.78E-06, avg # of iterations = 4.9 total cpu time spent up to now is 9.9 secs total energy = -340.73675460 Ry Harris-Foulkes estimate = -340.73675500 Ry estimated scf accuracy < 0.00000899 Ry iteration # 6 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.04E-08, avg # of iterations = 3.6 total cpu time spent up to now is 11.3 secs total energy = -340.73675809 Ry Harris-Foulkes estimate = -340.73675783 Ry estimated scf accuracy < 0.00000092 Ry iteration # 7 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.10E-09, avg # of iterations = 2.1 total cpu time spent up to now is 12.6 secs total energy = -340.73675826 Ry Harris-Foulkes estimate = -340.73675828 Ry estimated scf accuracy < 0.00000005 Ry iteration # 8 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.03E-10, avg # of iterations = 2.7 total cpu time spent up to now is 13.9 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 4385 PWs) bands (ev): -21.7548 -21.7548 -17.6291 -17.6291 -17.4910 -17.4910 -17.4910 -17.4910 -12.2471 -12.2471 -12.2471 -12.2471 -9.7919 -9.7919 -9.6255 -9.6255 -9.6255 -9.6255 -6.4122 -6.4122 -5.5342 -5.5342 -5.5342 -5.5342 -0.7504 -0.7504 -0.6582 -0.6582 -0.3066 -0.3066 -0.3066 -0.3066 1.2614 1.2614 1.2731 1.2731 1.2731 1.2731 1.3020 1.3020 1.3278 1.3278 1.3278 1.3278 8.7923 8.7923 10.1504 10.1504 10.4113 10.4113 10.4113 10.4113 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 ( 4428 PWs) bands (ev): -21.7533 -21.7533 -17.6114 -17.6114 -17.4907 -17.4907 -17.4787 -17.4787 -12.2792 -12.2792 -12.2477 -12.2477 -9.7950 -9.7950 -9.6925 -9.6925 -9.6261 -9.6261 -6.4124 -6.4124 -5.5419 -5.5419 -5.5261 -5.5261 -1.6831 -1.6831 -0.5254 -0.5254 -0.2841 -0.2841 0.5992 0.5992 1.2485 1.2485 1.2560 1.2560 1.2684 1.2684 1.2846 1.2846 1.3141 1.3141 1.6328 1.6328 8.7511 8.7511 9.9664 9.9664 10.6121 10.6121 10.6432 10.6432 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 ( 4432 PWs) bands (ev): -21.7508 -21.7508 -17.5909 -17.5909 -17.4901 -17.4901 -17.4492 -17.4492 -12.3342 -12.3342 -12.2485 -12.2485 -9.7999 -9.7999 -9.7976 -9.7976 -9.6268 -9.6268 -6.4132 -6.4132 -5.5537 -5.5537 -5.5130 -5.5130 -2.2148 -2.2148 -0.4804 -0.4804 -0.2476 -0.2476 0.7917 0.7917 1.2056 1.2056 1.2324 1.2324 1.2452 1.2452 1.2666 1.2666 1.2862 1.2862 2.8017 2.8017 8.0459 8.0459 10.0824 10.0824 10.9726 10.9726 11.0385 11.0385 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 ( 4419 PWs) bands (ev): -21.7518 -21.7518 -17.5915 -17.5915 -17.4832 -17.4832 -17.4752 -17.4752 -12.2951 -12.2951 -12.2640 -12.2640 -9.7979 -9.7979 -9.7233 -9.7233 -9.6608 -9.6608 -6.4151 -6.4151 -5.5469 -5.5469 -5.5271 -5.5271 -2.2168 -2.2168 -0.4000 -0.4000 0.0831 0.0831 0.7135 0.7135 1.1665 1.1665 1.2518 1.2518 1.2657 1.2657 1.2897 1.2897 1.3809 1.3809 1.9980 1.9980 8.8202 8.8202 9.6444 9.6444 10.6525 10.6525 10.8548 10.8548 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 ( 4419 PWs) bands (ev): -21.7495 -21.7495 -17.5687 -17.5687 -17.4802 -17.4802 -17.4495 -17.4495 -12.3443 -12.3443 -12.2700 -12.2700 -9.8145 -9.8145 -9.8023 -9.8023 -9.6728 -9.6728 -6.4197 -6.4197 -5.5618 -5.5618 -5.5212 -5.5212 -2.6809 -2.6809 -0.3390 -0.3390 0.3110 0.3110 0.7844 0.7844 1.1128 1.1128 1.2307 1.2307 1.2413 1.2413 1.2695 1.2695 1.4084 1.4084 3.0606 3.0606 8.1502 8.1502 9.7069 9.7069 10.7929 10.7929 11.1981 11.1981 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 ( 4426 PWs) bands (ev): -21.7475 -21.7475 -17.5434 -17.5434 -17.4556 -17.4556 -17.4466 -17.4466 -12.3758 -12.3758 -12.2930 -12.2930 -9.8661 -9.8661 -9.8067 -9.8067 -9.7204 -9.7204 -6.4306 -6.4306 -5.5822 -5.5822 -5.5260 -5.5260 -3.1291 -3.1291 -0.2614 -0.2614 0.7532 0.7532 0.9470 0.9470 1.0925 1.0925 1.2016 1.2016 1.2124 1.2124 1.2810 1.2810 1.5283 1.5283 3.7997 3.7997 7.9985 7.9985 9.1539 9.1539 10.7068 10.7068 11.0381 11.0381 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 ( 4419 PWs) bands (ev): -21.7505 -21.7505 -17.5686 -17.5686 -17.4765 -17.4765 -17.4733 -17.4733 -12.2957 -12.2957 -12.2953 -12.2953 -9.8008 -9.8008 -9.7238 -9.7238 -9.7237 -9.7237 -6.4196 -6.4196 -5.5496 -5.5496 -5.5364 -5.5364 -2.6410 -2.6410 -0.0065 -0.0065 0.1866 0.1866 0.7475 0.7475 1.1254 1.1254 1.1429 1.1429 1.3450 1.3450 1.3811 1.3811 1.3825 1.3825 2.3387 2.3387 8.8376 8.8376 9.4952 9.4952 10.1923 10.1923 11.2275 11.2275 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 ( 4413 PWs) bands (ev): -21.7484 -21.7484 -17.5423 -17.5423 -17.4744 -17.4744 -17.4494 -17.4494 -12.3497 -12.3497 -12.2962 -12.2962 -9.8232 -9.8232 -9.8056 -9.8056 -9.7245 -9.7245 -6.4273 -6.4273 -5.5671 -5.5671 -5.5372 -5.5372 -3.0629 -3.0629 0.1247 0.1247 0.4746 0.4746 0.7968 0.7968 1.0410 1.0410 1.1389 1.1389 1.3407 1.3407 1.3649 1.3649 1.4322 1.4322 3.3003 3.3003 8.1247 8.1247 9.5110 9.5110 10.1979 10.1979 11.4687 11.4687 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 ( 4447 PWs) bands (ev): -21.7465 -21.7465 -17.5118 -17.5118 -17.4547 -17.4547 -17.4459 -17.4459 -12.3763 -12.3763 -12.3240 -12.3240 -9.8667 -9.8667 -9.8105 -9.8105 -9.7772 -9.7772 -6.4402 -6.4402 -5.5860 -5.5860 -5.5512 -5.5512 -3.4686 -3.4686 0.3302 0.3302 0.8021 0.8021 0.9421 0.9421 1.0011 1.0011 1.1466 1.1466 1.3661 1.3661 1.4026 1.4026 1.5307 1.5307 3.9992 3.9992 7.7796 7.7796 8.9975 8.9975 9.8636 9.8636 11.1342 11.1342 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 ( 4428 PWs) bands (ev): -21.7449 -21.7449 -17.4649 -17.4649 -17.4491 -17.4491 -17.4433 -17.4433 -12.3772 -12.3772 -12.3771 -12.3771 -9.8679 -9.8679 -9.8677 -9.8677 -9.8164 -9.8164 -6.4561 -6.4561 -5.5934 -5.5934 -5.5891 -5.5891 -3.8504 -3.8504 0.8730 0.8730 0.9278 0.9278 0.9384 0.9384 1.1065 1.1065 1.1325 1.1325 1.5015 1.5015 1.5320 1.5320 1.5451 1.5451 4.5978 4.5978 7.2751 7.2751 8.7082 8.7082 8.8099 8.8099 11.6567 11.6567 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 5.4091 ev ! total energy = -340.73675827 Ry Harris-Foulkes estimate = -340.73675827 Ry estimated scf accuracy < 1.8E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -160.60388941 Ry hartree contribution = 105.11634426 Ry xc contribution = -80.54216903 Ry ewald contribution = -204.70704409 Ry smearing contrib. (-TS) = 0.00000000 Ry convergence has been achieved in 8 iterations Writing output data file RbPbF3.save init_run : 0.58s CPU 0.68s WALL ( 1 calls) electrons : 10.42s CPU 10.66s WALL ( 1 calls) Called by init_run: wfcinit : 0.33s CPU 0.36s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 8.22s CPU 8.40s WALL ( 8 calls) sum_band : 1.63s CPU 1.67s 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 : 0.52s CPU 0.55s WALL ( 9 calls) mix_rho : 0.02s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.03s CPU 0.02s WALL ( 170 calls) cegterg : 7.96s CPU 8.10s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.17s CPU 0.19s WALL ( 80 calls) addusdens : 0.41s CPU 0.42s WALL ( 8 calls) Called by *egterg: h_psi : 5.36s CPU 5.47s WALL ( 366 calls) s_psi : 0.19s CPU 0.22s WALL ( 366 calls) g_psi : 0.02s CPU 0.01s WALL ( 276 calls) cdiaghg : 1.91s CPU 1.99s WALL ( 356 calls) cegterg:over : 0.26s CPU 0.23s WALL ( 276 calls) cegterg:upda : 0.23s CPU 0.21s WALL ( 276 calls) cegterg:last : 0.05s CPU 0.07s WALL ( 80 calls) cdiaghg:chol : 0.11s CPU 0.12s WALL ( 356 calls) cdiaghg:inve : 0.06s CPU 0.07s WALL ( 356 calls) cdiaghg:para : 0.09s CPU 0.12s WALL ( 712 calls) Called by h_psi: h_psi:vloc : 4.80s CPU 4.93s WALL ( 366 calls) h_psi:vnl : 0.54s CPU 0.53s WALL ( 366 calls) add_vuspsi : 0.28s CPU 0.25s WALL ( 366 calls) General routines calbec : 0.34s CPU 0.35s WALL ( 446 calls) fft : 0.08s CPU 0.08s WALL ( 263 calls) ffts : 0.02s CPU 0.01s WALL ( 68 calls) fftw : 5.38s CPU 5.47s WALL ( 54792 calls) interpolate : 0.04s CPU 0.04s WALL ( 68 calls) Parallel routines fft_scatter : 1.86s CPU 1.86s WALL ( 55123 calls) PWSCF : 13.63s CPU 15.12s WALL This run was terminated on: 17:48:35 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=