Program PWSCF v.6.0 starts on 5Feb2017 at 20:45: 0 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 Serial multi-threaded version, running on 1 processor cores 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 G-vector sticks info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Sum 703 703 211 12387 12387 1989 bravais-lattice index = 2 lattice parameter (alat) = 9.2700 a.u. unit-cell volume = 199.1495 (a.u.)^3 number of atoms/cell = 1 number of atomic types = 1 number of electrons = 14.00 number of Kohn-Sham states= 11 kinetic-energy cutoff = 60.0000 Ry charge density cutoff = 240.0000 Ry convergence threshold = 1.0E-10 mixing beta = 0.7000 number of iterations used = 8 plain mixing Exchange-correlation = PZ ( 1 1 0 0 0 0) celldm(1)= 9.270000 celldm(2)= 0.000000 celldm(3)= 0.000000 celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000 crystal axes: (cart. coord. in units of alat) a(1) = ( -0.500000 0.000000 0.500000 ) a(2) = ( 0.000000 0.500000 0.500000 ) a(3) = ( -0.500000 0.500000 0.000000 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( -1.000000 -1.000000 1.000000 ) b(2) = ( 1.000000 1.000000 1.000000 ) b(3) = ( -1.000000 1.000000 -1.000000 ) PseudoPot. # 1 for Pb read from file: ./pb_s.UPF MD5 check sum: 7774f596943a2055f0d8d19eea6f2ae8 Pseudo is Norm-conserving, Zval = 14.0 Generated using "atomic" code by A. Dal Corso (Quantum ESPRESSO distribution) Using radial grid of 1281 points, 2 beta functions with: l(1) = 2 l(2) = 1 atomic species valence mass pseudopotential Pb 14.00 207.20000 Pb( 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) = ( 0 1 -1 ) ( 1 0 -1 ) ( 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 ) ( -1 0 1 ) ( -1 1 0 ) 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) = ( 0 -1 1 ) ( 0 -1 0 ) ( 1 -1 0 ) 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) = ( -1 0 1 ) ( 0 -1 1 ) ( 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 ) ( 0 1 -1 ) ( -1 1 0 ) 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) = ( 1 0 -1 ) ( 1 0 0 ) ( 1 -1 0 ) 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) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 -1 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) = ( 1 0 0 ) ( 1 -1 0 ) ( 1 0 -1 ) 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 1 -1 ) ( -1 1 0 ) ( 0 1 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 -1 1 ) ( 0 0 1 ) ( -1 0 1 ) 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 1 0 ) ( 0 1 0 ) ( 0 1 -1 ) 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) = ( 0 0 -1 ) ( 0 -1 0 ) ( -1 0 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) = ( 0 0 1 ) ( -1 0 1 ) ( 0 -1 1 ) 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 -1 0 ) ( 1 0 -1 ) ( 1 0 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) = ( -1 0 1 ) ( -1 1 0 ) ( -1 0 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 1 0 ) ( 0 0 1 ) ( 1 0 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) = ( 1 0 -1 ) ( 0 0 -1 ) ( 0 1 -1 ) 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 -1 0 ) ( 1 -1 0 ) ( 0 -1 1 ) 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 0 -1 ) ( 0 1 -1 ) ( 1 0 -1 ) 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) = ( -1 1 0 ) ( -1 0 0 ) ( -1 0 1 ) 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 0 1 ) ( 1 0 0 ) ( 0 1 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) = ( 1 -1 0 ) ( 0 -1 1 ) ( 0 -1 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) = ( 0 -1 1 ) ( -1 0 1 ) ( 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 ) ( 1 0 -1 ) ( 1 -1 0 ) 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) = ( 0 1 -1 ) ( 0 1 0 ) ( -1 1 0 ) 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) = ( 1 0 -1 ) ( 0 1 -1 ) ( 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 ) ( 0 -1 1 ) ( 1 -1 0 ) 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) = ( -1 0 1 ) ( -1 0 0 ) ( -1 1 0 ) 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) = ( 1 0 0 ) ( 0 0 1 ) ( 0 1 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) = ( -1 0 0 ) ( -1 1 0 ) ( -1 0 1 ) 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 -1 1 ) ( 1 -1 0 ) ( 0 -1 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 1 -1 ) ( 0 0 -1 ) ( 1 0 -1 ) 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 -1 0 ) ( 0 -1 0 ) ( 0 -1 1 ) 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) = ( 0 0 1 ) ( 0 1 0 ) ( 1 0 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) = ( 0 0 -1 ) ( 1 0 -1 ) ( 0 1 -1 ) 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 1 0 ) ( -1 0 1 ) ( -1 0 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) = ( 1 0 -1 ) ( 1 -1 0 ) ( 1 0 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 -1 0 ) ( 0 0 -1 ) ( -1 0 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) = ( -1 0 1 ) ( 0 0 1 ) ( 0 -1 1 ) 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 1 0 ) ( -1 1 0 ) ( 0 1 -1 ) 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 0 1 ) ( 0 -1 1 ) ( -1 0 1 ) 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) = ( 1 -1 0 ) ( 1 0 0 ) ( 1 0 -1 ) 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 0 -1 ) ( -1 0 0 ) ( 0 -1 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) = ( -1 1 0 ) ( 0 1 -1 ) ( 0 1 0 ) cart. s(48) = ( 0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) point group O_h (m-3m) there are 10 classes the character table: E 8C3 3C2 6C4 6C2' i 8S6 3s_h 6S4 6s_d A_1g 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 A_2g 1.00 1.00 1.00 -1.00 -1.00 1.00 1.00 1.00 -1.00 -1.00 E_g 2.00 -1.00 2.00 0.00 0.00 2.00 -1.00 2.00 0.00 0.00 T_1g 3.00 0.00 -1.00 1.00 -1.00 3.00 0.00 -1.00 1.00 -1.00 T_2g 3.00 0.00 -1.00 -1.00 1.00 3.00 0.00 -1.00 -1.00 1.00 A_1u 1.00 1.00 1.00 1.00 1.00 -1.00 -1.00 -1.00 -1.00 -1.00 A_2u 1.00 1.00 1.00 -1.00 -1.00 -1.00 -1.00 -1.00 1.00 1.00 E_u 2.00 -1.00 2.00 0.00 0.00 -2.00 1.00 -2.00 0.00 0.00 T_1u 3.00 0.00 -1.00 1.00 -1.00 -3.00 0.00 1.00 -1.00 1.00 T_2u 3.00 0.00 -1.00 -1.00 1.00 -3.00 0.00 1.00 1.00 -1.00 the symmetry operations in each class and the name of the first element: E 1 identity 8C3 17 19 20 18 24 21 22 23 120 deg rotation - cart. axis [-1,-1,-1] 3C2 2 4 3 180 deg rotation - cart. axis [0,0,1] 6C4 7 8 15 16 12 11 90 deg rotation - cart. axis [0,0,-1] 6C2' 5 6 14 13 10 9 180 deg rotation - cart. axis [1,1,0] i 25 inversion 8S6 41 43 44 42 48 45 46 47 inv. 120 deg rotation - cart. axis [-1,-1,-1] 3s_h 26 28 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] 6s_d 29 30 38 37 34 33 inv. 180 deg rotation - cart. axis [1,1,0] Cartesian axes site n. atom positions (alat units) 1 Pb tau( 1) = ( 0.0000000 0.0000000 0.0000000 ) Crystallographic axes site n. atom positions (cryst. coord.) 1 Pb tau( 1) = ( 0.0000000 0.0000000 0.0000000 ) number of k points= 4 Methfessel-Paxton smearing, width (Ry)= 0.0250 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0740741 k( 2) = ( -0.3333333 0.3333333 -0.3333333), wk = 0.5925926 k( 3) = ( 0.0000000 0.6666667 0.0000000), wk = 0.4444444 k( 4) = ( 0.6666667 -0.0000000 0.6666667), wk = 0.8888889 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0740741 k( 2) = ( 0.0000000 0.0000000 0.3333333), wk = 0.5925926 k( 3) = ( 0.0000000 0.3333333 0.3333333), wk = 0.4444444 k( 4) = ( 0.0000000 0.3333333 -0.3333333), wk = 0.8888889 Dense grid: 12387 G-vectors FFT dimensions: ( 36, 36, 36) Estimated max dynamical RAM per process > 10.78MB Initial potential from superposition of free atoms starting charge 13.99923, renormalised to 14.00000 Starting wfc are 9 randomized atomic wfcs + 2 random wfc total cpu time spent up to now is 0.4 secs per-process dynamical memory: 8.1 Mb Self-consistent Calculation iteration # 1 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 4.2 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 6.74E-05, avg # of iterations = 1.8 total cpu time spent up to now is 0.9 secs total energy = -119.53911991 Ry Harris-Foulkes estimate = -119.54461162 Ry estimated scf accuracy < 0.00956363 Ry iteration # 2 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.83E-05, avg # of iterations = 3.2 total cpu time spent up to now is 1.1 secs total energy = -119.54067050 Ry Harris-Foulkes estimate = -119.54090760 Ry estimated scf accuracy < 0.00048537 Ry iteration # 3 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.47E-06, avg # of iterations = 2.2 total cpu time spent up to now is 1.3 secs total energy = -119.54073374 Ry Harris-Foulkes estimate = -119.54073978 Ry estimated scf accuracy < 0.00002247 Ry iteration # 4 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.60E-07, avg # of iterations = 1.0 total cpu time spent up to now is 1.5 secs total energy = -119.54073497 Ry Harris-Foulkes estimate = -119.54073496 Ry estimated scf accuracy < 0.00000011 Ry iteration # 5 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.56E-10, avg # of iterations = 2.2 total cpu time spent up to now is 1.8 secs total energy = -119.54073501 Ry Harris-Foulkes estimate = -119.54073503 Ry estimated scf accuracy < 0.00000002 Ry iteration # 6 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.27E-10, avg # of iterations = 2.0 total cpu time spent up to now is 2.0 secs total energy = -119.54073502 Ry Harris-Foulkes estimate = -119.54073502 Ry estimated scf accuracy < 1.8E-09 Ry iteration # 7 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.32E-11, avg # of iterations = 2.0 total cpu time spent up to now is 2.2 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1471 PWs) bands (ev): -5.5635 -5.5635 -5.5635 -5.4022 -5.4022 -0.5737 18.6327 18.6327 18.6327 21.9477 22.9972 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 k =-0.3333 0.3333-0.3333 ( 1540 PWs) bands (ev): -5.5864 -5.5317 -5.5317 -5.3154 -5.3154 1.5754 9.7460 17.2082 17.2082 20.8948 22.5624 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.6763 -0.0000 -0.0000 -0.0000 -0.0000 k = 0.0000 0.6667 0.0000 ( 1580 PWs) bands (ev): -5.7000 -5.5748 -5.3124 -5.3124 -5.3002 2.2921 11.2650 14.2103 14.2103 18.1260 18.9957 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 k = 0.6667-0.0000 0.6667 ( 1568 PWs) bands (ev): -5.6238 -5.4488 -5.4004 -5.3676 -5.3022 4.1691 9.4624 9.7032 14.9408 18.7583 25.8659 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0294 0.7698 -0.0000 -0.0000 -0.0000 the Fermi energy is 9.8187 ev ! total energy = -119.54073502 Ry Harris-Foulkes estimate = -119.54073502 Ry estimated scf accuracy < 5.0E-12 Ry The total energy is the sum of the following terms: one-electron contribution = -32.99451180 Ry hartree contribution = 24.93109769 Ry xc contribution = -14.53469691 Ry ewald contribution = -96.93991011 Ry smearing contrib. (-TS) = -0.00271389 Ry convergence has been achieved in 7 iterations Forces acting on atoms (cartesian axes, Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 The non-local contrib. to forces atom 1 type 1 force = -0.00000000 -0.00000000 0.00000000 The ionic contribution to forces atom 1 type 1 force = -0.00000000 -0.00000000 -0.00000000 The local contribution to forces atom 1 type 1 force = 0.00000000 -0.00000000 -0.00000000 The core correction contribution to forces atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 The Hubbard contrib. to forces atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 The SCF correction term to forces atom 1 type 1 force = 0.00000000 -0.00000000 -0.00000000 Total force = 0.000000 Total SCF correction = 0.000000 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= -20.78 -0.00014127 -0.00000000 -0.00000000 -20.78 -0.00 -0.00 -0.00000000 -0.00014127 0.00000000 -0.00 -20.78 0.00 -0.00000000 0.00000000 -0.00014127 -0.00 0.00 -20.78 kinetic stress (kbar) 34581.92 -0.00 0.00 -0.00 34581.92 -0.00 0.00 -0.00 34581.92 local stress (kbar) 32811.08 0.00 -0.00 0.00 32811.08 -0.00 -0.00 -0.00 32811.08 nonloc. stress (kbar) -46353.32 0.00 0.00 0.00 -46353.32 -0.00 0.00 -0.00 -46353.32 hartree stress (kbar) 6138.59 -0.00 0.00 -0.00 6138.59 0.00 0.00 0.00 6138.59 exc-cor stress (kbar) -3330.29 0.00 0.00 0.00 -3330.29 0.00 0.00 0.00 -3330.29 corecor stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ewald stress (kbar) -23868.75 -0.00 0.00 -0.00 -23868.75 0.00 0.00 0.00 -23868.75 hubbard stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 london stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 XDM stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 dft-nl stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 TS-vdW stress (kbar) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 Writing output data file pb.save init_run : 0.27s CPU 0.28s WALL ( 1 calls) electrons : 1.82s CPU 1.85s WALL ( 1 calls) forces : 0.02s CPU 0.02s WALL ( 1 calls) stress : 0.16s CPU 0.16s WALL ( 1 calls) Called by init_run: wfcinit : 0.07s CPU 0.07s WALL ( 1 calls) wfcinit:atom : 0.00s CPU 0.00s WALL ( 4 calls) wfcinit:wfcr : 0.07s CPU 0.07s WALL ( 4 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 1.46s CPU 1.48s WALL ( 8 calls) sum_band : 0.29s CPU 0.30s WALL ( 8 calls) v_of_rho : 0.06s CPU 0.06s WALL ( 8 calls) v_h : 0.01s CPU 0.01s WALL ( 8 calls) v_xc : 0.05s CPU 0.05s WALL ( 8 calls) mix_rho : 0.01s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.05s CPU 0.05s WALL ( 76 calls) cegterg : 1.43s CPU 1.44s WALL ( 32 calls) Called by sum_band: Called by *egterg: h_psi : 1.28s CPU 1.29s WALL ( 111 calls) g_psi : 0.01s CPU 0.01s WALL ( 75 calls) cdiaghg : 0.02s CPU 0.02s WALL ( 103 calls) cegterg:over : 0.07s CPU 0.07s WALL ( 75 calls) cegterg:upda : 0.05s CPU 0.05s WALL ( 75 calls) cegterg:last : 0.03s CPU 0.03s WALL ( 32 calls) Called by h_psi: h_psi:pot : 1.28s CPU 1.28s WALL ( 111 calls) h_psi:calbec : 0.02s CPU 0.02s WALL ( 111 calls) vloc_psi : 1.23s CPU 1.24s WALL ( 111 calls) add_vuspsi : 0.02s CPU 0.02s WALL ( 111 calls) General routines calbec : 0.03s CPU 0.03s WALL ( 131 calls) fft : 0.04s CPU 0.04s WALL ( 38 calls) fftw : 1.31s CPU 1.32s WALL ( 2170 calls) PWSCF : 2.38s CPU 2.50s WALL This run was terminated on: 20:45: 2 5Feb2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=